KR20110085785A - 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치 - Google Patents

거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치 Download PDF

Info

Publication number
KR20110085785A
KR20110085785A KR1020100005753A KR20100005753A KR20110085785A KR 20110085785 A KR20110085785 A KR 20110085785A KR 1020100005753 A KR1020100005753 A KR 1020100005753A KR 20100005753 A KR20100005753 A KR 20100005753A KR 20110085785 A KR20110085785 A KR 20110085785A
Authority
KR
South Korea
Prior art keywords
light
phase delay
polynomial
gain
projection
Prior art date
Application number
KR1020100005753A
Other languages
English (en)
Other versions
KR101675112B1 (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 KR1020100005753A priority Critical patent/KR101675112B1/ko
Priority to US12/837,814 priority patent/US8611610B2/en
Priority to JP2011011423A priority patent/JP5739174B2/ja
Publication of KR20110085785A publication Critical patent/KR20110085785A/ko
Application granted granted Critical
Publication of KR101675112B1 publication Critical patent/KR101675112B1/ko

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/521Depth or shape recovery from laser ranging, e.g. using interferometry; from the projection of structured light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • G01S13/42Simultaneous measurement of distance and other co-ordinates
    • G01S13/426Scanning radar, e.g. 3D radar
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20056Discrete and fast Fourier transform, [DFT, FFT]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Optics & Photonics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Measurement Of Optical Distance (AREA)

Abstract

거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치가 개시된다. 개시된 방법에 따르면, 수학적으로 이상적인 파형이 아닌 실제의 비선형적 파형을 이용하여 광학 장치와 피사체 사이의 거리 정보를 추출한다. 따라서, 사용 파형에 제한을 받지 않고 정확한 거리 정보를 추출할 수 있으며, 왜곡 및 비선형성이 거의 없는 고가의 광원이나 광변조 소자를 사용할 필요가 없고, 복잡한 오차 보상 수단이 요구되지 않는다. 따라서, 기존의 광원, 광변조 소자 및 광학 장치를 그대로 이용할 수 있어서 추가 비용이 소요되지 않는다. 더욱이, 미리 계산된 거리 정보가 저장되어 있는 룩업 테이블을 사용하기 때문에, 거리 정보를 추출하는데 있어서 연산량이 매우 작으므로 실시간 거리 정보 영상의 촬영이 가능하다.

Description

거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치{Method of extractig depth information and optical apparatus employing the method}
개시된 내용은 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치에 관한 것으로, 더욱 상세하게는 수학적으로 이상적인 파형이 아닌 실제의 비선형적 파형을 이용하여 광학 장치와 피사체 사이의 거리 정보를 추출하는 방법 및 상기 방법을 채용한 광학 장치에 관한 것이다.
피사체(object)와의 거리 정보를 획득할 수 있는 3D 카메라나 레이저 레이더(Laser Radar; LADAR)에 관한 연구가 최근 증가하는 추세이다. 일반적으로, 피사체에 대한 거리 정보는, 두 대의 카메라를 이용한 양안 입체시(Stereo Vision) 방법이나 구조광(Structured Light)과 카메라를 이용한 삼각 측량법(Triangulation)을 이용하여 얻을 수 있다. 그러나 이러한 방법은 피사체의 거리가 멀어질수록 거리 정보에 대한 정확도가 급격히 저하되고 피사체의 표면 상태에 의존적이어서 정밀한 거리 정보를 얻기 어렵다.
이러한 문제를 개선하기 위하여 광시간비행법(Time-of-Flight; TOF)이 도입되었다. TOF 방법은 레이저 빔을 피사체에 조사한 후, 피사체로부터 반사되는 광이 수광부에서 수광되기까지의 광 비행시간을 측정하는 방법이다. TOF 방법에 따르면, 기본적으로 특정 파장의 빛(예컨대, 850nm의 근적외선)을 발광 다이오드(LED) 또는 레이저 다이오드(LD)를 이용하여 피사체에 투사하고, 피사체로부터 반사된 동일한 파장의 빛을 수광부에서 수광한 후, 거리 정보를 추출하기 위한 특별한 처리 과정을 거치게 된다. 이러한 일련의 광 처리 과정에 따라 다양한 TOF 방법이 소개되었다.
예를 들어, SLP(Shuttered Light Pulse) 방법은 피사체로부터 반사된 영상을 영상증배관(Image Intensifier) 또는 다른 특별한 고체 변조기 소자를 이용하여 광변조한 후, 광변조된 영상을 이미지 센서로 촬영하여 그 세기 값으로부터 거리 정보를 얻는 방법이다. 이 방법의 경우, 거리에 따른 빛의 위상차 또는 이동 시간을 식별하기 위해 수십~수백 MHz의 초고속의 광변조 속도가 필요하다. 이를 위해, 예를 들어 MCP(Multi-Channel Plate)를 구비한 영상증배관을 사용하거나, GaAs 계열의 고체 변조기 소자를 사용한다. 최근에는, 특성이 개선된 GaAs 기반의 변조기 소자 및 전광(Electro-Optic) 물질을 이용한 박형의 변조기 소자가 제시되었다.
한편, 거리 정보 추출을 위한 광처리 과정에 있어서, 광원과 광변조 소자를 펄스 구동하는 방법, 삼각형파(예컨대, 램프(ramp) 파형) 등의 특별한 파형을 이용하는 방법, 또는 사인파를 사용하는 방법 등이 소개되었다. 또한 각각의 사용 파형에 따라 광원과 광변조 소자의 다양한 구동 방법이 소개되어 있으며, 촬영된 세기 값으로부터 거리 정보를 추출하기 위한 다양한 알고리즘이 소개되어 있다.
그런데, 상술한 TOF 방법에서는, 사용된 특정한 파형 즉, 펄스, 삼각형파, 사인파 등이 이상적인 파형을 갖는 것으로 가정하고 거리 정보를 추출하였다. 그러나 실제로는, 광원 또는 광변조 소자의 작동 오차 및 비선형성 등으로 인하여 정확한 이상적인 파형의 구현이 실질적으로 어려울 수가 있다. 예를 들어, LED 광원은 문턱 전류 이상에서 작동하며, 작동 구간 내에서도 입력 전류에 대한 출력 광파워의 비선형성 및 포화(saturation) 등이 존재 한다. 광변조 소자도 역시 작동 구간 내에서 광변조의 이득(gain)이 선형적이지 않은 경우가 많다. 특히, 고전압 또는 고전류를 사용하는 경우에, 구동 회로의 설계에 따라 비선형성 또는 작동 오차가 커질 수 있다. 이러한 작동 오차는 거리 정보 추출 결과에 직접적인 영향을 주며, 따라서 오차의 요인으로 작용한다. 이에 따라, 매우 정확한 이상적인 파형을 발생시킬 수 있는 고가의 파형 발생기를 사용하거나, 또는 오차를 보상하기 위한 복잡한 추가적인 알고리즘 및 구동 회로가 요구될 수 있다.
수학적으로 이상적인 파형이 아닌 실제의 비선형적 파형을 이용하여 광학 장치와 피사체 사이의 거리 정보를 추출하는 방법을 제공한다.
또한, 상기 방법을 채용하여 거리 정보를 추출할 수 있는 광학 장치를 제공한다.
일 유형에 따른 거리 정보 추출 방법은, 적어도 3개의 상이한 투사광을 순차적으로 피사체에 각각 투사하는 단계; 피사체에서 반사된 반사광을 이득 파형을 갖는 광변조기로 변조하는 단계; 상기 변조된 반사광을 광검출기에서 검출하여 상기 적어도 3개의 상이한 투사광에 각각 대응하는 적어도 3개의 세기 영상을 얻는 단계; 상기 적어도 3개의 세기 영상을 이용하여 반사광의 위상 지연을 구하는 단계; 및 상기 위상 지연으로부터 피사체까지의 거리를 계산하는 단계를 포함할 수 있다. 여기서, 상기 위상 지연을 구하는 단계는 상기 적어도 3개의 상이한 투사광과 이득 파형을 미리 실측하여 구성한 룩업 테이블을 참조하여 이루어질 수 있다.
여기서, 상기 적어도 3개의 상이한 투사광과 상기 광변조기의 이득 파형은 동일한 주기를 가질 수 있다.
예컨대, 상기 적어도 3개의 상이한 투사광은 파형이 동일하고 위상이 서로 다르거나 또는 파형이 서로 다를 수 있다.
상기 적어도 3개의 상이한 투사광은 그 평균값이 서로 같도록 제어될 수 있다.
상기 광검출기는 상기 적어도 3개의 상이한 투사광과 각각 대응하는 적어도 3개의 변조된 반사광을 소정의 노출 시간 동안 각각 누적하여 상기 적어도 3개의 상이한 투사광에 각각 대응하는 적어도 3개의 세기 영상을 형성할 수 있다.
상술한 거리 정보 추출 방법은 상기 위상 지연으로부터 피사체의 반사도, 및 투과광/반사광 이외의 성분인 외광의 크기를 구하는 단계를 더 포함할 수 있다.
상기 위상 지연을 구하는 단계는, 상기 적어도 3개의 세기 영상을 이용하여 3차원 공간 상의 세기 영상 벡터를 구성하는 단계; 상기 세기 영상 벡터가 3차원 공간 상의 기준 벡터로부터 기울어진 각도를 구하는 단계; 및 상기 기울어진 각도로 상기 룩업 테이블을 참조하여 반사광의 위상 지연을 구하는 단계를 포함할 수 있다.
이 경우, 상기 룩업 테이블은 상기 기울어진 각도와 반사광의 위상 지연 사이의 미리 결정된 관계를 포함할 수 있다.
또한, 상기 세기 영상 벡터는 상기 적어도 3개의 세기 영상들의 크기들 사이의 차들의 조합으로 구성될 수 있다.
예를 들어, 상기 적어도 3개의 투사광과 이득 파형의 상이한 조합별로 각각 별도의 룩업 테이블이 존재할 수 있으며, 상기 적어도 3개의 투사광과 이득 파형의 특정한 조합에 따라 그에 해당하는 룩업 테이블을 참조할 수 있다.
한편, 상기 룩업 테이블의 구성 방법은, 상기 적어도 3개의 상이한 투사광과 상기 광변조기의 이득 파형을 측정하여 각각의 투사광과 이득 파형을 나타내는 다항식을 각각 구하는 단계; 상기 적어도 3개의 상이한 투사광과 상기 광변조기의 이득 파형에 관한 다항식을 이용하여 적어도 3개의 세기 영상에 관한 다항식을 구하는 단계; 상기 적어도 3개의 세기 영상에 관한 다항식을 각각 위상 지연에 관한 함수와 위상 지연과 관련 없는 상수항의 합으로 정리하는 단계; 상기 위상 지연에 관한 적어도 3개의 함수들 사이의 차들의 조합을 이용하여 3차원 공간 상에 위상 지연에 관한 벡터 함수를 정의하고, 특정한 벡터를 기준 벡터로서 정의하는 단계; 상기 적어도 3개의 세기 영상에 관한 다항식들 사이의 차들의 조합으로 3차원 공간 상에 세기 영상 벡터를 정의하는 단계; 세기 영상 벡터가 기준 벡터로부터 기울어진 각도와 상기 각도에 대응하는 위상 지연 사이의 관계를 해석적으로 결정하는 단계; 및 상기 해석적으로 결정된 관계를 표 또는 함수로 나타내는 단계를 포함할 수 있다.
여기서, 상기 세기 영상에 관한 다항식을 구하는 단계는, 상기 적어도 3개의 상이한 투사광이 피사체에서 반사된 후의 반사광에 관한 일반적인 다항식을 상기 적어도 3개의 상이한 투사광에 관한 다항식으로부터 유도하는 단계; 상기 반사광에 관한 다항식과 상기 이득 파형에 관한 다항식을 곱하여 광변조된 광에 대한 다항식을 구하는 단계; 및 상기 광변조된 광에 대한 다항식을 광검출기에서의 노출 시간으로 적분하여 상기 적어도 3개의 상이한 투사광에 각각 대응하는 적어도 3개의 세기 영상에 관한 다항식을 구하는 단계를 포함할 수 있다.
또한, 상기 적어도 3개의 상이한 투사광과 상기 광변조기의 이득 파형을 측정하여 각각의 투사광과 이득 파형을 나타내는 다항식을 각각 구하는 단계는, 상기 적어도 3개의 상이한 투사광을 각각 측정하는 단계; 상기 각각의 투사광을 측정한 결과를 시간 축상 직교하는 함수로 나타내는 단계; 상기 광변조기의 이득 파형을 측정하는 단계; 및 상기 이득 파형을 측정한 결과를 시간 축상 직교하는 함수로 나타내는 단계를 포함할 수 있다.
여기서, 상기 각각의 투사광을 측정한 결과를 시간 축상 직교하는 함수로 나타내는 단계는, 푸리에 급수의 사인항의 계수와 코사인항의 계수 또는 체비셰프 다항식의 계수를 결정함으로써 각각의 투사광에 대한 다항식을 구하는 단계를 포함할 수 있다.
또한, 상기 이득 파형을 측정한 결과를 시간 축상 직교하는 함수로 나타내는 단계는, 푸리에 급수의 사인항의 계수와 코사인항의 계수 또는 체비셰프 다항식의 계수를 결정함으로써 이득 파형에 대한 다항식을 구하는 단계를 포함할 수 있다.
또한, 상기 룩업 테이블의 구성 방법은, 상기 각각의 투사광과 이득 파형을 측정한 결과로부터 각각의 투사광의 평균값과 이득 파형의 평균값을 결정하는 단계를 더 포함할 수 있다.
상기 적어도 3개의 세기 영상에 관한 다항식을 각각 위상 지연에 관한 함수와 위상 지연과 관련 없는 상수항의 합으로 정리하는 단계는 상기 적어도 3개의 세기 영상에 관한 다항식에 존재하는 직류 항만을 고려할 수 있다.
예를 들어, 상기 차원 공간 상의 위상 지연에 관한 벡터 함수는
Figure pat00001
와 같이 정의될 수 있으며, 여기서
Figure pat00002
이고(i,j=1,2,3), F(i)TOF)는 i번째 투사광에 대한 위상 지연에 관한 함수를 나타낸다.
예를 들어, 상기 위상 지연에 관한 벡터 함수 중에서 위상 지연 ΦTOF=0인 벡터
Figure pat00003
를 기준 벡터로 정의할 수 있다.
또한, 상기 3차원 공간 상의 세기 영상 벡터는
Figure pat00004
와 같이 정의될 수 있으며, 여기서
Figure pat00005
이고(i,j=1,2,3), I(i) CCD는 i번째 투사광에 대한 세기 영상을 나타낸다.
한편, 다른 유형에 따른 거리 정보 추출 방법은, 하나의 주기적인 파형을 갖는 투사광을 피사체에 투사하는 단계; 피사체에서 반사된 반사광을 적어도 3개의 상이한 이득 파형을 갖는 광변조기로 순차적으로 변조하는 단계; 상기 변조된 반사광을 광검출기에서 검출하여 상기 적어도 3개의 상이한 이득 파형에 각각 대응하는 적어도 3개의 세기 영상을 얻는 단계; 상기 적어도 3개의 세기 영상을 이용하여 반사광의 위상 지연을 구하는 단계; 및 상기 위상 지연으로부터 피사체까지의 거리를 계산하는 단계를 포함할 수 있다. 여기서, 상기 위상 지연을 구하는 단계는 상기 투사광과 상기 적어도 3개의 상이한 이득 파형을 미리 실측하여 구성한 룩업 테이블을 참조할 수 있다.
또한, 상술한 거리 정보 추출 방법을 수행하기 위한 일 유형에 따른 광학 장치는, 주기적인 투사광을 발생시키는 광원; 피사체로부터 반사된 광을 주기적인 이득 파형으로 광변조하기 위한 광변조기; 상기 광변조기에 의해 광변조된 영상을 검출하기 위한 광검출기; 및 상기 광검출기의 출력을 기초로 피사체와의 거리를 계산하기 위한 거리정보영상 처리기를 포함할 수 있다.
또한 상기 광학 장치는, 상기 광원을 구동시켜 투사광의 파형을 제어하는 광원 구동부; 상기 광변조기를 구동시켜 이득 파형을 제어하는 광변조기 구동부; 및 상기 광원 구동부, 광변조기 구동부 및 광검출기의 동작을 제어하기 위한 제어부를 더 포함할 수 있다.
또한 상기 광학 장치는, 상기 광변조기의 광입사면에서 반사광을 광변조기의 영역 내에 집광하는 제 1 렌즈; 상기 제 1 렌즈와 광변조기 사이에서 소정의 파장을 갖는 광만을 투과시키는 필터; 및 상기 광변조기와 광검출기 사이에서 상기 광변조된 영상을 광검출기의 영역 내에 집광하는 제 2 렌즈를 더 포함할 수 있다.
또한, 상기 광검출기는 2차원 또는 1차원 어레이를 갖는 CCD 이미지 센서, CMOS 이미지 센서, 포토다이오드 어레이를 사용하거나 또는 단일 지점의 거리 측정을 위한 1개의 포토다이오드를 사용할 수 있다.
이러한 구성의 광학 장치는 예를 들어 3D 카메라나 레이저 레이더(LADAR)에서 사용될 수 있다.
한편, 또 다른 유형의 룩업 테이블 구성 방법은, 하나의 투사광 및 광변조기의 적어도 3개의 상이한 이득 파형을 측정하여 각각의 투사광과 이득 파형을 나타내는 다항식을 각각 구하는 단계; 상기 투사광과 상기 광변조기의 적어도 3개의 상이한 이득 파형에 관한 다항식을 이용하여 적어도 3개의 세기 영상에 관한 다항식을 구하는 단계; 상기 적어도 3개의 세기 영상에 관한 다항식을 각각 위상 지연에 관한 함수와 위상 지연과 관련 없는 상수항의 합으로 정리하는 단계; 상기 위상 지연에 관한 적어도 3개의 함수들 사이의 차들의 조합을 이용하여 3차원 공간 상에 위상 지연에 관한 벡터 함수를 정의하고, 특정한 벡터를 기준 벡터로서 정의하는 단계; 상기 적어도 3개의 세기 영상에 관한 다항식들 사이의 차들의 조합으로 3차원 공간 상에 세기 영상 벡터를 정의하는 단계; 세기 영상 벡터가 기준 벡터로부터 기울어진 각도와 상기 각도에 대응하는 위상 지연 사이의 관계를 해석적으로 결정하는 단계; 및 상기 해석적으로 결정된 관계를 표 또는 함수의 형태로 나타내는 단계를 포함할 수 있다.
또한, 일 유형에 따른 저장 매체는 상술한 방법으로 구성된 룩업 테이블을 저장할 수 있다.
개시된 방법에 따르면, 수학적으로 이상적인 파형이 아닌 실제의 비선형적 파형을 고려하여 거리 정보를 추출하기 때문에, 사용 파형에 제한을 받지 않고 정확한 거리 정보를 추출할 수 있다. 그 결과, 왜곡 및 비선형성이 거의 없는 고가의 광원이나 광변조 소자를 사용할 필요가 없으며, 복잡한 오차 보상 수단이 요구되지 않는다. 따라서, 기존의 광원, 광변조 소자 및 광학 장치를 그대로 이용할 수 있어서 추가 비용이 소요되지 않는다. 더욱이, 미리 계산된 거리 정보가 저장되어 있는 룩업 테이블을 사용하기 때문에, 거리 정보를 추출하는데 있어서 연산량이 매우 작으므로 실시간 거리 정보 영상(depth image)의 촬영이 가능하다.
도 1은 광시간비행법(TOF)을 이용하여 거리 정보를 추출할 수 있는 광학 장치의 예시적인 구조를 개략적으로 도시한다.
도 2는 광원에서 실제로 출력된 출력광과 이상적인 사인파 출력광을 예시한다.
도 3은 광변조기의 투과율의 시간축 파형을 예시적으로 도시한다.
도 4는 광원에서 발생하는 3개의 투사광의 형태를 예시적으로 도시한다.
도 5는 3개의 상이한 반사광을 변조한 후 광검출기에서 3개의 상이한 세기 영상을 생성하는 과정을 도시한다.
도 6은 3개의 세기 영상으로 거리 정보를 갖는 영상을 생성하는 과정을 개략적으로 도시한다.
도 7은 하나의 동일한 투사광과 3개의 상이한 변조로 3개의 상이한 세기 영상을 생성하는 과정을 도시한다.
도 8은 룩업 테이블의 구성을 위해 위상 지연을 변수로서 갖는 벡터 함수와 세기 영상에 관한 벡터를 3차원 공간 상에 표시한 예를 도시한다.
도 9는 룩업 테이블을 예시적으로 도시하는 그래프이다.
도 10은 룩업 테이블을 구성하는 과정을 개략적으로 보이는 흐름도이다.
도 11은 룩업 테이블을 이용하여 거리 정보를 추출하는 과정을 개략적으로 보이는 흐름도이다.
도 12는 시뮬레이션을 통해 거리 정보를 추출한 결과를 참값과 비교한 그래프이다.
도 13은 도 12의 거리 정보 추출 결과와 참값과의 오차를 나타낸 그래프이다.
도 14는 실험을 위해 레이저 다이오드의 광출력을 측정한 결과이다.
도 15는 도 14의 레이저 다이오드의 광출력 스펙트럼을 고속 푸리에 변환으로 분석한 결과이다.
도 16은 광변조기의 변조 이득 파형을 측정한 결과이다.
도 17은 도 16의 이득 파형 스펙트럼을 고속 푸리에 변환을 통해 분석한 결과이다.
도 18에는 개시된 거리 정보 추출 알고리즘에 따라 파형의 왜곡을 고려한 결과와 이상적인 파형만을 고려한 결과를 비교하는 그래프이다.
이하, 첨부된 도면들을 참조하여, 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치에 대해 상세하게 설명한다. 이하의 도면들에서 동일한 참조부호는 동일한 구성요소를 지칭하며, 도면상에서 각 구성요소의 크기는 설명의 명료성과 편의상 과장되어 있을 수 있다.
1. 광학 장치의 구조 및 동작
먼저, 도 1은 광시간비행법(TOF)을 이용하여 거리 정보를 추출할 수 있는 광학 장치(100)의 예시적인 구조를 개략적으로 도시하고 있다. 도 1을 참조하면, 상기 광학 장치(100)는 소정의 파장을 갖는 광을 발생시키는 광원(101), 상기 광원(101)을 구동시키기 위한 광원 구동부(102), 피사체(200)로부터 반사된 광을 광변조하기 위한 광변조기(103), 광변조기(103)를 구동시키기 위한 광변조기 구동부(104), 광변조기(103)에 의해 광변조된 영상을 검출하기 위한 광검출기(105), 상기 광원 구동부(102), 광변조기 구동부(104) 및 광검출기(105)의 동작을 제어하기 위한 제어부(106), 및 광검출기(105)의 출력을 기초로 거리 정보를 계산하기 위한 거리정보영상 처리기(depth image processor)(107)를 포함할 수 있다. 또한, 상기 광변조기(103)의 광입사면에는 반사광을 광변조기(103)의 영역 내에 집광하기 위한 제 1 렌즈(108)와 소정의 파장을 갖는 광만을 투과시키는 필터(109)가 배치될 수 있다. 광변조기(103)와 광검출기(105) 사이에는 광변조된 영상을 광검출기(105)의 영역 내에 집광하기 위한 제 2 렌즈(110)가 더 배치될 수 있다.
광원(101)은, 예를 들어, 안전을 위해 인간의 눈에는 보이지 않는 약 850nm의 근적외선(NIR) 파장을 갖는 광을 방출시킬 수 있는 발광 다이오드(LED) 또는 레이저 다이오드(LD)일 수 있지만, 파장의 대역과 광원의 종류는 제한을 받지 않는다. 광원 구동부(102)는 제어부(106)로부터 수신된 제어 신호에 따라 광원(101)을 예를 들어 크기 변조(Amplitude Modulation) 방식으로 구동할 수 있다. 이에 따라, 광원(101)으로부터 피사체(200)로 투사되는 투사광은 소정의 주기(Te)를 갖는 주기적인 연속 함수의 형태를 가질 수 있다. 예를 들어, 투사광은 사인파, 램프파, 사각파 등과 같이 특수하게 정의된 파형을 가질 수도 있지만, 정의되지 않은 일반적인 형태의 파형을 가질 수도 있다.
앞서 설명한 바와 같이, 지금까지 제안된 거리 정보 추출 방법들의 경우, 광원이 이상적인 형태의 사인파, 램프파, 사각파 등을 투사광으로서 출력하는 것으로 가정하였다. 그러나 실제의 광원은, 광원 자체의 동작 비선형성 및 신호 왜곡 등에 의해, 이상적인 파형을 구동신호로 입력하는 경우에도 이상적인 파형을 그대로 출력하지 않는다. 이는 거리 정보를 추출하는데 있어서 오차의 요소로서 작용한다. 도 2는 광원에서 실제로 출력된 출력광과 이상적인 사인파 출력광을 예시하고 있다. 도 2에서 점선으로 표시된 것은 이상적인 사인파 출력광을 나타내며, 실선으로 표시된 것은 실제의 출력광을 나타낸다. 도 2에 도시된 실제 출력광을 보면, 이상적인 사인파에 비하여, 실제 출력광은 왜곡과 DC-오프셋을 갖는다는 것을 알 수 있다. 이에 따라, 측정 정밀도를 높이기 위하여, 이상적인 파형을 발생시킬 수 있는 고가의 파형 발생기나 고가의 광원이 요구될 수 있다.
도 1에 도시된 광학 장치(100)의 경우, 후술하는 바와 같이 실제의 출력 파형을 측정하고 그 측정 결과를 거리 정보 계산에 적용한다. 따라서, 상기 광학 장치(100)는 왜곡과 DC-오프셋이 있는 출력광을 그대로 이용할 수 있으므로, 상술한 고가의 파형 발생기나 고가의 광원을 사용하지 않는다. 또한, 상기 광학 장치(100)의 광원(101)은 단일한 파형이 아닌 사인파, 램프파, 사각파 등을 조합한 일반적인 파형을 갖는 광을 투사광으로서 출력할 수도 있다.
한편 광변조기(103)는 피사체(200)로부터 반사된 광을 광변조기 구동부(104)의 제어에 따라 광변조한다. 광변조기 구동부(104)는 제어부(106)로부터 수신된 제어 신호에 따라 광변조기(103)를 구동시킨다. 예를 들어, 광변조기(103)는 광변조기 구동부(104)에 의해 제공된 소정의 파형을 갖는 광변조 신호에 따라 이득을 변화시켜 반사광의 크기를 변조시킬 수 있다. 이를 위하여, 광변조기(103)는 가변 이득(gain)을 갖는 투과 필터의 역할을 할 수 있다. 광변조기(103)는 거리에 따른 빛의 위상차 또는 이동 시간을 식별하기 위해 수십~수백 MHz의 높은 광변조 속도로 동작할 수 있다. 이에 부합하는 광변조기(103)로서, 예를 들어 MCP(Multi-Channel Plate)를 구비한 영상증배관, GaAs 계열의 고체 변조기 소자, 전광(Electro-Optic) 물질을 이용한 박형의 변조기 소자 등이 사용될 수 있다. 도 1에는 광변조기(103)가 투과형인 것으로 도시되어 있지만, 반사형 광변조기를 사용하는 것도 가능하다.
도 3은 광변조기(103)의 주기적인 이득 변화 즉, 투과율의 시간축 파형을 예시적으로 도시하고 있다. 여기서, 광변조기(103)의 동작 주파수는 광원(101)의 동작 주파수(fe)와 같도록 구성된다. 도 3에서 점선으로 표시된 것은 광변조기(103)의 이상적인 이득 변화를 나타내며, 실선으로 표시된 것은 광변조기(103)의 실제 이득 변화를 나타낸다. 광원(101)과 마찬가지로, 광변조기(103)도 역시 소자 구조에 따라 비선형성과 DC-오프셋을 갖는다. 따라서, 도 3에 도시된 바와 같이, 광변조기(103)의 실제 이득 변화에는 이상적인 광변조기의 이득 변화와 달리 왜곡이 존재하게 된다. 지금까지 제안된 거리 정보 추출 방법들의 경우, 광변조기의 이득 변화가 이상적이라고 가정하였기 때문에, 이러한 광변조기(103)의 비선형성과 DC-오프셋도 역시 오차의 원인이 된다. 도 1에 도시된 광학 장치(100)의 경우, 후술하는 바와 같이 광변조기(103)의 실제 이득 변화를 측정하고, 그 측정 결과를 거리 정보 계산에 반영한다. 따라서, 상기 광학 장치(100)는 왜곡이 존재하는 일반적인 광변조기를 그대로 이용할 수 있다. 또한, 광변조기(103)의 이득 파형은 사인파 형태뿐만 아니라 다양한 형태의 파형을 가질 수도 있다.
광검출기(105)는 광변조기(103)에 의해 광변조된 영상을 제어부(106)의 제어에 따라 검출하는 역할을 한다. 만약 피사체(200)의 어느 한 점까지의 거리만을 측정하고자 하는 경우, 광검출기(105)는 예를 들어 포토다이오드나 적분기와 같은 하나의 단일한 광센서를 사용할 수도 있다. 그러나 피사체(200) 상의 다수의 점들까지의 거리들을 동시에 측정하고자 하는 경우에는, 광검출기(105)는 다수의 포토다이오드들의 2차원 또는 1차원 어레이 또는 다른 1차원이나 2차원 광검출 어레이를 가질 수도 있다. 예를 들어, 광검출기(105)는 2차원 또는 1차원 어레이를 갖는 CCD 이미지 센서 또는 CMOS 이미지 센서일 수도 있다. 특히, CCD 이미지 센서 또는 CMOS 이미지 센서와 같은 광검출 어레이를 사용할 경우, 광학 장치(100)는 깊이 정보를 갖는 3D 영상을 촬영할 수 있는 3D 카메라에 적용될 수 있다.
거리정보영상 처리기(107)는 상기 광검출기(105)의 출력을 기초로, 후술할 거리 정보 추출 알고리즘에 따라 거리 정보를 계산하는 역할을 한다. 거리정보영상 처리기(107)는 예를 들어 전용의 집적회로(IC)으로 구현될 수도 있으며, 또는 개인용 컴퓨터와 같은 일반적인 컴퓨터 장치로도 구현될 수도 있다.
이하에서는, 상술한 구조를 갖는 광학 장치(100)의 동작을 개략적으로 설명한다.
먼저, 광원(101)은 제어부(106) 및 광원 구동부(102)의 제어에 따라 소정의 주기(Te) 및 소정의 파형을 갖는 3종류 이상의 상이한 투사광을 서로 다른 시간 동안 차례로 발생시킨다. 이렇게 발생한 3종류 이상의 투사광은 피사체(200)에 차례로 투사된다. 예를 들어 3개의 상이한 투사광을 사용하는 경우, 시간 T1 동안 제 1 투사광을 발생시켜 피사체(200)에 투사시키고, 다음의 시간 T2 동안 제 2 투사광을 발생시켜 피사체(200)에 투사시키고, 이어서 시간 T3 동안 제 3 투사광을 발생시켜 피사체(200)에 투사시킬 수 있다. 3종류 이상의 상이한 투사광은 다양한 방식으로 구현될 수 있다. 예를 들어, 3종류 이상의 투사광은 동일한 파형을 가지면서 단순히 위상만이 다를 수도 있다. 그러나 더욱 정확한 거리 측정을 위하여 파형 자체가 다른 3종류 이상의 투사광을 발생시킬 수도 있다. 파형의 형태 및 주기와 위상의 크기에는 특별한 제한이 없지만, 모든 투사광들은 동일한 주기의 파형을 갖는다.
도 4는 광원(101)에서 발생할 수 있는 3개의 투사광의 형태를 예시적으로 도시하고 있다. 도 4를 참조하면, 제 1 투사광과 제 3 투사광은 비선형적으로 변형된 삼각파의 파형에 가까우며 그 위상이 다르다. 또한, 제 2 투사광은 왜곡된 사각파의 파형에 가깝다. 비교를 위하여, 도 4에는 이상적인 형태의 삼각파가 점선으로 함께 표시되어 있다. 그 외에도, 예를 들어, 사인파, 삼각파, 사각파 및 이들의 조합파와 더불어서 왜곡이 존재하는 파형이 모두 사용 가능하다. 이러한 상이한 파형들의 사용은, 후술할 거리 정보 추출 알고리즘에서 설명하는 바와 같이, 실제의 출력 파형을 측정하고 그 측정 결과를 거리 정보 계산에 적용하기 때문에 가능한 것이다.
광원(101)으로부터 방출된 투사광은 피사체(200)의 표면에서 반사되어, 제 1 렌즈(108)로 입사한다. 일반적으로 실제의 피사체(200)는 광학 장치(100)로부터 촬영되는 표면까지의 거리, 즉 깊이(depth)가 서로 다른 다수의 표면들이 2차원 어레이를 이룰 것이지만, 도 1에는 설명의 단순화를 위하여 깊이(depth)가 서로 다른 5개의 표면(P1~P5)을 갖는 피사체(200)가 예시적으로 도시되어 있다. 투사광이 각각의 표면(P1~P5)에서 반사되면서, 상이하게 시간 지연된(즉, 위상이 상이한) 5개의 반사광이 발생한다. 예를 들어, 제 1 투사광이 피사체(200)의 5개의 표면(P1~P5)에서 각각 반사되면서 위상이 상이한 5개의 제 1 반사광이 발생하며, 제 2 투사광이 동일한 피사체(200)의 5개의 표면(P1~P5)에서 각각 반사되면서 위상이 상이한 5개의 제 2 반사광이 발생하고, 마찬가지로 제 3 투사광도 피사체(200)의 5개의 표면(P1~P5)에서 각각 반사되면서 위상이 상이한 5개의 제 3 반사광이 발생한다. 이때, 광학 장치(100)로부터의 거리가 가장 먼 표면(P1)에서 반사된 반사광은 TOF1 만큼의 시간 지연 후 제 1 렌즈(108)에 도달하고, 광학 장치(100)로부터의 거리가 가장 가까운 표면(P5)에서 반사된 반사광은 TOF1보다 작은 TOF5 만큼의 시간 지연 후 제 1 렌즈(108)에 도달할 것이다.
제 1 렌즈(108)는 반사광을 광변조기(103)의 영역 내에 포커싱한다. 제 1 렌즈(108)와 광변조기(103) 사이에는 사용 파장 이외의 배경광이나 잡광을 제거하기 위하여 소정의 파장을 갖는 광만을 투과시키는 필터(109)가 배치될 수 있다. 예를 들어, 광원(101)이 약 850nm의 근적외선(NIR) 파장을 갖는 광을 방출하는 경우, 상기 필터(109)는 약 850nm의 근적외선 파장 대역을 통과시키는 근적외선 대역 통과 필터(IR band pass Filter)일 수 있다. 따라서, 광변조기(103)에 입사하는 광은 광원(101)으로부터 방출되어 피사체(200)에서 반사된 광이 지배적일 수 있다. 도 1에는 제 1 렌즈(108)와 광변조기(103) 사이에 필터(109)가 배치된 것으로 도시되어 있지만, 제 1 렌즈(108)와 필터(109)의 위치는 서로 바뀔 수도 있다. 예를 들어, 필터(109)를 먼저 통과한 근적외선광이 제 1 렌즈(108)에 의해 광변조기(103)로 포커싱될 수도 있다.
광변조기(103)는 소정의 파형을 갖는 광변조 신호로 반사광의 크기를 변조한다. 예를 들어, 반사광은 도 3에서 예시한 이득을 곱한 양만큼 광변조기(103)에 의해 크기가 변조될 수 있다. 여기서, 이득 파형의 주기는 투사광의 주기(Te)와 같다. 도 1에 도시된 예에서, 광변조기(103)는 피사체(200)의 5개의 표면(P1~P5)에서 각각 반사된 5개의 제 1 반사광을 광변조하여 광검출기(105)에 제공하고, 이어서 5개의 제 2 반사광과 5개의 제 3 반사광을 차례로 광변조하여 광검출기(105)에 제공할 수 있다.
광변조기(103)에 의해 크기가 변조된 광은 제 2 렌즈(110)를 통과하면서 배율 조정 및 재포커싱된 후 광검출기(105)에 도달한다. 따라서, 변조된 광은 제 2 렌즈(110)에 의해 광검출기(105)의 영역 내에 집광된다. 광검출기(105)는 상기 변조된 광을 소정의 노출 시간 동안 수광하여 세기 영상(intensity image)을 생성한다. 예를 들어 광검출기(105)는, 도 5에서 (A)로 표시된 바와 같이, 피사체(200)의 5개의 표면(P1~P5)에서 각각 반사된 5개의 제 1 반사광에 대해 변조된 광을 소정의 노출 시간 동안 수광하여 제 1 세기 영상(IIM1)을 생성한다. 이어서, 도 5에서 (B)로 표시된 바와 같이, 광검출기(105)는 피사체(200)의 5개의 표면(P1~P5)에서 각각 반사된 5개의 제 2 반사광에 대해 변조된 광을 소정의 노출 시간 동안 수광하여 제 2 세기 영상(IIM2)을 생성한다. 다음으로, 도 5에서 (C)로 표시된 바와 같이, 광검출기(105)는 피사체(200)의 5개의 표면(P1~P5)에서 각각 반사된 5개의 제 3 반사광에 대해 변조된 광을 소정의 노출 시간 동안 수광하여 제 3 세기 영상(IIM3)을 생성한다. 이러한 방식으로, 도 5에서 (D)로 표시된 바와 같이, 세 개의 상이한 세기 영상(IIM1,IIM2,IIM3)을 얻을 수 있다. 여기서, 각각의 세기 영상(IIM1,IIM2,IIM3)은 피사체(200)의 5개의 표면(P1~P5)에서 각각 반사된 5개의 반사광에 대한 세기 정보(I1~I5)를 포함하고 있다.
거리 정보를 갖는 한 프레임의 영상을 얻기 위해서는 상술한 3개 이상의 세기 영상들이 필요하다. 도 1 및 도 5에서는, 3개의 투사광을 이용하여 3개의 세기 영상(IIM1,IIM2,IIM3)을 얻는 예가 도시되어 있다. 그러나, 4개 또는 그 이상의 투사광을 이용하여 4개 또는 그 이상의 세기 영상을 얻는 것도 가능하다. 각각의 세기 영상(IIM1,IIM2,IIM3)은 거리 정보를 갖는 한 프레임의 영상을 얻기 위한 서브 프레임(sub-frame)일 수 있다. 만약 한 프레임의 주기를 Td라고 한다면, 3개의 세기 영상(IIM1,IIM2,IIM3) 각각을 얻기 위한 광검출기(105)의 노출 시간(T)은 대략적으로 1/3*Td일 수 있다. 그러나 상기 세기 영상(IIM1,IIM2,IIM3)들을 얻기 위한 노출 시간들이 모두 동일할 필요는 없으며, 설계에 따라 임의로 선택될 수도 있다. 예를 들어, 3개의 세기 영상(IIM1,IIM2,IIM3)에 대해 서로 다른 노출 시간 T1, T2 및 T3가 각각 할당될 수 있으며, 단지 T1 + T2 + T3 ≤ Td 만을 만족하면 된다.
다시 도 5의 (A)를 참조하면, 광원(101)으로부터 방출된 제 1 투사광이 피사체(200)의 5개의 표면(P1~P5)에서 각각 반사되어 5개의 제 1 반사광이 생성된다. 5개의 제 1 반사광은 광변조기(103)에 의해 변조된 후 광검출기(105)에 도달한다. 도 5에는, 설명의 용이성을 위해 광검출기(105)가 5개의 표면(P1~P5)에 각각 대응하는 5개의 화소를 갖는 것으로 도시되어 있다. 5개의 제 1 반사광은 대응하는 5개의 화소에 각각 입사할 수 있다. 도 5의 (A)에 도시된 바와 같이, 표면(P1~P5)에서 각각 반사된 5개의 제 1 반사광은 광학 장치(100)로부터의 거리에 따른 상이한 위상 지연(ΦP1P5)을 각각 갖는다. 예를 들어, 가장 먼 표면(P1)으로부터의 반사광이 가장 큰 위상 지연(ΦP1)을 갖게 된다. 첫번째 서브 프레임에서, 광검출기(105)는 예컨대 T1의 노출 시간 동안 상기 변조된 제 1 반사광을 촬영하여, T1의 시간 동안 누적된 제 1 세기 영상(IIM1)을 생성할 수 있다.
두 번째 서브 프레임은 첫번째 서브 프레임의 촬영 후 T2의 시간 동안 촬영될 수 있다. 도 5의 (B)를 참조하면, 제 2 투사광이 피사체(200)에서 반사되어 생성된 제 2 반사광이 광변조기(103)에 의해 변조된 후 광검출기(105)에 도달한다. 5개의 표면(P1~P5)에서 각각 반사된 5개의 제 2 반사광은 상술한 5개의 제 1 반사광과 동일한 위상 지연(ΦP1P5)들을 각각 갖는다. 광검출기(105)는 예컨대 T2의 노출 시간 동안 상기 변조된 제 2 반사광을 촬영하여, T2의 시간 동안 누적된 제 2 세기 영상(IIM2)을 생성할 수 있다.
세 번째 서브 프레임은 두 번째 서브 프레임의 촬영 후 T3의 시간 동안 촬영될 수 있다. 도 5의 (C)를 참조하면, 제 3 투사광이 피사체(200)에서 반사되어 생성된 제 3 반사광이 광변조기(103)에 의해 변조된 후 광검출기(105)에 도달한다. 5개의 표면(P1~P5)에서 각각 반사된 5개의 제 3 반사광은 상술한 5개의 제 1 반사광과 동일한 위상 지연(ΦP1P5)들을 각각 갖는다. 광검출기(105)는 예컨대 T3의 노출 시간 동안, 상기 변조된 제 3 반사광을 촬영하여, T3의 시간 동안 누적된 제 3 세기 영상(IIM3)을 생성할 수 있다.
상술한 방식으로 순차적으로 얻은 3개의 세기 영상(IIM1,IIM2,IIM3)은 거리정보영상 처리기(107)로 전달된다. 거리정보영상 처리기(107)는 상기 3개의 세기 영상(IIM1,IIM2,IIM3)을 기초로 후술할 알고리즘에 따라 거리 정보를 갖는 영상을 생성할 수 있다. 도 6은 3개의 세기 영상(IIM1,IIM2,IIM3)으로 거리 정보를 갖는 영상을 생성하는 과정을 개략적으로 도시하고 있다. 예를 들어, 첫 번째 프레임에서 상술한 방식으로 생성된 3개의 세기 영상(IIM1,IIM2,IIM3)으로 제 1 거리 정보 영상(depth image)(IMG1)을 얻을 수 있다. 다음으로, 두 번째 프레임에서도 동일한 방식으로 생성된 3개의 세기 영상(IIM1,IIM2,IIM3)으로 제 2 거리 정보 영상(IMG2)을 얻을 수 있으며, 세 번째 프레임에서 같은 방식으로 제 3 거리 정보 영상(IMG3)을 얻을 수 있다. 이러한 동작을 반복함으로써, 소정의 프레임 속도(fd = 1/Td)로 피사체(200)의 다수의 표면까지의 동작을 실시간으로 관찰할 수 있다.
지금까지는 3개 이상의 세기 영상(IIM1,IIM2,IIM3)을 얻기 위하여 3개 이상의 상이한 투사광을 사용하는 경우에 대해 설명하였다. 그러나, 3개 이상의 서브 프레임 동안, 광원(101)은 동일한 한 파형의 투사광을 출력하고, 광변조기(103)가 3개 이상의 상이한 이득 파형으로 반사광을 변조하는 것도 가능하다. 도 7은 하나의 동일한 투사광과 3개의 상이한 변조로 3개의 상이한 세기 영상(IIM1,IIM2,IIM3)을 생성하는 과정을 도시하고 있다. 도 7을 참조하면, 제 1 내지 제 3 반사광은 모두 동일한 파형과 위상을 가지며, 단지 반사된 피사체(200)의 표면(P1~P5)에 따라 상이한 위상 지연(ΦP1P5)들이 존재한다. 또한, 도 7의 (A)~(C)에 각각 도시된 바와 같이, 첫 번째 서브 프레임에서 광변조기(103)는 제 1 광변조 신호로 반사광을 변조하며, 두 번째 서브 프레임에서 광변조기(103)는 제 1 광변조 신호와 다른 제 2 광변조 신호로 반사광을 변조하고, 세 번째 서브 프레임에서 광변조기(103)는 제 1 및 제 2 광변조 신호와 다른 제 3 광변조 신호로 반사광을 변조한다. 이에 따라, 도 7의 (D)에 각각 도시된 바와 같이, 서로 다른 3개의 세기 영상(IIM1,IIM2,IIM3)을 얻을 수 있다. 여기서, 3개의 광변조 신호는 위상만이 서로 다를 수도 있고, 또는 파형 자체가 서로 다를 수도 있다. 이후에 거리 정보 영상을 생성하는 방식은 도 6을 참조로 설명한 것과 같다.
2. 거리 정보 추출 알고리즘
이하에서는, 상술한 3개 이상의 세기 영상(IIM1,IIM2,IIM3)을 이용하여 거리정보영상 처리기(107)가 거리 정보를 추출하는 알고리즘에 대해 설명한다.
이하에서 설명하는 거리 정보 추출 알고리즘의 개요는, 광원(101)에서 방출된 투사광의 실제 파형과 광변조기(103)의 실제 이득 파형을 측정하여 푸리에 급수 등으로 수학적 모델을 세운 다음, 이를 기초로 거리와 세기 영상 사이의 이론적인 관계를 결정하여 룩업 테이블을 구성하고, 실제 촬영하여 얻은 세기 영상을 룩업 테이블과 비교하여 거리 정보를 실시간으로 추출하는 것이다.
이하의 이론적 설명은, 도 5에 도시된 바와 같이, 3종류의 투사광을 사용하고 광변조기(103)의 이득 파형은 한 가지를 사용하는 방법을 위주로 기술한다. 또한, 광검출기(105)의 출력이 2차원 어레이 영상인 경우라도, 각각의 화소에 적용되는 거리 정보 추출 방법은 동일하므로, 하나의 화소에 적용되는 방법만을 설명한다. 다만, 2차원 어레이 영상에서 다수의 화소로부터 거리 정보를 동시에 추출하는 경우에는, 데이터 관리 및 메모리 할당 등을 효율적으로 처리하여 중복 처리되는 부분을 생략함으로써 계산양을 줄일 수도 있을 것이다.
2-1. 투사광과 반사광의 파형, 광변조기의 이득 파형에 대한 모델링
먼저, 주기 Te를 갖는 일반적인 투사광의 파형은 시간 축상 직교하는 함수로 나타낼 수 있다. 예를 들어, 아래의 수학식(1)과 같이, 푸리에 급수의 사인항과 코사인항으로 나타낼 수도 있으며, 또는 체비셰프 다항식(Chebyshev polynomial)과 같은 일반적인 다항식으로 나타낼 수도 있다. 이하에서는 투사광의 파형을 푸리에 급수로 나타내는 경우를 중심으로 설명한다.
Figure pat00006
여기서, s는 투사광의 파형의 종류를 구분하는 식별자이다. 예를 들어, 3개의 투사광을 사용하는 경우에 s = 1,2,3이다. ω는 투사광의 파형의 구동 각주파수로서 ω=2πfe (fe는 투사광의 파형의 구동 주파수, fe = 1/Te)의 관계를 갖는다. 또한, ak (s)는 투사광 (s)에서 k차 사인항의 계수이며, bk (s)는 투사광 (s)에서 k차 코사인항의 계수이다. 상기 계수 ak (s)와 bk (s)는 광원(101)의 광출력을 직접 계측하여 얻을 수 있다. 예를 들어, 광원(101)에서 방출된 투사광의 파형을 측정하고, 측정된 파형에 대한 고속 푸리에 변환(FFT) 등과 같은 푸리에 변환을 통해 주파수 성분을 파악하여 상기 계수들을 구할 수 있다. 이때 차수 m은 광원(101)에서 출력된 투사광의 파형 계측치를 적절히 표현할 수 있도록 충분히 크게 잡는다. 한편, 출력된 투사광의 평균값은 아래의 수학식(2)와 같이 정의된다.
Figure pat00007
여기서, T는 3개의 투사광 중에서 어느 하나에 대응하는 하나의 세기 영상을 얻기 위한 광검출기(105)의 노출 시간, 즉 서브 프레임의 주기이다. 만약 3개의 투사광에 대한 노출 시간이 서로 다르게 설정된다면, 예를 들어 s = 1에 대해 T를 T1으로 놓고, s = 2에 대해 T를 T2로 놓고, s = 3에 대해 T를 T3로 놓을 수 있다. 이하에서는 3개의 투사광에 대해 동일한 노출 시간 T가 설정된 것으로 가정한다.
3개 투사광의 파형을 구현하는 경우에, 제어부(106)는 각각의 투사광에 대한 평균 출력이 같은 값이 되도록 제어 신호를 적절히 조절한다. 도 2와 도 4는 수학식(1)에서 표현된 제 1 내지 제 3 투사광을 예시하는 것으로, 예를 들어 Te=25ns, m=5차항을 사용한 경우이고 투사광의 평균 출력값은 0.5(W)이다. 투사광의 파형으로는 상기 수학식(1)로 표현 가능한 모든 형태의 파형이 제한 없이 사용될 수 있다. 예를 들어, 사인파, 삼각형파, 사각파 또는 이들의 조합, 그리고 왜곡이 존재하는 경우까지 모두 사용 가능하다.
수학식(1)로 표현된 투사광이 피사체(200)에서 반사되어 광학 장치(100)로 다시 입사할 때, 그 반사광은 아래의 수학식(3)과 같은 다항식으로 나타낼 수 있다.
Figure pat00008
광학 장치(100)로 되돌아오는 반사광(s)는 거리에 따른 시간 지연 tTOF 및 피사체(200)의 반사도 r이 투사광(s)에 작용한 것이다. 여기서 반사도 r은 투사광이 피사체(200)에서 반사된 후 광학 장치(100)로 입사되는 비율로 정의될 수 있으며, 피사체(200)와의 거리, 피사체(200)의 표면의 재질 및 기울기, 제 1 렌즈(108)의 크기 등 다양한 요소와 연관되어 있는 변수이다.
Figure pat00009
는 투사광과 관계 없이 광학 장치(100)에 유입되는 외광(ambient light) 성분이다. 거리에 따른 반사광의 위상 지연 ΦTOF는 다음의 수학식(4)와 같이 시간 지연 TOF 및 투사광의 구동 파형의 주기 Te로 나타낼 수 있다.
Figure pat00010
한편, 광변조기(103)의 일반적인 이득 파형도 역시, 아래의 수학식(5)와 같이, 푸리에 급수와 같은 일반적인 다항식으로 나타낼 수 있다.
Figure pat00011
여기서, 사인항과 코사인항의 각 차수의 계수 ck와 dk는 광변조기(103)의 이득 파형을 직접 계측하고, 계측된 이득 파형에 대한 고속 푸리에 변환(FFT) 등과 같은 푸리에 변환을 통해 실험적으로 구한 값을 사용할 수 있다. 이때, 차수 n은 계측치를 적절히 표현할 수 있도록 충분히 크게 잡는다. 또한 광변조기(103)의 이득 파형의 평균값은 아래의 수학식(6)과 같이 정의될 수 있다.
Figure pat00012
여기서, T는 앞서 설명한 바와 같이 한 서브 프레임의 주기이다. 도 3은 수학식(5)로 표현된 광변조기(103)의 일반적인 이득 파형을 예시하는 것으로, 예를 들어 Te=25ns, n=5차항을 사용한 경우이고 이득 파형의 평균값은 0.4(W)이다.
2-2. 광변조기에 의한 광변조 및 광검출기에서의 세기 영상 생성
수학식(3)으로 표현된 반사광은 수학식(5)로 표현된 광변조기(103)의 이득 파형과 곱해져서 광검출기(105)에 도달한다. 상기 광변조기(103)를 투과한 후의 변조된 광의 광량 I(s)(t)(s=1,2,3)는 아래의 수학식(7)과 같은 다항식으로 나타낼 수 있다.
Figure pat00013
광변조기(103)를 투과한 광 영상은 광검출기(105)에서 소정의 노출 시간(T) 동안 누적되며, 이로부터 세기 영상이 형성된다. 3개의 세기 영상을 이용하여 거리 정보를 추출하는 경우에, 노출 시간 T는 거리 정보 영상(depth image)의 출력 주기의 1/3로 할 수 있다. 예를 들어, 초당 30프레임의 거리 정보 영상을 출력하는 경우(즉, fd=30, fd는 거리 정보 영상의 프레임 속도), 노출 시간 T는 약 0.011초가 된다. 광검출기(105)의 감도에 따라, 광검출기(105)에 도달하는 광량과 세기 영상의 크기 사이에 소정의 변환 비율이 존재할 수 있지만, 편의상 이를 단순화하여 세기 영상의 크기를 아래의 수학식(8)과 같이 정의할 수 있다.
Figure pat00014
상기 수학식(8)에서 세기 영상을 정의할 때, 광검출기(105)로의 입력 광 세기와 광검출기(105)로부터의 출력 세기 사이의 스케일 계수(scale factor)를 1로 가정하였다. 따라서, 실제 계산시에는 광검출기(105)에 대한 보정(calibration)을 통해 스케일 계수를 결정하여 수학식(8)에 적용할 필요가 있을 것이다.
일반적으로, 투사광의 주기 Te(예컨대, 25ns)는 광검출기(105)의 노출 시간 T(예컨대, 0.011s)에 비해서 매우 작다. 따라서, 수학식(8)의 적분시에, 수학식(7)에 존재하는 직류 항(DC term)이 주요한 기여를 하며, 교류 항(AC term)은 서로 음수/양수로 상쇄되어 그 적분 결과에 기여하는 양이 매우 적다. 따라서, 수학식(7)을 수학식(8)에 대입하여 그 직류 항만을 고려하면, 아래와 같은 세기 영상에 관한 관계식(9)를 얻을 수 있다. 이에 대한 상세한 유도 과정은 생략한다.
Figure pat00015
여기서,
Figure pat00016
이며,
Figure pat00017
이다.
수학식(10)에서 차수 l은 차수 m과 차수 n중 작은 값이다. F(s)TOF)는 미지수인 위상 지연 ΦTOF의 함수로서 상기의 수학식(10)과 같이 수학적인 표현이 가능한 항이다. 따라서, 수학식(9) 내지 수학식(11)로부터 알 수 있듯이, 세기 영상의 크기는 미지수인 위상 지연 ΦTOF에 관한 함수 F(s)TOF)와 위상 지연과 관련 없는 상수항
Figure pat00018
의 합으로서 정리될 수 있다.
2-3. 거리 정보 추출을 위한 문제의 정의 및 추출 방법
거리 정보를 추출하기 위한 문제는 아래의 식(12)과 같이 정의될 수 있다.
[수학식 12]
ㆍ 풀이할 대상:
Figure pat00019
ㆍ 측정으로부터 구할 수 있는 기지의 파라미터
Figure pat00020
: 광출력을 측정하여 구함
Figure pat00021
: 광변조기(103)의 이득을 측정하여 구함
Figure pat00022
: 광검출기(105)의 검출 결과
위의 문제를 풀기 위해서는 3개 이상의 방정식이 필요하다. 이를 위해서, 수학식(9)로 표현된 세기 영상에 관한 관계식으로부터 3개의 방정식을 구성할 수 있다. 즉, 3개의 투사광에 대응하는 3개의 세기 영상을 광검출기(105)로 순차적으로 검출하고, 수학식(9)를 이용하면 다음과 같은 연립 방정식을 구성할 수 있다.
Figure pat00023
Figure pat00024
Figure pat00025
따라서, 상기 수학식(12)의 문제 정의에서 언급된 TOF 위상 지연(ΦTOF), 반사도(r) 및 외광 크기(
Figure pat00026
)의 3개 미지수를 위의 수학식(13) 내지 수학식(15)의 3개의 방정식으로 구성된 연립방정식을 이용하여 추출할 수 있다.
그런데, 수학식(13) 내지 수학식(15)로 구성된 연립방정식은 미지수에 대해서 비선형 함수이므로, 수치해석을 통해 직접 미지수를 추출하는 것은 많은 계산 시간이 소요된다. 실시간으로 거리 정보 영상을 추출하기 위해서는, 비교적 짧은 계산 시간으로 정확한 결과를 도출할 수 있는 다른 방법이 필요하다. 따라서, 수학식(13) 내지 수학식(15)의 연립방정식을 푸는데 있어서 벡터 공간의 해법을 활용하여 룩업 테이블(Look-up Table)을 구성하는 방법을 제안한다.
이를 위하여, 먼저 아래의 수학식(16) 및 수학식(17)과 같은 3차원 공간 상의 벡터를 정의하면, 이를 기초로 룩업 테이블을 구성할 수 있다.
Figure pat00027
Figure pat00028
여기서,
Figure pat00029
Figure pat00030
수학식(16) 및 수학식(18)을 통해 알 수 있듯이, 세기 영상에 관한 3차원 공간 상의 벡터는 3개의 투사광과 각각 관련된 세기 영상들의 크기들 사이의 차들의 조합으로 구성될 수 있다. 마찬가지로, 수학식(17) 및 수학식(19)를 통해 알 수 있듯이, 위상 지연에 관한 3차원 공간 상의 벡터 함수는 3개의 투사광과 각각 관련된 위상 지연에 관한 함수들 사이의 차들의 조합으로 구성될 수 있다. 차들의 조합 대신에 차들의 비율들의 조합을 사용하는 것도 역시 가능하다. 도 8은 위상 지연에 관한 벡터 함수
Figure pat00031
와 세기 영상에 관한 세기 영상 벡터
Figure pat00032
을 3차원 공간 상에 표시한 예를 도시하고 있다.
Figure pat00033
는 수학적으로 알고 있는 벡터 함수로서 위상 지연 ΦTOF을 변수로서 갖는다. 세기 영상 벡터
Figure pat00034
는 측정을 통해서 결정되는 벡터이다. 한편, 수학식(16) 내지 수학식(19)로 정의된 벡터는 위수(rank)가 2이며, 도 8에 도시된 바와 같이, 법선 벡터가 (1,1,1)인 디스크를 형성한다. 따라서, 세기 영상 벡터
Figure pat00035
가 기준 벡터, 예를 들어
Figure pat00036
로부터 기울어진 회전 각도 α는 위상 지연 ΦTOF와 일대일 대응 관계를 가지고 있다. 따라서, 각도 α와 그에 대응하는 위상 지연 ΦTOF을 수치적으로 미리 계산하면, 그 관계를 표(table) 또는 아래의 수학식(20)과 같은 룩업 테이블 함수 LUT(α)의 형태로 정리할 수 있다. 이것이 본 발명에서 활용하는 룩업 테이블이다.
Figure pat00037
도 9는 이렇게 얻은 룩업 테이블을 예시적으로 도시하는 그래프이다. 도 9의 그래프는 각도 α와 위상 지연 ΦTOF 사이의 관계가 표로 정리되어 있는 것을 그래프의 형태로 그린 것이다. 도 9로부터 알 수 있듯이, 광검출기(105)의 측정을 통해 세기 영상 벡터
Figure pat00038
를 구할 수 있으면, 세기 영상 벡터
Figure pat00039
가 기준 벡터로부터 기울어진 각도 α를 구할 수 있으며, 상기 각도 α를 이용하여 룩업 테이블을 참조함으로써 위상 지연 ΦTOF을 알 수 있게 된다.
한편, 수학식(13) 내지 수학식(15)의 연립방정식에서
Figure pat00040
를 소거하면, 다음과 같은 해(solution)를 위한 조건을 얻을 수 있다. 이에 대한 자세한 유도 과정은 생략한다.
Figure pat00041
상기 수학식(21)은, 전술한 수학식(16)과 수학식(17)에서 정의된 두 벡터가 아래의 수학식(22)와 같이 3차원 공간 상에서 평행이기 위한 조건이다.
Figure pat00042
즉, 미지수 ΦTOF에 관한 조건은, 해석치
Figure pat00043
와 측정치
Figure pat00044
가 도 8에 도시된 3차원 공간의 디스크 상에서, 기준 벡터로부터 기울어진 각도 α 값이 동일하다는 조건이다. 따라서, 측정된
Figure pat00045
가 기준 벡터로부터 기울어진 각도 α에 해당하는 위상 지연 ΦTOF이 해가 된다. 이 해는 α를 먼저 계산한 다음, 룩업 테이블에서 해당하는 ΦTOF를 읽어서 찾을 수 있다. 각도 α를 구하는 방법은 여러 가지가 있을 수 있다. 예를 들어, 기준 벡터
Figure pat00046
로부터 기울어진 각도 α는 다음의 수학식(23)과 같이 계산될 수 있다.
Figure pat00047
여기서,
Figure pat00048
Figure pat00049
수학식(23)으로 계산된 각도 α와 미리 계산된 룩업 테이블을 이용해 위상 지연 ΦTOF를 구한다. 예를 들어, 도 9의 그래프에서 α=250도인 경우에 해당되는 위상 지연은 ΦTOF=270도이다. 이렇게 구한 위상 지연을 이용하여 광검출기(105)의 해당 화소에 대한 광학 장치(100)와 피사체(200) 사이의 거리(즉, depth) d를 아래의 수학식(26)으로 구할 수 있다.
Figure pat00050
여기서, c는 빛의 속도이다.
또한, 반사도 r은 수학식(21)로부터 아래와 같이 구할 수 있다.
Figure pat00051
또한, 외광 크기
Figure pat00052
는 수학식(11) 및 수학식(13)을 이용하여 아래와 같이 구할 수 있다.
Figure pat00053
수학식(8)에서 세기 영상을 정의할 때, 광검출기(105)로의 입력 광 세기와 광검출기(105)로부터의 출력 세기 사이의 스케일 계수(scale factor)를 1로 가정하였으나, 실제로는 특정한 비율로 크기 변환(scaling)된다. 따라서, 반사도 r과 외광 크기
Figure pat00054
의 물리적 수치도 역시 위의 특정한 비율로 수학식(27)과 수학식(28)의 결과를 크기 변환한 값이어야 하지만, 여기서는 편의상 스케일 계수를 1로 가정하여 그러한 크기 변환을 하지 않았다.
3. 흐름도(Flow Chart)
도 10은 상술한 거리 정보 추출 알고리즘에서 룩업 테이블을 구성하는 과정만을 요약하여 설명하는 흐름도를 도시하고 있다. 도 10을 참조하면, 단계(S10)에서, 광원(101)으로부터 방출될 적어도 3개의 상이한 투사광의 파형을 측정한다. 그러면 그 측정된 결과로부터 수학식(1)과 같이 투사광을 나타내는 일반적인 다항식의 계수들을 결정할 수 있다. 예를 들어, 수학식(1)에서 투사광 (s)에서의 k차 사인항의 계수 ak (s)와 투사광 (s)에서의 k차 코사인항의 계수 bk (s)를 결정할 수 있다. 만약 3개의 투사광을 사용하는 경우에, 3개의 투사광 각각에 대한 계수들을 측정을 통해 차례로 결정할 수 있다. 또한, 투사광의 평균값
Figure pat00055
도 역시 결정한다. 앞서 설명한 바와 같이, 3개의 투사광에 대한 평균값들이 같도록 제어될 수 있다.
다음으로 단계(S11)에서, 광변조기(103)의 이득 파형을 측정한다. 그러면 그 측정된 결과로부터 수학식(5)와 같이 이득 파형을 나타내는 일반적인 다항식의 계수들을 결정할 수 있다. 예를 들어, 수학식(5)에서 사인항과 코사인항의 각 차수의 계수 ck와 dk를 상기 측정을 통해 결정할 수 있다. 또한, 이득 파형의 평균값
Figure pat00056
도 역시 결정할 수 있다. 만약, 3개의 상이한 투사광과 1개의 이득 파형을 사용하는 대신, 1개의 투사광과 3개의 상이한 이득 파형을 사용하는 경우에는, 3개의 상이한 이득 파형을 모두 측정하여 그 각각에 대한 다항식의 계수들을 결정할 수 있다. 이하에서는 편의상 3개의 상이한 투사광과 1개의 이득 파형을 사용하는 경우에 대해 설명한다.
이어서 단계(S12)에서, 3개의 투사광이 피사체(200)에서 반사된 후의 반사광에 관한 일반적인 다항식을 수학식(3)과 같이 구한다. 그러면 수학식(3)으로부터, 상기 반사광이 광변조기(103)에서 변조된 후의 광량 I(s)(t)(s=1,2,3)을 수학식(7)과 같이 구할 수 있다. 즉, 수학식(7)은 반사광에 관한 수학식(3)과 이득 파형에 관한 수학식(5)를 곱하여 얻어질 수 있다. 그런 후, 변조된 광에 관한 수학식(7)을 수학식(8)과 같이 광검출기(105)에서의 노출 시간 T로 적분하면, 광변조기(103)에서 변조된 반사광을 광검출기(105)에서 소정의 노출 시간 T 동안 누적하여 얻은 세기 영상(또는 세기 영상의 크기)에 관한 다항식을 유도할 수 있다. 여기서, 3개의 투사광의 각각에 대해 세기 영상이 얻어지므로, 3개의 세기 영상에 관하여 다항식이 각각 존재하게 된다.
단계(S13)에서, 앞서 설명한 바와 같이, 수학식(7)에 존재하는 직류 항만이 적분 결과에 주요한 기여를 하고 교류 항은 기여가 적기 때문에, 직류 항만을 고려하여 세기 영상의 크기에 대한 수학식(7) 및 수학식(8)의 다항식을 정리할 수가 있다. 그러면, 수학식(9)와 같이 미지수인 위상 지연 ΦTOF에 관한 함수 F(s)TOF)(수학식(10) 참조)와 위상 지연과 관련 없는 상수항
Figure pat00057
(수학식(11) 참조)의 합으로서 상기 세기 영상의 크기에 대한 다항식이 정리될 수 있다. 수학식(11)과 같이 상수항
Figure pat00058
는 투사광의 평균값, 이득 파형의 평균값 및 투과광/반사광 이외의 성분인 외광 성분으로 구성된다.
이렇게 정리된 세기 영상의 크기에 관한 다항식과 위상 지연에 관한 함수는 수학식(9)와 수학식(10)으로 표현된 바와 같이, 3개의 투사광에 대해 각각 하나씩 존재한다. 단계(S14)에서, 수학식(17) 및 수학식(19)로 표현되는 바와 같이, 위상 지연에 관한 3개의 함수 F(1)TOF), F(2)TOF), F(3)TOF)들 사이의 차들의 조합을 이용하여 3차원 공간 상에 위상 지연에 관한 벡터 함수
Figure pat00059
를 정의할 수 있다. 상술한 바와 같이, 상기 벡터 함수는 법선 벡터가 (1,1,1)인 디스크를 형성할 수 있다(도 8 참조). 그리고, 이러한 3차원 공간 내에서 좌표 기준을 설정하기 위하여, 특정한 값의 벡터를 기준 벡터로서 정의할 수 있다. 예를 들어, 위상 지연에 관한 벡터 함수 중에서 위상 지연 ΦTOF=0인 벡터
Figure pat00060
를 x축으로 설정할 경우, 기준 벡터는 수학식(24) 및 수학식(25)와 같이 정의될 수 있다.
단계(S15)에서, 세기 영상의 크기에 관한 3개의 정리된 수학식(13) 내지 수학식(15)을 이용하여, 수학식(16) 및 수학식(18)로 표현되는 바와 같이, 3개의 투사광과 각각 관련된 세기 영상들의 크기들 사이의 차들의 조합으로 3차원 공간 상에 세기 영상 벡터
Figure pat00061
를 정의할 수 있다. 상술한 바와 같이, 상기 세기 영상 벡터는 측정을 통해 결정되는 벡터로서, 위상 지연에 관한 벡터 함수와 마찬가지로 3차원 공간 내에서 법선 벡터가 (1,1,1)인 디스크를 형성할 수 있다. 이러한 세기 영상 벡터가 기준 벡터로부터 기울어진 각도 α는 위상 지연 ΦTOF와 일대일 대응 관계를 가지고 있다. 따라서, 단계(S16)에서 세기 영상 벡터가 기준 벡터와 이루는 각도 α와 그에 대응하는 위상 지연 사이의 관계를, 예를 들어, 수치해석과 같은 해석적인 방식으로 미리 계산하여 결정해 둘 수 있다. 그러면 단계(S17)에서, 수학식(20)과 같이 상기 결정된 관계를 나타내는 룩업 테이블을 표 또는 함수의 형태로 구성하는 것이 가능하게 된다.
이렇게 구성된 룩업 테이블은 3개의 특정 투사광들과 1개의 특정 이득 파형의 조합에 관하여 미리 측정된 결과로 얻은 것이다. 만약 3개의 투사광들과 1개의 이득 파형 중에서 어느 하나의 위상 또는 파형 종류가 바뀐다면 그에 대응하는 새로운 룩업 테이블을 구성하여 둘 수 있다. 예를 들어, 광학 장치(100)에서 사용할 수 있는 투사광과 이득 파형의 다양한 조합별로 각각 별도의 룩업 테이블을 구성할 수 있다. 그러면, 광학 장치(100)의 동작시, 사용된 투사광과 이득 파형의 특정한 조합에 따라 그에 해당하는 룩업 테이블을 참조하는 것이 가능하다. 이러한 투사광과 이득 파형의 다양한 조합은 3개의 투사광들과 1개의 이득 파형의 조합뿐만 아니라, 4개 이상의 투사광과 1개의 이득 파형의 조합, 1개의 투사광과 3개 이상의 이득 파형의 조합도 포함할 수 있다. 어느 경우이든, 상술한 룩업 테이블의 구성 방법은 동일하게 적용될 수 있다. 이러한 룩업 테이블은, 예를 들어, 도 1에 도시된 거리정보영상 처리기(107) 내에 함께 통합되어 있을 수도 있으며, 또는 도시되지 않은 별도의 저장 매체에 저장되어 있을 수 있다.
도 11은 상술한 방식으로 얻은 룩업 테이블을 이용하여 실제로 거리 정보를 추출하는 과정을 요약하여 설명하는 흐름도를 도시하고 있다. 다른 방식의 투사광과 이득 파형의 조합도 가능하지만, 여기서도 편의상 3개의 상이한 투사광과 1개의 이득 파형을 사용하는 경우에 대해 예시적으로 설명한다.
먼저, 단계(S20)에서, 광원(101)이 3개의 상이한 투사광을 시간 순서에 따라 차례로 발생시켜 피사체(200)에 투사한다. 앞서 설명한 바와 같이, 3개의 투사광은 동일한 주기를 가지며, 단지 위상만이 서로 다를 수도 있고, 또는 파형 자체가 서로 다를 수도 있다. 상기 투사광은 예를 들어 약 850nm의 파장을 갖는 근적외선 광일 수 있다.
그런 후, 3개의 투사광은 피사체(200)에서 반사되어 차례로 광학 장치(100)로 되돌아온다. 이때, 필터(109)를 이용하여 반사광 중에서 예를 들어 약 850nm의 근적외선 광만을 투과시킨다. 그러면 단계(S21)에서, 광변조기(103)는 피사체(200)에서 반사된 3개의 반사광을 소정의 이득 파형으로 광변조하고, 변조된 광을 광검출기(105)에 제공한다. 광변조기(103)의 이득 파형은 어떠한 것이라도 사용될 수 있지만, 그 주기는 투사광의 주기와 동일하다.
단계(S22)에서, 광검출기(105)는 상기 광변조기(103)에서 변조된 반사광을 소정의 노출 시간 T 동안 누적하여 출력한다. 그러면, 3개의 투사광에 각각 대응하는 3개의 세기 영상을 얻을 수 있다. 이렇게 얻은 3개의 세기 영상은 거리정보영상 처리기(107)에 제공될 수 있다. 단계(S23)에서, 거리정보영상 처리기(107)는 3개의 세기 영상의 크기를 이용하여, 수학식(16) 및 수학식(18)로 설명된 바와 같이, 3차원 공간 상의 세기 영상 벡터를 구성한다.
다음 단계(S24)에서, 거리정보영상 처리기(107)는 수학식(23)을 이용하여 세기 영상 벡터가 3차원 공간 상의 기준 벡터로부터 기울어진 각도 α를 구한다. 상술한 바와 같이, 상기 기준 벡터는 수학식(24) 및 수학식(25)로부터 미리 정의되어 있다. 각도 α를 구한 후에, 거리정보영상 처리기(107)는 미리 저장되어 있는 룩업 테이블을 참조하여 각도 α에 대응하는 위상 지연을 구할 수 있다. 여기서 룩업 테이블은 사용 가능한 투사광과 이득 파형의 다양한 조합별로 다양하게 마련되어 있을 수 있다. 따라서, 거리정보영상 처리기(107)는, 광원(101)과 광변조기(103)에서 사용된 투사광과 이득 파형의 특정한 조합에 따라 그에 해당하는 룩업 테이블을 참조할 수 있다.
마지막으로 단계(S25)에서, 거리정보영상 처리기(107)는 상술한 방식으로 구한 위상 지연을 이용하여 피사체로부터의 거리, 피사체의 반사도 및 외광의 크기를 계산할 수 있다. 예를 들어, 수학식(26) 내지 수학식(28)에 기재된 공식에 위상 지연을 대입함으로써 피사체로부터의 거리, 피사체의 반사도 및 외광의 크기를 계산하는 것이 가능하다.
지금까지 광검출기(105)의 어느 한 화소에 대해 거리 정보를 추출하는 방식을 설명하였지만, 광검출기(105) 내의 다수의 화소들의 각각에 대해 동일한 방식이 적용될 수 있다. 또한, CCD와 같이 2차원 어레이로 배열된 다수의 화소를 갖는 광검출기(105)뿐만 아니라, 1차원 어레이로 배열된 다수의 화소를 갖는 광검출기 또는 단지 하나의 포토다이오드만으로 구성된 광검출기에 대해서도 위에서 설명한 방법이 그대로 적용될 수 있다.
도 10 및 도 11에 도시된 작업을 수행하는 데 필요한 연산량을 살펴보면, 피사체까지의 거리 정보를 추출하는데 필요한 계산의 대부분은 도 10의 룩업 테이블 구성 과정에서 이루어진다. 룩업 테이블은, 예를 들어, 광학 장치(100)의 제조시에 제조업자가 미리 계산하여 만들어 놓은 것일 수 있다. 따라서, 거리 정보를 추출하는데 있어서 필요한 실시간 연산은, 도 11에 도시된 세기 영상 벡터를 구성하기 위한 계산, 룩업 테이블을 참조하여 위상 지연을 구하기 위한 계산, 및 위상 지연으로부터 피사체로부터의 거리, 피사체의 반사도 및 외광의 크기 정보를 추출하기 위한 계산 정도이다. 이러한 점에서 광학 장치(100)가 직접 수행하는 연산량이 비교적 적기 때문에, 광학 장치(100)를 이용하여 실시간으로 거리 정보를 추출하는 것이 가능하다. 따라서, 다수의 화소를 갖는 화상 정보와 함께 화상 내의 각 화소별 거리 정보를 동시에 획득할 수 있는 3D 카메라의 구현이 가능하다.
4. 실험 결과
상술한 방식의 거리 정보 추출 방법을 확인하기 위하여 시뮬레이션을 수행하였다. 투사광으로는 도 2에 도시된 25ns의 주기를 갖는 3종의 파형(실선)을 사용하였으며, 이득 파형으로는 도 3에 도시된 25ns의 주기를 갖는 파형(실선)을 사용하였다. 0.5~3.5m의 거리에 있는 피사체에 대해 상술한 방식으로 거리를 측정하였으며, 광검출기(예컨대, CCD)의 노이즈를 고려하기 위해 세기 영상에 그 크기의 0.1% 정도인 랜덤 노이즈(random noise)를 의도적으로 인가하였다.
도 12는 시뮬레이션을 통해 거리 정보를 추출한 결과를 참값과 비교한 그래프이며, 도 13은 도 12의 거리 정보 추출 결과와 참값과의 오차를 나타낸 그래프이다. 도 12 및 도 13을 참조하면, 약 3m 거리 내에서 2cm 가량의 오차를 나타낸다는 것을 알 수 있다. 이는 광검출기에 의도적으로 인가한 노이즈에 기인한다. 거리가 멀수록 반사광의 크기가 작아져 신호대잡음비가 작아지기 때문에, 거리가 멀수록 오차가 커질 수 있다. 3m 내에서 약 2cm 이내의 오차는 상용의 거리 정보 카메라로서 사용 가능한 수준이라고 볼 수 있다.
또한, 상술한 알고리즘을 수행하기 위한 하드웨어를 제작하여 실험을 하였다. 이 실험에서, 레이저 다이오드(LD)와 광변조기를 약 20MHz로 구동하여 20cm 간격으로 20cm~180cm의 거리를 측정하였다. 도 14는 레이저 다이오드의 광출력을 측정한 결과이고, 도 15는 레이저 다이오드의 광출력 스펙트럼을 고속 푸리에 변환으로 분석한 결과이다. 도 15를 보면, 레이저 다이오드의 광출력에서 20MHz의 기본 주파수 이외에 고차항이 존재함을 알 수 있다. 이러한 고차항을 상술한 거리 정보 추출 알고리즘에 반영하였다. 또한, 도 16은 광변조기의 변조 이득 파형을 측정한 결과로서, DC-오프셋 및 상당량의 파형 왜곡을 관찰할 수 있다. 도 17은 상기 이득 파형 스펙트럼을 고속 푸리에 변환을 통해 분석한 결과이다. 도 17을 보면, 광변조기의 변조 이득 파형은 20MHz의 기본 주파수 이외에 고차항이 상당히 존재함을 알 수 있다. 이러한 고차항을 상술한 거리 정보 추출 알고리즘에 반영하였다.
도 18은 상술한 알고리즘에 따라 실제로 거리 정보를 추출한 결과이다. 도 18에는 상술한 알고리즘에 따라 파형의 왜곡을 실측하여 고려한 결과와 이상적인 파형만을 고려한 결과가 비교되어 있다. 도 18에서 실선으로 표시한 것은 실측된 파형 왜곡을 고려한 경우이고 점선으로 표시한 것은 파형 왜곡을 고려하지 않는 경우이다. 도 18의 결과를 살펴보면, 파형 왜곡을 실측하여 고려한 경우에 전체적으로 약 0.7cm의 평균적인 오차가 있었지만, 파형 왜곡을 고려하지 않은 경우에는 전체적으로 약 2.1cm의 평균적인 오차가 발생하였다.
상술한 바와 같이, 개시된 알고리즘은 수학적으로 이상적인 파형이 아닌 실제의 실측된 비선형적 파형을 고려하여 거리 정보를 추출하기 때문에, 사용 파형에 제한을 받지 않고 정확한 거리 정보를 추출할 수 있다. 따라서, 왜곡 및 비선형성이 거의 없는 고가의 광원이나 광변조 소자를 사용할 필요가 없으며, 복잡한 오차 보상 수단이 요구되지 않는다. 또한, 기존의 광원, 광변조 소자 및 광학 장치를 그대로 이용할 수 있어서 추가 비용이 소요되지 않는다. 더욱이, 미리 계산된 거리 정보가 저장되어 있는 룩업 테이블을 사용하기 때문에, 거리 정보를 추출하는데 있어서 연산량이 매우 작으므로 실시간 거리 정보 영상의 촬영이 가능하다. 따라서, 상술한 알고리즘이 채용된 광학 장치(100)는, 예를 들어 3D 카메라나 레이저 레이더(Laser Radar; LADAR)와 같은 분야에서 유용하게 사용될 수 있다.
지금까지, 본 발명의 이해를 돕기 위하여 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치에 대한 예시적인 실시예가 설명되고 첨부된 도면에 도시되었다. 그러나, 이러한 실시예는 단지 본 발명을 예시하기 위한 것이고 이를 제한하지 않는다는 점이 이해되어야 할 것이다. 그리고 본 발명은 도시되고 설명된 설명에 국한되지 않는다는 점이 이해되어야 할 것이다. 이는 다양한 다른 변형이 본 기술분야에서 통상의 지식을 가진 자에게 일어날 수 있기 때문이다.
100....광학 장치 101....광원
102....광원 구동부 103....광변조기
104....광변조기 구동부 105....광검출기
106....제어부 107....거리정보영상 처리기
108....제 1 렌즈 109....필터
110....제 2 렌즈 200....피사체

Claims (37)

  1. 적어도 3개의 상이한 투사광을 순차적으로 피사체에 각각 투사하는 단계;
    피사체에서 반사된 반사광을 이득 파형을 갖는 광변조기로 변조하는 단계;
    상기 변조된 반사광을 광검출기에서 검출하여 상기 적어도 3개의 상이한 투사광에 각각 대응하는 적어도 3개의 세기 영상을 얻는 단계;
    상기 적어도 3개의 세기 영상을 이용하여 반사광의 위상 지연을 구하는 단계; 및
    상기 위상 지연으로부터 피사체까지의 거리를 계산하는 단계를 포함하며,
    상기 위상 지연을 구하는 단계는 상기 적어도 3개의 상이한 투사광과 이득 파형을 미리 실측하여 구성한 룩업 테이블을 참조하는 거리 정보 추출 방법.
  2. 제 1 항에 있어서,
    상기 적어도 3개의 상이한 투사광과 상기 광변조기의 이득 파형은 동일한 주기를 갖는 거리 정보 추출 방법.
  3. 제 2 항에 있어서,
    상기 적어도 3개의 상이한 투사광은 파형이 동일하고 위상이 서로 다르거나 또는 파형이 서로 다른 거리 정보 추출 방법.
  4. 제 2 항에 있어서,
    상기 적어도 3개의 상이한 투사광은 그 평균값이 서로 같도록 제어되는 거리 정보 추출 방법.
  5. 제 1 항에 있어서,
    상기 광검출기는 상기 적어도 3개의 상이한 투사광과 각각 대응하는 적어도 3개의 변조된 반사광을 소정의 노출 시간 동안 각각 누적하여 상기 적어도 3개의 상이한 투사광에 각각 대응하는 적어도 3개의 세기 영상을 형성하는 거리 정보 추출 방법.
  6. 제 1 항에 있어서,
    상기 위상 지연으로부터 피사체의 반사도 및 외광의 크기를 구하는 단계를 더 포함하는 거리 정보 추출 방법.
  7. 제 1 항에 있어서,
    상기 위상 지연을 구하는 단계는:
    상기 적어도 3개의 세기 영상을 이용하여 3차원 공간 상의 세기 영상 벡터를 구성하는 단계;
    상기 세기 영상 벡터가 3차원 공간 상의 기준 벡터로부터 기울어진 각도를 구하는 단계; 및
    상기 기울어진 각도로 상기 룩업 테이블을 참조하여 반사광의 위상 지연을 구하는 단계를 포함하는 거리 정보 추출 방법.
  8. 제 7 항에 있어서,
    상기 룩업 테이블은 상기 기울어진 각도와 반사광의 위상 지연 사이의 미리 결정된 관계를 포함하는 거리 정보 추출 방법.
  9. 제 7 항에 있어서,
    상기 세기 영상 벡터는 상기 적어도 3개의 세기 영상들의 크기들 사이의 차들의 조합으로 구성되는 거리 정보 추출 방법.
  10. 제 1 항에 있어서,
    상기 적어도 3개의 투사광과 이득 파형의 상이한 조합별로 각각 별도의 룩업 테이블이 존재하며, 상기 적어도 3개의 투사광과 이득 파형의 특정한 조합에 따라 그에 해당하는 룩업 테이블을 참조하는 거리 정보 추출 방법.
  11. 제 1 항에 있어서,
    상기 룩업 테이블은:
    상기 적어도 3개의 상이한 투사광과 상기 광변조기의 이득 파형을 측정하여 각각의 투사광과 이득 파형을 나타내는 다항식을 각각 구하는 단계;
    상기 적어도 3개의 상이한 투사광과 상기 광변조기의 이득 파형에 관한 다항식을 이용하여 적어도 3개의 세기 영상에 관한 다항식을 구하는 단계;
    상기 적어도 3개의 세기 영상에 관한 다항식을 각각 위상 지연에 관한 함수와 위상 지연과 관련 없는 상수항의 합으로 정리하는 단계;
    상기 위상 지연에 관한 적어도 3개의 함수들 사이의 차들의 조합을 이용하여 3차원 공간 상에 위상 지연에 관한 벡터 함수를 정의하고, 특정한 벡터를 기준 벡터로서 정의하는 단계;
    상기 적어도 3개의 세기 영상에 관한 다항식들 사이의 차들의 조합으로 3차원 공간 상에 세기 영상 벡터를 정의하는 단계;
    세기 영상 벡터가 기준 벡터로부터 기울어진 각도와 상기 각도에 대응하는 위상 지연 사이의 관계를 해석적으로 결정하는 단계; 및
    상기 해석적으로 결정된 관계를 나타내는 룩업 테이블을 구성하는 단계를 포함하는 방식으로 구성되는 거리 정보 추출 방법.
  12. 제 11 항에 있어서,
    상기 세기 영상에 관한 다항식을 구하는 단계는:
    상기 적어도 3개의 상이한 투사광이 피사체에서 반사된 후의 반사광에 관한 일반적인 다항식을 상기 적어도 3개의 상이한 투사광에 관한 다항식으로부터 유도하는 단계;
    상기 반사광에 관한 다항식과 상기 이득 파형에 관한 다항식을 곱하여 광변조된 광에 대한 다항식을 구하는 단계; 및
    상기 광변조된 광에 대한 다항식을 광검출기에서의 노출 시간으로 적분하여 상기 적어도 3개의 상이한 투사광에 각각 대응하는 적어도 3개의 세기 영상에 관한 다항식을 구하는 단계를 포함하는 거리 정보 추출 방법.
  13. 제 11 항에 있어서,
    상기 적어도 3개의 상이한 투사광과 상기 광변조기의 이득 파형을 측정하여 각각의 투사광과 이득 파형을 나타내는 다항식을 각각 구하는 단계는:
    상기 적어도 3개의 상이한 투사광을 각각 측정하는 단계;
    상기 각각의 투사광을 측정한 결과를 시간 축상 직교하는 함수로 나타내는 단계;
    상기 광변조기의 이득 파형을 측정하는 단계; 및
    상기 이득 파형을 측정한 결과를 시간 축상 직교하는 함수로 나타내는 단계를 포함하는 거리 정보 추출 방법.
  14. 제 13 항에 있어서,
    상기 각각의 투사광을 측정한 결과를 시간 축상 직교하는 함수로 나타내는 단계는, 푸리에 급수의 사인항의 계수와 코사인항의 계수 또는 체비셰프 다항식의 계수를 결정함으로써 각각의 투사광에 대한 다항식을 구하는 단계를 포함하는 거리 정보 추출 방법.
  15. 제 13 항에 있어서,
    상기 이득 파형을 측정한 결과를 시간 축상 직교하는 함수로 나타내는 단계는, 푸리에 급수의 사인항의 계수와 코사인항의 계수 또는 체비셰프 다항식의 계수를 결정함으로써 이득 파형에 대한 다항식을 구하는 단계를 포함하는 거리 정보 추출 방법.
  16. 제 13 항에 있어서,
    상기 각각의 투사광과 이득 파형을 측정한 결과로부터 각각의 투사광의 평균값과 이득 파형의 평균값을 결정하는 단계를 더 포함하는 거리 정보 추출 방법.
  17. 제 11 항에 있어서,
    상기 적어도 3개의 세기 영상에 관한 다항식을 각각 위상 지연에 관한 함수와 위상 지연과 관련 없는 상수항의 합으로 정리하는 단계는 상기 적어도 3개의 세기 영상에 관한 다항식에 존재하는 직류 항만을 고려하는 거리 정보 추출 방법.
  18. 제 11 항에 있어서,
    상기 차원 공간 상의 위상 지연에 관한 벡터 함수는
    Figure pat00062
    와 같이 정의되며, 여기서
    Figure pat00063
    이고(i,j=1,2,3),
    F(i)TOF)는 i번째 투사광에 대한 위상 지연에 관한 함수를 나타내는 거리 정보 추출 방법.
  19. 제 18 항에 있어서,
    상기 위상 지연에 관한 벡터 함수 중에서 위상 지연 ΦTOF=0인 벡터
    Figure pat00064
    를 기준 벡터로 정의하는 거리 정보 추출 방법.
  20. 제 11 항에 있어서,
    상기 3차원 공간 상의 세기 영상 벡터는,
    Figure pat00065
    와 같이 정의되며, 여기서
    Figure pat00066
    이고(i,j=1,2,3),
    I(i) CCD는 i번째 투사광에 대한 세기 영상을 나타내는 거리 정보 추출 방법.
  21. 제 11 항에 있어서,
    상기 적어도 3개의 투사광과 이득 파형의 상이한 조합별로 각각 별도의 룩업 테이블을 구성하는 거리 정보 추출 방법.
  22. 하나의 주기적인 파형을 갖는 투사광을 피사체에 투사하는 단계;
    피사체에서 반사된 반사광을 적어도 3개의 상이한 이득 파형을 갖는 광변조기로 순차적으로 변조하는 단계;
    상기 변조된 반사광을 광검출기에서 검출하여 상기 적어도 3개의 상이한 이득 파형에 각각 대응하는 적어도 3개의 세기 영상을 얻는 단계;
    상기 적어도 3개의 세기 영상을 이용하여 반사광의 위상 지연을 구하는 단계; 및
    상기 위상 지연으로부터 피사체까지의 거리를 계산하는 단계를 포함하며,
    상기 위상 지연을 구하는 단계는 상기 투사광과 상기 적어도 3개의 상이한 이득 파형을 미리 실측하여 구성한 룩업 테이블을 참조하는 거리 정보 추출 방법.
  23. 제 22 항에 있어서,
    상기 투사광과 상기 광변조기의 적어도 3개의 상이한 이득 파형은 동일한 주기를 갖는 거리 정보 추출 방법.
  24. 제 22 항에 있어서,
    상기 광검출기는 상기 적어도 3개의 상이한 이득 파형과 각각 대응하는 적어도 3개의 변조된 반사광을 소정의 노출 시간 동안 각각 누적하여 상기 적어도 3개의 상이한 이득 파형에 각각 대응하는 적어도 3개의 세기 영상을 형성하는 거리 정보 추출 방법.
  25. 제 22 항에 있어서,
    상기 위상 지연을 구하는 단계는:
    상기 적어도 3개의 세기 영상을 이용하여 3차원 공간 상의 세기 영상 벡터를 구성하는 단계;
    상기 세기 영상 벡터가 3차원 공간 상의 기준 벡터로부터 기울어진 각도를 구하는 단계; 및
    상기 기울어진 각도로 상기 룩업 테이블을 참조하여 반사광의 위상 지연을 구하는 단계를 포함하는 거리 정보 추출 방법.
  26. 제 25 항에 있어서,
    상기 룩업 테이블은 상기 기울어진 각도와 반사광의 위상 지연 사이의 미리 결정된 관계를 포함하는 거리 정보 추출 방법.
  27. 제 25 항에 있어서,
    상기 세기 영상 벡터는 상기 적어도 3개의 세기 영상들의 크기들 사이의 차들의 조합으로 구성되는 거리 정보 추출 방법.
  28. 제 22 항에 있어서,
    상기 투사광과 상기 적어도 3개의 이득 파형의 상이한 조합별로 각각 별도의 룩업 테이블이 존재하며, 상기 투사광과 상기 적어도 3개의 이득 파형의 특정한 조합에 따라 그에 해당하는 룩업 테이블을 참조하는 거리 정보 추출 방법.
  29. 주기적인 투사광을 발생시키는 광원;
    피사체로부터 반사된 광을 주기적인 이득 파형으로 광변조하기 위한 광변조기;
    상기 광변조기에 의해 광변조된 영상을 검출하기 위한 광검출기; 및
    상기 광검출기의 출력을 기초로 피사체와의 거리를 계산하기 위한 거리정보영상 처리기를 포함하며,
    상기 거리정보영상 처리기는 제 1 항 내지 제 28 항 중 어느 한 항에 따른 방법으로 거리 정보를 추출하는 광학 장치.
  30. 제 29 항에 있어서,
    상기 광원을 구동시켜 투사광의 파형을 제어하는 광원 구동부;
    상기 광변조기를 구동시켜 이득 파형을 제어하는 광변조기 구동부; 및
    상기 광원 구동부, 광변조기 구동부 및 광검출기의 동작을 제어하기 위한 제어부를 더 포함하는 광학 장치.
  31. 제 29 항에 있어서,
    상기 광변조기의 광입사면에서 반사광을 광변조기의 영역 내에 집광하는 제 1 렌즈;
    상기 제 1 렌즈와 광변조기 사이에서 소정의 파장을 갖는 광만을 투과시키는 필터; 및
    상기 광변조기와 광검출기 사이에서 상기 광변조된 영상을 광검출기의 영역 내에 집광하는 제 2 렌즈를 더 포함하는 광학 장치.
  32. 제 29 항에 있어서,
    상기 광검출기는 2차원 또는 1차원 어레이를 갖는 CCD 이미지 센서, CMOS 이미지 센서, 포토다이오드 어레이를 사용하거나 또는 단일 지점의 거리 측정을 위한 1개의 포토다이오드를 사용하는 광학 장치.
  33. 제 29 항에 따른 광학 장치를 포함하는 3D 카메라.
  34. 제 29 항에 따른 광학 장치를 포함하는 레이저 레이더.
  35. 적어도 3개의 상이한 투사광 및 광변조기의 이득 파형을 측정하여 각각의 투사광과 이득 파형을 나타내는 다항식을 각각 구하는 단계;
    상기 적어도 3개의 상이한 투사광과 상기 광변조기의 이득 파형에 관한 다항식을 이용하여 적어도 3개의 세기 영상에 관한 다항식을 구하는 단계;
    상기 적어도 3개의 세기 영상에 관한 다항식을 각각 위상 지연에 관한 함수와 위상 지연과 관련 없는 상수항의 합으로 정리하는 단계;
    상기 위상 지연에 관한 적어도 3개의 함수들 사이의 차들의 조합을 이용하여 3차원 공간 상에 위상 지연에 관한 벡터 함수를 정의하고, 특정한 벡터를 기준 벡터로서 정의하는 단계;
    상기 적어도 3개의 세기 영상에 관한 다항식들 사이의 차들의 조합으로 3차원 공간 상에 세기 영상 벡터를 정의하는 단계;
    세기 영상 벡터가 기준 벡터로부터 기울어진 각도와 상기 각도에 대응하는 위상 지연 사이의 관계를 해석적으로 결정하는 단계; 및
    상기 해석적으로 결정된 관계를 표 또는 함수의 형태로 나타내는 단계를 포함하는 거리 정보 추출용 룩업 테이블을 구성하는 방법.
  36. 하나의 투사광 및 광변조기의 적어도 3개의 상이한 이득 파형을 측정하여 각각의 투사광과 이득 파형을 나타내는 다항식을 각각 구하는 단계;
    상기 투사광과 상기 광변조기의 적어도 3개의 상이한 이득 파형에 관한 다항식을 이용하여 적어도 3개의 세기 영상에 관한 다항식을 구하는 단계;
    상기 적어도 3개의 세기 영상에 관한 다항식을 각각 위상 지연에 관한 함수와 위상 지연과 관련 없는 상수항의 합으로 정리하는 단계;
    상기 위상 지연에 관한 적어도 3개의 함수들 사이의 차들의 조합을 이용하여 3차원 공간 상에 위상 지연에 관한 벡터 함수를 정의하고, 특정한 벡터를 기준 벡터로서 정의하는 단계;
    상기 적어도 3개의 세기 영상에 관한 다항식들 사이의 차들의 조합으로 3차원 공간 상에 세기 영상 벡터를 정의하는 단계;
    세기 영상 벡터가 기준 벡터로부터 기울어진 각도와 상기 각도에 대응하는 위상 지연 사이의 관계를 해석적으로 결정하는 단계; 및
    상기 해석적으로 결정된 관계를 표 또는 함수의 형태로 나타내는 단계를 포함하는 거리 정보 추출용 룩업 테이블을 구성하는 방법.
  37. 제 35 항 또는 제 36 항에 따른 방법으로 구성된 룩업 테이블을 저장하고 있는 저장 매체.
KR1020100005753A 2010-01-21 2010-01-21 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치 KR101675112B1 (ko)

Priority Applications (3)

Application Number Priority Date Filing Date Title
KR1020100005753A KR101675112B1 (ko) 2010-01-21 2010-01-21 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치
US12/837,814 US8611610B2 (en) 2010-01-21 2010-07-16 Method and apparatus for calculating a distance between an optical apparatus and an object
JP2011011423A JP5739174B2 (ja) 2010-01-21 2011-01-21 距離情報抽出方法及び該方法を採用した光学装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020100005753A KR101675112B1 (ko) 2010-01-21 2010-01-21 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치

Publications (2)

Publication Number Publication Date
KR20110085785A true KR20110085785A (ko) 2011-07-27
KR101675112B1 KR101675112B1 (ko) 2016-11-22

Family

ID=44277614

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020100005753A KR101675112B1 (ko) 2010-01-21 2010-01-21 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치

Country Status (3)

Country Link
US (1) US8611610B2 (ko)
JP (1) JP5739174B2 (ko)
KR (1) KR101675112B1 (ko)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130101176A1 (en) * 2011-10-25 2013-04-25 Samsung Electronic Co.,Ltd. 3d image acquisition apparatus and method of calculating depth information in the 3d image acquisition apparatus
KR20130091194A (ko) * 2012-02-07 2013-08-16 삼성전자주식회사 3차원 영상 획득 장치 및 상기 3차원 영상 획득 장치에서 깊이 정보를 추출하는 방법
KR20140121710A (ko) * 2013-04-08 2014-10-16 삼성전자주식회사 3차원 영상 획득 장치 및 3차원 영상 획득 장치에서의 깊이 영상 생성 방법
WO2015102212A1 (ko) * 2013-12-31 2015-07-09 고려대학교 산학협력단 레이저 거리 센서에 의해 측정된 측정 거리의 거리 유형을 평가하는 방법 및 이를 이용한 이동 로봇의 위치 추정 방법
US9310488B2 (en) 2010-09-17 2016-04-12 Samsung Electronics Co., Ltd. Apparatus and method for generating depth image
US9350925B2 (en) 2011-11-02 2016-05-24 Samsung Electronics Co., Ltd. Image processing apparatus and method
WO2016200166A1 (ko) * 2015-06-09 2016-12-15 (주) 솔 특성 평가를 통한 광센서 어레이 모듈 보정방법
KR20160144148A (ko) * 2015-06-08 2016-12-16 삼성전자주식회사 광학 장치 및 깊이 정보 생성 방법
KR20160144692A (ko) * 2015-06-09 2016-12-19 (주) 솔 광센서 어레이의 특성평가를 통한 보정방법
US9667944B2 (en) 2012-10-10 2017-05-30 Samsung Electronics Co., Ltd. Imaging optical system and 3D image acquisition apparatus including the imaging optical system
US10101453B2 (en) 2014-08-11 2018-10-16 Infineon Technologies Ag Time of flight apparatuses and an illumination source
US10205933B2 (en) 2014-08-29 2019-02-12 Samsung Electronics Co., Ltd. Depth image acquisition apparatus and method of acquiring depth information
US10430956B2 (en) 2016-09-23 2019-10-01 Samsung Electronics Co., Ltd. Time-of-flight (TOF) capturing apparatus and image processing method of reducing distortion of depth caused by multiple reflection
US10509125B2 (en) 2015-12-24 2019-12-17 Samsung Electronics Co., Ltd. Method and device for acquiring distance information
US10545237B2 (en) 2016-08-22 2020-01-28 Samsung Electronics Co., Ltd. Method and device for acquiring distance information
US10802119B2 (en) 2016-07-26 2020-10-13 Samsung Electronics Co., Ltd. Lidar device and method of measuring distance using the same
WO2021086031A1 (ko) * 2019-10-31 2021-05-06 부산대학교 산학협력단 이중 가변 기반의 발진 광주파수 스캐닝 레이저 광원 그리고 그를 이용하는 측정 장치 및 중심파장 별 전파 각도 전환을 이용하는 대상물의 각도에 따른 거리 측정 장치
KR102302424B1 (ko) * 2020-10-26 2021-09-16 (주)시그봇 오브젝트 인식을 위한 센서 장치
KR102353365B1 (ko) * 2020-08-31 2022-01-19 부산대학교 산학협력단 컬러천이 레이저의 중심파장 별 전파 각도 전환을 이용하는 대상물의 각도에 따른 거리 측정 장치

Families Citing this family (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8294809B2 (en) 2005-05-10 2012-10-23 Advanced Scientific Concepts, Inc. Dimensioning system
US9753128B2 (en) * 2010-07-23 2017-09-05 Heptagon Micro Optics Pte. Ltd. Multi-path compensation using multiple modulation frequencies in time of flight sensor
US9857469B2 (en) * 2010-10-22 2018-01-02 Heptagon Micro Optics Pte. Ltd. System and method for multi TOF camera operation using phase hopping
WO2013012335A1 (en) * 2011-07-21 2013-01-24 Ziv Attar Imaging device for motion detection of objects in a scene, and method for motion detection of objects in a scene
JP6047318B2 (ja) 2012-07-06 2016-12-21 矢崎総業株式会社 ドライバ状態検出装置及びドライバ状態報知装置
WO2014041350A1 (en) * 2012-09-13 2014-03-20 Mbda Uk Limited Room occupancy sensing apparatus and method
KR102056904B1 (ko) * 2013-05-22 2019-12-18 삼성전자주식회사 3차원 영상 획득 장치 및 그 구동 방법
US10230934B2 (en) * 2013-06-14 2019-03-12 Microsoft Tehcnology Licensing, Llc Depth map correction using lookup tables
US9542749B2 (en) 2014-01-06 2017-01-10 Microsoft Technology Licensing, Llc Fast general multipath correction in time-of-flight imaging
US9360554B2 (en) * 2014-04-11 2016-06-07 Facet Technology Corp. Methods and apparatus for object detection and identification in a multiple detector lidar array
CN106464858B (zh) * 2014-04-26 2019-08-20 泰特拉维公司 用于3d成像中的深度感测的耐用且延伸的照明波形的方法和系统
JP6507529B2 (ja) 2014-08-29 2019-05-08 株式会社デンソー 光飛行型測距装置
US9958383B2 (en) * 2014-12-18 2018-05-01 Microsoft Technology Licensing, Llc. Range camera
JP6734292B2 (ja) * 2015-03-02 2020-08-05 ノースウェスタン ユニバーシティ Tof深度画像化システム及び方法
EP3308365B1 (en) * 2015-06-12 2023-10-25 Google LLC Using infrared images of a monitored scene to identify false alert regions
CN105004321B (zh) * 2015-07-17 2017-05-10 湖北省电力勘测设计院 顾及曝光不同步的无人机gps辅助光束法平差方法
CN107407560B (zh) * 2016-02-29 2019-09-24 松下电器产业株式会社 图像处理装置以及图像处理方法
CN105846919B (zh) * 2016-03-10 2018-03-27 中国计量科学研究院 一种调制度定标方法及装置
US10620297B2 (en) * 2016-12-22 2020-04-14 Apple Inc. Radar methods and apparatus using in phased array communication systems
CN106908041B (zh) * 2017-03-20 2019-05-21 成都通甲优博科技有限责任公司 一种近场标定相机实施远场测量的方法和装置
US11752635B2 (en) * 2017-06-02 2023-09-12 Pixart Imaging Inc. Mobile robot performing multiple detections using image frames of same optical sensor
US11691264B2 (en) * 2017-06-02 2023-07-04 Pixart Imaging Inc. Mobile robot performing multiple detections using image frames of same optical sensor
US10627518B2 (en) 2017-06-02 2020-04-21 Pixart Imaging Inc Tracking device with improved work surface adaptability
US11821985B2 (en) 2017-06-02 2023-11-21 Pixart Imaging Inc. Mobile robot performing multiple detections using image frames of same optical sensor
CN107656238B (zh) * 2017-09-12 2020-10-30 中国海洋大学 一种基于双查找表方法的高频地波雷达测向新方法
US11212512B2 (en) 2017-12-28 2021-12-28 Nlight, Inc. System and method of imaging using multiple illumination pulses
US10785422B2 (en) 2018-05-29 2020-09-22 Microsoft Technology Licensing, Llc Face recognition using depth and multi-spectral camera
US11245875B2 (en) 2019-01-15 2022-02-08 Microsoft Technology Licensing, Llc Monitoring activity with depth and multi-spectral camera
JPWO2020218283A1 (ko) * 2019-04-22 2020-10-29
US10812125B1 (en) * 2019-05-31 2020-10-20 Intel Corporation Radiation exposure control for beamforming technologies
JP7401211B2 (ja) * 2019-06-25 2023-12-19 ファナック株式会社 外光照度測定機能を備えた測距装置及び外光照度測定方法
KR102090789B1 (ko) * 2019-11-06 2020-03-19 주식회사 뉴이스트원테크 레이더 수신기의 성능 측정을 위한 모사 신호 발생 장치 및 그 방법
KR102332511B1 (ko) 2019-11-12 2021-11-29 현대모비스 주식회사 Tdc를 구비하는 라이다 시스템 및 그것의 멀티 위상 클럭 생성 방법
EP3835819A1 (en) * 2019-12-10 2021-06-16 Melexis Technologies NV Optical range calculation apparatus and method of range calculation
US11922606B2 (en) 2021-10-04 2024-03-05 Samsung Electronics Co., Ltd. Multipass interference correction and material recognition based on patterned illumination without frame rate loss
CN114721006B (zh) * 2022-06-08 2022-09-02 中国人民解放军国防科技大学 高精度单像素成像方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006201037A (ja) * 2005-01-20 2006-08-03 Ricoh Co Ltd 形状測定装置
JP2006308522A (ja) * 2005-05-02 2006-11-09 Matsushita Electric Works Ltd 空間情報検出装置
JP2008241695A (ja) * 2007-02-27 2008-10-09 Fujifilm Corp 測距装置及び測距方法
JP2009079987A (ja) * 2007-09-26 2009-04-16 Fujifilm Corp 測距方法および装置

Family Cites Families (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5081530A (en) * 1987-06-26 1992-01-14 Antonio Medina Three dimensional camera and range finder
JP2665231B2 (ja) * 1988-05-13 1997-10-22 浜松ホトニクス株式会社 光波形測定装置
US4935616A (en) * 1989-08-14 1990-06-19 The United States Of America As Represented By The Department Of Energy Range imaging laser radar
JP3305083B2 (ja) * 1993-12-22 2002-07-22 キヤノン株式会社 光レーダ
KR0134309B1 (ko) * 1994-03-11 1998-04-23 김광호 디지탈 주파수 자동조절회로
DE69635891T2 (de) * 1995-06-22 2006-12-14 3Dv Systems Ltd. Verbesserte optische kamera zur entfernungsmessung
IL114278A (en) * 1995-06-22 2010-06-16 Microsoft Internat Holdings B Camera and method
US6088086A (en) * 1995-09-11 2000-07-11 Sandia Corporation Range determination for scannerless imaging
ES2194302T3 (es) * 1998-02-08 2003-11-16 3Dv Systems Ltd Obturador optico de imagen de gran apertura.
JP2000055616A (ja) * 1998-08-12 2000-02-25 Fuji Xerox Co Ltd 3次元画像撮影方法および装置
US6118946A (en) * 1999-06-29 2000-09-12 Eastman Kodak Company Method and apparatus for scannerless range image capture using photographic film
US6856355B1 (en) * 1999-11-30 2005-02-15 Eastman Kodak Company Method and apparatus for a color scannerless range image system
US6794628B2 (en) 2000-01-03 2004-09-21 3Dv Systems, Ltd. Solid state optical shutter
US6881986B1 (en) * 2001-11-07 2005-04-19 Taiwan Semiconductor Manufacturing Company Design and fabrication method for finger n-type doped photodiodes with high sensitivity for CIS products
US6856382B2 (en) * 2003-02-06 2005-02-15 Eastman Kodak Company Formation of three-dimensional video sequences with a scannerless range imaging system
GB2404013B (en) * 2003-07-17 2006-05-31 Isis Innovation Apparatus for and method of measuring fluorescence lifetime
WO2005036372A2 (en) * 2003-10-09 2005-04-21 Honda Motor Co., Ltd. Systems and methods for determining depth using shuttered light pulses
JP4161910B2 (ja) * 2004-01-28 2008-10-08 株式会社デンソー 距離画像データ生成装置及び生成方法,プログラム
KR101159495B1 (ko) 2004-03-11 2012-06-22 이코스비젼 시스팀스 엔.브이. 파면 조정 및 향상된 3?d 측정을 위한 방법 및 장치
JP4590914B2 (ja) 2004-04-21 2010-12-01 ブラザー工業株式会社 画像表示装置
JP4293194B2 (ja) 2005-09-02 2009-07-08 財団法人雑賀技術研究所 距離測定装置、及び距離測定方法
JP2007078821A (ja) 2005-09-12 2007-03-29 Casio Comput Co Ltd 投影装置、投影方法及びプログラム
KR100842617B1 (ko) * 2007-05-29 2008-06-30 삼성전자주식회사 프로젝터
EP2026246A1 (en) * 2007-08-03 2009-02-18 Harman/Becker Automotive Systems GmbH Method and apparatus for evaluating an image
JP2009281999A (ja) * 2007-08-07 2009-12-03 Sanyo Electric Co Ltd レーザレーダ用光学部材駆動機構
WO2009073652A1 (en) * 2007-12-04 2009-06-11 Particle Measuring Systems, Inc. Two-dimensional optical imaging methods and systems for particle detection
KR101520029B1 (ko) * 2008-12-31 2015-05-15 삼성전자주식회사 고정세화 패턴을 갖는 광 변조기
KR101547327B1 (ko) * 2009-01-15 2015-09-07 삼성전자주식회사 광 이미지 변조기와 이를 포함하는 광학장치와 광 이미지 변조기의 제조 및 구동방법
KR101603778B1 (ko) * 2009-01-19 2016-03-25 삼성전자주식회사 광 이미지 셔터
KR101623960B1 (ko) 2009-06-04 2016-05-25 삼성전자주식회사 광전자 셔터, 이의 동작 방법 및 광전자 셔터를 채용한 광학 장치
KR101638974B1 (ko) 2009-06-17 2016-07-13 삼성전자주식회사 광 변조기와 그 제조 및 동작방법과 광 변조기를 포함하는 광학장치
KR101706354B1 (ko) 2009-06-24 2017-02-13 삼성전자주식회사 고속 광 변조기 및 이를 이용한 광 변조방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006201037A (ja) * 2005-01-20 2006-08-03 Ricoh Co Ltd 形状測定装置
JP2006308522A (ja) * 2005-05-02 2006-11-09 Matsushita Electric Works Ltd 空間情報検出装置
JP2008241695A (ja) * 2007-02-27 2008-10-09 Fujifilm Corp 測距装置及び測距方法
JP2009079987A (ja) * 2007-09-26 2009-04-16 Fujifilm Corp 測距方法および装置

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9310488B2 (en) 2010-09-17 2016-04-12 Samsung Electronics Co., Ltd. Apparatus and method for generating depth image
US9418425B2 (en) * 2011-10-25 2016-08-16 Samsung Electronic Co., Ltd. 3D image acquisition apparatus and method of calculating depth information in the 3D image acquisition apparatus
US20130101176A1 (en) * 2011-10-25 2013-04-25 Samsung Electronic Co.,Ltd. 3d image acquisition apparatus and method of calculating depth information in the 3d image acquisition apparatus
US9350925B2 (en) 2011-11-02 2016-05-24 Samsung Electronics Co., Ltd. Image processing apparatus and method
KR20130091194A (ko) * 2012-02-07 2013-08-16 삼성전자주식회사 3차원 영상 획득 장치 및 상기 3차원 영상 획득 장치에서 깊이 정보를 추출하는 방법
US9667944B2 (en) 2012-10-10 2017-05-30 Samsung Electronics Co., Ltd. Imaging optical system and 3D image acquisition apparatus including the imaging optical system
US9998730B2 (en) 2012-10-10 2018-06-12 Samsung Electronics Co., Ltd. Imaging optical system and 3D image acquisition apparatus including the imaging optical system
KR20140121710A (ko) * 2013-04-08 2014-10-16 삼성전자주식회사 3차원 영상 획득 장치 및 3차원 영상 획득 장치에서의 깊이 영상 생성 방법
WO2015102212A1 (ko) * 2013-12-31 2015-07-09 고려대학교 산학협력단 레이저 거리 센서에 의해 측정된 측정 거리의 거리 유형을 평가하는 방법 및 이를 이용한 이동 로봇의 위치 추정 방법
US10107897B2 (en) 2013-12-31 2018-10-23 Korea University Research And Business Foundation Method for evaluating type of distance measured by laser range finder and method for estimating position of mobile robot by using same
US10101453B2 (en) 2014-08-11 2018-10-16 Infineon Technologies Ag Time of flight apparatuses and an illumination source
US10205933B2 (en) 2014-08-29 2019-02-12 Samsung Electronics Co., Ltd. Depth image acquisition apparatus and method of acquiring depth information
US11079475B2 (en) 2015-06-08 2021-08-03 Samsung Electronics Co., Ltd. Optical device and depth information generating method
KR20160144148A (ko) * 2015-06-08 2016-12-16 삼성전자주식회사 광학 장치 및 깊이 정보 생성 방법
WO2016200166A1 (ko) * 2015-06-09 2016-12-15 (주) 솔 특성 평가를 통한 광센서 어레이 모듈 보정방법
KR20160144692A (ko) * 2015-06-09 2016-12-19 (주) 솔 광센서 어레이의 특성평가를 통한 보정방법
US11112301B2 (en) 2015-06-09 2021-09-07 Sol Inc. Method for correcting optical sensor array module through characteristic evaluation
US11680850B2 (en) 2015-06-09 2023-06-20 Sol Inc. Method for correcting optical sensor array module through characteristic evaluation
US10509125B2 (en) 2015-12-24 2019-12-17 Samsung Electronics Co., Ltd. Method and device for acquiring distance information
US10802119B2 (en) 2016-07-26 2020-10-13 Samsung Electronics Co., Ltd. Lidar device and method of measuring distance using the same
US10545237B2 (en) 2016-08-22 2020-01-28 Samsung Electronics Co., Ltd. Method and device for acquiring distance information
US10430956B2 (en) 2016-09-23 2019-10-01 Samsung Electronics Co., Ltd. Time-of-flight (TOF) capturing apparatus and image processing method of reducing distortion of depth caused by multiple reflection
WO2021086031A1 (ko) * 2019-10-31 2021-05-06 부산대학교 산학협력단 이중 가변 기반의 발진 광주파수 스캐닝 레이저 광원 그리고 그를 이용하는 측정 장치 및 중심파장 별 전파 각도 전환을 이용하는 대상물의 각도에 따른 거리 측정 장치
KR102353365B1 (ko) * 2020-08-31 2022-01-19 부산대학교 산학협력단 컬러천이 레이저의 중심파장 별 전파 각도 전환을 이용하는 대상물의 각도에 따른 거리 측정 장치
KR102302424B1 (ko) * 2020-10-26 2021-09-16 (주)시그봇 오브젝트 인식을 위한 센서 장치

Also Published As

Publication number Publication date
JP5739174B2 (ja) 2015-06-24
KR101675112B1 (ko) 2016-11-22
US20110176709A1 (en) 2011-07-21
US8611610B2 (en) 2013-12-17
JP2011149942A (ja) 2011-08-04

Similar Documents

Publication Publication Date Title
KR101675112B1 (ko) 거리 정보 추출 방법 및 상기 방법을 채용한 광학 장치
KR101722641B1 (ko) 3차원 영상 획득 장치 및 상기 3차원 영상 획득 장치에서 깊이 정보를 추출하는 방법
KR101955334B1 (ko) 3차원 영상 획득 장치 및 상기 3차원 영상 획득 장치에서 깊이 정보를 추출하는 방법
US9749613B2 (en) 3D image acquisition apparatus and method of generating depth image in the 3D image acquisition apparatus
US9310488B2 (en) Apparatus and method for generating depth image
KR101854188B1 (ko) 3차원 영상 획득 장치 및 3차원 영상 획득 장치에서 깊이 정보 산출 방법
US7379163B2 (en) Method and system for automatic gain control of sensors in time-of-flight systems
EP2477043A1 (en) 3D time-of-flight camera and method
KR102056904B1 (ko) 3차원 영상 획득 장치 및 그 구동 방법
KR102194233B1 (ko) 깊이 영상 생성 장치 및 방법
US10264240B2 (en) Method and apparatus for generating depth image
KR102194237B1 (ko) 깊이 영상 촬영 장치 및 깊이 정보 획득 방법
US20120162370A1 (en) Apparatus and method for generating depth image
US20220252730A1 (en) Time-of-flight imaging apparatus and time-of-flight imaging method
KR20170076477A (ko) 거리 정보를 획득하는 방법 및 디바이스
JP4955507B2 (ja) 測距装置および測距方法並びにプログラム
van Frank Investigation of a 3D camera

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