Thank you Moritz,
i'll try to follow your suggestions
Il 14/06/2016 17:50, Moritz Lennert ha scritto:
> On 14/06/16 09:45, roy roy wrote:
>> I'm reading this wiki page:
>>
>> https://grasswiki.osgeo.org/wiki/Vector_Database_Management#The_concept_of_layers
>>
>>
>>
>> about the GRASS vector mod
Cool!
Thanks my friend!
A slightly modified version did the trick! My previous attempt included a
bunch of nested for loops that took forever :-/
map=`g.list rast pat=MOD13Q1*_NDVI`
for m in ${map}; do
i=`echo $m | cut -c 1-23`
echo $m $i
r.mapcalc --o \
expression="${m} = if("${i
Hi Vero
If I got you right, may be this could be a better approach,
#Assuming you are using MOD13Q1 product, and for example maps from
MOD13Q1.A2015353.h18v05.006.2016007180711.hdf;
map=`g.list rast pat=MOD13*_EVI`
for tile in ${map}; do
## Put the common part (In this case - MOD13Q1.A2015353
Hello everybody,
Is it possible to use r.mapcalc with `g.list` as part of the
expression??? Something like:
tiles=(h11v11 h12v11 h12v12)
for tile in ${tiles[*]}
do
r.mapcalc --o expression="`g.list rast pat=MOD13*$tile_EVI` =
if(`g.list rast pat=MOD13*$tile_pixel_reliability` == 0, `g.list rast
p