[GRASSLIST:7488] Re: need help in concatenating commands

Massimiliano Cannata massimiliano.cannata at supsi.ch
Thu Jul 7 09:22:40 EDT 2005


Markus Neteler wrote:

>On Wed, Jul 06, 2005 at 04:38:52PM +0200, Massimiliano Cannata wrote:
>  
>
>>Hello,
>>I just need help in concatenating commands...
>>I build a command that, given a East and North coordinate scan a raster 
>>finding the nearest cell with value majon than a treshold...
>>    
>>
>>>r.nearest.coor input=accum100a threshold=500 radius=20 
>>>      
>>>
>>easting=692019.11712163 northing=119228.853597
>>this is the output of the command....
>>    
>>
>>>easting=690796.625000 northing=119081.002604
>>>      
>>>
>>I want now call the r.watershed.outlet.... concatenating this output...
>>    
>>
>>>r.water.outlet drainage=drainage100 basin=mybasin (easting=...... 
>>>      
>>>
>>northing=........)
>>
>>if I run:
>>    
>>
>>>r.nearest.coor input=accum100a threshold=500 radius=20 
>>>      
>>>
>>easting=692019.11712163 northing=119228.853597 | r.water.outlet 
>>drainage=drainage100 basin=a7642899097zzoo
>>
>>I have the message that r.water.outlet need easting and northing 
>>parameters......
>>
>>What to do?
>>    
>>
>
>Provide easting and northing to r.water.outlet :-)
>
>These are parameters, I don't think that r.water.outlet reads from
>stdin (pipe). So you have to slightly modify the last command.
>
>Markus
>
>  
>
I've found....

 r.water.outlet drainage=drainage100 basin=a7642899097zzooAAZ123 
`r.nearest.coor input=accum100a threshold=500 radius=20 
easting=694109.756531 northing=113416.282572`

that's it........
Ciao

-- 
-----------------------------------------------------
Massimiliano Cannata
Environmental & Geomatic Engineer

Institute of Earth Scences - SUPSI
Trevano, C.P. 72, CH-6952 Canobbio, SWITZERLAND
Tel: +41 (0)58 / 666 62 18  Fax: +41 (0)58 / 666 62 09
massimiliano.cannata at supsi.ch, www.supsi.ch
-------------------------------------------------------




More information about the grass-user mailing list