AIDT 온디바이스 STEM 초음파 전단파 탄성 영상(SWE) 점탄성 전파 서보 엔진

1. 개요

전단파 탄성초음파 영상(Shear Wave Elastography, SWE)은 음향 복사력 임펄스(Acoustic Radiation Force Impulse, ARFI)를 생체 조직 심부에 밀리초 단위로 집속 조사하여 조직 국소 변위를 일으킨 후, 방사형으로 퍼져나가는 횡파(전단파, Shear Wave)의 전파 궤적을 초당 수천 프레임(> 5,000 fps)의 초고속 평면파(Plane Wave) 초음파로 추적하여 조직의 탄성 계수(Young's Modulus, \(E\))와 점성 감쇄율을 비침습적으로 정량 계측하는 의료 진단 영상 기술입니다. 간경변, 유방암 악성 결절 등 병변 부위의 조직 섬유화는 전단파 전파 속도를 급격히 가속시킵니다. 본 고에서는 10,000Hz 0-Alloc C# 엔진으로 2차원 교차 상관 변위 신호로부터 전단파 도달 시간(Time-of-Flight)을 역해석하고 음향 푸시 빔 초점을 능동 보정하는 서보 아키텍처를 구현합니다.

2. 수학적 모델링 및 지배 방정식

등방성 점탄성 매질에서 전단파 전파 속도 \(c_s\)와 영률(Young's Modulus, \(E\)), 전단 탄성률(\(\mu\)) 간의 관계식은 다음과 같습니다:

c_s = \sqrt{\frac{\mu}{\rho}} \approx \sqrt{\frac{E}{3\rho}} \implies E = 3 \rho c_s^2

생체 조직 밀도 \(\rho \approx 1000\,\text{kg/m}^3\)에서 정상 간(\(c_s \sim 1.0\text{--}1.2\,\text{m/s}\))은 약 \(3\text{--}4\,\text{kPa}\)이나, 심각한 간 섬유화 4단계(\(c_s > 2.5\,\text{m/s}\))에서는 \(E > 20\,\text{kPa}\)로 급증합니다. 포이트(Voigt) 점탄성 감쇄 모델에 따른 주파수 분산 관계는 다음과 같습니다:

c_s(\omega) = \sqrt{\frac{2(\mu^2 + \omega^2 \eta^2)}{\rho(\mu + \sqrt{\mu^2 + \omega^2 \eta^2})}}

환자의 호흡 및 맥박에 의한 모션 아티팩트를 제거하기 위해 조직 변위 상관 계수 \(r_{cross} > 0.95\) 조건에서만 10,000Hz 피드백으로 전단파 속도를 샘플링합니다.

3. 실시간 피드백 제어 아키텍처

graph TD A[128채널 선형 초음파 프로브 및 인체 간/유방 조직 표면 접촉] --> B[집속 음향 복사력 ARFI 푸시 빔 인가로 수 마이크로미터 수직 변위 발생] B --> C[측면 수평 방향으로 1~5 m/s 속도의 초음파 전단파 전파 개시] C --> D[5,000fps 초고속 평면파 연속 송수신으로 전단 변위 파면 추적] D --> E[IQ 빔포밍 데이터로부터 1D/2D 상호 상관 기반 전단파 도달 시간 ToF 계측] E --> F[10,000Hz 0-Alloc C# 전단파 속도 cs 및 영률 E = 3*rho*cs^2 실시간 산출] F --> G[조직 점탄성 경계선 추출 및 ARFI 푸시 빔 초점 깊이 최적화 지령 송출] G --> H[조직 강도 탄성 맵 컬러 오버레이 실시간 렌더링 및 간 섬유화 진단] H --> B

4. 고성능 0-Alloc C# 서보 엔진 구현

가비지 컬렉터 부하 없이 10,000Hz 주기로 전단파 속도를 계산하고 탄성 계수를 산출하는 의료 진단급 C# 엔진입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.ShearWaveElastography
{
    public readonly record struct SweConfig(
        double TissueDensityKgPerM3,    // 조직 밀도 ρ (예: 1000.0 kg/m^3)
        double LateralSpacingMm,        // 래터럴 추적점 간격 (예: 2.5 mm)
        double MinValidCorrelation,     // 최소 유효 상호 상관 계수 (예: 0.92)
        double FibrosisThresholdKpa     // 간 섬유화 판정 임계치 (예: 12.5 kPa)
    );

    public ref struct SweServoState
    {
        public double MeasuredShearSpeedMps;
        public double YoungsModulusKpa;
        public double MotionArtifactIndex;
        public double RecommendedPushDepthMm;
        public bool IsTissueHardened;
    }

    public static class ShearWaveElastographyEngine
    {
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void UpdateSweServo(
            in SweConfig cfg,
            ref SweServoState state,
            double measuredTimeOfFlightMs,
            double correlationScore,
            double currentPushDepthMm)
        {
            state.RecommendedPushDepthMm = currentPushDepthMm;

            // 1. 상관 계수 기반 모션 아티팩트 판정
            state.MotionArtifactIndex = 1.0 - Math.Clamp(correlationScore, 0.0, 1.0);

            // 2. 도달 시간(Time-of-Flight)으로부터 전단파 속도 산출
            // cs = Δx / Δt
            double dtSec = Math.Max(1e-6, measuredTimeOfFlightMs * 1e-3);
            double dxM = cfg.LateralSpacingMm * 1e-3;
            state.MeasuredShearSpeedMps = dxM / dtSec;

            // 3. 영률(Young's Modulus) 계산: E = 3 * ρ * cs^2 (kPa 단위)
            double cs = state.MeasuredShearSpeedMps;
            double modulusPa = 3.0 * cfg.TissueDensityKgPerM3 * (cs * cs);
            state.YoungsModulusKpa = modulusPa * 1e-3;

            // 4. 병변 경화(섬유화/악성 종양) 판정
            state.IsTissueHardened = (correlationScore >= cfg.MinValidCorrelation) && 
                                    (state.YoungsModulusKpa >= cfg.FibrosisThresholdKpa);
        }
    }
}

5. 교육 및 산업 응용 가치

본 서보 시스템은 대학병원 프리미엄 임상 초음파 장비의 핵심 소프트웨어 백본인 정량 전단파 탄성 영상의 전 파이프라인을 온디바이스 C#으로 재현합니다. 의공학 및 물리 교육과정에서는 고전 파동역학의 횡파/종파 분리 전파 원리와 생체 조직의 점탄성 유변학(Rheology)이 현대 첨단 의료 진단에 어떻게 통합되는지 명확한 디지털 시뮬레이션을 제공합니다.