Criterion B Details

# Default ecosystem code for template development.
# This line is replaced by build_ecosystem_pages.py for each ecosystem.
ecosystem_code = 'F2.1.1'

Import Python modules.

import os
import yaml
from pathlib import Path
from lonboard import Map
from rle_python_gee.ecosystems import Ecosystems
from rle_python_gee.eoo import make_eoo
from rle_python_gee.aoo import make_aoo_grid

Load the country config file.

project_root = os.environ.get('PIXI_PROJECT_ROOT', str(Path('..').resolve()))
config_path = Path(project_root) / 'config' / 'country_config.yaml'
with open(config_path) as f:
    config = yaml.safe_load(f)

Load & Filter Ecosystem Data

Load data for all the ecosystems.

source = config['ecosystem_source']
ecosystems = Ecosystems.from_file(
    source['data'],
    ecosystem_column=source['ecosystem_code_column'],
    ecosystem_name_column=source['ecosystem_name_column'],
    functional_group_column=source['functional_group_column']
)

Filter by the F2.1.1 and check the number of features.

ecosystem = ecosystems.filter(ecosystem_code)
has_data = ecosystem.size() > 0
print(f'{ecosystem.size() = }')
if not has_data:
    from IPython.display import Markdown, display
    display(Markdown(
        f'**No spatial data found for {ecosystem_code}.** '
        f'Criterion B calculations are skipped.'
    ))

ecosystem.size() = 14

Extent of Occurrence (EOO) (subcriterion B1)

Extent of occurrence (EOO). The EOO of an ecosystem is the area (km2) of a minimum convex polygon – the smallest polygon in which no internal angle exceeds 180° that encompasses all known current spatial occurrences of the ecosystem type.

The minimum convex polygon (also known as a convex hull) must not exclude any areas, discontinuities or disjunctions, regardless of whether the ecosystem can occur in those areas or not. Regions such as oceans (for terrestrial ecosystems), land (for coastal or marine ecosystems), or areas outside the study area (such as in a different country) must remain included within the minimum convex polygon to ensure that this standardised method is comparable across ecosystem types. In addition, these features contribute to spreading risks across the distribution of the ecosystem by making different parts of its distribution more spatially independent.

Calculate EOO

Start by calculating the convex hull of the ecosystem’s distribution.

import geopandas as gpd

if has_data:
    ecosystem_geometry = ecosystem.geometry.union_all()
    gdf_ecosystem_polygons = gpd.GeoDataFrame(geometry=[ecosystem_geometry], crs=ecosystem.geometry.crs)
    hull = ecosystem_geometry.convex_hull
    gdf_hull = gpd.GeoDataFrame(geometry=[hull], crs=ecosystem.geometry.crs)

Display the ecosystem’s distribution and the convex hull.

from lonboard import Map, PolygonLayer
from rle_python_gee.viz import smart_map

if has_data:
    eoo_hull = make_eoo(ecosystem).compute()
    display(smart_map([eoo_hull, ecosystem]))

/home/runner/work/rle-tyler-20/rle-tyler-20/.pixi/envs/default/lib/python3.11/site-packages/lonboard/_geoarrow/ops/reproject.py:116: UserWarning: Input being reprojected to EPSG:4326 CRS.
Lonboard is only able to render data in EPSG:4326 projection.
  warnings.warn(
/home/runner/work/rle-tyler-20/rle-tyler-20/.pixi/envs/default/lib/python3.11/site-packages/lonboard/_geoarrow/ops/reproject.py:116: UserWarning: Input being reprojected to EPSG:4326 CRS.
Lonboard is only able to render data in EPSG:4326 projection.
  warnings.warn(

if has_data:
    hull_ea = gdf_hull.to_crs("ESRI:54034")
    eoo = hull_ea.geometry.iloc[0].area / 1e6
    print(f'EOO is {eoo:.1f} km2')

EOO is 198756.9 km2

Then calculate the area of the convex hull polygon.

Direct calculation of EOO

EOO can also be calculated directly using …

if has_data:
    ecosystem.eoo

Verify that the area returned by calling make_eoo(ecosystem).compute().area_km2 is the same as the area of the convex hull polygon.

if has_data:
    assert ecosystem.eoo == eoo

Area of Occupancy (AOO) (subcriterion B2)

The protocol for this adjustment includes the following steps:

Intersect AOO grid with the ecosystem’s distribution map.

Calculate extent of the ecosystem type in each grid cell (area) and sum these areas to obtain the total ecosystem area (total area).

Arrange grid cells in ascending order based on their area (smaller first). Calculate accumulated sum of area per cell (cumulative area).

Calculate cumulative proportion by dividing cumulative area by total area (cumulative proportion takes values between 0 and 1)

Calculate AOO by counting the number of cells with a cumulative proportion greater than 0.01 (i.e. exclude cells that in combination account for up to 1% of the total mapped extent of the ecosystem type).

AOO Calculation Details

Intersect AOO grid and ecosystem map

Intersect AOO grid with the ecosystem’s distribution map

from pathlib import Path
from rle_python_gee.aoo import make_aoo_grid_cached

if has_data:
    cache_path = Path(project_root) / '.cache' / 'aoo_grid.parquet'
    aoo_grid = make_aoo_grid_cached(ecosystems, cache_path=cache_path)
    aoo_grid_filtered = aoo_grid.filter_by_ecosystem(ecosystem_code)

Visualize variations in the AOO grid.

from matplotlib.colors import LinearSegmentedColormap
from lonboard.colormap import apply_continuous_cmap
from rle_python_gee.aoo import slugify_ecosystem_name

ecosystem_column = slugify_ecosystem_name(ecosystem_code)
if has_data:
    cmap = LinearSegmentedColormap.from_list("white_red", ["white", "red"])
    values = aoo_grid_filtered.grid_cells[ecosystem_column].values
    normalized = (values - values.min()) / (values.max() - values.min())
    colors = apply_continuous_cmap(normalized, cmap)
    display(smart_map([(aoo_grid_filtered, {"get_fill_color": colors}), ecosystem]))

/home/runner/work/rle-tyler-20/rle-tyler-20/.pixi/envs/default/lib/python3.11/site-packages/lonboard/_geoarrow/ops/reproject.py:116: UserWarning: Input being reprojected to EPSG:4326 CRS.
Lonboard is only able to render data in EPSG:4326 projection.
  warnings.warn(

Calculate grid cell area and total area

Calculate extent of the ecosystem type in each grid cell (area) and sum these areas to obtain the total ecosystem area (total area).

if has_data:
    keep = ['geometry', 'grid_col', 'grid_row', ecosystem_column]
    gdf = aoo_grid_filtered.grid_cells[keep]
    display(gdf)

	geometry	grid_col	grid_row	F2_1_1
0	POLYGON ((-77.16528 0.99486, -77.16528 1.08531...	-860	11	0.044107
1	POLYGON ((-77.16528 1.08531, -77.16528 1.17576...	-860	12	0.012848
2	POLYGON ((-77.07545 0.99486, -77.07545 1.08531...	-859	11	0.127875
3	POLYGON ((-77.07545 1.08531, -77.07545 1.17576...	-859	12	0.235943
4	POLYGON ((-75.09916 10.36476, -75.09916 10.456...	-837	114	0.135626
5	POLYGON ((-75.09916 10.45667, -75.09916 10.548...	-837	115	0.096738
6	POLYGON ((-75.09916 10.54861, -75.09916 10.640...	-837	116	0.007393
7	POLYGON ((-75.00933 10.36476, -75.00933 10.456...	-836	114	0.203752
8	POLYGON ((-75.00933 10.45667, -75.00933 10.548...	-836	115	0.558854
9	POLYGON ((-75.00933 10.54861, -75.00933 10.640...	-836	116	0.350296
10	POLYGON ((-74.91949 10.45667, -74.91949 10.548...	-835	115	0.006663
11	POLYGON ((-74.91949 10.54861, -74.91949 10.640...	-835	116	0.000125
12	POLYGON ((-73.75168 5.34332, -73.75168 5.43414...	-822	59	0.000006
13	POLYGON ((-73.75168 5.43414, -73.75168 5.52498...	-822	60	0.087679
14	POLYGON ((-73.66185 5.34332, -73.66185 5.43414...	-821	59	0.003016
15	POLYGON ((-73.66185 5.43414, -73.66185 5.52498...	-821	60	0.171675
16	POLYGON ((-72.9432 5.43414, -72.9432 5.52498, ...	-813	60	0.081912
17	POLYGON ((-72.9432 5.52498, -72.9432 5.61584, ...	-813	61	0.026972
18	POLYGON ((-72.85337 5.43414, -72.85337 5.52498...	-812	60	0.094895
19	POLYGON ((-72.85337 5.52498, -72.85337 5.61584...	-812	61	0.352469

The column F2_1_1 contains the (fractional) area of the ecosystem in each grid cell.

Sum up the areas of each grid cell to get the total area.

if has_data:
    total_area = gdf[ecosystem_column].sum()
    display(total_area)

np.float64(2.5988422769822472)

Calculate cumulative area

Arrange grid cells in ascending order based on their area (smaller first). Calculate accumulated sum of area per cell (cumulative area).

if has_data:
    gdf = gdf.sort_values(by=ecosystem_column)
    gdf["cumulative_area"] = gdf[ecosystem_column].cumsum()
    display(gdf)

	geometry	grid_col	grid_row	F2_1_1	cumulative_area
12	POLYGON ((-73.75168 5.34332, -73.75168 5.43414...	-822	59	0.000006	0.000006
11	POLYGON ((-74.91949 10.54861, -74.91949 10.640...	-835	116	0.000125	0.000131
14	POLYGON ((-73.66185 5.34332, -73.66185 5.43414...	-821	59	0.003016	0.003147
10	POLYGON ((-74.91949 10.45667, -74.91949 10.548...	-835	115	0.006663	0.009810
6	POLYGON ((-75.09916 10.54861, -75.09916 10.640...	-837	116	0.007393	0.017202
1	POLYGON ((-77.16528 1.08531, -77.16528 1.17576...	-860	12	0.012848	0.030050
17	POLYGON ((-72.9432 5.52498, -72.9432 5.61584, ...	-813	61	0.026972	0.057022
0	POLYGON ((-77.16528 0.99486, -77.16528 1.08531...	-860	11	0.044107	0.101129
16	POLYGON ((-72.9432 5.43414, -72.9432 5.52498, ...	-813	60	0.081912	0.183041
13	POLYGON ((-73.75168 5.43414, -73.75168 5.52498...	-822	60	0.087679	0.270720
18	POLYGON ((-72.85337 5.43414, -72.85337 5.52498...	-812	60	0.094895	0.365615
5	POLYGON ((-75.09916 10.45667, -75.09916 10.548...	-837	115	0.096738	0.462353
2	POLYGON ((-77.07545 0.99486, -77.07545 1.08531...	-859	11	0.127875	0.590228
4	POLYGON ((-75.09916 10.36476, -75.09916 10.456...	-837	114	0.135626	0.725854
15	POLYGON ((-73.66185 5.43414, -73.66185 5.52498...	-821	60	0.171675	0.897529
7	POLYGON ((-75.00933 10.36476, -75.00933 10.456...	-836	114	0.203752	1.101281
3	POLYGON ((-77.07545 1.08531, -77.07545 1.17576...	-859	12	0.235943	1.337224
9	POLYGON ((-75.00933 10.54861, -75.00933 10.640...	-836	116	0.350296	1.687520
19	POLYGON ((-72.85337 5.52498, -72.85337 5.61584...	-812	61	0.352469	2.039988
8	POLYGON ((-75.00933 10.45667, -75.00933 10.548...	-836	115	0.558854	2.598842

Calculate cumulative proportion

Calculate cumulative proportion by dividing cumulative area by total area (cumulative proportion takes values between 0 and 1)

if has_data:
    gdf["cumulative_proportion"] = gdf["cumulative_area"] / total_area
    display(gdf)

	geometry	grid_col	grid_row	F2_1_1	cumulative_area	cumulative_proportion
12	POLYGON ((-73.75168 5.34332, -73.75168 5.43414...	-822	59	0.000006	0.000006	0.000002
11	POLYGON ((-74.91949 10.54861, -74.91949 10.640...	-835	116	0.000125	0.000131	0.000050
14	POLYGON ((-73.66185 5.34332, -73.66185 5.43414...	-821	59	0.003016	0.003147	0.001211
10	POLYGON ((-74.91949 10.45667, -74.91949 10.548...	-835	115	0.006663	0.009810	0.003775
6	POLYGON ((-75.09916 10.54861, -75.09916 10.640...	-837	116	0.007393	0.017202	0.006619
1	POLYGON ((-77.16528 1.08531, -77.16528 1.17576...	-860	12	0.012848	0.030050	0.011563
17	POLYGON ((-72.9432 5.52498, -72.9432 5.61584, ...	-813	61	0.026972	0.057022	0.021941
0	POLYGON ((-77.16528 0.99486, -77.16528 1.08531...	-860	11	0.044107	0.101129	0.038913
16	POLYGON ((-72.9432 5.43414, -72.9432 5.52498, ...	-813	60	0.081912	0.183041	0.070432
13	POLYGON ((-73.75168 5.43414, -73.75168 5.52498...	-822	60	0.087679	0.270720	0.104169
18	POLYGON ((-72.85337 5.43414, -72.85337 5.52498...	-812	60	0.094895	0.365615	0.140684
5	POLYGON ((-75.09916 10.45667, -75.09916 10.548...	-837	115	0.096738	0.462353	0.177907
2	POLYGON ((-77.07545 0.99486, -77.07545 1.08531...	-859	11	0.127875	0.590228	0.227112
4	POLYGON ((-75.09916 10.36476, -75.09916 10.456...	-837	114	0.135626	0.725854	0.279299
15	POLYGON ((-73.66185 5.43414, -73.66185 5.52498...	-821	60	0.171675	0.897529	0.345357
7	POLYGON ((-75.00933 10.36476, -75.00933 10.456...	-836	114	0.203752	1.101281	0.423758
3	POLYGON ((-77.07545 1.08531, -77.07545 1.17576...	-859	12	0.235943	1.337224	0.514546
9	POLYGON ((-75.00933 10.54861, -75.00933 10.640...	-836	116	0.350296	1.687520	0.649335
19	POLYGON ((-72.85337 5.52498, -72.85337 5.61584...	-812	61	0.352469	2.039988	0.784960
8	POLYGON ((-75.00933 10.45667, -75.00933 10.548...	-836	115	0.558854	2.598842	1.000000

Count AOO cells

Calculate AOO by counting the number of cells with a cumulative proportion greater than 0.01 (i.e. exclude cells that in combination account for up to 1% of the total mapped extent of the ecosystem type).

if has_data:
    aoo = len(gdf[gdf["cumulative_proportion"] > 0.01])
    print(f'AOO is {aoo} cells')

AOO is 15 cells

AOO Calculation (direct call)

if has_data:
    aoo_count = ecosystem.aoo
    print(f'AOO: {aoo_count} grid cells')

AOO: 15 grid cells

if has_data:
    display(smart_map([aoo_grid_filtered, ecosystem]))

/home/runner/work/rle-tyler-20/rle-tyler-20/.pixi/envs/default/lib/python3.11/site-packages/lonboard/_geoarrow/ops/reproject.py:116: UserWarning: Input being reprojected to EPSG:4326 CRS.
Lonboard is only able to render data in EPSG:4326 projection.
  warnings.warn(

Criterion B Summary

Ecosystem Code	Ecosystem Name	EOO	AOO
F2.1.1	INDEF	198757 km²	15 cells