文档

Geoskill: Land Subsidence InSAR

试用

Analyze land subsidence from InSAR displacement data. Use when the user wants to analyze changes, detect hazards, or generate assessment reports.

它能做什么

Analyze land subsidence from InSAR displacement data. Use when the user wants to analyze changes, detect hazards, or generate assessment reports.

技能文档

Land Subsidence InSAR

Analyze land subsidence from InSAR displacement data.

CLI Usage

python scripts/land_subsidence_insar.py --displacement disp.tif
python scripts/land_subsidence_insar.py --displacement disp.tif --reference baseline.tif
python scripts/land_subsidence_insar.py --displacement disp.tif --output-dir my_output

Data Download

This skill can auto-fetch a Sentinel-1 GRD scene from the Microsoft Planetary Computer when given a bounding box + date range. The downloaded asset is then passed through --displacement (note: real InSAR displacement products require additional SAR processing; the raw backscatter/intensity is only a rough proxy for subsidence hotspots — see Limitations below).

# Fetch a S1 GRD scene for the Beijing area and analyze
python scripts/land_subsidence_insar.py \
  --bbox 116.4,39.9,116.42,39.92 \
  --date-range 2024-06-01,2024-06-30 \
  --output-dir ./subsidence-output

# Or via an AOI polygon file
python scripts/land_subsidence_insar.py \
  --aoi-file ./aoi.geojson \
  --output-dir ./subsidence-output

The PYTHONPATH must include the parent of _geoskill_data_fetcher/ (the same directory the 50 skills live in). Set it once:

export PYTHONPATH="/path/to/行业Skill创意-20260727"

Parameters

ArgumentRequiredDefaultDescription
--displacementYes*Path to InSAR displacement raster (mm, positive = uplift, negative = subsidence)
--referenceNoOptional reference displacement raster for differential analysis
--bboxNoBounding box W,S,E,N — auto-downloads S1 GRD from MPC
--date-rangeNoDate range START,END (ISO-8601) for auto-download
--aoi-fileNoPath to AOI GeoJSON polygon for auto-download
--output-dir, -oNosubsidence-outputDirectory to write outputs

*Required unless --synthetic or --bbox+--date-range is supplied.

Output

FileDescription
subsidence-report.jsonMachine-readable subsidence stats (max/min/mean rate, affected area)
report.htmlHuman-readable HTML report
output-manifest.jsonRun metadata + result summary

Exit Codes

CodeMeaning
0Success
2Argument error
7Processing failure

相关技能

Extract flood extent from SAR backscatter imagery. Use when the user wants to analyze changes, detect hazards, or generate assessment reports.

1 次安装

Multi-criteria suitability analysis for infrastructure site selection. Use when the user wants to analyze changes, compare multi-temporal rasters, compute indices, or generate assessment reports.

1 次安装

融合InSAR形变速率、后向散射变化与DEM坡度综合加权评分,双门限提取疑似滑坡连通域并矢量化分级(high/medium/low),输出滑坡GeoJSON、形变速率/风险评分GeoTIFF与风险汇总JSON。SAR landslide detection fusing InSAR deformation, backscatter change and slope.

1 次安装

Detect forest disturbance from multi-temporal NDVI. Use when the user wants to analyze changes, detect hazards, or generate assessment reports.

1 次安装

Estimate impervious surface fraction from multi-band satellite imagery (Sentinel-2) using spectral indices (NDBI, NDVI, MNDWI). Supports binary classification and continuous fraction estimation, with zone-level aggregation and change detection. Use when mapping urban impervious surfaces, computing impervious ratios by watershed/admin unit, or analyzing temporal changes in built-up areas.

1 次安装