Повеќе

GRASS скрипти за производство на мулти-бенд растер?

GRASS скрипти за производство на мулти-бенд растер?


Се обидувам да го заокружам р.сонце функционира во скрипта на Пајтон.

Како и да е, ми се чини дека не разбирам како можам да ги соберам растерските мапи што секоја јамка ќе ги произведе во еден растер со повеќе опсези (по можност во геотифен формат).

Мислам дека еден начин ќе биде користење на листата произведена со е.листа или (уште подобро) g.mlist како влезен параметар за i.група функција, а потоа извезете ја креираната група како датотека со повеќе опсези (како што е прикажано на примерот овде). Еве го мојот едноставен код, изменет од оваа страница.

# увезете ги потребните модули од osgeo import gdal увезете numpy како np увезете трева.скрипт како трева # подгответе го потребниот влез dem = "[email protected]_test" doy = np.аранжман (1,366) # ден од годината # пресметајте наклони и аспекти grass.run_command ("r.slope.aspect", elevation = dem, slope = "[email protected]_test", aspect = "[email protected]_test") # i е денот од годината (DOY) што треба да се користи за трчање r.sun за d in doy: # сакаат да добијат глобално зрачење radname = "globrad" + str (d) # run the solar model model.run_command ("r.sun", flags = "s", elevin = dem, aspin = "[email protected]_test", slopein = "[email protected]_test", insol_time = insolname, glob_rad = radname, day = d, step = 1) # извезувај ги излезите во датотеки GTiff grass.run_command ("r.out .gdal ", влез = radname +" @globrad_test ", излез = radname)

Последниот ред ќе извезува секоја карта во датотека со геотиф по секоја јамка. Наместо тоа, би ги ставил сите во иста датотека со повеќе опсези соn*m*pдимензии, кадеминсе редови и колони на моите слики (секоја со иста големина, се разбира) истрбројот на излезни карти по јамките (односно 365, бројот на денови во една година).


Потребна е само мала измена:

Треба да ја генерирате групата пред да извезувате (т.е. да користите i.group), потоа да извезувате во повеќеслојна датотека со доделување на името на групата како влез во повикот r.out.gdal.


Погледнете го видеото: Ткачество поясов на бёрдо. Заправочные узоры.