生态系统服务评估 | 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 / 输出

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

输出

数据源 / Source

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

隐私声明 / Privacy

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

License

MIT