Hi, EPSG:4326 bounds do work all the time but in order to see your image the epsg:4326 extents must be the same than the epsg:26719 extents of the image projected to epsg:4326.
This may go wrong, I have not used cs2cs program much, but let's hope I did this correctly So I tried to take the 26716 coordinates and reproject them to 4321 d:\FWTools2.4.6>cs2cs +init=epsg:26716 +to +init=epsg:4326 #reproject from 26716 -> to 4326 295125.999 4761593.049 #input coordinates 89d30'45.961"W 42d58'52.872"N 0.000 #cs2cs output 308264.099 4778305.843 #input coordinates 89d21'27.097"W 43d8'6.525"N 0.000 #cs2cs output If this want OK it means that your image should be inside an epsg:4326 box defined by EXTENT -90 42 -89 44 If it went wrong some more experienced cs2cs user will certainly help you further. -Jukka- > -----Alkuperäinen viesti----- > Lähettäjä: [email protected] > [mailto:[email protected]] Puolesta Jon Miner > Lähetetty: 1. helmikuuta 2010 14:46 > Vastaanottaja: [email protected] > Aihe: Re: [mapserver-users] New to mapserver;question about > USGS DRG/GeoTIFF files > > > > Thanks Jukka, > > That did it. I guess I need to better understand the > relationship between EXTENT and UNITS. I thought that my > original bounds, in degrees, would have worked? > > EXTENT -72 42 -71 44 > UNITS DD > > - Jon > > > > > > ----- Original Message ---- > From: Rahkonen Jukka <[email protected]> > To: Jon Miner <[email protected]>; [email protected] > Sent: Mon, February 1, 2010 12:28:09 AM > Subject: Re: [mapserver-users] New to mapserver;question > about USGS DRG/GeoTIFF files > > Hi, > > One possibility is that your image is not inside the spatial > window of your request. Extents or the image in epsg:26719 > system seem to be > Lower Left ( 295125.999, 4761593.049) ( 71d30'45.52"W, > 42d58'52.88"N) > Upper Right ( 308264.099, 4778305.843) ( 71d21'26.67"W, 43d 8'6.54"N) > > Extents of the MAP are in epsg:4326 system EXTENT -72 42 -71 > 44. If the image is not inside this box you won't see it > which is natural. > > Make a test by changing the MAP estents and projection > temporary to suit the image > EXTENT 295125.999 4761593.049 308264.099 4778305.843 > PROJECTION > "init=epsg:26719" > END > > If you get the image now with shp2img you know that > everything is OK with the geotiff. > > -Jukka Rahkonen- > > Jon Miner wrote: > > Hi Roger - > > Thanks for the advice, but no change. I did stumble across > the advice to try to run "shp2img" and even with debug set to > "5", get very little debugging information. > > # shp2img -m test.map -all_debug 5 -o test.png > GDAL: > GDALOpen(/opt/fgs/apps/gmap-demo-cvs_MS_VERSION_54/htdocs/./.. > /data/o43071a4.tif, this=0x9bd61b0) succeeds as GTiff. > GDAL: GDALDeregister_GTiff() called. > > The resulting test.png ends up being all white. > > The actual TIFF file (10MB) can be found here: > > http://trainguy.dyn.dhs.org/~jminer/o43071a4.tif > > Thanks > > > > > ________________________________ > From: Roger André <[email protected]> > To: Jon Miner <[email protected]> > Cc: [email protected] > Sent: Fri, January 29, 2010 4:16:08 PM > Subject: Re: [mapserver-users] New to mapserver; question > about USGS DRG/GeoTIFF files > > Hi Jon, > > You don't appear to have any projection info in your map > file. Since you have "UNITS dd" in it, I assume you want to > display the map in lat/lon WGS84? If so, your MAP file needs > an entry that looks like this: > > PROJECTION > "init=epsg:4326" > END > > and your usgs LAYER (the name should be quoted btw) should > have the UTM projection defined in it, like this: > > LAYER > NAME "usgs" > PROJECTION "init=epsg:26719" END > > Best of luck, > > Roger > -- > > > > > > On Fri, Jan 29, 2010 at 4:38 AM, Jon Miner <[email protected]> wrote: > > > > > > > >>Hi - > > > >>I'm new to Mapserver, and so far impressed by it's > capabilities. I've installed the FGS package via > fgs-mapserver_basic_5.6.0-fgs_9.5-linux-i386.bin and played > with their demo program. > > > >>I'm having a problem with a USGS DRG file that is encoded > as a GeoTIFF file. It doesn't display on the webpage, and > I'm not sure where to start debugging this. I've checked the > file with gdalinfo: > > > >>Driver: GTiff/GeoTIFF > >>Files: o43071a4.tif > >>Size is 5388, 6854 > >>Coordinate System is: > >>PROJCS["NAD27 / UTM zone 19N", > >> GEOGCS["NAD27", > >> DATUM["North_American_Datum_1927", > >> SPHEROID["Clarke 1866",6378206.4,294.9786982138982, > >> AUTHORITY["EPSG","7008"]], > >> AUTHORITY["EPSG","6267"]], > >> PRIMEM["Greenwich",0], > >> UNIT["degree",0.0174532925199433], > >> AUTHORITY["EPSG","4267"]], > >> PROJECTION["Transverse_Mercator"], > >> PARAMETER["latitude_of_origin",0], > >> PARAMETER["central_meridian",-69], > >> PARAMETER["scale_factor",0.9996], > >> PARAMETER["false_easting",500000], > >> PARAMETER["false_northing",0], > >> UNIT["metre",1, > >> AUTHORITY["EPSG","9001"]], > >> AUTHORITY["EPSG","26719"]] > >>Origin = (295125.999493999988772,4778305.842659999616444) > >>Pixel Size = (2.438400000000000,-2.438400000000000) > >>Metadata: > >> AREA_OR_POINT=Area > >> TIFFTAG_IMAGEDESCRIPTION=USGS GeoTIFF DRG 1:24000 Quad of > Manchester North. Product:522644 > >> TIFFTAG_SOFTWARE=USGS CD Archiver program tif2usgsdrg v.1.0 > >> TIFFTAG_DATETIME=1997:01:10 07:16:28 > >> TIFFTAG_XRESOLUTION=72 > >> TIFFTAG_YRESOLUTION=72 > >> TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) > >>Image Structure Metadata: > >> COMPRESSION=PACKBITS > >> INTERLEAVE=BAND > >>Corner Coordinates: > >>Upper Left ( 295125.999, 4778305.843) ( 71d31'7.66"W, 43d > 7'54.16"N) > >>Lower Left ( 295125.999, 4761593.049) ( 71d30'45.52"W, > 42d58'52.88"N) > >>Upper Right ( 308264.099, 4778305.843) ( 71d21'26.67"W, > 43d 8'6.54"N) > >>Lower Right ( 308264.099, 4761593.049) ( 71d21'5.95"W, > 42d59'5.20"N) > >>Center ( 301695.049, 4769949.446) ( 71d26'6.44"W, 43d > 3'29.80"N) > >>Band 1 Block=5388x1 Type=Byte, ColorInterp=Palette > >> Color Table (RGB with 256 entries) > >> 0: 0,0,0,255 > >> 1: 255,255,255,255 > >> 2: 0,151,164,255 > >> 3: 203,0,23,255 > >> 4: 131,66,37,255 > >> 5: 201,234,157,255 > >> 6: 137,51,128,255 > >> 7: 255,234,0,255 > >> 8: 167,226,226,255 > >> 9: 255,184,184,255 > >> 10: 218,179,214,255 > >> 11: 209,209,209,255 > >> 12: 207,164,142,255 > >>(etc...) > > > >>And my map file looks like this: > > > >>MAP > >>NAME DEMO > >>STATUS ON > >>SIZE 400 300 > >>SYMBOLSET ../etc/symbols.sym > >>EXTENT -72 42 -71 44 > >>UNITS dd > >>SCALEDENOM 24000 > >>SHAPEPATH "../data" > >>IMAGECOLOR 255 255 255 > >>FONTSET ../etc/fonts.txt > >>CONFIG "MS_ERRORFILE" "/tmp/ms_error.txt" > >># > >># Start of web interface definition > >># > >>WEB > >> MINSCALEDENOM 24000 > >> MAXSCALEDENOM 24000 > >># > >># On Windows systems, /tmp and /opt/fgs/tmp/ms_tmp/ should > be created at the root > >># of the drive where the .MAP file resides. > >># > >> IMAGEPATH "/opt/fgs/tmp/ms_tmp/" > >> IMAGEURL "/ms_tmp/" > >> LOG "/tmp/gmap.log" > >>END > > > >>LAYER > >> NAME usgs > >> METADATA > >> "DESCRIPTION" "USGS GeoTIFF" > >> END > >> TYPE RASTER > >> STATUS ON > >> DATA o43071a4.tif > >>END > > > >>END # Map File > > > > > > > >>_______________________________________________ > >>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
