[Mapserver-users] geotiff output problem

[Aaron D. Hunt]

I worked when I used shp2img. Is there any reason it is not working in the
script?


----- Original Message ----- 
From: "Frank Warmerdam" <warmerdam at pobox.com>
To: "Aaron D. Hunt" <hunt at zedxinc.com>
Cc: <mapserver-users at lists.gis.umn.edu>
Sent: Tuesday, January 27, 2004 1:30 PM
Subject: Re: [Mapserver-users] geotiff output problem


> 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
>