Integrate terrain, geology, rainfall, land cover, roads, and historical landslide data to produce interpretable susceptibility zoning with spatial cross-validation.
1 installs
Documents
Geoskill: Land Subsidence InSAR
Try itAnalyze land subsidence from InSAR displacement data. Use when the user wants to analyze changes, detect hazards, or generate assessment reports.
What it does
Analyze land subsidence from InSAR displacement data. Use when the user wants to analyze changes, detect hazards, or generate assessment reports.
The skill document
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
| Argument | Required | Default | Description |
|---|---|---|---|
--displacement | Yes* | — | Path to InSAR displacement raster (mm, positive = uplift, negative = subsidence) |
--reference | No | — | Optional reference displacement raster for differential analysis |
--bbox | No | — | Bounding box W,S,E,N — auto-downloads S1 GRD from MPC |
--date-range | No | — | Date range START,END (ISO-8601) for auto-download |
--aoi-file | No | — | Path to AOI GeoJSON polygon for auto-download |
--output-dir, -o | No | subsidence-output | Directory to write outputs |
*Required unless --synthetic or --bbox+--date-range is supplied.
Output
| File | Description |
|---|---|
subsidence-report.json | Machine-readable subsidence stats (max/min/mean rate, affected area) |
report.html | Human-readable HTML report |
output-manifest.json | Run metadata + result summary |
Exit Codes
| Code | Meaning |
|---|---|
| 0 | Success |
| 2 | Argument error |
| 7 | Processing failure |
Related skills
Extract flood extent from SAR backscatter imagery. Use when the user wants to analyze changes, detect hazards, or generate assessment reports.
1 installs
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 installs
融合InSAR形变速率、后向散射变化与DEM坡度综合加权评分，双门限提取疑似滑坡连通域并矢量化分级（high/medium/low），输出滑坡GeoJSON、形变速率/风险评分GeoTIFF与风险汇总JSON。SAR landslide detection fusing InSAR deformation, backscatter change and slope.
1 installs
Detect forest disturbance from multi-temporal NDVI. Use when the user wants to analyze changes, detect hazards, or generate assessment reports.
1 installs
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 installs