classify raster through WMS

Question

I try to use WMS to dispaly raster map in geotiff format on web. I want to classify the raster file. how can I do this? I use mapserver for windows. The following is my .map file. what I get through this is

MAP
NAME PM10
IMAGECOLOR 255 255 255
SIZE 600 800
IMAGETYPE PNG24 ## use AGG to for anti-aliassing
OUTPUTFORMAT
NAME 'AGG'
DRIVER AGG/PNG
MIMETYPE "image/png"
IMAGEMODE RGB
EXTENSION "png"
END # outputformat
PROJECTION
"init=epsg:3035" #latlon on etrs 1989 laea
 END
 EXTENT 3487500 2297500 4402500 3202500 # meters extents of region2
 WEB
 IMAGEPATH "c:/tmp/ms_tmp/"
 IMAGEURL "/ms_tmp/"
 METADATA
  "ows_enable_request"   "*"
  "map" "C:/ms4w/apps/airpollution/config.map"
  "ows_schemas_location" "http://schemas.opengeospatial.net"
  "ows_title" "Sample WMS"
  "ows_enable_request"   "*"
  "ows_onlineresource" "http://localhost:7070/cgi-bin/mapserv.exe?       map=C:/ms4w/apps/airpollution/config.map&"
  "ows_srs" "EPSG:3035 " #latlon      
  "wms_feature_info_mime_type" "text/plain"
  "wms_feature_info_mime_type" "text/html"
  "wms_server_version" "1.1.1"
  "wms_formatlist" "image/png,image/gif,image/jpeg, image/geotiff"
  "wms_format" "image/geotiff"
 END #metadata
 END #web
 LAYER
 NAME "pm10"
DATA "pm10.tif"
 TYPE RASTER
 STATUS ON
    METADATA
  "ows_title" "pollution"
 END #metadata
 PROJECTION
  "init=epsg:3035"
END #projection
CLASSITEM "[pixel]"
 # class using simple string comparison, equivelent to ([pixel] = 0)
 CLASS
 EXPRESSION "0"
 STYLE
  COLOR 20 20 20
 END
 END
# class using an EXPRESSION using only [pixel].
CLASS
EXPRESSION ([pixel] >= 0AND [pixel] < 7)
STYLE
  COLOR 255 0 0
END
CLASS
EXPRESSION ([pixel] >= 7AND [pixel] < 20)
STYLE
  COLOR 0 255 0
END
END
CLASS
EXPRESSION ([pixel] >= 7AND [pixel] < 50)
STYLE
  COLOR 0 0 255
END
END


END #layer pm10
END #map

and what I get as a responce is this image

it seems that it only reads line 3

IMAGECOLOR 255 255 255

Answer 1

I don't know map server very well, but I know that you can style Rasters from a WMS using SLD (Styled Layer Descriptor), which is simply a xml file that you can pass in a WMS request according to OGC standards.

In other words: You can specify the classification in a XML schema. Following is an example of a simple SLD that styles everything in a raster black except white pixels, which are styled opaque.

    <?xml version="1.0" encoding="ISO-8859-1"?>
<StyledLayerDescriptor version="1.0.0" xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.opengis.net/sld http://schemas.opengis.net/sld/1.0.0/StyledLayerDescriptor.xsd">
    <NamedLayer>
        <Name>undefined</Name>
        <UserStyle>
            <Name>rasterr</Name>
            <Title>Rasterr</Title>
            <Abstract>A simple raaster style</Abstract>
            <FeatureTypeStyle>
                <FeatureTypeName>Feature</FeatureTypeName>
                <Rule>
                    <RasterSymbolizer>
                        <Opacity>1.0</Opacity>
                        <ColorMap>
                            <ColorMapEntry color="#ffffff" quantity="0" opacity="0.0" />
                            <ColorMapEntry color="#000000" quantity="1" />
                            <ColorMapEntry color="#000000" quantity="2" />
                            <ColorMapEntry color="#000000" quantity="3" />
                            <ColorMapEntry color="#000000" quantity="4" />
                            <ColorMapEntry color="#000000" quantity="5" />
                            <ColorMapEntry color="#000000" quantity="6" />
                            <ColorMapEntry color="#000000" quantity="7" />
                            <ColorMapEntry color="#000000" quantity="8" />
                            <ColorMapEntry color="#000000" quantity="9" />
                            <ColorMapEntry color="#000000" quantity="10" />
                        </ColorMap>
                    </RasterSymbolizer>
                </Rule>
            </FeatureTypeStyle>
        </UserStyle>
    </NamedLayer>
</StyledLayerDescriptor>

Pass the SLD like this:

http://demo.mapserver.org/cgi-bin/wms?SERVICE=wms&VERSION=1.1.1&REQUEST=GetMap&LAYERS=country_bounds&SLD=http://demo.mapserver.org/ogc-demos/map/sld/sld_line_simple.xml

Read more here:

http://mapserver.org/ogc/sld.html - This is for mapserver, use the RasterSymbolizer and a ColorMap to do your classification. On that page it is also described how the colormap works.

http://www.opengeospatial.org/standards/sld Wikipedia