Daugiau

GRASS scenarijus kelių juostų rastrų gamybai?


Aš stengiuosi įjungti kilpą r.saule funkcija „Python“ scenarijuje.

Bet kokiu atveju, atrodo, nesuprantu, kaip galiu sukrauti rastrinius žemėlapius, kuriuos kiekviena kilpa sudarytų į vieną kelių juostų rastrą (pageidautina geotifo formatu).

Manau, kad vienas iš būdų būtų naudoti sąrašą, sukurtą kartu su g.sąrašas arba (dar geriau) g.mlist kaip įvesties parametras i.grupe funkcija ir tada eksportuokite sukurtą grupę kaip kelių juostų failą (kaip parodyta čia esančiame pavyzdyje). Čia yra mano paprastas kodas, pakeistas iš šio puslapio.

# importuokite reikiamus modulius iš osgeo importo gdal import numpy kaip np importuoti žolę. scenarijus kaip žolė # paruoškite reikiamą įvestį dem = "[email protected]_test" doy = np.arange (1 366) # metų diena # apskaičiuokite nuolydžio ir kraštinių žemėlapius grass.run_command ("r.slope.aspect", elevation = dem, slope = "[email protected]_test", aspekt = "[email protected]_test") # i yra metų diena (DOY), kurią reikia naudoti bėgimui r.saulė d = "[email protected]_test", slopein = "[email protected]_test", insol_time = insolname, glob_rad = radname, day = d, step = 1) # eksportuokite išvestis į GTiff failus grass.run_command ("r.out .gdal ", input = radname +" @globrad_test ", output = radname)

Paskutinė eilutė kiekvieną žemėlapį eksportuos į geotifo failą po kiekvienos kilpos. Atvirkščiai, aš juos visus sudėčiau į tą patį kelių juostų failąn*m*pmatmenys, kurmirnyra mano vaizdų eilutės ir stulpeliai (žinoma, vienodo dydžio) irpišvesties žemėlapių skaičius po kilpų (tai yra 365, dienų skaičius per metus).


Reikia tik nedidelių pakeitimų:

Prieš eksportuodami turite sukurti grupę (t. Y. Naudoti i.group), tada eksportuoti į daugiasluoksnį failą, priskirdami grupės pavadinimą kaip įvestį skambutyje r.out.gdal.


Žiūrėti video įrašą: Ткачество поясов на бёрдо. Заправочные узоры. (Spalio Mėn 2021).