[GRASS-dev] r.stream.* to trunk - what about r.stream.basins?

Margherita Di Leo dileomargherita at gmail.com
Thu Mar 27 06:38:30 PDT 2014


Hi,


On Thu, Mar 27, 2014 at 2:20 PM, Moritz Lennert <
mlennert at club.worldonline.be> wrote:

> On 27/03/14 14:08, Newcomb, Doug wrote:
>
>>
>>
>>
>>     1) to keep r.stream.basins as an add-on, demanding the delineation of
>>     multiple subbasins from coordinates to a user´s script looping
>>     r.water.outlet;
>>     2) to include r.stream.basins and keep also r.water.outlet;
>>     3) to include r.stream.basins in the core and remove r.water.outlet as
>>     obsolete.
>>     4) ?
>>
>>     regards from Vienna
>>
>>     Madi, Helmut
>>
>>
>> If r.stream.basins can replicate the functionality of r.water.outlet, I
>> would prefer 3  . Some argument might be made for a conservative
>> approach for option 2 for version 7.0 and drop r.water.outlet in version
>> 7.1, but 7.0 is in development status. Why not make the change now?
>>
>
> As long as results are identical, +1. If not, this should be explored
> further.
>
>
I reproduced the example given in the manual page (NC dataset) of
r.water.outlet using r.stream.basin:

g.region rast=elev_lid792_1m -p

r.watershed elev_lid792_1m thresh=5000 accum=accum_5K
drain=draindir_5K basin=basin_5K
r.water.outlet input=draindir_5K output=basin_A30
coord=638740.423248,220271.519225

r.stream.basins dirs=draindir_5K at testing
coors=638740.423248,220271.519225 basins=basin_rstreambasins

r.mapcalc expression="diff=basin_rstreambasins-basin_A30"

r.info diff
 +----------------------------------------------------------------------------+
 | Layer:    diff                           Date: Thu Mar 27 14:34:25 2014    |
 | Mapset:   testing                        Login of Creator: madi            |
 | Location: nc_spm_08                                                        |
 | DataBase: /home/madi/grassdata                                             |
 | Title:     ( diff )                                                        |
 | Timestamp: none                                                            |
 |----------------------------------------------------------------------------|
 |                                                                            |
 |   Type of Map:  raster               Number of Categories: 0               |
 |   Data Type:    CELL                                                       |
 |   Rows:         750                                                        |
 |   Columns:      700                                                        |
 |   Total Cells:  525000                                                     |
 |        Projection: Lambert Conformal Conic                                 |
 |            N:     220750    S:     220000   Res:     1                     |
 |            E:     639000    W:     638300   Res:     1                     |
 |   Range of data:    min = 0  max = 0                                       |
 |                                                                            |
 |   Data Description:                                                        |
 |    generated by r.mapcalc                                                  |
 |                                                                            |
 |   Comments:                                                                |
 |    basin_rstreambasins - basin_A30                                         |
 |                                                                            |
 +----------------------------------------------------------------------------+

Best regards,
madi







-- 
Best regards,

Dr. Margherita DI LEO
Scientific / technical project officer

European Commission - DG JRC
Institute for Environment and Sustainability (IES)
Via Fermi, 2749
I-21027 Ispra (VA) - Italy - TP 261

Tel. +39 0332 78 3600
margherita.di-leo at jrc.ec.europa.eu

Disclaimer: The views expressed are purely those of the writer and may not
in any circumstance be regarded as stating an official position of the
European Commission.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20140327/72413437/attachment.html>


More information about the grass-dev mailing list