[mapserver-users] Problems getting a JP2 or TIFF raster layer drawen

Frank Warmerdam warmerdam at pobox.com
Wed May 27 12:56:50 EDT 2009


Michael Da Silva Pereira wrote:
> It displays everything normally, except the raster, which is just not drawn
> 
> GDAL: Auto register /usr/lib/gdal15plugins/gdal_GRASS.so using 
> GDALRegister_GRASS.
> [Wed May 27 18:04:41 2009].197691 msLoadMap(): 0.570s
> [Wed May 27 18:04:41 2009].414187 msDrawMap(): Layer 0 
> (vereenigingroads), 0.194s
> [Wed May 27 18:04:41 2009].414307 msDrawRasterLayerLow(2627D): entering.
> JPEG2000: IHDR box found. Dump: width=20018, height=22269, numcmpts=3, bpp=8
> JPEG2000: Component 0: bpp=8, signedness=0
> JPEG2000: Component 1: bpp=8, signedness=0
> JPEG2000: Component 2: bpp=8, signedness=0
> GDALJP2Metadata: Got projection from GeoJP2 (geotiff) box: PROJCS["WGS 
> 84 / UTM zone 35S",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 
> 84",6378137,298.2572235629972,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",27],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","32735"]]
> GDAL: GDALOpen(/var/www/ms-tmp/./2627D.JP2) succeeds as JPEG2000.
> GDAL: GDALClose(/var/www/ms-tmp/./2627D.JP2)

Michael,

The file is opening fine, so I would suggest very carefully considering
the georeferenced bounds reported by gdalinfo for the file, and the
map region you are drawing.

In fact, the coordinate system of file this is UTM 35S, but your map is
in lat/long degrees.  If you want MapServer to reproject a UTM image
to lat/long on the fly you will need to specify appropriate PROJECTION
blocks for the map and the raster layer - which I observe you have not done.

For the MAP perhaps something like:

PROJECTION
   "+proj=latlong +datum=WGS84"
END

and for the raster layer something like:

PROJECTION
   "+proj=utm +zone=35 +south +datum=WGS84"
END

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



More information about the mapserver-users mailing list