Re: [mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files
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ä: mapserver-users-boun...@lists.osgeo.org > [mailto:mapserver-users-boun...@lists.osgeo.org] Puolesta Jon Miner > Lähetetty: 1. helmikuuta 2010 14:46 > Vastaanottaja: mapserver-users@lists.osgeo.org > 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 > To: Jon Miner ; mapserver-users@lists.osgeo.org > 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é > To: Jon Miner > Cc: mapserver-users@lists.osgeo.org > 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 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 no
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 To: Jon Miner ; mapserver-users@lists.osgeo.org 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é To: Jon Miner Cc: mapserver-users@lists.osgeo.org 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 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",50], >>PARAMETER["false_northing",0], >>UNIT["metre",1, >>AUTHORITY["EPSG","9001"]], >>AUTHORITY["EPSG","26719"]] >>Origin = (295125.99949388772,4778305.84265616444) >>Pixel Size = (2.4384000,-2.4384000) >>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"
Re: [mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files
Hi, Another thing, your mapfile has: MINSCALEDENOM 24000 MAXSCALEDENOM 24000 I means that the image is shown only at a scale 1:24000, exactly. Comment out these lines with #. The demo application is obviously asking maps always in this scale but it is difficult for you do the same with direct requests. And learn to send request with a browser by using CGI calls or WMS, then you will have full control on all the parameters. -Jukka- ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
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é To: Jon Miner Cc: mapserver-users@lists.osgeo.org 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 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",50], >>PARAMETER["false_northing",0], >>UNIT["metre",1, >>AUTHORITY["EPSG","9001"]], >>AUTHORITY["EPSG","26719"]] >>Origin = (295125.99949388772,4778305.84265616444) >>Pixel Size = (2.4384000,-2.4384000) >>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, ColorInt
Re: [mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files
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é To: Jon Miner Cc: mapserver-users@lists.osgeo.org 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 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",50], >>PARAMETER["false_northing",0], >>UNIT["metre",1, >>AUTHORITY["EPSG","9001"]], >>AUTHORITY["EPSG","26719"]] >>Origin = (295125.99949388772,4778305.84265616444) >>Pixel Size = (2.4384000,-2.4384000) >>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 >mapserver-users@lists.osgeo.org >http://lists.osgeo.org/mailman/listinfo/mapserver-users > ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
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 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",50], >PARAMETER["false_northing",0], >UNIT["metre",1, >AUTHORITY["EPSG","9001"]], >AUTHORITY["EPSG","26719"]] > Origin = (295125.99949388772,4778305.84265616444) > Pixel Size = (2.4384000,-2.4384000) > 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 > mapserver-users@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/mapserver-users > ___ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users
[mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files
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",50], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AUTHORITY["EPSG","26719"]] Origin = (295125.99949388772,4778305.84265616444) Pixel Size = (2.4384000,-2.4384000) 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 mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users