Het verhaal¶
Met enig proberen en fouten maken was ik in staat een uitvoer, die leek op arcHYDRO, uit QGIS te produceren. Beginnend met een USGS Hydrosheds DEM die werd geconverteerd naar een Regularized Spline with Tension met behulp van de functies voor GRASS in QGIS en training door Helena Mitsova. Met oppervlakte RST werd dit oppervlak geanalyseerd om de opvanggebieden en stroompaden te zoeken, en resulterende vectortopologiën werden geschoond in QGIS om het hydrologisch netwerk te genereren en ook om de attributentabellen op te bouwen.
Dit netwerk werd gerangschikt door de extensie r.stream in GRASS om de volgorde voor stromen volgens build Hortons en Hacks te bouwen. Daarna werd een ratser voor de oppervlakte van afstand tot uitvoer gegenereerd en geanalyseerd door gebiedsstatistieken om het verste punt te zoeken. Daarna r.drain in GRASS om een lijn voor tijd van concentratie (TOC) te traceren. Tenslotte werd v.net in GRASS gebruikt om een netwerk van gietpunten, stroompaden, TOC-paden en HDMP te bouwen voor elke waterberging. In deze werkstroom was het gereedschap van QGIS Attributen per plaats van onschatbare waarde.
Nadat het basis hydrologische netwerk was voltooid, gebruikte ik UN FAO soils database met samenvoegingen van tabellen voor het classificeren van de bodem. Daarna werden met behulp van functies voor query in QGIS de FAO-bodems geconverteerd naar USDA-bodems. Een eerder gedownload raster van USGS ETM7+ werd geclassificeerd door i.cluser en r.maxlik in GRASS naar een raster voor NLCD. Door tabellen voor hereclassificatie te gebruiken converteerde ik het landgebruik naar het raster Manning’s n-value grid voor de stroming over land. Tenslotte schreef ik een opzoektabel in r.mapcalc en voegde de rasters voor bodem en landgebruik samen naar NRCS TR-55 boognummers gevolgd door de methodologie arcHYDRO
Ik voerde een statistische analyse uit op GHCN-tabellen om de mate van regenval voor terugkerende ontworpen stormen voor 1, 2,5 ,10, 50, 100 jaar te zoeken om de raster voor de neerslag te maken. Daarna downloadde ik een gratis voorbeeldraster van PRISM voor China en, uitgaande van een lineaire relatie,, converteerde gemiddelde celwaarden naar maximale celwaarden om de rasters voor maximale intensiteit van neerslag te bouwen.
