KR100620807B1 - 분해능을 향상시킨 파면왜곡 측정 방법 - Google Patents

분해능을 향상시킨 파면왜곡 측정 방법 Download PDF

Info

Publication number
KR100620807B1
KR100620807B1 KR1020040105272A KR20040105272A KR100620807B1 KR 100620807 B1 KR100620807 B1 KR 100620807B1 KR 1020040105272 A KR1020040105272 A KR 1020040105272A KR 20040105272 A KR20040105272 A KR 20040105272A KR 100620807 B1 KR100620807 B1 KR 100620807B1
Authority
KR
South Korea
Prior art keywords
image
center
point
center point
measuring
Prior art date
Application number
KR1020040105272A
Other languages
English (en)
Other versions
KR20060066802A (ko
Inventor
백성훈
박승규
차병헌
Original Assignee
한국원자력연구소
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 한국원자력연구소 filed Critical 한국원자력연구소
Priority to KR1020040105272A priority Critical patent/KR100620807B1/ko
Publication of KR20060066802A publication Critical patent/KR20060066802A/ko
Application granted granted Critical
Publication of KR100620807B1 publication Critical patent/KR100620807B1/ko

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
    • G01J9/04Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength by beating two waves of a same source but of different frequency and measuring the phase shift of the lower frequency obtained
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00Photometry, e.g. photographic exposure meter
    • G01J1/42Photometry, e.g. photographic exposure meter using electric radiation detectors
    • G01J1/44Electric circuits
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M11/00Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
    • G01M11/02Testing optical properties
    • G01M11/0242Testing optical properties by measuring geometrical properties or aberrations
    • G01M11/025Testing optical properties by measuring geometrical properties or aberrations by determining the shape of the object to be tested
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00Photometry, e.g. photographic exposure meter
    • G01J1/42Photometry, e.g. photographic exposure meter using electric radiation detectors
    • G01J1/44Electric circuits
    • G01J2001/4446Type of detector
    • G01J2001/448Array [CCD]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J9/00Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
    • G01J2009/002Wavefront phase distribution

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Geometry (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

본 발명은 분해능을 향상시킨 파면왜곡 측정 방법에 관한 것으로, 그 목적은 점 영상의 중심점 추출 정확도와 재현성을 향상시키기 위하여 기존의 무게중심법에 영상성분의 위치에 따른 가중치를 이용하여 신호 성분의 가중치는 크게 하고, 노이즈 성분의 가중치는 작게 하는 원리를 이용하는 방법을 제공하는 데 있다.
본 발명의 구성은 측정 대상이 되는 빛이 빔 확대기를 통해 측정기기에 적절한 크기로 입사시킨 후 배열렌즈에 의해 빛이 여러 개의 작은 구경으로 나뉘어 각각의 구경에 대한 초점을 맺도록 하는 단계와; 각각의 초점 영상이 렌즈를 통하여 카메라 등의 영상획득 장치로 입사되어 영상을 획득하는 단계와; 획득된 영상에 있는 점 영상의 중심점들의 중심 위치를 정확히 판단하기 위하여 각 점영상의 세기에 가중치를 주어 가중치를 준 무게 중심을 구함으로써 점 영상의 주변 노이즈를 제거한 위치를 구하는 단계와; 이렇게 구한 중심점의 위치를 기준이 되는 중심점의 위치와 비교함으로써 그 이동량으로부터 파면왜곡 정도를 측정하는 단계로 이루어진 방법을 특징으로 한다.
파면왜곡, 레이저빔, 하트만센서, 중심점추출, 파면측정장치

Description

분해능을 향상시킨 파면왜곡 측정 방법{High-resolution wavefront measurement method}
도 1은 파면측정장치 구성도
도 2는 하트만 센서에서의 배열 렌즈에 의한 점영상 형태도,
도 3은 문턱치 필터와 창이 없을 때의 중심점 위치 측정 결과도,
도 4는 문턱치 필터와 창이 없을 때의 중심점 위치 값들의 표준편차도,
도 5는 문턱치 필터는 없고 창의 이동이 있을 때의 중심점 위치 측정도,
도 6은 문턱치 필터와 창이 있을 때의 중심점 위치값들의 표준편차도이다.
<도면의 주요부분에 대한 부호의 설명>
(2, 3) : 렌즈 (4) : 배열렌즈
(5, 6) : 거울 (7) : 중계 렌즈
(11) : 레이저빔 (12) : 빔 확대기
(13) : 선형 이송기 (14) : CCD 카메라
본 발명은 분해능을 향상시킨 파면왜곡 측정 방법에 관한 것으로, 특히 레이저빔 파면보정과 광학계 수차 측정 등에 사용되는 파면측정장치인 하트만 센서와 이 센서의 분해능 향상 방법에 대한 것으로 자세하게는 하트만센서의 점영상의 중심점을 추출함에 있어서 가중치를 선별적으로 주어 주변 잡음의 효과를 줄이는 방법에 대한 것이다.
변형된 파면을 보정하는 적응광학 장치나 광학부품 검사 등에 사용되는 파면측정 장치는 여러 가지 종류가 있으나, 하트만 센서가 가장 일반적으로 사용된다.
상기 하트만 센서는 측정 대상이 되는 빔을 어레이 렌즈 등을 사용하여 여러 개의 부분 구경으로 나눈 후 각각의 구경에서의 집속 점의 변위를 측정하여 전체 파면을 측정하는 방식이다. 따라서 하트만 센서의 분해능은 각각의 집속 점의 중심점의 위치 측정 분해능과 직접 연관되어 있다. 또한 이러한 패턴의 중심점 측정은 하트만 센서 이외에도, 레이저 삼각측정 방법 등 레이저 계측장치에도 사용되고 있다.
상기 파면왜곡을 측정하는 하트만 센서는 배열렌즈 각각이 맺은 점영상의 위 치가 변하는 것을 측정하여 파면왜곡을 측정하게 된다. 따라서 측정 분해능을 결정짓는 가장 중요한 요소는 점영상의 중심점 추출 정확도이다. Roy Irwan, Richard G. Lane, "Analysis of optimal centroid estimation applied to Shack-Hartmann sensing," Applied Optics, Vol. 38, 6737 (1999)나 Daoshan Yang, Jihong Chen, Huichen Zhou, Shawn Buckley, "New algorithm to calculate the center of laser reflections," Proceedings of SPIE, Vol. 3306, 54 (1998) 등 몇몇 논문에서 이러한 점 영상의 중심을 추출하는 방법을 다루고 있으나, 대부분 무게중심법을 근간으로 하는 방법이며, 무게중심법은 아직까지도 파면측정 센서 신호처리의 가장 대표적인 방법이다.
또한 무게중심법은 하트만 센서 뿐 아니라 레이저 삼각측정법 등에서도 가장 많이 사용하는 방법이다.
그러나 점 영상에 많은 노이즈가 있거나, 파면 왜곡에 의한 점 영상의 왜곡 정도가 클 때에는 정확한 중심점을 추출하기 어려운 단점이 있다.
최근에는 일종의 필터링 방법(문턱치 노이즈 제거 필터, 고주파 노이즈 성분 제거 필터 등)을 사용하여 기존의 무게중심법의 단점을 보완하는 알고리즘을 많이 사용하기도 한다.
하지만 이러한 필터는 그 자체가 정보를 왜곡시킬 수 있으며, 추가의 계산 시간을 요하기 때문에 고속 적응광학 장치에는 약점이 된다
상기와 같은 문제점을 해결하기 위한 본 발명의 목적은 점 영상의 중심점 추출 정확도와 재현성을 향상시키기 위하여 기존의 무게중심법에 영상성분의 위치에 따른 가중치를 이용하여 신호 성분의 가중치는 크게 하고, 노이즈 성분의 가중치는 작게 하는 원리를 이용하는데, 계산 대상이 되는 점영상 신호의 그레이 레벨 값을 그 대로 사용하는 것이 아니라, 제곱, 세제곱 등의 지수 승을 사용함으로써, 큰 신호의 기여도는 살리고 작은 노이즈 성분의 기여도는 저하시키는 방법을 제공하는 데 있다.
상기한 바와 같은 목적을 달성하고 종래의 결점을 제거하기 위한 과제를 수행하는 본 발명의 구성은 다음과 같다.
영상의 중심점을 계산하는데 가장 일반적으로 사용되는 무게중심법에 의한 x-축 방향의 중심점은 다음 식과 같이 구한다.
Figure 112004058772816-pat00001
-----------------------------------(1)
상기 식(1)에서 Iij는 i-번째 행과 j-번째 열의 픽셀에서의 밝기를 나타내며, xi는 행과 열의 위치를 나타내는 좌표이고, x0y0 는 대략적인 중심의 위치이며, Wx x Wy 는 정밀한 중심 위치 추출의 계산을 위한 작은 창(Window)의 크기이다. 이 창의 크기가 크면 서브-픽셀 분해능의 정밀한 중심점을 계산하는데 시간이 많이 걸린다. 그러므로 대부분의 중심점 추출 방법으로는, 일단 큰 창을 대상으로 픽셀 분해능 정도의 대략의 중심점 위치(x0, y0) 계산하고, 그렇게 찾은 대략의 중심점 위치를 기준으로 작은 창(Wx x Wy)을 설정하여 그 영역 내에서 정밀한 중심점 계산을 하는 방식을 사용한다.
본 발명에서 새로 발명한 중심 위치 추출 알고리즘은 중심점 근처의 강한 빔세기(Iij) 값에는 가중치를 주어 중심점 추출에 기여를 많이 하도록 하고, 중심점에서 멀리 떨어진 세기가 약한 빔 세기는 중심점 계산시 기여를 작게 하도록 하는 것이다. 즉, 중심을 추출하는 과정에서 빔 세기 대신에 빔 세기의 제곱이나 세제곱 등의 값을 사용함으로써 이러한 효과를 얻을 수 있다. 이 경우 중심점에서 어느 정도 떨어진 위치에서의 각종 노이즈 성분을 작게 할 수 있으며, 전자적인 기저 노이즈도 효율적으로 줄일 수 있다. 이러한 개념을 더한 새로운 중심 위치 추출 알고리즘은 다음과 같이 나타낼 수 있다.
Figure 112004058772816-pat00002
------------------------------------(2)
상기 (2)식에서 α는 빔세기에 적용한 차수로, 가중치의 역할을 한다.
이하 본 발명의 실시예인 구성과 그 작용을 첨부도면에 연계시켜 상세히 설명하면 다음과 같다. 이하의 설명은 파장이 일반 CCD 카메라를 사용하는 파면측정 장치에 관한 예를 든 것이나, 다른 파면측정 장치에서도 동일한 원리의 적용이 가 능하다.
본 발명에서 발명한 파면측정 장치의 구성은 첨부한 도 1과 같다. 레이저빔(11)이 입사동으로 입사하면 렌즈(2)와 렌즈(3)으로 구성된 빔 확대기(12)에서 측정에 적절한 크기로 확대/축소 된다. 빔확대기를 거친 레이저빔은 배열렌즈(4)(여기서는 19x19의 배열렌즈를 사용함)에서 작은 여러 개의 구경으로 나뉘어지고 각각의 구경에 해당하는 빔은 초점면에 점영상을 맺는다. 이 점영상은 거울(5)와 거울(6)을 거쳐서 중계 렌즈(7)에 의해 CCD 카메라(14)의 센서 면에 영상을 맺게 된다. 정확한 영상을 얻기 위해서 선형 이송기(13)을 사용하게 된다.
이렇게 하여 획득된 점영상은 도 2와 같은 형태를 갖는다. 이러한 영상은 컴퓨터와 영상처리 장치에서 처리되어 파면을 측정하게 된다. 도 2에서 각각의 국부적인 영역에서의 점 영상은 어느 정도의 기저 노이즈를 가지며, 완전한 대칭 꼴이 아님을 알 수 있다.
본 발명 개발된 장치를 이용하여 앞의 식(1)의 방법과 식(2)의 방법으로 중심점을 추출하고, 그 결과를 비교하여 분석하였다. 각각의 경우에 점영상을 계산하는 창의 위치를 점 영상을 따라 이동시키는 경우와 이동시키지 않은 경우, 문턱치를 설정하여 어느 정도의 잡음을 제거한 경우와 전혀 잡음을 제거하지 않은 경우에 대하여 각각 계산하여 비교하였다.
도 3은 세기의 가중치인 지수 승을 변화시킬 때의 중심점의 위치를 계산한 결과이다. 이 때는 문턱치를 제거하는 필터가 없고, 중심점을 계산하는 창(여기서는 11x11 픽셀을 설정)의 이동이 없는 경우로 많은 오차가 발생할 수 있는 조건이다. 도면에서 알 수 있듯이 세기 가중치를 주어서 지수를 3으로 한 경우에는 점 들이 거의 선형적인 그래프를 보이며 직선을 보이지만, 가중치가 없는 경우에는 많은 오차를 가지며 그래프가 휘는 것을 알 수 있다.
도 4는 각각의 위치에서 측정한 50개 점 들의 측정 결과의 표준편차를 계산한 것이다. 세기의 지수 값에 따라서 표준편차도 영향을 받음을 도면에서 알 수 있다. 지수 값이 3일 때의 표준편차는 지수 값이 1일 때에 비해 반 이하인 것을 알 수 있다. 이 측정 결과에서 측정의 재현성에서도 지수의 값을 바꾸어 가중치를 준 것이 기존의 무게중심법에 비해 더 우수한 결과를 나타냄을 알 수 있다.
도 5는 창을 이동하면서 중심점을 추출할 때, 가중치를 준 경우와 가중치를 주지 않는 경우의 데이터(data)를 비교한 것이다. 기저 노이즈를 없애기 위한 문턱치 제거는 하지 않았다.
도 6은 문턱치를 30으로 설정하여 밝기의 30 이하는 전부 제거하고, 창의 위치도 변화시키면서 중심점을 추출했을 때의 각각의 표준편차를 계산한 것이다. 도 면에서 표준편차의 값이 1/5 이하로 줄어드는 것을 알 수 있다. 이것은 그만큼 측정의 재현성이 좋다는 것을 의미하며, 따라서 측정 정밀도를 높인다는 것을 알 수 있다.
본 발명은 특정 점영상의 중심점을 추출할 때의 정확도와 정밀도 향상에 대해 기술하였지만, 파면왜곡 장치 이외에도 레이저 패턴 영사 장치나, 삼각측정 장치 등 거의 모든 레이저 패턴 처리 방식에 적용하여 최적화함으로써 측정 분해능 향상에 좋은 결과를 얻을 수 있다.
본 발명은 상술한 특정의 바람직한 실시예에 한정되지 아니하며, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 당해 발명이 속하는 기술분야에서 통상의 지식을 가진 자라면 누구든지 다양한 변형실시가 가능한 것은 물론이고, 그와 같은 변경은 청구범위 기재의 범위 내에 있게 된다.
상기와 같은 본 발명은 기존의 무게중심법에 영상성분의 위치에 따른 가중치를 이용하여 신호 성분의 가중치는 크게 하고, 노이즈 성분의 가중치는 작게 하는 방법으로, 구체적으로 계산 대상이 되는 점영상 신호의 그레이 레벨 값을 그 대로 사용하는 것이 아니라, 제곱, 세제곱 등의 지수 승을 사용함으로써, 큰 신호의 기여도는 살리고 작은 노이즈 성분의 기여도는 저하시키는 방법과 대략적으로 구한 중심의 위치에서 가까운 정보는 가중치를 주어 주변의 노이즈 성분을 죽이는 효과 를 갖고 있어서, 점 영상의 중심점 추출 정확도와 재현성을 향상되었다는 장점을 가진 유용한 발명으로 산업상 그 이용이 크게 기대되는 발명인 것이다.

















Claims (4)

  1. 삭제
  2. 레이저빔의 파면왜곡이나 광학계의 수차 등을 측정하기 위한 하트만 방식의 측정방법에 있어서,
    측정 대상이 되는 빛이 빔 확대기를 통해 측정기기에 적절한 크기로 입사시킨 후 배열렌즈에 의해 빛이 여러 개의 작은 구경으로 나뉘어 각각의 구경에 대한 초점을 맺도록 하는 단계와;
    각각의 초점 영상은 렌즈를 통하여 카메라 등의 영상획득 장치로 입사되어 영상을 획득하는 단계와;
    획득된 영상에 있는 점 영상의 중심점들의 중심 위치를 정확히 판단하기 위하여 각 점영상의 세기에 가중치를 주어 가중치를 준 중심을 구함으로써 점 영상의 주변 노이즈를 제거한 위치를 구하는 단계와;
    이렇게 구한 중심점의 위치를 기준이 되는 중심점의 위치와 비교함으로써 그 이동량으로부터 파면왜곡 정도를 측정하는 단계로 이루어진 것을 특징으로 하는 분해능을 향상시킨 파면왜곡 측정 방법.
  3. 제 2항에 있어서,
    획득된 점영상의 중심점을 추출할 때 각영상의 세기의 제곱 또는 세제곱이 되는 양의 무게 중심을 구함으로써 큰 세기 값의 효과를 더욱 크게 하여 분해능을 향상시키는 단계로 구성되는 것을 특징으로 하는 분해능을 향상시킨 파면왜곡 측정 방법.
  4. 제 2항에 있어서,
    획득된 점영상의 중심점을 무게중심법으로 추출할 때 일차적으로 구한 중심점에 가까운 위치에 있는 세기값은 2 또는 3의 가중치를 주어서 더 큰 효과를 내도 록 함으로써 파면측정 장치의 분해능을 높이는 단계로 구성되는 것을 특징으로 하는 분해능을 향상시킨 파면왜곡 측정 방법.
KR1020040105272A 2004-12-14 2004-12-14 분해능을 향상시킨 파면왜곡 측정 방법 KR100620807B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020040105272A KR100620807B1 (ko) 2004-12-14 2004-12-14 분해능을 향상시킨 파면왜곡 측정 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020040105272A KR100620807B1 (ko) 2004-12-14 2004-12-14 분해능을 향상시킨 파면왜곡 측정 방법

Publications (2)

Publication Number Publication Date
KR20060066802A KR20060066802A (ko) 2006-06-19
KR100620807B1 true KR100620807B1 (ko) 2006-09-12

Family

ID=37161482

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020040105272A KR100620807B1 (ko) 2004-12-14 2004-12-14 분해능을 향상시킨 파면왜곡 측정 방법

Country Status (1)

Country Link
KR (1) KR100620807B1 (ko)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101404690B1 (ko) * 2014-04-30 2014-06-09 국방과학연구소 레이저 파면센서의 기준파면 점영상 위치 측정 방법

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07332952A (ja) * 1994-06-03 1995-12-22 Fuji Photo Optical Co Ltd 干渉計による球面測定解析方法
JP2000097609A (ja) 1998-09-21 2000-04-07 Nikon Corp 波面干渉計
JP2000275005A (ja) 1999-03-24 2000-10-06 Fuji Photo Optical Co Ltd 波面形状測定方法
US6578963B2 (en) 2000-04-19 2003-06-17 Alcon Universal Ltd. Wavefront sensor for objective measurement of an optical system and associated methods
KR20030080765A (ko) * 2002-04-10 2003-10-17 한국원자력연구소 다중 입사광의 파면왜곡 정밀측정 장치 및 그 방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07332952A (ja) * 1994-06-03 1995-12-22 Fuji Photo Optical Co Ltd 干渉計による球面測定解析方法
JP2000097609A (ja) 1998-09-21 2000-04-07 Nikon Corp 波面干渉計
JP2000275005A (ja) 1999-03-24 2000-10-06 Fuji Photo Optical Co Ltd 波面形状測定方法
US6578963B2 (en) 2000-04-19 2003-06-17 Alcon Universal Ltd. Wavefront sensor for objective measurement of an optical system and associated methods
KR20030080765A (ko) * 2002-04-10 2003-10-17 한국원자력연구소 다중 입사광의 파면왜곡 정밀측정 장치 및 그 방법

Also Published As

Publication number Publication date
KR20060066802A (ko) 2006-06-19

Similar Documents

Publication Publication Date Title
US8432479B2 (en) Range measurement using a zoom camera
US8773550B2 (en) Range measurement using multiple coded apertures
CN110431404A (zh) 表面缺陷检查方法及表面缺陷检查装置
CN111083365B (zh) 一种最佳焦平面位置快速检测方法及装置
US8260074B2 (en) Apparatus and method for measuring depth and method for computing image defocus and blur status
CN108401114A (zh) 焦点检测装置
US20110267485A1 (en) Range measurement using a coded aperture
CN110520768B (zh) 高光谱光场成像方法和系统
CN110261069B (zh) 一种用于光学镜头的检测方法
Flesia et al. Sub-pixel straight lines detection for measuring through machine vision
KR20160034928A (ko) 주요지점 식별
CN115375924A (zh) 一种基于图像识别的桥梁健康监测方法与系统
CN114565565A (zh) 一种视觉测量目标中心亚像素定位方法
CN104813217A (zh) 用于设计能够估计景深的无源单通道成像器的方法
CN117073988B (zh) 抬头显示虚像距离测量系统及方法、电子设备
CN107478174B (zh) 一种针对暗弱信号的夏克哈特曼探测器质心探测方法
CN113298700A (zh) 一种在散射场景中的高分辨图像重构方法
CN113188671A (zh) 一种基于交叉迭代自动位置矫正的波前检测方法
KR100620807B1 (ko) 분해능을 향상시킨 파면왜곡 측정 방법
RU2680657C1 (ru) Способ определения волновых аберраций оптической системы
JP2006528499A (ja) オンライン波面測定および表示
CN116952357A (zh) 基于线-面相机结合的分光成像视觉振动测量系统及方法
CN111222544B (zh) 一种卫星颤振对相机成像影响的地面模拟测试系统
Baer Circular-edge spatial frequency response test
Hong et al. Phase retrieval based on transport of intensity equation and image interpolation

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20110711

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20120711

Year of fee payment: 7

LAPS Lapse due to unpaid annual fee