[mapserver-users] New to mapserver; question about USGS DRG/GeoTIFF files

Sun Jan 31 21:28:09 PST 2010

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
>