Associare il nome delle regioni ISTAT a dei punti demanio

  • autore: Totò Fiandaca

Procedimento

Dopo aver scaricato i due dataset occorre convertirli in shapefile oppure importarli in un database sqlite (con estensione spaziale), verificare che I DUE STRATI abbiano stesso EPSG (sistema di riferimento delle coordinate)

Query spaziale

CREATE TABLE opendemanio_3857_reg AS
SELECT a.*, r.DEN_REG
FROM opendemanio_3857 a, regioni_3857 r
WHERE ST_Intersects(a.geometry,r.geometry)
AND a.rowid IN (SELECT rowid 
			FROM SpatialIndex 
			WHERE f_table_name = 'opendemanio_3857' AND search_frame = r.geometry);

Mapshaper

mapshaper opendemanio_3857.shp -join regioni_3857.shp fields=COD_REG,DEN_REG -o opendemanio_v01.csv

Installare mapshaper

## per installarlo

sudo npm install -g mapshapersudo apt-get update
sudo apt-get install nodejs
sudo apt-get install npm
sudo npm install -g mapshaper

## per aggiornarlo:

sudo npm update -g mapshaper

Cosa fa questa ricetta

Associa il nome delle regioni italiane ISTAT ai punti (demanio) che vi ricadono dentro; la query spaziale si puo” usare sia in spatialite_gui che nella riga di comando, mentre la ricetta di mapshaper solo da riga di comando.