Forum QGIS

Witam,

chcę napisać skrypt w Pythonie, który ze wszystkich plików w formacie NetCDF, zapisze wybrane pasma do formatu ASCII GRID - nie znam jednak potrzebnych mi poleceń...
Czy ktoś może wrzucić przykładowy fragment skryptu wykonującego taką czynność

Chciałbym również, aby pliki były zapisywane w postaci trzech kolumn - współrzędne i wartość - ze względu na to, że nie znam odpowiedniej opcji w QGIS, napisałem skrypt Pythona przekształcający ASCII GRID do wymaganej postaci. Może istnieje jednak prostszy sposób

To akurat nie ma nic do QGISa, tylko do GDALa. Można to zrobić pythonowymi łączami gdala, ale chyba szybciej i prościej gotowym programem gdal_translate [1]. Uwagi dotyczące obsługi tych dwóch formatów przez gdala są pod linkiem [2].

[1] http://www.gdal.org/gdal_translate.html
[2] http://www.gdal.org/frmt_various.html

Kwestia wywołania gdal_translate w pętli to już w zasadzie zadanie systemu operacyjnego, ale w razie kłopotów zawsze można wołać go z Pythona

(19-03-2011, 22:36)borys napisał(a): [ -> ]Kwestia wywołania gdal_translate w pętli to już w zasadzie zadanie systemu operacyjnego, ale w razie kłopotów zawsze można wołać go z Pythona

Niestety wywoływanie daje różne efekty...
Zawsze używałem do tego celu w Pythonie modułu subprocess - niestety nie byłem w stanie odwołać się do wszystkich wtyczek GLDAS-a (mimo podawania ich bezwzględnej lokalizacji na dysku).

Rozwiązanie - użyłem polecenia print, po czym skopiowałem wyplute komendy i wkleiłem do konsoli