Coding

geoskill-light-pollution-assessment

Try it

由 VIIRS 夜光辐射值按生态阈值分 0-5 级光污染等级,对数响应模型估算生态影响指数,附天空辉光代理。Assesses light pollution grades from VIIRS night-time lights and ecological thresholds. 输出等级/生态影响/天空辉光三张 GeoTIFF。

What it does

由 VIIRS 夜光辐射值按生态阈值分 0-5 级光污染等级,对数响应模型估算生态影响指数,附天空辉光代理。Assesses light pollution grades from VIIRS night-time lights and ecological thresholds. 输出等级/生态影响/天空辉光三张 GeoTIFF。

The skill document

光污染评估 | Light Pollution Assessment

Light pollution grade thresholds (0.25/1/4/15/50 nW·cm⁻²·sr⁻¹) follow the global light pollution classes of Falchi et al. (2016): 0 = pristine dark sky, 5 = extreme light pollution. The ecological impact index uses a logarithmic response I = log10(1+k·R)/log10(1+k·Rmax) normalized to [0,1], characterizing the disruption of artificial light on nocturnal organism rhythms; the skyglow proxy = radiance × scattering coefficient.

Use cases: dark-sky reserve delineation, ecological light-environment assessment, and urban lighting planning.

Dependencies / 依赖

pip install numpy rasterio

Usage / 使用方法

Example 1: synthetic urban radiance gradient

python geoskill-light-pollution-assessment.py --bbox 116.0 39.0 117.0 40.0 --synthetic --output-dir ./output

Example 2: real VIIRS annual mean radiance raster

python geoskill-light-pollution-assessment.py --input viirs_annual.tif --output-dir ./real

Example 3: a different region

python geoskill-light-pollution-assessment.py --bbox 121 31 122 32 --synthetic --output-dir ./shanghai

Example 4: quick validation on a tiny area

python geoskill-light-pollution-assessment.py --bbox 116.39 39.90 116.40 39.91 --synthetic --output-dir ./tiny

Example 5: silent batch run

python geoskill-light-pollution-assessment.py --bbox 113 23 114 24 --synthetic --quiet --output-dir ./batch

Output / 输出

FileFormatDescription
light_pollution_grade.tifGeoTIFF (float32)Light pollution grade 0-5
ecological_impact.tifGeoTIFF (float32)Ecological impact index ∈ [0,1]
skyglow_proxy.tifGeoTIFF (float32)Skyglow proxy
light_pollution_params.jsonJSONThreshold table, per-grade pixel counts, radiance statistics
output-manifest.jsonJSONRun manifest (input/output/QA/software versions)

Data Source / 数据源 / Source

Local VIIRS night-time lights GeoTIFF (nW·cm⁻²·sr⁻¹); grading thresholds follow Falchi et al. 2016 (Science Advances, open access); synthetic mode generates an urban radiance gradient locally with no external data source.

Privacy / 隐私声明 / Privacy

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

License / License

MIT



name: geoskill-light-pollution-assessment description: '由 VIIRS 夜光辐射值按生态阈值分 0-5 级光污染等级,对数响应模型估算生态影响指数,附天空辉光代理。Assesses light pollution grades from VIIRS night-time lights and ecological thresholds. 输出等级/生态影响/天空辉光三张 GeoTIFF。'

光污染评估 | Light Pollution Assessment

光污染等级阈值(0.25/1/4/15/50 nW·cm⁻²·sr⁻¹)参考 Falchi et al. (2016) 全球光污染分级:0=原始暗夜,5=极端光污染。生态影响指数用对数响应 I = log10(1+k·R)/log10(1+k·Rmax) 归一化到 [0,1],刻画人造光对夜行生物节律的干扰;天空辉光代理 = 辐射值×散射系数。

适用场景:暗夜保护区划定、生态光环境评估、城市照明规划。

依赖

pip install numpy rasterio

使用方法

示例 1:合成城市梯度夜光

python geoskill-light-pollution-assessment.py --bbox 116.0 39.0 117.0 40.0 --synthetic --output-dir ./output

示例 2:真实 VIIRS 年平均辐射栅格

python geoskill-light-pollution-assessment.py --input viirs_annual.tif --output-dir ./real

示例 3:不同区域

python geoskill-light-pollution-assessment.py --bbox 121 31 122 32 --synthetic --output-dir ./shanghai

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

python geoskill-light-pollution-assessment.py --bbox 116.39 39.90 116.40 39.91 --synthetic --output-dir ./tiny

示例 5:静默批量

python geoskill-light-pollution-assessment.py --bbox 113 23 114 24 --synthetic --quiet --output-dir ./batch

输出

文件格式说明
light_pollution_grade.tifGeoTIFF (float32)光污染等级 0-5
ecological_impact.tifGeoTIFF (float32)生态影响指数 ∈ [0,1]
skyglow_proxy.tifGeoTIFF (float32)天空辉光代理
light_pollution_params.jsonJSON阈值表、等级像元计数、辐射统计
output-manifest.jsonJSON运行清单(输入/输出/QA/软件版本)

数据源 / Source

本地 VIIRS 夜光 GeoTIFF(nW·cm⁻²·sr⁻¹);分级阈值参考 Falchi et al. 2016(Science Advances, 开放获取);合成模式本地生成城市梯度夜光,无外部数据源。

隐私声明 / Privacy

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

License

MIT

Related skills

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

评估多光谱影像辐射质量(SNR/条纹/坏线)与几何质量(云量/清晰度),输出质量评分 JSON + HTML 报告。Assess radiometric (SNR, striping, dead lines) and geometric (cloud, sharpness) quality of multispectral imagery; outputs a scored JSON and HTML report.

1 installs

交通噪声衰减模型制图 dB(A) 等级:FHWA 简化声源级、点/线源几何发散、建筑屏障衰减与地面吸收。Maps traffic noise levels with distance attenuation, building shielding and terrain. 输出噪声级 GeoTIFF + 参数 JSON。

1 installs

融合 NDVI 趋势(Sen 斜率/线性回归)、反照率与植被稀缺度,综合评分分级荒漠化(稳定/轻/中/重度),输出等级与趋势栅格、面积统计。Fuses NDVI trend (Sen/linear slope), albedo and vegetation scarcity to grade desertification.

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

Calculate solar PV energy potential from NASA POWER solar radiation data. Computes annual GHI, optimal tilt angle, estimated PV output, and economic analysis.