Történet¶
Próbálkozásaim eredményeképp sikerült ArcHydro szerű outputot elérnem QGIS-szel. Az USGS Hydroshed magasságmodelljeivel kezdtem a műveletet, melyet feszített spline módszerrel (regularized spline with tension) interpoláltam a QGIS-ben található GRASS eszközök segítségével, Helena Mitsova útmutatásai alapján. A spline felszínen végeztem a hidrológiai műveleket, mint vízgyűjtőlehatárolás, folyásirány meghatározása, majd az eredmény vektor topológiákat QGIS-ben tisztítottam a helyes hidrológiai hálózat létrehozásához, már leíró adatokkal egyűtt.
A folyók Horton és Hacks rendűségét a GRASS r.stream eszközével határoztam meg. Ezután egy szabályos rácshálót készítettem a torkolattávolsággal, majd területi statisztika eszközzel elemeztem a legtávolabbi pont megtalálása érdekében. Következett az r.drain GRASS eszköz az összefolyások meghatározásához. Végezetül pedig a v.net eszközt használva elkészülhetett a hálózat a torkolatokból, összefolyásokból, stb. Ebben a munkafolyamatban a QGIS attribute by location eszköze volt rendkívül hasznos.
Miután az alapvető hidrológiai hálózat térképezése elkészült a UN FAO talaj adatbázisának táblát csatoltam a talajok osztályozása érdekében. Utána QGIS lekérdezési funkciókat használtam fel a FAO talajok konvertálásra kerültek USDA talajok osztályaivá. A korábban letöltött USGS ETM7+ grideket osztályoztam i.cluser és r.maxlik GRASS utasítások segítségével NLCD grid formába. Újraosztályoztam a táblákat és konvertáltam a felszínborítási kategóriákat Manning-féle n-value grid formába a szárazföldi lefolyás meghatározására. Végül írtam egy kereső táblát r.mapcalc segítségével és egyesítettem a talaj és a felszínborítás grideket a NRCS (Natural Resources Conservation Service) TR-55 görbe, a lefolyás mennyiségének meghatározására az arcHYDRO eljárását követve.
A csapadékmennyiség gridek előállításához statisztikai elemzéseket végeztem a GHCN táblákon, hogy meghatározzam a csapadékmennyiségének alakulását a viharok tervezésének 1,2,5,10,50,100 éves előrejelzésében. Ezután letöltöttem egy ingyenes mintát a PRISM adatbázisból Kína területére gridben, feltételezve a lineáris kapcsolatot létrehoztam az átlagos mező értékek alapján maximum cella értékeket annak érdekében, hogy előállítsak egy maximum csapadék intenzitás grid térképet.