New notebook for loading data from Microsoft Planetary Computer

robbibt commented 1 year ago

Proposed changes

This PR adds a new notebook for loading data from Microsoft Planetary Computer and combining it with DEA data, including:

How to use pystac_client to search for ESA WorldCover data for a study area and time period
Load ESA WorldCover from Planetary Computer into the notebook using odc.stac
Compare ESA WorldCover outputs against DEA's Land Cover product
Load a time series of USGS Collection 2 Landsat surface temperature from Planetary Computer
Add this data to a time series of DEA Landsat ARD so they can be analysed together

Closes issues (optional)

Closes Issue #1088

Checklist (replace `[ ]` with `[x]` to check off)

[x] Notebook created using the DEA-notebooks template
[x] Remove any unused Python packages from Load packages
[x] Remove any unused/empty code cells
[x] Remove any guidance cells (e.g. General advice)
[x] Ensure that all code cells follow the PEP8 standard for code. The jupyterlab_code_formatter tool can be used to format code cells to a consistent style: select each code cell, then click Edit and then one of the Apply X Formatter options (YAPF or Black are recommended).
[x] Include relevant tags in the final notebook cell (refer to the DEA Tags Index, and re-use tags if possible)
[x] Clear all outputs, run notebook from start to finish, and save the notebook in the state where all cells have been sequentially evaluated
[ ] Test notebook on both the NCI and DEA Sandbox (flag if not working as part of PR and ask for help to solve if needed)
[ ] If applicable, update the Notebook currently compatible with the NCI|DEA Sandbox environment only line below the notebook title to reflect the environments the notebook is compatible with

review-notebook-app[bot] commented 1 year ago

Check out this pull request on

See visual diffs & provide feedback on Jupyter Notebooks.

Powered by ReviewNB

review-notebook-app[bot] commented 1 year ago

View / edit / reply to this conversation on ReviewNB

geoscience-aman commented on 2023-07-04T05:21:21Z ----------------------------------------------------------------

Markdown for a warning box maybe?

<div class="alert alert-block alert-warning">

<b>Warning:</b> Keep these extents as small as possible for reasonable loading times!

</div>

_robbibt commented on 2023-07-04T05:29:20Z_ ----------------------------------------------------------------

Great idea, will try it out

_robbibt commented on 2023-07-04T05:41:20Z_ ----------------------------------------------------------------

Will do a teal "info" box instead, feels friendlier!

review-notebook-app[bot] commented 1 year ago

View / edit / reply to this conversation on ReviewNB

geoscience-aman commented on 2023-07-04T05:21:22Z ----------------------------------------------------------------

Is it worth adding a note to say that not all datasets in the catalogue are currently available in the Planetary Computer API?

E.g. https://planetarycomputer.microsoft.com/dataset/storage/ghe

review-notebook-app[bot] commented 1 year ago

View / edit / reply to this conversation on ReviewNB

geoscience-aman commented on 2023-07-04T05:21:23Z ----------------------------------------------------------------

Line #9.    axes[1].set_title("DEA Land Cover\nLevel 3 class 124 ('Natural Aquatic Vegetation')");

Maybe the following code under this line to add some space between the legend label of the left plot and the ylabel of the right plot?

plt.subplots_adjust(wspace=0.35)

_robbibt commented on 2023-07-04T05:29:46Z_ ----------------------------------------------------------------

Neat, will do

_robbibt commented on 2023-07-04T05:34:30Z_ ----------------------------------------------------------------

Actually, I might switch to a more intuitive cmap and strip out the legends completely, they don't really do much here

review-notebook-app[bot] commented 1 year ago

View / edit / reply to this conversation on ReviewNB

geoscience-aman commented on 2023-07-04T05:21:24Z ----------------------------------------------------------------

Line #4.    # Apply USGS scaling factors to obtain temperature in Kelvins

Is it worth saying 'Collection 2 USGS scaling ...' in case someone accidentally applies the same scaling to Landsat Collection 1?

_robbibt commented on 2023-07-04T05:30:10Z_ ----------------------------------------------------------------

Yes, absolulety - will update

review-notebook-app[bot] commented 1 year ago

View / edit / reply to this conversation on ReviewNB

geoscience-aman commented on 2023-07-04T05:21:25Z ----------------------------------------------------------------

Line #9.    axes[1].set_title("USGS Landsat surface temperature");

Same as before under this line too to prevent overlapping text maybe?

plt.subplots_adjust(wspace=0.35)

robbibt commented 1 year ago

Great idea, will try it out

View entire conversation on ReviewNB

robbibt commented 1 year ago

Neat, will do

View entire conversation on ReviewNB

robbibt commented 1 year ago

Yes, absolulety - will update

View entire conversation on ReviewNB

robbibt commented 1 year ago

Actually, I might switch to a more intuitive cmap and strip out the legends completely, they don't really do much here

View entire conversation on ReviewNB

robbibt commented 1 year ago

Will do a teal "info" box instead, feels friendlier!

View entire conversation on ReviewNB

GeoscienceAustralia / dea-notebooks