由 NDVI 幂律异速生长方程估算地上生物量碳,叠加根茎比地下碳与类型化土壤碳密度。Estimates carbon stocks from biomass allometry and soil carbon density. 输出地上碳/土壤碳/总碳三张 GeoTIFF 与汇总 JSON。
编程
geoskill-biodiversity-mapping
试用基于生境异质性假说,用 NDVI 生产力、纹理结构异质性与地形粗糙度三类代理估算物种丰富度空间分布。Maps species richness proxies from NDVI, texture and terrain heterogeneity. 输出物种丰富度与生境质量 GeoTIFF + 参数 JSON。
它能做什么
基于生境异质性假说,用 NDVI 生产力、纹理结构异质性与地形粗糙度三类代理估算物种丰富度空间分布。Maps species richness proxies from NDVI, texture and terrain heterogeneity. 输出物种丰富度与生境质量 GeoTIFF + 参数 JSON。
技能文档
生物多样性制图 | Biodiversity Mapping
Vegetation productivity (NDVI), habitat structural diversity (local standard deviation of NDVI), and terrain heterogeneity (gradient magnitude of the DEM) are the most commonly used remote sensing proxies for biodiversity. This skill normalizes each of the three proxies individually and then fuses them with weights into a "habitat quality" index, which is mapped to relative species richness via the saturating curve S = Smax·(1 − exp(−k·q)) to avoid linear extrapolation.
Applicable scenarios: protected area siting, ecological baseline surveys, and identification of priority zones for biodiversity conservation. Supports two weighting schemes: heterogeneity (emphasizing structural/terrain diversity) and productivity (emphasizing energy availability).
Dependencies / 依赖
pip install numpy rasterio scipy
Usage / 使用方法
Example 1: Synthetic Data, Offline Run
python geoskill-biodiversity-mapping.py --bbox 116.0 39.0 117.0 40.0 --synthetic --output-dir ./output
Example 2: Real Multispectral Imagery (band1=Red, band2=NIR, band3=DEM)
python geoskill-biodiversity-mapping.py --input scene.tif --method heterogeneity --output-dir ./real
Example 3: Productivity-Dominant Mode
python geoskill-biodiversity-mapping.py --bbox 121.0 31.0 122.0 32.0 --synthetic --method productivity --output-dir ./prod
Example 4: Larger Texture Window + Higher Saturation Cap
python geoskill-biodiversity-mapping.py --bbox 116 39 117 40 --synthetic --window 9 --s-max 300 --k 2.5 --output-dir ./tuned
Example 5: Quiet Batch Mode
python geoskill-biodiversity-mapping.py --bbox 113 23 114 24 --synthetic --quiet --output-dir ./batch
Output / 输出
| File | Format | Description |
|---|---|---|
species_richness.tif | GeoTIFF (float32) | Relative species richness [0, Smax], EPSG:4326 |
habitat_quality.tif | GeoTIFF (float32) | Habitat quality index q ∈ [0, 1] |
richness_params.json | JSON | Parameters such as weights/window/proxy means |
output-manifest.json | JSON | Run manifest (inputs/outputs/QA/software versions) |
Data Source / 数据源 / Source
Local multispectral GeoTIFF (band1=red, band2=NIR, band3=optional DEM); synthetic mode generates vegetation/water/bare-soil scenes and a rugged DEM locally, with no external data source.
Privacy / 隐私声明 / Privacy
- Runs fully offline by default; makes no network requests
--syntheticmode reads no external data- All computation is performed locally; no user data is uploaded
License / License
MIT
name: geoskill-biodiversity-mapping description: '基于生境异质性假说,用 NDVI 生产力、纹理结构异质性与地形粗糙度三类代理估算物种丰富度空间分布。Maps species richness proxies from NDVI, texture and terrain heterogeneity. 输出物种丰富度与生境质量 GeoTIFF + 参数 JSON。'
生物多样性制图 | Biodiversity Mapping
植被生产力(NDVI)、生境结构多样性(NDVI 局部标准差)与地形异质性(DEM 梯度模)是生物多样性最常用的遥感代理量。本 skill 将三者各自归一化后加权融合为「生境质量」,再经饱和曲线 S = Smax·(1 - exp(-k·q)) 映射为相对物种丰富度,避免线性外推。
适用场景:保护区选址、生态本底调查、生物多样性保护优先区识别。支持 heterogeneity(强调结构/地形多样性)与 productivity(强调能量可得性)两种加权方案。
依赖
pip install numpy rasterio scipy
使用方法
示例 1:合成数据离线运行
python geoskill-biodiversity-mapping.py --bbox 116.0 39.0 117.0 40.0 --synthetic --output-dir ./output
示例 2:真实多光谱影像(band1=红, band2=近红外, band3=DEM)
python geoskill-biodiversity-mapping.py --input scene.tif --method heterogeneity --output-dir ./real
示例 3:生产力主导模式
python geoskill-biodiversity-mapping.py --bbox 121.0 31.0 122.0 32.0 --synthetic --method productivity --output-dir ./prod
示例 4:加大纹理窗口 + 调高饱和上限
python geoskill-biodiversity-mapping.py --bbox 116 39 117 40 --synthetic --window 9 --s-max 300 --k 2.5 --output-dir ./tuned
示例 5:静默批量模式
python geoskill-biodiversity-mapping.py --bbox 113 23 114 24 --synthetic --quiet --output-dir ./batch
输出
| 文件 | 格式 | 说明 |
|---|---|---|
species_richness.tif | GeoTIFF (float32) | 相对物种丰富度 [0, Smax],EPSG:4326 |
habitat_quality.tif | GeoTIFF (float32) | 生境质量指数 q ∈ [0,1] |
richness_params.json | JSON | 权重/窗口/代理量均值等参数 |
output-manifest.json | JSON | 运行清单(输入/输出/QA/软件版本) |
数据源 / Source
本地多光谱 GeoTIFF(band1=红、band2=近红外、band3 可选 DEM);合成模式本地生成植被/水体/裸土场景与起伏 DEM,无外部数据源。
隐私声明 / Privacy
- 默认完全离线运行,不发起任何网络请求
--synthetic模式不读取任何外部数据- 所有计算在本地完成,不上传用户数据
License
MIT
相关技能
由多变量环境栅格训练物种分布模型,输出 0-1 栖息地适宜性概率与变量贡献度。Models habitat suitability probability from environmental rasters with RF or logistic regression. 输出适宜性 GeoTIFF + 模型参数 JSON(含交叉验证 AUC)。
Map green infrastructure from high-resolution NDVI, tree crown detection, green-space classification and patch connectivity.
基于 NDVI 时序积分估算累积生物量,叠加气候校正得到生产力指数。Estimates cumulative biomass from NDVI time-series integration with climate correction to yield a productivity index.
用当量因子法评估供给、调节、支持、文化四类生态系统服务价值。Estimates four ecosystem service values with a simplified InVEST plus equivalent-factor method. 输出四类服务价值 GeoTIFF 与总量 JSON。
多时相指数阈值分类检测入侵新增像元,计算面积相对扩散速率,用环境适宜性×距离衰减预测入侵风险。Monitors invasive species spread from multi-temporal classification and risk prediction. 输出新增入侵与风险 GeoTIFF。