[논문 리뷰] Inversion of the Spherical Mean Transform with Sources on a Hyperplane

Aleksei Beltukov|ArXiv.org|2009. 10. 08.

Numerical methods in inverse problems참고 문헌 12인용 수 19

한 줄 요약

이 논문은 n차원 상반공간의 초평면 위에 있는 원천을 가진 구 평균 변환에 대한 새로운 역공식을 제시한다. 이는 구 평균 변환과 포물면 커널에 대한 컨볼루션을 연결하는 상호연결 연산자를 사용한다. 주요 기여는 푸리에 분석과 고전적 라돈 이론을 바탕으로 한 명시적이고 안정적인 재구성 공식으로, 모든 빠르게 감소하는 함수에 대해 유효하며 임의의 차원에서 적용 가능하다.

ABSTRACT

The object of this study is an integral operator $\mathcal{S}$ which averages functions in the Euclidean upper half-space $\mathbb{R}_{+}^{n}$ over the half-spheres centered on the topological boundary $\partial \mathbb{R}_{+}^{n}$. By generalizing Norton's approach to the inversion of arc means in the upper half-plane, we intertwine $\mathcal{S}$ with a convolution operator $\mathcal{P}$. The latter integrates functions in $\mathbb{R}^{n}$ over the translates of a paraboloid of revolution. Our main result is a set of inversion formulas for $\mathcal{P}$ and $\mathcal{S}$ derived using a combination of Fourier analysis and classical Radon theory. These formulas appear to be new and are suitable for practical reconstructions.

연구 동기 및 목표

• 상반공간 ℝⁿ⁺¹에서 초평면 위에 있는 원천을 가진 구 평균 변환에 대한 새로운 명시적 역공식을 개발한다.
• 2차원 경우의 Norton의 접근을 임의의 차원 n ≥ 1으로 일반화한다.
• 구 평균 변환과 고전적 라돈 변환 사이에 컨볼루션 연산자를 통해 직접적인 연결을 수립한다.
• 단층촬영 응용에 적합한 안정적인 백프로젝션 유형의 재구성 알고리즘을 제공한다.

제안 방법

• ℝⁿ⁺¹의 함수에 작용하는 컨볼루션 연산자 𝒫와 구 평균 변환 𝒮를 연결하는 두 개의 사상 쌍을 도입한다.
• 𝒫를 회전 포물면의 이동 병진에 대한 적분을 포함하는 컨볼루션 연산자로 정의한다.
• 푸리에 분석과 분포 이론을 사용하여 𝒫의 역함수 𝒫⁻¹을 유도한다.
• 𝒮⁻¹ = 𝒮* ∘ 𝒦_y ∘ 𝒮의 조합을 통해 구 평균 변환의 역함수를 수립한다. 여기서 𝒦_y는 고전적 Λ-연산자와 관련된 일변수 연산자이다.
• 고전적 라돈 역함수 이론을 활용하여 역공식의 구조와 타당성을 정당화한다.
• 미분기하학적 기법을 통해 공식을 검증한다. 이는 체적 형식의 당김과 야코비안 행렬의 행렬식 계산을 포함한다.

실험 결과

연구 질문

• RQ1 원천이 초평면 위에 있는 구 평균 변환은 함수 정의역에 강력한 해석적 제약 조건을 요구하지 않는 방법으로 역행할 수 있는가?
• RQ2 어떻게 구 평균 변환을 컨볼루션 연산자를 통해 고전적 라돈 변환과 연결할 수 있는가?
• RQ3 모든 차원 n ≥ 1에 대해 적용 가능한 명시적이고 안정적인 구 평균 변환의 역공식을 유도할 수 있는가?
• RQ4 역공식은 수치적 안정성을 향상시키기 위해 수반 연산자 𝒮*를 포함하는 백프로젝션 형태로 표현될 수 있는가?
• RQ5 특히 포물면 컨볼루션 커널 측면에서 역연산자의 정확한 기하학적 및 해석적 구조는 무엇인가?

주요 결과

• 논문은 스웨츠 공간 𝒮(ℝⁿ⁺¹)에 속하는 모든 빠르게 감소하는 함수에 대해 유효한 새로운 구 평균 변환 𝒮의 역공식을 도출한다. 이는 제한적인 해석적 가정 없이도 가능하다.
• 역함수는 𝒮⁻¹(f) = 𝒮* ∘ 𝒦_y ∘ 𝒮(f)로 표현되며, 여기서 𝒦_y는 고전적 Λ-연산자의 일변수 당김으로서 단순성과 안정성을 보장한다.
• 𝒮와 상호연결되는 컨볼루션 연산자 𝒫는 분포 이론을 통해 가역성이 입증되어, 𝒫⁻¹에 대한 명시적 공식을 도출할 수 있다.
• 체적 형식의 당김에 사용된 변환의 야코비안 행렬식은 1/(1+|z|²)⁽ⁿ⁺¹⁾/²로 계산되며, 이는 역함수에 대한 정확한 스케일링을 확인한다.
• 이 방법은 구 평균 변환과 고전적 라돈 변환 사이에 직접적이고 자연스러운 연결을 수립하며, 후자는 역함수의 이론적 기초를 제공한다.
• 구성적 구조와 푸리에 기반 수치 알고리즘과의 호환성 덕분에 제안된 역공식은 실용적인 단층촬영 재구성에 적합하다.