rigacci.org

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
====== 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**).

<code>
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
</code>

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:

<code>
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
</code>

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:

<code xml>
<?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>
</code>

===== Query su raster =====

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

<code xml>
<?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>
</code>

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

<code xml>
<?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>
</code>