KR102142616B1 - Method for utilizing of hyperspectral image using optimal band ratio - Google Patents

Method for utilizing of hyperspectral image using optimal band ratio Download PDF

Info

Publication number
KR102142616B1
KR102142616B1 KR1020190000699A KR20190000699A KR102142616B1 KR 102142616 B1 KR102142616 B1 KR 102142616B1 KR 1020190000699 A KR1020190000699 A KR 1020190000699A KR 20190000699 A KR20190000699 A KR 20190000699A KR 102142616 B1 KR102142616 B1 KR 102142616B1
Authority
KR
South Korea
Prior art keywords
hyperspectral
band ratio
wavelengths
value
module
Prior art date
Application number
KR1020190000699A
Other languages
Korean (ko)
Other versions
KR20200084973A (en
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 KR1020190000699A priority Critical patent/KR102142616B1/en
Publication of KR20200084973A publication Critical patent/KR20200084973A/en
Application granted granted Critical
Publication of KR102142616B1 publication Critical patent/KR102142616B1/en

Links

Images

Classifications

    • 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
    • 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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10032—Satellite or aerial image; Remote sensing
    • G06T2207/10036—Multispectral image; Hyperspectral image

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Image Processing (AREA)

Abstract

본 발명은 (a) 초분광 관측값 입력모듈(310)에 미리 결정된 다수의 지점 각각의 초분광 관측값(Y)이 입력되는 단계; (b) 촬영영상 입력모듈(210)에 초분광정보를 갖는 촬영영상이 입력되고, 상기 미리 결정된 다수의 지점 각각에서 다수의 파장(x)에 대한 화소값인 초분광 촬영값(I(x))이 확인되는 단계; (c) 밴드비 연산모듈(320)이, 상기 초분광 촬영값에 포함된 다수의 파장(x) 중 임의의 2개의 파장의 쌍(λ1, λ2)을 선택하고, 각각에 대한 초분광 촬영값(I(λ1), I(λ2))을 확인하는 단계; (d) 상기 밴드비 연산모듈(320)이, 상기 확인된 2개의 초분광 촬영값(I(λ1), I(λ2))에서 어느 하나를 다른 하나로 나누어 밴드비(X)를 연산하는 단계; 및 (e) 회귀분석모듈(330)이, 상기 초분광 관측값(Y)과 상기 연산된 밴드비(X)가 선형관계에 있음을 가정하고, 상기 다수의 지점에서의 상기 초분광 관측값(Y)과 상기 연산된 밴드비(X)를 이용하여, 회귀식(Y=aX+b)의 a와 b를 연산함으로써, 회귀식(Y=aX+b)을 생성하는 단계를 포함하는, 밴드비를 이용한 초분광영상 활용 방법을 제공한다.The present invention includes the steps of: (a) inputting the hyperspectral observation value Y of each of a plurality of predetermined points into the hyperspectral observation value input module 310; (b) A photographed image having hyperspectral information is input to the photographed image input module 210, and a hyperspectral photographing value (I(x)) that is a pixel value for a plurality of wavelengths (x) at each of the plurality of predetermined points ) Is confirmed; (c) The band ratio calculation module 320 selects a pair of arbitrary two wavelengths (λ 1 , λ 2 ) among a plurality of wavelengths (x) included in the hyperspectral photographing value, and hyperspectral for each Checking the photographed values I(λ 1 ) and I(λ 2 ); (d) The band ratio calculation module 320 calculates a band ratio (X) by dividing any one from the identified two hyperspectral photographing values (I(λ 1 ) and I(λ 2 )) into the other. step; And (e) the regression analysis module 330 assumes that the hyperspectral observation value (Y) and the calculated band ratio (X) are in a linear relationship, and the hyperspectral observation values at the plurality of points ( Y) and the calculated band ratio (X) by calculating a and b of a regression equation (Y=aX+b), thereby generating a regression equation (Y=aX+b). Provides a method of utilizing hyperspectral images using rain.

Description

밴드비를 이용한 초분광영상 활용 방법{Method for utilizing of hyperspectral image using optimal band ratio}Method for utilizing of hyperspectral image using optimal band ratio}

본 발명은 밴드비(band ratio)를 이용한 초분광영상 활용 방법에 관한 것으로, 구체적으로, 관측값과 초분광영상을 이용하여 실재 물체를 가장 잘 대변하는 밴드비를 확인함으로써, 관측값이 없는 지점에서도 초분광영상만을 이용하여 높은 정확도로 물체를 식별할 수 있는 방법에 관한 것이다. The present invention relates to a method of utilizing a hyperspectral image using a band ratio, and specifically, a point without an observed value by confirming the band ratio that best represents an actual object using an observation value and a hyperspectral image. Also, it relates to a method that can identify objects with high accuracy using only hyperspectral images.

카메라에 포함되는 일반적인 이미지 센서(image sensor)는 가시광영역의 정보만을 감지하는데 반하여, 초분광 센서(hyperspectral sensor)를 사용할 경우 자외선영역, 적외선영역을 더 포함하여 수백배에 이르는 많은 정보를 확보할 수 있다. While a general image sensor included in a camera detects only information in the visible light region, using a hyperspectral sensor can secure hundreds of times more information including the ultraviolet region and the infrared region. have.

수많은 정보를 포함한 초분광정보를 이용할 경우, 해당 영상만으로도 촬영한 지점의 물체, 수심, 환경 등을 정밀하게 식별할 수 있다. 이를 위해, 실재의 파장별 밴드가 무엇을 의미하는지를 매칭시킨 데이터베이스가 필요하다. 예를 들어, 2010년에 발사된 천리안 해양관측 위성(Goestationary Ocean Color Imager, GOGI)에 의해여, 8개 밴드(412, 490, 512, 555, 660, 680, 765, 875nm)의 값을 알 경우, 각각의 값이 어떠한 해색(海色)을 의미하는지 여부에 대한 데이터베이스는 이미 확보되어 있는 상태이다. 이를 이용하여, 해양의 초분광영상을 촬영함으로써 지점별 초분광정보가 확인되면, 데이터베이스로부터 해당 파장에서의 값과 이에 매칭된 해색을 확인함으로써, 직접 해색을 관측하지 않아도, 초분광영상만으로 해색을 추정하게 된다.In the case of using hyperspectral information including a large number of information, the object, depth, and environment of the captured point can be accurately identified with only the corresponding image. To this end, a database is needed to match what the actual bands by wavelength mean. For example, if you know the value of 8 bands (412, 490, 512, 555, 660, 680, 765, 875nm) by the Goestationary Ocean Color Imager (GOGI) launched in 2010 The database on whether or not each value represents a sea color (海色) is already secured. Using this, when the hyperspectral information for each point is confirmed by photographing a hyperspectral image of the ocean, the value at the corresponding wavelength and the sea color matched thereto are checked from the database. Is estimated.

또한, 모든 지점에서의 파장별 관측값을 확인하고 데이터베이스화할 수 없기에, 회귀식이 사용된다. 즉, 다수의 데이터를 이용하여 회귀식을 구축하여 두면, 데이터베이스화되어 있지 않은 위치에서도, 회귀식을 이용하여, 초분광영상에서 확인되는 파장별 밴드를 이용하여 실재의 값을 추정할 수 있다. In addition, since observation values for each wavelength at all points cannot be confirmed and converted into a database, a regression equation is used. That is, if a regression equation is constructed using a plurality of data, the actual value can be estimated by using the band for each wavelength identified in the hyperspectral image by using the regression equation even at a location that is not in a database.

하지만, 이러한 방식으로 특정 파장의 촬영값(즉, 파장별 밴드)을 이용하여 정보를 추정하는 방식은 일정 수준의 오차가 있는 것으로 알려져 있다. 이러한 오차의 원인으로서, 초분광정보 촬영 단계에서 발생하는 오차(예를 들어, 반사광의 영향으로 촬영값이 실재 위치의 값을 대변하지 못하는 상황), 영상 수정 단계에서 발생하는 오차(예를 들어, 기하 보정의 오류, 영상 정합의 오류 등), 데이터베이스의 오류(예를 들어, 촬영값은 정확하였지만 식별 정보가 잘못된 경우), 기타 외부요인 등이 제시된다. However, it is known that the method of estimating information using photographing values of a specific wavelength (ie, bands for each wavelength) in this manner has a certain level of error. As a cause of such errors, errors occurring in the hyperspectral information photographing step (for example, a situation in which the photographed value does not represent the actual position value due to the influence of reflected light), and errors occurring in the image correction step (for example, Errors in geometric correction, errors in image registration, etc.), errors in the database (for example, when the photographed value is correct but the identification information is incorrect), and other external factors are presented.

더욱이, 정확도가 낮은 값을 이용하여 구축한 회귀식을 통해 값을 추정할 경우, 이러한 오차는 점차 커지게 된다.Moreover, when a value is estimated through a regression equation constructed using a value with low accuracy, this error gradually increases.

이에, 오차를 줄여서 정확도를 높일 수 있는 새로운 방식의 식별 방법에 대한 요구가 증가하는 추세이다. Accordingly, there is a growing demand for a new method of identification that can improve accuracy by reducing errors.

관련된 특허 문헌을 검토한다. Review related patent documents.

한국등록특허 제10-1463351호는 초분광영상에서 밴드비를 활용하여 특이영역을 검출하고, 초분광영상에서 검출된 영역을 제외함으로써 정보 연산의 속도를 상승시키는 기술을 개시한다. 후술하는 본 발명과 밴드비를 사용한다는 점만 공통될 뿐, 전체적인 발명의 목적, 구성 및 효과가 상이하다. Korean Patent Registration No. 10-1463351 discloses a technique for increasing the speed of information calculation by detecting a specific region in a hyperspectral image by using a band ratio and excluding the detected region in a hyperspectral image. The present invention described later and the use of a band ratio are common only, and the objects, configurations, and effects of the overall invention are different.

한국등록특허 제10-1619836호는 초분광영상 촬영을 위한 카메라를 짐벌과 완충장치를 이용한 드론을 개시한다. 드론의 흔들림을 보상한 영상이 확보된다는 점에서 장점이나, 라인 스캐닝에 의한 영상 정합 기술을 개시하지 못한다.Korean Patent Registration No. 10-1619836 discloses a drone using a gimbal and a shock absorber for a camera for shooting hyperspectral images. It is advantageous in that an image that compensates for the shaking of a drone is secured, but the image matching technology by line scanning cannot be disclosed.

한국등록특허 제10-1806488호는 초분광영상 센서를 드론에 탑재하여 암체를 식별하는 기술을 개시한다. 한국등록특허 제10-1744662호는 초분광영상을 이용하여 식물을 분류하는 기술을 개시한다. 이를 위해, 다양한 암체 내지 식물에 대한 초분광영상 표준자료가 미리 확보되어 있어야 하며, 촬영된 영상을 이와 비교함으로써, 피촬영 오브젝트를 식별한다. 하지만, 드론의 흔들림으로 인한 정확도 저하 문제, 영상 정합의 효율성을 상승시키는 문제 등 촬영 단계에서 확인되는 원천적인 문제를 해결하지 못한다. Korean Patent Registration No. 10-1806488 discloses a technology for identifying a dark body by mounting a hyperspectral image sensor on a drone. Korean Patent Registration No. 10-1744662 discloses a technology for classifying plants using hyperspectral images. To this end, standard data of hyperspectral images for various dark bodies or plants must be secured in advance, and the object to be photographed is identified by comparing the photographed images with them. However, it does not solve the fundamental problems found in the shooting stage, such as a problem of lowering accuracy due to the shaking of a drone, and a problem of increasing the efficiency of image matching.

한국등록특허 제10-1689197호는 드론에 부착될 수 있는 짐벌 어셈블리를 개시한다. 여기에 개시된 어셈블리는 회전/역회전이 가능한 제1유닛과, 승강이 가능한 제2유닛과, 이들을 회동시키는 제3유닛으로 구성되어, 비행체의 흔들림을 보상한다. 하지만, 이와 같이 구성될 경우, 짐벌의 무게가 과다하게 증가하여, 드론의 비행 성능이 낮아진다는 문제가 발생할 것이다.Korean Patent Registration No. 10-1689197 discloses a gimbal assembly that can be attached to a drone. The assembly disclosed herein is composed of a first unit capable of rotation/reverse rotation, a second unit capable of elevating and descending, and a third unit that rotates them, thereby compensating for shaking of the aircraft. However, when configured in this way, the weight of the gimbal will increase excessively, resulting in a problem that the flying performance of the drone is degraded.

한국공개특허 제10-2014-0014819호는 초분광영상에서 확인된 데이터에 대하여 색인/질의 기법을 적용함으로써, 마치 인공지능 알고리즘과 같은 형태로, 데이터를 식별하는 방법을 개시한다. 이러한 방법 활용의 전제조건은 초분광영상의 정확도인데, 무인비행체와 같은 환경에서 촬영된 영상의 정확도를 상승시키는 방법까지 제시하지 못한다. Korean Patent Publication No. 10-2014-0014819 discloses a method of identifying data in the form of an artificial intelligence algorithm by applying an index/query technique to data identified in a hyperspectral image. The prerequisite for using this method is the accuracy of hyperspectral images, but it cannot even suggest a method to increase the accuracy of images captured in an environment such as an unmanned aerial vehicle.

미국특허 제7,181,055호는 초분광영상을 정합하는 알고리즘을 제시한다. 여기에서는, 반사도 이미지를 활용하여 정합 정확도를 상승시킨다. US Patent No. 7,181,055 proposes an algorithm for matching hyperspectral images. Here, the reflectivity image is utilized to increase the matching accuracy.

(특허문헌 1) 한국등록특허 제10-1463351호 (Patent Document 1) Korean Patent Registration No. 10-1463351

(특허문헌 2) 한국등록특허 제10-1619836호 (Patent Document 2) Korean Patent Registration No. 10-1619836

(특허문헌 3) 한국등록특허 제10-1806488호(Patent Document 3) Korean Patent Registration No. 10-1806488

(특허문헌 4) 한국등록특허 제10-1744662호(Patent Document 4) Korean Patent Registration No. 10-1744662

(특허문헌 5) 한국등록특허 제10-1689197호(Patent Document 5) Korean Patent Registration No. 10-1689197

(특허문헌 6) 한국공개특허 제10-2014-0014819호(Patent Document 6) Korean Patent Publication No. 10-2014-0014819

(특허문헌 7) 미국특허 제7,181,055호(Patent Document 7) U.S. Patent No. 7,181,055

본 발명은 상기와 같은 문제점을 해결하기 위하여 안출된 것이다. The present invention was devised to solve the above problems.

구체적으로, 정확한 초분광정보를 획득하여 특정 파장이 갖고 있는 촬영값이 정확함에도, 해당 촬영값에 대한 데이터베이스가 정확한 식별 정보를 반영하지 못하는 상황을 방지하기 위한, 신규한 초분광정보 식별 방법을 제안하고자 한다. Specifically, a novel hyperspectral information identification method is proposed to prevent the situation in which the database for the corresponding photographed value does not reflect accurate identification information even though the photographed value of a specific wavelength is accurate by obtaining accurate hyperspectral information. I want to.

특히, 특정 파장의 촬영값을 사용하는 것이 아니라, 2개의 파장의 쌍의 촬영값인 밴드비를 이용하는 신규한 초분광영상 활용 방법을 제안하고자 한다.In particular, a novel hyperspectral image utilization method using a band ratio, which is an image value of a pair of two wavelengths, is proposed instead of using a photographing value of a specific wavelength.

상기와 같은 과제를 해결하기 위한 본 발명의 일 실시예는, (a) 초분광 관측값 입력모듈(310)에 미리 결정된 다수의 지점 각각의 초분광 관측값(Y)이 입력되는 단계; (b) 촬영영상 입력모듈(210)에 초분광정보를 갖는 촬영영상이 입력되고, 상기 미리 결정된 다수의 지점 각각에서 다수의 파장(x)에 대한 화소값인 초분광 촬영값(I(x))이 확인되는 단계; (c) 밴드비 연산모듈(320)이, 상기 초분광 촬영값에 포함된 다수의 파장(x) 중 임의의 2개의 파장의 쌍(λ1, λ2)을 선택하고, 각각에 대한 초분광 촬영값(I(λ1), I(λ2))을 확인하는 단계; (d) 상기 밴드비 연산모듈(320)이, 상기 확인된 2개의 초분광 촬영값(I(λ1), I(λ2))에서 어느 하나를 다른 하나로 나누어 밴드비(X)를 연산하는 단계; 및 (e) 회귀분석모듈(330)이, 상기 초분광 관측값(Y)과 상기 연산된 밴드비(X)가 선형관계에 있음을 가정하고, 상기 다수의 지점에서의 상기 초분광 관측값(Y)과 상기 연산된 밴드비(X)를 이용하여, 회귀식(Y=aX+b)의 a와 b를 연산함으로써, 회귀식(Y=aX+b)을 생성하는 단계를 포함하는, 밴드비를 이용한 초분광영상 활용 방법을 제공한다.An embodiment of the present invention for solving the above problems includes the steps of: (a) inputting a hyperspectral observation value (Y) of each of a plurality of predetermined points into the hyperspectral observation value input module 310; (b) A photographed image having hyperspectral information is input to the photographed image input module 210, and a hyperspectral photographing value (I(x)) that is a pixel value for a plurality of wavelengths (x) at each of the plurality of predetermined points ) Is confirmed; (c) The band ratio calculation module 320 selects a pair of arbitrary two wavelengths (λ 1 , λ 2 ) among a plurality of wavelengths (x) included in the hyperspectral photographing value, and hyperspectral for each Checking the photographed values I(λ 1 ) and I(λ 2 ); (d) The band ratio calculation module 320 calculates a band ratio (X) by dividing any one from the identified two hyperspectral photographing values (I(λ 1 ) and I(λ 2 )) into the other. step; And (e) the regression analysis module 330 assumes that the hyperspectral observation value (Y) and the calculated band ratio (X) are in a linear relationship, and the hyperspectral observation values at the plurality of points ( Y) and the calculated band ratio (X) by calculating a and b of a regression equation (Y=aX+b), thereby generating a regression equation (Y=aX+b). Provides a method of utilizing hyperspectral images using rain.

또한, 상기 (e) 단계 이후, (f) 상기 회귀분석모듈(330)이, 상기 초분광 촬영값에 포함된 상기 다수의 파장 중 상기 (c) 단계에서의 2개의 파장의 쌍을 제외한 다른 임의의 파장을 선택하여, 이에 대한 초분광 촬영값을 확인하고, 상기 (d) 내지 (e) 단계를 반복함으로써, 다수의 2개의 파장의 쌍마다 각각 회귀식을 생성하는 단계; (g) 회귀식 결정모듈(340)이, 상기 생성된 다수의 회귀식 중 상관계수(correlation)가 가장 높은 회귀식을 선택하는 단계; 및 (h) 상기 회귀식 결정모듈(340)이, 상기 (g) 단계에서 선택된 회귀식에 사용된 밴드비(X)에 적용된 2개의 파장을, 최적 밴드비 파장의 쌍으로 결정하는 단계를 포함하는 것이 바람직하다.In addition, after the step (e), (f) the regression analysis module 330, of the plurality of wavelengths included in the hyperspectral photographing value, other arbitrary other than the pair of two wavelengths in step (c) Generating a regression equation for each of a plurality of pairs of two wavelengths by selecting a wavelength of, confirming the hyperspectral photographing value for this, and repeating steps (d) to (e); (g) selecting, by the regression equation determination module 340, a regression equation having the highest correlation among the plurality of generated regression equations; And (h) determining, by the regression equation determination module 340, two wavelengths applied to the band ratio (X) used in the regression equation selected in step (g) as a pair of optimal band ratio wavelengths. It is desirable to do it.

또한, 상기 (h) 단계 이후, (i) 관측값 추정모듈(350)이, 상기 다수의 지점 이외의 지점의 초분광 촬영값에서, 상기 (h) 단계에서 결정된 최적 밴드비 파장의 쌍에 해당하는 2개의 초분광 촬영값을 확인하고, 어느 하나를 다른 하나로 나누어 해당 지점의 밴드비를 연산하는 단계; 및 (j) 상기 관측값 추정모듈(350)이, 상기 (g) 단계에서 선택된 회귀식에 상기 연산된 밴드비를 입력하여 해당 지점의 관측값을 추정하는 단계를 더 포함하는 것이 바람직하다.In addition, after step (h), (i) the observation value estimation module 350 corresponds to the pair of optimal band ratio wavelengths determined in step (h) from the hyperspectral photographing values of points other than the plurality of points. Checking the two hyperspectral photographing values, dividing one of them into the other, and calculating a band ratio of the corresponding point; And (j) the observation value estimating module 350, inputting the calculated band ratio to the regression equation selected in step (g) to estimate the observation value of the corresponding point.

또한, 상기 (g) 단계는, 상기 회귀식 결정모듈(340)이, 상기 생성된 다수의 회귀식에 대한 상관계수를 연산하여, 2개의 파장의 쌍에서 각각의 파장을 X축 및 Y축으로 하는 상관계수 맵을 생성하는 단계를 포함하고, 상기 (h) 단계는, 상기 회귀식 결정모듈(340)이, 상기 (g) 단계에서 생성된 상관계수 맵을 이용하여 상관계수가 가장 높은 지점을 결정하고, 해당 지점을 구성하는 각각의 파장을 상기 최적 밴드비 파장의 쌍으로 결정하는 단계를 포함하는 것이 바람직하다.In addition, in the step (g), the regression equation determination module 340 calculates correlation coefficients for the plurality of generated regression equations, and converts each wavelength from the pair of two wavelengths to the X axis and the Y axis. And generating a correlation coefficient map, wherein in step (h), the regression equation determination module 340 uses the correlation coefficient map generated in step (g) to determine the point with the highest correlation coefficient. It is preferable to include the step of determining and determining each wavelength constituting the corresponding point as a pair of the optimal band ratio wavelengths.

또한, 상기 (a) 단계는, (a1) 비행경로 설정모듈(220)이, 수집영역을 이용하여 무인비행체의 비행경로를 설정하는 단계; 및 (a2) 무인비행체가, 상기 설정된 비행경로를 따라 비행하면서, 초분광 센서를 이용하여, 라인 스캐닝을 통해, 기 설정된 영상 단위마다 상기 초분광정보를 획득하는 단계를 포함하는 것이 바람직하다.In addition, the (a) step, (a1) the flight path setting module 220, by using the collection area to set the flight path of the unmanned aerial vehicle; And (a2) acquiring the hyperspectral information for each preset image unit through line scanning by using a hyperspectral sensor while the unmanned aerial vehicle is flying along the set flight path.

상기와 같은 과제를 해결하기 위한 본 발명의 다른 실시예는, 컴퓨터 판독 가능한 기록 매체에 저장되어, 전술한 방법을 실행시키는 프로그램과, 그 코드가 기록된 컴퓨터 판독 가능한 기록 매체를 제공한다.Another embodiment of the present invention for solving the above problems provides a program that is stored in a computer-readable recording medium and executes the above-described method, and a computer-readable recording medium in which the code is recorded.

본 발명에 따른 방법에 의해, 밴드비를 활용함으로써, 초분광정보에서의 촬영값을 통하여, 실재 물체를 높은 정확도로 추정할 수 있다. According to the method according to the present invention, by utilizing the band ratio, the real object can be estimated with high accuracy through the photographed value in the hyperspectral information.

종래에는 실재 물체의 파장의 값과 가장 가까운 것으로 식별되는 값을 갖는 오직 1개의 파장만을 사용하기에, 해당 매칭이 부정확하거나 측정값에 오류가 있는 경우 실재 물체와 다른 물체를 추정하는 오류가 발생할 수 있다는 단점이 있었다. 본 발명은, 실재 물체의 파장의 값과 가장 가까운 것으로 식별되는 파장의 촬영값이 분자로 사용되고, 실재 물체의 파장의 값과 가장 먼 것으로 식별의 촬영값이 분모로 사용함으로써, 이중의 매칭 과정을 통해 실재 물체의 값을 추정하기에, 정확도를 높일 수 있는 것이다. Conventionally, only one wavelength with a value identified as the closest to the actual object's wavelength value is used, so if the matching is incorrect or there is an error in the measured value, an error in estimating an object different from the real object may occur. There was a drawback. In the present invention, the photographed value of the wavelength identified as the closest to the value of the wavelength of the real object is used as a numerator, and the photographed value of the identification as the farthest from the value of the wavelength of the real object is used as the denominator. Through the estimation of the value of the real object, the accuracy can be improved.

이러한 방법은, 컴퓨터 프로그램화되어, 다수의 지점의 관측값과 촬영된 초분광영상이 입력되면, 자동으로 회귀식이 결정될 수 있어서, 높은 정확도를 유지하면서도 신속한 회귀식 연산 및 이를 기반으로 한 대상 물체의 식별이 가능하다는 장점을 갖는다. This method is computer-programmed, and when observation values of multiple points and captured hyperspectral images are input, a regression equation can be automatically determined, so that a fast regression equation is calculated while maintaining high accuracy and It has the advantage of being identifiable.

도 1은 본 발명에 따른 방법을 수행하기 위한 시스템을 도시한다.
도 2는 본 발명에 따른 방법을 설명하기 위한 순서도이다.
도 3은 본 발명에 따른 방법의 밴드비를 설명하기 위한 도면이다.
도 4는 본 발명에 따른 방법에서 밴드비와 관측값을 이용하여 회귀분석이 수행되는 과정을 설명하기 위한 도면이다.
도 5는 본 발명에 따른 방법에 의하여 설정된 상관관계 맵을 도시한다.
도 6은 본 발명에 따른 방법에 의하여 상관관계 맵을 이용하여 회귀식을 선택하는 방법을 설명한다.
도 7은 본 발명에 따른 방법이 구현된 프로그램의 화면을 도시한다.1 shows a system for carrying out the method according to the invention.
2 is a flow chart for explaining the method according to the present invention.
3 is a diagram for explaining the band ratio of the method according to the present invention.
4 is a diagram illustrating a process of performing regression analysis using a band ratio and an observed value in the method according to the present invention.
5 shows a correlation map established by the method according to the present invention.
6 illustrates a method of selecting a regression equation using a correlation map by the method according to the present invention.
7 shows a screen of a program in which the method according to the present invention is implemented.

이하, 도면을 참조하여 본 발명을 보다 상세히 설명한다. Hereinafter, the present invention will be described in more detail with reference to the drawings.

시스템의 설명Description of the system

본 발명에 따른 방법이 수행되는 시스템은, 영상처리부(200)와 영상분석부(300)로 구분될 수 있다. 영상처리부(200)는 드론과 같은 무인비행체(미도시)에서 촬영한 촬영영상을 처리하는 기능을 수행하며, 영상분석부(300)는 여기에서 전처리된 영상을 이용하여 밴드비를 활용해 회귀분석을 수행하며, 관측값을 추정하는 기능을 수행한다.The system in which the method according to the present invention is performed may be divided into an image processing unit 200 and an image analysis unit 300. The image processing unit 200 performs a function of processing a photographed image captured by an unmanned aerial vehicle (not shown) such as a drone, and the image analysis unit 300 performs regression analysis using the band ratio using the image preprocessed here. And performs the function of estimating the observed value.

영상처리부(200)는 통신모듈(210), 비행경로 설정모듈(220), 촬영영상 입력모듈(230), 정규화 수행모듈(240), 기하보정 수행모듈(250), 영상정합 수행모듈(260) 및 초분광영상 출력모듈(270)을 포함한다. The image processing unit 200 includes a communication module 210, a flight path setting module 220, a photographed image input module 230, a normalization module 240, a geometric correction module 250, and an image matching module 260. And a hyperspectral image output module 270.

통신모듈(210)은 무인비행체(미도시) 와 유선 또는 무선으로 통신하여, 촬영영상 입력모듈(230)이 초분광센서(미도시)에서 획득되는 초분광정보를 수신하도록 한다. 수신하는 정보에는, 위치 정보가 포함될 수도 있다. The communication module 210 communicates with an unmanned aerial vehicle (not shown) by wire or wirelessly so that the photographed image input module 230 receives hyperspectral information obtained from a hyperspectral sensor (not shown). The information to be received may include location information.

비행경로 설정모듈(220)은 초분광정보를 수집하고자 하는 수집영역을 기초로 하여, 그 영상을 확보하기 위한 무인비행체(미도시)의 비행경로를 설정하는 기능을 수행한다. The flight path setting module 220 performs a function of setting a flight path of an unmanned aerial vehicle (not shown) to secure the image based on the collection area to which the hyperspectral information is to be collected.

정규화 수행모듈(240)은, 확보된 초분광정보의 정규화 기능을 수행한다. 모든 기하학적 변환에서, 회전, 크기, 평행성 보존 등은 원점(0, 0)을 기준으로 발생하기에, 각 영상에 존재하는 XY 좌표의 각 순서쌍의 정규화가 필요하다. The normalization execution module 240 performs a normalization function of the secured hyperspectral information. In all geometric transformations, rotation, size, and preservation of parallelism occur based on the origin (0, 0), so it is necessary to normalize each ordered pair of XY coordinates present in each image.

기하보정 수행모듈(250)은, 확보된 초분광정보에서 지면의 굴곡 등 기하학적 조건에 의한 왜곡을 보상하기 위한 기하보정을 수행한다.The geometric correction module 250 performs geometric correction for compensating for distortion caused by geometric conditions such as curvature of the ground from the obtained hyperspectral information.

영상정합 수행모듈(260)은, 확보된 다수의 초분광정보를 정합하는 기능을 수행한다. 이를 위하여, 적분영상 연산모듈(261) 및 대응점 확인모듈(262)이 구비된다.The image matching module 260 performs a function of matching a plurality of secured hyperspectral information. To this end, an integral image calculation module 261 and a correspondence point confirmation module 262 are provided.

초분광영상 출력모듈(270)은, 정규화되고, 기하보정되고, 영상정합된 최종 결과물로서 획득된 초분광영상을 출력하는 기능을 수행한다.The hyperspectral image output module 270 performs a function of outputting a hyperspectral image obtained as a normalized, geometrically corrected, and image-matched final result.

한편, 이와 같은 영상 처리 방식은 일례에 불과하며, 종래에 알려진 어떠한 처리 방식으로 초분광영상이 처리되어도 무방하다.Meanwhile, such an image processing method is only an example, and the hyperspectral image may be processed by any conventionally known processing method.

영상분석부(300)는, 초분광 관측값 입력모듈(310), 밴드비 연산모듈(320), 회귀분석모듈(330), 회귀식 결정모듈(340) 및 관측값 추정모듈(350)을 포함한다. The image analysis unit 300 includes a hyperspectral observation value input module 310, a band ratio calculation module 320, a regression analysis module 330, a regression equation determination module 340, and an observation value estimation module 350. do.

초분광 관측값 입력모듈(310)은, 미리 결정된 다수의 지점에 대한 실재 초분광 관측값을 입력하는 기능을 수행한다. 즉, 초분광영상에서 확보되는 촬영값과 상관관계를 확인하기 위한 실재의 관측값이 입력되는 것이다. The hyperspectral observation value input module 310 performs a function of inputting actual hyperspectral observation values for a plurality of predetermined points. That is, an actual observation value for confirming the correlation with the photographed value secured in the hyperspectral image is input.

밴드비 연산모듈(320)은, 영상처리부(200)로부터의 초분광영상에서 확보되는 촬영값(예를 들어, 화소값)을 통해 밴드비를 연산한다. The band ratio calculation module 320 calculates a band ratio using a photographing value (eg, a pixel value) secured from the hyperspectral image from the image processing unit 200.

회귀분석모듈(330)은, 관측값과 밴드비의 회귀분석을 수행하여 회귀식을 연산한다. The regression analysis module 330 calculates a regression equation by performing regression analysis of the observed value and the band ratio.

회귀식 결정모듈(340)은, 다수의 회귀식이 상관관계를 연산하여, 가장 상관성이 높은 회귀식을 결정하고, 이를 통해 최적의 밴드비를 결정하는 기능을 수행한다. The regression equation determination module 340 performs a function of determining a regression equation having the highest correlation by calculating a correlation between a plurality of regression equations, and determining an optimal band ratio through this.

관측값 추정모듈(350)은, 회귀식 결정모듈(340)에서 결정된 회귀식을 이용하여, 특정 초분광영상이 입력된 경우(즉, 임의의 밴드비가 입력된 경우), 이를 기반으로 관측값을 추정하고, 이를 토대로 실재 물체를 식별하는 기능을 수행한다.Observation value estimation module 350, using the regression equation determined in the regression equation determination module 340, when a specific hyperspectral image is input (that is, when an arbitrary band ratio is input), based on the observation value It performs the function of estimating and identifying real objects based on this.

각 모듈의 구체적인 기능 수행 방법은 후술한다.A method of performing specific functions of each module will be described later.

방법의 설명Description of the method

도 2 및 도 4 내지 도 7을 참조하여, 본 발명에 따른 초분광영상 활용 방법을 설명한다. A method of utilizing hyperspectral images according to the present invention will be described with reference to FIGS. 2 and 4 to 7.

먼저, 초분광영상이 확보되어 있어야 한다. 예를 들어, 영상처리부(200)의 비행 경로 설정모듈(220)이 수집영역을 이용하여 무인비행체(미도시)의 비행경로를 설정한 후, 무인비행체(미도시)가, 설정된 비행경로를 따라 비행하면서, 초분광 센서를 이용하여, 라인 스캐닝을 통해, 기 설정된 영상 단위마다 획득한 초분광영상일 수 있다. 이러한 초분광영상에는 초분광정보가 포함되어 있다.First, hyperspectral images must be secured. For example, after the flight path setting module 220 of the image processing unit 200 sets the flight path of the unmanned aerial vehicle (not shown) by using the collection area, the unmanned aerial vehicle (not shown) follows the set flight path. It may be a hyperspectral image acquired for each preset image unit through line scanning using a hyperspectral sensor while flying. This hyperspectral image includes hyperspectral information.

또한, 미리 결정된 다수의 지점마다 각각의 초분광 관측값(Y)이 확인되어 있어야 한다(S210). 즉, 초분광 관측값 입력모듈(310)에 미리 결정된 다수의 지점 각각의 초분광 관측값(Y)이 입력된 것이다. 여기서, 회귀분석의 정확도를 위하여, 관측값(Y)이 확인되는 지점은 다수일수록 바람직하다.In addition, each hyperspectral observation value (Y) must be confirmed for each of a plurality of predetermined points (S210). That is, the hyperspectral observation value Y of each of a plurality of predetermined points is input to the hyperspectral observation value input module 310. Here, for the accuracy of the regression analysis, it is preferable that the number of points at which the observed value Y is confirmed is larger.

다음, 촬영영상 입력모듈(210)에는 초분광정보를 갖는 촬영영상이 입력되고, 관측값이 입력되었던 미리 결정된 다수의 지점 각각에서 다수의 파장(x)에 대한 화소값인 초분광 촬영값(I(x))이 확인된다(S220). Next, a photographed image having hyperspectral information is input to the photographed image input module 210, and a hyperspectral photographing value (I), which is a pixel value for a plurality of wavelengths (x) at each of a plurality of predetermined points at which the observation value was input. (x)) is confirmed (S220).

다음, 밴드비 연산모듈(320)이, 초분광 촬영값에 포함된 다수의 파장(x) 중 임의의 2개의 파장의 쌍(λ1, λ2)을 선택하고, 각각에 대한 초분광 촬영값(I(λ1), I(λ2))을 확인한다(S230). 이 때에, 확인된 2개의 초분광 촬영값(I(λ1), I(λ2))에서 어느 하나를 다른 하나로 나누어 밴드비(X)를 연산한다(S240)(도 3 참조). 이를 수식으로 표현하면, 아래의 수학식 1과 같다. Next, the band ratio calculation module 320 selects a pair of arbitrary two wavelengths (λ 1 , λ 2 ) among a plurality of wavelengths (x) included in the hyperspectral photographing value, and hyperspectral photographing values for each (I(λ 1 ), I(λ 2 )) is checked (S230). At this time, the band ratio X is calculated by dividing any one of the identified two hyperspectral photographing values I(λ 1 ) and I(λ 2 ) by the other (S240) (see FIG. 3). If this is expressed by an equation, it is as shown in Equation 1 below.

Figure 112019000704320-pat00001
Figure 112019000704320-pat00001

달리 표현하면, 어느 하나의 지점마다 밴드비가 연산되어 초분광 관측값(Y)과 밴드비(X)가 확인되는데, 다수의 지점마다 각각 개별적으로 수행되는 것이다.In other words, the band ratio is calculated for each point to confirm the hyperspectral observation value (Y) and the band ratio (X), which are performed individually for each of a plurality of points.

다음, 회귀분석모듈(330)이, 초분광 관측값(Y)과 연산된 밴드비(X)가 선형관계에 있음을 가정하고, 다수의 지점 각각에서 S210 단계에서 확인된 초분광 관측값(Y)과 S240 단계에서 연산된 밴드비(X)를 이용함으로써, 회귀식(Y=aX+b)의 a와 b를 연산하여, 회귀식(Y=aX+b)을 생성한다(S250). 여기서 a를 연산하는 구체적인 수식은 아래의 수학식 2와 같으며, a의 값이 연산되면 b의 값은 당연히 연산된다. n은 관측값의 개수이며, X는 밴드비이고, Y는 관측값이다. 상부 바는 평균값을 의미한다(이는 이하의 모든 수식에서 동일하다).Next, the regression analysis module 330 assumes that the hyperspectral observation value (Y) and the calculated band ratio (X) are in a linear relationship, and the hyperspectral observation value (Y) confirmed in step S210 at each of a plurality of points ) And the band ratio (X) calculated in step S240, a and b of the regression equation (Y=aX+b) are calculated to generate a regression equation (Y=aX+b) (S250). Here, a specific formula for calculating a is the same as in Equation 2 below, and when the value of a is calculated, the value of b is naturally calculated. n is the number of observations, X is the band ratio, and Y is the observation. The top bar means the average value (this is the same for all formulas below).

Figure 112019000704320-pat00002
Figure 112019000704320-pat00002

도 4는 그 예시를 도시한다. 여기에는, 5개의 지점(P1, P2, P3, P4, P5)에서의 밴드비(X)와 관측값(Y)을 이용하여 연산된 회귀식이 우측 하단에 도시된다.Fig. 4 shows an example. Here, a regression equation calculated using the band ratio (X) and observed value (Y) at five points (P1, P2, P3, P4, P5) is shown in the lower right.

한편, 이러한 과정은 S230 단계에서 선택된 2개의 파장을 변경하면서 각각의 쌍에 대하여 이루어질 수 있다. 미리 결정된 지점이 동일하고, 각 지점에서의 관측값(Y)이 동일하더라도, 2개의 파장을 다르게 선택하면, 최종 연산되는 회귀식은 다르게 된다. Meanwhile, this process may be performed for each pair while changing the two wavelengths selected in step S230. Even if the predetermined points are the same and the observed values Y at each point are the same, if two wavelengths are selected differently, the final calculated regression equation is different.

구체적으로, 회귀분석모듈(330)이, 먼저 선택된 2개의 파장의 쌍을 제외한 다른 임의의 파장을 선택하여, 이에 대한 초분광 촬영값을 확인하고, S240 내지 S250 단계를 반복하여, 다수의 회귀식을 생성한다(S260).Specifically, the regression analysis module 330 first selects an arbitrary wavelength other than the selected pair of two wavelengths, checks the hyperspectral photographing value for this, and repeats steps S240 to S250, and a plurality of regression equations To generate (S260).

다음, 회귀식 결정모듈(340)이, 생성된 다수의 회귀식의 상관계수를 각각 연산하고, 그 중에서 상관계수가 가장 높은 회귀식을 선택한다(S270). 이 때에 선택된 회귀식에 사용된 밴드비(X)에 적용된 2개의 파장을, 최적 밴드비 파장의 쌍으로 결정된다. Next, the regression equation determination module 340 calculates correlation coefficients of the plurality of generated regression equations, respectively, and selects a regression equation having the highest correlation coefficient among them (S270). At this time, two wavelengths applied to the band ratio (X) used in the selected regression equation are determined as a pair of optimal band ratio wavelengths.

예를 들어, 초분광정보에 포함된 파장의 개수가 N개라면, 2개의 파장의 쌍의 개수는 그 조합(combination)에 의하여 NC2개이므로, 해당 개수만큼 회귀식과 상관계수가 연산되어, 최적의 회귀식과 밴드비와 그 파장이 연산될 것이다.For example, if the number of wavelengths included in the hyperspectral information is N, the number of pairs of two wavelengths is N C 2 due to the combination, so the regression equation and the correlation coefficient are calculated as much as the number, The optimal regression equation and band ratio and their wavelength will be calculated.

한편, 이러한 방식은 상관계수 맵을 이용하여 이루어질 수 있다. 상관계수 맵의 예시가 도 5에 도시된다. 여기에는, 회귀식 결정모듈(340)이, 생성된 다수의 회귀식에 대한 상관계수를 연산하되, 2개의 파장의 쌍에서 어느 하나의 파장을 X축으로 하고 다른 하나의 파장을 Y축으로 하여 생성한 상관계수 맵이 도시된다.Meanwhile, this method can be achieved using a correlation coefficient map. An example of a correlation coefficient map is shown in FIG. 5. Here, the regression equation determination module 340 calculates a correlation coefficient for a plurality of generated regression equations, but in the pair of two wavelengths, one wavelength is the X axis and the other wavelength is the Y axis. The generated correlation coefficient map is shown.

다음, 회귀식 결정모듈(340)은 상관계수 맵을 이용하여 상관계수가 가장 높은 지점을 결정하고, 해당 지점을 구성하는 각각의 파장을 상기 최적 밴드비 파장의 쌍으로 결정하게 된다.Next, the regression equation determination module 340 determines a point having the highest correlation coefficient using a correlation coefficient map, and determines each wavelength constituting the point as a pair of the optimal band ratio wavelengths.

도 6에는 이러한 방식을 자동화하는 개념이 도시된다. 분광특성인 파장마다의 촬영값을 추출하고, 이를 통해 각각의 밴드비를 산정하고, 관측값을 이용하여 회귀분석하면, 상관계수 맵이 생성되며, 각 파장의 쌍마다 상관계수를 연산함으로써, 최적 밴드비 파장과 최적 회귀식이 결정된다.In Figure 6 a concept of automating this approach is shown. By extracting the photographed value for each wavelength, which is a spectral characteristic, and calculating each band ratio through it, and by regression analysis using the observed value, a correlation coefficient map is generated, and by calculating the correlation coefficient for each pair of wavelengths, the optimum The band ratio wavelength and the optimal regression equation are determined.

도 7은 이를 구현한 프로그램의 출력 화면을 도시한다. 우측의 상관계수 맵이 도시되며, 그 아래에는 각 파장에 대한 상관계수 및 회귀식이 자동으로 연산됨을 확인할 수 있다. 7 shows an output screen of a program implementing this. The correlation coefficient map on the right is shown, and it can be seen that the correlation coefficient and regression equation for each wavelength are automatically calculated below it.

이와 같은 과정을 통하여, 최적의 회귀식과 최적의 밴드비가 결정되면, 관측값이 없는 위치에서도 초분광영상에 포함된 초분광정보의 파장과 촬영값만으로도 관측값을 추정할 수 있다. Through such a process, when the optimum regression equation and the optimum band ratio are determined, the observation value can be estimated using only the wavelength and photographing value of the hyperspectral information included in the hyperspectral image even at a location where there is no observation value.

예를 들어, 관측값 추정모듈(350)이, 초분광영상에 포함된 촬영값 중에서 최적 밴드비 파장의 쌍에 해당하는 2개의 초분광 촬영값을 확인하여 밴드비를 연산하면, 이를 최적의 회귀식에 입력하여 해당 지점의 관측값을 높은 정확도로 추정할 수 있는 것이다. For example, when the observation value estimation module 350 calculates a band ratio by checking two hyperspectral photographing values corresponding to a pair of optimal band ratio wavelengths among photographed values included in the hyperspectral image, it is optimally regressed. By entering the equation, you can estimate the observed value at that point with high accuracy.

이상, 본 명세서에는 본 발명을 당업자가 용이하게 이해하고 재현할 수 있도록 도면에 도시한 실시예를 참고로 설명되었으나 이는 예시적인 것에 불과하며, 당업자라면 본 발명의 실시예로부터 다양한 변형 및 균등한 타 실시예가 가능하다는 점을 이해할 것이다. 따라서 본 발명의 보호범위는 특허청구범위에 의해서 정해져야 할 것이다.As described above, the present specification has been described with reference to the embodiments shown in the drawings so that those skilled in the art can easily understand and reproduce the present invention, but these are merely exemplary, and those skilled in the art can make various modifications and equivalents from the embodiments of the present invention. It will be understood that embodiments are possible. Therefore, the protection scope of the present invention should be defined by the claims.

200: 영상처리부
210: 통신모듈
220: 비행경로 설정모듈
230: 촬영영상 입력모듈
240: 정규화 수행모듈
250: 기하보정 수행모듈
260: 영상정합 수행모듈
261: 적분영상 연산모듈
262: 대응점 확인모듈
270: 초분광영상 출력모듈
300: 영상분석부
310: 초분광 관측값 입력모듈
320: 밴드비 연산모듈
330: 회귀분석모듈
340: 회귀식 결정모듈
350: 관측값 추정모듈200: image processing unit
210: communication module
220: flight path setting module
230: photographed image input module
240: normalization execution module
250: geometric correction performance module
260: image matching module
261: integral image operation module
262: correspondence point confirmation module
270: hyperspectral image output module
300: image analysis unit
310: Hyperspectral observation value input module
320: band ratio calculation module
330: regression analysis module
340: regression determination module
350: observation value estimation module

Claims (7)

(a) 초분광 관측값 입력모듈(310)에 미리 결정된 다수의 지점 각각의 초분광 관측값(Y)이 입력되는 단계;
(b) 촬영영상 입력모듈(210)에 초분광정보를 갖는 촬영영상이 입력되고, 상기 미리 결정된 다수의 지점 각각에서 다수의 파장(x)에 대한 화소값인 초분광 촬영값(I(x))이 확인되는 단계;
(c) 밴드비 연산모듈(320)이, 상기 초분광 촬영값에 포함된 다수의 파장(x) 중 임의의 2개의 파장의 쌍(λ1, λ2)을 선택하고, 각각에 대한 초분광 촬영값(I(λ1), I(λ2))을 확인하는 단계;
(d) 상기 밴드비 연산모듈(320)이, 상기 확인된 2개의 초분광 촬영값(I(λ1), I(λ2))에서 어느 하나를 다른 하나로 나누어 밴드비(band ratio)(X)를 연산하는 단계; 및
(e) 회귀분석모듈(330)이, 상기 초분광 관측값(Y)과 상기 연산된 밴드비(X)가 선형관계에 있음을 가정하고, 상기 다수의 지점에서의 상기 초분광 관측값(Y)과 상기 연산된 밴드비(X)를 이용하여, 회귀식(Y=aX+b)의 a와 b를 연산함으로써, 회귀식(Y=aX+b)을 생성하는 단계를 포함하고,
상기 (e) 단계 이후,
(f) 상기 회귀분석모듈(330)이, 상기 초분광 촬영값에 포함된 상기 다수의 파장 중 상기 (c) 단계에서의 2개의 파장의 쌍을 제외한 다른 임의의 파장을 선택하여, 이에 대한 초분광 촬영값을 확인하고, 상기 (d) 내지 (e) 단계를 반복함으로써, 다수의 2개의 파장의 쌍마다 각각 회귀식을 생성하는 단계;
(g) 회귀식 결정모듈(340)이, 상기 생성된 다수의 회귀식 중 상관계수(correlation)가 가장 높은 회귀식을 선택하는 단계; 및
(h) 상기 회귀식 결정모듈(340)이, 상기 (g) 단계에서 선택된 회귀식에 사용된 밴드비(X)에 적용된 2개의 파장을, 최적 밴드비 파장의 쌍으로 결정하는 단계를 더 포함하는,
밴드비를 이용한 초분광영상 활용 방법.
(a) inputting a hyperspectral observation value Y of each of a plurality of predetermined points into the hyperspectral observation value input module 310;
(b) A photographed image having hyperspectral information is input to the photographed image input module 210, and a hyperspectral photographing value (I(x)) that is a pixel value for a plurality of wavelengths (x) at each of the plurality of predetermined points ) Is confirmed;
(c) The band ratio calculation module 320 selects a pair of arbitrary two wavelengths (λ 1 , λ 2 ) among a plurality of wavelengths (x) included in the hyperspectral photographing value, and hyperspectral for each Checking the photographed values I(λ 1 ) and I(λ 2 );
(d) The band ratio calculation module 320 divides any one of the identified two hyperspectral photographing values (I(λ 1 ) and I(λ 2 )) into the other and divides the band ratio (X Calculating ); And
(e) The regression analysis module 330 assumes that the hyperspectral observation value (Y) and the calculated band ratio (X) are in a linear relationship, and the hyperspectral observation value (Y) at the plurality of points. ) And the calculated band ratio (X) by calculating a and b of the regression equation (Y=aX+b), thereby generating a regression equation (Y=aX+b),
After step (e),
(f) The regression analysis module 330 selects an arbitrary wavelength other than the pair of two wavelengths in step (c) among the plurality of wavelengths included in the hyperspectral photographing value, and Generating a regression equation for each of a plurality of pairs of two wavelengths by checking the spectroscopic image value and repeating steps (d) to (e);
(g) selecting, by the regression equation determination module 340, a regression equation having the highest correlation among the plurality of generated regression equations; And
(h) further comprising the step of determining, by the regression equation determination module 340, two wavelengths applied to the band ratio (X) used in the regression equation selected in step (g) as a pair of optimal band ratio wavelengths doing,
How to use hyperspectral images using the band ratio.
삭제delete 제 1 항에 있어서,
상기 (h) 단계 이후,
(i) 관측값 추정모듈(350)이, 상기 다수의 지점 이외의 지점의 초분광 촬영값에서, 상기 (h) 단계에서 결정된 최적 밴드비 파장의 쌍에 해당하는 2개의 초분광 촬영값을 확인하고, 어느 하나를 다른 하나로 나누어 해당 지점의 밴드비를 연산하는 단계; 및
(j) 상기 관측값 추정모듈(350)이, 상기 (g) 단계에서 선택된 회귀식에 상기 연산된 밴드비를 입력하여 해당 지점의 관측값을 추정하는 단계를 더 포함하는,
밴드비를 이용한 초분광영상 활용 방법.
The method of claim 1,
After step (h),
(i) The observation value estimation module 350 checks two hyperspectral photographing values corresponding to the pair of optimal band ratio wavelengths determined in step (h) from the hyperspectral photographing values of points other than the plurality of points. And calculating a band ratio of a corresponding point by dividing one of them into the other; And
(j) the observation value estimation module 350, the step of estimating the observation value of the corresponding point by inputting the calculated band ratio to the regression equation selected in step (g),
How to use hyperspectral images using the band ratio.
제 3 항에 있어서,
상기 (g) 단계는,
상기 회귀식 결정모듈(340)이, 상기 생성된 다수의 회귀식에 대한 상관계수를 연산하여, 2개의 파장의 쌍에서 각각의 파장을 X축 및 Y축으로 하는 상관계수 맵을 생성하는 단계를 포함하고,
상기 (h) 단계는,
상기 회귀식 결정모듈(340)이, 상기 (g) 단계에서 생성된 상관계수 맵을 이용하여 상관계수가 가장 높은 지점을 결정하고, 해당 지점을 구성하는 각각의 파장을 상기 최적 밴드비 파장의 쌍으로 결정하는 단계를 포함하는,
밴드비를 이용한 초분광영상 활용 방법.
The method of claim 3,
The step (g),
The regression equation determination module 340 calculates correlation coefficients for the plurality of generated regression equations, and generates a correlation coefficient map in which each wavelength is an X axis and a Y axis in the pair of two wavelengths. Including,
The step (h),
The regression equation determination module 340 determines a point with the highest correlation coefficient using the correlation coefficient map generated in step (g), and sets each wavelength constituting the point to the optimal band ratio wavelength pair. Including the step of determining as,
How to use hyperspectral images using the band ratio.
제 1 항에 있어서,
상기 (a) 단계는,
(a1) 비행경로 설정모듈(220)이, 수집영역을 이용하여 무인비행체의 비행경로를 설정하는 단계; 및
(a2) 무인비행체가, 상기 설정된 비행경로를 따라 비행하면서, 초분광 센서를 이용하여, 라인 스캐닝을 통해, 기 설정된 영상 단위마다 상기 초분광정보를 획득하는 단계를 포함하는,
밴드비를 이용한 초분광영상 활용 방법.
The method of claim 1,
The step (a),
(a1) the flight path setting module 220, setting the flight path of the unmanned aerial vehicle using the collection area; And
(a2) including the step of acquiring the hyperspectral information for each preset image unit through line scanning by using a hyperspectral sensor while flying by the unmanned aerial vehicle along the set flight path,
How to use hyperspectral images using the band ratio.
컴퓨터 판독 가능한 기록 매체에 저장되어
제 1 항 및 제 3 항 내지 제 5 항 중 어느 한 항에 따른 방법을 실행시키는 프로그램.
Stored on a computer-readable recording medium
A program for executing the method according to any one of claims 1 and 3 to 5.
제 1 항 및 제 3 항 내지 제 5 항 중 어느 한 항에 따른 방법을 실행하기 위한 프로그램 코드가 기록된 컴퓨터 판독 가능한 기록 매체. A computer-readable recording medium on which a program code for executing the method according to any one of claims 1 and 3 to 5 is recorded.
KR1020190000699A 2019-01-03 2019-01-03 Method for utilizing of hyperspectral image using optimal band ratio KR102142616B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020190000699A KR102142616B1 (en) 2019-01-03 2019-01-03 Method for utilizing of hyperspectral image using optimal band ratio

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020190000699A KR102142616B1 (en) 2019-01-03 2019-01-03 Method for utilizing of hyperspectral image using optimal band ratio

Publications (2)

Publication Number Publication Date
KR20200084973A KR20200084973A (en) 2020-07-14
KR102142616B1 true KR102142616B1 (en) 2020-08-10

Family

ID=71526579

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020190000699A KR102142616B1 (en) 2019-01-03 2019-01-03 Method for utilizing of hyperspectral image using optimal band ratio

Country Status (1)

Country Link
KR (1) KR102142616B1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102504280B1 (en) * 2020-11-13 2023-02-27 단국대학교 산학협력단 Image Registration Technique of Hyperspectral Image using Optical Flow Algorithm
KR102302745B1 (en) * 2021-07-23 2021-09-15 국방과학연구소 Method of acquiring image data, apparatus for acquiring image data, and computer program for the method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101463354B1 (en) 2013-10-22 2014-11-18 영남대학교 산학협력단 Apparatus for target detection of hyperspectral image and method thereof
CN106053370A (en) 2016-08-02 2016-10-26 山东省科学院海洋仪器仪表研究所 Inversion method for offshore secchi disk depth based on HICO simulation
KR101756960B1 (en) * 2016-05-31 2017-07-13 (주)아세아항측 Monitoring method of heavy metals in soil using the aerial hyperspectral images

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101463354B1 (en) 2013-10-22 2014-11-18 영남대학교 산학협력단 Apparatus for target detection of hyperspectral image and method thereof
KR101756960B1 (en) * 2016-05-31 2017-07-13 (주)아세아항측 Monitoring method of heavy metals in soil using the aerial hyperspectral images
CN106053370A (en) 2016-08-02 2016-10-26 山东省科学院海洋仪器仪表研究所 Inversion method for offshore secchi disk depth based on HICO simulation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
김태우 외 2, ‘수질 매개변수 추정에 있어서 항공 초분광영상의 가용성 고찰’ (대한원격탐사학회지, 제30권, 2013.11.20.) pp. 61~73.

Also Published As

Publication number Publication date
KR20200084973A (en) 2020-07-14

Similar Documents

Publication Publication Date Title
KR102313438B1 (en) Method for estimation of river bed change rate using hyperspectral image
KR102169687B1 (en) Method for acquisition of hyperspectral image using an unmanned aerial vehicle
JP6729791B2 (en) Vegetation index calculation device, vegetation index calculation method, and program
KR102142616B1 (en) Method for utilizing of hyperspectral image using optimal band ratio
WO2017149869A1 (en) Information processing device, method, program, and multi-camera system
CN110998657B (en) Image processing apparatus, image processing method, and program
KR102481896B1 (en) System and method for establishing structural exterior map using image stitching
JP2017527785A (en) Method and system for photogrammetric processing of images
JP7379299B2 (en) Position and orientation estimation device, position and orientation estimation method, and program
US20200051261A1 (en) Image processing device and image capture apparatus
WO2021108626A1 (en) System and method for correspondence map determination
WO2022213218A1 (en) System and method for vegetation detection from aerial photogrammetric multispectral data
Wierzbicki et al. Method of radiometric quality assessment of NIR images acquired with a custom sensor mounted on an unmanned aerial vehicle
KR102164522B1 (en) Method for estimation of depth of water using hyperspectral image
CN113326807A (en) Black and odorous water body identification method and device and electronic equipment
JP6201507B2 (en) Image processing apparatus, image processing method, and image processing program
KR102177985B1 (en) Apparatus and Method for Automatically Detecting Landslide Area
KR102295630B1 (en) Method for utilizing of hyperspectral image using optimal band index
JP5807111B2 (en) Image generation system and image generation method
WO2020179276A1 (en) Image processing device, image processing method, and program
KR102295631B1 (en) Method for utilizing of hyperspectral image using optimal spectral characteristic index
JP7451465B2 (en) Learning method, program and image processing device
CN112785693B (en) Method, system and device for generating intrinsic hyperspectral point cloud
WO2016189853A1 (en) Image processing device, image processing system, image processing method, and program recording medium
CN114511781A (en) Method, device and medium for identifying disguise personnel by carrying multispectral camera by unmanned aerial vehicle

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant