2022 개정 교육과정 연계 온디바이스 기상관측·원격탐사 STEM 시리즈 #75

이중편파 기상 레이더 수상체(Hydrometeor) 위상차 및 강수 분류 서보 엔진

비차등 위상차($K_{dp}$) 미분 필터링 및 0-Alloc C# 우박·돌발홍수 실시간 수상체 식별 제어 루프

ZeliDesk 기상레이더·원격탐사팀 국방기상정보·악기상조기경보 분과

2026년 9월 12일

1. 이중편파 레이더 관측 변수 및 전파 위상 지연 모델

수평 반사도 인자 $Z_h$와 수직 반사도 인자 $Z_v$의 로그비인 차등 반사도 $Z_{dr}$는 강수 입자의 평균 편평도(축비)를 나타냅니다:

편평한 강수 입자를 통과할 때 수평 편파가 더 큰 유전체 매질 경로를 겪어 수직 편파 대비 위상 지연이 발생하며, 경로 거리 $r$에 따른 비차등 위상차 $K_{dp}$는 총 위상차 $\Psi_{dp}$의 공간 도함수로 정의됩니다:

$K_{dp}$는 부분 빔 차폐(Beam Blockage)나 레이더 감쇠의 영향을 받지 않아 집중호우 시 강우 강도($R$)를 정확히 추정($R = a K_{dp}^b$)하는 핵심 지표이며, 교차상관계수 $\rho_{hv}$는 입자의 균일성을 평가하여 비기상 에코(새, 곤충, 지형 클러터)를 $0.8$ 이하로 즉각 필터링합니다.

2. 0-Alloc C# 스택 기반 기상 레이더 수상체 분류 서보

게이트 단위 편파 관측치로부터 노이즈를 제거하고 퍼지 소속도 함수 기반으로 수상체 유형(비, 눈, 우박, 클러터)을 실시간 판별하는 ref struct 엔진입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.MeteorologicalRadar
{
    public enum HydrometeorType : byte
    {
        ClearOrClutter = 0,
        LightRain = 1,
        HeavyRain = 2,
        Hail = 3,
        WetSnow = 4
    }

    public readonly ref struct DualPolGateTelemetry
    {
        public readonly double ReflectivityZhDbz;
        public readonly double DifferentialReflectivityZdrDb;
        public readonly double DifferentialPhasePsiDpDeg;
        public readonly double CorrelationRhoHv;
        public readonly double GateRangeKm;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public DualPolGateTelemetry(double zh, double zdr, double psiDp, double rhoHv, double rangeKm)
        {
            ReflectivityZhDbz = zh;
            DifferentialReflectivityZdrDb = zdr;
            DifferentialPhasePsiDpDeg = psiDp;
            CorrelationRhoHv = rhoHv;
            GateRangeKm = rangeKm;
        }
    }

    public readonly ref struct RadarClassificationResult
    {
        public readonly double FilteredKdpDegPerKm;
        public readonly double RainRateMmPerHour;
        public readonly HydrometeorType IdentifiedType;
        public readonly bool IsSevereWeatherAlert;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public RadarClassificationResult(double kdp, double rainRate, HydrometeorType type, bool alert)
        {
            FilteredKdpDegPerKm = kdp;
            RainRateMmPerHour = rainRate;
            IdentifiedType = type;
            IsSevereWeatherAlert = alert;
        }
    }

    public ref struct DualPolClassificationServo
    {
        private double _prevPsiDp;
        private double _prevRangeKm;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public DualPolClassificationServo(double initialRangeKm)
        {
            _prevPsiDp = 0.0;
            _prevRangeKm = initialRangeKm;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public RadarClassificationResult Step(in DualPolGateTelemetry telemetry)
        {
            // 1. 공간 유한 차분을 통한 Kdp(비차등 위상차) 수치 미분
            double deltaR = telemetry.GateRangeKm - _prevRangeKm;
            double rawKdp = 0.0;
            if (deltaR > 1e-4)
            {
                rawKdp = 0.5 * (telemetry.DifferentialPhasePsiDpDeg - _prevPsiDp) / deltaR;
            }
            _prevPsiDp = telemetry.DifferentialPhasePsiDpDeg;
            _prevRangeKm = telemetry.GateRangeKm;

            // Kdp 음수 클리핑 (물리적 증가량 보존)
            double filteredKdp = Math.Max(0.0, rawKdp);

            // 2. 비기상 에코(Non-meteorological Clutter) 필터링
            if (telemetry.CorrelationRhoHv < 0.80 || telemetry.ReflectivityZhDbz < 5.0)
            {
                return new RadarClassificationResult(filteredKdp, 0.0, HydrometeorType.ClearOrClutter, false);
            }

            // 3. 퍼지 수상체 판별 로직
            HydrometeorType type;
            bool alert = false;

            if (telemetry.ReflectivityZhDbz > 52.0 && telemetry.DifferentialReflectivityZdrDb < 1.2)
            {
                // 강한 반사도이나 입자가 둥글어 Zdr이 낮음 -> 우박(Hail)
                type = HydrometeorType.Hail;
                alert = true;
            }
            else if (telemetry.ReflectivityZhDbz > 40.0 || filteredKdp > 1.5)
            {
                // 높은 반사도 및 큰 Kdp -> 집중호우(Heavy Rain)
                type = HydrometeorType.HeavyRain;
                alert = true;
            }
            else if (telemetry.ReflectivityZhDbz > 20.0 && telemetry.DifferentialReflectivityZdrDb > 0.5)
            {
                type = HydrometeorType.LightRain;
            }
            else
            {
                type = HydrometeorType.WetSnow;
            }

            // 4. 강우 강도 추정 (Kdp 기반 관계식: R = 40.5 * Kdp^0.85)
            double rainRate = 0.0;
            if (filteredKdp > 0.1)
            {
                rainRate = 40.5 * Math.Pow(filteredKdp, 0.85);
            }
            else if (telemetry.ReflectivityZhDbz > 10.0)
            {
                // 마셜-파머 Z-R 관계식: Z = 200 * R^1.6
                double zLinear = Math.Pow(10.0, telemetry.ReflectivityZhDbz / 10.0);
                rainRate = Math.Pow(zLinear / 200.0, 1.0 / 1.6);
            }

            return new RadarClassificationResult(filteredKdp, rainRate, type, alert);
        }
    }
}

3. 교육적 가치 및 중고교 실험 시뮬레이션 적용