An angular distortion matrix approach for joint wave-speed tomography and aberration correction in scattering media

기술 요약: 공동 파동 속도 토모그래피 및 수차 교정을 위한 각도 왜곡 행렬 접근법

문제 정의
반사 이미징(초음파, 광학, 지진파)에서, 정확한 매질 반사율 재구성은 파동 속도 분포에 의해 제어되는 정밀한 파동 전파 모델에 의존한다. 그러나 복잡한 산란 매질에서 파동 속도는 사전에 알려진 경우가 드물다. 가정된 파동 속도와 실제 파동 속도 사이의 불일치는 위상 수차(phase aberrations)를 유발하여 초점을 저하시키고, 점 퍼짐 함수(point-spread function)를 넓히며, 산란체의 겉보기 위치를 왜곡시킨다. 이러한 수차는 근저에 있는 파동 속도 불균질성에 대한 정보를 인코딩하지만, 이를 매질의 반사율과 분리해내는 것은 특히 스펙클(speckle)이 지배적인 매질에서 매우 어렵다. 기존 방법들은 근본적인 결합 문제에 직면해 있다: 수차를 교정하려면 정확한 파동 속도 모델이 필요하지만, 이 모델을 추출하려면 수차에 의해 편향되지 않은 관측값이 필요하다. 또한, 전통적인 왜곡 행렬 접근법은 수차가 입사 경로 또는 반사 경로 중 하나에서 발생한다고 가정하지만, 실제로는 두 기여분이 본질적으로 혼합되어 있어, 속도 모델을 업데이트하지 않고서는 해결할 수 없는 깊이 이동(depth shifts)을 초래한다.

방법론
저자들은 파동 속도 분포를 반복적으로 추정하고 수차를 교정하기 위해 새로운 각도 왜곡 행렬($D_{\theta\theta}$)을 사용하는 통합 프레임워크를 제안한다. 방법론은 다음과 같이 진행된다:

1. 행렬 이미징 및 주파수-파수 표현:
   과정은 반사 행렬 $R_{uu}$(트랜스듀서-대-트랜스듀서 응답)를 획득하는 것으로 시작된다. 이는 초점 기저(focused basis)로 투영되어 초점 반사 행렬 $R_{xx}$를 형성한다. 왜곡을 격리하기 위해, 저자들은 주파수-파수 표현의 파동 왜곡 $D_{kk}$를 도입한다. 이 연산자는 임의의 지점 $x$로부터 산란이 발생한 것처럼 가정하여 기대 기하학적 위상 성분을 제거함으로써 구축된다. 이를 통해 사후 처리 과정에서 임의의 가상 실험에 대한 왜곡을 합성할 수 있다.

2. 각도 왜곡 행렬 구축:
   $D_{kk}$ 연산자는 **가우시안 빔 기저(Gaussian-beam basis)**로 투영되어 각도 왜곡 행렬 $D_{\theta\theta}$를 형성한다. 이 기저는 입사각($\theta_{in}$)과 반사각($\theta_{out}$)으로 정의된 조향된 가우시안 빔을 활용한다. 이 투영은 파동이 무질서함에도 불구하고 특정 각도 방향을 따라 강한 상관관계를 보이는 **각도 메모리 효과(angular memory effect)**를 활용한다. $D_{\theta\theta}$의 안티대각선(constant mid-angle $\theta_m = (\theta_{in} + \theta_{out})/2$)을 분석함으로써, 이 방법은 잘 정의된 전파 경로와 관련된 위상 수차를 격리한다.

3. 국소 수차 추정:
   시야(FOV)는 겹치는 "등방성 패치(isoplanatic patches)"로 나뉜다. 각 패치 내에서는 각도 왜곡 행렬 요소로부터 국소 상관 행렬이 계산된다. 파동 왜곡이 작은 영역 내에서 불변한다는 가정하에, 이 상관 행렬의 첫 번째 고유벡터는 국소 위상 수차 법칙($u_m$)의 추정치를 제공한다. 위상 에일리어싱(aliasing)을 방지하기 위해, 이 방법은 연속적인 반-차이 각도($\theta_d$) 사이의 상대적 위상 차이를 활용한다.

4. 파동 속도 토모그래피를 위한 반복 역전(Iterative Inversion):
   추정된 위상 수차는 음속 분포를 재구성하는 데 사용된다. 문제는 측정된 위상 차와 전파 모델에 의해 예측된 위상 차 사이의 미스핏(misfit)을 최소화하는 최적화 과제로 공식화된다. 역전은 다음과 같이 반복적으로 수행된다:

   • 초기 균질 모델(예: $c_0 = 1540$ m/s)을 가정한다.
   • 수차를 추정하고 교정한다.
   • 직선 경로 근사(Jacobian을 위한)와 티코노프 규제화(Tikhonov regularization)를 사용하는 선형 역문제(linearized inverse problem)를 통해 슬로니스 섭동($\Delta s$)을 구한다.
   • 파동 전파 모델을 업데이트하고, 잔류 수차가 무시할 수 있을 정도로 작아질 때까지 이 과정을 반복한다.

5. 다중 산란 처리:
   메모리 효과를 저하시키는 고도로 불균질한 매질(예: 유방 조직)의 경우, 왜곡 행렬을 계산하기 전에 반사 행렬에 **적응형 컨포컬 필터(adaptive confocal filter)**를 적용한다. 이는 단일 산란과 관련된 대각 에너지(diagonal energy)는 보존하면서 비간섭적 배경 노이즈를 억제한다.

주요 결과
이 접근법은 세 가지 시나리오에서 초음파 실험을 통해 검증되었다:

• 조직 모사 팬텀: 알려진 층상 구조(지방, 근육, 간)와 다양한 음속(1420–1550 m/s)을 가진 팬텀이 사용되었다. 반복법은 평균 오차가 참값의 1% 이내인 음속 지도를 성공적으로 복원하였다. 스트렐 비율(Strehl ratio, 초점 품질의 척도)은 0.34에서 0.96으로 향상되었으며, 반사율 이미지는 날카로운 계면과 향상된 스펙클 대비를 보여주었다.
• 생체 내 유방 이미징: 침습성 암종 환자에게 적용되었다. 이 방법은 주변 지방 및 섬유선 조직과 구별되는, 높아진 음속(~1574 m/s)이 나타나는 국소 영역을 식별하였다. 수차 교정은 교정되지 않은 이미지가 해부학적 특징의 겉보기 깊이를 편향시켰음을 드러냈다.
• 생체 내 간 이미징: 건강한 피험자에게 적용되었다. 이 방법은 높은 축 방향 해상도로 복벽(피부, 지방, 근육)과 간 실질의 층상 구조를 분해하였다. 추출된 음속은 건강한 조직의 문헌 값과 일치하였으며, 스트렐 비율은 0.61에서 0.93으로 증가하여 반사 이미지를 크게 날카롭게 만들고 무에코(anechoic) 혈관과 같이 이전에 가려졌던 구조를 드러냈다.

의의 및 주장
본 논문은 각도 왜곡 행렬 접근법이 파동 속도 토모그래피와 수차 교정의 결합된 문제에 대한 강력한 솔루션을 제공한다고 주장한다. 주요 의의는 다음과 같다:

• 통합: 이 방법은 (일반적으로 수차 교정에 사용되는) 왜곡 행렬 개념을 단일 반복 행렬 이미징 프레임워크 내에서 에코 모드 기반의 계산 초음파 토모그래피(CUTE)와 통합한다.
• 모델 독립성: 각도 왜곡 행위의 공식화는 초점을 위한 특정 파동 전파 모델에 독립적이므로, 방법 자체가 모델을 점진적으로 정교화할 수 있게 한다.
• 고해상도: 가우시안 빔 기저로 투영함으로써, 이 방법은 고전적인 왜곡 행렬 접근법보다 더 작은 등방성 패치(약 1 mm)에 대해 수차를 추정할 수 있어 토모그래피에서 더 높은 공간 해상도를 가능하게 한다.
• 임상적 유용성: 이 방법은 양적인 음속 지도(병리학적 지표가 될 수 있는)를 제공하는 동시에, 복잡하고 불균질한 조직에서 구조적 이미지 품질을 개선할 수 있음을 입증하였다.
• 범용성: 초음파에서 시연되었으나, 반사 행렬 이미징을 구현할 수 있는 모든 파동 양식(광학, 지진파)에 적용 가능하다.