다중에코 방사형 FLASH와 모델 기반 재구성을 이용한 고해상도 무왜곡 태아 뇌 R2 매핑

초록

본 연구는 3T MRI에서 모성 자유 호흡 중에도 2 초/슬라이스로 1.1 × 1.1 × 3 mm³ 해상도의 무왜곡 R2* 지도와 동시에 물‑지방 분리 및 B0 맵을 제공하는 2D 다중에코 방사형 FLASH 시퀀스와 보정 없는 모델 기반 재구성 방법을 제안한다. 수치·NIST 팬텀 및 26–36주 태아 9명을 대상으로 검증했으며, 기존 PICS‑Graph Cut 및 다중에코 EPI와 비교해 높은 정확도(ICC > 0.96), 재현성(ICC ≈ 0.92)와 왜곡 감소를 확인하였다.

상세 분석

이 논문은 태아 뇌의 정량적 R2* 매핑을 위해 두 가지 핵심 기술을 결합하였다. 첫 번째는 블립 그라디언트를 삽입해 에코 차원에서 회전하는 2D 다중에코 방사형 FLASH 시퀀스로, 이는 전통적인 EPI에 비해 읽기 시간(읽기 라인)이 짧아 B0 비균일성에 의한 기하학적 왜곡을 크게 감소시킨다. 또한 골든 앵글(≈68.75°) 기반의 스포크 회전 설계와 35개의 에코를 활용해 충분한 T2* 감쇠 정보를 확보함으로써, 태아 뇌와 같이 T2가 짧은 조직에서도 안정적인 피팅이 가능하도록 하였다. 두 번째는 물·지방, R2, B0를 동시에 추정하는 비보정 모델 기반 재구성이다. 신호 모델식 M(TE)= (W+F·z)·exp(i2πfB0·TE)·exp(−R2*·TE) 를 기반으로, 코일 감도와 샘플링 마스크를 포함한 비선형 연산자를 정의하고, ℓ1‑Wavelet 스파시티와 Sobolev 정규화를 결합한 비용 함수를 최소화한다. IRGNM‑FISTA 알고리즘을 GPU 가속으로 구현해 5–10분 내에 수렴하도록 설계하였다. 정규화 파라미터는 가우스‑뉴턴 반복마다 1/3씩 감소시켜 초기 강한 억제에서 점진적 세부 복원으로 전환한다. 이러한 접근은 전통적인 2단계(이미지 재구성 → 픽셀별 피팅) 방식보다 파라미터 공간을 직접 제약함으로써 undersampling에 대한 강인성을 높이고, 물·지방 신호가 겹치는 주변 조직에서도 정확한 분리를 가능하게 한다. 수치 시뮬레이션에서는 10 %~95 % 지방 함량, −50 ~ +50 Hz 오프레조넌스, R2* 10–200 s⁻¹ 범위 전반에 걸쳐 평균 오차가 0.03 % 이하(FF)와 0.17 s⁻¹ 이하(R2*)에 머물렀으며, 노이즈 레벨이 증가해도 Bland‑Altman 분석에서 편향이 거의 없었다. 실험 팬텀(NIST)과 인체 태아 데이터에서도 PICS‑Graph Cut과의 ICC가 0.96, 반복성 ICC가 0.92에 달해 재현성이 뛰어남을 확인했다. 특히 EPI와 비교했을 때 해상도는 1.1 mm에서 2–3 mm로 향상되고, 왜곡이 현저히 감소했음에도 R2* 값 자체는 ICC 0.80으로 통계적으로 일치하였다. 이는 모델 기반 재구성이 고해상도와 정량적 정확도 사이의 전통적 트레이드오프를 효과적으로 해소함을 의미한다.