Associare il nome delle regioni ISTAT a dei punti demanio

A cura di Totò Fiandaca  | guida/e Andrea Borruso

Dataset

• 01_demanio.csv
• regioni.shp

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.