InSAR 변위 탐지

PSI · MiaplPy

PSphase linkingwiki/psinsar-miaplpy-workflow

한 줄 요약

PSInSAR(PSI)는 시간이 지나도 위상이 안정적인 영구산란체(Persistent Scatterer)만 점으로 골라 풀해상도로 mm/yr 변위를 뽑는 기법으로, 인접 픽셀을 멀티룩 평균하는 SBAS와는 정반대 철학이다. MiaplPy(ISCE2 stack 위에서 동작)로 죽전 도심 PSI를 끝까지 완주해 108,690점을 대시보드에 적재했다. 가장 값비싼 phase_linking 단계는 BLAS 스레드 오버서브스크립션을 잡아 5시간에서 6분으로(~15배) 줄였다.

PSI vs SBAS — 정반대 철학

둘 다 같은 coregistered SLC 스택에서 갈리는데, SBAS는 인접 픽셀을 멀티룩 평균해 면적 변위를 얻는 분산 산란체(Distributed Scatterer) 접근이다.
PSI는 위상이 안정적인 점만 선별해 풀해상도를 유지하므로 건물·인공구조물처럼 도심에서 강하다.
도심 정밀은 PSI, 광역 면적은 거친 SBAS가 정석이며, 매립지·수변 간척지도 PSI 영역이다 — 송도에서 PSI 29,929점 성공 vs SBAS는 temporal coherence 0.17로 적재 불가였다.

	SBAS	PSI
산란체	Distributed (멀티룩 평균)	Persistent (점 선별)
해상도	멀티룩만큼 희생 (10–55m)	풀해상도 유지
강한 지형	면적·평지·(거칠게) 산지	도심·인공구조물

같은 스택에서 대상·해상도·철학이 갈린다

MiaplPy 파이프라인

입력은 SBAS와 달리 멀티룩하지 않은 풀해상도 SLC이며, r6 도심 coregistration 결과(coreg_secondarys 87장)를 재사용했다.
merged의 SLC/geom이 VRT(가상)만 존재해, MiaplPy가 요구하는 물리 바이너리(.full)로 materialize하는 전처리가 필요했다(~74GB).
마지막 timeseries_correction에서 ERA5·demErr 보정 후 velocity로 변환하고 geocode한다.

📥load_dataAOI subset 크롭

🔗phase_linking최대 비용 단계

🧩ifgram·unwrap

🔄invert_network

🗺timeseries_correctionERA5·demErr → velocity → geocode

load_data부터 timeseries_correction까지 — phase_linking이 최대 비용

phase_linking 스레드 함정 — 5시간→6분

PSI 처리에서 가장 값비싼 단계로, 기본 num_worker=4에 워커당 ~16 BLAS 스레드가 64코어에 과다 배정되면 thread-barrier spin으로 us=100%지만 헛돈다(총 ~5시간).
87×87 EMI 소행렬은 BLAS 멀티스레드 스케일링이 나빠, 스레드를 3으로 묶고 프로세스 병렬 numProcessor를 패치 수 이상으로 올리는 게 정답이다.
그 결과 21개 패치가 1-wave로 병렬 처리되어 patch당 3~5분·총 ~6분이 됐고, flag.npy가 있는 patch는 자동 skip이라 중단·재시작도 안전하다.

QA 결과 — 죽전 도심 완주

레이더좌표 609×1393에서 진짜 PS(maskPS)는 848,337점 중 56,975점(6.7%)이 선별됐다.
temporal coherence는 mean 0.485, ≥0.7이 25.6%로 도시 영역에서 정상 수준이며, 보정 LOS velocity는 mean +0.20·std 0.94 mm/yr였다.
최종적으로 stack_id=3 'jukjeon_psi_r6'으로 108,690점(10m 재-geocode)이 대시보드에 적재됐다.

주의점 · 함정

대시보드는 SBAS의 unwrap 검증 지표(conncomp_ok)로 '신뢰점'을 세는데, PSI엔 connComp가 없어 NULL이 되므로 '신뢰 점만' 토글을 켜면 0점이 된다 — 거짓으로 conncomp_ok=true를 박는 건 부정직하고, 정직한 해법은 대시보드가 temporalCoherence 신뢰등급을 PSI stack에 별도 지원하는 것이다. 또한 108,690 적재점도 단일 PS로 인용하지 말고 건물 단위(n_PS≥10)·단지·AOI 평균으로만 봐야 한다.