[논문 리뷰] A dissipative time reversal technique for photo-acoustic tomography in a cavity

Linh V. Nguyen, Leonid Kunyansky|arXiv (Cornell University)|2015. 07. 06.

Photoacoustic and Ultrasonic Imaging참고 문헌 29인용 수 5

한 줄 요약

이 논문은 에너지 보존으로 인해 기존 시간 역전 기법이 실패하는 완전하게 반사하는 벽을 가진 캐비티에서 광음향 단층촬영(photo-acoustic tomography)을 위한 소산성 시간 역전 기법을 제안한다. 딜리클레 경계 조건을 정규 및 시간 도함수를 포함하는 비표준 조건으로 대체함으로써 오차 방정식은 소산성 경계 조건을 만족하게 되어, 측정 시간 T가 증가함에 따라 초기 압력 재구성의 지수적 수렴을 가능하게 한다. 이는 비반복적 및 반복적(네우만 급수) 구현 방식을 통해 전체 및 부분 데이터에 모두 적용 가능하다.

ABSTRACT

We consider the inverse source problem arising in thermo- and photo-acoustic tomography. It consists in reconstructing the initial pressure from the boundary measurements of the acoustic wave. Our goal is to extend versatile time reversal techniques to the case of perfectly reflecting boundary of the domain. Standard time reversal works only if the solution of the direct problem decays in time, which does not happen in the setup we consider. We thus propose a novel time reversal technique with a non-standard boundary condition. The error induced by this time reversal technique satisfies the wave equation with a dissipative boundary condition and, therefore, decays in time. For larger measurement times, this method yields a close approximation; for smaller times, the first approximation can be iteratively refined, resulting in a convergent Neumann series for the approximation.

연구 동기 및 목표

• 에너지가 보존되는 경계 조건을 가진 공진 캐비티 내에서 열 및 광음향 단층촬영(TAT/PAT)의 역원천 문제를 다루며, 기존 시간 역전 기법이 에너지 감쇠가 없어 실패하는 이유를 해결한다.
• 해결책의 감쇠에 의존하는 표준 시간 역전 기법이 에너지가 보존되는 환경에서는 적용 불가능한 한계를 극복한다.
• 측정 시간이 짧을 경우에도 안정적이고 수치적으로 실현 가능한 시간 역전 방법을 개발하여 오차 감쇠 및 수렴을 보장한다.
• 엄격한 고유값 조건이나 조화 확장을 요구하는 이전 방법의 제약을 넘어 전체 및 부분 경계 데이터 모두에서 정확한 재구성을 가능하게 한다.

제안 방법

• 해결책의 법선 및 시간 도함수를 조합한 비표준 경계 조건을 사용하는 새로운 시간 역전 알고리즘을 제안하여 오차가 소산성 파동 방정식을 만족하도록 보장한다.
• 오차 방정식이 소산성 경계 조건으로 인해 시간이 지남에 따라 감쇠되도록 설정함으로써 측정 시간 T가 증가함에 따라 수렴이 가능해진다.
• 비반복적 근사와 반복적(네우만 급수) 전략을 포함한 두 가지 재구성 전략을 도입한다.
• 이전 연구에서 고유값 기반 조건보다 더 덜 제약적인 일반화된 제어 조건(GCC) 하에서 이론적 수렴을 확립한다.
• 경계 데이터의 조화 확장을 계산할 필요를 피함으로써 이전 평균 시간 역전 방법에 비해 수치적 구현을 단순화한다.
• 원형 도메인에서 2차원으로 수치적 검증을 수행하여 전체 및 부분 데이터, 노이즈 유무 상황에서의 성능을 입증한다.

실험 결과

연구 질문

• RQ1 에너지가 보존되는 완전하게 반사하는 벽을 가진 캐비티에서 시간 역전 기법을 개발할 수 있는가? 이 경우 직접 문제의 해는 시간이 지남에 따라 감쇠되지 않는다.
• RQ2 노이즈 경계 조건으로 인해 파동 에너지가 보존되는 상황에서 시간 역전의 오차를 어떻게 제어하고 감쇠시킬 수 있는가?
• RQ3 오차 방정식이 소산 조건을 만족하고 시간이 지남에 따라 감쇠되도록 하기 위해 시간 역전에서 어떤 경계 조건을 사용할 수 있는가?
• RQ4 제안된 방법은 도메인의 제한적인 스펙트럼 조건 없이 전체 및 부분 데이터 문제 모두에서 수렴을 달성하는가?
• RQ5 일부 이전 방법에서 요구하는 바와 달리, 경계 데이터의 조화 확장을 계산하지 않아도 효율적으로 구현 가능한가?

주요 결과

• 제안된 소산성 시간 역전 방법은 에너지 감쇠가 없더라도 측정 시간 T가 증가함에 따라 재구성 오차가 지수적으로 0에 수렴함을 보였다.
• 전체 데이터의 경우 첫 번째 시뮬레이션에서 $L^2(Ω)$ 상대 오차는 1.1%였고, $T = 1.6$일 때 다섯 번의 반복 후 3.3%로 감소하였다 (비교 기준: $T(Ω,\Gamma) \approx 2.828$).
• 부분 데이터의 경우 첫 번째 시뮬레이션에서 $L^2(Ω)$ 상대 오차는 6.8%였고, $T = 3$일 때 다섯 번의 반복 후 5.4%로 감소하였다 (비교 기준: $T(Ω,\Gamma) \approx 5.6569$).
• 50%의 화이트 노이즈가 첨가된 상황에서도 방법은 강인했으며, 첫 번째 근사에서 전체 데이터의 경우 상대 오차는 19%, 부분 데이터의 경우 22%로 나타나 데이터 노이즈에 대해 낮은 민감도를 보였다.
• 일반화된 제어 조건(GCC) 하에서 이론적 수렴이 입증되었으며, 이는 이전 연구에서 사용된 고유값 기반 조건보다 더 덜 제약적인 것으로 확인되었다.
• 수치 결과는 이론적 하한값에 가까운 측정 시간에서도 알고리즘이 잘 작동함을 보여주어, 날카로운 시간 요구 조건이 실질적으로 너무 엄격하지 않음을 시사한다.