Exact reconstruction formulas for a Radon transform over cones

기술 요약: 원뿔에 대한 라돈 변환의 정확한 재구성 공식

문제 정의
본 논문은 $R^{(p)}$로 표기되는 특정 변형된 라돈 변환인 원뿔형 라돈 변환(conical Radon transform)의 수학적 역변환을 다룬다. 이 변환은 상반 평면 $\mathbb{R}^{d-1} \times (0, \infty)$ 내에 지지 집합을 갖는 매끄러운 함수 $f$를 일방향 원뿔 곡면(one-sided conical surfaces)에 대해 적분한다. 이 원뿔들은 하이퍼플레인 $\mathbb{R}^{d-1} \times \{0\}$에 위치한 꼭짓점, 이 평면에 수직인 대칭축(y축과 정렬됨), 그리고 반개각(half-opening angle) $\theta \in (0, \pi/2)$를 특징으로 한다. 적분에는 원뿔의 생성선(generator)을 따른 거리 $s$에 대한 방사형 가중치 $1/s^p$가 포함된다.

이 변환을 연구하는 주요 동기는 컴프턴 산란(Compton scattering)을 고려하는 컴퓨터 단층 촬영(CT)에 있다. 이러한 영상 기법에서 산란된 광자의 검출은 광원의 위치를 선이나 평면이 아닌 원뿔의 표면으로 추적할 수 있게 한다. 3차원($p \in \{0, 2\}$)의 경우 특정 푸리에 재구성 공식이 존재하고 2차원의 경우 역투영(back-projection) 공식이 존재하지만, 임의의 공간 차원 $d \geq 2$ 및 임의의 가중치 매개변수 $p \in \mathbb{R}$에 대한 일반적인 사례는 기존에 연구된 바 없다.

방법론
저자 Markus Haltmeier는 푸리에 변환과 한켈(Hankel) 변환을 활용하는 조화 해석학 기법을 사용하여 역변환 공식을 도출한다. 방법론의 핵심은 고전적인 라돈 변환의 잘 알려진 푸리에 슬라이스 정리와 유사하게, 원뿔형 라돈 변환에 대한 "푸리에 슬라이스 항등식(Fourier slice identity)"을 확립하는 것이다.

도출 과정은 세 가지 논리적 단계로 진행된다:

1. 가중치가 없는 경우로의 축소: 임의의 가중치 $p$에 대한 변환 $R^{(p)}$가 $\cos(\theta)$와 $y^{-p}$를 포함하는 단순 곱셈 연산자를 통해 가중치가 없는 변환 $R^{(0)}$와 연관될 수 있음을 보여주는 보조정리를 확립한다. 이를 통해 저자는 $p=0$인 경우의 도출에 집중하고 이후 결과를 일반화할 수 있다.
2. 푸리에 슬라이스 항등식의 도출: 공간 변수 $u \in \mathbb{R}^{d-1}$에 대한 푸리에 변환을 적용하고 구 $S^{d-2}$ 상에서의 베셀 함수(Bessel function)의 적분 표현을 활용하여, 원래 함수 $f$의 푸리에-한켈 변환과 원뿔형 라돈 변환 데이터의 푸리에 변환을 연결하는 항등식을 도출한다.
3. 역투영을 통한 역변환: 도출된 푸리에 슬라이스 항등식을 사용하여 역 푸리에 변환과 역 한켈 변환을 적용함으로써 $f$를 역변환한다. 이 과정은 첫 번째 인수에 대한 리에스 포텐셜(Riesz potential) $I^{(1-d)}$의 적용으로 식별되는 필터링 과정을 포함하는 필터링된 역투영(filtered back-projection) 형태의 명시적 재구성 공식을 이끈다.

주요 기여 및 결과
본 논문은 $R^{(p)}f$로부터 $f(x, y)$를 복구하기 위한 두 가지 주요 명시적 재구성 공식을 제공한다:

1. 각도 매개변수화 (정리 1.1, 식 1.3):
   $$f(x, y) = \frac{y^p}{(2\pi)^{d-1}} \int_0^{\pi/2} \frac{1}{\cos(\theta)^{1+p}} \int_{S^{d-2}} \left( I^{(1-d)} R^{(p)} f \right)(x + y \tan(\theta) n, \theta) \, dn \, d\theta$$
   이 공식은 개각 $\theta$와 단위 구 위의 방향 $n$에 의해 매개변수화된, 재구성 점 $(x, y)$를 통과하는 모든 원뿔들을 적분한다.

2. 공간 매개변수화 (정리 1.1, 식 1.4):
   $$f(x, y) = \frac{1}{(2\pi)^{d-1}} \int_{\mathbb{R}^{d-1}} \frac{(|u-x|^2 + y^2)^{\frac{p-1}{2}}}{|u-x|^{d-2}} \left( I^{(1-d)} R^{(p)} f \right)\left(u, \arctan\frac{|u-x|}{y}\right) \, du$$
   이 형태는 밑면 평면 $\mathbb{R}^{d-1}$에 대해 적분하며, 꼭짓점이 $u$에 위치하는 원뿔들의 기여를 합산한다.

특정 차원:

• $d=2$인 경우: 리에스 포텐셜 $I^{-1}$은 미분과 힐베르트 변환($\partial_u H_u$)의 합성체에 해당한다.
• $d=3$인 경우: 리에스 포텐셜 $I^{-2}$는 음의 라플라시안($-\Delta_u$)에 해당한다. 결과적으로 3D 재구성 공식은 데이터에 라플라시안을 적용한 후 원뿔형 역투영을 수행하는 과정을 포함한다.

의의 및 주장
본 논문은 임의의 차원 $d \geq 2$ 및 임의의 실수 가중치 $p$에 대해 원뿔형 라돈 변환에 대한 최초의 명시적인 필터링된 역투영 형태의 재구성 공식을 제공한다고 주장한다.

• 국소성(Locality): 저자는 차원 $d$의 기우성(parity)에 따른 구분을 언급한다. 홀수 차원의 경우, 재구성은 **국소적(local)**이다. 즉, $f$의 값은 $(x, y)$를 통과하는 임의의 작은 이웃 영역 내의 원뿔들에 대한 적분에만 의존한다. 짝수 차원의 경우, 재구성은 **비국소적(non-local)**이며, 모든 원뿔 곡면에 대한 적분 정보를 필요로 한다. 이러한 거동은 고전적인 라돈 변환과 유사하다.
• 신규성: $d=2$에 대한 동등한 공식들은 다른 연구들에서 이미 도출되었으나, $d=3$ 및 그 이상의 차원에 대한 공식은 새로운 기여로 제시되었다.
• 실용성: 저자는 3D 재구성 공식이 특히 단순하며, 라플라시안 적용과 원뿔형 역투영으로 구성되어 있다고 주장한다. 저자는 수치적 구현이 고전적 및 구형 라돈 변환에서 확립된 방법들을 따라 직관적일 것이라고 제언하지만, 수치적 연구는 향후 과제로 남겨두었다.

이 연구는 컴프턴 카메라를 사용하는 방출 단층 촬영의 이론적 토대를 마련하며, 이러한 장치에서 생성되는 데이터를 역변환하기 위한 필요한 수학적 도구를 제공한다.