Thanks, that solved my problem - the map projects perfectly now. I feel foolish. ________________________________________ From: Smith, Michael ERDC-CRREL-NH [[email protected]] Sent: Saturday, August 20, 2011 2:25 AM To: Ed Resnick; [email protected] Subject: <Ext> Re: [mapserver-users] Problem displaying a NAD83 tile in WGS84 (or via WMS)
Ed, Your mrsid image is in UTM Zone 10, NAD 1983. You should set PROJECTION "init=epsg:26910" END For your layer (http://spatialreference.org/ref/epsg/26910/) -- Mike Michael Smith Remote Sensing/GIS Center US Army Corps of Engineers On 8/19/11 5:23 PM, "Ed Resnick" <[email protected]> wrote: >Here's a WMS request, broken up for readability: > >------------------------------------------------------- >http://10.10.10.206/cgi-bin/mapserv.exe?map=E:/osm4326/mrsid.map > &service=WMS > &version=1.1.1 > &request=GetMap > &layers=default > &format=image/png > &width=256 > &height=256 > &srs=EPSG:4326 > &bbox=-122.875,45.25,-122.75,45.375 > &transparent=TRUE >------------------------------------------------------- > >Here's the gdalinfo dump: > >------------------------------------------------------- >Driver: MrSID/Multi-resolution Seamless Image Database (MrSID) >Files: wa.sid.sid > wa.sid.sdw > wa.sid.sid.aux.xml >Size is 1742, 1458 >Coordinate System is: >PROJCS["NAD_1983_UTM_Zone_10N", > GEOGCS["GCS_North_American_1983", > DATUM["North_American_Datum_1983", > SPHEROID["GRS_1980",6378137.0,298.257222101]], > PRIMEM["Greenwich",0.0], > UNIT["Degree",0.0174532925199433]], > PROJECTION["Transverse_Mercator"], > PARAMETER["false_easting",500000.0], > PARAMETER["false_northing",0.0], > PARAMETER["central_meridian",-123.0], > PARAMETER["scale_factor",0.9996], > PARAMETER["latitude_of_origin",0.0], > UNIT["Meter",1.0]] >Origin = (555230.000000000000000,5048039.000000000000000) >Pixel Size = (1.000009239300614,-1.000000000000000) >Metadata: > IMAGE__INPUT_NAME=C:\DOCUME~1\PAULW~2.ATT\LOCALS~1\Temp\img1AE.tmp > IMAGE__INPUT_FILE_SIZE=7619508.000000 > GEO__spref_horizsys_cordsysn_geogcsn=GCS_North_American_1983 > GEO__spref_horizsys_cordsysn_projcsn=NAD_1983_UTM_Zone_10N > GEO__spref_horizsys_planar_gridsys_gridsysn=Universal Transverse >Mercator > GEO__spref_horizsys_planar_gridsys_utm_utmzone=10 > GEO__spref_horizsys_planar_gridsys_utm_transmer_sfctrmer=0.999600 > GEO__spref_horizsys_planar_gridsys_utm_transmer_longcm=-123.000000 > GEO__spref_horizsys_planar_gridsys_utm_transmer_latprjo=0.000000 > GEO__spref_horizsys_planar_gridsys_utm_transmer_feast=500000.000000 > GEO__spref_horizsys_planar_gridsys_utm_transmer_fnorth=0.000000 > GEO__spref_horizsys_planar_planci_plance=row and column > GEO__spref_horizsys_planar_planci_plandu=meters > GEO__spref_horizsys_planar_planci_coordrep_absres=1.000009 > GEO__spref_horizsys_planar_planci_coordrep_ordres=1.000000 > GEO__spref_horizsys_geodetic_horizdn=North American Datum of 1983 > GEO__spref_horizsys_geodetic_ellips=Geodetic Reference System 80 > GEO__spref_horizsys_geodetic_semiaxis=6378137.000000 > GEO__spref_horizsys_geodetic_denflat=298.257222 > IMAGE__COMPRESSION_VERSION=1,6,3 > IMAGE__TARGET_COMPRESSION_RATIO=20.000000 > IMAGE__COMPRESSION_NLEV=6 > IMAGE__COMPRESSION_WEIGHT=4.000000 > IMAGE__COMPRESSION_GAMMA=2.000000 > IMAGE__COMPRESSION_BLOCK_SIZE=512 > IMAGE__CREATION_DATE=Wed Dec 02 14:18:16 2009 > > IMAGE__NO_DITHER=0 > VERSION=MG2 > IMAGE__WIDTH=1742 > IMAGE__HEIGHT=1458 > IMAGE__NO_DATA_VALUE=15,172,40 > IMAGE__TRANSPARENT_DATA_VALUE=0,0,176 > IMAGE__COLOR_SCHEME=0 > IMAGE__DATA_TYPE=0 > IMAGE__XY_ORIGIN=555230.500005,5048038.500000 > IMAGE__X_RESOLUTION=1.000009 > IMAGE__Y_RESOLUTION=1.000000 >Image Structure Metadata: > INTERLEAVE=PIXEL >Corner Coordinates: >Upper Left ( 555230.000, 5048039.000) (122d17'31.28"W, 45d35' 1.32"N) >Lower Left ( 555230.000, 5046581.000) (122d17'31.87"W, 45d34'14.07"N) >Upper Right ( 556972.016, 5048039.000) (122d16'10.90"W, 45d35' 0.81"N) >Lower Right ( 556972.016, 5046581.000) (122d16'11.51"W, 45d34'13.57"N) >Center ( 556101.008, 5047310.000) (122d16'51.39"W, 45d34'37.44"N) >Band 1 Block=1742x128 Type=Byte, ColorInterp=Red > Minimum=60.000, Maximum=222.000, Mean=119.606, StdDev=36.132 > Overviews: 871x729, 436x365, 218x183, 109x92, 55x46, 28x23 >Band 2 Block=1742x128 Type=Byte, ColorInterp=Green > Minimum=84.000, Maximum=220.000, Mean=133.956, StdDev=21.319 > Overviews: 871x729, 436x365, 218x183, 109x92, 55x46, 28x23 >Band 3 Block=1742x128 Type=Byte, ColorInterp=Blue > Minimum=95.000, Maximum=205.000, Mean=119.015, StdDev=12.620 > Overviews: 871x729, 436x365, 218x183, 109x92, 55x46, 28x23 >------------------------------------------------------- > >And finally here is the map file: > >------------------------------------------------------- >MAP >OUTPUTFORMAT > NAME agg > DRIVER AGG/PNG > IMAGEMODE RGB > FORMATOPTION "INTERLACE=false" > MIMETYPE "image/png" >END >NAME "MrSID" >IMAGETYPE agg >MAXSIZE 4000 >SIZE 800 800 >IMAGECOLOR "#B3C6D4" >WEB > METADATA > ows_enable_request "*" > labelcache_map_edge_buffer "-10" > wms_srs "EPSG:4326 EPSG:4269" > wms_title "MrSID" > wms_onlineresource >"http://10.10.10.206/cgi-bin/mapserv.exe?map=E:/osm4326/mrsid.map"; > END >END >CONFIG "PROJ_LIB" "C:\ms4w\proj\nad" >CONFIG "MS_ERRORFILE" "E:\errors.txt" >DEBUG 5 >SHAPEPATH "E:/osm4326/MrSID/" >EXTENT -123.0 45.0 -122.0 46.0 >UNITS dd >PROJECTION >"init=epsg:4326" >END >LAYER > METADATA > wms_title "A" > wms_group_title "default" > wms_srs "EPSG:4269" > wms_extent "555230 5046581 556972 5048039" > END > TYPE RASTER > STATUS ON > GROUP "default" > NAME "A" > PROJECTION > "init=epsg:4269" > END > EXTENT 555230 5046581 556972 5048039 > DATA "wa.sid.sid" >END >END >------------------------------------------------------- > > >________________________________________ >From: Rahkonen Jukka [[email protected]] >Sent: Thursday, August 18, 2011 9:51 PM >To: Ed Resnick; '[email protected]' >Subject: <Ext> Re: [mapserver-users] Problem displaying a NAD83 tile in >WGS84 (or via WMS) > >Hi, > >Could you give a bit more information? Your use case is the most common >that exists and all instructions should be at >http://mapserver.org/ogc/wms_server.html >It would be good to see the gdalinfo report of your image and your >mapfile and the exact WMS request you were using including the WMS >version. At least these parts of mapfile are important > >MAP level: extent, projection and metadata (wms_srs and >wms_inable_request at least) >LAYER level projection. > >I can see some really odd LatLonBoundingBoxes in my own MS 6.0 WMS >service too. I will investigate them a bit further. If it is just that >the Mapserver is creating the GetCapabilities document wrong then the WMS >service should still be usable. > >-Jukka Rahkonen- > > > > >> -----Alkuperäinen viesti----- >> Lähettäjä: [email protected] >> [mailto:[email protected]] Puolesta Ed Resnick >> Lähetetty: 19. elokuuta 2011 2:47 >> Vastaanottaja: [email protected] >> Aihe: [mapserver-users] Problem displaying a NAD83 tile in >> WGS84 (or via WMS) >> >> I have a MrSID file in NAD83 projection (EPGS 4269) which I >> want to display via WMS in WGS84 projection (EPGS 4326). I >> cannot get this to work. >> >> I have tried setting the entire map-file's projection to >> 4269, with the extent specified in that projection (in >> meters), and shp2img shows the image, but attempting to view >> it via WMS (requesting the 4326 projection) does not work, >> and the GetCapabilities response file is completely wrong - >> the LatLonBoundingBox tags actually list the 4269 bounding >> box in meters. Excerpt: >> >> <LatLonBoundingBox minx="555230" miny="5.04658e+006" >> maxx="556972" maxy="5.04804e+006" /> >> <BoundingBox SRS="EPSG:4269" >> minx="555230" miny="5.04658e+006" >> maxx="556972" maxy="5.04804e+006" /> >> >> The BoundingBox is correct; the LatLonBoundingBox obviously is not. >> >> I have also tried specifying 4326 as the map projection and >> 4269 as the layer projection, with the map's extent in >> lat/lon and the layer's extent as the above value (meters). >> Not only did this not work, but even shp2img just shows a >> blank image. If I change the map's extent to the value in >> meters (even though its projection is specified as 4326), >> then shp2img works, but the WMS feed still does not; the >> GetCapabilities still is wrong. >> >> Is MapServer simply buggy or is there a >> solution?_______________________________________________ >> mapserver-users mailing list >> [email protected] >> http://lists.osgeo.org/mailman/listinfo/mapserver-users >> >_______________________________________________ >mapserver-users mailing list >[email protected] >http://lists.osgeo.org/mailman/listinfo/mapserver-users _______________________________________________ mapserver-users mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/mapserver-users
