AIDT 온디바이스 STEM 위성 영상 레이더 간섭계(InSAR) 위상 언래핑 및 대기 지연 보정 서보 엔진

1. 개요: 우주 레이더 간섭계(InSAR)와 지각 변동의 극초정밀 계측

합성개구레이더 간섭계(Interferometric Synthetic Aperture Radar, InSAR)는 동일한 지표 영역을 서로 다른 시각 또는 궤도에서 관측한 두 복소 SAR 영상의 위상차(\(\Delta \phi\))를 계산하여 밀리미터(\(\text{mm}\)) 단위의 지표 변위와 지형 고도를 원격 탐사하는 우주 지구물리 기술입니다. 화산 마그마 융기, 지진 단층 활동 및 지반 침하를 모니터링할 때 획득되는 위상은 삼각함수 주기성에 의해 \([-\pi, \pi)\) 구간으로 랩핑(Phase Wrapping)되어 있으므로, 인접 픽셀 간의 정수 \(2\pi k\) 불연속을 복원하는 2차원 위상 언래핑(Phase Unwrapping)과 기상 조건에 따른 대류권 수증기 지연(\(\phi_{atm}\)) 필터링이 필수적입니다.

2. 간섭 위상 성분 분해 및 위상 언래핑 수식 모델

두 관측 영상 사이의 총 간섭 위상 \(\Delta \phi\)는 지형 위상, 지표 변위, 대기 전파 지연 및 노이즈로 분해됩니다:

\Delta \phi = \phi_{topo} + \phi_{def} + \phi_{atm} + \phi_{noise} \pmod{2\pi}

외부 수치표고모델(DEM)로 지형 위상을 제거한 차분 간섭계(D-InSAR)에서 레이더 시선 방향(LOS, Line-of-Sight) 지표 변위 \(\Delta r\)은 다음과 같이 변환됩니다:

\Delta r = \frac{\lambda}{4\pi} \phi_{def}

여기서 \(\lambda\)는 마이크로파 파장(예: C-band 5.6 cm, X-band 3.1 cm)입니다. 이토(Itoh) 조건(\(|\nabla \phi_{true}| < \pi\)) 하에서 비랩핑 위상 구배는 랩핑 연산자 \(\mathcal{W}\{\cdot\}\)를 통해 연속 복원됩니다:

\nabla \psi = \mathcal{W}\{ \nabla \phi_{wrapped} \} \implies \nabla^2 \psi = \nabla \cdot \mathcal{W}\{ \nabla \phi_{wrapped} \}

3. 0-Alloc C# InSAR 위상 언래핑 서보 엔진 구현

국소 잔차(Residue)를 감지하고 대기 위상 지터를 10,000Hz 주기로 제거하는 ref struct 아키텍처입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.Geophysics
{
    public readonly ref struct InSarDeformationState
    {
        public readonly float TrueDisplacementMm;
        public readonly float UnwrappedPhaseRad;
        public readonly float AtmosphericDelayMm;
        public readonly float Coherence;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public InSarDeformationState(float dispMm, float unwrappedRad, float atmMm, float coh)
        {
            TrueDisplacementMm = dispMm;
            UnwrappedPhaseRad = unwrappedRad;
            AtmosphericDelayMm = atmMm;
            Coherence = coh;
        }
    }

    public ref struct SatelliteInSarServo
    {
        private const float RadarWavelengthMm = 55.46f; // Sentinel-1 C-band (5.55 cm)
        private const float RadToMm = RadarWavelengthMm / (4.0f * MathF.PI);

        private float _previousWrappedPhase;
        private int _accumulatedCycleK;
        private float _estimatedAtmosphericMm;

        public SatelliteInSarServo(float initialPhase = 0.0f)
        {
            _previousWrappedPhase = initialPhase;
            _accumulatedCycleK = 0;
            _estimatedAtmosphericMm = 0.0f;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public InSarDeformationState UpdateServo(float currentWrappedPhase, float atmTempWaterVaporIndex, float dt = 0.0001f)
        {
            // 1D Itoh unwrapping differentiator
            float phaseDiff = currentWrappedPhase - _previousWrappedPhase;
            if (phaseDiff > MathF.PI)
            {
                _accumulatedCycleK--;
            }
            else if (phaseDiff < -MathF.PI)
            {
                _accumulatedCycleK++;
            }
            _previousWrappedPhase = currentWrappedPhase;

            float unwrappedPhase = currentWrappedPhase + (_accumulatedCycleK * 2.0f * MathF.PI);

            // Tropospheric water vapor delay estimation (low-pass spatial filter proxy)
            float targetAtmMm = atmTempWaterVaporIndex * 2.8f;
            _estimatedAtmosphericMm += (targetAtmMm - _estimatedAtmosphericMm) * (dt * 20.0f);

            // Corrected displacement in mm
            float rawDisplacementMm = unwrappedPhase * RadToMm;
            float correctedDisplacementMm = rawDisplacementMm - _estimatedAtmosphericMm;

            float coherence = MathF.Clamp(1.0f - MathF.Abs(phaseDiff) / (2.0f * MathF.PI), 0.1f, 1.0f);

            return new InSarDeformationState(correctedDisplacementMm, unwrappedPhase, _estimatedAtmosphericMm, coherence);
        }
    }
}

4. 실시간 위상 랩핑/언래핑 및 지반 침하 계측 인터랙티브 시뮬레이터

지표 침하량과 대기 수증기 지연을 조절하며 무지개빛 랩핑 간섭 줄무늬(Fringes)가 연속 변위 곡선으로 복원되는 과정을 확인하십시오.

지표 침하 중심 변위 (Peak Subsidence): 28 mm 대기 수증기 요동 (Atmospheric Noise): 5 mm

5. 교육 연계 및 국방/재난 방재 실무적 의의

2022 개정 교육과정 지구과학 I/II: 판 구조론과 지각 변동, 화산 폭발 전조 분화구 융기, 지반 침하 및 지구 관측 위성 기술.
물리학 II: 전자기파 마이크로파 간섭, 위상차와 경로차의 기하학적 관계, 위상 모듈로(Modulo) 불연속 복원.
국방 위성 정찰 / 국토 안전: 적 군사 기지 지하 갱도 굴착 탐지, 댐/교량/활주로 밀리미터 안전 진단.