KR101055313B1 - Image Restoration System and Method in 3D Space Using Terahertz Pulse Reflection - Google Patents
Image Restoration System and Method in 3D Space Using Terahertz Pulse Reflection Download PDFInfo
- Publication number
- KR101055313B1 KR101055313B1 KR1020090042488A KR20090042488A KR101055313B1 KR 101055313 B1 KR101055313 B1 KR 101055313B1 KR 1020090042488 A KR1020090042488 A KR 1020090042488A KR 20090042488 A KR20090042488 A KR 20090042488A KR 101055313 B1 KR101055313 B1 KR 101055313B1
- Authority
- KR
- South Korea
- Prior art keywords
- reflected wave
- pulse
- dimensional space
- algorithm
- wave
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000006835 compression Effects 0.000 claims abstract description 29
- 238000007906 compression Methods 0.000 claims abstract description 29
- 230000003287 optical effect Effects 0.000 claims abstract description 26
- 230000001678 irradiating effect Effects 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 17
- 238000004613 tight binding model Methods 0.000 claims description 14
- 238000001228 spectrum Methods 0.000 claims description 10
- 238000010187 selection method Methods 0.000 claims description 9
- 230000003595 spectral effect Effects 0.000 claims description 7
- 238000000605 extraction Methods 0.000 claims description 4
- 230000005693 optoelectronics Effects 0.000 claims description 4
- 230000006870 function Effects 0.000 description 9
- 238000005259 measurement Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 239000013598 vector Substances 0.000 description 8
- 230000005684 electric field Effects 0.000 description 7
- 238000005457 optimization Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000001328 terahertz time-domain spectroscopy Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
-
- 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/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
- G01N21/3586—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation by Terahertz time domain spectroscopy [THz-TDS]
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N13/00—Stereoscopic video systems; Multi-view video systems; Details thereof
- H04N13/10—Processing, recording or transmission of stereoscopic or multi-view image signals
- H04N13/106—Processing image signals
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03B—GENERATION OF OSCILLATIONS, DIRECTLY OR BY FREQUENCY-CHANGING, BY CIRCUITS EMPLOYING ACTIVE ELEMENTS WHICH OPERATE IN A NON-SWITCHING MANNER; GENERATION OF NOISE BY SUCH CIRCUITS
- H03B2200/00—Indexing scheme relating to details of oscillators covered by H03B
- H03B2200/006—Functional aspects of oscillators
- H03B2200/0084—Functional aspects of oscillators dedicated to Terahertz frequencies
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Signal Processing (AREA)
- Multimedia (AREA)
- Toxicology (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
테라헤르츠 펄스 반사를 이용한 3차원 공간에서의 영상 복원 시스템 및 그 방법을 개시한다. 영상 복원 시스템은 3차원 공간에 테라헤르츠 영역의 펄스파를 조사하여 상기 3차원 공간 내의 산란 물체에 의해 반사되는 반사파를 수신하는 광학계; 및, 상기 수신된 반사파로부터 압축 센싱 알고리즘을 적용하여 상기 산란 물체의 위치를 복원하는 영상 복원부를 포함한다.Disclosed is a system and method for image reconstruction in three-dimensional space using terahertz pulse reflection. The image reconstruction system includes: an optical system for irradiating a pulse wave of a terahertz region to a three-dimensional space to receive a reflected wave reflected by a scattering object in the three-dimensional space; And an image reconstruction unit reconstructing the position of the scattering object by applying a compression sensing algorithm from the received reflected wave.
테라헤르츠, 압축 센싱, 반사형 펄스, 테라헤르츠 TDS, 영상 복원 Terahertz, Compression Sensing, Reflective Pulse, Terahertz TDS, Image Restoration
Description
본 발명의 실시예들은 테라헤르츠 파장 영역의 펄스를 이용하여 3차원 공간 내에 존재하는 물체의 위치를 복원하는 영상 복원 시스템 및 그 방법에 관한 것이다.Embodiments of the present invention relate to an image restoration system and method for restoring a position of an object existing in a three-dimensional space using pulses in the terahertz wavelength region.
테라헤르츠 영역에서 넓은 스펙트럼을 가지는 펄스파를 생성하고 수신하는 방법으로는 기계적인 스캐닝을 이용하는 테라헤르츠 TDS(terahertz time domain spectroscopy) 시스템과 AOS(asynchronous optical sampling)에 기반을 둔 테라헤르츠 TDS(terahertz time domain spectroscopy) 시스템이 있다. PCA(photo conductive antenna)는 NIR(near infrared ray) 영역의 펨토초 펄스(femtosecond pulse)와 반응하여 테라헤르츠 펄스파를 생성한다. 펄스파는 보통의 경우 수 테라헤르츠(THz)까지의 넓은 스펙트럼으로 구성되어 있다.A method of generating and receiving a wide spectrum pulse wave in the terahertz region includes a terahertz time domain spectroscopy (Terahertz time domain spectroscopy) system using mechanical scanning and a terahertz terahertz time based on asynchronous optical sampling (AOS). domain spectroscopy). A photo conductive antenna (PCA) generates terahertz pulse waves in response to femtosecond pulses in the near infrared ray (NIR) region. Pulse waves usually consist of a broad spectrum up to several terahertz (THz).
3차원 공간상에서의 물체의 모형을 복원하는 방법에는 여러 가지가 있다. 그 중에서도 펄스파를 이용하여 복원하는 방법으로는 키르히호프 방법(Kirchhoff's migration Method)이 존재한다. 키르히호프 방법은 반사되어 나온 펄스파의 시간지연에 해당하는 구형을 측정 평면에서부터 같은 거리에 존재하는 3차원 복셀(voxel)에 값들을 모두 더하는 과정으로 CT(tomography)에서의 역전사 과정(Back-projection)과 유사하다. 이러한 기술의 단점은 측정되어야 할 데이터의 개수가 많고, 복원하는 과정에 많은 시간이 필요하다는 것이다.There are many ways to restore the model of an object in three-dimensional space. Among them, Kirchhoff's migration method exists as a method of restoring using pulse waves. Kirchhoff's method adds all the values of the spherical time delay of the reflected pulse wave to the three-dimensional voxels that exist at the same distance from the measurement plane and back-projection in CT (tomography). Similar to). The disadvantages of this technique are that the number of data to be measured is large, and the time required for the restoration process is large.
이와는 다른 방식으로 2차원 영상을 복원하는 방법이 제시되었는데, 이는 퓨리에 평면에서 랜덤 패턴 마스크(Random pattern mask)를 위치하여 데이터를 한 점에서 얻은 후 이를 압축 센싱(compressed sensing) 이론에 따라 적은 수의 데이터로 2차원 영상을 얻어내는 방법이다. 이는 빠른 시간 안에 본래의 영상에 가까운 결과를 보여준다. 그러나, 압축 센싱을 이용한 영상 복원 방법은 아직 3차원 공간 상에서 영상을 복원하는 과정에 적용되지 못하고 있는 실정이다.A method of reconstructing two-dimensional images in a different manner has been proposed, in which a random pattern mask is located in a Fourier plane to obtain data at one point, and then a small number of data are obtained according to the compressed sensing theory. A method of obtaining a two-dimensional image from data. This shows a result close to the original image in a short time. However, an image restoration method using compression sensing has not yet been applied to a process of restoring an image in a three-dimensional space.
본 발명의 일실시예는 압축 센싱 이론을 바탕으로 한 알고리즘을 이용하여 3차원 공간에 존재하는 물체의 위치를 보다 정확하고 빠르게 복원하기 위한 영상 복원 시스템 및 그 방법을 제공한다.One embodiment of the present invention provides an image restoration system and method for more accurately and quickly reconstructing the position of an object in three-dimensional space using an algorithm based on compression sensing theory.
본 발명의 일실시예에 따른 영상 복원 시스템은 3차원 공간에 테라헤르츠 영역의 펄스파를 조사하여 상기 3차원 공간 내의 산란 물체에 의해 반사되는 반사파를 수신하는 광학계; 및, 상기 수신된 반사파로부터 압축 센싱 알고리즘을 적용하여 상기 산란 물체의 위치를 복원하는 영상 복원부를 포함한다.An image reconstruction system according to an embodiment of the present invention includes an optical system for receiving a reflected wave reflected by a scattering object in the three-dimensional space by irradiating a pulse wave of the terahertz region in the three-dimensional space; And an image reconstruction unit reconstructing the position of the scattering object by applying a compression sensing algorithm from the received reflected wave.
본 발명의 일실시예에서 상기 광학계는 상기 펄스파를 조사하는 조사 위치와 상기 반사파를 수신하는 수신 위치가 동일한 위치 또는 일정 범위 이내의 위치에 존재하거나, 상기 반사파를 수신하는 수신 위치가 상기 펄스파를 조사하는 조사 위치와 반대편 혹은 다른 위치에 존재하도록 한다.In one embodiment of the present invention, the optical system is located at the same position or within a predetermined range of the irradiation position for irradiating the pulse wave and the receiving position for receiving the reflected wave, or the receiving position for receiving the reflected wave is the pulse wave Be on the other side of the survey or in a different location.
본 발명의 일실시예에서 상기 광학계는 상기 펄스파를 생성하는 펄스 생성부와 상기 반사파를 수신하는 펄스 수신부를 포함한 구조의 광전자 안테나로 구성된다.In an embodiment of the present invention, the optical system includes an optoelectronic antenna having a structure including a pulse generator for generating the pulse wave and a pulse receiver for receiving the reflected wave.
본 발명의 일실시예에서 상기 영상 복원부는 상기 반사파를 스펙트럼 정보로 변환한 후 상기 스펙트럼 정보와 일정 조건의 그린 함수를 상기 압축 센싱 알고리즘에 적용하여 상기 산란 물체의 위치를 산출한다.In one embodiment of the present invention, the image reconstruction unit converts the reflected wave into spectral information and calculates the position of the scattering object by applying the spectral information and the green function of a predetermined condition to the compression sensing algorithm.
본 발명의 일실시예에서 상기 영상 복원부는 상기 반사파를 퓨리에 변환을 통해 시간 도메인 스펙트럼으로 변환하고 상기 변환된 시간 도메인 스펙트럼으로부터 상기 반사파에 대한 선형적인 행렬을 구한 후 상기 선형적인 행렬에 상기 그린 함수의 행렬을 적용한다.In one embodiment of the present invention, the image reconstructor converts the reflected wave into a time domain spectrum through a Fourier transform, obtains a linear matrix for the reflected wave from the converted time domain spectrum, and then applies the green function to the linear matrix. Apply the matrix.
본 발명의 일실시예에서 상기 압축 센싱 알고리즘은 압축 센싱 이론에 적용되는 알고리즘으로 선형 프로그래밍, 포커스 알고리즘, 추출 알고리즘, 기저함수 선택 방법, 직교 기저함수 선택 방법, 약한 직교 기저함수 선택 방법 중 적어도 하나를 사용할 수 있다.In one embodiment of the present invention, the compression sensing algorithm is an algorithm applied to a compression sensing theory, and includes at least one of linear programming, focus algorithm, extraction algorithm, basis function selection method, orthogonal basis function selection method, and weak orthogonal basis function selection method. Can be used.
본 발명의 일실시예에 따른 영상 복원 방법은 광학계와 영상 복원부를 포함한 영상 복원 시스템의 영상 복원 방법에 있어서, 상기 광학계에서 3차원 공간에 테라헤르츠 영역의 펄스파를 조사하여 상기 3차원 공간 내의 산란 물체에 의해 반사되는 반사파를 수신하는 단계; 및, 상기 영상 복원부에서 상기 수신된 반사파로부터 압축 센싱 알고리즘을 적용하여 상기 산란 물체의 위치를 복원하는 단계를 포함한다.An image restoration method according to an embodiment of the present invention is an image restoration method of an image restoration system including an optical system and an image restoration unit, wherein the optical system irradiates a pulse wave of a terahertz region to a three-dimensional space and scatters the three-dimensional space. Receiving a reflected wave reflected by an object; And reconstructing the position of the scattering object by applying a compression sensing algorithm from the received reflected wave in the image reconstruction unit.
본 발명의 일실시예에 따르면, 전자기파 중 테라헤르츠에서 넓은 영역을 갖는 펄스파를 이용하여 3차원 공간 상에서의 물체의 위치를 복원함으로써 보다 향상된 해상도의 영상을 얻을 수 있다.According to an embodiment of the present invention, an image having a higher resolution may be obtained by restoring the position of an object in a three-dimensional space by using a pulse wave having a large area in terahertz among electromagnetic waves.
본 발명의 일실시예에 따르면, 압축 센싱 이론을 바탕으로 한 알고리즘을 이용하여 테라헤르츠 펄스에 대한 반사파(pulse echo)를 적용함으로써 적은 수의 데이터로 3차원 공간의 물체를 복원할 수 있어 영상 복원에 소요되는 시간을 단축시킬 수 있다.According to an embodiment of the present invention, by applying a reflected echo to a terahertz pulse using an algorithm based on a compression sensing theory, it is possible to recover an object in a three-dimensional space with a small number of data. This can shorten the time required.
이하에서, 본 발명에 따른 실시예들을 첨부된 도면을 참조하여 상세하게 설명한다. 그러나, 본 발명이 실시예들에 의해 제한되거나 한정되는 것은 아니다. 각 도면에 제시된 동일한 참조 부호는 동일한 부재를 나타낸다.Hereinafter, embodiments of the present invention will be described in detail with reference to the accompanying drawings. However, the present invention is not limited or limited by the embodiments. Like reference numerals in the drawings denote like elements.
도 1은 본 발명의 일실시예에 따른 3차원 공간에서의 물체의 위치를 복원하는 영상 복원 시스템을 도시한 것이다. 도 1을 참조하면, 본 발명의 일실시예에 따른 영상 복원 시스템은 광학계(110)와, 영상 복원부(130)를 포함한다.1 illustrates an image reconstruction system for reconstructing a position of an object in a three-dimensional space according to an embodiment of the present invention. Referring to FIG. 1, an image restoration system according to an embodiment of the present invention includes an
도 2는 본 발명의 일실시예에 따른 영상 복원 시스템의 영상 복원 과정을 도시한 것이다. 여기서, 본 발명의 일실시예에 따른 영상 복원 과정은 도 1의 영상 복원 시스템에 의해 실행될 수 있다.2 illustrates an image restoration process of an image restoration system according to an embodiment of the present invention. Here, the image restoration process according to an embodiment of the present invention may be performed by the image restoration system of FIG. 1.
도 1과 도 2를 참조하면, 광학계(110)는 테라헤르츠 영역의 펄스파(이하, '테라헤르츠 펄스'라 약칭함)를 3차원 공간으로 방출한 후(S201) 상기 3차원 공간 내에 존재하는 물체(이하, '산란 물체'라 칭함)에 의해 반사되는 반사파를 수신하는(S202) 역할을 수행한다. 광학계(110)는 근적외선 영역의 펨토초 펄스와 반응하여 광원인 테라헤르츠 펄스를 생성하는 광전자 안테나(PCA)로 구성된다. 즉, 광학계(110)는 테라헤르츠 펄스가 생성되는 펄스 생성부(111)와, 테라헤르츠 펄스가 산란 물체에 의해서 반사되는 반사파를 수신하는 펄스 수신부(113)를 포함한 구조의 광전자 안테나를 사용할 수 있다. 이때, 펄스 생성부(111)와 펄스 수신부(113)는 하나의 구조물로 구성되어 테라헤르츠 펄스를 조사하는 펄스 생성부(111)의 조사 위치와 반사파를 수신하는 펄스 수신부(113)의 수신 위치가 동일한 위치 또는 일정 범위 이내의 위치에서 이루어지도록 한다.Referring to FIGS. 1 and 2, the
영상 복원부(130)는 광학계(110)를 통해 수신되는 반사파로부터 압축 센싱 이론을 바탕으로 하는 알고리즘을 적용하여 3차원 공간 상의 산란 물체의 위치를 복원하는(S203) 역할을 수행한다. 영상 복원부(130)는 테라헤르츠 TDS(time domain spectroscopy) 시스템을 포함할 수 있다. 이때, TDS 시스템은 AOS(asynchronous optical sampling)에 기반을 둔 테라헤르츠 TDS 시스템을 이용할 수 있다. 영상 복원부(130)는 수신되는 반사파를 스펙트럼 정보로 변환한 후 상기 변환된 스펙트럼 정보와 일정 조건의 그린 함수를 압축 센싱 이론에 적용하여 상기 산란 물체의 위치를 산출할 수 있다.The
본 발명의 일실시예는 전자기파 중 테라헤르츠에서 넓은 빔 영역을 갖는 펄스파를 이용함으로써 테라헤르츠 펄스의 스펙트럼을 해석할 때 보다 많은 정보를 포함할 수 있다. 또한, 본 발명의 일실시예는 테라헤르츠 펄스에 대한 반사파(pulse echo)의 스펙트럼을 압축 센싱 이론에 적용함으로써 3차원 공간 내 물체의 위치를 정확하게 복원할 수 있다.One embodiment of the present invention may include more information when analyzing the spectrum of the terahertz pulse by using a pulse wave having a wide beam area in the terahertz of the electromagnetic waves. In addition, one embodiment of the present invention can accurately restore the position of the object in the three-dimensional space by applying the spectrum of the pulse (pulse echo) for the terahertz pulse to the compression sensing theory.
도 3 내지 도 5를 참조하여, 3차원 공간 내 산란 물체의 위치를 복원하는 영상 복원 과정을 상세하게 설명한다.3 to 5, the image restoration process of restoring the position of the scattering object in the 3D space will be described in detail.
도 3을 참조하면, 광학계(110)의 펄스 생성부(111)는 광원인 테라헤르츠 펄스를 생성한다. 펄스 생성부(111)는 영상 복원부(130)의 제어 하에 측정면(301)의 임의 조사 위치(301a)에서 구형의 테라헤르츠 펄스를 생성하여 방출한다. 본 발명의 일실시예는 3차원 공간을 복원하기 위하여 펄스 수신부(113)에서 구형의 반사파를 수신하여 이용하기 때문에 펄스 생성부(111)에서 생성된 테라헤르츠 펄스를 렌즈 등에 의해 평행빔으로 변환하지 않아도 된다. 또한, 본 발명의 일실시예는 측정면(301)에서 조사 위치(301a) 및 수신 위치(301b)의 정확한 위치를 알 수 있는 조건이라면 3차원 공간(303)에 대한 측정면(301)이 반드시 평면일 필요는 없다.Referring to FIG. 3, the
자유 공간인 3차원 공간(303) 상에서 전파되는 테라헤르츠 펄스는 전자기장의 한 종류로써 멕스웰 방정식(Maxwell's equation)을 만족시키고, 또한 스칼라 헬름홀츠 방정식(scalar Helmhortz equation)을 만족시킨다. 전자기파는 자유 공간에서 전파되는 도중 이동경로에 존재하는 산란 물체(305)에 의해서 산란된다. 이때, 산란 물체(305)를 지난 후의 전자기파는 산란에 의한 전기장과 입사되는 전기장이 합쳐져서 표현될 수 있으며 수학식 1과 같이 정의할 수 있다.The terahertz pulse propagating in the three-
여기서, ψT는 산란 후의 측정되는 전기장, ?娶?i는 조사 위치(301a)로부터 입사되는 전기장, ?娶?s는 산란 물체(305)에 의해 산란된 전기장이다. 산란된 전기장(ψs)은 본 근사식(Born approximation)과 헬름홀츠 방정식에 의해서 수학식 2와 같이 정의할 수 있다.Where ψ T is the electric field measured after scattering,? i is the electric field incident from the
수학식 2의 k(k=w/c)는 파장벡터(wavenumber vector), gk(r|r')는 그린함수(Green's function)이다. 산란된 전기장(ψs)에 의한 반사파는 산란 물체의 특성인 f(r)에 따라 영향을 받는다. 산란 물체(305)이 흡수율이 크면 반사파는 세기가 매우 약해져서 검출이 안될 것이다.K (k = w / c) in Equation 2 is a wavenumber vector, and g k (r | r ') is a Green's function. The reflected wave due to the scattered electric field ψ s is affected by f (r), which is a characteristic of the scattering object. If the
본 발명의 일실시예는 테라헤르츠 펄스에 대한 반사파를 이용하므로 수학식2를 이용하여 반사파를 측정할 수 있다. 광학계(110)의 펄스 수신부(113)는 테라헤르츠 펄스를 조사한 조사 위치(301a)와 동일한 위치의 수신 위치(301b)에서 산란 물체(305)에 의해 산란되어 반사되는 반사파를 측정한다. 광학계(110)는 영상 복원부(130)의 제어 하에 측정면(301) 내에서 위치를 이동해가면서 반사파를 측정 및 수집한다.Since the exemplary embodiment of the present invention uses the reflected wave for the terahertz pulse, the reflected wave can be measured using Equation 2. The
그린 함수는 일반적인 선형 편미분 방정식에서 델타 함수를 도출하는 기본함수로써 보통 해를 구할 때 사용된다. 수학식 3은 헬름홀츠 방정식에 해당하는 그린 함수이다.Green's function is a basic function that derives the delta function from a general linear partial differential equation and is usually used to solve the solution. Equation 3 is a green function corresponding to the Helmholtz equation.
여기서, r, r'는 광학계(110)의 조사 위치(301a)와 수신 위치(301b)를 각각 나타내며 직교 좌표계로 생각해서 이용할 수 있다. 본 발명에서 조사 위치(301a)와 수신 위치(301b)는 동일한 위치 또는 일정 범위 이내의 위치에 놓임으로써 수학식 1을 단순화 시키고 단순화된 식을 수신 위치(301b)에서 측정된 반사파에 대한 데이터 값(y)으로 나타내면 수학식 4와 같이 정의할 수 있다.Here, r and r 'represent the
여기서, S(k)는 광원인 테라헤르츠 펄스가 갖는 파형의 스펙트럼을 나타낸다. 수학식 4는 파장과 반비례하는 파장벡터(wavenumber) 도메인을 가진다. 파형에 따라서 측정하게 되는 반사파의 데이터 크기나 노이즈 정도가 변할 수 있다. 수학식 4와 같이 3차원 공간(303) 상에 존재하는 산란 물체(305)에 대한 정보를 측정면(301)에서 수신되는 반사파를 통해서 얻을 수 있다.Here, S (k) represents the spectrum of the waveform of the terahertz pulse which is a light source. Equation 4 has a wavenumber domain in inverse proportion to the wavelength. Depending on the waveform, the data size and noise level of the reflected wave to be measured may change. As shown in Equation 4, information about the
또한, 도 3에서, 광학계(110)의 조사 위치(301a)는 그대로 유지하고 수신 위치(301b)를 조사 위치(301a)가 위치한 측정면(301)의 반대편 즉, 3차원 공간(303)에서 조사 위치(301a)와 서로 대응되는 위치에 위치시킬 수 있다. 조사 위치(301a)의 반대편에 위치한 수신 위치(301b)는 수학식 5와 같이 정의할 수 있다.3, the
이와 같은 수신 위치(301b)의 변경은 실제 알고리즘에 적용되는 선형 행렬 계산의 그린함수 행렬을 변경시키나, 똑같은 알고리즘이 적용되어 복원될 수 있다.Such a change in the
영상 복원부(130)는 펄스 수신부(113)에 의해서 측정되는 반사파를 디지털화된 데이터로 저장한다. 이때, 반사파에 대한 데이터(401)는 도 4와 같이 모형화 될 수 있다. 영상 복원부(130)는 반사파의 데이터(401)를 수학식 4에 맞게 퓨리에 변환을 통해 스펙트럼 정보로 변환할 수 있다. 또한, 수학식 4에서 알 수 있듯이 반사파에 대한 데이터는 연속적인 시간 도메인에 존재하는데, 영상 복원부(130)는 반사파를 이산화하여 실제 알고리즘에 적합한 형태로 바꿀 수 있고 선형적인 행렬 계산으로 변환할 수 있다.The image reconstructor 130 stores the reflected wave measured by the
여기서, n은 측정하는 위치 즉, 광학계(110)를 이동하면서 반사파의 수신이 이루어진 수신 위치(301b)에 대한 개수를 의미한다. n이 개수를 의미하더라도 여기에 해당하는 그린 함수의 행렬은 각각의 펄스 수신부(113)에 대한 위치정보를 그린 함수의 제곱 꼴로 가지고 있다.Here, n denotes the number of the measurement position, that is, the
따라서, 수학식 6은 수학식 4와 동치된다. 여기서, 그린 함수의 행렬은 수학식 7과 같이 구성된다. 지수 n은 각각의 수신 위치(301b)를 나타낸다. 기본적으로, 반사파에 대한 데이터는 y이며 이들은 구하고자 하는 f 벡터와 그린 함수 행렬과의 행렬 곱이 된다. f 벡터가 듬성하게(sparse) 분포되어 있는 산란 물체라고 가정할 경우 압축 센싱 이론에 따라 적은 수의 y로도 본래의 물체를 완벽에 가깝게 복원할 수 있다.Therefore, equation (6) is equivalent to equation (4). Here, the matrix of the green function is configured as shown in Equation (7). The exponent n represents each receiving
압축 센싱 이론을 적용하기 위해서는 한가지 조건이 필요하다. 선형 연산자에서 각각의 컬럼 벡터(column vector)가 서로의 스칼라 곱셈의 값이 작아야 한다. 컬럼 벡터간의 유사성을 상호 유사성(mutual coherence)이라는 인자를 기준으로 해석적인 해를 구한다. 상호 유사성 인자로 어느 정도의 범위까지 복원이 가능한지 알아낼 수 있다. 상호 유사성을 구하는 식은 수학식 8과 같다.In order to apply the compression sensing theory, one condition is required. In a linear operator, each column vector must have a small value of each other's scalar multiplication. Analyze the similarity between column vectors based on the factor of mutual coherence. The mutual similarity factor can be used to determine the extent to which restoration is possible. The equation for obtaining mutual similarity is shown in Equation 8.
여기서, rn은 펄스 수신부(113)의 위치를, r'p, r'q은 산란 물체의 위치를 의미한다. 도 5를 참조하면, 상호 유사성은 전체의 면적과 두 개의 산란 물체(r'p, r'q) 사이에서 같은 거리에 떨어진 평면(503)과 측정면(501)이 만나는 선 위의 임의의 두 점(rn)이 산란 물체(r'p, r'q)와 떨어진 거리의 비율이 상호 유사성이 된다. 이에 따른 상호 유사성은 매우 작은 값이 되므로 본 발명의 일실시예는 압축 센싱 이론을 효과적으로 적용할 수 있다.Here, r n means the position of the
여기서, ∥f∥0는 산란 물체의 개수이다. 복원 가능한 산란 물체의 개수가 상호 유사성과 연관이 있으나, 상호 유사성이 매주 작으므로 많은 수의 산란 물체를 복원할 수 있다. 실제로 상호 유사성은 수학식 6에서 측정되는 벡터의 차원(dimension)과 관련이 있으므로, 복원 정도의 범위는 변수들이 변함에 따라 변할 수 있다.Where ∥f∥ 0 is the number of scattering objects. Although the number of recoverable scattering objects is related to each other, the number of scattering objects can be restored because the mutual similarity is very small. In fact, since the mutual similarity is related to the dimension of the vector measured in Equation 6, the range of the degree of reconstruction may change as the variables change.
상기한 과정에서 얻어진 반사파에 대한 데이터(즉, 수학식 4)를 미리 계산해 놓은 그림 함수의 행렬을 통해 정해진 위치의 그린 함수를 선택하여 압축 센싱 이론을 바탕으로 한 알고리즘으로 산란물체를 복원한다.The scattering object is reconstructed by an algorithm based on compression sensing theory by selecting a green function at a predetermined position through a matrix of a picture function previously calculated for the reflected wave data (ie, Equation 4) obtained in the above process.
여기서, f는 복원할 산란 물체의 1차원 벡터를 뜻한다. 수학식 10은 0이 아닌 값(Nonzero value)을 최소로 하는 L0 문제이다. L0 최소화는 복원할 산란 물체에 대한 벡터의 위치 특성이 고려되지 않으므로 유일한 해를 구하기 어렵다.Here, f means a one-dimensional vector of the scattering object to be restored. Equation 10 is a L0 problem of minimizing a nonzero value. The L0 minimization is difficult to find the only solution because the positional characteristics of the vector relative to the scattering object to be restored are not taken into account.
그러나, 3차원 공간 상의 산란 물체가 y(i)=G(i)f를 만족하는 범위에서 최대한 듬성하게 분포되어 있을 경우 압축 센싱 이론을 적용하면, 수학식 10이 수학식 11과 같이 변환된다.However, when the scattering object in the three-dimensional space is distributed as maximally as possible in the range satisfying y (i) = G (i) f, if the compression sensing theory is applied, Equation 10 is converted to Equation 11.
본래의 0이 아닌 값(Nonzero value)만 복원하는 최적화 문제였으나 압축 센싱 이론에 의해 L1 최적화 문제로 변화되며 이는 L0 문제보다 특정 해를 구하기가 쉬운 문제이다.It was an optimization problem that restores the original nonzero value, but it is changed to the L1 optimization problem by compression sensing theory, which is easier to find a specific solution than the L0 problem.
본 발명의 일실시예에서 L1 최적화 알고리즘으로는 L1 최적화에 쓰이는 선형 프로그래밍(Linear programming), 포커스 알고리즘(FOCUSS Algorithm) 등을 포함하는 방법, 필요한 성분들만 분리하는 추출 알고리즘(Greedy Algorithm), 기저함수 선택 방법(Matching Pursuit Method), 직교 기저함수 선택 방법(Orthogonal Matching Pursuit Method), 약한 직교 기저함수 선택 방법(Weak Orthogonal Matching Pursuit Method) 등을 사용할 수 있다.In one embodiment of the present invention, the L1 optimization algorithm includes a linear programming method used for L1 optimization, a focus algorithm (FOCUSS Algorithm), etc., an extraction algorithm (Greedy Algorithm) for separating only necessary components, and a basis function selection. Matching Pursuit Method, Orthogonal Matching Pursuit Method, Weak Orthogonal Matching Pursuit Method, etc. may be used.
본 발명의 일실시예는 압축 센싱 이론을 바탕으로 한 L1 최적화 알고리즘을 사용함으로써 적은 수의 데이터로 3차원 공간 내 산란 물체의 위치를 복원할 수 있다.One embodiment of the present invention can restore the position of the scattering object in the three-dimensional space with a small number of data by using the L1 optimization algorithm based on the compression sensing theory.
본 발명의 실시예들은 다양한 컴퓨터로 구현되는 동작을 수행하기 위한 프로그램 명령을 포함하는 컴퓨터 판독 가능 매체를 포함한다. 상기 컴퓨터 판독 가능 매체는 프로그램 명령, 데이터 파일, 데이터 구조 등을 단독으로 또는 조합하여 포함할 수 있다. 상기 매체는 프로그램 명령은 본 발명을 위하여 특별히 설계되고 구성된 것들이거나 컴퓨터 소프트웨어 당업자에게 공지되어 사용 가능한 것일 수도 있다. 컴퓨터 판독 가능 기록 매체의 예에는 하드 디스크, 플로피 디스크 및 자기 테이프와 같은 자기 매체(magnetic media), CD-ROM, DVD와 같은 광기록 매체(optical media), 플롭티컬 디스크(floptical disk)와 같은 자기-광 매체(magneto-optical media), 및 롬(ROM), 램(RAM), 플래시 메모리 등과 같은 프로그램 명령을 저장하고 수행하도록 특별히 구성된 하드웨어 장치가 포함된다. 프로그램 명령의 예에는 컴파일러에 의해 만들어지는 것과 같은 기계어 코드뿐만 아니라 인터프리터 등을 사용해서 컴퓨터에 의해서 실행될 수 있는 고급 언어 코드를 포함한다.Embodiments of the present invention include computer readable media including program instructions for performing various computer implemented operations. The computer readable medium may include program instructions, data files, data structures, etc. alone or in combination. The media may be program instructions that are specially designed and constructed for the present invention or may be available to those skilled in the art of computer software. Examples of computer-readable recording media include magnetic media such as hard disks, floppy disks, and magnetic tape, optical media such as CD-ROMs, DVDs, and magnetic disks, such as floppy disks. Magneto-optical media, and hardware devices specifically configured to store and execute program instructions, such as ROM, RAM, flash memory, and the like. Examples of program instructions include not only machine code generated by a compiler, but also high-level language code that can be executed by a computer using an interpreter or the like.
이상과 같이 본 발명에서는 구체적인 구성 요소 등과 같은 특정 사항들과 한정된 실시예 및 도면에 의해 설명되었으나 이는 본 발명의 보다 전반적인 이해를 돕기 위해서 제공된 것일 뿐, 본 발명은 상기의 실시예에 한정되는 것은 아니며, 본 발명이 속하는 분야에서 통상적인 지식을 가진 자라면 이러한 기재로부터 다양한 수정 및 변형이 가능하다.As described above, the present invention has been described by specific embodiments such as specific components and the like. For those skilled in the art to which the present invention pertains, various modifications and variations are possible.
따라서, 본 발명의 사상은 설명된 실시예에 국한되어 정해져서는 아니되며, 후술하는 특허청구범위뿐 아니라 이 특허청구범위와 균등하거나 등가적 변형이 있 는 모든 것들은 본 발명 사상의 범주에 속한다고 할 것이다.Therefore, the spirit of the present invention should not be limited to the described embodiments, and all the things that are equivalent to or equivalent to the scope of the claims as well as the claims to be described later belong to the scope of the present invention. will be.
도 1은 본 발명의 일실시예에 따른 3차원 공간에서의 물체의 위치를 복원하는 영상 복원 시스템의 구성을 도시한 도면이다.1 is a diagram illustrating a configuration of an image restoration system for restoring a position of an object in a three-dimensional space according to an embodiment of the present invention.
도 2는 본 발명의 일실시예에서 3차원 공간에서의 물체의 위치를 복원하는 과정을 설명하기 위한 도면이다.2 is a view for explaining a process of restoring the position of an object in a three-dimensional space in an embodiment of the present invention.
도 3은 테라헤르츠 영역의 펄스를 이용하여 3차원 공간 내 물체의 위치를 산출하는 과정을 설명하기 위한 도면이다.FIG. 3 is a diagram for describing a process of calculating a position of an object in a three-dimensional space using pulses in the terahertz region.
도 4는 3차원 공간으로부터 측정된 반사파를 설명하기 위한 도면이다.4 is a diagram for explaining reflected waves measured from a three-dimensional space.
도 5는 압축 센싱 이론에 적용하여 3차원 공간 내 물체의 위치를 산출하는 과정을 설명하기 위한 도면이다.5 is a view for explaining a process of calculating the position of the object in the three-dimensional space by applying the compression sensing theory.
<도면의 주요 부분에 대한 부호의 설명><Explanation of symbols for the main parts of the drawings>
110: 광학계110: optical system
111: 펄스 생성부111: pulse generator
113: 펄스 수신부113: pulse receiving unit
130: 영상 복원부130: image restoration unit
Claims (14)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020090042488A KR101055313B1 (en) | 2009-05-15 | 2009-05-15 | Image Restoration System and Method in 3D Space Using Terahertz Pulse Reflection |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020090042488A KR101055313B1 (en) | 2009-05-15 | 2009-05-15 | Image Restoration System and Method in 3D Space Using Terahertz Pulse Reflection |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20100123322A KR20100123322A (en) | 2010-11-24 |
KR101055313B1 true KR101055313B1 (en) | 2011-08-09 |
Family
ID=43408018
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020090042488A KR101055313B1 (en) | 2009-05-15 | 2009-05-15 | Image Restoration System and Method in 3D Space Using Terahertz Pulse Reflection |
Country Status (1)
Country | Link |
---|---|
KR (1) | KR101055313B1 (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10153547A (en) * | 1996-11-06 | 1998-06-09 | Lucent Technol Inc | Analysis processing method for medium |
US20050100866A1 (en) * | 1999-07-23 | 2005-05-12 | Teraview Limited | Radiation probe and detecting tooth decay |
JP2005265793A (en) * | 2004-03-22 | 2005-09-29 | Pentax Corp | Apparatus for tracking specified substance |
JP2008046574A (en) * | 2006-08-18 | 2008-02-28 | Junichi Nishizawa | Frequency variable terahertz wave endoscope |
-
2009
- 2009-05-15 KR KR1020090042488A patent/KR101055313B1/en active IP Right Grant
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH10153547A (en) * | 1996-11-06 | 1998-06-09 | Lucent Technol Inc | Analysis processing method for medium |
US20050100866A1 (en) * | 1999-07-23 | 2005-05-12 | Teraview Limited | Radiation probe and detecting tooth decay |
JP2005265793A (en) * | 2004-03-22 | 2005-09-29 | Pentax Corp | Apparatus for tracking specified substance |
JP2008046574A (en) * | 2006-08-18 | 2008-02-28 | Junichi Nishizawa | Frequency variable terahertz wave endoscope |
Also Published As
Publication number | Publication date |
---|---|
KR20100123322A (en) | 2010-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lindell et al. | Wave-based non-line-of-sight imaging using fast fk migration | |
Iseringhausen et al. | Non-line-of-sight reconstruction using efficient transient rendering | |
Haltmeier et al. | Compressed sensing and sparsity in photoacoustic tomography | |
Lobel et al. | A new regularization scheme for inverse scattering | |
Bao et al. | Imaging of local surface displacement on an infinite ground plane: the multiple frequency case | |
US11412118B2 (en) | Electromagnetic wave phase/amplitude generation device, electromagnetic wave phase/amplitude generation method, and electromagnetic wave phase/amplitude generation program | |
Wang et al. | Data-driven multichannel seismic impedance inversion with anisotropic total variation regularization | |
Recur et al. | Terahertz radiation for tomographic inspection | |
Aldawood et al. | The possibilities of compressed-sensing-based Kirchhoff prestack migration | |
Carin et al. | In situ compressive sensing | |
Rigaud et al. | Reflective imaging solved by the radon transform | |
KR101055313B1 (en) | Image Restoration System and Method in 3D Space Using Terahertz Pulse Reflection | |
Galdi et al. | Moderately rough surface underground imaging via short-pulse quasi-ray Gaussian beams | |
Zhu et al. | A three-step reconstruction method for fluorescence molecular tomography based on compressive sensing | |
Galdi et al. | Multifrequency subsurface sensing in the presence of a moderately rough air–soil interface via quasi-ray Gaussian beams | |
JP7101811B2 (en) | Methods for image reconstruction of scanners and structure of objects | |
KR102180569B1 (en) | Apparatus for Enhancing Contrast and Resolution in a Grating Interferometer by Machine Learning | |
Shustak et al. | Time reversal based detection of subsurface scatterers | |
KR20130008739A (en) | Method for compressive diffuse optical tomography which is non-iterative and exact reconstruction using joint sparsity and apparatus thereof | |
Baumgart et al. | Procedure enabling simulation and in-depth analysis of optical effects in camera-based time-of-flight sensors | |
Elbau et al. | Reconstructing the optical parameters of a layered medium with optical coherence elastography | |
Chen et al. | Robust speckle-autocorrelation non-line-of-sight imaging with generative adversarial networks | |
Dourthe et al. | Regularized bi-conjugate gradient algorithm for tomographic reconstruction of buried objects | |
Gaullier et al. | Introducing shape constraints into object-based traveltime tomography | |
LaHaie et al. | 3D conformal inverse synthetic aperture radar (ISAR) imaging using basis pursuit |
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: 20140724 Year of fee payment: 4 |
|
FPAY | Annual fee payment |
Payment date: 20150729 Year of fee payment: 5 |
|
FPAY | Annual fee payment |
Payment date: 20160727 Year of fee payment: 6 |
|
FPAY | Annual fee payment |
Payment date: 20180123 Year of fee payment: 7 |
|
FPAY | Annual fee payment |
Payment date: 20180725 Year of fee payment: 8 |