강도 간섭계와 프티코그래피를 결합한 새로운 영상 복원 기법

본 논문은 1950년대부터 별의 각지름을 측정해 온 강도 간섭계(Intensity Interferometry, II)의 한계를 극복하고, 복잡한 물체를 고해상도로 이미징할 수 있는 새로운 프레임워크를 제시한다. 전통적인 II는 두 검출기에서 측정된 순간 강도 곱의 평균, 즉 2차 상관 함수를 이용해 물체의 공간 주파수 정보를 얻는다. 이 상관 함수는 물체의 푸리에 변환의 절댓값 제곱에 해당하지만, 위상 정보가 소실돼 직접적인 역푸리에 변환으로는 원본 이미지를 복원할 수 없다. 기존에는 G‑S, ER, HIO와 같은 위상 회복 알고리즘을 적용했으나, 복잡한 구조물에 대해서는 수렴이 느리거나 지역 최소값에 머무르는 문제가 있었다. 저자들은 이러한 문제를 해결하기 위해 프티코그래피(ptychography) 반복 알고리즘을 II와 결합한다. 프티코그래피는 겹치는 조명을 이용해 물체의 여러 부분을 동시에 측정함으로써 데이터의 중복성을 확보하고, 이를 통해 위상 정보를 효율적으로 추정한다. 특히, 이 연구는 기존 프티코그래피가 요구하던 일관된 코히런트 광원 대신, 비상관성(thermal) 광원을 사용한다. 회전하는 지면 유리(GG)를 통과한 레이저 빔은 의사열광을 생성하고, 이는 Rayleigh 분포와 균일 위상(−π~π) 특성을 가진다. 이러한 광원은 스펙클 패턴을 형성하고, 각 스캔 위치에서 기록된 스펙클을 통해 자동상관 함수를 계산한다. 수학적 전개는 두 검출기 I₁(t₁,x₁), I₂(t₂,x₂)의 상관 ⟨I₁I₂⟩가 물체 O(α,β)의 푸리에 변환 |F{O}|²와 비례함을 보이며, 이는 Van Cittert‑Zernike 정리와 일치한다. 따라서 자동상관 함수의 제곱근을 모듈러스로 사용하고, 위상은 반복적으로 업데이트한다. 구체적인 절차는 다음과 같다. (1) 초기 객체 O(r)와 프로브 P(r−R)을 곱해 빛장 X(r,R)를 만든다. (2) X를 수치적으로 파동 전파시켜 CCD 평면의 푸리에 영역 F(k)로 변환하고, 그 절댓값을 자동상관 함수의 제곱근으로 교체한다. (3) 수정된 F(k)를 역전파해 객체 평면으로 되돌리고, 객체 추정값 Oₙ(r)를 업데이트한다. (4) 모든 스캔 위치에 대해 위 과정을 반복하고, 사전 정의된 반복 횟수에 도달하면 종료한다. 프티코그래피의 핵심 제약 중 하나는 프로브의 정확한 위치와 크기를 알아야 한다는 점이다. 실제 실험에서는 스캔 단계에서 10 %~50 % 수준의 이동 오차가 발생할 수 있다. 이를 보완하기 위해 저자들은 ‘느슨한 지원(loose support)’ 개념을 도입한다. 지원 영역을 의도적으로 넓혀(예: 0, 5, 15, 20 픽셀) 프로브 이동 오류에 대한 내성을 확보한다. 시뮬레이션과 실험 결과, 25 %·50 % 이동 오류에서도 지원을 확대하면 물체 형태가 명확히 복원되며, 상대 잔차는 0.03 수준으로 유지된다. 실험 설정은 다음과 같다. 457 nm 파장의 고체 레이저를 10배 증폭한 뒤, 회전 속도 0.6 rpm인 지면 유리 위에 비추어 의사열광을 만든다. 프로브는 직경 5 mm의 원형 구멍이며, X축으로 0.5 mm 간격, 총 16단계 스캔한다. 각 단계마다 1000장의 스펙클을 20 ms 노출로 기록한다. CCD는 물체 뒤 250 mm 지점에 배치되어 스펙클 패턴을 수집한다. 실험에서는 프로브를 물체 앞·뒤에 배치하는 두 가지 구성을 검증했으며, 장거리(수 미터) 관측에도 적용 가능함을 보였다. 알고리즘 성능을 기존 ER, HIO와 비교하면, ER은 1000회 반복 후에도 세부 구조가 흐릿하고, HIO는 전반적으로 블러된 결과를 보인다. 반면 프티코그래피 기반 방법은 10~20회 반복만으로도 물체의 세부가 뚜렷하게 복원된다. 또한, 연산 비용이 낮아 5~20 루프 내에 수렴하며, 실험에서는 11~21 루프로 충분히 재구성하였다. 결론적으로, 이 연구는 (1) 복잡한 물체에 대한 높은 복원 정확도, (2) 프로브 위치·크기 불확실성에 대한 강인성, (3) 낮은 연산 비용이라는 세 가지 장점을 제공한다. 이는 기존 체렌코프 텔레스코프 배열을 보완하거나, 실험실 수준에서 저비용 고해상도 이미징을 구현하는 데 유용할 것으로 기대된다. 향후 연구에서는 더 큰 스케일의 천문 관측, 다중 파장 확장, 그리고 실시간 복원 알고리즘 개발이 과제로 남는다.