Michael Barton wrote:
On Jan 15, 2009, at 5:22 AM, <grass-dev-requ...@lists.osgeo.org> wrote:
[...]
Comment (by mmetz):
I have submitted a new version of r.watershed to trunk with various
changes
[...]
Markus,
Can this go into develbranch_6 too?
I hope so:-) But I would like to get the new functionality tested by
others before it is backported in case I missed something. You mentioned
that you are testing the new functionality, and you are using Mac while
I use Linux 64bit, so if the new version runs fine and produces proper
results on your system that would be encouraging!
Below are test scripts for the North Carolina sample dataset nc_spm_08
and the Spearfish dataset spearfish60. You can copy the commands for
each dataset into a script and then run the script, after updating to
today's version of r.watershed:-)
Markus M
Testing commands for the elevation map elev_lid792_1m in the North
Carolina sample dataset nc_spm_08. SFD and MFD output is strikingly
different.
g.region rast=elev_lid792_1m
# multiply elev_lid792_1m to get millimetres as vertical units
r.mapcalc "elev_lid792_1m.1mm = round(elev_lid792_1m * 1000.0)"
# SFD (D8) flow distribution
r.watershed elev=elev_lid792_1m.1mm threshold=3000 \
accum=elev_lid792_1m.1mm.acc.sfd drain=elev_lid792_1m.1mm.drain.sfd \
basin=elev_lid792_1m.1mm.basin.sfd stream=elev_lid792_1m.1mm.stream.sfd \
half_b=elev_lid792_1m.1mm.halfb.sfd vis=elev_lid792_1m.1mm.viz.sfd \
length_sl=elev_lid792_1m.1mm.ls.sfd slope_st=elev_lid792_1m.1mm.sls.sfd -s
# MFD flow distribution
r.watershed elev=elev_lid792_1m.1mm threshold=3000 \
accum=elev_lid792_1m.1mm.acc.mfd drain=elev_lid792_1m.1mm.drain.mfd \
basin=elev_lid792_1m.1mm.basin.mfd stream=elev_lid792_1m.1mm.stream.mfd \
half_b=elev_lid792_1m.1mm.halfb.mfd vis=elev_lid792_1m.1mm.viz.mfd \
length_sl=elev_lid792_1m.1mm.ls.mfd slope_st=elev_lid792_1m.1mm.sls.mfd \
convergence=5
Testing commands for the elevation map elevation.10m in the Spearfish
sample dataset spearfish60. Differences between SFD and MFD are mainly
in the northern half.
g.region rast=elevation.10m
# multiply elevation.10m to get centimetres as vertical units
r.mapcalc "elevation.10m.1cm = round(elevation.10m * 100.0)"
# SFD (D8) flow distribution
r.watershed elev=elevation.10m.1cm threshold=3000 \
accum=elevation.10m.1cm.acc.sfd drain=elevation.10m.1cm.drain.sfd \
basin=elevation.10m.1cm.basin.sfd stream=elevation.10m.1cm.stream.sfd \
half_b=elevation.10m.1cm.halfb.sfd vis=elevation.10m.1cm.viz.sfd \
length_sl=elevation.10m.1cm.ls.sfd slope_st=elevation.10m.1cm.sls.sfd -s
# MFD flow distribution
r.watershed elev=elevation.10m.1cm threshold=3000 \
accum=elevation.10m.1cm.acc.mfd drain=elevation.10m.1cm.drain.mfd \
basin=elevation.10m.1cm.basin.mfd stream=elevation.10m.1cm.stream.mfd \
half_b=elevation.10m.1cm.halfb.mfd vis=elevation.10m.1cm.viz.mfd \
length_sl=elevation.10m.1cm.ls.mfd slope_st=elevation.10m.1cm.sls.mfd \
convergence=5
Michael
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev
