J'ai une énorme quantité de données sous la forme de CSV contenant les coordonnées UTM en Xet Yet une valeur d'élévation comme Zinformation. J'ai besoin de convertir ces données en DEM en tant que GeoTIFF pour une analyse plus approfondie. Dans ce cas, une énorme quantité signifie 16 m. lignes, avec un point dans X, Yet Zpar ligne. Les points sont également répartis, donc une interpolation n'est pas nécessaire; chaque point doit simplement être converti en cellule raster.
Les données d'origine étaient fournies sans séparateur, avec des largeurs de colonne fixes. J'ai déjà compris comment convertir la syntaxe du fichier pour utiliser un séparateur au lieu de largeurs fixes et éliminer tous les caractères d'espace, en utilisant l'éditeur de texte de flux sed . A partir de là, normalement mon flux de travail serait d'importer les données dans ArcGIS en créant une classe d'entités des X, Yet des Zdonnées et en tant que deuxième étape, la conversion du point shapefile en GeoTIFF, en utilisant le point à Raster outil. Cependant, le fichier que j'ai actuellement est trop volumineux pour ce processus.
Au lieu du workflow décrit ci-dessus, je cherchais une alternative efficace et j'ai découvert GDAL. Cependant, dans gdal_translate, le format pris en charge le plus proche que je peux trouver dans la liste des types de fichiers pris en charge, est la grille ASCII mais pas de XYZ séparé par des virgules. Une autre difficulté est que j'ai des coordonnées UTM , alors que la plupart des exemples semblent utiliser des coordonnées en degrés décimaux. Cependant, je dois rester dans le système UTM (ou au moins, ma sortie GeoTIFF doit être dans un système de coordonnées UTM).
Je cherche donc un moyen de convertir le CSV XYZ en GeoTIFF, en utilisant GDAL , mais jusqu'à présent je n'ai pas pu trouver d'exemples traitant de ce problème exact. Je serais très heureux pour quelques conseils ou même des exemples de code.