Gfoss -- Geographic Free and Open Source Software - Italian mailing list

raster e mapserver

‹ Previous Topic Next Topic ›
 
classic Classic list List threaded Threaded
♦
♦
 4 messages Options
Luca Delucchi
Reply | Threaded
Open this post in threaded view
♦
♦
|

raster e mapserver

Luca Delucchi
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
_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
Non inviate messaggi commerciali.
I messaggi di questa lista non rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
496 iscritti al 24.1.2011
luca.casagrande@gmail.com
Reply | Threaded
Open this post in threaded view
♦
♦
|

Re: raster e mapserver

luca.casagrande@gmail.com
Ciao Luca,
puoi postare il layer del mapfile?

L.

Il 08 febbraio 2011 11:39, Luca Delucchi <[hidden email]> 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
_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
Non inviate messaggi commerciali.
I messaggi di questa lista non rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
496 iscritti al 24.1.2011
Luca Delucchi
Reply | Threaded
Open this post in threaded view
♦
♦
|

Re: raster e mapserver

Luca Delucchi
Il 08 febbraio 2011 11:45, [hidden email]
<[hidden email]> 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
_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
Non inviate messaggi commerciali.
I messaggi di questa lista non rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
496 iscritti al 24.1.2011
Marco Stelluti
Reply | Threaded
Open this post in threaded view
♦
♦
|

Re: raster e mapserver

Marco Stelluti
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

_______________________________________________
Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
[hidden email]
http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
Questa e' una lista di discussione pubblica aperta a tutti.
Non inviate messaggi commerciali.
I messaggi di questa lista non rispecchiano necessariamente
le posizioni dell'Associazione GFOSS.it.
496 iscritti al 24.1.2011
Free forum by Nabble Edit this page