add shapefile to dictionary: plantations and NCRA_Marine_region

AusClimateService · Aug 22, 2024 · e79c94d · e79c94d
1 parent fd72856
commit e79c94d
Showing 1 changed file with 8 additions and 3 deletions.
diff --git a/acs_area_statistics.py b/acs_area_statistics.py
@@ -5,7 +5,6 @@
 
 This method has used guidance from
 [https://github.com/aus-ref-clim-data-nci/shapefiles/blob/master/python_tutorial.ipynb]"""
-
 import xarray as xr
 import geopandas as gpd
 import regionmask
@@ -17,27 +16,34 @@
                   "aus_states_territories",
                   "australia", 
                   "broadacre_regions", 
+                  "NCRA_Marine_region",
                   "ncra_regions", 
                   "NCRA_regions_coastal_waters_GDA94", 
-                  "nrm_regions", ]
+                  "nrm_regions",
+                  "plantations"]
+
 
 name_dict = {"aus_local_gov":"LGA_NAME22", 
              "aus_states_territories":"STE_NAME21",
              "australia": "AUS_NAME21",
              "broadacre_regions": "name",
+             "NCRA_Marine_region":"Label",
              "ncra_regions": "regionname", 
              "NCRA_regions_coastal_waters_GDA94": "regionname",
              "nrm_regions":"SubClusNm",
+             "plantations":"CodeName",
             }
 
 
 abbr_dict = {"aus_local_gov":"LGA_CODE22", 
              "aus_states_territories":"ABBREV",
              "australia": "AUS_CODE21",
              "broadacre_regions": "aagis",
+             "NCRA_Marine_region":"RegionID",
              "ncra_regions": "short_name", 
              "NCRA_regions_coastal_waters_GDA94": "NCRA",
              "nrm_regions": "SubClusAb",
+             "plantations":"REGCODE",
             }
 
 def get_regions(shapefiles):
@@ -67,7 +73,6 @@ def get_regions(shapefiles):
     return regionmask.from_geopandas(gdf, names="NAME", abbrevs="abbrevs", name= "-".join(shapefiles), overlap=True) 
 
 
-
 # read in the shapefile with regions you will use Australia and NCRA regions
 # from acs_plotting_maps import regions_dict
 # regions = get_regions(["ncra_regions", "australia"])