Ho un'enorme quantità di dati sotto forma di CSV contenente coordinate UTM come Xe Ye un valore di elevazione come Zinformazioni. Devo convertire questi dati in un DEM come GeoTIFF per ulteriori analisi. In questo caso, una quantità enorme significa 16 m. linee, con un punto X, Ye Zper linea. I punti sono equamente distribuiti, pertanto non è necessaria un'interpolazione; ogni punto deve solo essere convertito in una cella raster.
I dati originali sono arrivati senza separatore, con larghezze di colonna fisse. Ho già capito come convertire la sintassi del file per utilizzare un separatore invece di larghezze fisse ed eliminare tutti i caratteri dello spazio, utilizzando l'editor di testo di flusso sed . Da qui in poi, normalmente il mio flusso di lavoro consisterebbe nell'importare i dati in ArcGIS creando una classe di funzionalità da X, Ye Zdati e come secondo passo, convertendo il file shapefile in un GeoTIFF, usando lo strumento Point to Raster . Tuttavia, il file che ho attualmente è troppo grande per questo processo.
Invece del flusso di lavoro sopra descritto, stavo cercando un'alternativa efficiente e ho scoperto GDAL. Tuttavia, in gdal_translate, il formato supportato più vicino che posso trovare nell'elenco dei tipi di file supportati, è la griglia ASCII ma nessuna XYZ separata da virgola. Un'altra difficoltà è che ho le coordinate UTM , mentre la maggior parte degli esempi sembra utilizzare coordinate in gradi decimali. Tuttavia, devo rimanere all'interno del sistema UTM (o almeno, il mio output GeoTIFF deve essere in un sistema di coordinate UTM).
Quindi sto cercando un modo per convertire il CSV XYZ in un GeoTIFF, usando GDAL , ma finora non sono riuscito a trovare esempi che affrontino questo esatto problema. Sarei molto felice per alcuni suggerimenti o anche esempi di codice.