[논문 리뷰] Full-Wave Iterative Image Reconstruction in Photoacoustic Tomography with Acoustically Inhomogeneous Media

Chao Huang, Kun Wang외 3명|arXiv (Cornell University)|2013. 03. 22.

Photoacoustic and Ultrasonic Imaging참고 문헌 39인용 수 7

한 줄 요약

이 논문은 정확한 광음파 파동 방정식에 기반한 이산 이미징 모델을 사용하여 음향적으로 비균일한 매질에서 광음파 단층촬영(PACT)을 위한 전파형 반복 이미지 복원 방법을 제안한다. k-space 가짜스펙트럼 방법을 통해 일치하는 정방향 및 역방향 프로젝션 연산자를 구현함으로써, 노이즈, 불완전한 데이터, 음향 비균일성으로 인한 잡음 감소와 함께 강건한 복원을 가능하게 하였으며, 시뮬레이션 및 실험 데이터를 통해 검증되었다.

ABSTRACT

Existing approaches to image reconstruction in photoacoustic computed tomography (PACT) with acoustically heterogeneous media are limited to weakly varying media, are computationally burdensome, and/or cannot effectively mitigate the effects of measurement data incompleteness and noise. In this work, we develop and investigate a discrete imaging model for PACT that is based on the exact photoacoustic (PA) wave equation and facilitates the circumvention of these limitations. A key contribution of the work is the establishment of a procedure to implement a matched forward and backprojection operator pair associated with the discrete imaging model, which permits application of a wide-range of modern image reconstruction algorithms that can mitigate the effects of data incompleteness and noise. The forward and backprojection operators are based on the k-space pseudospectral method for computing numerical solutions to the PA wave equation in the time domain. The developed reconstruction methodology is investigated by use of both computer-simulated and experimental PACT measurement data.

연구 동기 및 목표

• 음향적으로 비균일한 매질에서 기존 PACT 복원 방법의 한계를 해결하기 위해, 약하게 변화하는 매질에 대한 가정과 높은 계산 비용을 줄이기 위한 목표.
• 실제 PACT 응용에서 데이터의 불완전성, 노이즈, 모델링 오차로 인한 잡음 문제를 해결하기 위한 목표.
• 변동하는 음속, 밀도, 감쇠를 고려한 비균일 매질에서의 파동 전파를 정확히 표현하기 위해 정확한 광음파 파동 방정식에 기반한 이산 이미징 모델을 개발하기 위한 목표.
• 복잡한 매질에서 반복 복원 알고리즘을 효율적으로 구현하기 위해 정방향 및 역방향 프로젝션 연산자 쌍을 확립하기 위한 목표.
• 최신 반복 복원 기법(예: 원시-쌍대 하이브리드 경사법 기반 기법)을 현실적인 음향 비균일성을 포함한 전파형 PACT 모델에 적용할 수 있도록 하기 위한 목표.

제안 방법

• 정확한 광음파 파동 방정식에 기반한 이산 이미징 모델을 수립하여, 공간적으로 변화하는 음속, 질량 밀도, 주파수-거듭제곱법 감쇠를 포함하기.
• PA 파동 방정식의 시간 도메인 수치적 해를 정확하게 구하기 위해 k-space 가짜스펙트럼 방법을 사용하여 정방향 및 역방향 프로젝션 연산자를 구현하기.
• 반복 복원 알고리즘의 효율적 적용을 가능하게 하기 위해 정방향, 역방향 및 그 수반 연산자를 포함한 일치하는 연산자 쌍을 개발하기.
• 역문제를 노이즈 및 불완전한 데이터 조건에서 안정적으로 해결하기 위해 연산자를 원시-쌍대 하이브리드 경사법(PDHG) 알고리즘에 통합하기.
• 핵심 연산자 정의: $\mathcal{P}_l$ 은 이산 발산 유사 연산을 계산하기 위한 것이며, $\mathcal{P}_c$ 는 비음수성 강제를 위한 것이며, $\mathcal{P}_l^T$ 는 수반 연산을 위한 것이며, $\mathcal{P}_p$ 는 경사 업데이트에서 벡터 정규화를 위한 것이다.
• 알고리즘을 컴퓨터 시뮬레이션 및 실험적 PACT 데이터에 적용하여 실제 비균일 매질에서의 성능을 검증하기.

실험 결과

연구 질문

• RQ1 정확한 광음파 파동 방정식에 기반한 전파형 반복 복원 프레임워크는 음향 비균일성이 존재할 경우 PACT에서 발생하는 잡음 요소를 효과적으로 줄일 수 있는가?
• RQ2 제안된 정방향 및 역방향 프로젝션 연산자 쌍은 데이터의 불완전성과 노이즈 존재 조건에서 복원 정확도와 수렴성을 어떻게 향상시키는가?
• RQ3 k-space 가짜스펙트럼 방법은 3차원 비균일 매질에서 PACT를 위한 효율적이고 정확한 파동 전파 모델링을 얼마나 잘 가능하게 하는가?
• RQ4 강한 음향 비편향(예: 두개골에 의한 군집파면 왜곡)이 존재하는 실험 데이터에서 반복 복원 방법은 어떻게 성능을 발휘하는가?
• RQ5 실제 측정 조건에서 전통적인 시간 역전 또는 레이 기반 접근법에 비해 제안된 방법이 잡음 억제 및 공간 해상도 측면에서 우수한 성능을 보일 수 있는가?

주요 결과

• 제안된 전파형 반복 복원 방법은 두개골에 의한 군집파면 왜곡 등 음향 비균일성으로 인한 영상 잡음 요소를 크게 감소시키며, 이는 뇌를 통과하는 PACT 실험 데이터를 통해 입증되었다.
• 특히 데이터가 불완전하거나 노이즈가 많은 조건에서 기존의 시간 역전 및 레이 기반 접근법에 비해 더 높은 영상 품질과 공간 해상도를 달성하였다.
• k-space 가짜스펙트럼 방법의 사용은 강한 음향 비균일성 조건에서도 PA 파동 방정식의 정확하고 안정적인 수치적 해를 가능하게 하였다.
• 일치하는 정방향 및 역방향 프로젝션 연산자 쌍은 반복 알고리즘의 효율적 구현을 가능하게 하여 계산 부담을 줄이면서도 높은 복원 정밀도를 유지하였다.
• PDHG 기반 알고리즘은 희소하거나 노이즈가 많은 측정 데이터 조건에서도 안정적으로 수렴하며 정확한 복원을 생성하였으며, 이는 시뮬레이션 및 실험 데이터셋을 통해 검증되었다.
• 정확한 파동 방정식에 기반한 방법적 기반 덕분에 음속 및 감쇠의 변화와 같은 모델링 오차와 물리적 매개변수의 불확실성에 대해서도 강건성을 보였다.