[논문 리뷰] Numerical methods for phase retrieval

Eliyahu Osherovich|arXiv (Cornell University)|2012. 03. 11.

Advanced X-ray Imaging Techniques참고 문헌 52인용 수 14

한 줄 요약

이 논문은 단일 촬영 강도 측정에서 비편광 회절 영상(CDI)에서 미세해상도를 달성하기 위해 희박성 제약 조건과 근사 푸리에 위상 지식을 활용하는 새로운 최적화 기반 위상 복원 방법을 제안한다. 알려진 기저(예: 웨이브릿 또는 기울기 영역)에서 신호의 희박성에 대한 사전 지식을 활용함으로써, 기존의 회절 한계를 초월하는 고해상도 영상을 복원한다. 이는 Fienup의 HIO나 Gerchberg-Saxton과 같은 전통적 알고리즘보다 뚜렷한 성능 향상을 보인다.

ABSTRACT

In this work we consider the problem of reconstruction of a signal from the magnitude of its Fourier transform, also known as phase retrieval. The problem arises in many areas of astronomy, crystallography, optics, and coherent diffraction imaging (CDI). Our main goal is to develop an efficient reconstruction method based on continuous optimization techniques. Unlike current reconstruction methods, which are based on alternating projections, our approach leads to a much faster and more robust method. However, all previous attempts to employ continuous optimization methods, such as Newton-type algorithms, to the phase retrieval problem failed. In this work we provide an explanation for this failure, and based on this explanation we devise a sufficient condition that allows development of new reconstruction methods---approximately known Fourier phase. We demonstrate that a rough (up to $π/2$ radians) Fourier phase estimate practically guarantees successful reconstruction by any reasonable method. We also present a new reconstruction method whose reconstruction time is orders of magnitude faster than that of the current method-of-choice in phase retrieval---Hybrid Input-Output (HIO). Moreover, our method is capable of successful reconstruction even in the situations where HIO is known to fail. We also extended our method to other applications: Fourier domain holography, and interferometry. Additionally we developed a new sparsity-based method for sub-wavelength CDI. Using this method we demonstrated experimental resolution exceeding several times the physical limit imposed by the diffraction light properties (so called diffraction limit).

연구 동기 및 목표

• 단일 강도 측정만 가능할 때 발생하는 비편광 회절 영상(CDI)에서의 위상 복원 문제에 대한 근본적 과제를 해결하기 위해.
• 신호의 희박성에 대한 사전 지식을 활용하여 단일 촬영 실험에서 초해상도 영상 촬영을 가능하게 하기 위해.
• 하드웨어 수정 없이도 광학, 전자 및 기타 현미경 시스템에 일반적으로 적용 가능한 알고리즘을 개발하기 위해.
• 희박성이 알려진 기저에서만 초해상도를 가능하게 하는 핵심 요소임을 입증하기 위해, 비음성 또는 유한한 지지 영역 조건보다 더 높은 성능을 달성함을 보여주기 위해.
• 기존 영상 하드웨어와 호환되는 실시간이고 계산 효율적인 복원 프레임워크를 제공하기 위해.

제안 방법

• 복소수 기울기와 헤시안을 사용하여 위상 복원 문제를 제약 조건이 있는 최적화 문제로 수식화한다.
• 목적 함수에 희박성 촉진 정규화 항을 도입하며, l0 유사 페널티 또는 다른 노름을 사용한다.
• 푸리에 도메인에서의 데이터 일치성과 신호 도메인에서의 희박성 강화를 번갈아 적용하는 반복 최적화 기법을 사용한다.
• 수렴 안정성과 속도 향상을 위해 근사 푸리에 위상 지식을 사전 지식으로 통합한다.
• 신호에 독립적인 기저 함수 선택 전략과 자동 그리드 결정 방법을 사용하여 복원 품질을 최적화한다.
• 기준 광선 위상 불확실성을 모델링하고 반복적으로 보정함으로써 허브로그래픽 영상에 적용한다.

실험 결과

연구 질문

• RQ1 알려진 기저에서의 희박성이 단일 촬영 강도 측정에서의 위상 복원에서 초해상도를 가능하게 할 수 있는가?
• RQ2 근사 푸리에 위상 지식은 위상 복원에서 복원 안정성과 정확도를 어떻게 향상시키는가?
• RQ3 기저 함수 선택과 그리드 해상도는 복원 품질과 수렴 속도에 어떤 영향을 미치는가?
• RQ4 제안된 방법이 Fienup의 하이브리드 입력-출력(HIO) 및 Gerchberg-Saxton 알고리즘보다 해상도와 내성 면에서 뛰어나게 성능을 발휘할 수 있는가?
• RQ5 이 방법은 원자력 또는 자기력 현미경과 같은 비광학적 영상 시스템으로까지 일반화될 수 있는가?

주요 결과

• 이 방법은 단지 강도 측정과 희박성에 대한 사전 지식만을 사용하여 CDI에서 초해상도를 달성하며, 기존의 디컨볼루션 및 위상 복원 기법보다 뚜렷한 성능 향상을 보였다.
• 희박성과 근사 위상 지식을 함께 사용할 경우 복원 품질이 크게 향상되며, 희박성만으로는 고해상도 복원이 불가능함을 확인했다.
• 처리 시간이 짧아 실시간 응용에 적합하며, Lena와 Phantom 등의 테스트 이미지에서의 속도 비교를 통해 높은 복원 속도를 확보했다.
• 다양한 그리드 피치 크기에서 목적 함수 값이 안정적으로 수렴하며, 정확한 그리드 해상도(예: 21 픽셀)에서 최적 성능를 보여, 파rameter 조정에 대한 강건성을 입증했다.
• 실험 결과는 소음 조건과 기준 광선의 위상 오차가 존재하는 상황에서도 합성 및 실제 CDI 데이터에서 초해상도 특징을 성공적으로 복원함을 보였다.
• 특히 저신호 대비 잡음 비율 및 좁은 지지 영역 조건에서 Fienup의 HIO 및 Gerchberg-Saxton 알고리즘보다 해상도와 내성 면에서 뛰어난 성능을 보였다.