Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-12 Thread Markus Neteler
On Mon, Dec 10, 2012 at 1:17 PM, Glynn Clements gl...@gclements.plus.com There isn't currently an efficient way to join slices either using commands or the public C API. I wonder if the segment lib could be (ab)used for this? Markus ___ grass-user

Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-12 Thread Markus Metz
On Wed, Dec 12, 2012 at 7:37 PM, Markus Neteler nete...@osgeo.org wrote: On Mon, Dec 10, 2012 at 1:17 PM, Glynn Clements gl...@gclements.plus.com There isn't currently an efficient way to join slices either using commands or the public C API. I wonder if the segment lib could be (ab)used for

Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-11 Thread Andranik Hayrapetyan
Ok, thank you. I will try it and inform about results here. On Mon, Dec 10, 2012 at 4:17 PM, Glynn Clements gl...@gclements.plus.comwrote: Andranik Hayrapetyan wrote: r.mapcalc? I don't know how to overlay rasters with r.mapcalc. can you, please be a bit detailed?

Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-10 Thread Andranik Hayrapetyan
So, here are two versions of script, of parallel GRASS jobs with *split map into spatial chunks (possibly with overlap to gain smooth results)* approach. First script ( script1.bash ) is splitting the region in specified number of smaller regions, do the calculation ( calculating NDVI for example

Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-10 Thread Glynn Clements
Andranik Hayrapetyan wrote: r.mapcalc? I don't know how to overlay rasters with r.mapcalc. can you, please be a bit detailed? r.mapcalc 'outmap = if(isnull(map1),map2,map1)' But this is just as inefficient as r.patch. There isn't currently an efficient way to join slices

Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-09 Thread Andranik Hayrapetyan
r.mapcalc? I don't know how to overlay rasters with r.mapcalc. can you, please be a bit detailed? On Thu, Dec 6, 2012 at 4:53 PM, Helmut Kudrnovsky hel...@web.de wrote: So as I understand the only way to concatenate chunks of compressed image is r.patch r.mapcalc?

Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-06 Thread Andranik Hayrapetyan
Thank you for your response, I found an alternative way to avoid r.patch. But it is only for special cases when you need to output image after calculation at once. I put all chunks into a group with i.group and after it export the group with r.out.gdal . It is much more faster then when I

Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-06 Thread Andranik Hayrapetyan
Hi again, In order not to confuse you, I'd like to inform that my last post was wrong, I have some mistakes in my script and came across wrong results. Actually making group of chunks and exporting it making multiband image. It is working slow, as it is filling nodata of each chunk with 0, and

Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-06 Thread Helmut Kudrnovsky
So as I understand the only way to concatenate chunks of compressed image is r.patch r.mapcalc? http://grass.osgeo.org/grass64/manuals/r.mapcalc.html - best regards Helmut -- View this message in context:

[GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-05 Thread Andranik Hayrapetyan
Hi all, I am trying to make a script which will split a raster map into chunks, do the calculation on every chunk separately ( in parallel way ), and gather together all the peaces. So I am cutting the region into peaces, doing the calculation, restoring region to default region and finally get

Re: [GRASS-user] g.copy several map layers into one layer without overwriting .....

2012-12-05 Thread Glynn Clements
Andranik Hayrapetyan wrote: I am trying to make a script which will split a raster map into chunks, do the calculation on every chunk separately ( in parallel way ), and gather together all the peaces. So I am cutting the region into peaces, doing the calculation, restoring region to