-
Notifications
You must be signed in to change notification settings - Fork 2
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
review time series part 2 #11
base: main
Are you sure you want to change the base?
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I added something about temporal topology and relations. I prefer keeping the example removed, it's just a different way of trying to catch and maintain reader's attention
# Import the GRASS GIS packages we need | ||
import grass.script as gs | ||
import grass.jupyter as gj | ||
|
||
path_to_project = "eu_laea/italy_LST_daily" |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
path_to_project = "eu_laea/italy_LST_daily" | |
path_to_project = "italy_eu_laea/italy_LST_daily" |
@@ -61,9 +61,9 @@ scheme. In this way, we have: | |||
To run this tutorial locally or in Google Colab, you should install | |||
GRASS GIS 8.4+, and download the | |||
[daily MODIS LST project](https://zenodo.org/records/3564515). |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
[daily MODIS LST project](https://zenodo.org/records/3564515). | |
[daily MODIS LST project](https://zenodo.org/doi/10.5281/zenodo.3564514). |
tmin = gs.list_grouped(type="raster", pattern="lst_minimum_*")["italy_LST_daily"] | ||
tmax = gs.list_grouped(type="raster", pattern="lst_maximum_*")["italy_LST_daily"] | ||
tavg = gs.list_grouped(type="raster", pattern="lst_average_*")["italy_LST_daily"] |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I needed numbers, that's why I used the ??
. Does it do the right selection with the *?
:::{.callout-caution title="Question" collapse="true"} | ||
How would you obtain the average spring warming over the whole time series? | ||
|
||
```{python} | ||
# Average spring warming | ||
gs.run_command("t.rast.series", | ||
input="annual_spring_warming", | ||
output="avg_spring_warming", | ||
method="average") | ||
|
||
# Display raster map with interactive class | ||
spw_map = gj.InteractiveMap(width = 500, use_region=True, tiles="CartoDark") | ||
spw_map.add_raster("avg_spring_warming") | ||
spw_map.add_layer_control(position = "bottomright") | ||
spw_map.show() | ||
``` | ||
::: |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I would not remove this example. The question is just a trigger for the next topic. It's a different way of introducing a topic. And in any case, it is collapsed so it is up the the reader to unfold it.
::: {.callout-note title="Temporal relations"} | ||
Something about temporal relations... |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
::: {.callout-note title="Temporal relations"} | |
Something about temporal relations... | |
::: {.callout-note title="Temporal topology"} | |
Temporal topology refers to the temporal relations among time intervals or instances in different time series. They are mostly useful when sampling one STRDS with another (see [t.sample](https://grass.osgeo.org/grass-stable/manuals/t.sample.html) manual). Some examples of temporal relations are *start, equal, during, contain, overlap,* etc. See below a graphic representation of relations and sampling taken from Gebbert and Pebesma (2014). Note that A1 starts in S1, B3 is during S3, B4 equals S4, B1 overlaps S1. | |
![temporal_sampling_and_relations](https://github.com/user-attachments/assets/a0ffbf4b-2e0e-4a39-a392-8a70df17887b) | |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Dunno if attaching img like this really work...
I completely removed one of the examples within the callout questions, it was placed before the tool was introduced, which seemed problematic. It could be integrated later if needed.
I replaced the knitr engine for jupyter in the heading, my understanding is knitr is for R.