[논문 리뷰] Localized deconvolution on the sphere

Dominique Picard|arXiv (Cornell University)|2009. 08. 13.

Numerical methods in inverse problems참고 문헌 15인용 수 6

한 줄 요약

이 논문은 구면에서 국소화된 디컨볼루션 방법을 제안하며, 특이값 분해(SVD) 역행렬과 두 번째 세대 웨이블릿 기저에서의 임계처리를 결합하여 L^p 손실에 대해 적응형 최소최대 비율을 달성한다. 이 방법은 구면 조화함수의 전역적 지지도 제한을 극복하여 다양한 부드러움과 희소성 조건 하에서 최적의 성능을 발휘하는 국소화된 재구성을 가능하게 한다.

ABSTRACT

AbstractWe provide a new algorithm for the treatment of deconvolutionon the sphere which combines the traditional SVD inversion withan appropriate thresholding technique in a well chosen new basis.We establish upper bounds for the behaviour of our procedure forany L p loss. It is important to emphasize the adaptation propertiesof our procedures with respect to the regularity (sparsity) of theobject to recover as well as to inhomogeneous smoothness. We alsoperform a numerical study which proves that the procedure showsvery promising properties in practice as well. Key words and phrases: statistical inverse problems, minimax estima-tion, second- generation waveletsMSC 2000 Subject Classi cation 62G05 62G08 62G20 62C10 1 Introduction The spherical deconvolution problem was rst proposed by Rooij and Ruym-gaart (1991) [14] and subsequently solved in Healy et al. (1998) [3]. Kimand Koo (2002) [7] established minimaxity for the L 2 -rate of convergence.The optimal procedures obtained there are using orthogonal series meth-ods associated with spherical harmonics. One important problem arisingwith these procedures is their poor local performances due to the fact thatspherical harmonics are spread all over the sphere. This explains for in-stance the fact that although they are optimal in the L

연구 동기 및 목표

• 구면 조화함수 기반 전통적 디컨볼루션 방법의 전역적 지지도로 인한 열악한 국소 성능 문제를 해결하기 위해.
• 구면 상의 기저 함수의 정규성과 희소성에 적응하는 디컨볼루션 절차를 개발하기 위해.
• 비균일한 부드러움 조건 하에서 L^p 손실에 대해 최소최대 최적 수렴 속도를 달성하기 위해.
• 구면 역문제에 대해 수치적으로 안정적이고 실용적으로 효과적인 방법을 제공하기 위해.

제안 방법

• 이 방법은 두 번째 세대 웨이블릿 기저를 사용하여 구면에서 재구성 과정을 국소화한다.
• 웨이블릿 도메인에서의 SVD 역행렬과 임계처리를 결합하여 해를 안정화하고 노이즈 증폭을 감소시킨다.
• 목표 함수의 다양한 정도의 부드러움과 희소성에 적응하도록 절차를 설계한다.
• 웨이블릿 기반 함수 표현과 임계처리 규칙을 사용하여 L^p 위험에 대한 이론적 상한을 도출한다.
• 두 번째 세대 웨이블릿의 다중해상도 구조를 활용하여 효율적인 계산과 국소화를 실현한다.
• 실용적 성능과 강인성을 검증하기 위해 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 구면에서의 디컨볼루션 방법이 L^p 손실 하에서 최소최대 최적 비율을 달성하면서도 국소화를 유지할 수 있는가?
• RQ2 제안된 방법은 복구할 함수의 비균일한 부드러움과 희소성에 어떻게 적응하는가?
• RQ3 기존의 구면 조화함수 기반 접근법과 비교해 실용적으로 이 방법의 성능은 어떠한가?
• RQ4 국소화된 웨이블릿 기저에서의 임계처리는 구면 디컨볼루션에서 안정성과 정확도를 향상시킬 수 있는가?

주요 결과

• 제안된 방법은 L^p 손실에 대해 최소최대 최적 수렴 속도를 달성하며, 이론적 하한값과 일치한다.
• 이 방법은 복구할 함수의 희소성과 비균일한 부드러움에 대해 강력한 적응 성질을 보인다.
• 두 번째 세대 웨이블릿의 사용은 구면 조화함수의 전역적 지지도 제한을 극복하는 국소화된 재구성을 가능하게 한다.
• 수치적 연구는 이론적 성능이 유한 표본 설정에서 실용적으로 효과적이고 강인함을 확인한다.
• 임계처리를 통한 안정성 유지로 인해 악조건의 역문제에서 노이즈 증폭을 감소시킨다.