<p>Dear all,</p>
<p>I have a large raster project (~60,000 x 48,000 pixels)<br>and I extracted the drainage from DEM for ~1,000 subbasin.<br>Now I would like to join all the pieces of drainage, but<br>whem I tryed :<br> g.region  rast=`g.mlist pat=ANA_COD5_drenbin_* sep=,`<br>
 r.series input=`g.mlist pat=ANA_COD5_drenbin_* sep=,` output=ANA_COD5_drenbin_join method=maximum --o</p>
<p>I get error because the (I suppose) the shell not support<br>my large number of raster names on rast= and input=<br>at g.region and r.series, respectivelly.</p>
<p>As my base map is very large, if I try something like</p>
<p>r.mapcalc &quot;ANA_COD5_drenbin_joinAUX=if(AnyMAP&gt;0, AnyMAP, ANA_COD5_drenbin_join)&quot;<br>g.rename rast=ANA_COD5_drenbin_joinAUX,ANA_COD5_drenbin_join</p>
<p>it is very time expensive, because of the size of my base/large<br>map. Is there a way of I update only part of a map, using<br>g.region or r.region for each subbasing, update my<br>base map without loss other information already stored<br>
on base map?</p>
<p>A little bit late, but still wishing Merry Xtmas.</p>
<p>miltinho<br></p>