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
