Sunday, August 21, 2011
"Using transport bicycles such as this 1950 Schwinn cruiser instead of a car saved me enough money in nine years to pay for three years of overseas travel by bike, boat, and plane. It carries a comfortable payload of 80 lbs in lockable, waterproof panniers." from http://oasisdesign.net/design/principles.htm
Thursday, August 18, 2011
Sunday, August 14, 2011
Thursday, July 14, 2011
For those learning the Soil and Water Assessment Tool (SWAT)...
I suggest checking out this site.
http://www.smm.org/scwrs/tapwaters
http://www.smm.org/scwrs/tapwaters
Wednesday, June 29, 2011
Using r.buffer for watershed analysis
A problem that I ran into when using r.watershed with a masked basin was that the streams were being cut off, and the resulting basin was incorrect.
![]() |
Close up of the HUC 10 watershed boundary in black and the streams derived from r.watershed in blue and pink |
![]() |
basin resulting from original mask |
To solve this I created a buffer raster of the original raster being used as a mask.
r.buffer input=rioGrandeHUC10mask output=rioGrandeBuffer distances=100
Then, following the grass book
r.watershed -f elev_cm_rioGrande2 thresh=10000 accum=accum_10K drain=draindir_10K basin=basin_10K stream=rivers_10K --o
r.mapcalc "streams_der_30m=if(abs(accum_10K)>100,1,null())"
r.thin streams_der_30m out=streams_der_30m_tmp
r.to.vect -s streams_der_30m_tmp out=streams_der_30m
r.to.vect -s basin_10K out=basin_10K feature=area
Here is the resulting stream
![]() |
Streams in green after using a buffer with 100m distance |
r.water.outlet drainage=draindir_10K basin=outlet_O8 easting=170670 northing=270651 --o
r.mapcalc 'outlet_O8=if(outlet_O8==1,1,null())'
![]() |
Resulting basin with buffered mask |
Subscribe to:
Posts (Atom)