Wydaje się, że powinno to być proste, ale nie mogę tego zrozumieć w życiu.
Biorąc pod uwagę istniejący GeoTiff zawierający jeden zespół, jak możesz dodać kolejny zespół do tego GeoTiff? Czy jest to jedyny sposób, aby to zrobić przy użyciu GDALCreate (), aby utworzyć nowy plik wyjściowy z dodatkowym pasmem, a następnie skopiować istniejące pasmo z oryginalnego pliku do nowego pliku? To jedyne rozwiązanie, jakie do tej pory wymyśliłem, ale wydaje się, że powinien istnieć prosty sposób na dodanie pasma do istniejącego zestawu danych.
Odpowiedzi:
Rozwiązaniem, jeśli sterownik je obsługuje, jest wywołanie GDALOpen () z dostępem GA_Update, a następnie użycie GDALAddBand lub GDALDataset :: AddBand . Jednak sterownik geotiff nie obsługuje AddBand.
źródło
Aby rozwinąć odpowiedź Luke'a i podać konkretny przykład w Pythonie, oto fragment, który dodaje pasmo alfa do źródłowego rastra i zapisuje go jako PNG.
Użyłem
MEM
sterownika,VRT
ponieważ ten ostatni nie obsługujeWriteRaster()
iWriteArray()
(błąd „Zapis przez VRTSourcedRasterBand nie jest obsługiwany.”). Podejrzewam, że użycie sterownika vrt może być nadal możliwe za pomocą innych metod.źródło