[mapserver-users] many raster projections

[Andy Colson]

To take OpenLayers out of the equation, I made myself a test.map file 
that just included two layers and used shp2img.

If I switch on the 2007 (raster07) then I get no arial (the result image 
is just a gray slate).

If I switch on raster09_chil, the result image has the arial's as expected.


Hopefully the attachment will make it through.

-Andy



On 2/8/2010 11:05 AM, Andy Colson wrote:
> Thank you, but that didnt seem to help.
>
> I left the main PROJECTION as is, and I changed my raster to this:
>
> LAYER
> NAME "raster07"
> STATUS off
> TILEINDEX "raster07.shp"
> TILEITEM "location"
> TYPE RASTER
> PROJECTION
> "+proj=utm +zone=15 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"
> END
> END
>
>
> Now, here's another problem (to add to the layers of projections). I'm
> using OpenLayers on the front end (for the website), and I have no idea
> what all this projection stuff means, and wasnt sure how to set it in
> OpenLayers.
>
> I create the map like this:
>
> map = new OpenLayers.Map("map", {
> controls: [
> new OpenLayers.Control.Navigation( {zoomWheelEnabled: false} ),
> new OpenLayers.Control.PanZoomBar(),
> ls
> ],
> numZoomLevels: 12,
> maxExtent: new OpenLayers.Bounds([%minx%], [%miny%], [%maxx%], [%maxy%]),
> maxResolution: 250,
> unit: 'm',
> projection: "EPSG:2816"}
> );
>
> It seems to work. I totally guessed at EPSG:2816, I have no idea if it
> matches:
> "+proj=tmerc +lat_0=35.83333333333334 +lon_0=-92.5 +k=0.999933333
> +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs"
>
> The shapefiles show up, as do most of the arials.
>
>  > Are these all publicly available datasets?
>
> I dont know, they are from the Livingston county Missouri Assessors
> office. Hum... did I say Minnesota before? That was wrong, its MO not
> MN. Its about 2 gig, and I could email and ask them if it really comes
> down to it.
>
> -Andy
>
> On 2/8/2010 10:36 AM, Fawcett, David (MPCA) wrote:
>> Andy,
>>
>> It looks like you have data in at least two or three different
>> projections. If you want to display them in the same map using
>> MapServer, you need:
>>
>> - a projection block at the MAP level. This defines your output
>> projection. The map UNITS and EXTENT need tobe consistent with this
>> projection.
>>
>> - a projection block for each layer. This is used to tell MapServer
>> what projection the data for that layer is currently in.
>>
>> MapServer then uses this info to warp all of the layers into the same
>> output projection.
>>
>> Does this help?
>>
>> Are these all publicly available datasets?
>>
>> David.
>>
>>
>>
>> -----Original Message----- From:
>> mapserver-users-bounces at lists.osgeo.org
>> [mailto:mapserver-users-bounces at lists.osgeo.org] On Behalf Of Andy
>> Colson Sent: Monday, February 08, 2010 10:11 AM To: mapserver-users
>> Subject: [mapserver-users] many raster projections
>>
>> Ok, I'm confused. Most of the problem is probably because I'm a
>> programmer, not a mapper.
>>
>> I have arial imagery of Minnesota, multiple sets, each seeming to be
>> in a different projection.
>>
>> At the top of my mapfile I have this:
>>
>> MAP EXTENT 1285920 1377857 1426390 1506589
>>
>> PROJECTION "+proj=tmerc +lat_0=35.83333333333334 +lon_0=-92.5
>> +k=0.999933333 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs"
>> END
>>
>> There are no other PROJECTION settings. Most of the images seem to
>> work. I got the projtext above from loading one of the shapefiles
>> into qgis and looking at its properties. I just kinda hoped they'd
>> all be the same. I did try adding a PROJECTION to the raster layer,
>> but it didnt seem to work. So I thought I'd "just" re-project the
>> imagery.
>>
>>
>> Three sets of imagery seem to work (I didnt touch them, they always
>> worked), they have (again copied from qgis) proj text's of:
>>
>> +proj=tmerc +lat_0=35.83333333333334 +lon_0=-92.5 +k=0.999933333
>> +x_0=500000.0000000001 +y_0=0 +ellps=GRS80 +datum=NAD83
>> +to_meter=0.3048006096012192 +no_defs
>>
>> +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
>>
>> and
>>
>> +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
>>
>> One set is displayed way off the map, way far away from where its
>> supposed to be. Its projtext is:
>>
>> +proj=utm +zone=15 +ellps=GRS80 +datum=NAD83 +units=m +no_defs
>>
>> I attempted to use gdal_translate and gdal_warp.
>>
>> gdal_translate made a new image, with projtext: +proj=tmerc
>> +lat_0=35.83333333333334 +lon_0=-92.5 +k=0.999933333 +x_0=500000
>> +y_0=0 +ellps=GRS80 +units=m +no_defs
>>
>> but its in the exact same spot, it didnt re-position to where images
>> are supposed to be.
>>
>> A gdalwarp moved it, but not to the right place.
>>
>> So I'm confused, and not sure what I'm doing. Hopefully there is
>> some info in here someone can use to set me on the right path. Any
>> help would be appreciated.
>>
>> -Andy _______________________________________________ mapserver-users
>> mailing list mapserver-users at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/mapserver-users
>
> _______________________________________________
> mapserver-users mailing list
> mapserver-users at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/mapserver-users

-------------- next part --------------
MAP
	NAME LivingstonMO
	SIZE 800 800
	UNITS feet
	SHAPEPATH '/pub/maps/LivingstonMO/shape'
	FONTSET "/pub/www/fonts/fonts.txt"

	EXTENT 1285920 1377857 1426390 1506589

	PROJECTION
		"+proj=tmerc +lat_0=35.83333333333334 +lon_0=-92.5 +k=0.999933333 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs"
	END

	IMAGECOLOR 192 192 192
	IMAGEQUALITY 95
	IMAGETYPE jpeg
	OUTPUTFORMAT
		NAME jpeg
		DRIVER 'GD/JPEG'
		MIMETYPE 'image/jpeg'
		#IMAGEMODE PC256
		EXTENSION 'jpeg'
	END


	LAYER
		NAME "raster09_chil"
		STATUS off
		TILEINDEX "raster09_chil.shp"
		TILEITEM "location"
		TYPE RASTER
	END

	LAYER
		NAME "raster07"
		STATUS on
		TILEINDEX "raster07.shp"
		TILEITEM "location"
		TYPE RASTER
		PROJECTION
		"+proj=utm +zone=15 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"
		END
	END
END