Treffer: Rouse-alpha Inversion & Transect Vertical Profile Toolkit

This repository provides two Python tools for suspended-sediment analysis using the Rouse framework. RouseAlphaSolver.py — Rouse-α Inversion (Map & Calculator modes) Map mode: Inverts the Rouse number (P) pixelwise from a depth raster (H) using a target α (alpha), producing a georeferenced P map (GeoTIFF). Calc mode: Computes P directly for user-supplied depth values h (list or CSV). Important: target_alpha = 1.27 is a global-dataset-derived default (prior). α varies by river/reach, sediment (incl. flocculation), and hydraulics. For scientific use, calibrate/validate α for each study area. RouseAlphaVerticalEstimation.py — Vertical SSC along a Transect from a Mapped Rouse Number Samples P (GeoTIFF), reference/surface SSC (GeoTIFF), and depth H (GeoTIFF) along a user line (shp/geojson/gpkg), then reconstructs vertical SSC profiles using the classical Rouse power-law:C(z) = C_ref * ((z_ref / z) * ((h − z) / (h − z_ref)))^Pwith configurable reference level (default z_ref = h − 0.5 m). Outputs point samples (distance, P, SSC_ref, h) and a vertical SSC matrix (layers × points) as .npy and optional CSV, plus an optional quicklook plot. What’s included RouseAlphaSolver.py — CLI (mode=map|calc), ENVI-based georeferencing, masks, robust solver (scipy.optimize.fsolve), plotting. RouseAlphaVerticalEstimation.py — CRS-aware transect sampling (rasterio), optional median filtering, vertical grid control, standardized outputs. Quick start Build a Rouse-number map (GeoTIFF)RouseAlphaSolver.py --mode map --ndwi NDWI.npy --depth H_map.npy --envi_dat subset_hsi --out_tif r_map_alpha127.tif --reduction_factor 3 --plot Direct calculatorpython RouseAlphaSolver.py --mode calc --h 0.8 1.0 2.5 5.0 --out_csv P_values.csv Transect-based vertical SSC from mapped Ppython RouseAlphaVerticalEstimation.py --line transect.shp --spacing 2 --rmap r_map_alpha127.tif --ssc TSS_map.tif --depth h_map.tif --layers 22 --ref_offset 0.5 --out_points transect_points.gpkg --out_npy vertical_ssc.npy --out_csv vertical_ssc.csv --plot Inputs & outputs ...