Pasakojimas¶
Po kelių bandymų man pavyko su QGIS sukurti rezultatą, panašų į arcHYDRO. Pradžioje USGS Hydroshed DEM buvo konvertuotas į Reguliarizuotas juosteles su įtampa naudojant QGIS GRASS funkcijas pagal Helenos Mitsovos pamokas. RST paviršius buvo išanalizuotas ieškant telkinių ir tėkmės kelių. Gauta vektorinė topologija buvo išvalyta QGIS’u generuojant Hydrologinį tinklą ir taip pat kuriant atributų lenteles.
Šis tinklas buvo reitinguotas naudojant GRASS plėtinį r.strem, kad būtų sukurti Hortono ir Hacko tėkmių rikiavimai. Tada atstumas iki paviršiaus gardelės buvo sugeneruotas ir išanalizuotas naudojant zoninę statistiką, ieškant tolimiausio taško. Tada, GRASS r.drain panaudotas atsekti koncentracijos linijos laiką. Ir pabaigoje GRASS v.net buvo panaudotas pylimo taškų tinklo sukūrimui, tekėjimo keliams, TOC keliams ir HMDP kiekvienam vandens telkiniui. Tokia darbų seka QGIS įrankio atributai pagal vietą buvo neįkainojami.
Pabaigus paprastąjį hidrologinį tinklą, aš naudojau UN FAO dirvožemio duomenų bazę jungiant lenteles su dirvožemio klasifikacija. Tada QGIS užklausų funkcijos FAO dirvožemius konvertavo į USDA dirvožemius. Anksčiau atsiųsta USGS ETM7+ gardelė buvo klasifikuota naudojant i.cluster ir r.maxlik į NLCD gardelę. Naudojant perklasifikavimo lenteles aš konvertavau žemės paviršiaus informaciją į Manning’o n-reikšmių gardelę viršžeminei tėkmei. Pabaigoje parašiau lentelę r.mapcalc ir suliejau dirvožemio ir žemės paviršiaus gardeles į NRCS TR-55 kreivių skaičius naudojant archHYDRO metodologiją.
Kad sukurčiau lietaus gardeles, vykdžiau GHCN lentelių statistinę analizę, radau lietaus gylius ir dizaino audras 1,2,5,10,50,100 metams. Tada atsisiunčiau nemokamą Kinijos PRISM gardelę ir darant prielaidą, kad priklausomybė yra linijinė, konvertavau vidutines celės reikšmes į maksimalias celės reikšmes, kad sukurčiau maksimalaus kritulių intensyvumo gardeles.