All-optical photoacoustic tomography via beam deflection

논문 개요

이 논문은 기존 광음향 영상 (Photoacoustic Imaging, PAI) 의 한계를 극복하기 위해 제안된 새로운 전광학 (All-optical) 광음향 단층촬영 (Tomography) 모달리티를 소개합니다. 기존 방식은 압력 변환기 (Transducer) 를 사용하여 음파를 감지하지만, 이 연구에서는 광학 빔의 편향 (Beam Deflection) 현상을 이용하여 음압 구배 (Pressure Gradient) 를 측정하고, 이를 통해 초기 압력 분포를 재구성하는 알고리즘을 개발했습니다.

1. 문제 제기 (Problem)

기존 광음향 영상 기술은 높은 광학적 대비와 깊은 조직 침투력을 장점으로 하지만, 다음과 같은 근본적인 한계가 존재합니다.

• 압력 변환기의 한계: 기존 시스템은 압전 소자 (PZT, PVDF 등) 어레이를 사용하여 음파를 감지합니다.
  • 저감도 및 왜곡: 박막 변환기는 감도가 낮고, 결정체나 고체 변환기는 기계적 공진으로 인해 실제 광음향 파형이 왜곡되는 문제가 있습니다.
  • 대역폭 제한: 변환기의 주파수 응답 특성이 영상의 해상도를 제한합니다.
• 측정 정확도: 기존 센서들의 물리적 한계로 인해 측정 오차가 발생하며, 이는 영상의 정밀도를 떨어뜨립니다.

2. 방법론 (Methodology)

가. 물리적 모델링 (Physical Modeling)

• 원리: 광음향 효과로 생성된 음파는 매질의 밀도 변동을 일으키며, 이는 굴절률 변화를 유발합니다. 탐사 빔 (Probe beam) 이 이 음장을 통과할 때 빔이 편향됩니다.
• 측정 신호: 편향된 빔의 각도는 광음향 파 내의 **밀도 구배 (즉, 압력 구배, $\nabla p$)**에 비례합니다.
• 수학적 표현: 광검출기 (Photodetector) 가 측정하는 신호 $\beta$는 압력 구배의 **라돈 변환 (Radon Transform)**으로 모델링됩니다.
  $$\beta = \int_{z_L}^{z_U} \mathbf{n}' \cdot \nabla p(\mathbf{x}, \tau) \, dz$$
  여기서 $\mathbf{n}'$는 측정 평면의 법선 벡터이며, 적분은 빔 경로 ($z_L$에서 $z_U$) 를 따라 수행됩니다.

나. 역문제 해결 알고리즘 (Inverse Problem Algorithm)

측정된 데이터로부터 초기 압력 분포 ($p_0$) 를 복원하는 역문제는 불완전하고 불안정 (Ill-posed) 하므로, 최적화 기반 접근법을 사용합니다.

1. 방향성 전처리 (Directional Preprocessing):
   • 측정된 라돈 변환 데이터를 역라돈 변환 (Inverse Radon Transform) 하여 각 방향 ($x_1, x_2, x_3$) 에 대한 **압력 구배 필드 ($\nabla p$)**를 먼저 복원합니다.
   • 이 과정은 고주파 성분을 증폭시키고 저주파 배경 노이즈를 억제하여 구조적 정보를 강조합니다.
2. 최적화 기반 복원 (Optimization-based Reconstruction):
   • **ISTA (Iterative Shrinkage-Thresholding Algorithm)**를 사용하여 압력 구배 필드를 반복적으로 업데이트합니다.
   • 목적 함수: 측정 데이터와 시뮬레이션 데이터 간의 오차 (Sensor Loss) 를 최소화하면서, $L_1$ 정규화 항을 추가하여 에지 (Edge) 보존과 희소성 (Sparsity) 을 유도합니다.
   • 어드join트 연산자 (Adjoint Operator): 오차 항을 역전파하여 초기 추정을 수정합니다.
3. 갈레르킨 방법 (Galerkin Method) 을 통한 압력 복원:
   • 복원된 3 차원 압력 구배 필드 ($\nabla p$) 를 입력으로 받아 **푸아송 방정식 (Poisson Equation)**을 수치적으로 적분 (유한 요소법, FEM) 하여 최종 **초기 압력 필드 ($p_0$)**를 재구성합니다.

다. 실험 설정 (Numerical Simulation)

• 시뮬레이션 도구: k-Wave 패키지를 사용하여 파동 전파를 모사하고, Python-Scikit-Image 를 사용하여 라돈 변환을 수행했습니다.
• 센서 배치: 입방체 영역의 3 개 이상 서로 직교하는 면에 센서를 배치하여 모든 방향의 구배 정보를 수집하도록 설계되었습니다.

3. 주요 기여 (Key Contributions)

1. 새로운 영상 모달리티 제안: 압력 변환기 대신 광학 빔 편향을 이용한 전광학 광음향 단층촬영 방식을 제안했습니다.
2. 다중 물리 현상 모델링: 광학 빔 편향과 음파 전파를 결합한 정교한 물리 모델을 수립했습니다.
3. 강건한 3D 역문제 해결 알고리즘: 최적화 (ISTA) 와 갈레르킨 방법을 결합하여 3 차원 공간에서 압력 구배와 압력 필드를 안정적으로 복원하는 알고리즘을 개발했습니다.
4. 성능 평가: 단순한 구형 소스부터 복잡한 Shepp-Logan 팬텀, 기계적 메타물질까지 다양한 시나리오에서 알고리즘의 유효성을 입증했습니다.

4. 결과 (Results)

• 수렴성 (Convergence): Shepp-Logan 팬텀 (65³ 격자) 시뮬레이션에서 50 회 반복 후 압력 구배 필드의 RMSE 는 약 60%, 압력 필드는 약 32% 로 수렴했습니다. PSNR(피크 신호 대 잡음비) 은 초기 반복부터 빠르게 향상되었습니다.
• 노이즈 내성 (Noise Robustness): 센서 측정 오차 (가우시안 노이즈) 가 최대 1.0 에 달하는 극한 조건에서도 제안된 방식은 기존 압력 기반 방식보다 **SSIM(구조적 유사성)**과 PSNR에서 우수한 성능을 보였습니다. 이는 갈레르킨 기반의 푸아송 방정식 해법이 무작위 노이즈를 효과적으로 억제하기 때문입니다.
• 세부 구조 복원:
  • 큐브 팬텀: 날카로운 모서리와 작은 돌출부 (bulge) 를 정확하게 복원하여 미세 구조 검출 능력을 입증했습니다.
  • 스피노달 메타물질 (Spinodal Metamaterial): 나노/마이크로 스케일의 복잡한 층상 구조를 광음향 신호로 비파괴적으로 재구성하는 데 성공했습니다.

5. 의의 및 결론 (Significance)

• 해상도 및 감도 향상: 광학 센서의 넓은 대역폭과 높은 감도 덕분에 기존 압전 센서의 기계적 왜곡 없이 더 선명한 영상을 얻을 수 있는 잠재력을 가집니다.
• 비파괴 검사 가능성: 미세 기계적 메타물질과 같은 복잡한 구조물의 동적 특성을 비파괴적으로 평가할 수 있는 새로운 도구를 제공합니다.
• 임상 및 연구 적용: LED 기반 광원 및 소형 광학 시스템과 결합하여 휴대용 및 저비용 광음향 영상 시스템 개발의 길을 열 수 있습니다.
• 한계 및 향후 과제: 현재는 수치 시뮬레이션 단계이며, 실제 실험을 통한 검증이 필요합니다. 또한, 굴절률 구배가 큰 이질 매질에서는 광선 추적 (Ray tracing) 등 더 정교한 광학 모델링이 필요하며, 계산 비용이 높다는 점은 개선이 필요한 부분입니다.

요약하자면, 이 연구는 광학 빔 편향 원리를 활용하여 광음향 신호를 측정하고, 최적화 알고리즘과 수학적 적분을 결합하여 고해상도 3D 영상을 재구성하는 혁신적인 전광학 단층촬영 시스템을 제안했습니다. 이는 기존 광음향 영상의 물리적 한계를 극복하고, 미세 구조 이미징 및 비파괴 검사 분야에서 중요한 진전을 이루는 계기가 될 것으로 기대됩니다.