Trying to display a geotiff DEM within MapServer
Brent Fraser
bfraser at GEOANALYTIC.COM
Thu Apr 26 09:50:58 PDT 2007
John,
The three most common problems when trying to display raster
using a TILEINDEX:
Coordinate System Problems:
1. EXTENT given in MAP not suitable for data in LAYER
- change EXTENT to match the bounds of the data
Data Access:
2. Pathing from TILEINDEX to raster files wrong
- put the full paths to the data in the LOCATION
attribute of the TILEINDEX
Data Organization:
3. Data representation not suitable
- use the LAYER's PROCESSING object to scale the
32bit float height values to greyscale.
- it could be that Mapserver using GDAL for input
will not read the data from a 32bit tiff.
See http://mapserver.gis.umn.edu/docs/howto/raster_data for
more details.
Brent Fraser
GeoAnalytic Inc.
Calgary, Alberta
----- Original Message -----
From: "John Mitchell" <mitchelljj98 at GMAIL.COM>
To: <MAPSERVER-USERS at LISTS.UMN.EDU>
Sent: Thursday, April 26, 2007 10:21 AM
Subject: [UMN_MAPSERVER-USERS] Trying to display a geotiff
DEM within MapServer
Hi,
I am trying to display a DEM which is in geotiff format
(gdalinfo
listing below) within MapServer (mapfile layer
listing below) and the image is all white.
When trying to display this same image within the desktop
application
called "global mapper" it first asks
whether you want to display this data as an image or as a
DEM. If
select as an image then is displays an image
that is all the same color similar to mapserver and if
select as a DEM
then it displays a color coded elevated image.
Is their a way to display a DEM within mapserver that
displays the
elevated image instead of the all white image
that I am currently displaying?
Thanks,
John J. Mitchell
------------------------------------------------------------
-----------
MAP
# all temp files get prefixed with this string
NAME All_Raster_Layers_ignore
# always returns a map
STATUS ON
# image format options
OUTPUTFORMAT
NAME png
DRIVER "GD/PNG"
MIMETYPE "image/png"
IMAGEMODE RGB
EXTENSION "png"
END
# projection of the INTERFACE that this mapfile supports
PROJECTION
"init=epsg:4326"
END
# minx miny maxx maxy
EXTENT -180 -90 180 90 # World
# units of map/app
UNITS DD
WEB
# this is set in the form which initializes the page, but
can be set
here if you want
IMAGEPATH "/ms4w/tmp/ms_tmp/"
# this is set in the form which initializes the page, but
can be set
here if you want
IMAGEURL "/ms_tmp/"
END
# start of layer definitions
LAYER
# name of layer
NAME Columbia_1000000012_utm_19_fullResolution_ignore
GROUP "Columbia_1000000012_utm_19"
PROJECTION
"init=epsg:32619"
END
TYPE RASTER
STATUS ON
TILEINDEX
"Z:\GlobalView\NGA\Columbia\1000000012\utm_19\fullResolution
.shp"
TILEITEM "Location"
# OFFSITE 0 0 0
MAXSCALE 10000
END
------------------------------------------------------------
-----------
Driver: GTiff/GeoTIFF
Size is 2792, 2794
Coordinate System is:
PROJCS["WGS 84 / UTM zone 19N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235630016,
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",-69],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32619"]]
Origin = (168297.500000000000000,761002.500000000000000)
Pixel Size = (5.000000000000000,-5.000000000000000)
Metadata:
AREA_OR_POINT=Area
Corner Coordinates:
Upper Left ( 168297.500, 761002.500) ( 72d 0'3.71"W,
6d52'31.02"N)
Lower Left ( 168297.500, 747032.500) ( 72d 0'0.88"W,
6d44'56.70"N)
Upper Right ( 182257.500, 761002.500) ( 71d52'29.42"W,
6d52'33.81"N)
Lower Right ( 182257.500, 747032.500) ( 71d52'26.72"W,
6d44'59.44"N)
Center ( 175277.500, 754017.500) ( 71d56'15.18"W,
6d48'45.26"N)
Band 1 Block=256x256 Type=Float32, ColorInterp=Gray
Minimum=238.364, Maximum=1652.747, Mean=684.315,
StdDev=268.118
Metadata:
STATISTICS_MINIMUM=238.36408996582
STATISTICS_MAXIMUM=1652.7473144531
STATISTICS_MEAN=684.31510635898
STATISTICS_MEDIAN=628.00527274609
STATISTICS_MODE=506.19808101654
STATISTICS_STDDEV=268.11784359311
STATISTICS_HISTONUMBINS=256
STATISTICS_HISTOMIN=240.4369354248
STATISTICS_HISTOMAX=1657.8297119141
LAYER_TYPE=athematic
STATISTICS_HISTOBINVALUES=5200|20332|20380|24760|23240|16028
|7444|4828|4632|7912|18480|38392|44548|44292|52980|55312|530
24|60224|56604|51976|49164|50864|45892|47932|47932|48004|458
88|46656|49592|50648|49180|51948|53668|53692|51920|54680|595
28|65148|73368|85132|88648|79764|73992|73216|72108|74412|764
88|99608|103600|95308|78448|73684|68140|70020|72932|75032|79
512|82636|83640|77584|75844|77228|66504|59480|58160|55356|55
064|54572|53656|52912|51696|51936|52816|52924|53996|52560|52
192|51984|51144|50164|48892|48960|49412|49352|48220|49576|50
632|49420|48748|48368|46484|47664|46132|47004|46636|45124|43
416|41900|40692|40824|39976|40188|41860|40652|40368|39692|39
740|40220|38624|37064|37240|37096|37844|39032|39880|40800|40
172|38876|37964|39788|38152|37716|37064|34944|34712|33448|33
896|33800|33420|33012|32080|30868|30976|29740|30700|29972|29
664|29712|29248|29088|27892|27036|26092|24756|24736|25132|24
360|24300|23608|22616|21032|20160|19760|18748|18616|18964|18
400|17624|16692|16348|15628|15376|14644|14520|14064|13992|14
160|15048|16600|15976|14892|13428|12740|12360|11868|11520|10
672|9772|9344|9000|8672|8732|8520|9020|8604|8052|8136|8240|8
036|7824|7844|8280|7816|7624|7180|7308|6744|7024|6880|6800|6
904|6336|6056|6120|5968|6008|5936|5388|5176|4996|5020|4520|4
232|4288|4184|3856|3820|3752|3420|3632|3472|3316|3448|3316|3
120|2924|2816|2652|2512|2296|1940|1624|1336|1260|1252|936|66
8|568|576|548|484|512|456|412|376|360|264|288|168|112|104|76
|72|12|16|4|
--
John J. Mitchell
More information about the MapServer-users
mailing list