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

[Jon Miner]

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