Nel dibattito acceso sulla lunghezza della costa italiana, una questione cruciale emerge con forza: la discrepanza tra i dati forniti dal governo italiano e quelli proposti da varie associazioni ambientaliste / enti ed istituzioni politiche dello stato. La chiave per comprendere questa discrepanza è nel sistema SID, uno strumento di raccolta dati governativo attivo da numerosi anni. Mentre diverse associazioni hanno presentato stime inferiori riguardo alla lunghezza della costa italiana, non hanno preso in considerazione le informazioni dettagliate e aggiornate raccolte dal sistema SID.
Il governo italiano, insieme al tavolo tecnico consultivo sulle concessioni balneari, ha effettuato un’analisi approfondita basandosi sui dati raccolti dal SID. Questo sistema ha una lunga storia nella raccolta di informazioni geografiche e topografiche sull’intero territorio italiano, inclusa la linea costiera. A differenza di altre fonti, il SID contiene dati che sono stati continuamente aggiornati e arricchiti nel corso degli anni, fornendo un quadro più completo e attuale della situazione.
Per chi è interessato ad analizzare autonomamente la lunghezza della costa, il SID offre la possibilità di scaricare questi dati tramite il suo GeoServer, seguendo le linee guida disponibili sul sito ufficiale GeoServer Documentation. Il file relativo alla linea costiera, di dimensioni superiori a 700 MB e contenente oltre 655 Multiline Features, rappresenta una fonte di dati preziosa e accurata.
Tuttavia, è importante notare che i dati sul GeoServer del SID sono forniti nel formato “EPSG:3857“. Per una corretta analisi, è necessario convertire questi dati nel formato “EPSG:4326“. E’ importante sapere che EPSG:3857 è una proiezione cilindrica che distorce le dimensioni e le distanze, specialmente vicino ai poli, in particolare per calcoli su lunghe distanze, su scale globali, la distorsione può portare a risultati imprecisi.
Nelle prossime sezioni, forniremo una guida su come effettuare questa conversione e calcolare autonomamente la lunghezza della linea costiera italiana, utilizzando i dati forniti dal SID. Questo processo permetterà a chiunque di verificare in modo indipendente la lunghezza della costa italiana, offrendo una maggiore trasparenza e precisione nelle discussioni pubbliche riguardanti la gestione delle concessioni balneari e la conformità alle direttive europee.
1. Download della Linea di Costa in Formato GeoJSON
Il primo passo per analizzare la linea costiera italiana consiste nel scaricare i dati dal GeoServer del SID. È essenziale selezionare il Layer corretto per il download: il Layer “LINEA_DI_COSTA” contiene le informazioni necessarie per questa analisi. Questo Layer in formato GeoJSON fornisce una rappresentazione della linea costiera italiana.
Non forniamo direttamente i link per il download di questi dati. Tuttavia, chi è interessato può procedere autonomamente, con l’opzione di affidarsi a figure professionali nel campo della geografia o della cartografia per un’assistenza qualificata. Inoltre, il team di https://mappaturaspiagge.it/ è disponibile per fornire assistenza nel download dei dati e degli script di calcolo.
2. Divisione del File in Singole Feature
Una volta scaricato il file, il passo successivo è dividerlo in singole feature. Tale processo è particolarmente importante, dato che il file della linea costiera è estremamente grande, con dimensioni superiori a 700 MB. Questa grandezza può rendere difficoltoso aprire e gestire il file su molti computer. Pertanto, scomporre il file in singole features non solo consente di esaminare e analizzare in dettaglio ogni segmento della linea costiera, ma rende anche il processo più gestibile e meno oneroso in termini di risorse computazionali.
Questo approccio semplificato ci permette di affrontare l’analisi in modo più efficiente, facilitando la comprensione di ogni aspetto della complessa geografia costiera italiana.
Esempio (divisione features):
import geojson
from pyproj import Transformer
import os
def transform_feature_geometry(geometry, transformer):
if geometry['type'] == 'Point':
return {'type': 'Point', 'coordinates': transformer.transform(*geometry['coordinates'])}
elif geometry['type'] in ['LineString', 'MultiPoint']:
return {'type': geometry['type'], 'coordinates': [transformer.transform(*coord) for coord in geometry['coordinates']]}
elif geometry['type'] in ['Polygon', 'MultiLineString']:
return {'type': geometry['type'], 'coordinates': [[transformer.transform(*coord) for coord in ring] for ring in geometry['coordinates']]}
elif geometry['type'] == 'MultiPolygon':
return {'type': 'MultiPolygon', 'coordinates': [[[transformer.transform(*coord) for coord in ring] for ring in polygon] for polygon in geometry['coordinates']]}
def convert_and_save_features(geojson_file, output_folder):
try:
with open(geojson_file, 'r') as file:
data = geojson.load(file)
transformer = Transformer.from_crs("EPSG:3857", "EPSG:4326", always_xy=True)
if not os.path.exists(output_folder):
os.makedirs(output_folder)
for idx, feature in enumerate(data['features']):
transformed_geometry = transform_feature_geometry(feature['geometry'], transformer)
transformed_feature = geojson.Feature(geometry=transformed_geometry, properties=feature['properties'])
output_filename = os.path.join(output_folder, f'feature_{idx}.geojson')
with open(output_filename, 'w') as f:
geojson.dump(transformed_feature, f)
return f"Conversion complete. Files saved in {output_folder}"
except Exception as e:
return f"An error occurred: {e}"
geojson_file = 'geo_json-LINEA_DI_COSTA.json'
output_folder = 'features'
print(convert_and_save_features(geojson_file, output_folder))
3. Analisi delle Singole Feature
L’ultimo passo è l’analisi di ogni singola feature che ci permetterà, di calcolare con precisione la lunghezza totale della costa italiana, fornendo una base solida per eventuali discussioni politiche o ambientali relative alle concessioni balneari.
Esempio Calcolo:
import geojson
import os
from geopy.distance import geodesic
def calculate_length(geometry):
if geometry['type'] == 'LineString':
return sum(geodesic(geometry['coordinates'][i], geometry['coordinates'][i+1]).meters for i in range(len(geometry['coordinates']) - 1))
elif geometry['type'] == 'MultiLineString':
return sum(calculate_length({'type': 'LineString', 'coordinates': line}) for line in geometry['coordinates'])
def process_folder(folder_path):
total_length = 0
for filename in os.listdir(folder_path):
if filename.endswith('.geojson'):
with open(os.path.join(folder_path, filename), 'r') as file:
data = geojson.load(file)
if data['geometry']['type'] in ['LineString', 'MultiLineString']:
length = calculate_length(data['geometry'])
total_length += length
print(f"{filename}: Length = {length} meters, Properties = {data['properties']}")
return total_length
folder_path = 'features'
total_length = process_folder(folder_path)
print("Total length of all features (in meters):", total_length)
Questi passaggi, se eseguiti correttamente, offrono a ricercatori e al pubblico in generale gli strumenti per comprendere e valutare in modo indipendente la lunghezza della linea costiera italiana adottata dal SID.
500+ .....
feature_500.geojson: Length = 53221.69567599052 meters, Properties = {'id': 544}
feature_478.geojson: Length = 6747.467253626024 meters, Properties = {'id': 587}
feature_345.geojson: Length = 2253.1096425180676 meters, Properties = {'id': 334}
feature_355.geojson: Length = 6221.821254350592 meters, Properties = {'id': 354}
feature_327.geojson: Length = 7682.529944461105 meters, Properties = {'id': 340}
feature_337.geojson: Length = 35798.59247802161 meters, Properties = {'id': 435}
feature_572.geojson: Length = 11103.16533568514 meters, Properties = {'id': 573}
feature_562.geojson: Length = 5083.714916052935 meters, Properties = {'id': 633}
feature_614.geojson: Length = 11806.774408349127 meters, Properties = {'id': 228}
feature_604.geojson: Length = 19605.606135198734 meters, Properties = {'id': 403}
feature_139.geojson: Length = 42757.38472699991 meters, Properties = {'id': 157}
feature_129.geojson: Length = 70853.07589252737 meters, Properties = {'id': 71}
feature_180.geojson: Length = 11536.011893451025 meters, Properties = {'id': 110}
feature_64.geojson: Length = 2375.7287512498037 meters, Properties = {'id': 225}
feature_74.geojson: Length = 44679.59965741988 meters, Properties = {'id': 56}
feature_190.geojson: Length = 7722.530326315258 meters, Properties = {'id': 111}
Total length of all features (in meters): 12841807.659889873Conclusione: La Verità nei Numeri
Dopo un’analisi approfondita e tecnica della linea costiera italiana, utilizzando il file GeoJSON fornito dal SID, emergono dati che gettano una nuova luce sulla controversia della lunghezza della costa.
Il calcolo effettuato dal file della linea di costa del SID rivela una lunghezza totale di 12.841.807,66 metri, equivalenti a 12.841,80 chilometri. Questo dato è significativamente maggiore rispetto alle stime proposte dalle varie associazioni di categoria. È plausibile che il tavolo tecnico, nel suo impegno per fornire una stima accurata, abbia ulteriormente affinato questo calcolo, sicuramente analizzando i dati con maggiore dettaglio.
La nostra comprensione è che questi dati siano stati valutati con precisione. Il tavolo tecnico ha condotto le sue analisi basandosi sui dati forniti dal SID, considerato da tempo la fonte ufficiale per tali informazioni. Ciononostante, è importante notare che potrebbero esistere altri modelli o fonti di dati, che si basano su linee di costa non aggiornate o che adottano criteri di valutazione differenti.
Siamo aperti al dialogo e al confronto con altre fonti di dati per approfondire ulteriormente questa importante questione. Il nostro obiettivo principale è promuovere una discussione basata su dati accurati e trasparenti, al fine di guidare decisioni informate per la gestione delle nostre coste.
Per chi desidera eseguire analisi indipendenti o ulteriori studi, il pacchetto dati è disponibile per il download contattando il servizio di https://mappaturaspiagge.it/contatti
Massimo Chiarillo
Founder MappaturaSpiagge.it
[email protected]
Angelo Giulio Guerrieri
Referente Relazioni MappaturaSpiagge.it
Consulenza Richiesta Da:
Mauro Della Valle
Presidente Confimprese Demaniali Italia
[email protected]
Lascia un commento