[GRASS-user] Maps in a defined region

Ondřej Pešek pesej.ondrek at gmail.com
Thu Oct 3 06:48:12 PDT 2019


čt 3. 10. 2019 v 15:36 odesílatel Ken Mankoff <mankoff at gmail.com> napsal:

>
> On 2019-10-03 at 15:25 +02, Ondřej Pešek <pesej.ondrek at gmail.com> wrote...
> > but r.out.gdal would export "something" from all the maps in the loop,
> > wouldn't it? I thought that if I would loop through all my rasters and
> > pass them to r.out.gdal, it would create an output for every of them,
> > just the ones outside the computational region will be full of no data
> > values. I wanted to avoid this export of empty rasters.
>
> In the loop you cant test each raster for some non-null values:
>
> for r in $(g.list...); do
>   eval $(r.univar -g ...)
>   if [[ ${sum} != 0 ]]; then...
>
>
Oh, right. That seems like the smoothest solution if there is no
specialized tool. Thanks to all of you.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20191003/705a19af/attachment.html>


More information about the grass-user mailing list