Mam ogromną ilość danych w postaci CSV zawierających współrzędne UTM jako X
i Y
oraz wartość wysokości jako Z
informację. Muszę przekonwertować te dane do DEM jako GeoTIFF do dalszej analizy. W tym przypadku ogromna ilość oznacza 16 m. linie, z jednym punkcie X
, Y
a Z
na linię. Punkty są równomiernie rozmieszczone, dlatego interpolacja nie jest potrzebna; każdy punkt wystarczy przekształcić w komórkę rastrową.
Oryginalne dane były dostarczane bez separatora, ze stałą szerokością kolumn. Już wymyśliłem, jak przekonwertować składnię pliku, aby użyć separatora zamiast stałych szerokości i wyeliminować wszystkie znaki spacji, używając edytora tekstu strumieniowego sed . Odtąd, normalnie mój workflow byłoby importować dane do ArcGIS tworząc klasę elementów z X
, Y
i Z
danych, jak i drugim etapie przekształcania punkt shapefile do GeoTIFF, korzystając z punktu do Raster narzędzia. Jednak plik, który mam obecnie, jest zdecydowanie za duży dla tego procesu.
Zamiast wyżej opisanego przepływu pracy szukałem wydajnej alternatywy i odkryłem GDAL. Jednak gdal_translate
najbliższym obsługiwanym formatem, jaki mogę znaleźć na liście obsługiwanych typów plików, jest siatka ASCII, ale XYZ nie jest oddzielony przecinkami. Inną trudnością jest to, że mam współrzędne UTM , podczas gdy większość przykładów wydaje się używać współrzędnych dziesiętnych. Muszę jednak pozostać w systemie UTM (a przynajmniej mój wynik GeoTIFF musi znajdować się w układzie współrzędnych UTM).
Dlatego szukam sposobu na konwersję CSV XYZ na GeoTIFF, używając GDAL , ale jak dotąd nie byłem w stanie znaleźć przykładów dotyczących tego dokładnie problemu. Byłbym bardzo szczęśliwy z powodu niektórych wskazówek, a nawet przykładów kodu.