희소 회절 데이터에서 위상 복원을 위한 볼록 최적화와 윈도잉 기법

** 본 논문은 회절 기반 이미징에서 위상 정보를 잃은 채 강도만 측정되는 전형적인 문제를 다루며, 특히 포아송 잡음이 지배적인 희소 데이터 상황에 초점을 맞춘다. 전통적인 위상 복원 기법은 과잉 샘플링(over‑sampling)과 물체 지원 영역(support) 제약을 결합해 비선형 최적화 문제를 풀지만, 이는 고신호‑대‑노이즈 비율에서만 안정적으로 동작한다. 현대 X‑ray 자유 전자 레이저(FEL) 실험에서는 단일 광자 검출기와 결함이 있는 검출기 배열 때문에 데이터가 매우 희소하고 잡음이 크게 섞여 있다. 이러한 환경에서는 기존 방법이 지역 최소에 머무르거나 전역 최적을 찾지 못한다는 한계가 있다. 저자들은 이 문제를 해결하기 위해 두 가지 주요 전략을 제시한다. 첫 번째는 지원 제약을 완화하고, 강도 자체를 변수로 삼아 자가상관(autocorrelation) 형태인 Patterson 함수를 직접 최적화한다는 점이다. 구체적으로, 원본 물체 P의 Fourier 변환 X와 강도 Y=|X|² 사이의 관계를 이용한다. Y의 Fourier 변환은 ℱ(Y)=ℱ(X)∗ℱ(X)ᶜ= P∗P 로 표현되며, 이는 물체와 그 자체의 컨볼루션이다. 기존의 지원 제약은 X에 직접 적용되었지만, 여기서는 Y에 대해 “지원의 자기 컨볼루션”을 적용함으로써 제약을 약화한다. 이는 선형 연산(ℱ)과 비선형 제약(지원) 사이의 복합성을 감소시켜, 최적화 문제를 볼록(convex) 형태로 변환한다. 두 번째 전략은 포아송 잡음 모델을 최대우도(MLE) 방식으로 포함하는 것이다. 실제 검출기에서 측정된 데이터 B는 포아송 분포 B_{ij}∼Poisson(r|X_{ij}|²) 로 모델링된다. 강도 Y를 복원할 때, 로그우도 −∑_{ij}