Bringing this back onto the mailing list....
Re: [mapserver-users] INSIDE Idaho... From: "greg.lu...@scu.edu.au" <greg.lu...@scu.edu.au>...Add to Contacts To: Jon Miner <cp_...@yahoo.com> ______________________________ Hi Jon, I'm very much NOT a mapserver expert (a beginner in fact), but one thing I noticed (it might mean nothing) is: UNIT["Degree",0.0174532925199433]] I realise that ogr is saying this, but I thought UTM projection unit is metre not degree, and I'm guessing that number (0.0174532925199433) is how many degrees in a metre at that latitude. Also the gdalinfo is listing the origin and pixel size in metres. One thing you could do is just put one layer in a mapfile and set the mapfile projection to NAD83 (or whatever the epsg code is for that) to see if you can actually see the data at all, and then project it into what you need, using "metres" (or maybe "meters") as the unit, somehow. Sorry I can't help more, but maybe this might point you in a useful direction. Thanks, Greg. ----- Forwarded Message ---- From: Jon Miner <cp_...@yahoo.com> To: glu...@scu.edu.au Sent: Mon, February 8, 2010 8:16:22 PM Subject: Re: [mapserver-users] INSIDE Idaho Hi Greg - I've looked at the parameters... Here is the DRG file: http://data.insideidaho.org/data2/24k/usgs/drg/idtm/quad/nocollar/archive/43112/cirbu3se_24k_usgs_drg.tgz So I extract that file and then run gdalinfo on the tif file: gdalinfo cirbu3se_24k_usgs_drg.tif Driver: GTiff/GeoTIFF Files: cirbu3se_24k_usgs_drg.tif Size is 4221, 5752 Coordinate System is `' Origin = (590767.311137848766521,281203.236009103886317) Pixel Size = (2.438400000000000,-2.438400000000000) Metadata: TIFFTAG_SOFTWARE=Arc/Info TIFFTAG_XRESOLUTION=0.41 TIFFTAG_YRESOLUTION=0.41 TIFFTAG_RESOLUTIONUNIT=1 (unitless) Image Structure Metadata: COMPRESSION=PACKBITS INTERLEAVE=BAND Corner Coordinates: Upper Left ( 590767.311, 281203.236) Lower Left ( 590767.311, 267177.559) Upper Right ( 601059.798, 281203.236) Lower Right ( 601059.798, 267177.559) Center ( 595913.554, 274190.398) ...etc... Note the blank "coordinate system" entry... Using values from http://spatialreference.org, I came up with the PROJECTION values in my map file, which seems to work, since mapserver renders the DRG. So far so good... Next I download some vector data: http://data.insideidaho.org/data/USCB/TIGER2000/archive/butte/roadsl_butte_uscb.tgz and extract it, then check it with ogrinfo: $ ogrinfo roadsl_butte_uscb.shp roadsl_butte_uscb -so INFO: Open of `roadsl_butte_uscb.shp' using driver `ESRI Shapefile' successful. Layer name: roadsl_butte_uscb Geometry: Line String Feature Count: 2596 Extent: (-113.641719, 43.284631) - (-112.694959, 44.233133) Layer SRS WKT: 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]] ..etc.. If I'm not mistaken, I need to translate this to match the projection of the DRG file... I've been doing this with ogr2ogr: ogr2ogr -t_srs idtm83.wkt -f "ESRI Shapefile" temp roadsl_butte_uscb.shp where the idtm83.wkt contains: PROJCS["NAD_1983_Transverse_Mercator", GEOGCS["GCS_North_American_1983", DATUM["D_North_American_1983", SPHEROID["GRS_1980",6378137.0,298.257222101]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["False_Easting",2500000.0], PARAMETER["False_Northing",1200000.0], PARAMETER["Central_Meridian",-114.0], PARAMETER["Scale_Factor",0.9996], PARAMETER["Latitude_Of_Origin",42.0], UNIT["Meter",1.0]] ogrinfo on the new file reports: Layer name: roadsl_butte_uscb Geometry: Line String Feature Count: 2596 Extent: (2528908.642544, 1342938.885743) - (2605077.684651, 1448186.034775) Layer SRS WKT: PROJCS["NAD_1983_Transverse_Mercator", GEOGCS["GCS_North_American_1983", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980",6378137.0,298.257222101]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["False_Easting",2500000.0], PARAMETER["False_Northing",1200000.0], PARAMETER["Central_Meridian",-114.0], PARAMETER["Scale_Factor",0.9996], PARAMETER["Latitude_Of_Origin",42.0], UNIT["Meter",1.0]] To me (as a programmer, not a map guru) everything seems to match up... Can you spot something that is not right? Thanks - Jon ________________________________ From: "greg.lu...@scu.edu.au" <greg.lu...@scu.edu.au> To: Jon Miner <cp_...@yahoo.com> Sent: Mon, February 8, 2010 5:50:47 PM Subject: Re: [mapserver-users] INSIDE Idaho Hi John, According to this page: http://inside.uidaho.edu/tutorial/data/data.htm all INSIDE data: "are available in Idaho Transverse Mercator (IDTM) projection. Additionally, some data are available in Universal Transverse Mercator (UTM) projection." (UTM zone 11 or zone 12). Parameters are here: http://inside.uidaho.edu/geodata/projparam.htm Is this any help? Thanks, Greg. Greg Luker GIS Lab Manager Southern Cross University 61 2 66203026 _______________________________________________ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
