[Gfoss] raster e mapserver

Local Chapters GFOSS.it
1

Stamattina ho messo su un servizio wms composto da un singolo raster,
ma non riesco ad ottenere una color table corretta. ho provato tutte
le soluzioni proposte qui [0] ma al massimo riesco ad ottenere una
scala di grigi (con PROCESSING "SCALE=AUTO")
il geotiff è il seguente

gis@srvcarto:~$ gdalinfo winkler2003.tif
Driver: GTiff/GeoTIFF
Size is 468, 393
Coordinate System is:
PROJCS["WGS 84 / UTM zone 32N",
    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",9],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32632"]]
Origin = (612000.000000000000000,5157500.000000000000000)
Pixel Size = (250.000000000000000,-250.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Corner Coordinates:
Upper Left ( 612000.000, 5157500.000) ( 10d27'40.65"E, 46d33'42.23"N)
Lower Left ( 612000.000, 5059250.000) ( 10d26'17.17"E, 45d40'39.90"N)
Upper Right ( 729000.000, 5157500.000) ( 11d59'10.48"E, 46d31'55.45"N)
Lower Right ( 729000.000, 5059250.000) ( 11d56'20.12"E, 45d38'56.34"N)
Center ( 670500.000, 5108375.000) ( 11d12'22.41"E, 46d 6'27.56"N)
Band 1 Block=468x2 Type=Float64, ColorInterp=Gray
  Metadata:
    COLOR_TABLE_RULES_COUNT=4
    COLOR_TABLE_RULE_RGB_0=0.000000e+00 5.000000e+02 191 191 191 255 255 0
    COLOR_TABLE_RULE_RGB_1=5.000000e+02 1.000000e+03 255 255 0 0 255 0
    COLOR_TABLE_RULE_RGB_2=1.000000e+03 1.500000e+03 0 255 0 255 127 0
    COLOR_TABLE_RULE_RGB_3=1.500000e+03 2.500000e+03 255 127 0 255 0 0

avete qualche suggerimento?

[0] http://mapserver.org/input/raster.html#special-processing-directives

--
ciao
Luca

http://gis.cri.fmach.it/delucchi/
www.lucadelu.org

2

Ciao Luca,
puoi postare il layer del mapfile?

L.

Il 08 febbraio 2011 11:39, Luca Delucchi <lucadeluge@gmail.com> ha scritto:

Stamattina ho messo su un servizio wms composto da un singolo raster,
ma non riesco ad ottenere una color table corretta. ho provato tutte
le soluzioni proposte qui [0] ma al massimo riesco ad ottenere una
scala di grigi (con PROCESSING "SCALE=AUTO")
il geotiff è il seguente

gis@srvcarto:~$ gdalinfo winkler2003.tif
Driver: GTiff/GeoTIFF
Size is 468, 393
Coordinate System is:
PROJCS["WGS 84 / UTM zone 32N",
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",9],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32632"]]
Origin = (612000.000000000000000,5157500.000000000000000)
Pixel Size = (250.000000000000000,-250.000000000000000)
Metadata:
AREA_OR_POINT=Area
Corner Coordinates:
Upper Left ( 612000.000, 5157500.000) ( 10d27'40.65"E, 46d33'42.23"N)
Lower Left ( 612000.000, 5059250.000) ( 10d26'17.17"E, 45d40'39.90"N)
Upper Right ( 729000.000, 5157500.000) ( 11d59'10.48"E, 46d31'55.45"N)
Lower Right ( 729000.000, 5059250.000) ( 11d56'20.12"E, 45d38'56.34"N)
Center ( 670500.000, 5108375.000) ( 11d12'22.41"E, 46d 6'27.56"N)
Band 1 Block=468x2 Type=Float64, ColorInterp=Gray
Metadata:
COLOR_TABLE_RULES_COUNT=4
COLOR_TABLE_RULE_RGB_0=0.000000e+00 5.000000e+02 191 191 191 255 255 0
COLOR_TABLE_RULE_RGB_1=5.000000e+02 1.000000e+03 255 255 0 0 255 0
COLOR_TABLE_RULE_RGB_2=1.000000e+03 1.500000e+03 0 255 0 255 127 0
COLOR_TABLE_RULE_RGB_3=1.500000e+03 2.500000e+03 255 127 0 255 0 0

avete qualche suggerimento?

[0] http://mapserver.org/input/raster.html#special-processing-directives

--
ciao
Luca

--
Luca Casagrande
http://www.lucacasagrande.net
twitter: lucacasagrande

3

Il 08 febbraio 2011 11:45, luca.casagrande@gmail.com
<luca.casagrande@gmail.com> ha scritto:

Ciao Luca,
puoi postare il layer del mapfile?

        LAYER
                NAME "winkler"
                STATUS ON
                PROJECTION
                        "init=epsg:32632"
                END
                DUMP TRUE
                TYPE RASTER
# PROCESSING "BANDS=1,2,3"
# PROCESSING "COLOR_MATCH_THRESHOLD=4"
# PROCESSING "DITHER=YES"
                PROCESSING "SCALE=AUTO"
# PROCESSING "LUT=/home/gis/winkler_curve"
                DATA "/home/gis/winkler2003.tif"
                METADATA
                        "wms_title" "winkler2003"
                        "wms_srs" "EPSG:32632"
                END
        END

L.

--
ciao
Luca

http://gis.cri.fmach.it/delucchi/
www.lucadelu.org

4

Ciao,
io avevo risolto utilizzando un vrt file:
http://osgeo-org.1803224.n2.nabble.com/gdal-dev-Modify-an-image-color-table-td5221037.html

oppure se le classi dei colori all’interno dell’immagine non sono moltissime, puoi descriverle una per una sul map file tipo:

CLASS
NAME “NOME CLASSE”
EXPRESSION ([pixel]=4) #QUI POUI INSERIRE UNA QUALUNQUE ESPRESSIONE
COLOR 188 0 0 #COLORE RGB
END

Ciao
m.s