[논문 리뷰] Hamilton-Green solver for the forward and adjoint problems in photoacoustic tomography

Francesc Rul⋅lan, Marta M. Betcke|arXiv (Cornell University)|2018. 10. 31.

Photoacoustic and Ultrasonic Imaging참고 문헌 69인용 수 4

한 줄 요약

이 논문은 고주파 근사법을 활용해 광음향단층촬영(PAT)에서 정방향 및 수반 파동 방정식을 푸는 데 사용되는 해밀턴-그린(HG) 솔버를 소개한다. 이 방법은 해밀턴 경로를 따라 그린 함수를 근사함으로써 고주파 영역에서 효율적인, 센서에 종속되지 않는 계산을 가능하게 하며, 전체 파동장 해법 대비 빠른 속도를 제공한다. 이는 음속의 비균일성과 고주파 영역에서의 초점화 및 공명면(카우스틱)이 존재하는 2차원 비균질 영역에서 k-Wave 시뮬레이션과의 비교 검증을 통해 정확한 결과를 도출함을 입증한다.

ABSTRACT

The majority of the solvers for the acoustic problem in Photoacoustic Tomography (PAT) rely on full solution of the wave equation which makes them less suitable for real-time and dynamic applications where only partial data is available. This is in contrast to other tomographic modalities, e.g. X-ray tomography, where partial data implies partial cost for the application of the forward and adjoint operators. In this work we present a novel solver for the forward and adjoint wave equations for the acoustic problem in PAT. We term the proposed solver Hamilton-Green as it approximates the fundamental solution to the respective wave equation along the trajectories of the Hamiltonian system resulting from the high frequency asymptotics for the wave equation. This approach is fast and scalable in the sense that it allows computing the solution for each sensor independently at a fraction of the cost of the full wave solution. The theoretical foundations of our approach are rooted in results available in seismics and ocean acoustics. We present results for 2D numerical phantom with heterogeneous sound speed which we evaluate against a full wave solution obtained with a pseudospectral method implemented in k-Wave toolbox.

연구 동기 및 목표

• 실시간 및 동적 응용에서 부분적 또는 희소한 데이터를 가진 경우에 전체 파동 해법의 계산 비효율성을 해결하기 위해.
• 전체 파동장 계산을 피하고, 정방향 및 수반 파동 방정식을 위한 확장 가능하고 센서에 종속되지 않는 솔버를 개발하기 위해.
• 레이어 추적과 해밀턴 역학을 기반으로 한 고주파 근사 근사를 활용해 PAT의 음향 문제를 정확하게 해결하기 위해.
• 표준 레이 방법이 진폭 폭주로 인해 실패하는 곳인 비균질 매질 내에서의 공명면과 초점 렌즈와 같은 도전적인 특징을 다루기 위해.
• 반복 재구성 및 변분 방법에 활용 가능한 수치적으로 안정적이고 효율적인 전체 파동 해법의 대안을 제공하기 위해.

제안 방법

• 해밀턴-그린 솔버는 고주파 근사 해의 해밀턴 시스템에서 유도된 궤적을 따라 레이를 추적함으로써 파동 방정식의 기본 해를 근사한다.
• 레이의 진폭이 레이 맵의 자코비안 행렬식에 의해 조절되는 그린 함수를 계산하며, 공명면에서의 특이성을 피하기 위해 정규화 조치를 취한다.
• 수반 문제에서는 시간 역전 소스를 사용하며, 센서 데이터가 레이를 따라 역방향 전파되어 初기 압력 분포를 재구성한다.
• 실제 측정 제약 조건과의 호환성을 확보하기 위해 통합 시간 창을 제한하기 위해 부드러운 커팅 함수를 사용한다.
• 공명면에서는 레이 튜브 행렬식 |q(t; x₀)| < ε 인 경우 레이 진폭을 0으로 설정하고, 좌표계 반전에 대비해 q의 절대값을 사용하여 보정한다.
• 2차원에서 구현되었으며, 각 센서당 고정된 수의 레이를 사용하여, 센서 간 독립적인 계산이 가능해 확장성이 확보된다.

실험 결과

연구 질문

• RQ1 레이어 추적 기반의 고주파 근사법이 비균질 매질에서 PAT의 정방향 및 수반 파동 방정식에 대해 정확하고 효율적인 해법을 제공할 수 있는가?
• RQ2 레이 진폭이 발산하고 표준 레이 방법이 실패하는 공명면이 존재하는 상황에서 해밀턴-그린 솔버의 성능은 어떠한가?
• RQ3 한정된 시야 및 동적 PAT 환경에서 k-Wave와 같은 전체 파동 해법 대비 계산 비용을 얼마나 줄일 수 있으며, 정확도를 유지할 수 있는가?
• RQ4 공명면에서 레이 진폭을 정규화함으로써 센서에서 재구성된 압력 신호의 정확도에 어떤 영향을 미치는가?
• RQ5 효율적인 센서 기반 계산과 정확한 수반 맵핑을 고려할 때, 반복 재구성 프레임워크에서 이 솔버를 효과적으로 사용할 수 있는가?

주요 결과

• 해밀턴-그린 솔버는 음속 비균일성과 렌즈에 의한 레이 집중 및 공명면 형성으로 인해 발생하는 2차원 영역에서 정확한 정방향 및 수반 해를 도출한다.
• 공명면 이전에 위치한 초기 압력에 대해서는 k-Wave 해와 매우 유사한 결과를 보이며, 최소한의 오차를 기록한다.
• 공명면에서 가장 큰 오차가 발생하는데, 이는 |q(t; x₀)| < ε 인 경우 진폭을 0으로 설정함으로써 위상 연속성이 깨지기 때문이다.
• 공명면 이후에는 오차가 감소하지만 여전히 존재하며, 이는 공명면 형성과 관련된 복잡한 위상 이동 exp(−im(t)π/2)를 완전히 반영하지 못하기 때문이다.
• 레이의 집중과 렌즈 주변 굴절을 성공적으로 처리하며, 외곽 레이들은 공명면 영향을 받지 않으며, 진폭 정규화를 통해 수치적 안정성을 유지한다.
• 전체 파동 해법 대비 뚜렷한 계산적 이점과 함께, 센서 기반 독립 계산을 가능케 하여 비용의 일부분만으로도 수행할 수 있다.