[Mapserver-users] geotiff output problem

[Frank Warmerdam]

Aaron D. Hunt wrote:
> I am trying to get mapserver to output geotiffs. When I do so the tiff 
> is made fine but the geo in-coding is blank (see listgeo out put below). 
> I have mapserver 4.0.1 with gdal-cvs-031112, tiff-v3.6.0 and 
> libgeotiff-1.2.1 on redhat 9.0 (configure call below). The perl script 
> and Map file I am using is also bellow. Does anyone have any ideas.

Aaron,

What you did seems fine.  It is possible that the "init=epsg:26915" isn't
being expanded properly though even then the coordinates should have been
written.  I tried running shp2img on the following map file and produced a
GeoTIFF file with coordinates and a coordinate system.

MAP
NAME  "w2qpgi5qsg4k"
DEBUG on
SIZE  400 400
EXTENT  -6500000 -6500000 6500000 6500000

OUTPUTFORMAT
   NAME GT
   DRIVER "GDAL/GTiff"
   IMAGEMODE RGB
END

LAYER
NAME  "world"
TYPE  Raster
STATUS  default
DATA "/home/warmerda/public_html/world.tif"
PROCESSING "DITHER=YES"

PROJECTION
"proj=latlong"
END
END


#PROJECTION
#"proj=ortho"
#"lat_0=0.0"
#"lon_0=0.0"
#"x_0=0.0"
#"y_0=0.0"
#END
PROJECTION
   "init=epsg:26915"
END

END
warmerda at gdal2200[242]% gdalinfo out.tif
GDAL: GDALOpen(out.tif) succeeds as GTiff.

Driver: GTiff/GeoTIFF
Size is 400, 400
Coordinate System is:
PROJCS["NAD83 / UTM zone 15N",
     GEOGCS["NAD83",
         DATUM["North_American_Datum_1983",
             SPHEROID["GRS 1980",6378137,298.2572221010042,
                 AUTHORITY["EPSG","7019"]],
             AUTHORITY["EPSG","6269"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4269"]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",-93],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","26915"]]
Origin = (-6500000.000000,6500000.000000)
Pixel Size = (32500.00000000,-32500.00000000)
Corner Coordinates:
Upper Left  (-6500000.000, 6500000.000) ( 95d47'40.98"W,280d 4'9.61"N)
Lower Left  (-6500000.000,-6500000.000) ( 95d47'40.98"W,280d 4'9.61"S)
Upper Right ( 6500000.000, 6500000.000) ( 86d 6'45.06"E,101d28'3.37"N)
Lower Right ( 6500000.000,-6500000.000) ( 86d 6'45.06"E,101d28'3.37"S)
Center      (   0.0000000,   0.0000000) ( 97d29'19.48"W,  0d 0'0.01"N)
Band 1 Block=400x20 Type=Byte, ColorInterp=Red
Band 2 Block=400x20 Type=Byte, ColorInterp=Green
Band 3 Block=400x20 Type=Byte, ColorInterp=Blue
GDAL: GDALClose(out.tif)

GDAL: GDALDeregister_GTiff() called.

I really don't know what is going wrong at your end.   You could try
updating to more recent Mapserver versions, but I don't really see that
this should be needed.  Could you try running shp2img on your .map and
verify that this also exibits the same problem?

Best regards,

-- 
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent