[논문 리뷰] Mathematical Challenges Arising in Thermoacoustic Tomography with Line Detectors

Markus Haltmeier, T C FIDLER|ArXiv.org|2006. 10. 04.

Photoacoustic and Ultrasonic Imaging참고 문헌 63인용 수 10

한 줄 요약

이 논문은 선형 검출기를 사용한 열음향 단층촬영에 대한 수학적 프레임워크를 수립하여 3D 이미지 재구성 문제를 선 위에서 측정된 경계 데이터를 가진 2차원 파동 방정식으로 환원한다. 데이터의 적분 가능성 문제에도 불구하고 빠르고 안정적인 수치 구현이 가능한 푸리에 기반의 해석적 재구성 공식을 유도하고 철저히 분석한다.

ABSTRACT

Thermoacoustic computed tomography (thermoacoustic CT) has the potential to become a mayor non-invasive medical imaging method. In this paper we derive a general mathematical framework of a novel measuring setup introduced in [P. Burgholzer, C. Hofer, G. Paltauf, M. Haltmeier, and O. Scherzer, "Thermoacoustic tomography with integrating area and line detectors", IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 52 (2005)], that uses line shaped detectors instead of the usual point like ones. We show that the three dimensional thermoacoustic imaging problem reduces to the mathematical problem of reconstructing the initial data of the two dimensional wave equation from boundary measurements of its solution. We derive and analyze an analytic reconstruction formula which allows for fast numerical implementation.

연구 동기 및 목표

• 점 검출기 대신 선형 검출기를 사용한 열음향 단층촬영을 위한 엄밀한 수학적 모델을 개발하기 위해.
• 고대역폭, 소형 구현에서의 기술적 제약으로 인해 기존의 점 검출기 근사 모델이 약한 점검을 겪는 문제를 해결하기 위해.
• 선 검출기 측정치를 이용하여 3차원 열음향 영상 문제를 선 위에서의 초기 데이터 재구성 문제로 환원하기 위해.
• 실제 구현에 적합한 수학적으로 타당하고 수치적으로 효율적인 재구성 공식을 제공하기 위해.
• 비절대 적분 가능 측정치를 포함한 현실적인 데이터 조건 하에서 재구성 공식의 수렴성과 안정성을 분석하기 위해.

제안 방법

• 열음향 과정을 초기 압력 분포 f(x)가 흡수된 에너지를 나타내는 3차원 파동 방정식으로 모델링한다.
• 선 검출기를 선분 위에서 압력장의 적분으로 모델링하여 기록 곡선 C_rec 상의 경계 측정치를 도출한다.
• 검출기 선에 수직인 평면으로 도메인을 잘라내어 3차원 전진 연산자를 선형 검출기 측정치를 기반으로 한 2차원 파동 방정식 문제의 집합으로 분해한다.
• 선 측정치로부터 2차원 파동 방정식의 초기 데이터를 위한 푸리에 기반 해석적 재구성 공식을 유도한다.
• 측정 데이터를 초기 압력 분포로 연결하기 위해 C-변환과 그 역변환을 적용하여 푸리에 분석을 수행한다.
• C-변환이 대칭적 성분을 가진 토폴로지적 등장사상으로서의 성질을 증명함으로써 재구성 공식의 수렴성과 안정성, 가역성을 보장한다.

실험 결과

연구 질문

• RQ1 선 검출기를 사용할 경우 3차원 열음향 영상 문제는 어떻게 2차원 파동 방정식 문제로 환원될 수 있는가?
• RQ2 선 검출기 열음향 단층촬영에서 전진 연산자의 수학적 구조는 무엇이며, 점 검출기 모델과 어떻게 다를 수 있는가?
• RQ3 측정 데이터가 절대 적분 가능하지 않은 경우에도 안정적이고 빠른 해석적 재구성 공식을 도출할 수 있는가?
• RQ4 2차원 파동 방정식 프레임워크 내에서 C-변환은 측정 데이터와 초기 압력 분포를 연결하는 데 어떤 역할을 하는가?
• RQ5 제안된 방법은 현실적인 물리적 측정 조건 하에서 수치적 안정성과 수렴성을 어떻게 확보하는가?

주요 결과

• 선 검출기를 사용한 3차원 열음향 영상 문제는 2차원 파동 방정식의 해를 선 위에 제한한 것에서 초기 데이터를 재구성하는 문제로 환원된다.
• 푸리에 기반의 해석적 재구성 공식이 도출되고 엄밀히 분석되어, 측정 데이터가 절대 적분 가능하지 않은 경우에도 수학적으로 타당함이 보장된다.
• C-변환이 대칭적 성분을 가진 토폴로지적 등장사상으로서의 성질을 증명하여 측정 연산자의 안정적 역행렬화가 가능함이 입증된다.
• 역 C-변환은 측정 데이터의 푸리에 변환과余弦 조절을 포함하는 이중 적분을 통해 명시적으로 기술된다.
• 재구성 공식은 주파수 도메인과 공간 도메인에서의 분리 가능성과 해석적 구조 덕분에 빠른 수치적 구현이 가능하다.
• 단일 축을 중심으로 선 검출기 어레이를 회전시킴으로써 3차원 이미지 재구성을 달성하며, 3차원 직접 역행렬화에 비해 감소된 수치적 복잡도를 유지한다.