[Mapserver-users] MS4.0 and Raster (GeoTiff) assistance needed

humberto at pastoraldacrianca.org.br humberto at pastoraldacrianca.org.br
Thu Nov 20 09:38:04 PST 2003


Phil,

if you have gdal instaled, try to run gdalinfo:

gdalinfo name.tif

here is a example from Markus Neteler

A good idea is always to use 'gdalinfo' to get information about Projection,
coordinate system, datum etc.:

gdalinfo shdrlfi020l.tif
Driver: GTiff/GeoTIFF
Size is 10366, 7273
Coordinate System is:
PROJCS["unnamed",
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["Lambert_Azimuthal_Equal_Area"],
PARAMETER["latitude_of_center",45],
PARAMETER["longitude_of_center",-100],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-6086629.000000,4488761.000000)
Pixel Size = (1000.000000,-1000.000000)
Metadata:
TIFFTAG_IMAGEDESCRIPTION=Color shaded relief of North America
TIFFTAG_SOFTWARE=USGS DRG production software
TIFFTAG_DATETIME=2003:02:06 13:19:34
TIFFTAG_XRESOLUTION=254
TIFFTAG_YRESOLUTION=254
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Corner Coordinates:
Upper Left (-6086629.000, 4488761.000) (156d 3'31.10"E, 37d45'55.65"N)
Lower Left (-6086629.000,-2784239.000) (154d27'44.36"W, 3d 2'12.16"N)
Upper Right ( 4279371.000, 4488761.000) ( 4d 0'23.08"W, 53d59'22.66"N)
Lower Right ( 4279371.000,-2784239.000) ( 61d 9'32.19"W, 11d18'12.56"N)
Center ( -903629.000, 852261.000) (113d12'19.80"W, 51d59'8.84"N)
Band 1 Block=10366x1 Type=Byte, ColorInterp=Red
Band 2 Block=10366x1 Type=Byte, ColorInterp=Green
Band 3 Block=10366x1 Type=Byte, ColorInterp=Blue

-> This is "Lambert_Azimuthal_Equal_Area".

Cheers

Markus

Humberto Cereser Ibanez
Pastoral da Criança - Brazil

Citando Frank Warmerdam <warmerdam at pobox.com>:

> Phil Shevlin wrote:
> > I have a map with a number of layers working.  I am trying to include a
> > raster layer (.tif).  When I include the layer:
> > 
> >   LAYER
> >     NAME "Test"
> >     STATUS ON
> >     DATA "E:/Inetpub/wwwroot/mapserver/raster/4424.tif"
> >     TYPE RASTER
> >     TEMPLATE "template_hi.html"
> >   END
> > 
> > The map shows but the layer does not appear.  I am assuming that it is a
> > projection problem.  My extent is in degrees:
> > 
> > EXTENT -81.917574 26.569844 -81.755026 26.691756
> > 
> > But the .tfw file appears otherwise:
> > 
> > 4.00000000000000
> > 0.00000000000000
> > 0.00000000000000
> > -4.00000000000000
> > 669000.00000000000000
> > 855000.00000000000000
> > 
> > How do I go about including this layer?  Or, can someone point me to some
> > tutorials (I have found one that helps me).
> > 
> 
> Phil,
> 
> The raster is clearly in a different projection.  You will need to establish
> what coordinate system it is, and properly describe it with a PROJECTION
> block
> in the layer.  You will also need an appropriate lat/long projectio block
> for
> the map level if your generated map is currently in lat/long.
> 
> To know the coordinate system you will have to look at your dataset and
> talk to whoever you got it from.
> 
> Good luck,
> 
> -- 
> ---------------------------------------+--------------------------------------
> 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
> 
> 
> _______________________________________________
> Mapserver-users mailing list
> Mapserver-users at lists.gis.umn.edu
> http://lists.gis.umn.edu/mailman/listinfo/mapserver-users
> 




-------------------------------------------------
E-Mail enviado pelo Webmail Pastoral da Criança
usando IMP: http://horde.org/imp/



More information about the MapServer-users mailing list