[논문 리뷰] Regularized Newton methods for simultaneous Radon inversion and phase retrieval in phase contrast tomography

Simon Maretzke|arXiv (Cornell University)|2015. 02. 17.

Advanced X-ray Imaging Techniques인용 수 5

한 줄 요약

이 논문은 근거리 및 원거리 전파 모델을 활용하여 단일 거리에서의 복소 굴절률과 투과도를 동시에 복원하는 정규화된 뉴턴 방법을 제안한다. 이는 생물학적 세포와 같은 약간 투과성이 낮은 샘플의 3D 納米 스케일 영상에서 안정성과 정확도를 크게 향상시키는 데 기여한다.

ABSTRACT

Promoted by the advent of coherent synchrotron light sources, phase contrast tomography allows to resolve three-dimensional variations of an unknown sample's complex refractive index from scattering intensities recorded at different incident angles of an X-ray beam. By diffractive free-space propagation of the transmitted wave field, this method is sensitive not only to absorption but also to refractive phase shifts induced by the specimen, permitting three-dimensional nanoscale imaging of quasi-transparent samples such as biological cells. However, the reconstruction of the specimen structure from the observed data constitutes an algorithmically challenging nonlinear ill-posed inverse problem, mainly due to the characteristic loss of phase information in the detection of the wave field. In this work, regularized Newton methods are developed for the solution of this tomographic phase retrieval problem, based on a detailed analysis of its mathematical structure. We consider both the near-field- or Fresnel regime characterized by a moderate propagation length between sample and detector and the far-field limit of large detector distances, where propagation is governed by the Fourier transform. In the former setting, excellent numerical reconstructions are obtained via the chosen Newton-type approach, supplemented by novel theoretical results stating that measurements from a single detector distance are sufficient to uniquely recover both refraction and absorption of a sample. The proposed algorithm simultaneously performs tomographic- and phase reconstruction, which is found to stabilize the latter by exploiting correlations between the diffraction patterns recorded under different incident angles.

연구 동기 및 목표

• 강도 측정에서 위상 정보가 누락되어 발생하는 비선형적이고 불안정한 역문제를 해결하기 위함.
• quasi-투명한 샘플의 복소 굴절률과 단층 영상 구조를 동시에 회복하는 안정적이고 효율적인 복원 방법을 개발하기 위함.
• 단일 검출기 거리로도 굴절성과 투과성 성분을 유일하게 복원할 수 있음을 증명하여, 이전에 다중 거리가 필요하다는 가정을 도전하기 위함.
• 다양한 입사 각도 간의 각도 상관관계를 활용하여 위상 복원을 안정화하고 복원 품질을 향상시키기 위함.
• 프레넬 및 원거리 전파 모델의 맥락에서 정규화된 뉴턴 방법에 대한 수학적으로 엄밀한 프레임워크를 제공하기 위함.

제안 방법

• 근거리 조건에서 파라엑셜 파동 방정식과 프레넬 전파를 사용하여 정방 모델을 수립하고, 자유 공간 전파 후 검출기에서 측정된 강도를 기반으로 한다.
• 뉴턴 유형의 반복적 역복원을 가능하게 하기 위해 정방 연산자의 프레셰 도함수를 유도하며, 국소 2차 수렴 보장.
• 불안정한 역문제의 해를 안정화하기 위해 티호노프 유형의 정규화를 적용하고, 굴절률의 매끄러움을 촉진하는 펜alty 항을 도입.
• 원통형 푸리에 및 라돈 변환을 사용하여 근거리 및 원거리 영역에서 정방 연산자의 수학적 구조를 분석.
• 첨부 연산자와 효율적인 수치 해석기법을 사용한 이산화된 뉴턴 방법을 구현하여 대규모 3D 복원을 수행.
• 다양한 입사 각도의 데이터를 통합하여 공동 복원 프레임워크를 구성함으로써 각도 상관관계를 활용하고 위상 복원의 안정성을 향상시킴.

실험 결과

연구 질문

• RQ1 단일 검출기 거리로도 복소 굴절률의 굴절성 및 투과성 성분을 유일하게 결정할 수 있는가?
• RQ2 순차적 접근 방식과 비교해 복합적인 단층 영상 및 위상 복원이 위상 복원의 안정성과 정확도를 어떻게 향상시키는가?
• RQ3 근거리 및 원거리 영역에서 정방 연산자의 수학적 구조는 어떠한가? 이는 뉴턴 유형의 수렴에 어떤 영향을 미치는가?
• RQ4 다른 입사 각도 간의 각도 상관관계는 위상 성분 복원에 얼마나 기여하는가?
• RQ5 정규화를 어떻게 최적의 방식으로 뉴턴 반복과 조합하여 비선형 역문제의 해를 안정화할 수 있는가?

주요 결과

• 정방 연산자의 단사성에 대한 이론적 분석을 통해 단일 검출기 거리로도 복소 굴절률의 굴절성 및 투과성 성분을 유일하게 복원할 수 있음이 입증됨.
• 제안된 정규화된 뉴턴 방법은 근거리 영역에서 뛰어난 수치적 복원 성능을 보이며, 약간 투과성이 낮은 샘플의 3D 영상에서 높은 정밀도와 안정성을 확보함.
• 다양한 입사 각도에서의 衍射 패atters 간의 상관관계를 활용함으로써 동시에 단층 영상 및 위상 복원이 이루어질 경우 위상 복원의 안정성이 크게 향상됨.
• 실제 노이즈 수준에서도 수렴이 이루어지며, 정규화가 잡음과 잡음 유도 잡음의 억제에 효과적으로 기여하고 해상도 향상에 기여함.
• 이론적 분석을 통해 근거리 영역에서 미약한 조건 하에서도 정방 연산자가 단사적임을 확인하여 해의 유일성에 대한 정당성 제공.
• 원거리 근사는 푸리에 변환에 의해 지배됨을 입증하여, 큰 검출기 거리에서 단순화된 정확한 모델을 제공함.