Tengo una gran cantidad de datos en forma de CSV que contiene coordenadas UTM como Xy Yy un valor de elevación como Zinformación. Necesito convertir estos datos en un DEM como GeoTIFF para su posterior análisis. En este caso, una gran cantidad significa 16 m. líneas, con uno de cada punto X, Yy Zpor línea. Los puntos están distribuidos equitativamente, por lo tanto, no se necesita una interpolación; cada punto solo necesita convertirse en una celda ráster.
Los datos originales llegaron sin separador, con anchos de columna fijos. Ya descubrí cómo convertir la sintaxis del archivo para usar un separador en lugar de anchos fijos y eliminar todos los caracteres de espacio, usando el editor de texto de flujo sed . A partir de aquí, normalmente mi flujo de trabajo sería para importar los datos en ArcGIS mediante la creación de una clase de entidad de los X, Yy Zlos datos y como segundo paso, convirtiendo el shapefile de puntos en un GeoTIFF, utilizando el punto de la trama herramienta. Sin embargo, el archivo que tengo actualmente es demasiado grande para este proceso.
En lugar del flujo de trabajo descrito anteriormente, estaba buscando una alternativa eficiente y descubrí GDAL. Sin embargo, en gdal_translate, el formato compatible más cercano que puedo encontrar en la lista de tipos de archivos admitidos es la cuadrícula ASCII pero no XYZ separados por comas. Otra dificultad es que tengo coordenadas UTM , mientras que la mayoría de los ejemplos parecen usar coordenadas de grados decimales. Sin embargo, necesito permanecer dentro del sistema UTM (o al menos, mi salida GeoTIFF debe estar en un sistema de coordenadas UTM).
Entonces, estoy buscando una manera de convertir el CSV XYZ en un GeoTIFF, usando GDAL , pero hasta ahora no pude encontrar ejemplos que trataran este problema exacto. Estaría muy feliz por algunos consejos o incluso ejemplos de código.