+- Forum QGIS (http://forum.quantum-gis.pl)
+-- Dział: Desktop GIS (http://forum.quantum-gis.pl/forum-4.html)
+--- Dział: QGIS (http://forum.quantum-gis.pl/forum-5.html)
+--- Wątek: skrypt służący do konwersji plików (/thread-413.html)
skrpyt służący do konwersji plików - mateokkk - 19-03-2011
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
RE: skrpyt służący do konwersji plików - borys - 19-03-2011
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
RE: skrpyt służący do konwersji plików - matooook - 09-04-2011
(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
