[논문 리뷰] A family of inversion formulas in Thermoacoustic Tomography

Linh V. Nguyen|ArXiv.org|2009. 02. 15.

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

한 줄 요약

이 논문은 일정한 음속을 가진 열음향단층촬영(TAT)에 대해 시간 영역 및 주파수 영역 공식 모두에 적용 가능한 폐쇄형 역공식의 통합된 가족을 제시한다. 하이켈 함수와 상호연결 연산자를 활용하여 기존의 필터링 배경화형 공식(피ン치, 패치, 레이크시, 쿤얀스키, 주-왕 포함)을 일반화하고 통합하며, 차원과 데이터 유형 간의 동치성에 대한 모호함을 해결한다.

ABSTRACT

We present a family of closed form inversion formulas in thermoacoustic tomography in the case of a constant sound speed. The formulas are presented in both time-domain and frequency-domain versions. As special cases, they imply most of the previously known filtered backprojection type formulas.

연구 동기 및 목표

• 열음향단층촬영에서 이전에 유도된 역공식 간의 동치성과 관계에 대해 해결되지 않은 질문을 해결하기 위해.
• 특히 피닉-패치-레이크시, 쿤얀스키, 주-왕의 공식을 포함한 기존의 폐쇄형 역공식들을 하나의 일관된 프레임워크로 통합하기 위해.
• 역공식의 적용 범위를 홀수 차원과 짝수 차원 모두로 확장하고, 유한 시간 데이터 복원까지 포함하기 위해.
• 푸리에 및 하이켈 변환을 통해 시간 영역 공식과 주파수 영역 공식 간의 연결 고리를 설정하기 위해.
• 기존에 딜리클레 경계 데이터에 국한된 연구를 넘어서, 노이만 형 경계 데이터(압력의 법선 도함수)로부터의 역공식 유도를 위해.

제안 방법

• 주파수 영역에서 열음향 데이터의 짝수 연장의 푸리에 변환인 $\hat{g}_0(y,\lambda)$를 사용하여 역공식의 가족을 유도한다.
• 첫 번째 종류의 하이켈 함수를 기반으로 한 커널 함수 $G(s,\lambda) = \frac{i}{4}\left(\frac{\lambda}{2\pi s}\right)^{\frac{n-2}{2}}H^{(1)}_{\frac{n-2}{2}}(\lambda s)$를 도입한다.
• 파동 방정식을 다르쿠스 방정식으로 변환하는 상호연결 연산자 $\mathcal{W}$를 적용하여 구면 평균 연산자로의 변환을 가능하게 한다.
• 측정된 데이터와 初기 압력 분포 간의 관계를 위해 구면 라돈 변환 $\mathcal{R}_S$ 와 구면 평균 연산자 $\mathcal{M}_S$ 를 사용한다.
• 커널 $K(x,z,\lambda)$의 대칭성과 $\mathcal{W}$ 및 그 수반 연산자 $\mathcal{W}^*$ 간의 이중성 관계를 활용하여 역공식 항등식을 도출한다.
• 부분 적분과 푸리에-베셀 변환의 역공식 기법을 적용하여 데이터 $g$ 로부터 초기 함수 $f(x)$ 를 복원한다.

실험 결과

연구 질문

• RQ1 짝수 차원 열음향단층촬영에서 유한 시간 공식과 무한 시간 공식 간의 관계는 무엇인가?
• RQ2 쿠난스키 및 주-왕의 공식과 같은 기존의 역공식들은 어떻게 하나의 수학적 프레임워크로 통합될 수 있는가?
• RQ3 동일한 공식 가족이 딜리클레 데이터 외에도 노이만 형 경계 데이터(압력의 법선 도함수)로부터도 $f(x)$ 를 복원할 수 있는가?
• RQ4 하이켈 함수와 상호연결 연산자 $\mathcal{W}$ 는 모든 차원에서 통합된 역공식을 구성하는 데 어떤 역할을 하는가?
• RQ5 연산자 $\mathcal{W}\frac{d}{ds}\mathcal{W}^*$ 가 항등원과 동치가 되는 조건은 무엇이며, 이는 역공식과 어떻게 관련되는가?

주요 결과

• 제안된 역공식 가족은 기존에 알려진 모든 필터링 배경화형 공식을 특수 케이스로 포함하며, 피닉, 패치, 레이크시(홀수 차원), 쿤얀스키(모든 차원), 주-왕(n=3)의 공식을 포함한다.
• 짝수 차원에서는 유한 시간 역공식을 제공하며, 적절한 조건 하에서 피닉, 할트마이어, 레이크시의 무한 시간 공식과 동치가 된다.
• 커널 $G(s,\lambda)$ 는 모든 차원 $n \geq 2$ 에서 유효한 역공식을 보장하며, 짝수/홀수에 따라 부드러움 특성이 다름: 홀수 $n$ 에서는 $C^\infty$, 짝수 $n$ 에서는 영점 외에서 $C^\infty$ 이고 영점에서 로그 특이성을 가짐.
• 상호연결 연산자 $\mathcal{W}$ 는 2차 도함수를 베셀 유형의 연산자로 변환하여, 구면 평균과 푸리에-베셀 변환을 통한 역공식 유도를 가능하게 한다.
• 커널 $K(x,z,\lambda)$ 의 대칭성은 $\int_{\mathbb{R}} K(x,z,\lambda) d\lambda = \int_{\mathbb{R}} K(z,x,\lambda) d\lambda$ 를 암시하며, 이는 역공식에서 사용된 이중성의 기초가 된다.
• 노이만 데이터에 대한 역공식이 도출됨: $f(x) = 2\int_S \mathcal{W}(\partial_{\nu_y}u)(y, |x-y|) d\sigma(y)$ 로서, 이는 이전 결과를 일반 $n$ 과 노이만 형 측정까지 확장한다.