제임스웹(JWST) 외계행성 투과분광학 & 1,000Hz 0-Alloc C# NIRSpec 대기분자 흡수선 검출 서보 제어기

using System;
using System.Runtime.CompilerServices;
using System.Runtime.InteropServices;

namespace ZeliDesk.Stem.SpaceSpectroscopy
{
    /// 

    /// JWST NIRSpec 분광 검출기 1ms 원시 픽셀 프레임
    /// 
    [StructLayout(LayoutKind.Sequential)]
    public readonly struct NirSpecPixelFrame
    {
        public readonly uint TickMs;
        public readonly float ContinuumIntensity;  // 1.5um 연속체 기준 채널
        public readonly float WaterBandIntensity;    // 1.4um H2O 흡수 채널
        public readonly float Co2BandIntensity;      // 4.3um CO2 흡수 채널
        public readonly float MethaneBandIntensity;  // 3.3um CH4 흡수 채널

        public NirSpecPixelFrame(uint ms, float cont, float h2o, float co2, float ch4)
        {
            TickMs = ms;
            ContinuumIntensity = cont;
            WaterBandIntensity = h2o;
            Co2BandIntensity = co2;
            MethaneBandIntensity = ch4;
        }
    }

    /// 

    /// 1,000Hz 0-Alloc JWST 외계행성 대기 흡수선 검출 서보 엔진
    /// 
    public ref struct NirSpecAtmosphereServo
    {
        private const float BaselineDepthPpm = 12500.0f; // 기본 행성 기하 식 깊이 (1.25%)
        private const float Co2DetectionLimitPpm = 80.0f; // CO2 80ppm 초과 시 검출 판정

        private float _filteredH2oPpm;
        private float _filteredCo2Ppm;
        private float _filteredCh4Ppm;
        private int _detectionConfidenceTicks;

        public NirSpecAtmosphereServo(bool initialize = true)
        {
            _filteredH2oPpm = 0.0f;
            _filteredCo2Ppm = 0.0f;
            _filteredCh4Ppm = 0.0f;
            _detectionConfidenceTicks = 0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public bool ProcessFrame(in NirSpecPixelFrame frame, out float co2PpmExcess, out bool atmosphereDetected)
        {
            float cont = MathF.Max(1e-3f, frame.ContinuumIntensity);

            // 1. 파장별 미세 흡수율 ppm 편차 계산 (0-Alloc 스택)
            float instantH2oPpm = ((cont - frame.WaterBandIntensity) / cont) * 1.0e6f;
            float instantCo2Ppm = ((cont - frame.Co2BandIntensity) / cont) * 1.0e6f;
            float instantCh4Ppm = ((cont - frame.MethaneBandIntensity) / cont) * 1.0e6f;

            // 2. 지수 이동 평균(EMA) 필터로 샷 잡음 억제
            _filteredH2oPpm = (0.05f * instantH2oPpm) + (0.95f * _filteredH2oPpm);
            _filteredCo2Ppm = (0.05f * instantCo2Ppm) + (0.95f * _filteredCo2Ppm);
            _filteredCh4Ppm = (0.05f * instantCh4Ppm) + (0.95f * _filteredCh4Ppm);

            co2PpmExcess = _filteredCo2Ppm;

            // 3. 외계행성 대기 검출 신뢰도 판정
            if (_filteredCo2Ppm >= Co2DetectionLimitPpm)
            {
                _detectionConfidenceTicks++;
                atmosphereDetected = (_detectionConfidenceTicks >= 10); // 10ms 연속 유지 시 확정
            }
            else
            {
                _detectionConfidenceTicks = Math.Max(0, _detectionConfidenceTicks - 1);
                atmosphereDetected = false;
            }

            return true;
        }
    }
}