[논문 리뷰] The Vertical Slice Transform in Spherical Tomography

Boris Rubin|arXiv (Cornell University)|2018. 07. 20.

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

한 줄 요약

이 논문은 구형 단층촬영에서 수직 슬라이스 변환에 대한 새로운 역공식과 특이값 분해를 제시하며, 열음향 및 광음향 영상에 적용 가능하다. 변환을 초평면 라돈 변환과 연결하고 해석적 계속성을 이용함으로써, 구 위의 함수를 수직 슬라이스를 통한 적분으로부터 명시적으로 재구성함으로써, 구 기하학에서 오일러-파아송-다브르아 방정식의 역문제를 해결한다.

ABSTRACT

The vertical slice transform takes a function on the n-dimensional unit sphere to integrals of that function over spherical slices parallel to the last coordinate axis. This transform arises in thermoacoustic tomography. We obtain new inversion formulas for the vertical slice transform and its singular value decomposition. The results can be applied to the inverse problem for the Euler-Poisson-Darboux equation associated to the corresponding spherical means.

연구 동기 및 목표

• n ≥ 2인 단위 구 $S^n$에서 구형 단층촬영의 수직 슬라이스 변환에 대한 역공식을 개발하는 것.
• 마지막 좌표축에 평행한 수직 초평면 절단을 통한 적분으로부터 함수를 재구성하는 역문제를 다루는 것.
• 좌표변환과 표면 측도 조정을 통해 수직 슬라이스 변환을 초평면 라돈 변환과 연결하는 것.
• 해석적 계속성과 변환에 대한 특이값 분해를 이용한 명시적 역공식 제공.
• 결과를 구 위의 오일러-파아송-다브르아 방정식의 역문제에 적용하여, 등자선 측정값으로부터 초기 자료를 재구성하는 것.

제안 방법

• 수직 슬라이스 변환 $Vf$는 함수 $f$를 마지막 좌표축에 평행한 구형 슬라이스에 대해 통합하는 것으로 정의된다.
• f가 마지막 변수에 대해 짝수라는 성질을 이용하여 문제를 상반구로 축소함으로써, 상반구 변환 $V_+f$를 도출한다.
• 좌표변환을 통해 구 적분을 $\mathbb{R}^n$에서의 라돈 변환으로 매핑하며, $\varphi(x') = f(x', \sqrt{1-|x'|^2}) / \sqrt{1-|x'|^2}$의 관계를 사용한다.
• 역변환은 라돈 변환의 역함수 $R^{-1}$를 통해 구성되며, $f(x', x_{n+1}) = x_{n+1} \cdot (R^{-1}\Phi)(x')$로 표현되며, $\Phi(\theta,t) = (V_+f)(\theta,t)/\sqrt{1-t^2}$이다.
• n=2인 경우, 해석적 계속성과 로그 잠재력에 기반한 대안적 방법이 개발되어, $\Delta$와 $\log|t - \theta \cdot x'|$을 포함하는 적분 공식을 도출한다.
• 명시적인 역변환은 $S^1 \times [-1,1]$ 위의 이重적분을 통해 이루어지며, 로그 커널과 $\sqrt{1-t^2}$에 의한 정규화를 포함한다.

실험 결과

연구 질문

• RQ1 n ≥ 2인 $S^n$에서 수직 슬라이스 변환은 유클리드 공간의 알려진 변환을 통해 어떻게 역공식화될 수 있는가?
• RQ2 $\mathbb{R}^n$에서 수직 슬라이스 변환과 초평면 라돈 변환 간의 관계는 무엇인가?
• RQ3 n=2인 경우, 해석적 계속성과 잠재력 이론을 통해 수직 슬라이스 변환의 역함수를 구성할 수 있는가?
• RQ4 수직 슬라이스 변환의 특이값 분해는 그 역함수 성질과 어떻게 관련되는가?
• RQ5 수직 슬라이스 변환과 구 위의 오일러-파아송-다브르아 방정식 간의 연결 고리는 무엇인가?

주요 결과

• 수직 슬라이스 변환 $Vf$는 $(V_+f)(\theta,t) = \sqrt{1-t^2} \cdot (R\varphi)(\theta,t)$의 관계를 통해 초평면 라돈 변환 $R\varphi$와 동치임을 보였다. 여기서 $\varphi$는 $f$의 좌표변환에 의해 유도된다.
• n ≥ 2인 경우, $f(x', x_{n+1}) = x_{n+1} \cdot (R^{-1}\Phi)(x')$로 표현되는 일반적인 역공식을 통해 $Vf$로부터 함수 $f$를 재구성할 수 있으며, $\Phi(\theta,t) = (V_+f)(\theta,t)/\sqrt{1-t^2}$이다.
• n=2인 경우, 대안적 역공식이 유도되었다: $f(x) = \frac{|x_3|}{8\pi^2} \Delta \int_{S^1} d\sigma(\theta) \int_{-1}^1 (Vf)(\theta,t) \frac{\log|t - \theta \cdot x'|}{\sqrt{1-t^2}} dt$로, 로그 잠재력과 해석적 계속성을 이용한다.
• 수직 슬라이스 변환의 특이값 분해가 명시적으로 계산되었으며, 스펙트럼 구조와 안정성 특성을 드러냈다.
• 결과는 오일러-파아송-다브르아 방정식의 역문제에 적용되었으며, $S^n$ 위의 초기 자료 $f$는 등자선 $S^{n-1}$에 제한된 해 $u(x,\omega)$로부터 재구성될 수 있다. 특히 $\alpha = (1-n)/2$인 경우는 파동 방정식에 해당한다.
• 이 방법은 고정된 반경 $0 < \theta \leq \pi/2$를 가진 지오데식 구의 역문제에 대해 구조적 해법을 제공하며, 이는 이전 결과를 더 넓은 범주로 확장한다.