[mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files
Jon Miner
cp_n18 at yahoo.com
Mon Feb 1 04:45:47 PST 2010
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 <Jukka.Rahkonen at mmmtike.fi>
To: Jon Miner <cp_n18 at yahoo.com>; mapserver-users at 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é <randre at gmail.com>
To: Jon Miner <cp_n18 at yahoo.com>
Cc: mapserver-users at 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 <cp_n18 at yahoo.com> 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
>mapserver-users at lists.osgeo.org
>http://lists.osgeo.org/mailman/listinfo/mapserver-users
>
More information about the MapServer-users
mailing list