I têm uma enorme quantidade de dados na forma de CSV contendo coordenadas UTM como X
e Y
e um valor de elevação como Z
informação. Preciso converter esses dados em um DEM como GeoTIFF para análise posterior. Nesse caso, uma quantidade enorme significa 16 m. linhas, com um ponto em X
, Y
e Z
por linha. Os pontos são igualmente distribuídos, portanto, não é necessária uma interpolação; cada ponto só precisa ser convertido em uma célula raster.
Os dados originais vieram sem separador, com larguras de coluna fixas. Eu já descobri como converter a sintaxe do arquivo para usar um separador em vez de larguras fixas e eliminar todos os caracteres de espaço, usando o editor de texto stream sed . A partir daqui, normalmente meu trabalho seria para importar os dados para o ArcGIS através da criação de uma classe de recurso dos X
, Y
e Z
de dados e, numa segunda etapa, a conversão do shapefile ponto em um GeoTIFF, usando o ponto de Raster ferramenta. No entanto, o arquivo que tenho atualmente é muito grande para esse processo.
Em vez do fluxo de trabalho descrito acima, eu estava procurando uma alternativa eficiente e descobri o GDAL. No entanto, em gdal_translate
, o formato suportado mais próximo que posso encontrar na lista de tipos de arquivos suportados é a grade ASCII, mas sem XYZ separado por vírgula. Outra dificuldade é que eu tenho coordenadas UTM , enquanto a maioria dos exemplos parece usar coordenadas de graus decimais. No entanto, preciso permanecer no sistema UTM (ou pelo menos, minha saída GeoTIFF precisa estar em um sistema de coordenadas UTM).
Portanto, estou procurando uma maneira de converter o CSV XYZ em um GeoTIFF, usando GDAL , mas até agora não consegui encontrar exemplos que tratassem desse problema exato. Eu ficaria muito feliz por algumas dicas ou mesmo exemplos de código.