# Default ecosystem code for template development.
# This line is replaced by build_ecosystem_pages.py for each ecosystem.
ecosystem_code = 'Zona Pantanosa Andina'Criterion B Details
Import Python modules.
import os
import yaml
from pathlib import Path
from lonboard import Map
from rle.core import Ecosystems, criterion_b_status, rle_category
from rle.core.eoo import make_eoo
from rle.core.aoo import make_aoo_gridLoad 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)
# Ecosystem index (the COG pixel value for this ecosystem), looked up by code
# from the canonical index table. Absent if the index has not been built yet.
import csv
ecosystem_index = None
index_csv = Path(project_root) / 'config' / 'ecosystems' / 'index.csv'
if index_csv.exists():
with open(index_csv) as f:
for row in csv.DictReader(f):
if row['code'] == ecosystem_code:
ecosystem_index = int(row['index'])
breakLoad & Filter Ecosystem Data
Load data for all the ecosystems.
import sys
# Shared config helpers live in scripts/; make them importable from here.
sys.path.insert(0, str(Path(project_root) / 'scripts'))
from _config import ensure_vector_source
source = config['ecosystem_source']
# ecosystem_code_column is optional: fall back to the name column.
ecosystem_column = source.get('ecosystem_code_column') or source.get('ecosystem_name_column')
# Prefer the ecosystem-sorted `optimized_data` copy when configured, so that
# filtering to one ecosystem uses parquet predicate pushdown instead of loading
# the whole national map into memory. Falls back to `data`.
ecosystems = Ecosystems.from_file(
ensure_vector_source(source.get('optimized_data') or source['data']),
ecosystem_column=ecosystem_column,
ecosystem_name_column=source.get('ecosystem_name_column'),
functional_group_column=source.get('functional_group_column')
)Filter by the Zona Pantanosa Andina 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() = 43Extent 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.core.viz import smart_map
if has_data:
eoo_hull = make_eoo(ecosystem).compute()
display(smart_map([eoo_hull, ecosystem]))/home/runner/work/rle-tyler-colombia/rle-tyler-colombia/.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-colombia/rle-tyler-colombia/.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 114484.9 km2Then calculate the area of the convex hull polygon.
Direct calculation of EOO
EOO can also be calculated directly using …
if has_data:
ecosystem.eooVerify 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 == eooArea 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 proportionby dividingcumulative areabytotal area(cumulative proportion takes values between 0 and 1)- Calculate AOO by counting the number of cells with a
cumulative proportiongreater 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.core.aoo import make_aoo_grid_cached
if has_data:
# Prefer a prebuilt grid cache (e.g. a gs:// URI) so the national AOO grid
# is not recomputed during CI renders — computing it from the full national
# ecosystem map peaks at many GB of RAM. Fall back to a local cache when no
# prebuilt cache is configured. Build one with `pixi run build-caches`.
cache_path = (
source.get('aoo_grid_cache_url')
or (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.core.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-colombia/rle-tyler-colombia/.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 | Zona_Pantanosa_Andina | |
|---|---|---|---|---|
| 0 | POLYGON ((-77.70427 0.90441, -77.70427 0.99486... | -866 | 10 | 0.001515 |
| 1 | POLYGON ((-77.61444 0.63307, -77.61444 0.72351... | -865 | 7 | 0.001161 |
| 2 | POLYGON ((-77.61444 0.72351, -77.61444 0.81396... | -865 | 8 | 0.010586 |
| 3 | POLYGON ((-77.61444 0.90441, -77.61444 0.99486... | -865 | 10 | 0.010972 |
| 4 | POLYGON ((-75.45848 6.6162, -75.45848 6.70724,... | -841 | 73 | 0.002899 |
| 5 | POLYGON ((-75.36865 4.70788, -75.36865 4.79862... | -840 | 52 | 0.013525 |
| 6 | POLYGON ((-75.27882 4.70788, -75.27882 4.79862... | -839 | 52 | 0.024420 |
| 7 | POLYGON ((-75.27882 4.79862, -75.27882 4.88937... | -839 | 53 | 0.006715 |
| 8 | POLYGON ((-74.29067 4.70788, -74.29067 4.79862... | -828 | 52 | 0.003951 |
| 9 | POLYGON ((-74.20084 4.16367, -74.20084 4.25434... | -827 | 46 | 0.005959 |
| 10 | POLYGON ((-74.20084 4.52643, -74.20084 4.61715... | -827 | 50 | 0.006292 |
| 11 | POLYGON ((-74.20084 4.61715, -74.20084 4.70788... | -827 | 51 | 0.035877 |
| 12 | POLYGON ((-74.20084 4.70788, -74.20084 4.79862... | -827 | 52 | 0.016856 |
| 13 | POLYGON ((-74.11101 4.52643, -74.11101 4.61715... | -826 | 50 | 0.006286 |
| 14 | POLYGON ((-74.11101 4.61715, -74.11101 4.70788... | -826 | 51 | 0.008861 |
| 15 | POLYGON ((-74.11101 4.70788, -74.11101 4.79862... | -826 | 52 | 0.053472 |
| 16 | POLYGON ((-74.02118 4.61715, -74.02118 4.70788... | -825 | 51 | 0.003801 |
| 17 | POLYGON ((-74.02118 4.70788, -74.02118 4.79862... | -825 | 52 | 0.020832 |
| 18 | POLYGON ((-73.93135 4.88937, -73.93135 4.98013... | -824 | 54 | 0.007423 |
| 19 | POLYGON ((-73.93135 4.98013, -73.93135 5.07091... | -824 | 55 | 0.006425 |
| 20 | POLYGON ((-73.75168 5.2525, -73.75168 5.34332,... | -822 | 58 | 0.023376 |
| 21 | POLYGON ((-73.75168 5.34332, -73.75168 5.43414... | -822 | 59 | 0.000452 |
| 22 | POLYGON ((-73.75168 5.43414, -73.75168 5.52498... | -822 | 60 | 0.082528 |
| 23 | POLYGON ((-73.66185 5.34332, -73.66185 5.43414... | -821 | 59 | 0.011092 |
| 24 | POLYGON ((-73.66185 5.43414, -73.66185 5.52498... | -821 | 60 | 0.122049 |
| 25 | POLYGON ((-73.57202 4.98013, -73.57202 5.07091... | -820 | 55 | 0.005972 |
| 26 | POLYGON ((-73.12286 5.70671, -73.12286 5.79759... | -815 | 63 | 0.004846 |
The column Zona_Pantanosa_Andina 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(0.49814296350714343)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 | Zona_Pantanosa_Andina | cumulative_area | |
|---|---|---|---|---|---|
| 21 | POLYGON ((-73.75168 5.34332, -73.75168 5.43414... | -822 | 59 | 0.000452 | 0.000452 |
| 1 | POLYGON ((-77.61444 0.63307, -77.61444 0.72351... | -865 | 7 | 0.001161 | 0.001613 |
| 0 | POLYGON ((-77.70427 0.90441, -77.70427 0.99486... | -866 | 10 | 0.001515 | 0.003128 |
| 4 | POLYGON ((-75.45848 6.6162, -75.45848 6.70724,... | -841 | 73 | 0.002899 | 0.006027 |
| 16 | POLYGON ((-74.02118 4.61715, -74.02118 4.70788... | -825 | 51 | 0.003801 | 0.009827 |
| 8 | POLYGON ((-74.29067 4.70788, -74.29067 4.79862... | -828 | 52 | 0.003951 | 0.013779 |
| 26 | POLYGON ((-73.12286 5.70671, -73.12286 5.79759... | -815 | 63 | 0.004846 | 0.018625 |
| 9 | POLYGON ((-74.20084 4.16367, -74.20084 4.25434... | -827 | 46 | 0.005959 | 0.024585 |
| 25 | POLYGON ((-73.57202 4.98013, -73.57202 5.07091... | -820 | 55 | 0.005972 | 0.030556 |
| 13 | POLYGON ((-74.11101 4.52643, -74.11101 4.61715... | -826 | 50 | 0.006286 | 0.036843 |
| 10 | POLYGON ((-74.20084 4.52643, -74.20084 4.61715... | -827 | 50 | 0.006292 | 0.043135 |
| 19 | POLYGON ((-73.93135 4.98013, -73.93135 5.07091... | -824 | 55 | 0.006425 | 0.049559 |
| 7 | POLYGON ((-75.27882 4.79862, -75.27882 4.88937... | -839 | 53 | 0.006715 | 0.056274 |
| 18 | POLYGON ((-73.93135 4.88937, -73.93135 4.98013... | -824 | 54 | 0.007423 | 0.063697 |
| 14 | POLYGON ((-74.11101 4.61715, -74.11101 4.70788... | -826 | 51 | 0.008861 | 0.072558 |
| 2 | POLYGON ((-77.61444 0.72351, -77.61444 0.81396... | -865 | 8 | 0.010586 | 0.083144 |
| 3 | POLYGON ((-77.61444 0.90441, -77.61444 0.99486... | -865 | 10 | 0.010972 | 0.094116 |
| 23 | POLYGON ((-73.66185 5.34332, -73.66185 5.43414... | -821 | 59 | 0.011092 | 0.105208 |
| 5 | POLYGON ((-75.36865 4.70788, -75.36865 4.79862... | -840 | 52 | 0.013525 | 0.118733 |
| 12 | POLYGON ((-74.20084 4.70788, -74.20084 4.79862... | -827 | 52 | 0.016856 | 0.135589 |
| 17 | POLYGON ((-74.02118 4.70788, -74.02118 4.79862... | -825 | 52 | 0.020832 | 0.156421 |
| 20 | POLYGON ((-73.75168 5.2525, -73.75168 5.34332,... | -822 | 58 | 0.023376 | 0.179797 |
| 6 | POLYGON ((-75.27882 4.70788, -75.27882 4.79862... | -839 | 52 | 0.024420 | 0.204217 |
| 11 | POLYGON ((-74.20084 4.61715, -74.20084 4.70788... | -827 | 51 | 0.035877 | 0.240094 |
| 15 | POLYGON ((-74.11101 4.70788, -74.11101 4.79862... | -826 | 52 | 0.053472 | 0.293566 |
| 22 | POLYGON ((-73.75168 5.43414, -73.75168 5.52498... | -822 | 60 | 0.082528 | 0.376094 |
| 24 | POLYGON ((-73.66185 5.43414, -73.66185 5.52498... | -821 | 60 | 0.122049 | 0.498143 |
Calculate cumulative proportion
- Calculate
cumulative proportionby dividingcumulative areabytotal 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 | Zona_Pantanosa_Andina | cumulative_area | cumulative_proportion | |
|---|---|---|---|---|---|---|
| 21 | POLYGON ((-73.75168 5.34332, -73.75168 5.43414... | -822 | 59 | 0.000452 | 0.000452 | 0.000907 |
| 1 | POLYGON ((-77.61444 0.63307, -77.61444 0.72351... | -865 | 7 | 0.001161 | 0.001613 | 0.003238 |
| 0 | POLYGON ((-77.70427 0.90441, -77.70427 0.99486... | -866 | 10 | 0.001515 | 0.003128 | 0.006280 |
| 4 | POLYGON ((-75.45848 6.6162, -75.45848 6.70724,... | -841 | 73 | 0.002899 | 0.006027 | 0.012099 |
| 16 | POLYGON ((-74.02118 4.61715, -74.02118 4.70788... | -825 | 51 | 0.003801 | 0.009827 | 0.019728 |
| 8 | POLYGON ((-74.29067 4.70788, -74.29067 4.79862... | -828 | 52 | 0.003951 | 0.013779 | 0.027660 |
| 26 | POLYGON ((-73.12286 5.70671, -73.12286 5.79759... | -815 | 63 | 0.004846 | 0.018625 | 0.037390 |
| 9 | POLYGON ((-74.20084 4.16367, -74.20084 4.25434... | -827 | 46 | 0.005959 | 0.024585 | 0.049353 |
| 25 | POLYGON ((-73.57202 4.98013, -73.57202 5.07091... | -820 | 55 | 0.005972 | 0.030556 | 0.061341 |
| 13 | POLYGON ((-74.11101 4.52643, -74.11101 4.61715... | -826 | 50 | 0.006286 | 0.036843 | 0.073960 |
| 10 | POLYGON ((-74.20084 4.52643, -74.20084 4.61715... | -827 | 50 | 0.006292 | 0.043135 | 0.086591 |
| 19 | POLYGON ((-73.93135 4.98013, -73.93135 5.07091... | -824 | 55 | 0.006425 | 0.049559 | 0.099488 |
| 7 | POLYGON ((-75.27882 4.79862, -75.27882 4.88937... | -839 | 53 | 0.006715 | 0.056274 | 0.112968 |
| 18 | POLYGON ((-73.93135 4.88937, -73.93135 4.98013... | -824 | 54 | 0.007423 | 0.063697 | 0.127868 |
| 14 | POLYGON ((-74.11101 4.61715, -74.11101 4.70788... | -826 | 51 | 0.008861 | 0.072558 | 0.145657 |
| 2 | POLYGON ((-77.61444 0.72351, -77.61444 0.81396... | -865 | 8 | 0.010586 | 0.083144 | 0.166909 |
| 3 | POLYGON ((-77.61444 0.90441, -77.61444 0.99486... | -865 | 10 | 0.010972 | 0.094116 | 0.188934 |
| 23 | POLYGON ((-73.66185 5.34332, -73.66185 5.43414... | -821 | 59 | 0.011092 | 0.105208 | 0.211199 |
| 5 | POLYGON ((-75.36865 4.70788, -75.36865 4.79862... | -840 | 52 | 0.013525 | 0.118733 | 0.238351 |
| 12 | POLYGON ((-74.20084 4.70788, -74.20084 4.79862... | -827 | 52 | 0.016856 | 0.135589 | 0.272189 |
| 17 | POLYGON ((-74.02118 4.70788, -74.02118 4.79862... | -825 | 52 | 0.020832 | 0.156421 | 0.314009 |
| 20 | POLYGON ((-73.75168 5.2525, -73.75168 5.34332,... | -822 | 58 | 0.023376 | 0.179797 | 0.360935 |
| 6 | POLYGON ((-75.27882 4.70788, -75.27882 4.79862... | -839 | 52 | 0.024420 | 0.204217 | 0.409957 |
| 11 | POLYGON ((-74.20084 4.61715, -74.20084 4.70788... | -827 | 51 | 0.035877 | 0.240094 | 0.481979 |
| 15 | POLYGON ((-74.11101 4.70788, -74.11101 4.79862... | -826 | 52 | 0.053472 | 0.293566 | 0.589321 |
| 22 | POLYGON ((-73.75168 5.43414, -73.75168 5.52498... | -822 | 60 | 0.082528 | 0.376094 | 0.754992 |
| 24 | POLYGON ((-73.66185 5.43414, -73.66185 5.52498... | -821 | 60 | 0.122049 | 0.498143 | 1.000000 |
Count AOO cells
- Calculate AOO by counting the number of cells with a
cumulative proportiongreater 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 24 cellsAOO Calculation (direct call)
if has_data:
aoo_count = ecosystem.aoo
print(f'AOO: {aoo_count} grid cells')AOO: 24 grid cellsCriterion B Summary
Criterion B status (spatial)
Vulnerable (VU) — Zona Pantanosa Andina (Zona Pantanosa Andina), index 82
Status reflects the spatial thresholds for EOO (B1) and AOO (B2) only. A final listing under B1/B2 additionally requires at least one of: (a) an observed or inferred continuing decline; (b) threatening processes likely to cause continuing decline within 20 years; or (c) few threat-defined locations — none of which are derived from the spatial metrics.
| Sub-criterion | Metric | Value | Category |
|---|---|---|---|
| B1 | EOO | 114485 km² | Least Concern (LC) |
| B2 | AOO | 24 cells | Vulnerable (VU) |
| Overall B | — | — | Vulnerable (VU) |