Coding

geoskill-ecosystem-services-valuation

Try it

用当量因子法评估供给、调节、支持、文化四类生态系统服务价值。Estimates four ecosystem service values with a simplified InVEST plus equivalent-factor method. 输出四类服务价值 GeoTIFF 与总量 JSON。

What it does

用当量因子法评估供给、调节、支持、文化四类生态系统服务价值。Estimates four ecosystem service values with a simplified InVEST plus equivalent-factor method. 输出四类服务价值 GeoTIFF 与总量 JSON。

The skill document

生态系统服务评估 | Ecosystem Services Valuation

Based on the equivalent-factor table of Xie Gaodi et al. (1 equivalent ≈ national mean grain production value of 3,000 CNY/ha/yr): five LULC classes — forest, grassland, cropland, water, and built-up land — are first derived from NDVI thresholds; the class × service equivalent coefficient is then multiplied by the pixel area to produce annual value-density rasters and regional totals for the four service categories: provisioning, regulating, supporting, and cultural.

Use cases: ecological asset accounting, Gross Ecosystem Product (GEP) estimation, and calculation of requisition–compensation balance and ecological compensation standards.

Dependencies / 依赖

pip install numpy rasterio geopandas shapely

Usage / 使用方法

Example 1: synthetic NDVI scenario

python geoskill-ecosystem-services-valuation.py --bbox 116.0 39.0 117.0 40.0 --synthetic --output-dir ./output

Example 2: real NDVI raster

python geoskill-ecosystem-services-valuation.py --input ndvi.tif --output-dir ./real

Example 3: different region (Shanghai)

python geoskill-ecosystem-services-valuation.py --bbox 121.0 31.0 122.0 32.0 --synthetic --output-dir ./shanghai

Example 4: silent batch run

python geoskill-ecosystem-services-valuation.py --bbox 113 23 114 24 --synthetic --quiet --output-dir ./batch

Example 5: tiny region for quick validation

python geoskill-ecosystem-services-valuation.py --bbox 116.39 39.90 116.40 39.91 --synthetic --output-dir ./tiny

Output / 输出

FileFormatDescription
value_provisioning.tifGeoTIFF (float32)Provisioning service value (CNY/yr/pixel)
value_regulating.tifGeoTIFF (float32)Regulating service value (CNY/yr/pixel)
value_supporting.tifGeoTIFF (float32)Supporting service value (CNY/yr/pixel)
value_cultural.tifGeoTIFF (float32)Cultural service value (CNY/yr/pixel)
service_value_params.jsonJSONPixel area, per-service totals, LULC pixel counts
output-manifest.jsonJSONRun manifest (inputs/outputs/QA/software versions)

Data Source / 数据源 / Source

Local GeoTIFF (band1 = NDVI); the equivalent-factor table follows the China ecosystem service value equivalent factors published by Xie Gaodi et al. (2015); synthetic mode generates data locally with no external data sources.

Privacy / 隐私声明 / Privacy

  • Runs fully offline by default and makes no network requests
  • --synthetic mode reads no external data
  • All computation is done locally; user data is never uploaded

License / License

MIT



name: geoskill-ecosystem-services-valuation description: '用当量因子法评估供给、调节、支持、文化四类生态系统服务价值。Estimates four ecosystem service values with a simplified InVEST plus equivalent-factor method. 输出四类服务价值 GeoTIFF 与总量 JSON。'

生态系统服务评估 | Ecosystem Services Valuation

基于谢高地等当量因子表(1 当量 ≈ 全国均值粮食产值 3000 元/ha/yr):先从 NDVI 阈值反演林/草/耕/水/建设用地五类 LULC,再按类别×服务的当量系数乘以像元面积,得到供给、调节、支持、文化四类服务的年价值密度栅格与区域总量。

适用场景:生态资产核算、GEP 估算、占补平衡与生态补偿标准测算。

依赖

pip install numpy rasterio geopandas shapely

使用方法

示例 1:合成 NDVI 场景

python geoskill-ecosystem-services-valuation.py --bbox 116.0 39.0 117.0 40.0 --synthetic --output-dir ./output

示例 2:真实 NDVI 栅格

python geoskill-ecosystem-services-valuation.py --input ndvi.tif --output-dir ./real

示例 3:不同区域(上海)

python geoskill-ecosystem-services-valuation.py --bbox 121.0 31.0 122.0 32.0 --synthetic --output-dir ./shanghai

示例 4:静默批量

python geoskill-ecosystem-services-valuation.py --bbox 113 23 114 24 --synthetic --quiet --output-dir ./batch

示例 5:极小区域快速验证

python geoskill-ecosystem-services-valuation.py --bbox 116.39 39.90 116.40 39.91 --synthetic --output-dir ./tiny

输出

文件格式说明
value_provisioning.tifGeoTIFF (float32)供给服务价值(元/yr/像元)
value_regulating.tifGeoTIFF (float32)调节服务价值(元/yr/像元)
value_supporting.tifGeoTIFF (float32)支持服务价值(元/yr/像元)
value_cultural.tifGeoTIFF (float32)文化服务价值(元/yr/像元)
service_value_params.jsonJSON像元面积、各服务总量、LULC 像元计数
output-manifest.jsonJSON运行清单(输入/输出/QA/软件版本)

数据源 / Source

本地 GeoTIFF(band1=NDVI);当量因子表参考谢高地等(2015)公开发表的中国生态系统服务价值当量因子;合成模式本地生成,无外部数据源。

隐私声明 / Privacy

  • 默认完全离线运行,不发起任何网络请求
  • --synthetic 模式不读取任何外部数据
  • 所有计算在本地完成,不上传用户数据

License

MIT

Related skills

多压力因子归一化加权叠加得综合影响指数,独立概率模型估算累积效应,按阈值划分5 级影响等级。Assesses environmental impact grades by multi-factor overlay and cumulative effects. 输出影响指数与等级 GeoTIFF。

基于生境异质性假说,用 NDVI 生产力、纹理结构异质性与地形粗糙度三类代理估算物种丰富度空间分布。Maps species richness proxies from NDVI, texture and terrain heterogeneity. 输出物种丰富度与生境质量 GeoTIFF + 参数 JSON。

由 NDVI 幂律异速生长方程估算地上生物量碳,叠加根茎比地下碳与类型化土壤碳密度。Estimates carbon stocks from biomass allometry and soil carbon density. 输出地上碳/土壤碳/总碳三张 GeoTIFF 与汇总 JSON。

简化 InVEST Budyko 水量平衡计算产水量,叠加 NDVI 调制的植被截留系数得水源涵养量,并估算养分截留净化量。Maps water retention and purification with a simplified InVEST water yield model. 输出产水/涵养/净化三张 GeoTIFF。

1 installs

Estimate forest carbon stock from remote sensing data using BEF, allometric equations, or IPCC Tier 1/2 methods. Includes Monte Carlo uncertainty analysis. Supports raster (GeoTIFF) and tabular (CSV) inputs.