Skip to content

Commit 56f10a6

Browse files
authored
Update README.md
1 parent 40495bc commit 56f10a6

File tree

1 file changed

+20
-4
lines changed

1 file changed

+20
-4
lines changed

README.md

Lines changed: 20 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -652,8 +652,27 @@ Yes, you can provide any shapefiles you like. Here is an example: [FAQ_example_c
652652
We have provided some helpful Australian regions from /g/data/ia39, but the functions are flexible to take custom regions. [See more about the provided shapefiles here](https://github.com/aus-ref-clim-data-nci/shapefiles/).
653653
You will need to define [regionmask regions](https://regionmask.readthedocs.io/en/stable/notebooks/mask_3D.html) with unique abbreviations and names
654654

655+
You may have region data in other formats. [area_statistics_example_basin_gpkg.ipynb](https://github.com/AusClimateService/plotting_maps/blob/main/example_notebooks/area_statistics_example_basin_gpkg.ipynb) is an example using custom regions defined by a GeoPackage (GPKG).
656+
657+
</details>
658+
659+
### Can I plot Antarctica or other non-Australian areas of the world?
660+
<details>
661+
<summary> Expand </summary>
662+
663+
Yes, although `acs_plotting_maps` is designed to plot Australian hazard data, the functions are flexible to plot data for any area of the world.
664+
665+
For example, to plot Antarctica, you must adjust `xlim`, `ylim` and `projection`. In [FAQ_example_antarctica.ipynb](https://github.com/AusClimateService/plotting_maps/blob/main/example_notebooks/FAQ_example_antarctica.ipynb), we use `projection=ccrs.SouthPolarStereo()` for a polar projection as used by the Bureau of Meteorology for Southern Hemisphere maps. Limit the longitude and latitude with `xlim=(-180, 180)` and `ylim=(-90, -60)`. To plot the outline of the Antarctic continent (and other coastlines), set `coastlines = True`.
666+
667+
![Antarctica_sst_climatology](https://github.com/user-attachments/assets/bbdbdf45-173f-4dee-9441-c842b0723dd7)
668+
669+
In a similar way you can Europe by setting `coastlines=True`, `xlim=(-15, 45)`, `ylim=(30, 70)`, and `projection=ccrs.AlbersEqualArea(15, 50)`. See [FAQ_example_antarctica.ipynb](https://github.com/AusClimateService/plotting_maps/blob/main/example_notebooks/FAQ_example_antarctica.ipynb) for the full code to recreate this plot.
670+
671+
![Europe-sst_climatology](https://github.com/user-attachments/assets/ae0e52dd-240f-4438-9dcb-26951ed51ea6)
672+
655673
</details>
656674

675+
657676
### Can I use any regions for the acs_regional_stats statistics function?
658677
<details>
659678
<summary> Expand </summary>
@@ -774,7 +793,7 @@ b) When you use acs_regional_stats to reduce the data with ```dims = ("lat", "lo
774793
da_summary = acs_regional_stats(ds=ds, var=var, mask = mask, dims = ("lat", "lon", "model",), how = ["mean"],)
775794
```
776795

777-
https://github.com/AusClimateService/plotting_maps/blob/main/example_notebooks/FAQ_example_timeseries_stats_for_ensemble_region.ipynb
796+
[FAQ_example_timeseries_stats_for_ensemble_region.ipynb](https://github.com/AusClimateService/plotting_maps/blob/main/example_notebooks/FAQ_example_timeseries_stats_for_ensemble_region.ipynb) shows example of calculating regional means over multidimensional data.
778797

779798
</details>
780799

@@ -798,17 +817,14 @@ A range of teams are actively using this code. Take care to maintain backward co
798817
<summary> Expand </summary>
799818

800819
**Figures to make:**
801-
- SLR observations with station data. For the Climate Hazards report, recreate the sea level observations of gridded ocean data and station data in the same plot.
802820
- Lightning plot. For the Climate Hazards report, recreate the lightning observations plot using the plot_acs_hazards function so that it is in the consistent format.
803821

804822
**Documentation:**
805-
- Clear tutorials for how to use different functionalities
806823
- “release” this version of the code as v1.0.0. It is in wide use and care is being taken to keep backward compatibility while adding features. Eg see https://github.com/AusClimateService/plotting_maps/releases/new
807824

808825
**Improve plotting function and axillaries:**
809826
- Improve the aesthetics and proportions of plotting, especially with dataset/date_range/baseline annotations. Design aesthetics were focused on vertical orientations for 4-panel plots without these annotations for a particular report.
810827
- Improve the aesthetics of plotting select_area. Eg remove boundaries of neighbouring regions (if desired)
811-
- Try plotting Antarctica and see what settings work. Especially xlim, ylim, projection. This code should work for any region of the world, with the right settings.
812828
- Forest mask for forested areas. For example, FFDI is not useful in places where there is not connected vegetation/fuel. This is probably particularly for arid desert areas of central Australia. Changes in climate and land use may cause changes over time.
813829
- Improve colormap for fire climate classes. This colour scheme is not completely colourblind-friendly. Perhaps modify the colours to increase the contrast.
814830
- Enable rotating all labels and tick labels so that they are horizontal (easier to read). We may need to reduce the labels to every second increment. Eg for temperature.

0 commit comments

Comments
 (0)