기술 요약: VITO – 역 위상 최적화를 이용한 혈관 기하학 및 혈류 추정

문제 정의

심혈관 질려(CVD) 평가는 환자 맞춤형 혈관 기하학으로부터 도출되는 벽 전단 응력(WSS) 및 분획 혈류 예비력(FFR)과 같은 혈역학적 기술자에 크게 의존한다. 표준 임상 워크플로우는 일반적으로 대비 증강 컴퓨터 단층 촬영(CE-CT)을 활용하여 이러한 기하학적 구조를 재구성한다. 그러나 일반적인 파이프라인(전형적으로 필터 역투영(FBP) 또는 반복 재구성 후 세그멘테이션을 수행하는 방식)은 방사선량 제한 및 생리적 움직임 아티팩트로 인한 희소하고 노이즈가 많은 시계열 획득 데이터에 취약하다. 이러한 한계는 언더샘플링 아티팩트를 유발하여, 누락된 분지나 인위적인 협착과 같은 기하학적 부정확성을 초래한다.

또한, 기존의 물리 법칙 제약 기반 역문제 방법들은 대개 고정된 사전 세그멘테이션 기하학을 가정한다. 잘못된 고정 기하학에 유동 물리학을 강제하는 것은 유효하지만 실제와는 다른 유동장을 생성하며, 이는 결정적인 혈역학적 파라미터를 포착하지 못하게 만든다. 최근 일부 연구에서 기하학적 구조와 유동을 공동으로 재구성하려는 시도가 있었으나, 이들은 원시 CT 투영 데이터(시노그램)로부터 직접 기하학을 복구하기보다는 측정된 속도 이미지(예: MRI)에 의존하는 경우가 많다. 따라서 사전 기하학적 가정 없이 원시 CE-CT 시노그램으로부터 혈관 위상과 유동 파라미터를 공동으로 복구할 수 있는 프레임워크에는 여전히 공백이 존재한다.

방법론

저자들은 혈관 기하학 복원을 밀도 기반의 위상 최적화(Topology Optimization, TO) 문제로 정식화한 물리 제약 재구성 프레임워크인 VITO를 제안한다. 기하학을 고정된 입력값으로 취급하는 대신, 혈관 루멘(lumen)을 미지의 설계 변수로 취급한다.

핵심 프레임워크

이 방법은 지배적인 물리 법칙을 따를 때 시뮬레이션된 투영 데이터가 획득된 희소 시노그램과 일치하도록 하는 함수를 찾는 역문제를 해결한다.

1. 설계 표현:

   • 계산 영역은 각 요소에 의사 밀도 변수 $\gamma_e \in [0, 1]$가 할당된 메쉬 요소로 이산화된다.
   • $\gamma = 0$은 유체(혈관 루멘)를 나타내고, $\gamma = 1$은 고체(조직/벽)를 나타낸다.
   • 재료 특성은 RAMP(Rational Approximation of Material Properties) 스킴을 사용하여 보간된다.
     • 브링크만 페널티 ($\alpha$): 고체 영역에서의 큰 페널티는 유체 속도를 억제하여, 명시적인 경계 추적 없이도 노슬립(no-slip) 조건을 효과적으로 강제한다.
     • 확산 계수 ($D$): 고체 영역 내에서 대비제 수송을 제한하여 인위적인 확산을 방지한다.

2. 지배 물리 (순방향 모델):

   • 정상 상태 유동: 브링크만 페널티 항을 사용하여 비압축성 나비에-스토크스 방정식을 풀고 유체 속도장 $\mathbf{u}$를 모델링한다.
   • 과도 상태 수송: 정체된 속도장 $\mathbf{u}$와 시변 입구 경계 조건을 기반으로 대비제 농도 $c$의 수송을 모델링하는 과도 이류-확산 방정식을 사용한다.
   • 순방향 투영: 결과물인 시계열 농도장은 비어-람베르트 법칙에 기반한 미분 가능한 평행 빔 투영 연산자를 사용하여 측정 공간으로 매핑된다. 이를 통해 예측 시노그램 $\hat{g}$가 생성된다.

3. 최적화 루프:

   • 목적 함수: 예측된 시노그램과 획득된 희소/노이즈 시노그램 데이터($g_{data}$) 사이의 평균 제곱 오차(MSE)이다.
   • 설계 변수: 의사 밀도 필드 $\gamma$(기하학)와 입구 속도 스케일링 파라미터 $a$(유동 조건)이다.
   • 민감도 분석: 프레임워크는 엔드 투 엔드 미분 가능한 파이프라인을 구현하기 위해 JAX 생태계를 활용한다.
     • 내재 함수 정리 (Implicit Function Theorem, IFT): 반복적인 뉴턴-라프슨 단계의 메모리 비용을 피하면서 정상 상태 나비에-스토크스 솔버의 민감도를 계산하는 데 사용된다.
     • 체크포인팅 (Checkpointing): 시간 단계에 따른 역방향 모드 자동 미분(AD) 중 메모리 사용량을 관리하기 위해 과도 이류-확산 솔버에 적용된다.
   • 최적화 도구: 설계 변수를 업데이트하기 위해 이동 점근선 방법(Method of Moving Asymptotes, MMA)이 사용된다. 이진 유체-고체 위상으로 해를 유도하고 지역 최솟값을 피하기 위해 연속 전략(페널티 및 투영 선명도를 점진적으로 증가시킴)이 채택된다.

주요 기여

• 결합 역재구성: 본 논문은 사전 세그멘테이션된 기하학 없이도 원시 CT 시노그램으로부터 혈관 위상과 입구 유동 파라미터를 동시에 복구하는 통합 프레임워크를 도입한다.
• 물리 기반 정규화: 혈역학적 물리학을 역문제와 결합함으로써, 단순한 이미지 사전 확률(prior)만이 아닌 물리 법칙을 통해 재구성을 정규화한다.
• 미분 가능한 다중 물리 파이프라인: 구현 과정에서 JAX를 활용하여 복잡한 결합 유동-수송-투영 시스템의 민감도를 처리하며, IFT와 체크포인팅을 사용하여 계산 가능성을 확보하였다.
• 데이터 제한에 대한 강건성: 저선량 또는 고속 획득 임상 CT에서 전형적인 희소 뷰(sparse-view) 및 노이즈가 많은 획득 조건에 대응하도록 설계되었다.

실험 결과

본 방법은 세 가지 혈관 시나리오를 나타내는 합성 팬텀을 통해 검증되었다:

1. 분지 동맥 (Bifurcated Artery): 희소 시노그램으로부터 분지 위상과 입구 유동 스케일링을 성공적으로 복구하였다.
2. 협착 동맥 (Stenosed Artery): 국소적인 50% 루멘 협착(stenosis)과 유동 파라미터를 정확하게 재구성하였다.
3. 협착 경동맥 분지 (Stenosed Carotid Bifurcation): 분지와 비대칭 협착이 결합된 복잡한 기하학을 복구하였다.

성능 지표:

• 정확도: 재구성된 기하학은 모든 벤치마크에서 참값(ground truth) 대비 높은 Dice 유사 계수(DSC > 0.99)와 낮은 정규화된 제곱 평균 제곱 오차(nRMSE)를 달로했다.
• 노이즈 및 희소성 강건성: 초희소 데이터(10개의 투영)와 높은 노이즈 수준(25% 상대 가우시안 노이즈)이 포함된 실험에서, 제안된 방법은 위상 일관성(DSC = 0.9830)을 유지하였다. 반면, 임계값 처리를 적용한 베이스라인 FBP 방식은 유사한 조건에서 경계 왜곡, 아티팩트 및 분지 연속성 상실를 겪었다.
• 유동 추정: 프레임워크는 입구 유동 스케일링 파라미터를 참값에 근접하게(1에 가깝게) 복구하였으며, 이는 기하학과 함께 경계 조건을 추론하는 능력을 입증한다.

의의 및 주장

본 논문은 재구성 과정에 기저의 혈역학을 통합하는 것이 강력한 기하학적 정규화를 제공하여, 희소 임상 투영 데이터로부터 복잡한 혈관 특징(분지, 협착)을 직접 복구할 수 있게 한다고 주장한다. 이 접근법은 현재의 워크플로우에서 발생하는 기하학적 오류가 혈역학 분석으로 전파되는 문제를 해결한다.

저자들은 이 연구가 희소 임상 데이터로부터 신뢰할 수 있는 혈관 분석을 위한 중요한 단계라고 위치시킨다. 저자들은 현재 연구가 2D 정상 유동 및 강체 경계로 제한되어 있음을 명시적으로 밝힌다. 향-후 연구로는 3D 과도 유동 확장, 탄성 벽을 포함한 유체-구조 상호작용(FSI) 도입, 입구/출구 위치를 미지수로 취급하는 방식, 불확실성 정량화, 그리고 특정 스캐너 기하학(예: 팬 빔 CT)에 대한 투영 연산자 적응 등이 필요하다고 언급하였다. 본 논문은 즉각적인 임상 배포를 주장하는 것이 아니라, 그러한 응용을 위한 수학적 및 계산적 토대를 구축하는 것을 목표로 한다.