KR101479249B1 - 간섭성 구조조명 이미징 방법 및 간섭성 구조조명 현미경 시스템 - Google Patents
간섭성 구조조명 이미징 방법 및 간섭성 구조조명 현미경 시스템 Download PDFInfo
- Publication number
- KR101479249B1 KR101479249B1 KR20130135555A KR20130135555A KR101479249B1 KR 101479249 B1 KR101479249 B1 KR 101479249B1 KR 20130135555 A KR20130135555 A KR 20130135555A KR 20130135555 A KR20130135555 A KR 20130135555A KR 101479249 B1 KR101479249 B1 KR 101479249B1
- Authority
- KR
- South Korea
- Prior art keywords
- image
- coherent
- hexagonal lattice
- illumination
- amplitude
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/02—Interferometers
- G01B9/02041—Interferometers characterised by particular imaging or detection techniques
- G01B9/02044—Imaging in the frequency domain, e.g. by using a spectrometer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/25—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures by projecting a pattern, e.g. one or more lines, moiré fringes on the object
- G01B11/2518—Projection by scanning of the object
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B9/00—Measuring instruments characterised by the use of optical techniques
- G01B9/04—Measuring microscopes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N21/6402—Atomic fluorescence; Laser induced fluorescence
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/32—Micromanipulators structurally combined with microscopes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Optics & Photonics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Microscoopes, Condenser (AREA)
Abstract
본 발명은 간섭성 구조조명 이미징 방법 및 간섭성 구조조명 현미경 시스템을 제공한다. 이 간섭성 구조조명 이미징 방법은 육각형 격자 패턴에 의하여 형성된 육각형 격자 진폭 구조조명을 측정 대상에 차례로 조사하는 단계; 상기 육각형 격자 진폭 구조조명들에 의하여 주파수 하향 변환된 원시 이미지들을 차례로 측정하는 단계; 및 상기 원시 이미지들을 처리하여 재구성된 물체 이미지를 산출하는 단계를 포함한다.
Description
본 발명은 초고분해능 구조조명 이미징 방법에 관한 것으로, 더 구체적으로 2차원공간분해능 증대 효과 및 등방성 향상을 위한 간섭성 구조조명 이미징 방법에 관한 것이다.
최근 고분해능 광학 이미징 기술의 유용성과 발전 가능성에 대한 인식이 점차 증대되고 있다. 초고분해능을 구현하는 방법으로 STED(stimulated emission depletion), PALM(photoactivation localization microscopy), 그리고 fPALM(fluorescence PALM) 등과 같이 형광 여기 검출법이 사용된다. 또는, SIM(structured illumination microscopy) 기술은 패턴 조명을 이용하여 높은 분해능을 제공할 수 있다.
광 회절한계를 뛰어넘는 초고분해능 광학 이미징의 실현은 종래의 광학 현미경으로 관찰할 수 없었던 100 nm 이하의 크기를 갖는 단백질 구조, 지질 분포, DNA, RNA, 그리고 염색체 등의 세포 구성물과 세포 소기관을 실시간으로 가시화할 수 있다.
현재까지 제안된 초고분해능 이미징 기법을 사용하면, 30 nm 수준의 공간 분해능은 원격장(far-field) 광학 이미징으로 가능하다. 그러나, 이러한 방법은 대부분 특수하게 고안된 형광체에 의존하여 이미지 컨트라스트(image contrast)를 구현하고, 과도한 측정시간(수 분 - 수 시간)을 요구한다.
한편, 구조조명 현미경(structured illumination microscopy: SIM)은 회절한계 공간주파수 대역폭 너머에 존재하는 시료의 높은 공간주파수(초고분해능) 성분을 낮은 공간 주파수 성분으로 변환한다. 이에 따라, 상기 구조조명 현미경은 일반 현미경으로 측정할 수 없는 초고분해능 시료를 측정 가능한 상태로 만든다. 이후, 이미지 복원처리과정(reconstruction process)을 통해 높은 분해능의 이미지가 얻어질 수 있다. 여기서, 무아레(Moire) 효과를 이용한 이미지 공간 주파수 하향변환 (image spatial frequency down conversion)을 위해서, 상기 구조조명 현미경에는 특정한 격자 패턴의 구조조명이 요구된다. 전시야(full-field) 방식의 구조조명 현미경(SIM) 기술은 레이저 이미지 스캔 또는 긴 시간이 소요되는 통계적 이미지 프로세싱 등의 과정을 요구하지 않는다. 따라서, 전시야 (full-field) 방식의 구조조명 현미경(SIM) 기술은 실시간 측정의 관점에서 장점을 가진다.
상기 구조조명을 이용한 방법은 원리적으로 형광 시료 뿐만 아니라 비형광 시료에 대해서도 초고분해능 이미징을 수행할 수 있다. 그러나, 비형광 시료에 대한 초고분해능 이미징은 장치 구현과 이미지 데이터 처리 방법이 용이하지 않아서 최근까지도 정립된 기술이 없는 상황이다. 간섭성이 없는 조명(illumination) 또는 형광체의 방출광(emission light)을 이용하는 기존의 구조조명 이미징 방법은 비교적 간단한 장치 구현과 광학 이론체계에 바탕을 둔다. 그러나, 레이저 광원으로 조명되는 투과/반사/산란체의 이미지 형성(image formation)은 수학적으로 보다 복잡한 간섭성 결상 이론(coherent imaging theory)에 의해 기술된다. 또한, 초고분해능 이미지 재구성에 필요한 구조조명 이미지 취득 방법과 데이터 처리 알고리즘이 종래의 방법으로는 적용이 불가능하다.
본 발명의 일 실시예에 따르면, 간섭성(coherent) 결상(image formation)원리에 기반한 비형광 2차원 구조조명 이미지 취득법 및 초고분해능 데이터 처리 알고리즘 체계가 제공된다. 본 발명의 일 실시예에 따르면, 종래에 제시된 직교 주기 격자 (orthogonal periodic pattern) 간섭성 구조조명 이미징 방법이 가지는 공간분해능 향상 효과의 비등방성(anisotropy) 문제를 개선할 수 있다.
본 발명의 해결하고자 하는 일 기술적 과제는 종래의 비간섭성(incoherent) 구조조명 이미징 방법을 적용할 수 없는 비형광(투과/반사/산란) 현미경 시료에 대한 2차원 초고분해능 간섭성 구조조명 이미징 기술을 제공한다. 준등방성(quasi-isotropic) 공간분해능 향상 효과가 실현될 수 있다.
본 발명의 일 실시예에 따른 간섭성 구조조명 이미징 방법은 육각형 격자 패턴에 의하여 형성된 육각형 격자 진폭 구조조명을 측정 대상에 차례로 조사하는 단계; 상기 육각형 격자 진폭 구조조명들에 의하여 주파수 하향 변환된 원시 이미지들을 차례로 측정하는 단계; 및 상기 원시 이미지들을 처리하여 재구성된 물체 이미지를 산출하는 단계를 포함한다.
본 발명의 일 실시예에 있어서, 상기 원시 이미지들을 처리하여 재구성된 물체 이미지를 산출하는 단계는 상기 육각형 격자 진폭 구조조명의 3 방위 진폭 조합에 의하여 정의되는 스펙트럼 혼합 행렬의 역행렬을 구하는 단계; 상기 원시 이미지들을 푸리에 변환한 원시 이미지 스펙트럼을 산출하는 단계; 상기 원시 이미지 스펙트럼으로 구성된 스펙트럼 열 벡터와 상기 스펙트럼 혼합 행열의 역행렬을 이용하여 상기 원시 이미지 스펙트럼을 구성하는 측정 스펙트럼 구성항을 추출하는 단계; 상기 측정 스펙트럼 구성항을 공간 주파수 되돌림하여 가상 스펙트럼 구성항을 산출하는 단계; 상기 가상 스펙트럼 구성항을 선형 결합하여 가상 이미지 스펙트럼을 형성하는 단계; 및 상기 가상 이미지 공간 주파수 스펙트럼을 푸리에 역변환하여 물체 이미지 또는 가상 이미지를 산출하는 단계를 포함할 수 있다.
본 발명의 일 실시예에 있어서, 상기 육각형 격자 패턴은 3 방위각 복소 진폭 간섭성 구조조명 패턴으로, 2차원 공간의 전체 방위각을 3 등분한 {0o, 120o, 240o} 방향으로 격자 벡터를 갖는 3 개의 선형 격자 패턴을 진폭 합산(amplitude summation)하여 얻어질 수 있다.
본 발명의 일 실시예에 있어서, 상기 육각형 격자 패턴은 3 방위 위상 조합에 의하여 형성되고, 상기 3 방위 위상 조합은 3축 위상 공간의 19 개 위상 벡터 점으로 표시될 수 있다.
본 발명의 일 실시예에 따른 간섭성 구조조명 현미경 시스템은 간섭성 광을 출력하는 광원부; 상기 광원부로부터 상기 간섭성 광을 제공받아 육각형 격자 패턴을 가지는 육각형 격자 진폭 구조조명를 생성하는 구조조명 패턴 형성부; 상기 육각형 격자 진폭 구조조명을 제공받아 측정 대상에 제공하는 현미경 광학계; 및 상기 육각형 격자 진폭 구조조명을 제공받은 측정 대상이 반사 또는 투과시키는 광을 측정하는 측정부를 포함한다.
본 발명의 일 실시예에 있어서, 상기 측정부는 상기 육각형 격자 진폭 구조조명에 의한 원시 이미지를 측정하고, 상기 원시 이미지를 처리하여 물체 이미지를 재구성하는 처리부를 더 포함할 수 있다.
본 발명의 일 실시예에 있어서, 상기 구조조명 패턴 형성부는 상기 측정 대상의 상 공액면(image conjugate plane)에 배치될 수 있다.
본 발명의 일 실시예에 있어서, 상기 구조조명 패턴 형성부는 2차원 회절 격자, 액정 디스플레이형 공간 광 변조기, 또는 디지털 거울 소자일 수 있다.
본 발명의 일 실시예에 있어서, 상기 측정부는 투과단 또는 반사단 방향에 설치된 경통 렌즈; 및 상기 현미경 광학계 및 상기 경통 렌즈를 통하여 형성된 확대 상(magnified image)을 2차원적으로 측정하는 광 검출기를 포함할 수 있다.
본 발명은 종래의 비간섭성(incoherent) 구조조명 이미징 측정법 및 데이터 처리 알고리즘의 적용이 불가능한 비형광(투과/반사/산란) 현미경 시료에 대한 2차원 초고분해능 간섭성 구조조명 이미징 기술을 제공한다. 본 발명의 일 실시예에 따른 간섭성 구조조명 이미징 방법은 2차원 공간분해능 향상 효과의 준등방성(quasi-isotropic)을 실현할 수 있다.
도 1은 비간섭성 구조조명에 의한 2차원 이미징 광학계의 광전달 함수(optical transfer function)(a)와 간섭성 전달 함수(coherent transfer function)(b)의 개념도이다.
도 2는 2차원 물체에 대한 광학 이미징 컴퓨터 시뮬레이션 결과를 나타내는 도면이다.
도 3은 비간섭성 선형격자패턴 구조조명(non-coherent linear grating pattern structured illumination)을 설명하는 도면이다.
도 4는 본 발명의 일 실시예에 따른 현미경 이미징 광학계를 설명하는 도면이다.
도 5는 본 발명의 일 실시예에 따른 등방성 2차원 간섭성 구조조명 이미징의 확장 간섭성 전달함수(extended coherent transfer function)를 설명하는 도면이다.
도 6은 본 발명의 일 실시예에 따른 3-방위각 간섭성 진폭 구조조명 육각형 격자 패턴을 설명하는 도면이다. (a) 선형 격자, (b) 육각형 격자, (c) 격자 위상 이동점.
도 7은 본 발명의 일 실시예에 따른 이미징 방법을 설명하는 흐름도이다.
도 8은 간섭성 구조조명 이미지로부터 스펙트럼 역혼합 계산을 통해 분리한 공간 주파수 하향 변환 측정 스펙트럼 구성항을 나타내는 도면이다.
도 9는 공간 주파수 하향 변환 측정 스펙트럼 구성항에 대한 공간 주파수 영역 재구성 복원 위치를 나타내는 도면이다.
도 10은 등방성 2차원 간섭성 구조조명 이미징에 의한 공간 주파수 스펙트럼 차단 주파수 확장 효과를 나타내는 도면이다.
도 11은 등방성 2차원 간섭성 구조조명 이미징에 실공간 이미지 공간 분해능 향상 효과를 나타내는 도면이다.
도 12는 종래의 직교 2-방위 구조조명 이미징과 본 발명의 등방성 3-방위 구조조명 이미징의 공간 분해능 특성 비교하는 도면이다.
도 2는 2차원 물체에 대한 광학 이미징 컴퓨터 시뮬레이션 결과를 나타내는 도면이다.
도 3은 비간섭성 선형격자패턴 구조조명(non-coherent linear grating pattern structured illumination)을 설명하는 도면이다.
도 4는 본 발명의 일 실시예에 따른 현미경 이미징 광학계를 설명하는 도면이다.
도 5는 본 발명의 일 실시예에 따른 등방성 2차원 간섭성 구조조명 이미징의 확장 간섭성 전달함수(extended coherent transfer function)를 설명하는 도면이다.
도 6은 본 발명의 일 실시예에 따른 3-방위각 간섭성 진폭 구조조명 육각형 격자 패턴을 설명하는 도면이다. (a) 선형 격자, (b) 육각형 격자, (c) 격자 위상 이동점.
도 7은 본 발명의 일 실시예에 따른 이미징 방법을 설명하는 흐름도이다.
도 8은 간섭성 구조조명 이미지로부터 스펙트럼 역혼합 계산을 통해 분리한 공간 주파수 하향 변환 측정 스펙트럼 구성항을 나타내는 도면이다.
도 9는 공간 주파수 하향 변환 측정 스펙트럼 구성항에 대한 공간 주파수 영역 재구성 복원 위치를 나타내는 도면이다.
도 10은 등방성 2차원 간섭성 구조조명 이미징에 의한 공간 주파수 스펙트럼 차단 주파수 확장 효과를 나타내는 도면이다.
도 11은 등방성 2차원 간섭성 구조조명 이미징에 실공간 이미지 공간 분해능 향상 효과를 나타내는 도면이다.
도 12는 종래의 직교 2-방위 구조조명 이미징과 본 발명의 등방성 3-방위 구조조명 이미징의 공간 분해능 특성 비교하는 도면이다.
이하, 첨부한 도면들을 참조하여 본 발명의 바람직한 실시예들을 상세히 설명하기로 한다. 그러나, 본 발명은 여기서 설명되어지는 실시예들에 한정되지 않고 다른 형태로 구체화될 수도 있다. 오히려, 여기서 소개되는 실시예는 개시된 내용이 철저하고 완전해질 수 있도록 그리고 당업자에게 본 발명의 사상이 충분히 전달될 수 있도록 하기 위해 제공되어지는 것이다. 도면들에 있어서, 구성요소는 명확성을 기하기 위하여 과장되어진 것이다. 명세서 전체에 걸쳐서 동일한 참조번호로 표시된 부분들은 동일한 구성요소들을 나타낸다.
[구조조명 이미징 기술의 원리]
광학 현미경에서의 공간분해능의 한계()는 이상적으로 무한히 작은 물체라 할지라도 회절현상에 의한 유한한 점퍼짐(point spreading)으로 나타난다. 파장()인 빛을 사용하는 수치구경(numerical aperture; )의 광학 시스템에서, 공간분해능의 한계()는 로 주어지는 아베 회절한계(Abbe diffraction limit) 값에 의해 결정된다.
물체의 광학적 강도 분포는 공간주파수(spatial frequency) 스펙트럼으로 변환될 수 있다. 이 경우, 이미징 시스템(또는 현미경 시스템)은 물체를 상면(image plane)으로 정보화하여 전달 가능한 차단주파수(cutoff frequency)를 가진다. 상기 차단 주파수는 전달 가능한 물체 공간주파수 대역의 상한선이다. 상기 차단 주파수(cutoff frequency)를 넘어서는 고주파수 성분의 물체 고유 정보는 이미징 과정에서 손실되는 것으로 이해할 수 있다.
도 1은 비간섭성 구조조명에 의한 2차원 이미징 광학계의 광전달 함수(optical transfer function)(a)와 간섭성 전달 함수(coherent transfer function)(b)의 개념도이다.
도 1의 (a)를 참조하면, 형광체와 같은 비간섭성 시료를 대상으로 하는 이미징 시스템의 공간주파수 응답(spatial frequency response)은광전달함수(optical transfer function: OTF)로 기술된다. 상기 광전달 함수(OTF)는 강도 차단 주파수 를 가진다. 즉, 강도 차단 주파수 를 반지름으로 하는 원 (점선으로 표시) 내부로 정의되는 2차원 공간주파수 영역이 비간섭성 이미징의 이론적 측정가능 대역(spatial frequency passband)이 된다. Fx는 x축 방향에 대응하는 제1 방향 공간 주파수이고, Fy는 y축 방향에 대응하는 제2 방향 공간 주파수이다.
도 1의 (b)를 참조하면, 레이저와 같이 간섭성을 갖는 광원을 사용하여 비형광(투과/반사/산란) 시료를 이미징하는 광학 시스템에서의 상 형성(image formation)은, 물체의 광 복소 진폭(complex amplitude) 분포와 이에 대한 공간주파수 응답(spatial frequency response)인 간섭성 전달 함수 (coherent transfer function: CTF)로 기술된다. 이 때, 상기 간섭성 전달 함수(CTF)는 진폭 차단 주파수 를 가진다. 따라서, 간섭성 이미징에 있어서의 시료 측정 가능 대역 또는 공간 주파수 투과대역(spatial frequency passband)은 진폭 차단 주파수 를 반지름으로 하는 원 (점선으로 표시) 내부로 정의되는 2차원 공간주파수 영역으로 제한된다. Fx는 제1 방향 공간 주파수이고, Fy는 제2 방향 공간 주파수이다.
광학 이미징 시스템은 물체 강도 분포 또는 물체 복소 진폭 분포의 공간주파수를 상면(image plane)으로 전달하는 과정에서 차단 주파수 한계 내의 물체 정보만을 전달한다. 따라서, 광학 이미징 시스템은 본래의 물체 강도 분포와 달리 선명도가 저하된 상(image)을 형성한다. 이 때, 동일한 광학 이미징 시스템이라 하더라도, 사용하는 광원과 관찰 대상 물체의 간섭성 유무에 따라 형성되는 상의 모양과 공간분해능 특성이 상이하게 나타난다.
도 2는 2차원 물체에 대한 광학 이미징 컴퓨터 시뮬레이션 결과를 나타내는 도면이다.
도 2를 참조하면, 동일한 시험 타켓에 대하여 비간섭성 이미징과 간섭성 이미징이 계산되었다.
도 2를 참조하면, 광원 및 물체의 간섭성 여부에 따라 이미징 시스템의 공간주파수 영역 특성이 각각 광전달 함수(OTF)또는 간섭성 전달 함수(CTF)로 다르게 표현된다. 뿐만 아니라, 결상 이론의 전개에 사용되는 수학적 체계가 각기 다르다. 따라서, 동일한 시험 타켓(중심 영상)에 대하여 비간섭성 이미징(좌측 영상)과 간섭성 이미징(우측 영상)이 서로 다르다.
먼저, 비간섭성 이미징의 공간분해능 향상을 위한 구조조명 방법이 도 1과 도 3을 참조하여 설명된다.
도 3은 비간섭성 선형격자패턴 구조조명(non-coherent linear grating pattern structured illumination)을 설명하는 도면이다.
도 3을 참조하면, 고주파수 성분을 포함하는 임의 형태의 강도 물체(intensity object)에 특정 주기(period) 의 선형격자패턴 구조조명(linear grating pattern structured illumination)을 도 3의 (a1) 또는 도 3의 (a2)와 같이 특정 방향으로 인가한다. 상기 선형격자패턴 구조조명은 물체의 공간주파수 스펙트럼을 상기 격자 패턴 주기의 역수 만큼 주파수 하향 변환(frequency down conversion)하여 유효 발광체 분포(effective emitter distribution)를 형성한다. 이에 따라, 구조조명 광학 시스템은 강도 차단 주파수 를 넘어서는 물체의 공간주파수 스펙트럼을 측정가능 대역 또는 공간주파수 투과대역(spatial frequency passband)으로 이동시키는 역할을 한다. 이후, 이미지 처리(image processing)과정은, 본래 물체의 투과대역 스펙트럼(passband spectrum)과 하향이동 스펙트럼(down-shifted spectrum)이 중첩된 상태로 측정된 이미지 데이터로부터 각 스펙트럼 성분을 분리하여 원상 복원한다. 결과적으로, 상기 광학 시스템의 본래의 차단 주파수(cutoff frequency)를 증가시키는 효과가 나타난다. 일례로, 도 3의 (a1)과 도 3의 (a2)에 나타낸 바와 같이, 직교한 2방향으로 선형격자 패턴조명을 인가하여 얻은 측정 이미지 데이터로부터 공간주파수 스펙트럼이 재구성되면, 도 1의 (a)에서 실선으로 도시한 외곽선까지 측정가능 대역(spatial frequency passband)이 확장된 물체 스펙트럼이 얻어진다. 따라서, 상기 확장된 물체 스펙트럼을 역 푸리에 변환(inverse Fourier transform)하면, 공간영역(spatial domain)에서의 이미지 공간분해능이 향상된다.
물체 스펙트럼에 대한 상기 비간섭성 구조조명에 의한 유효 발광체 스펙트럼의 공간주파수 이동량(shift)과 방향(direction)은, 격자 패턴 조명의 격자 주기(grating period)와 격자 벡터(grating vector) 방향에 의해 결정되며, 분해능 향상 정도와 유효 방향에 영향을 끼친다. 도 1의 (a) 참조를 참조하면, 비간섭성 선형광학 구조조명법에서 생성할 수 있는 광 격자 패턴의 최고 공간주파수 는 광전달 함수(OTF)의 상기 강도 차단 주파수 와 같다. 따라서, 공간분해능 향상의 효과는 최대 2배이다.
비간섭성 이미징의 경우, 측정 이미지 데이터의 강도 분포 함수는 선형적인 물리 법칙을 따른다. 따라서, 측정 이미지 데이터의 강도 분포 함수는 푸리에 변환(Fourier transform) 관계에 있는 공간주파수 스펙트럼으로 간편히 기술된다. 비간섭성 선형격자 패턴조명의 형성 및 이미지 데이터 측정은 구조조명 방위각 별로 독립적으로 시행된다. 이후, 개별 1차원 구조조명 이미지 처리 결과들은 모두 선형 중첩된다. 이에 따라, 2차원 초고분해능 구조조명 이미지가 복원될 수 있는 장점이 있다.
한편, 본 발명의 일 실시예에 따르면, 간섭성 구조조명 이미징에서의 초고분해능 구현 방법은 개념 상 앞서 기술한 비간섭성 구조조명 이미징 방법과 유사하게, 진폭 물체(amplitude object)의 높은 2차원 공간주파수 성분을 간섭성 구조조명에 의해 주파수 하향 변환(frequency down conversion)한다. 이에 따라, 진폭 차단 주파수 로 특정되는 측정가능 대역 또는 공간주파수 통과대역(spatial frequency passband)에서 이미지 데이터가 측정된다. 이후, 측정된 이미지 데이터는 물체 본래의 공간주파수 스펙트럼 성분을 분리/추출/재구성하는 초고분해능 영상 재구성 과정(super-resolution image reconstruction process)을 거쳐 초고분해능 공간 이미지로 복원된다.
간섭성 구조조명 이미징의 실제 구현은, 구조조명 패턴의 형성 방법과 초고분해능 공간 이미지 복원을 위한 광학 이론 및 계산 알고리즘에 있어서 비간섭성 구조조명 이미징과 매우 상이하다.
간섭성 이미징의 경우, 빛의 복소 진폭이 수학적 기술에 있어서 선형적인 물리량이다. 그러나, 일반적인 2차원 광 검출기에서 측정되는 광 강도는 간섭성 상 형성(coherent image formation)을 기술하는 수학적 체계에서 비선형적이다. 따라서, 선형광학에 근간한 기존의 이론적 취급에 수정이 가해져야 한다. 보다 실제적으로는, 독립적인 방위각 별 격자패턴으로 형성된 이미지 처리 결과는 2차원 구조조명 효과를 구현하기 위해 중첩되어 사용할 수 없는 난점이 존재한다.
간섭성 이미징에서의 2차원 구조조명 패턴 형성과 이미지 처리는 다음과 같다. 우선, 확장하고자 하는 공간주파수 스펙트럼 대역의 하향 변환 주파수에 부합하는 복수의 격자패턴 조명을 동시에 전사(projection)하여 원시 이미지 데이터가 얻어진다. 이후, 원시 이미지 데이터를 일괄적으로 취급하여, 목표로 하는 확장 간섭성 전달 함수(extended coherent transfer function)에 의한 이미징 결과와 동일한 효과를 복원해내는 해석적 이미지 처리 알고리즘을 탐색하여 적용해야 한다.
예를 들어, 도 1의 (b)에서 실선으로 도시한 외곽선까지 확장된 공간주파수 대역의 물체 스펙트럼이 포함된 구조조명 원시 이미지 데이터를 얻기 위한 2차원 간섭성 구조조명 이미징에 필요한 투사 조명 패턴이 요구된다. 상기 투사 조명 패턴은 도 3의 (b)에 나타낸 바와 같이 정방형 격자패턴(square lattice pattern)을 갖는 2차원 투사 조명일 수 있다. 상기 정방형 격자패턴 구조조명은 비간섭성 구조조명법에서 사용하는 직교한 2방위각으로의 개별 선형격자 패턴 조명을 동시에 중첩하였을 때 얻어지는 진폭 구조조명 형태이다.
도 3의 (b)를 참조하면, 정방형 격자 진폭 구조조명 이미징 기법에 의해 얻어지는 확장 간섭성 전달함수는, 도 1의 (b)에 나타낸 진폭 차단 주파수 경계의 형태에서 보듯이 상기 정방형 격자 진폭 구조조명을 구성하는 2개의 직교한 진폭선형격자(도 3의 (a1)과 도3의 (a2) 참조)의 벡터 방위각({0o, 90o, 180o, 270o})으로는 진폭 차단 주파수를 배가(doubling)시키는 효과가 있다. 그러나, 상기 선형격자 벡터의 방위각 중간에 위치한 진폭 물체의 공간주파수 성분은 진폭 차단 주파수 확장 효과가 감소한다. 예로써, 벡터 방위각({45o, 135o, 225o, 315o})에서 차단 주파수 확장 효과가 최소가 되며 약 1.4배에 그친다.
따라서, 모든 방위각에 대하여 분해능이 향상된 등방성 간섭성 구조조명 방법이 요구된다.
[2차원 등방성 초고분해능 간섭성 구조조명 이미징의 구현 방법]
본 발명의 일 실시예에 따른 2차원 등방성 초고분해능 간섭성 구조조명 이미징 기술이 이하 설명된다.
도 4는 본 발명의 일 실시예에 따른 현미경 이미징 광학계를 설명하는 도면이다.
도 5는 본 발명의 일 실시예에 따른 등방성 2차원 간섭성 구조조명 이미징의 확장 간섭성 전달함수(extended coherent transfer function)를 설명하는 도면이다.
도 6은 본 발명의 일 실시예에 따른 3-방위각 간섭성 진폭 구조조명 육각형 격자 패턴을 설명하는 도면이다. (a) 선형 격자, (b) 육각형 격자, (c) 격자 위상 이동점.
도 7은 본 발명의 일 실시예에 따른 이미징 방법을 설명하는 흐름도이다.
도 4를 참조하면, 상기 이미징 방법은 측정하고자 하는 진폭 물체(측정 대상)에 소정의 간섭성 구조조명 패턴을 투사하여 투과단 또는 반사단에서 전시야(full-field) 광 강도 분포를 이미징하는 현미경 이미징 광학계에서 구현될 수 있다.
본 발명의 일 실시예에 따른 간섭성 구조조명 현미경 시스템(100)은 간섭성 광을 출력하는 광원부(110), 상기 광원부(110)로부터 상기 간섭성 광을 제공받아 육각형 격자 패턴을 가지는 육각형 격자 진폭 구조조명를 생성하는 구조조명 패턴 형성부(126), 상기 육각형 격자 진폭 구조조명을 제공받아 측정 대상(20)에 제공하는 현미경 광학계(130), 및 상기 육각형 격자 진폭 구조조명을 제공받은 측정 대상이 반사 또는 투과시키는 광을 측정하는 측정부(140,150)를 포함한다.
상기 측정부(149,150)는 상기 육각형 격자 진폭 구조조명에 의한 원시 이미지를 측정할 수 있다. 처리부(160)는 상기 원시 이미지를 처리하여 상기 측정 대상의 물체 이미지를 재구성할 수 있다. 상기 구조조명 패턴 형성부(120)는 상기 측정 대상(20)의 상 공액면(image conjugate plane)에 배치될 수 있다. 상기 구조조명 패턴 형성부(120)는 2차원 회절 격자, 액정 디스플레이형 공간 광 변조기, 또는 디지털 거울 소자일 수 있다.
상기 측정부(140,150)는 투과단 또는 반사단 방향에 설치된 경통 렌즈(142,152), 및 상기 현미경 광학계(130) 및 상기 경통 렌즈(142,152)를 통하여 형성된 확대 상(magnified image)을 2차원적으로 측정하는 광 검출기(144,154)를 포함할 수 있다.
상기 육각형 격자 패턴을 가지는 육각형 격자 진폭 구조조명은 구조조명 패턴 형성부(120)에 의하여 형성될 수 있다. 상기 구조조명 패턴 형성부(120)는 2차원회절격자(2D diffraction grating), 액정 디스플레이(liquid crystal display; LCD)형 공간 광 변조기(spatial light modulator: SLM), 또는 디지털 거울 소자 (digital mirror device: DMD)일 수 있다. 상기 육각형 격자 패턴은 상기 측정 대상(20)의 상 공액면(image conjugate plane)에 균일하게 조사된 간섭성 조명광을 공간적으로 2차원 진폭 변조하여 얻을 수 있다. 상기 2차원 진폭 변조된 광 분포는 축소 광학계 (reduction optics, 127)를 통해 미세한 간섭성 구조조명 패턴으로 상기 현미경 광학계(130)를 통하여 상기 측정 대상(20)에 투사될 수 있다.
상기 광원부(110)는 간섭성이 우수한 레이저를 포함할 수 있다. 상기 광원부(110)는 레이저 광원(111), 전기광학 셔터(112), 광섬유(114), 그리고 집속 렌즈(116)를 포함할 수 있다. 상기 레이저 광원(111)은 간섭성이 우수한 적외선 대역 또는 가시광선 대역의 광을 출력할 수 있다. 상기 레이저 광원의 출력 광은 상기 전기광학 셔터(112)를 통하여 스위칭될 수 있다. 상기 전지광학 셔터(112)를 통과한 광은 광섬유(114)를 통하여 가이드될 수 있다. 상기 광섬유(114)를 통하여 가이드된 광은 집속 렌즈(collimating lens, 116)에 제공될 수 있다. 상기 집속 렌즈(116)는 평행광을 출력할 수 있다. 상기 평행광은 상기 구조조명 패턴 형성부(120)에 제공될 수 있다.
상기 구조조명 패턴 형성부(120)는 편광판(122), 빔스플릿터(beam splitter,124), 및 공간 광 변조부(spatial light modulator, 126)를 포함할 수 있다. 상기 평행광은 상기 편광판(122)에 제공될 수 있다. 이에 따라, 상기 편광판(122)은 선형 편광된 평행광을 출력할 수 있다. 상기 선형 편광된 평행광은 빔스플릿터(124)에 제공될 수 있다. 상기 빔스플릿터(124)는 50:50 광 분할기일 수 있다. 상기 빔스플릿터(124)는 입사하는 선형 편광된 평행광을 투과시킬 수 있다. 또한, 투과한 광은 상기 공간 광 변조부(126)에 의하여 공간 변조되고 반사될 수 있다. 상기 공간 변조된 광은 다시 상기 빔스플릿터(124)에 입사되어 90도 꺾여 진행할 수 있다. 상기 공간 광 변조부(126)는 2차원회절격자(2D diffraction grating), 액정 디스플레이(liquid crystal display; LCD)형 공간 광 변조기(spatial light modulator: SLM), 또는 디지털 거울 소자 (digital mirror device: DMD)일 수 있다. 상기 공간 광 변조부(124)는 시간에 따라 동일한 패턴을 가지나 서로 다른 공간 위상을 가지는 육각형 격자 패턴을 형성할 수 있다.
육각형 격자 진폭 구조조명은 육각형 격자 패턴을 제공할 수 있다. 육각형 격자 진폭 구조조명은 간섭성 구조조명 패턴의 형태에서는 동일하나, 서로 다른 공간 위상을 가질 수 있다. 구체적으로, 육각형 격자 패턴은 2차원 공간의 전체 방위각을 3 등분하여, {0o, 120o, 240o} 방향으로 격자 벡터를 갖는 3 개의 선형 격자 진폭 구조조명을 진폭 합산(amplitude summation)하여 얻어질 수 있다. 육각형 격자 패턴은 순차적으로 형성될 수 있다. 육각형 격자 패턴의 개수는 19개일 수 있다. 상기 공간 광 변조부는 현미경 시스템의 측정 대상의 상 공액면(image conjugate plane)에 배치될 수 있다.
축소광학계(127)는 상기 육각형 격자 패턴 또는 육각형 격자 진폭 구조조명을 축소시키어 현미경 광학계(130)에 제공할 수 있다.
빔스플릿터(128)는 상기 축소 광학계(127)가 제공한 축소된 육각형 격자 패턴을 현미경 광학계(130)에 제공할 수 있다. 상기 빔스플릿터(128)는 상기 현미경 광학계의 하부 대물렌즈부의 하부에 배치될 수 있다. 상기 빔스플릿터(128)는 상기 육각형 격자 진폭 구조조명을 반사시켜 상기 상기 현미경 광학계(130)에 제공하고, 상기 측정 대상(20)으로부터 오는 반사 광을 투과시킬 수 있다.
현미경 광학계(130)는 시료 지지대(136), 하부 대물 렌즈부(132), 그리고 상부 대물 렌즈부(134)를 포함할 수 있다. 상기 축소된 육각형 격자 패턴은 상기 하부 매물 렌즈부(132)를 통하여 상기 측정 대상(20)에 제공될 수 있다. 상기 측정 대상(20)은 시료 지지대(136)에 고정될 수 있다. 상기 하부 대물 렌즈부(132) 및 상부 대물 렌즈부(134)는 상기 시료 지지대(136)를 중심으로 서로 마주보도록 배치될 수 있다. 상기 시료 지지대(136)는 이동부(138)를 통하여 이동할 수 있다.
반사 측정부(140)는 상기 측정 대상(20)에서 반사된 광을 수집할 수 있다. 구체적으로, 상기 측정 대상(20)에서 반사된 광은 상기 하부 대물 렌즈부(132)를 통하여 상기 반사 측정부(140)에 제공될 수 있다. 상기 반사 측정부(140)는 경통 렌즈(142)와 2차원 광 검출기(144)를 포함할 수 있다. 상기 측정 대상에서 반사된 광은 상기 2차원 광 검출기(144)를 통하여 원시 이미지로 변환될 수 있다. 원시 이미지는 육각형 격자 진폭 구조조명의 3 방위 진폭 위상 조합에 따라 19개 일 수 있다.
투과 측정부(150)는 경통 렌즈(152)와 2차원 광 검출기(154)를 포함할 수 있다. 상기 측정 대상(20)을 투과한 광은 상기 2차원 광 검출기(154)를 통하여 원시 이미지로 변환될 수 있다. 상기 원시 이미지의 개수는 간섭성 구조조명 패턴의 개수와 동일할 수 있다.
처리부(160)는 상기 원시 이미지들을 처리하여 초고분해능의 공간 이미지를 생성할 수 있다.
본 발명에 이용되는 광원은 간섭성이 우수한 레이저인 것이 바람직하다. 현미경 광학계는 2차원 주기 패턴 (2D periodic pattern)으로 간섭성 구조조명된 상기 측정 대상의 확대상(magnified image)을 형성할 수 있다. 상기 확대상은 2차원 광 검출기(2D CCD camera)로 측정되고, 상기 2차원 광 검출기는 원시 이미지를 출력할 수 있다. 상기 현미경 광학계는 투과단 또는 반사단 쪽에 설치된 대물렌즈와 경통렌즈를 포함할 수 있다.
본 발명의 일 실시예에 따르면, 2차원 간섭성 구조조명 이미징에 적합한 광 진폭변조(amplitude modulation) 구조조명 패턴이 제안된다. 또한, 상기 간섭성 구조조명에 의한 측정 대상의 강도 이미지(원시 이미지)를 공간주파수 영역에서 분석하는 수학적으로 모델링이 제안된다. 상기 강도 이미지(원시 이미지)는 소정의 이미지 처리 알고리즘을 통하여 공간 이미지로 재구성된다.
보다 자세하게, 상기 이미지 처리 알고리즘은, 물체의 강도 이미지(원시 이미지)를 상기 간섭성 구조조명에 의해 유도된 주파수 하향변환 스펙트럼(down-converted spectrum) 간의 상관혼합 (correlation mixing)으로 모델링한다. 또한, 상기 이미지 처리 알고리즘은 물체의 초고분해능 공간주파수 스펙트럼 성분을 분리하는 역혼합 행렬(unmixing matrix) 도출과 물체의 확장된 공간주파수 스펙트럼(extended image spectrum)을 복원하는 수학적 체계에 근간하여 구현된다.
이하, 본 발명의 2차원 간섭성 구조조명 이미징의 구현에 필요한 광 진폭 변조(amplitude modulation) 구조조명 패턴 설계(A), 상기 광 진폭 변조 구조조명에 의한 물체 강도 이미지(원시 이미지) 형성에 관한 수학적 모델링(B), 원시 이미지 데이터(raw image data) 취득법(C), 그리고 이미지 공간주파수 영역에서 분석 및 처리를 통한 초고분해능 이미지 재구성 알고리즘을 제공하는 수학적 이론 체계(D)가 설명된다.
도 7을 참조하면, 본 발명의 일 실시예에 따른 간섭성 구조조명 이미지 처리 방법은 육각형 격자 패턴에 의하여 형성된 육각형 격자 진폭 구조조명을 측정 대상에 차례로 조사하는 단계(S110), 상기 육각형 격자 진폭 구조조명들에 의하여 주파수 하향 변환된 원시 이미지들을 차례로 측정하는 단계(S120), 및 상기 원시 이미지들을 처리하여 재구성된 물체 이미지를 산출하는 단계(S130)를 포함한다.
상기 원시 이미지들을 처리하여 재구성된 물체 이미지를 산출하는 단계(S130)는 상기 육각형 격자 진폭 구조조명의 3 방위 진폭 조합에 의하여 정의되는 스펙트럼 혼합 행렬의 역행렬을 구하는 단계(S131), 상기 원시 이미지들을 푸리에 변환한 원시 이미지 스펙트럼을 산출하는 단계(S132), 상기 원시 이미지 스펙트럼으로 구성된 스펙트럼 열 벡터와 상기 스펙트럼 혼합 행열의 역행렬을 이용하여 상기 원시 이미지 스펙트럼을 구성하는 측정 스펙트럼 구성항을 추출하는 단계(S133), 상기 측정 스펙트럼 구성항을 공간 주파수 되돌림하여 가상 스펙트럼 구성항을 산출하는 단계(S134), 상기 가상 스펙트럼 구성항을 선형 결합하여 가상 이미지 스펙트럼을 형성하는 단계(S135), 및 상기 가상 이미지 공간 주파수 스펙트럼을 푸리에 역변환하여 물체 이미지 또는 가상 이미지를 산출하는 단계(S136)를 포함한다.
상기 육각형 격자 패턴은 3 방위각 복소 진폭 간섭성 구조조명 패턴으로, 2차원 공간의 전체 방위각을 3 등분한 {0o, 120o, 240o} 방향으로 격자 벡터를 갖는 3 개의 선형 격자 패턴을 진폭 합산(amplitude summation)하여 얻어진다. 상기 육각형 격자 패턴은 3 방위 위상 조합에 의하여 형성되고, 상기 3 방위 위상 조합은 3축 위상 공간의 19 개 위상 벡터 점으로 표시될 수 있다.
본 발명의 일 실시예에 따르면, 간섭성 구조조명 이미지 처리 방법은 2차원 초고분해능 간섭성 구조조명 이미징 기술을 제공한다. 준등방성(quasi-isotropic) 공간분해능 향상을 구현하기 위하여, 본 발명에 따른 초고분해능 이미지 재구성결과로 얻어지는 확장 간섭성 전달함수(extended CTF)는 도 5에서 실선으로 나타낸 형태의 진폭 차단 주파수 경계를 가질 수 있다.
본 발명의 일 실시예에 따르면, 종래 간섭성 구조조명 이미징 기법에 비하여 진폭 차단 주파수 확장의 등방성이 향상될 수 있다. 이를 위하여, 육각형 격자 패턴(hexagonal lattice pattern) 형태의 2차원 진폭 구조조명이 제안된다. 상기 육각형 격자 진폭 구조조명은 2차원 공간의 전체 방위각을 3등분한 {0o, 120o, 240o}방향으로 격자 벡터(grating vector)를 갖는 3개의 선형격자 진폭 구조조명을 간섭성 중첩(amplitude summation)하여 얻을 수 있다. 상기 선형격자 진폭 구조조명의 주기는 진폭 차단 주파수 의 역수인 보다 큰 임의의 값이 허용되지만 최대한 근접한 작은 값으로 설정하는 것이 공간분해능 향상 효과 측면에서 바람직하다. 제1 격자 벡터(V1)는 직각 좌표계를 기준으로, 0도 방향의 단위 벡터이고, 제2 격자 벡터(V2)는 시계 반대 방향으로 120도 방향의 단위 벡터이고, 제3 격자 벡터(V3)는 시계 반대 방향으로 240도 방향의 단위 벡터일 수 있다. 육각형 격자 패턴의 단위 셀은 육각형의 중심에서 최대 강도를 가질 수 있다. 육각형 격자 패턴의 단위 셀은 3 방위 진폭 위상 조합에 의하여 변경될 수 있다. 구체적으로, 상기 육각형 격자 패턴의 단위 셀에 있어서, 3 방위 진폭 위상 조합에 의하여 최대 강도의 위치는 0번 내지 18번 지점으로 이동될 수 있다. 따라서, 19개의 육각형 격자 패턴이 형성될 수 있다.
[간섭성 이미징 광학 배경 이론]
간섭성 광원에 의해 조명된 진폭 물체의 강도 이미지 형성을 기술하기 위하여 관련 물리량을 다음과 같이 정의한다.
공간 변조된 복소 진폭을 갖는 조명 필드 패턴(또는 구조조명 광파) 에 의해 조명되는 진폭 물체(또는 시료 분포함수, 측정 대상에 대응) 에 의하여 형성되는 측정 이미지 강도 분포(원시 이미지에 대응) 는
으로 주어진다. 이 때, 은 2차원 공간 변위 벡터이고, 는 2차원 공간주파수 벡터이다. 수식에 사용된 표기인 와 는 각각 함수의 컨볼루션(convolution) 연산과 자기상관(auto-correlation) 연산을 의미한다. 대문자화된 함수 또는 ~ 표기가 된 함수는 푸리에 변환(Fourier transform )된 것을 의미한다. 특별히, 간섭성 점퍼짐 함수 (CSF)의 Fourier transform인 간섭성 전달함수 (CTF)는 로 표시된다. 비간섭성 점퍼짐 함수 (PSF)는 CSF와 의 관계가 성립하기 때문에, PSF의 푸리에 변환은 로 주어지며, 이를 광전달함수 (OTF)라 명하는 것이 관례이다. 공간 주파수 스펙트럼과 스펙트럼은 동일한 의미로 사용한다.
[목표로 하는 확장 간섭성 전달 함수에 의한 가상 초고분해능 이미지 수식화]
본 발명의 2차원 간섭성 구조조명법을 이용하여 도 5에 도시한 확장 간섭성 전달함수를 가지는 가상(virtual)의 초고분해능 광학계를 구현하고자 할 때, 목표로 하는 가상 CTF는
으로 정의할 수 있다. 이 때, 는 공간주파수 벡터이며, 가상 CTF는 이미징 광학계 본래의 간섭성 전달 함수 가 6개의 공간주파수 벡터 만큼 각각 변위되어 중첩된 형태이다. 간섭성 전달 함수의 공간주파수 변위 벡터 는
는 공간주파수 변위의 방위각(단위: 라디안)을 나타낸다.
으로 표현된다. 이 때, 로 표시한 연산은 두 함수의 교차상관(cross-correlation)이다. 간섭성 가상 이미지의 공간주파수 스펙트럼 은 모두 19개의 스펙트럼 구성항 으로 구성된다. 는 전체 확장 CTF를 구성하는 6개의 개별 요소 CTF에 의한 독립적 이미지 강도 스펙트럼의 합이다. 한편, 나머지 모든 항은 서로 다른 요소 CTF에 의한 이미지 진폭 스펙트럼 간의 교차상관 양상을 가진다.
[2차원 간섭성 구조조명에 의한 측정 이미지 수식화]
도 6의 (b) 참조하면, 상기 확장 CTF에 의한 간섭성 가상 이미지에 포함된 공간주파수 스펙트럼 구성항과 유사 내지 동등한 초고분해능 이미지 정보를 추출하기 위한 간섭성 구조조명은, 물리적 직관에 의해 육각형 격자(hexagonal lattice pattern) 형태의 2차원 진폭 변조 광파로 제공될 수 있다.
상기 정육각형 격자 형태의 진폭 변조 조명 필드 패턴 는 가상 확장 CTF를 구성하는 6개 요소 CTF의 공간주파수 변위 벡터 와 같은 공간주파수를 갖는 6개 광파의 중첩으로 이루어지며 다음과 같이 표현된다.
이므로, 수학식(7)을 이용하면 다음과 같이 전개된다.
[확장 CTF 초고분해능 가상 이미지와 간섭성 구조조명 이미지 간의 관계성]
확장 CTF를 갖는 가상의 초고분해능 이미징 광학계에서 얻을 수 있는 가상 이미지(수학식(6)으로 표현, 가상 이미지는 물체 이미지와 동일함)를 육각형 격자 진폭 구조조명에 의한 원시 이미지(수학식(9)로 표현)로부터 복원하는 본 발명의 이미지 처리 알고리듬에 대하여 설명한다.
상기 확장 CTF에 의한 가상 이미지 스펙트럼 과 상기 육각형 격자 진폭 구조조명에 의한 원시 이미지 스펙트럼 은 서로 같은 개수의 함수항으로 구성된다. 상기 두 이미지 스펙트럼을 각각 구성하는 동일 차수 의 함수항 과 함수항 는 다음과 같이 일대일 대응되는 해석식이 성립한다.
도 7을 참조하면, 적절한 측정을 통하여 상기 측정 스펙트럼 구성항을 독립적으로 구하면, 수학식(11)의 대응관계를 이용하여 모든 차수의 상기 가상 스펙트럼 구성항을 취득함으로써 목표로 하는 초고분해능의 가상 이미지 스펙트럼 을 재구성할 수 있다.
[구조조명 원시 이미지로부터의 측정 스펙트럼 구성항 획득 방법]
수학식(9) 참조하면, 본 발명의 구조조명 물체에 대해 측정한 원시 이미지를 푸리에 변환하여 얻은 원시 이미지 스펙트럼 은 구조조명의 3방위 진폭 위상 조합이 주어질 때, 다음과 같이 측정 스펙트럼 구성항 들의 선형적 합으로 기술된다.
따라서, 미지의 19개 측정 스펙트럼 구성항 를 구하기 위해서는, 적절히 선택된 19개의 3방위 진폭 위상조합 에 대하여 원시 이미지 스펙트럼을 구한다. 상기 원시 이미지 스펙트럼을 기술하는 연립 선형대수 방정식
으로부터 상기 측정 스펙트럼 구성항 의 해(solution)을 찾아야 한다. 이 때, 은 상기 3방위 진폭 위상 조합 의 설정에 따라 결정되는 19 x 19 스펙트럼 혼합 행렬(spectral mixing matrix)이다.
의 조건을 만족시키도록 해야 한다.
일례로써, 본 발명에서 제공하는 상기 3방위 진폭 위상 조합은,
수학식(14), (15), 및 (16)을 참조하면, 본 발명의 3방위 진폭 위상 조합은 도 6의(c)에 도시한 3축 위상 공간의 19개 위상 벡터 점으로 표시할 수 있다. 또한, 상기 3방위 진폭 위상 조합의 변경에 따른 육각형 격자 진폭 구조조명 패턴의 변화는 도 6의 (b)에 도시한 육각형 격자 구조조명 패턴이 상기 19개 위상 벡터에 대응하는 변위만큼 이동하는 효과로 나타난다.
상기 스펙트럼 혼합 행렬(spectral mixing matrix) 의 역행렬(matrix inverse)을 스펙트럼 역혼합 행렬(spectral unmixing matrix)이라 한다. 상기 3방위 진폭 위상 조합 에 대하여 구조조명에 대한 측정된 원시 이미지 스펙트럼으로 형성한 열 벡터(column vector)로부터
[초고분해능 이미지 스펙트럼 복원 및 실공간 이미지 재구성]
도 8은 간섭성 구조조명 이미지로부터 스펙트럼 역혼합 계산을 통해 분리한 공간 주파수 하향 변환 측정 스펙트럼 구성항을 나타내는 도면이다.
도 9는 공간 주파수 하향 변환 측정 스펙트럼 구성항에 대한 공간 주파수 영역 재구성 복원 위치를 나타내는 도면이다.
도 10은 등방성 2차원 간섭성 구조조명 이미징에 의한 공간 주파수 스펙트럼 차단 주파수 확장 효과를 나타내는 도면이다.
도 7을 참조하면, 상기 측정 스펙트럼 구성항은 각각 일대일 대응되는 가상 스펙트럼 구성항 이 공간주파수 원점 방향으로 특정 공간주파수 벡터만큼 이동된 함수이다. 상기 측정 스펙트럼 구성항은 확장 CTF에 의한 초고분해능의 가상 이미지를 구성하는 요소 스펙트럼이다.
따라서, 목표로 하는 초고분해능의 가상 이미지의 스펙트럼의 재구성은, 도 8에 도시한 바와 같은 원시 이미지 스펙트럼의 각 구성항을 상응하는 가상 스펙트럼 구성항 본래의 공간주파수 위치로 되돌림(spectral back-shift)으로써 이루어질 수 있다.
측정된 원시 이미지 스펙트럼 구성항 중, 항을 제외한 모든 항을 도 9에 도시한 바와 같이 대응되는 가상 스펙트럼 구성항 의 공간주파수 영역으로 재배치시켜서(spectral back-shift) 가상 스펙트럼 구성항 화시킨다.
도 10을 참조하면, 이후 가상 스펙트럼 구성항 은 전부 선형적으로 합쳐진다.이에 따라, 확장 CTF에 의한 가상 이미지 스펙트럼과 동등한 초고분해능의 이미지 스펙트럼으로 재구성된다.
또한, 상기 초고분해능의 가상 이미지 스펙트럼을 역 푸리에 변환(inverse Fourier transform)하면 실공간(real space)에서의 물체 공간 이미지(가상 이미지)로 복원된다.
[등방성 2차원 간섭성 구조조명 이미징에 의한 공간분해능 향상 특성]
위에서 설명한 등방성 2차원 간섭성 구조조명에 의한 측정과 공간주파수 영역 이미지 스펙트럼 재구성 계산을 적용해 복원시킨 진폭 물체(측정 대상)의 초고분해능 이미지 스펙트럼은 본 발명이 목표로 하는 이미징 광학계의 간섭성 전달 함수(CTF) 확장 효과를 포함한다.
도 11은 등방성 2차원 간섭성 구조조명 이미징에 실공간 이미지 공간 분해능 향상 효과를 나타내는 도면이다.
도 12는 종래의 직교 2-방위 구조조명 이미징과 본 발명의 등방성 3-방위 구조조명 이미징의 공간 분해능 특성 비교하는 도면이다.
도 2의 (0)에 도시한 진폭 물체 타겟을 대상으로 한 일반 간섭성 이미징과 본 발명의 등방성 2차원 간섭성 구조조명 이미징을 컴퓨터 시뮬레이션한 결과를 도 10과 도 11을 참조하여 설명한다.
도10의 (b) 참조하면, 본 발명의 등방성 2차원 간섭성 구조조명 이미징에 의한 진폭 물체의 이미지 스펙트럼이 표시된다. 도 10의 (a) 를 참조하면, 일반 간섭성 이미징에 의한 이미지 스펙트럼이 표시된다.
본 발명의 등방성 2차원 간섭성 구조조명 이미징에 의한 진폭 물체의 이미지 스펙트럼은 일반 간섭성 이미징에 의한 이미지 스펙트럼과 비교할 때, 진폭 물체의 공간주파수 스펙트럼 정보를 취득할 수 있는 영역이 확장된다. 즉, 2차원 진폭 물체의 공간주파수 스펙트럼을 취득할 수 있는 차단 주파수 한계가 본 발명에서 적용한 진폭 구조조명의 3개 방위각에 대하여 최대 2배 확장됨을 알 수 있다.
결과적으로, 본 발명의 등방성 2차원 간섭성 구조조명 이미징에 의해 차단 주파수 한계가 확장된 진폭 물체 공간주파수 스펙트럼을 역 푸리에 변환(inverse Fourier transform)하면 실공간(real space)에서 광 회절한계를 넘어서는 간섭성 이미지를 얻을 수 있다.
즉, 도 11의 (a)에 도시한 일반 간섭성 이미지에 비하여 2차원 공간분해능이 2배 향상된 초고분해능 간섭성 이미지(도 11(b) 참조)를 취득할 수 있다.
또한, 본 발명에 의한 등방성 2차원 간섭성 구조조명 이미징은 3-방위각 진폭 구조조명을 사용함으로써, 직교 2-방위각 진폭 구조조명을 사용한 종래의 간섭성 구조조명 이미징에 비해 차단 주파수가 확장되는 효과와 방향성이 공간주파수 영역에서 보다 균일해지는 장점이 있다.
도 11 참조하면, 보다 상세하게, 종래의 직교 2-방위각 간섭성 구조조명 이미징은 x축과 y축 방향으로 최대의 2배 의 공간분해능 향상 효과가 발생한다.
도 11(a)를 참조하면, x축 및 y축을 기준으로 한 대각선 방위에 대해서는 공간분해능 향상 정도가 약 30 퍼센트(%) 정도 저하된다.
반면, 본 발명이 제공하는 3-방위각 간섭성 구조조명 이미징은 최대 2배의 공간분해능 향상 효과가 전 방위에 걸쳐서 보다 등방적(isotropic)으로 나타나며, 공간분해능 향상 효과가 변화폭이 약 13 퍼센트(% )이하로 유지된다.
110: 광원부
120: 구조조명 패턴 형성부
130: 현미경 광학계
140,150: 측정부
160: 처리부
120: 구조조명 패턴 형성부
130: 현미경 광학계
140,150: 측정부
160: 처리부
Claims (8)
- 육각형 격자 패턴에 의하여 형성된 육각형 격자 진폭 구조조명을 측정 대상에 차례로 조사하는 단계;
상기 육각형 격자 진폭 구조조명들에 의하여 주파수 하향 변환된 원시 이미지들을 차례로 측정하는 단계; 및
상기 원시 이미지들을 처리하여 재구성된 물체 이미지를 산출하는 단계를 포함하는 것을 특징으로 하는 간섭성 구조조명 이미징 방법. - 제1 항에 있어서,
상기 원시 이미지들을 처리하여 재구성된 물체 이미지를 산출하는 단계는:
상기 육각형 격자 진폭 구조조명의 3 방위 진폭 조합에 의하여 정의되는 스펙트럼 혼합 행렬의 역행렬을 구하는 단계;
상기 원시 이미지들을 푸리에 변환한 원시 이미지 스펙트럼을 산출하는 단계;
상기 원시 이미지 스펙트럼으로 구성된 스펙트럼 열 벡터와 상기 스펙트럼 혼합 행열의 역행렬을 이용하여 상기 원시 이미지 스펙트럼을 구성하는 측정 스펙트럼 구성항을 추출하는 단계;
상기 측정 스펙트럼 구성항을 공간 주파수 되돌림하여 가상 스펙트럼 구성항을 산출하는 단계;
상기 가상 스펙트럼 구성항을 선형 결합하여 가상 이미지 스펙트럼을 형성하는 단계; 및
상기 가상 이미지 공간 주파수 스펙트럼을 푸리에 역변환하여 물체 이미지 또는 가상 이미지를 산출하는 단계;를 포함하는 것을 특징으로 하는 간섭성 구조조명 이미징 방법. - 제1 항에 있어서,
상기 육각형 격자 패턴은 3 방위각 복소 진폭 간섭성 구조조명 패턴으로, 2차원 공간의 전체 방위각을 3 등분한 {0o, 120o, 240o}방향으로 격자 벡터를 갖는 3 개의 선형 격자 패턴을 진폭 합산(amplitude summation)하여 얻어지는 것을 특징으로 하는 간섭성 구조조명 이미징 방법. - 제3 항에 있어서,
상기 육각형 격자 패턴은 3 방위 위상 조합에 의하여 형성되고, 상기 3 방위 위상 조합은 3축 위상 공간의 19 개 위상 벡터 점으로 표시되는 것을 특징으로 간섭성 구조조명 이미징 방법. - 간섭성 광을 출력하는 광원부;
상기 광원부로부터 상기 간섭성 광을 제공받아 육각형 격자 패턴을 가지는 육각형 격자 진폭 구조조명를 생성하는 구조조명 패턴 형성부;
상기 육각형 격자 진폭 구조조명을 제공받아 측정 대상에 제공하는 현미경 광학계; 및
상기 육각형 격자 진폭 구조조명을 제공받은 측정 대상이 반사 또는 투과시키는 광을 측정하는 측정부;를 포함하는 것을 특징으로 하는 간섭성 구조조명 현미경 시스템. - 제5 항에 있어서,
상기 측정부는 상기 육각형 격자 진폭 구조조명에 의한 원시 이미지를 측정하고,
상기 원시 이미지를 처리하여 물체 이미지를 재구성하는 처리부를 더 포함하는 것을 특징으로 하는 간섭성 구조조명 현미경 시스템. - 제5 항에 있어서,
상기 구조조명 패턴 형성부는 상기 측정 대상의 상 공액면(image conjugate plane)에 배치되는 것을 특징으로 하는 간섭성 구조조명 현미경 시스템. - 제5 항에 있어서,
상기 측정부는:
투과단 또는 반사단 방향에 설치된 경통 렌즈; 및
상기 현미경 광학계 및 상기 경통 렌즈를 통하여 형성된 확대 상(magnified image)을 2차원적으로 측정하는 광 검출기를 포함하는 것을 특징으로 하는 간섭성 구조조명 현미경 시스템.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR20130135555A KR101479249B1 (ko) | 2013-11-08 | 2013-11-08 | 간섭성 구조조명 이미징 방법 및 간섭성 구조조명 현미경 시스템 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR20130135555A KR101479249B1 (ko) | 2013-11-08 | 2013-11-08 | 간섭성 구조조명 이미징 방법 및 간섭성 구조조명 현미경 시스템 |
Publications (1)
Publication Number | Publication Date |
---|---|
KR101479249B1 true KR101479249B1 (ko) | 2015-01-05 |
Family
ID=52587895
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR20130135555A KR101479249B1 (ko) | 2013-11-08 | 2013-11-08 | 간섭성 구조조명 이미징 방법 및 간섭성 구조조명 현미경 시스템 |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101479249B1 (ko) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019078578A1 (ko) * | 2017-10-18 | 2019-04-25 | 고려대학교 산학협력단 | 구조 조명과 위상 검출을 이용한 단분자 중심위치 측정 장치 및 방법 |
KR20210003594A (ko) * | 2019-07-02 | 2021-01-12 | 한국과학기술연구원 | 체커보드 패턴 기반 구조 조명 현미경 시스템 |
WO2021020603A1 (ko) * | 2019-07-29 | 2021-02-04 | 한국표준과학연구원 | 복합 패턴을 이용한 초고속 편향 측정법을 이용한 자유곡면의 3차원 형상측정시스템 및 측정방법 |
KR20220167513A (ko) | 2021-06-14 | 2022-12-21 | 연세대학교 산학협력단 | 구조조명 현미경 시스템 및 그 제어방법 |
CN116380408A (zh) * | 2023-04-10 | 2023-07-04 | 南京航空航天大学 | 基于结构光和光场成像的三维超分辨流场测量方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006030178A (ja) | 2004-07-16 | 2006-02-02 | Carl Zeiss Jena Gmbh | 高度な分解能を持つ顕微鏡 |
KR100709922B1 (ko) | 1999-06-24 | 2007-04-24 | 에스브이지 리도그래피 시스템즈, 아이엔씨. | 광학 시스템의 특성화를 위한 방법 및 장치 |
KR20070104615A (ko) * | 2005-01-20 | 2007-10-26 | 지고 코포레이션 | 객체 표면의 특성을 결정하기 위한 간섭계 |
JP2011028208A (ja) | 2009-06-22 | 2011-02-10 | Yokogawa Electric Corp | 顕微鏡装置 |
-
2013
- 2013-11-08 KR KR20130135555A patent/KR101479249B1/ko active IP Right Grant
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100709922B1 (ko) | 1999-06-24 | 2007-04-24 | 에스브이지 리도그래피 시스템즈, 아이엔씨. | 광학 시스템의 특성화를 위한 방법 및 장치 |
JP2006030178A (ja) | 2004-07-16 | 2006-02-02 | Carl Zeiss Jena Gmbh | 高度な分解能を持つ顕微鏡 |
KR20070104615A (ko) * | 2005-01-20 | 2007-10-26 | 지고 코포레이션 | 객체 표면의 특성을 결정하기 위한 간섭계 |
JP2011028208A (ja) | 2009-06-22 | 2011-02-10 | Yokogawa Electric Corp | 顕微鏡装置 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019078578A1 (ko) * | 2017-10-18 | 2019-04-25 | 고려대학교 산학협력단 | 구조 조명과 위상 검출을 이용한 단분자 중심위치 측정 장치 및 방법 |
KR20210003594A (ko) * | 2019-07-02 | 2021-01-12 | 한국과학기술연구원 | 체커보드 패턴 기반 구조 조명 현미경 시스템 |
KR102260563B1 (ko) * | 2019-07-02 | 2021-06-07 | 한국과학기술연구원 | 체커보드 패턴 기반 구조 조명 현미경 시스템 |
WO2021020603A1 (ko) * | 2019-07-29 | 2021-02-04 | 한국표준과학연구원 | 복합 패턴을 이용한 초고속 편향 측정법을 이용한 자유곡면의 3차원 형상측정시스템 및 측정방법 |
KR20220167513A (ko) | 2021-06-14 | 2022-12-21 | 연세대학교 산학협력단 | 구조조명 현미경 시스템 및 그 제어방법 |
CN116380408A (zh) * | 2023-04-10 | 2023-07-04 | 南京航空航天大学 | 基于结构光和光场成像的三维超分辨流场测量方法及系统 |
CN116380408B (zh) * | 2023-04-10 | 2024-06-07 | 南京航空航天大学 | 基于结构光和光场成像的三维超分辨流场测量方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10606055B2 (en) | Aperture scanning Fourier ptychographic imaging | |
US8184279B2 (en) | Fourier domain sensing | |
Tahara et al. | Roadmap of incoherent digital holography | |
KR101479249B1 (ko) | 간섭성 구조조명 이미징 방법 및 간섭성 구조조명 현미경 시스템 | |
US6909105B1 (en) | Method and device for representing an object | |
US20120257197A1 (en) | Fourier domain sensing | |
CN101526400B (zh) | 阿达玛变换干涉光谱成像方法及设备 | |
US20110032586A1 (en) | Light microscope with novel digital method to achieve super-resolution | |
JP2014137558A (ja) | 撮像装置、撮像システム、および画像処理方法 | |
Wang et al. | Moving target tracking and imaging through scattering media via speckle-difference-combined bispectrum analysis | |
Cossairt et al. | Digital refocusing with incoherent holography | |
Shan et al. | Quasi-common-path off-axis interferometric quantitative phase microscopy based on amplitude-division | |
Nürnberg et al. | A simulation framework for the design and evaluation of computational cameras | |
O’holleran et al. | Methodology for imaging the 3D structure of singularities in scalar and vector optical fields | |
KR101501512B1 (ko) | 2차원 위상정합 비선형광학 구조조명 현미경 이미징 장치 및 그 방법 | |
Dumas et al. | A compressed sensing approach for resolution improvement in fiber-bundle based endomicroscopy | |
Zhang et al. | Bi-channel compressive hyperspectral imager based on polarization conversion metasurface | |
Hallada et al. | Fresnel zone plate light field spectral imaging simulation | |
Sola-Pikabea et al. | Optimal design of incoherent tunable-frequency structured illumination microscope scheme | |
Gureyev et al. | Spatial resolution and signal-to-noise ratio in x-ray imaging | |
Huang et al. | Non-interferometric accurate phase imaging via linear-convergence iterative optimization | |
Shabani et al. | Comparison of two structured illumination techniques based on different 3D illumination patterns | |
Classen | Enhancing x-ray diffractive imaging and optical microscopy by use of intensity correlations | |
Wang et al. | Super-resolution imaging lidar based on Fourier ptychography | |
Stokoe et al. | Theory and applications of structured light single pixel imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
E701 | Decision to grant or registration of patent right | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20181126 Year of fee payment: 5 |
|
FPAY | Annual fee payment |
Payment date: 20190926 Year of fee payment: 6 |