<div dir="ltr"><div>Hi Stefan</div><div><br></div><div>I'll test this weekend, would you mind providing us with a basic example to test? Thanks a lot!!</div><div><br></div><div>Vero<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El jue, 5 ene 2023 a las 11:01, Stefan Blumentrath (<<a href="mailto:Stefan.Blumentrath@gmx.de">Stefan.Blumentrath@gmx.de</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi again,<br>
<br>
Please see:<br>
<a href="https://github.com/OSGeo/grass/pull/2725" rel="noreferrer" target="_blank">https://github.com/OSGeo/grass/pull/2725</a><br>
<br>
It uses east, west, north and south in the where statement to select only maps that overlap with a given spatial extent (e.g. the computatinal region).<br>
I would be happy to get a review of the PR, before proceding to eventual changes of module code...<br>
<br>
<br>
Cheers<br>
Stefan<br>
<br>
<br>
<br>
<br>
Gesendet: Montag, 26. Dezember 2022 um 17:31 Uhr<br>
Von: "Sören Gebbert" <<a href="mailto:soerengebbert@gmail.com" target="_blank">soerengebbert@gmail.com</a>><br>
An: "Stefan Blumentrath" <<a href="mailto:Stefan.Blumentrath@gmx.de" target="_blank">Stefan.Blumentrath@gmx.de</a>><br>
Cc: <a href="mailto:grass-dev@lists.osgeo.org" target="_blank">grass-dev@lists.osgeo.org</a><br>
Betreff: Re: [GRASS-dev] Filter STDS by spatial index<br>
<br>
Hi,<br>
i am not sure if this helps: some t.rast.* commands should support a where statement. Maybe you can use east, west, north and south in the where statement to select only maps that are in a extent that is build in the where statement?<br>
<br>
Best regards<br>
Sören <br>
<br>
Am Mi., 14. Dez. 2022 um 15:24 Uhr schrieb Stefan Blumentrath <<a href="mailto:Stefan.Blumentrath@gmx.de" target="_blank">Stefan.Blumentrath@gmx.de</a>[mailto:<a href="mailto:Stefan.Blumentrath@gmx.de" target="_blank">Stefan.Blumentrath@gmx.de</a>]>:Hei,<br>
<br>
I do have a use case where I would like to filter maps in a STRDS by spatial extent before e.g. t.rast.univar processes all registered maps.<br>
In particular, I want to process only maps from a STRDS that intersect with the current computational region.<br>
<br>
The reason is that the maps registered in the STRDS are collections of single scenes of satellite images or collections of smaller mosaics that cover only parts of the total extent of the STRDS.<br>
<br>
Currently, that does not seem to be supported by TGIS, and I see two options to address that:<br>
a) patch the scenes together with a given granularity (e.g. daily mosaics), leave GRASS GIS functions untouched and hope that there is some coverage for reasonable region settings. or:<br>
b) add a spatial filter to the "get_registered_maps"[1] function in abstract_space_time_dataset.py, and then subsequently flags to relevant modules that make use modified function...<br>
<br>
For b) I would think that the most efficient way to handle it would be to add a where clause like it was done for the semantic labels? In that case, it may make sense to add an index for north, south, east, and west too, no? Maybe also top and bottom...<br>
<br>
My main question, however is if my case is just a corner case, or does some sort of spatial filteing possibility make sense in general / more broadly...<br>
<br>
Any feedback is much appreciated, esp. before I evtl. start working on b).<br>
<br>
Kind regards,<br>
Stefan<br>
<br>
1) <a href="https://github.com/OSGeo/grass/blob/main/python/grass/temporal/abstract_space_time_dataset.py#L1604[https://github.com/OSGeo/grass/blob/main/python/grass/temporal/abstract_space_time_dataset.py%23L1604]" rel="noreferrer" target="_blank">https://github.com/OSGeo/grass/blob/main/python/grass/temporal/abstract_space_time_dataset.py#L1604[https://github.com/OSGeo/grass/blob/main/python/grass/temporal/abstract_space_time_dataset.py#L1604]</a><br>
<br>
_______________________________________________<br>
grass-dev mailing list<br>
<a href="mailto:grass-dev@lists.osgeo.org" target="_blank">grass-dev@lists.osgeo.org</a>[mailto:<a href="mailto:grass-dev@lists.osgeo.org" target="_blank">grass-dev@lists.osgeo.org</a>]<br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-dev</a><br>
_______________________________________________<br>
grass-dev mailing list<br>
<a href="mailto:grass-dev@lists.osgeo.org" target="_blank">grass-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-dev</a><br>
</blockquote></div>