climaclass 0.1.0__py3-none-any.whl

climaclass/__init__.py

@@ -0,0 +1,66 @@
+# SPDX-License-Identifier: Apache-2.0
+# Copyright (c) 2024-2026 Darri Eythorsson
+
+"""climaclass - classical climate classification from monthly climatologies.
+
+Three classical schemes, one tiny input record, zero heavy dependencies:
+
+    >>> from climaclass import MonthlyClimate, classify
+    >>> rome = MonthlyClimate(
+    ...     temp=[8, 9, 11, 14, 18, 22, 25, 25, 21, 17, 12, 9],
+    ...     precip=[80, 75, 70, 65, 50, 35, 20, 30, 70, 110, 110, 95],
+    ... )
+    >>> classify(rome)["koppen"].code
+    'Csa'
+
+Each classifier consumes a :class:`MonthlyClimate` and returns a
+:class:`ClassificationResult`. The library is pure Python; the optional
+SYMFLUENCE integration lives in :mod:`climaclass.integrations.symfluence` and is
+only imported when you ask for it.
+"""
+
+from __future__ import annotations
+
+from typing import Dict
+
+from . import holdridge, koppen, thornthwaite
+from ._types import ClassificationResult, MonthlyClimate
+
+__version__ = "0.1.0"
+
+__all__ = [
+    "MonthlyClimate",
+    "ClassificationResult",
+    "classify",
+    "koppen",
+    "holdridge",
+    "thornthwaite",
+]
+
+_SCHEMES = {
+    "koppen": koppen.classify,
+    "holdridge": holdridge.classify,
+    "thornthwaite": thornthwaite.classify,
+}
+
+
+def classify(climate: MonthlyClimate, schemes=None) -> Dict[str, ClassificationResult]:
+    """Run one or more classification schemes on a climatology.
+
+    Args:
+        climate: The 12-month climatology to classify.
+        schemes: Iterable of scheme names to run. Defaults to all three
+            (``"koppen"``, ``"holdridge"``, ``"thornthwaite"``).
+
+    Returns:
+        Mapping of scheme name -> :class:`ClassificationResult`.
+    """
+    if schemes is None:
+        schemes = _SCHEMES.keys()
+    out: Dict[str, ClassificationResult] = {}
+    for name in schemes:
+        try:
+            out[name] = _SCHEMES[name](climate)
+        except KeyError as exc:
+            raise ValueError(f"Unknown scheme: {name!r}. Choose from {sorted(_SCHEMES)}") from exc
+    return out

climaclass/_types.py

@@ -0,0 +1,100 @@
+# SPDX-License-Identifier: Apache-2.0
+# Copyright (c) 2024-2026 Darri Eythorsson
+
+"""Shared input/output data structures for :mod:`climaclass`.
+
+The whole library is built around one small, explicit input record -
+:class:`MonthlyClimate` - a 12-month climatology of temperature and
+precipitation. Every classifier consumes that record and returns a typed
+result, so the schemes stay interchangeable and trivially testable.
+"""
+
+from __future__ import annotations
+
+from dataclasses import dataclass
+from typing import Optional, Sequence
+
+# Calendar constants. Index 0 == January, 11 == December throughout the package.
+DAYS_IN_MONTH = (31, 28.25, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31)
+# Day-of-year for the middle of each month (used for Thornthwaite day-length).
+MID_MONTH_DOY = (16, 46, 75, 106, 136, 167, 197, 228, 259, 289, 320, 350)
+
+
+def _validate_12(name: str, values: Sequence[float]) -> tuple:
+    vals = tuple(float(v) for v in values)
+    if len(vals) != 12:
+        raise ValueError(f"{name} must have 12 monthly values (Jan..Dec), got {len(vals)}")
+    return vals
+
+
+@dataclass(frozen=True)
+class MonthlyClimate:
+    """A 12-month climatology for a single location or catchment.
+
+    Args:
+        temp: Monthly mean air temperature [degC], January..December.
+        precip: Monthly total precipitation [mm], January..December.
+        tmin: Optional monthly mean daily-minimum temperature [degC]. When
+            absent it is approximated from ``temp`` where a scheme needs it.
+        tmax: Optional monthly mean daily-maximum temperature [degC].
+        latitude: Optional latitude [degrees, north positive]. Only used to
+            refine the Thornthwaite day-length correction and to disambiguate
+            the hemisphere; classification works without it.
+    """
+
+    temp: Sequence[float]
+    precip: Sequence[float]
+    tmin: Optional[Sequence[float]] = None
+    tmax: Optional[Sequence[float]] = None
+    latitude: Optional[float] = None
+
+    def __post_init__(self) -> None:
+        object.__setattr__(self, "temp", _validate_12("temp", self.temp))
+        object.__setattr__(self, "precip", _validate_12("precip", self.precip))
+        if self.tmin is not None:
+            object.__setattr__(self, "tmin", _validate_12("tmin", self.tmin))
+        if self.tmax is not None:
+            object.__setattr__(self, "tmax", _validate_12("tmax", self.tmax))
+
+    # --- convenient derived quantities (shared by several classifiers) ---
+
+    @property
+    def mat(self) -> float:
+        """Mean annual temperature [degC]."""
+        return sum(self.temp) / 12.0
+
+    @property
+    def map(self) -> float:
+        """Mean annual precipitation [mm]."""
+        return sum(self.precip)
+
+    @property
+    def hemisphere(self) -> str:
+        """``"north"`` or ``"south"``, inferred from the seasonal temperature cycle.
+
+        If ``latitude`` is given it wins; otherwise the warmer half-year decides.
+        """
+        if self.latitude is not None:
+            return "north" if self.latitude >= 0 else "south"
+        apr_sep = sum(self.temp[3:9]) / 6.0
+        oct_mar = sum(self.temp[9:12] + self.temp[0:3]) / 6.0
+        return "north" if apr_sep >= oct_mar else "south"
+
+
+@dataclass(frozen=True)
+class ClassificationResult:
+    """A single scheme's verdict.
+
+    Args:
+        scheme: ``"koppen"`` | ``"holdridge"`` | ``"thornthwaite"``.
+        code: Compact symbol (e.g. ``"Csb"``, ``"B"``).
+        zone: Integer index into the scheme's canonical legend (or ``None``).
+        name: Human-readable label (e.g. ``"Tropical rainforest"``).
+        details: Scheme-specific intermediate values, for transparency/QA.
+    """
+
+    scheme: str
+    code: str
+    zone: Optional[int]
+    name: str
+    details: dict

climaclass/holdridge.py

@@ -0,0 +1,140 @@
+# SPDX-License-Identifier: Apache-2.0
+# Copyright (c) 2024-2026 Darri Eythorsson
+
+"""Holdridge life-zone classification from a 12-month climatology.
+
+Follows Holdridge (1967), *Life Zone Ecology*. A life zone is the intersection
+of three logarithmic axes:
+
+* **Mean annual biotemperature** (degC) - mean of monthly temperatures each
+  clamped to [0, 30]. Sets the latitudinal/altitudinal belt.
+* **Mean annual precipitation** (mm).
+* **Potential evapotranspiration ratio**, ``PET_ratio = 58.93 * biotemp / precip``.
+  Sets the humidity province.
+
+This mirrors the variables computed in the original Earth-Engine notebook
+(``biotemperature``, ``pAnn``, ``PETR``) but evaluates them per-climatology in
+pure Python rather than over a global image.
+"""
+
+from __future__ import annotations
+
+from ._types import ClassificationResult, MonthlyClimate
+
+# Biotemperature breakpoints (degC) -> belt name, ordered cold..hot.
+# Each (upper_bound, belt) means: belt applies while biotemp < upper_bound.
+_BELTS = (
+    (1.5, "Polar"),
+    (3.0, "Subpolar"),
+    (6.0, "Boreal"),
+    (12.0, "Cool temperate"),
+    (18.0, "Warm temperate"),
+    (24.0, "Subtropical"),
+    (float("inf"), "Tropical"),
+)
+
+# Humidity provinces by PET ratio, ordered wet..dry. Index used to look up the
+# zone name within a belt. (upper_bound, province) means province applies while
+# PET_ratio < upper_bound.
+_PROVINCES = (
+    (0.25, "superhumid"),
+    (0.50, "perhumid"),
+    (1.00, "humid"),
+    (2.00, "subhumid"),
+    (4.00, "semiarid"),
+    (8.00, "arid"),
+    (16.0, "perarid"),
+    (float("inf"), "superarid"),
+)
+
+# Per-belt life-zone names, indexed by humidity province (0 = wettest).
+# Encodes the Holdridge hexagon; truncated rows repeat the driest realised zone.
+_ZONES = {
+    "Polar": [
+        "Polar desert", "Polar desert", "Polar desert", "Polar desert",
+        "Polar desert", "Polar desert", "Polar desert", "Polar desert",
+    ],
+    "Subpolar": [
+        "Subpolar rain tundra", "Subpolar wet tundra", "Subpolar moist tundra",
+        "Subpolar dry tundra", "Subpolar dry tundra", "Subpolar desert",
+        "Subpolar desert", "Subpolar desert",
+    ],
+    "Boreal": [
+        "Boreal rain forest", "Boreal wet forest", "Boreal moist forest",
+        "Boreal dry scrub", "Boreal desert", "Boreal desert",
+        "Boreal desert", "Boreal desert",
+    ],
+    "Cool temperate": [
+        "Cool temperate rain forest", "Cool temperate wet forest",
+        "Cool temperate moist forest", "Cool temperate steppe",
+        "Cool temperate desert scrub", "Cool temperate desert",
+        "Cool temperate desert", "Cool temperate desert",
+    ],
+    "Warm temperate": [
+        "Warm temperate rain forest", "Warm temperate wet forest",
+        "Warm temperate moist forest", "Warm temperate dry forest",
+        "Warm temperate thorn steppe", "Warm temperate desert scrub",
+        "Warm temperate desert", "Warm temperate desert",
+    ],
+    "Subtropical": [
+        "Subtropical rain forest", "Subtropical wet forest",
+        "Subtropical moist forest", "Subtropical dry forest",
+        "Subtropical thorn woodland", "Subtropical desert scrub",
+        "Subtropical desert", "Subtropical desert",
+    ],
+    "Tropical": [
+        "Tropical rain forest", "Tropical wet forest", "Tropical moist forest",
+        "Tropical dry forest", "Tropical very dry forest",
+        "Tropical thorn woodland", "Tropical desert scrub", "Tropical desert",
+    ],
+}
+
+# Stable integer codes: belt_index * 10 + province_index, plus 1 (1-based).
+_BELT_ORDER = [b for _, b in _BELTS]
+
+
+def biotemperature(climate: MonthlyClimate) -> float:
+    """Mean annual biotemperature [degC]: mean of monthly temps clamped to [0, 30]."""
+    clamped = [min(max(x, 0.0), 30.0) for x in climate.temp]
+    return sum(clamped) / 12.0
+
+
+def _belt(biotemp: float) -> tuple:
+    for idx, (upper, name) in enumerate(_BELTS):
+        if biotemp < upper:
+            return idx, name
+    return len(_BELTS) - 1, _BELTS[-1][1]
+
+
+def _province(pet_ratio: float) -> tuple:
+    for idx, (upper, name) in enumerate(_PROVINCES):
+        if pet_ratio < upper:
+            return idx, name
+    return len(_PROVINCES) - 1, _PROVINCES[-1][1]
+
+
+def classify(climate: MonthlyClimate) -> ClassificationResult:
+    """Classify a 12-month climatology into a Holdridge life zone."""
+    biotemp = biotemperature(climate)
+    precip = climate.map
+    # PET ratio is undefined for zero precipitation; treat as maximally arid.
+    pet_ratio = (58.93 * biotemp / precip) if precip > 0 else float("inf")
+
+    belt_idx, belt_name = _belt(biotemp)
+    prov_idx, prov_name = _province(pet_ratio)
+    name = _ZONES[belt_name][prov_idx]
+    code = belt_idx * 10 + prov_idx + 1
+
+    return ClassificationResult(
+        scheme="holdridge",
+        code=name,
+        zone=code,
+        name=name,
+        details={
+            "biotemperature": round(biotemp, 2),
+            "annual_precip": round(precip, 1),
+            "PET_ratio": round(pet_ratio, 3) if pet_ratio != float("inf") else None,
+            "belt": belt_name,
+            "humidity_province": prov_name,
+        },
+    )

climaclass/integrations/__init__.py

@@ -0,0 +1,9 @@
+# SPDX-License-Identifier: Apache-2.0
+# Copyright (c) 2024-2026 Darri Eythorsson
+
+"""Optional integrations for :mod:`climaclass`.
+
+Nothing here is imported by the core package. Each integration pulls in its
+host framework lazily, so installing ``climaclass`` never drags in SYMFLUENCE
+(and vice versa).
+"""

climaclass/integrations/symfluence.py

@@ -0,0 +1,155 @@
+# SPDX-License-Identifier: Apache-2.0
+# Copyright (c) 2024-2026 Darri Eythorsson
+
+"""SYMFLUENCE attribute-processor adapter for :mod:`climaclass`.
+
+This module lets SYMFLUENCE emit Koppen-Geiger, Holdridge and Thornthwaite
+classes as catchment attributes (``climate.koppen_*`` etc.), computed from the
+same WorldClim monthly rasters SYMFLUENCE already acquires for its climate
+attributes - no Earth Engine, no new data dependency.
+
+It is intentionally decoupled:
+
+* The pure mapping helper :func:`record_to_attributes` has no SYMFLUENCE
+  dependency and can be unit-tested standalone.
+* :class:`ClimateClassificationProcessor` only resolves the SYMFLUENCE base
+  class at import time; if SYMFLUENCE is absent the class still imports (its
+  base degrades to ``object``) so ``import climaclass`` never fails.
+
+Install for use inside SYMFLUENCE with::
+
+    pip install "climaclass[symfluence]"
+
+and register it as an attribute processor (see the project README).
+"""
+
+from __future__ import annotations
+
+from typing import Any, Dict, List, Optional, Sequence
+
+from .. import classify
+from .._types import MonthlyClimate
+
+# Resolve the SYMFLUENCE base class lazily/defensively so importing this module
+# never hard-fails when SYMFLUENCE is not installed.
+try:  # pragma: no cover - exercised only with SYMFLUENCE present
+    from symfluence.data.preprocessing.attribute_processors.base import (
+        BaseAttributeProcessor as _Base,
+    )
+
+    HAVE_SYMFLUENCE = True
+except Exception:  # noqa: BLE001 - any import failure means "not available"
+    _Base = object
+    HAVE_SYMFLUENCE = False
+
+
+def record_to_attributes(
+    temp: Sequence[float],
+    precip: Sequence[float],
+    *,
+    latitude: Optional[float] = None,
+    prefix: str = "",
+) -> Dict[str, Any]:
+    """Classify one HRU's climatology into flat ``climate.*`` attribute keys.
+
+    Args:
+        temp: 12 monthly mean temperatures [degC], Jan..Dec.
+        precip: 12 monthly total precipitation [mm], Jan..Dec.
+        latitude: Optional latitude to refine Thornthwaite PET.
+        prefix: Optional key prefix for distributed domains. SYMFLUENCE expects
+            ``"HRU_{id}_"`` so keys read ``HRU_3_climate.koppen_code``.
+
+    Returns:
+        Dict of attribute name -> value, ready to merge into a results dict.
+    """
+    climate = MonthlyClimate(temp=temp, precip=precip, latitude=latitude)
+    results = classify(climate)
+    k, h, t = results["koppen"], results["holdridge"], results["thornthwaite"]
+    return {
+        f"{prefix}climate.koppen_code": k.code,
+        f"{prefix}climate.koppen_zone": k.zone,
+        f"{prefix}climate.koppen_name": k.name,
+        f"{prefix}climate.holdridge_zone": h.code,
+        f"{prefix}climate.holdridge_code": h.zone,
+        f"{prefix}climate.holdridge_belt": h.details.get("belt"),
+        f"{prefix}climate.thornthwaite_code": t.code,
+        f"{prefix}climate.thornthwaite_moisture_index": t.details.get("moisture_index"),
+        f"{prefix}climate.thornthwaite_moisture_province": t.details.get("moisture_province"),
+    }
+
+
+class ClimateClassificationProcessor(_Base):
+    """SYMFLUENCE attribute processor that adds climate-classification attributes.
+
+    Reads 12 monthly WorldClim mean-temperature (``tavg``) and precipitation
+    (``prec``) rasters, reduces them over the catchment (lumped) or each HRU
+    (distributed) with zonal statistics, and classifies the resulting
+    climatology with all three schemes.
+    """
+
+    #: SYMFLUENCE auto-discovery hooks (kept here so the framework side stays thin).
+    name = "climate_classification"
+    provides = ("climate.koppen_code", "climate.holdridge_zone", "climate.thornthwaite_code")
+
+    def process(self) -> Dict[str, Any]:
+        if not HAVE_SYMFLUENCE:  # pragma: no cover - guard for standalone import
+            raise RuntimeError(
+                "ClimateClassificationProcessor requires SYMFLUENCE. "
+                "Install with: pip install 'climaclass[symfluence]'"
+            )
+
+        results: Dict[str, Any] = {}
+        worldclim_path = self._get_data_path(  # type: ignore[attr-defined]
+            "ATTRIBUTES_WORLDCLIM_PATH", "data/attributes/climate/worldclim"
+        )
+        if not worldclim_path.exists():
+            self.logger.warning(f"WorldClim path not found, skipping classification: {worldclim_path}")
+            return results
+
+        temp = self._monthly_zonal(worldclim_path, "tavg")
+        precip = self._monthly_zonal(worldclim_path, "prec")
+        if temp is None or precip is None:
+            self.logger.warning("Missing tavg/prec WorldClim rasters; skipping classification")
+            return results
+
+        # temp/precip are lists of 12-length climatologies, one per zone.
+        for zone_idx, (t12, p12) in enumerate(zip(temp, precip)):
+            if any(v is None for v in t12) or any(v is None for v in p12):
+                continue
+            prefix = "" if len(temp) == 1 else self._hru_prefix(zone_idx)
+            results.update(record_to_attributes(t12, p12, prefix=prefix))
+        return results
+
+    # --- helpers -----------------------------------------------------------
+
+    def _monthly_zonal(self, worldclim_path, var: str) -> Optional[List[List[Optional[float]]]]:
+        """Return per-zone 12-month means for ``var`` via zonal statistics."""
+        import geopandas as gpd  # lazy, per SYMFLUENCE convention
+        from rasterstats import zonal_stats
+
+        files = sorted(worldclim_path.glob(f"wc2.1_30s_{var}_*.tif"))
+        if len(files) < 12:
+            return None
+
+        catchment = gpd.read_file(self.catchment_path)  # type: ignore[attr-defined]
+        n_zones = len(catchment)
+        monthly: List[List[Optional[float]]] = [[None] * 12 for _ in range(n_zones)]
+        for month_idx, raster in enumerate(files[:12]):
+            stats = zonal_stats(str(self.catchment_path), str(raster), stats=["mean"])  # type: ignore[attr-defined]
+            for zone_idx, s in enumerate(stats):
+                monthly[zone_idx][month_idx] = s.get("mean")
+        return monthly
+
+    def _hru_prefix(self, zone_idx: int) -> str:
+        """Build the per-HRU attribute key prefix using the configured HRU id field."""
+        import geopandas as gpd
+
+        hru_field = self._get_config_value(  # type: ignore[attr-defined]
+            lambda: self.config.paths.catchment_hruid, default="HRU_ID", dict_key="CATCHMENT_SHP_HRUID"
+        )
+        catchment = gpd.read_file(self.catchment_path)  # type: ignore[attr-defined]
+        try:
+            hru_id = catchment.iloc[zone_idx][hru_field]
+        except Exception:  # noqa: BLE001
+            hru_id = zone_idx
+        return f"HRU_{hru_id}_"

climaclass/koppen.py

@@ -0,0 +1,162 @@
+# SPDX-License-Identifier: Apache-2.0
+# Copyright (c) 2024-2026 Darri Eythorsson
+
+"""Koppen-Geiger climate classification from a 12-month climatology.
+
+Implements the widely used formulation of Peel, Finlayson & McMahon (2007,
+*HESS*) and the numeric legend (1..30) of Beck et al. (2018, *Sci. Data*),
+which is the same legend rendered by the original Earth-Engine notebook this
+package grew out of. Pure-Python: no NumPy, no Earth Engine.
+"""
+
+from __future__ import annotations
+
+from typing import Sequence
+
+from ._types import ClassificationResult, MonthlyClimate
+
+# Beck et al. (2018) numeric legend: code string -> integer zone.
+KOPPEN_LEGEND = {
+    "Af": 1, "Am": 2, "Aw": 3,
+    "BWh": 4, "BWk": 5, "BSh": 6, "BSk": 7,
+    "Csa": 8, "Csb": 9, "Csc": 10,
+    "Cwa": 11, "Cwb": 12, "Cwc": 13,
+    "Cfa": 14, "Cfb": 15, "Cfc": 16,
+    "Dsa": 17, "Dsb": 18, "Dsc": 19, "Dsd": 20,
+    "Dwa": 21, "Dwb": 22, "Dwc": 23, "Dwd": 24,
+    "Dfa": 25, "Dfb": 26, "Dfc": 27, "Dfd": 28,
+    "ET": 29, "EF": 30,
+}
+
+KOPPEN_NAMES = {
+    "Af": "Tropical rainforest", "Am": "Tropical monsoon", "Aw": "Tropical savannah",
+    "BWh": "Hot desert", "BWk": "Cold desert", "BSh": "Hot steppe", "BSk": "Cold steppe",
+    "Csa": "Hot-summer Mediterranean", "Csb": "Warm-summer Mediterranean",
+    "Csc": "Cold-summer Mediterranean",
+    "Cwa": "Humid subtropical (dry winter)", "Cwb": "Subtropical highland (dry winter)",
+    "Cwc": "Cold subtropical highland (dry winter)",
+    "Cfa": "Humid subtropical", "Cfb": "Temperate oceanic", "Cfc": "Subpolar oceanic",
+    "Dsa": "Hot-summer Mediterranean continental", "Dsb": "Warm-summer Mediterranean continental",
+    "Dsc": "Mediterranean subarctic", "Dsd": "Mediterranean subarctic (severe winter)",
+    "Dwa": "Hot-summer humid continental (dry winter)",
+    "Dwb": "Warm-summer humid continental (dry winter)",
+    "Dwc": "Subarctic (dry winter)", "Dwd": "Severe-winter subarctic (dry winter)",
+    "Dfa": "Hot-summer humid continental", "Dfb": "Warm-summer humid continental",
+    "Dfc": "Subarctic", "Dfd": "Severe-winter subarctic",
+    "ET": "Tundra", "EF": "Ice cap",
+}
+
+# Northern-hemisphere summer = Apr..Sep (indices 3..8); winter = Oct..Mar.
+_NH_SUMMER = (3, 4, 5, 6, 7, 8)
+_NH_WINTER = (9, 10, 11, 0, 1, 2)
+
+
+def _seasonal_indices(hemisphere: str) -> tuple:
+    """Return (summer_months, winter_months) as 0-based index tuples."""
+    if hemisphere == "north":
+        return _NH_SUMMER, _NH_WINTER
+    return _NH_WINTER, _NH_SUMMER
+
+
+def classify(climate: MonthlyClimate) -> ClassificationResult:
+    """Classify a 12-month climatology into a Koppen-Geiger zone."""
+    t = climate.temp
+    p = climate.precip
+    mat = climate.mat
+    map_ = climate.map
+
+    t_cold = min(t)
+    t_hot = max(t)
+    t_mon10 = sum(1 for x in t if x >= 10.0)
+    p_dry = min(p)
+
+    summer, winter = _seasonal_indices(climate.hemisphere)
+    p_summer = sum(p[i] for i in summer)
+    p_winter = sum(p[i] for i in winter)
+    p_s_dry = min(p[i] for i in summer)
+    p_w_dry = min(p[i] for i in winter)
+    p_s_wet = max(p[i] for i in summer)
+    p_w_wet = max(p[i] for i in winter)
+
+    # Aridity threshold P_th (Peel et al. 2007).
+    if p_winter >= 0.7 * map_:
+        p_th = 2.0 * mat
+    elif p_summer >= 0.7 * map_:
+        p_th = 2.0 * mat + 28.0
+    else:
+        p_th = 2.0 * mat + 14.0
+
+    details = {
+        "MAT": round(mat, 2), "MAP": round(map_, 1),
+        "Tcold": round(t_cold, 2), "Thot": round(t_hot, 2),
+        "Tmon10": t_mon10, "Pdry": round(p_dry, 1),
+        "Pthreshold": round(p_th, 2), "hemisphere": climate.hemisphere,
+    }
+
+    # --- Main class. Aridity (B) is tested first and overrides A/C/D; E is the
+    #     cold residual. This ordering reproduces cold deserts (BWk) as mapped
+    #     by Beck et al. (2018). ---
+    if map_ < 10.0 * p_th:
+        code = _arid(map_, p_th, mat)
+    elif t_cold >= 18.0:
+        code = _tropical(p_dry, map_)
+    elif t_hot > 10.0 and t_cold > 0.0:
+        code = "C" + _cd_second(p_s_dry, p_w_wet, p_w_dry, p_s_wet) + _c_third(t_hot, t_mon10)
+    elif t_hot > 10.0 and t_cold <= 0.0:
+        code = "D" + _cd_second(p_s_dry, p_w_wet, p_w_dry, p_s_wet) + _d_third(t_hot, t_mon10, t_cold)
+    else:
+        code = "ET" if t_hot >= 0.0 else "EF"
+
+    return ClassificationResult(
+        scheme="koppen",
+        code=code,
+        zone=KOPPEN_LEGEND.get(code),
+        name=KOPPEN_NAMES.get(code, code),
+        details=details,
+    )
+
+
+def _arid(map_: float, p_th: float, mat: float) -> str:
+    second = "W" if map_ < 5.0 * p_th else "S"
+    third = "h" if mat >= 18.0 else "k"
+    return "B" + second + third
+
+
+def _tropical(p_dry: float, map_: float) -> str:
+    if p_dry >= 60.0:
+        return "Af"
+    if p_dry >= 100.0 - map_ / 25.0:
+        return "Am"
+    return "Aw"
+
+
+def _cd_second(p_s_dry: float, p_w_wet: float, p_w_dry: float, p_s_wet: float) -> str:
+    """Second letter for C/D climates: s (dry summer), w (dry winter) or f."""
+    dry_summer = p_s_dry < 40.0 and p_s_dry < p_w_wet / 3.0
+    dry_winter = p_w_dry < p_s_wet / 10.0
+    if dry_summer and not dry_winter:
+        return "s"
+    if dry_winter and not dry_summer:
+        return "w"
+    # If both criteria trip, assign by the wetter season (Peel et al. 2007).
+    if dry_summer and dry_winter:
+        return "s" if p_w_wet >= p_s_wet else "w"
+    return "f"
+
+
+def _c_third(t_hot: float, t_mon10: int) -> str:
+    if t_hot >= 22.0:
+        return "a"
+    if t_mon10 >= 4:
+        return "b"
+    return "c"
+
+
+def _d_third(t_hot: float, t_mon10: int, t_cold: float) -> str:
+    if t_hot >= 22.0:
+        return "a"
+    if t_mon10 >= 4:
+        return "b"
+    if t_cold <= -38.0:
+        return "d"
+    return "c"

climaclass/thornthwaite.py

@@ -0,0 +1,125 @@
+# SPDX-License-Identifier: Apache-2.0
+# Copyright (c) 2024-2026 Darri Eythorsson
+
+"""Thornthwaite climate classification from a 12-month climatology.
+
+Implements the Thornthwaite (1948) potential-evapotranspiration model and the
+moisture index. Potential evapotranspiration (PET) is computed from monthly
+mean temperature via the classic heat-index formula, with the high-temperature
+correction for T >= 26.5 degC and an optional latitude-based day-length
+adjustment. Classification uses:
+
+* **Moisture index** ``Im = 100 * (P - PET) / PET`` (the 1955 single-index form)
+  -> humidity province (A perhumid .. E arid).
+* **Thermal efficiency** = annual PET [mm] -> thermal province (A' .. E').
+
+Pure Python; ``latitude`` only refines the day-length correction.
+"""
+
+from __future__ import annotations
+
+import math
+
+from ._types import DAYS_IN_MONTH, MID_MONTH_DOY, ClassificationResult, MonthlyClimate
+
+# Moisture index thresholds (Im), wet..dry: (lower_bound, code, name).
+_MOISTURE = (
+    (100.0, "A", "Perhumid"),
+    (80.0, "B4", "Humid (B4)"),
+    (60.0, "B3", "Humid (B3)"),
+    (40.0, "B2", "Humid (B2)"),
+    (20.0, "B1", "Humid (B1)"),
+    (0.0, "C2", "Moist subhumid"),
+    (-33.3, "C1", "Dry subhumid"),
+    (-66.7, "D", "Semiarid"),
+    (float("-inf"), "E", "Arid"),
+)
+
+# Thermal-efficiency thresholds (annual PET, mm), warm..cold.
+_THERMAL = (
+    (1140.0, "A'", "Megathermal"),
+    (997.0, "B'4", "Mesothermal (B'4)"),
+    (855.0, "B'3", "Mesothermal (B'3)"),
+    (712.0, "B'2", "Mesothermal (B'2)"),
+    (570.0, "B'1", "Mesothermal (B'1)"),
+    (427.0, "C'2", "Microthermal (C'2)"),
+    (285.0, "C'1", "Microthermal (C'1)"),
+    (142.0, "D'", "Tundra"),
+    (float("-inf"), "E'", "Frost"),
+)
+
+
+def heat_index(temp: tuple) -> float:
+    """Annual heat index I = sum((T/5)^1.514) over months with T > 0."""
+    return sum((x / 5.0) ** 1.514 for x in temp if x > 0.0)
+
+
+def _alpha(i: float) -> float:
+    return 6.75e-7 * i**3 - 7.71e-5 * i**2 + 1.792e-2 * i + 0.49239
+
+
+def _daylength_factor(latitude: float, month: int) -> float:
+    """Correction factor (N/12)*(days/30) for the Thornthwaite PET adjustment."""
+    phi = math.radians(latitude)
+    decl = 0.409 * math.sin(2.0 * math.pi / 365.0 * MID_MONTH_DOY[month] - 1.39)
+    x = -math.tan(phi) * math.tan(decl)
+    x = max(-1.0, min(1.0, x))  # clamp for polar day/night
+    omega = math.acos(x)
+    daylight_hours = 24.0 / math.pi * omega
+    return (daylight_hours / 12.0) * (DAYS_IN_MONTH[month] / 30.0)
+
+
+def monthly_pet(climate: MonthlyClimate) -> list:
+    """Monthly potential evapotranspiration [mm], Thornthwaite (1948)."""
+    temp = climate.temp
+    i = heat_index(temp)
+    a = _alpha(i)
+    pet = []
+    for m, t in enumerate(temp):
+        if t <= 0.0:
+            unadjusted = 0.0
+        elif t < 26.5:
+            unadjusted = 16.0 * (10.0 * t / i) ** a if i > 0 else 0.0
+        else:
+            # High-temperature branch: PET saturates (Willmott et al. 1985).
+            unadjusted = -415.85 + 32.24 * t - 0.43 * t * t
+        if climate.latitude is not None:
+            unadjusted *= _daylength_factor(climate.latitude, m)
+        pet.append(max(unadjusted, 0.0))
+    return pet
+
+
+def _lookup(table: tuple, value: float) -> tuple:
+    for lower, code, name in table:
+        if value >= lower:
+            return code, name
+    return table[-1][1], table[-1][2]
+
+
+def classify(climate: MonthlyClimate) -> ClassificationResult:
+    """Classify a 12-month climatology by Thornthwaite moisture & thermal indices."""
+    pet = monthly_pet(climate)
+    annual_pet = sum(pet)
+    precip = climate.map
+
+    moisture_index = 100.0 * (precip - annual_pet) / annual_pet if annual_pet > 0 else float("inf")
+
+    m_code, m_name = _lookup(_MOISTURE, moisture_index)
+    t_code, t_name = _lookup(_THERMAL, annual_pet)
+
+    code = f"{m_code}{t_code}"
+    name = f"{m_name} / {t_name}"
+
+    return ClassificationResult(
+        scheme="thornthwaite",
+        code=code,
+        zone=None,
+        name=name,
+        details={
+            "moisture_index": round(moisture_index, 1) if moisture_index != float("inf") else None,
+            "annual_PET": round(annual_pet, 1),
+            "annual_precip": round(precip, 1),
+            "moisture_province": m_name,
+            "thermal_province": t_name,
+        },
+    )

climaclass-0.1.0.dist-info/METADATA

@@ -0,0 +1,111 @@
+Metadata-Version: 2.4
+Name: climaclass
+Version: 0.1.0
+Summary: Classical climate classification (Koppen-Geiger, Holdridge, Thornthwaite) from monthly climatologies - pure Python, no Earth Engine required.
+Project-URL: Homepage, https://github.com/DarriEy/Climate_Classifications
+Project-URL: Source, https://github.com/DarriEy/Climate_Classifications
+Author: Darri Eythorsson
+License-Expression: Apache-2.0
+License-File: LICENSE
+License-File: NOTICE
+Keywords: classification,climate,climatology,holdridge,hydrology,koppen,thornthwaite
+Classifier: Development Status :: 4 - Beta
+Classifier: Intended Audience :: Science/Research
+Classifier: Programming Language :: Python :: 3
+Classifier: Topic :: Scientific/Engineering :: Atmospheric Science
+Requires-Python: >=3.9
+Provides-Extra: dev
+Requires-Dist: numpy; extra == 'dev'
+Requires-Dist: pytest>=7; extra == 'dev'
+Requires-Dist: ruff>=0.1; extra == 'dev'
+Provides-Extra: symfluence
+Requires-Dist: geopandas>=0.12; extra == 'symfluence'
+Requires-Dist: rasterstats>=0.18; extra == 'symfluence'
+Description-Content-Type: text/markdown
+
+# climaclass
+
+**Classical climate classification from monthly climatologies — pure Python, no Earth Engine required.**
+
+`climaclass` turns a 12-month climatology of temperature and precipitation into a
+climate class under three classical schemes:
+
+| Scheme | Reference | Output |
+|--------|-----------|--------|
+| **Köppen–Geiger** | Peel et al. (2007); legend of Beck et al. (2018) | code (`Csb`) + zone 1–30 |
+| **Holdridge life zones** | Holdridge (1967) | life-zone name + belt/province |
+| **Thornthwaite** | Thornthwaite (1948) | moisture + thermal province |
+
+It grew out of the [`Climate_Classifications`](https://github.com/DarriEy/Climate_Classifications)
+Earth-Engine notebooks, but reimplements the *algorithms* as dependency-free
+Python so they run anywhere — a laptop, a CI job, or inside a hydrological model
+pipeline — without a Google Earth Engine account.
+
+## Install
+
+```bash
+pip install climaclass                 # core: zero dependencies
+pip install "climaclass[symfluence]"   # + SYMFLUENCE attribute-processor plugin
+```
+
+## Use
+
+```python
+from climaclass import MonthlyClimate, classify
+
+reykjavik = MonthlyClimate(
+    temp=[0.1, 0.4, 0.9, 3.3, 6.8, 9.4, 10.9, 10.5, 7.7, 4.5, 1.7, 0.4],   # °C, Jan..Dec
+    precip=[76, 72, 82, 58, 44, 50, 52, 62, 67, 86, 73, 79],               # mm,  Jan..Dec
+    latitude=64.1,
+)
+
+results = classify(reykjavik)
+print(results["koppen"].code)        # 'Cfc'  (subpolar oceanic)
+print(results["holdridge"].name)     # e.g. 'Boreal moist forest'
+print(results["thornthwaite"].code)  # moisture/thermal province code
+print(results["koppen"].details)     # MAT, MAP, Pthreshold, ... for QA
+```
+
+Run a single scheme directly:
+
+```python
+from climaclass import koppen
+koppen.classify(reykjavik).zone      # 16  (Beck et al. legend)
+```
+
+### Input
+
+Everything keys off one small record, `MonthlyClimate`:
+
+- `temp` — 12 monthly **mean** temperatures (°C), January → December *(required)*
+- `precip` — 12 monthly **total** precipitation (mm), January → December *(required)*
+- `latitude` — optional; refines Thornthwaite PET day-length and fixes hemisphere
+- `tmin` / `tmax` — optional, reserved for future scheme variants
+
+## SYMFLUENCE plugin
+
+The optional `symfluence` extra ships an attribute processor that emits climate
+classes as catchment attributes (`climate.koppen_code`, `climate.holdridge_zone`,
+`climate.thornthwaite_code`, …), computed from the **WorldClim** monthly rasters
+SYMFLUENCE already acquires — so there is no new data dependency and no Earth
+Engine in the loop. It is registered via the `symfluence.attribute_processors`
+entry point and can be used as a regionalization / PUB grouping variable.
+
+The integration is fully decoupled: importing `climaclass` never imports
+SYMFLUENCE, and the pure mapping helper `record_to_attributes(temp, precip)` is
+testable on its own.
+
+## Develop
+
+```bash
+pip install -e ".[dev]"
+pytest -q
+ruff check src/
+```
+
+## License
+
+Apache License 2.0 — see [LICENSE](LICENSE) and [NOTICE](NOTICE). The
+classification algorithms were prototyped in the author's
+`Climate_Classifications` Earth-Engine notebooks and relicensed by the
+copyright holder.

climaclass-0.1.0.dist-info/RECORD

@@ -0,0 +1,13 @@
+climaclass/__init__.py,sha256=ynQebYMVh72NbhKwuEAm8xfdsLXjKcp0VoGirvPg5Xk,1999
+climaclass/_types.py,sha256=FCS6nFb4hsEXhmKx9_2QlqrxDN0-mcST2gM-IfrJI-U,3739
+climaclass/holdridge.py,sha256=7PGNps0TjhDnQeebd_NwfL1wenwe5Pejyrm6F4dR6Vw,5102
+climaclass/koppen.py,sha256=e5lA1ZU79ZYTBlKVflrHfasbocLv4AXY9bx4h7Xsm90,5683
+climaclass/thornthwaite.py,sha256=u-7vJ80urNnKB4On6ByHSxz7TKSTcWNZ_pD5a58AGNU,4295
+climaclass/integrations/__init__.py,sha256=p6Olc_nkLJuhlgS8vWpiydu7erswCZ56ZV2He1LCVDA,307
+climaclass/integrations/symfluence.py,sha256=Kw3xmnzBpABbCOeRmc2HrP2Nj6R7j_t13-fqhhIOaJ8,6662
+climaclass-0.1.0.dist-info/METADATA,sha256=vFPDwm5o_1JQgj5oPFkS0WrNsgqFSb3gieTaN1D6jlM,4267
+climaclass-0.1.0.dist-info/WHEEL,sha256=QccIxa26bgl1E6uMy58deGWi-0aeIkkangHcxk2kWfw,87
+climaclass-0.1.0.dist-info/entry_points.txt,sha256=DBRgJhrjsRsg2o9FlXNjYFRX7phA6WcjFQ5Zofh9IUI,125
+climaclass-0.1.0.dist-info/licenses/LICENSE,sha256=Bb4OTVdNvJJ-JZSWB_tbpHnobJzSRUv-TugVl_E68vs,11351
+climaclass-0.1.0.dist-info/licenses/NOTICE,sha256=Y5zFcnG7UBwm_rNCL-H92cHWoRp-Ep1TEAuu41XBfqg,861
+climaclass-0.1.0.dist-info/RECORD,,

climaclass-0.1.0.dist-info/WHEEL

@@ -0,0 +1,4 @@
+Wheel-Version: 1.0
+Generator: hatchling 1.29.0
+Root-Is-Purelib: true
+Tag: py3-none-any

climaclass-0.1.0.dist-info/entry_points.txt

@@ -0,0 +1,2 @@
+[symfluence.attribute_processors]
+climate_classification = climaclass.integrations.symfluence:ClimateClassificationProcessor

climaclass-0.1.0.dist-info/licenses/LICENSE

@@ -0,0 +1,201 @@
+                                 Apache License
+                           Version 2.0, January 2004
+                        http://www.apache.org/licenses/
+
+   TERMS AND CONDITIONS FOR USE, REPRODUCTION, AND DISTRIBUTION
+
+   1. Definitions.
+
+      "License" shall mean the terms and conditions for use, reproduction,
+      and distribution as defined by Sections 1 through 9 of this document.
+
+      "Licensor" shall mean the copyright owner or entity authorized by
+      the copyright owner that is granting the License.
+
+      "Legal Entity" shall mean the union of the acting entity and all
+      other entities that control, are controlled by, or are under common
+      control with that entity. For the purposes of this definition,
+      "control" means (i) the power, direct or indirect, to cause the
+      direction or management of such entity, whether by contract or
+      otherwise, or (ii) ownership of fifty percent (50%) or more of the
+      outstanding shares, or (iii) beneficial ownership of such entity.
+
+      "You" (or "Your") shall mean an individual or Legal Entity
+      exercising permissions granted by this License.
+
+      "Source" form shall mean the preferred form for making modifications,
+      including but not limited to software source code, documentation
+      source, and configuration files.
+
+      "Object" form shall mean any form resulting from mechanical
+      transformation or translation of a Source form, including but
+      not limited to compiled object code, generated documentation,
+      and conversions to other media types.
+
+      "Work" shall mean the work of authorship, whether in Source or
+      Object form, made available under the License, as indicated by a
+      copyright notice that is included in or attached to the work
+      (an example is provided in the Appendix below).
+
+      "Derivative Works" shall mean any work, whether in Source or Object
+      form, that is based on (or derived from) the Work and for which the
+      editorial revisions, annotations, elaborations, or other modifications
+      represent, as a whole, an original work of authorship. For the purposes
+      of this License, Derivative Works shall not include works that remain
+      separable from, or merely link (or bind by name) to the interfaces of,
+      the Work and Derivative Works thereof.
+
+      "Contribution" shall mean any work of authorship, including
+      the original version of the Work and any modifications or additions
+      to that Work or Derivative Works thereof, that is intentionally
+      submitted to Licensor for inclusion in the Work by the copyright owner
+      or by an individual or Legal Entity authorized to submit on behalf of
+      the copyright owner. For the purposes of this definition, "submitted"
+      means any form of electronic, verbal, or written communication sent
+      to the Licensor or its representatives, including but not limited to
+      communication on electronic mailing lists, source code control systems,
+      and issue tracking systems that are managed by, or on behalf of, the
+      Licensor for the purpose of discussing and improving the Work, but
+      excluding communication that is conspicuously marked or otherwise
+      designated in writing by the copyright owner as "Not a Contribution."
+
+      "Contributor" shall mean Licensor and any individual or Legal Entity
+      on behalf of whom a Contribution has been received by Licensor and
+      subsequently incorporated within the Work.
+
+   2. Grant of Copyright License. Subject to the terms and conditions of
+      this License, each Contributor hereby grants to You a perpetual,
+      worldwide, non-exclusive, no-charge, royalty-free, irrevocable
+      copyright license to reproduce, prepare Derivative Works of,
+      publicly display, publicly perform, sublicense, and distribute the
+      Work and such Derivative Works in Source or Object form.
+
+   3. Grant of Patent License. Subject to the terms and conditions of
+      this License, each Contributor hereby grants to You a perpetual,
+      worldwide, non-exclusive, no-charge, royalty-free, irrevocable
+      (except as stated in this section) patent license to make, have made,
+      use, offer to sell, sell, import, and otherwise transfer the Work,
+      where such license applies only to those patent claims licensable
+      by such Contributor that are necessarily infringed by their
+      Contribution(s) alone or by combination of their Contribution(s)
+      with the Work to which such Contribution(s) was submitted. If You
+      institute patent litigation against any entity (including a
+      cross-claim or counterclaim in a lawsuit) alleging that the Work
+      or a Contribution incorporated within the Work constitutes direct
+      or contributory patent infringement, then any patent licenses
+      granted to You under this License for that Work shall terminate
+      as of the date such litigation is filed.
+
+   4. Redistribution. You may reproduce and distribute copies of the
+      Work or Derivative Works thereof in any medium, with or without
+      modifications, and in Source or Object form, provided that You
+      meet the following conditions:
+
+      (a) You must give any other recipients of the Work or Derivative
+          Works a copy of this License; and
+
+      (b) You must cause any modified files to carry prominent notices
+          stating that You changed the files; and
+
+      (c) You must retain, in the Source form of any Derivative Works
+          that You distribute, all copyright, patent, trademark, and
+          attribution notices from the Source form of the Work,
+          excluding those notices that do not pertain to any part of
+          the Derivative Works; and
+
+      (d) If the Work includes a "NOTICE" text file as part of its
+          distribution, then any Derivative Works that You distribute must
+          include a readable copy of the attribution notices contained
+          within such NOTICE file, excluding those notices that do not
+          pertain to any part of the Derivative Works, in at least one
+          of the following places: within a NOTICE text file distributed
+          as part of the Derivative Works; within the Source form or
+          documentation, if provided along with the Derivative Works; or,
+          within a display generated by the Derivative Works, if and
+          wherever such third-party notices normally appear. The contents
+          of the NOTICE file are for informational purposes only and
+          do not modify the License. You may add Your own attribution
+          notices within Derivative Works that You distribute, alongside
+          or as an addendum to the NOTICE text from the Work, provided
+          that such additional attribution notices cannot be construed
+          as modifying the License.
+
+      You may add Your own copyright statement to Your modifications and
+      may provide additional or different license terms and conditions
+      for use, reproduction, or distribution of Your modifications, or
+      for any such Derivative Works as a whole, provided Your use,
+      reproduction, and distribution of the Work otherwise complies with
+      the conditions stated in this License.
+
+   5. Submission of Contributions. Unless You explicitly state otherwise,
+      any Contribution intentionally submitted for inclusion in the Work
+      by You to the Licensor shall be under the terms and conditions of
+      this License, without any additional terms or conditions.
+      Notwithstanding the above, nothing herein shall supersede or modify
+      the terms of any separate license agreement you may have executed
+      with Licensor regarding such Contributions.
+
+   6. Trademarks. This License does not grant permission to use the trade
+      names, trademarks, service marks, or product names of the Licensor,
+      except as required for reasonable and customary use in describing the
+      origin of the Work and reproducing the content of the NOTICE file.
+
+   7. Disclaimer of Warranty. Unless required by applicable law or
+      agreed to in writing, Licensor provides the Work (and each
+      Contributor provides its Contributions) on an "AS IS" BASIS,
+      WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or
+      implied, including, without limitation, any warranties or conditions
+      of TITLE, NON-INFRINGEMENT, MERCHANTABILITY, or FITNESS FOR A
+      PARTICULAR PURPOSE. You are solely responsible for determining the
+      appropriateness of using or redistributing the Work and assume any
+      risks associated with Your exercise of permissions under this License.
+
+   8. Limitation of Liability. In no event and under no legal theory,
+      whether in tort (including negligence), contract, or otherwise,
+      unless required by applicable law (such as deliberate and grossly
+      negligent acts) or agreed to in writing, shall any Contributor be
+      liable to You for damages, including any direct, indirect, special,
+      incidental, or consequential damages of any character arising as a
+      result of this License or out of the use or inability to use the
+      Work (including but not limited to damages for loss of goodwill,
+      work stoppage, computer failure or malfunction, or any and all
+      other commercial damages or losses), even if such Contributor
+      has been advised of the possibility of such damages.
+
+   9. Accepting Warranty or Additional Liability. While redistributing
+      the Work or Derivative Works thereof, You may choose to offer,
+      and charge a fee for, acceptance of support, warranty, indemnity,
+      or other liability obligations and/or rights consistent with this
+      License. However, in accepting such obligations, You may act only
+      on Your own behalf and on Your sole responsibility, not on behalf
+      of any other Contributor, and only if You agree to indemnify,
+      defend, and hold each Contributor harmless for any liability
+      incurred by, or claims asserted against, such Contributor by reason
+      of your accepting any such warranty or additional liability.
+
+   END OF TERMS AND CONDITIONS
+
+   APPENDIX: How to apply the Apache License to your work.
+
+      To apply the Apache License to your work, attach the following
+      boilerplate notice, with the fields enclosed by brackets "[]"
+      replaced with your own identifying information. (Don't include
+      the brackets!)  The text should be enclosed in the appropriate
+      comment syntax for the file format. We also recommend that a
+      file or class name and description of purpose be included on the
+      same "printed page" as the copyright notice for easier
+      identification within third-party archives.
+
+   Copyright 2024-2026 Darri Eythorsson
+
+   Licensed under the Apache License, Version 2.0 (the "License");
+   you may not use this file except in compliance with the License.
+   You may obtain a copy of the License at
+
+       http://www.apache.org/licenses/LICENSE-2.0
+
+   Unless required by applicable law or agreed to in writing, software
+   distributed under the License is distributed on an "AS IS" BASIS,
+   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+   See the License for the specific language governing permissions and
+   limitations under the License.

climaclass-0.1.0.dist-info/licenses/NOTICE

@@ -0,0 +1,17 @@
+climaclass
+Copyright 2024-2026 Darri Eythorsson
+
+This product includes software developed by Darri Eythorsson.
+
+The classification algorithms reimplement the methods previously prototyped in
+the author's Climate_Classifications project (Google Earth Engine notebooks),
+relicensed by the copyright holder under the Apache License, Version 2.0.
+
+Scientific references:
+  - Peel, M. C., Finlayson, B. L., & McMahon, T. A. (2007). Updated world map
+    of the Koppen-Geiger climate classification. HESS, 11, 1633-1644.
+  - Beck, H. E., et al. (2018). Present and future Koppen-Geiger climate
+    classification maps at 1-km resolution. Scientific Data, 5, 180214.
+  - Holdridge, L. R. (1967). Life Zone Ecology. Tropical Science Center.
+  - Thornthwaite, C. W. (1948). An approach toward a rational classification of
+    climate. Geographical Review, 38(1), 55-94.