=== versions ===
python 3.12.0 | rasterio 1.5.2 | numpy 2.5.4 | scipy 1.18.1 | h5py 3.16.0 | matplotlib 3.11.2

=== LiCSAR interferograms ===
frame 143D_07132_131313: 4,434 interferograms online, 20141003 to 20260913
selected: 170 pairs of <= 18 days between 61 epochs, 20250911 to 20260913
temporal baselines (days): {5: 1, 6: 54, 7: 2, 12: 56, 13: 2, 18: 55}
heading -169.2°, mean incidence 33.9°

=== stack ===
170 interferograms x 400 x 500 pixels (3.6 arcsec, about 105 m x 111 m)
unwrapped in every pair: 94.3% of pixels; LOS up component 0.695..0.733 (local incidence 44.5°)

=== single pair 20260103_20260304 (60 days) ===
wrapped phase range -3.14..3.14 rad; unwrapped 1.2..16.8 rad
one fringe (2π) = 27.7 mm of LOS motion
LOS change relative to CU: 1st percentile -33 mm, 99th percentile 16 mm
median coherence 0.52; pixels with coherence < 0.3: 25.2%

=== coherence mask ===
mean coherence: median 0.59; reference area 0.61
kept 150,831 of 200,000 pixels (75.4%) with mean coherence >= 0.35 and unwrapped in all 170 pairs

=== SBAS inversion ===
design matrix (170, 60), rank 60 (needs 60: one connected network)
residual RMS per interferogram: median 0.15 mm, max 0.39 mm (20251011_20251029)
pixels with a residual > π in any pair (likely unwrapping errors): 0.01%

=== velocity ===
LOS velocity: 1st percentile -215 mm/yr, median -8 mm/yr, 99th percentile 9 mm/yr
vertical (LOS / up component): fastest -383 mm/yr; area subsiding faster than 100 mm/yr: 367 km², faster than 200 mm/yr: 136 km²
median 1-sigma velocity uncertainty (scatter about the line): 3.4 mm/yr
           place  height_m  coherence  los_mm_yr  vertical_mm_yr  los_total_mm
          Chalco      2242       0.70     -113.9          -156.9         -96.0
         Tláhuac      2238       0.72     -156.6          -218.0        -145.4
  Nezahualcóyotl      2236       0.74     -199.7          -278.7        -195.4
  Airport (AICM)      2226       0.35     -170.1          -238.8        -176.1
Centro Histórico      2242       0.72      -87.5          -123.7         -89.2
  Reference (CU)      2290       0.65        0.1             0.1           0.2

=== profile A-B ===
length 39.1 km; valid samples 81%; fastest -260 mm/yr at 28.6 km; terrain 2234-2563 m

=== outputs ===
wrote vertical_velocity_mm_yr.tif, los_timeseries_points.csv, data/ifgramStack.h5 and six figures

=== MintPy cross-check (python mintpy_crosscheck.py) ===
MintPy 1.6.4

MintPy time series: 61 dates, 20250911 to 20260913
LOS velocity, MintPy: 1st percentile -214 mm/yr, median -7 mm/yr
MintPy minus NumPy over 150,831 pixels: mean +0.62 mm/yr, max |difference| 0.83 mm/yr, correlation 1.0000