User Tools

Site Tools


Sidebar

No ai soldati italiani all'estero

Indice

Eventi

Energia

Rigacci.Org usa energia elettrica da fonti rinnovabili, grazie al gruppo di acquisto Merci Dolci.

Merci Dolci - Energia Rinnovabile

Software libero!

Petizione contro i brevetti software

Faunalia: Soluzioni GIS professionali

Debian

www.gnu.org www.kernel.org

Altre libertà

MoVimento 5 Stelle Campi Bisenzio

tecnica:gps_cartografia_gis:getfeatureinfo

GetFeatureinfo: query su WMS

Alcuni server WMS supportano le query tramite la richieste di tipo GetFeatureinfo.

La query viene eseguita rispetto al punto (X, Y), che sono le coordinate in pixel della mappa nelle dimensioni e bounding box specificate. La query viene fatta nel dintorno del punto, con la tolleranza indicata nel mapfile.

Il risultato può essere composto in formato testo, http oppure GML. Noi proveremo con quest'ultimo.

Con MapServer è necessario configurare alcuni parametri nel mapfile. Nel nostro caso aggiungiamo la complicazione dovuta al fatto che il layer è in un sistema di riferimento (EPSG:3003) diverso da quello che useremo nella query (EPSG:4326).

LAYER
  NAME "entita_lineari"
  DATA "275040el"
  STATUS ON
  TYPE LINE
  TEMPLATE "void"
  TOLERANCEUNITS pixels
  TOLERANCE 10
  DUMP TRUE
  PROJECTION
    "+init=epsg:3003 +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68"
  END
  METADATA
    "wms_title"         "Entità lineari"
    "wms_srs"           "EPSG:4326"
    "gml_include_items" "all"
    #"wms_include_items" "all"
  END
END

Se il client effettua richieste in formato text/plain (non GML), pare che per ottenere risposta sia necessario il parametro wms_include_items.

La richiesta WMS assume la seguente forma:

http://localhost/cgi-bin/mapserv?map=/var/www/wms_query.map&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetFeatureinfo&SRS=EPSG:4326&BBOX=11.16398206,43.74906364,11.33432601,43.79898334&LAYERS=entita_lineari&QUERY_LAYERS=entita_lineari&FEATURE_COUNT=100&X=150&Y=100&WIDTH=1024&HEIGHT=480&INFO_FORMAT=gml

Questi i parametri utilizzati:

SERVICE=WMS
VERSION=1.1.1
REQUEST=GetFeatureinfo
SRS=EPSG:4326
BBOX=11.16398206,43.74906364,11.33432601,43.79898334
LAYERS=entita_lineari
QUERY_LAYERS=entita_lineari
FEATURE_COUNT=100
X=150
Y=100
WIDTH=1024
HEIGHT=480
INFO_FORMAT=gml

Attenzione: il numero di feature (oggetti) ritornati dipende dalla dimensione del dintorno del punto (la TOLERANCE del mapfile) e dal parametro FEATURE_COUNT (numero massimo di feature ritornati) della query.

Il risultato è qualcosa del tipo:

<?xml version="1.0" encoding="UTF-8"?>
<msGMLOutput
    xmlns:gml="http://www.opengis.net/gml"
    xmlns:xlink="http://www.w3.org/1999/xlink"
    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
    <entita_lineari_layer>
        <entita_lineari_feature>
            <gml:boundedBy>
                <gml:Box srsName="EPSG:4326">
                    <gml:coordinates>11.173595,43.788665 11.197719,43.790197</gml:coordinates>
                </gml:Box>
            </gml:boundedBy>
            <FOGLIO>275030</FOGLIO>
            <CODICE>0301</CODICE>
            <RECORD>150872</RECORD>
            <IDENTIF></IDENTIF>
        </entita_lineari_feature>
        <entita_lineari_feature>
            <gml:boundedBy>
                <gml:Box srsName="EPSG:4326">
                    <gml:coordinates>11.180460,43.788754 11.189372,43.789513</gml:coordinates>
                </gml:Box>
            </gml:boundedBy>
            <FOGLIO>275030</FOGLIO>
            <CODICE>0102</CODICE>
            <RECORD>151559</RECORD>
            <IDENTIF></IDENTIF>
        </entita_lineari_feature>
    </entita_lineari_layer>
</msGMLOutput>

Query su raster

In modo analogo si possono effettuare query su layer raster, questo un esempio di risultato:

<?xml version="1.0" encoding="ISO-8859-1"?>
<msGMLOutput
    xmlns:gml="http://www.opengis.net/gml"
    xmlns:xlink="http://www.w3.org/1999/xlink"
    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
    <ctr25k_layer>
        <ctr25k_feature>
            <gml:boundedBy>
                <gml:Box srsName="EPSG:3003">
                    <gml:coordinates>1681609.346909,4707664.375214 1681609.346909,4707664.375214</gml:coordinates>
                </gml:Box>
            </gml:boundedBy>
            <x>1681609.3</x>
            <y>4707664.4</y>
            <value_0>0</value_0>
            <value_list>0</value_list>
            <red>255</red>
            <green>255</green>
            <blue>255</blue>
        </ctr25k_feature>
    </ctr25k_layer>
</msGMLOutput>

Questo il risultato di una query su un layer TIFF di tipo Int16, il colore (red, green, blue) deriva da una classificazione:

<?xml version="1.0" encoding="ISO-8859-1"?>
<msGMLOutput
    xmlns:gml="http://www.opengis.net/gml"
    xmlns:xlink="http://www.w3.org/1999/xlink"
    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
    <rainfall_layer>
        <rainfall_feature>
            <gml:boundedBy>
                <gml:Box srsName="EPSG:4326">
                    <gml:coordinates>-1.000000,-1.000000 -1.000000,-1.000000</gml:coordinates>
                </gml:Box>
            </gml:boundedBy>
            <x>0.2503605</x>
            <y>44.916937</y>
            <value_0>902</value_0>
            <value_list>902</value_list>
            <class>5</class>
            <red>247</red>
            <green>6</green>
            <blue>123</blue>
        </rainfall_feature>
    </rainfall_layer>
</msGMLOutput>
tecnica/gps_cartografia_gis/getfeatureinfo.txt · Last modified: 2010/01/28 15:53 by niccolo