Differences

This shows you the differences between two versions of the page.

--- doc:appunti:linux:lezioni:mapserver [2010/08/12 18:58] – niccolo
+++ doc:appunti:linux:lezioni:mapserver [2013/10/04 11:22] (current) – [gdaltindex] niccolo
@@ Line 75: / Line 75: @@
 MAP
   ...
-  IMAGETYPE agg_png
+  IMAGETYPE aggpng24
   ...
   OUTPUTFORMAT
-    NAME "agg_png"
+    NAME "aggpng24"
     DRIVER AGG/PNG
     IMAGEMODE RGB
@@ Line 235: / Line 235: @@
 </code>
-Il file indice è uno shapefile, analogamente a quanto visto con ''tile4ms''. Va indicato con la direttiva **''TILEINDEX''** che si intende **relativa a ''SHAPEPATH''**. L'indice deve contenere il percorso dei file raster relativi a ''SHAPEPATH'', non al file indice.
+Il file indice è uno shapefile, analogamente a quanto visto con ''tile4ms''. Va indicato con la direttiva **''TILEINDEX''** che si intende **relativa a ''SHAPEPATH''**. L'indice deve contenere il percorso dei file raster relativi a ''SHAPEPATH'', non al file indice. Con **''OFFSITE''** si indica quale colore deve essere considerato trasparente:
 <file>
@@ Line 246: / Line 246: @@
   TILEINDEX "toscana_ctr25k_index.dbf"
   TILEITEM "Location"
+  PROCESSING "RESAMPLE=AVERAGE"
+  OFFSITE 255 255 255
+END
+</file>
+==== gdalbuildvrt ====
+In alternativa all'indice di raster si può creare un **virtual raster** (che può essere agevolmente aperto con tutti i programmi che utilizzano una moderna libreria GDAL, come QGIS):
+<code>
+cd /usr/local/share/gisdata
+gdalbuildvrt toscana_ctr25k.vrt raster/toscana_ctr25k/*.tif
+</code>
+In questo caso il mapfile ha questa forma:
+<file>
+SHAPEPATH "/usr/local/share/gisdata"
+LAYER
+  NAME "ctr25k"
+  STATUS DEFAULT
+  TYPE RASTER
+  DATA toscana_ctr25k.vrt
   PROCESSING "RESAMPLE=AVERAGE"
   OFFSITE 255 255 255
@@ Line 281: / Line 305: @@
 gdaltindex 25k_retiled_index.shp 25k_retiled/*.tif
 </code>
+===== Ricetta per creazione di un servizio WMS =====
+  - Creare indice (quadro d'unione) dei raster:<code>gdaltindex 2002_index.shp 2002/*.tif</code>
+  - Creare le piramidi (overview). Per un singolo file .tif: <code>
+gdaladdo --config USE_RRD YES 287050.tif 2 4 8 16 32 64</code>
+  - Oppure per tutti i .tif in una directory: <code>ls *.tif | xargs -I nome_file gdaladdo --config USE_RRD YES nome_file 2 4 8 16 32 64</code>
+  - Con GDAL > 1.6.0 vedere ''man gdaladdo'' per il parametro **''%%-ro%%''** e **''%%-r gauss%%''** per creare file esterni e fare il resample di qualità migliore.
+  - Creare il file .map:
+    - Controllare l'estensione complessiva <code>ogrinfo 2002_index.shp -summary 2002_index</code>
+    - Verificare che esistano almeno i seguenti parametri <file>
+MAP.EXTENT                       1672800 4843801 1688801 4856418
+MAP.IMAGETYPE                    jpeg
+MAP.PROJECTION                   "+init=epsg:3003"
+WEB.METADATA.wms_title           "WMS Ortofoto"
+WEB.METADATA.wms_onlineresource  "http://192.168.2.10/cgi-bin/mapserv?map=..."
+WEB.METADATA.wms_srs             "EPSG:3003"
+LAYER.STATUS                     OFF
+LAYER.PROJECTION                 "+init=epsg:3003"
+LAYER.METADATA.wms_title         "Ortofoto 2001"
+LAYER.METADATA.wms_srs           "EPSG:3003"
+</file>
 ===== Legende =====