KR20110060626A - The method for measurign object's velocity using synthetic aperture radar image and the apparatus thereof - Google Patents

The method for measurign object's velocity using synthetic aperture radar image and the apparatus thereof Download PDF

Info

Publication number
KR20110060626A
KR20110060626A KR1020090117256A KR20090117256A KR20110060626A KR 20110060626 A KR20110060626 A KR 20110060626A KR 1020090117256 A KR1020090117256 A KR 1020090117256A KR 20090117256 A KR20090117256 A KR 20090117256A KR 20110060626 A KR20110060626 A KR 20110060626A
Authority
KR
South Korea
Prior art keywords
speed
interference
image
observation
mai
Prior art date
Application number
KR1020090117256A
Other languages
Korean (ko)
Other versions
KR101135070B1 (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 KR1020090117256A priority Critical patent/KR101135070B1/en
Publication of KR20110060626A publication Critical patent/KR20110060626A/en
Application granted granted Critical
Publication of KR101135070B1 publication Critical patent/KR101135070B1/en

Links

Images

Classifications

    • 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P3/00Measuring linear or angular speed; Measuring differences of linear or angular speeds
    • G01P3/42Devices characterised by the use of electric or magnetic means
    • G01P3/50Devices characterised by the use of electric or magnetic means for measuring linear speed
    • 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/50Systems of measurement based on relative movement of target
    • G01S13/52Discriminating between fixed and moving objects or between objects moving at different speeds
    • G01S13/522Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves
    • G01S13/524Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi
    • G01S13/53Discriminating between fixed and moving objects or between objects moving at different speeds using transmissions of interrupted pulse modulated waves based upon the phase or frequency shift resulting from movement of objects, with reference to the transmitted signals, e.g. coherent MTi performing filtering on a single spectral line and associated with one or more range gates with a phase detector or a frequency mixer to extract the Doppler information, e.g. pulse Doppler radar
    • 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/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

PURPOSE: A method and an apparatus for measuring the speed of an object using synthetic aperture radar images are provided to determine two dimensional speed of the object by observing the speeds of the surface observing direction and the aviating direction of the object. CONSTITUTION: A first interference diagram forming part(110) forms an along track SAR interferometry(ATI) interference diagram based on synthetic aperture radar images. A second interference diagram forming part(120) forms a multiple aperture SAR interferometry(MAI) interference diagram based on the synthetic aperture radar images. A first speed measuring part(130) measures the speed of the surface observing direction with respect to an object based on the ATI interference diagram. A second speed measuring part(140) measures the speed of the aviating direction with respect to the object based on the MAI interference diagram. A speed determining part(150) combines the speed of the surface observing direction and the aviating direction to determine the two dimensional speed of the object.

Description

합성 개구 레이더 영상을 이용한 물체의 속도측정 방법 및 그 장치{THE METHOD FOR MEASURIGN OBJECT'S VELOCITY USING SYNTHETIC APERTURE RADAR IMAGE AND THE APPARATUS THEREOF}TECHNICAL FIELD METHOD AND DEVICE MEASUREMENT USING A Synthetic Aperture Radar Image {THE METHOD FOR MEASURIGN OBJECT'S VELOCITY USING SYNTHETIC APERTURE RADAR IMAGE AND THE APPARATUS THEREOF}

본 발명은 물체의 속도 측정 방법 및 장치에 관한 것으로, 보다 상세하게는 합성 개구 레이더 영상을 이용하여 물체의 속도 측정 방법 및 장치에 관한 것이다. The present invention relates to a method and apparatus for measuring the speed of an object, and more particularly, to a method and apparatus for measuring the speed of an object using a synthetic aperture radar image.

ATI 기법은 비행방향으로 레이더 안테나 두 개를 탑재한 비행기 또는 위성으로부터 한 쌍의 합성 개구 레이더 영상 자료를 취득하고, 이들 자료로부터 레이더 간섭영상을 제작하며, 레이더 간섭영상의 위상값을 시선(ling-of-sight; LOS)방향으로의 속도로 변환한 후, 최종적으로 지표관측방향으로의 이동물체의 속도를 관측하는 방법이다. 여기서 두 개의 레이더 안테나 중 하나는 송수신기 기능을 지니고 있고, 다른 하나는 수신기의 기능만을 지니고 있다. The ATI technique acquires a pair of synthetic aperture radar image data from an airplane or satellite equipped with two radar antennas in the direction of flight, produces radar interference images from these data, and visualizes the phase values of radar interference images. After converting to the velocity in the direction of-sight (LOS), the velocity of the moving object in the surface observation direction is finally observed. Here, one of the two radar antennas has a transceiver function, and the other has only the function of a receiver.

이와 같은 ATI 기법은 지표관측방향으로의 이동물체의 속도를 매우 정밀하게 관측할 수 있어 자동차, 배의 속도 관측, 해파의 속도관측, 빙하속도 관측 등 매우 다양한 분야에서 이용되고 있다. The ATI technique can accurately observe the speed of a moving object in the surface observation direction and is used in various fields such as automobile, ship speed observation, sea wave speed observation and glacier speed observation.

그러나 이와 같은 ATI 기법은 레이더의 지표관측방향(비행방향의 수직방향) 으로의 이동속도만을 관측할 수 있는 단점을 지닌다. 그렇기 때문에 비행방향으로 움직이는 이동물체의 속도를 관측할 수 없었다. 이와 같은 최대 단점을 극복하기 위하여 자동차의 경우 도로를 따라 움직이기 때문에 도로의 방향을 이용하여 측정하고 있지만, 이 방법은 도로가 비행방향에 가까울수록 그 오차는 기하급수적으로 커지는 문제를 지니고 있다. 더욱이 자유로운 움직임을 보이는 해수의 이동, 배의 이동, 빙하의 이동 등에 관하여 풀리지 않는 숙제로 남아있다. However, this ATI technique has the disadvantage of observing only the moving speed of the radar in the surface observation direction (the vertical direction of the flight direction). Therefore, the speed of the moving object moving in the flight direction could not be observed. In order to overcome the biggest drawback, the vehicle is measured using the direction of the road because it moves along the road, but this method has a problem that the error increases exponentially as the road is closer to the flight direction. Moreover, it remains an unresolved assignment on the movement of seawater, the movement of ships, and the movement of glaciers, which show free movement.

본 발명은 상기한 문제점을 해결하기 위해 합성 개구 레이더 영상 자료를 이용하여 ATI 간섭도로부터 지표관측방향의 속도를 관측하고, MAI 간섭도로부터 비행방향의 속도를 관측하여 이들을 조합하여 2차원 속도를 추출하는 방법을 제공하는데 그 목적이 있다. In order to solve the above problems, the present invention observes the velocity of the surface observation direction from the ATI interference using synthetic aperture radar image data, extracts the two-dimensional velocity by observing the velocity of the flight direction from the MAI interference and combines them. The purpose is to provide a way to.

상기 목적을 달성하기 위한 본 발명에 따른, 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법은 합성 개구 레이더 영상을 이용하여 ATI(along-track SAR interferometry)간섭도를 제작하는 제1 제작 단계; 상기 ATI 간섭도로부터 이동물체에 대한 지표관측방향의 속도를 관측하는 제1 관측 단계; 상기 합성 개구 레이더 영상 자료를 이용하여 MAI(multiple aperture SAR interferometry)간섭도를 제작하는 제2 제작 단계; 상기 MAI 간섭도로부터 상기 이동물체에 대한 비행방향의 속도를 관측하는 제2 관측 단계; 및 상기 지표관측방향의 속도 및 상기 비행방향의 속도를 결합하여 상기 이동물체의 2차원 속도를 결정하는 단계;를 포함한다.According to an aspect of the present invention, there is provided a method of measuring velocity of an object using a synthetic aperture radar image, comprising: a first manufacturing step of producing an ATI (along-track SAR interferometry) interference using the synthetic aperture radar image; A first observation step of observing a velocity in a surface observation direction with respect to a moving object from the ATI interference degree; A second manufacturing step of manufacturing multiple aperture SAR interferometry (MAI) coherence using the synthesized aperture radar image data; A second observation step of observing a velocity of a flying direction with respect to the moving object from the MAI interference degree; And determining the two-dimensional velocity of the moving object by combining the velocity in the surface observation direction and the velocity in the flight direction.

그리고, 상기 제1 제작 단계는, 레이더 센서를 앞 및 뒤에 구비한 한 쌍의 센서탑재체에서 한 쌍의 전방관측 SLC 영상과 한 쌍의 후방관측 SLC 영상을 이용하여 ATI 간섭도를 제작하는 것이 바람직하다.And, in the first manufacturing step, it is preferable to produce the ATI interference degree using a pair of front observation SLC image and a pair of rear observation SLC image in a pair of sensor mounting body having a radar sensor in front and rear. .

또한, 상기 ATI 간섭도는, 상기 전방관측 SLC 영상과 상기 후방관측 SLC 영상의 컬레복소수를 곱함으로써 결정되는 것이 바람직하다.In addition, the ATI interference degree is preferably determined by multiplying the complex complex number of the front observation SLC image and the rear observation SLC image.

그리고, 상기 제1 관측 단계는, 상기 ATI 간섭도로부터 ATI 간섭도의 위상을 결정하는 단계; 상기 ATI 간섭도의 위상, 레이더의 파장, 상기 센서탑재체의 파장 및 상기 앞의 레이더 센서와 상기 뒤의 레이더 센서간의 거리를 이용하여 상기 이동물체에 대한 시선방향의 속도를 결정하는 단계; 및 상기 시선 방향의 속도 및 상기 레이더의 입사각을 이용하여 상기 이동물체에 대한 지표관측 방향의 속도를 관측하는 단계;를 포함하는 것이 바람직하다.The first observation may include determining a phase of an ATI interference degree from the ATI interference degree; Determining a speed in the visual direction with respect to the moving object using the phase of the ATI interference degree, the wavelength of the radar, the wavelength of the sensor mounting body, and the distance between the front radar sensor and the rear radar sensor; And observing the speed in the surface observation direction with respect to the moving object using the speed in the eye direction and the incident angle of the radar.

또한, 상기 제2 제작 단계는, 레이더 센서를 앞 및 뒤에 구비한 한 쌍의 센서탑재체에서 한 쌍의 전방관측(forward-looking) SLC영상과 한 쌍의 후방관측(backward-looking) SLC 영상을 제작하고, 상기 전방관측 SLC 영상 쌍으로부터 전방관측 InSAR 간섭도와 상기 후방관측 SLC 영상 쌍으로부터 후방관측 InSAR 간섭도를 제작하며, 상기 전방관측 InSAR 간섭도 및 상기 후방관측 InSAR 간섭도로부터 MAI 간섭도를 제작하는 것이 바람직하다.In the second manufacturing step, a pair of forward-looking SLC images and a pair of backward-looking SLC images are fabricated from a pair of sensor mounting bodies equipped with radar sensors in front and behind. And producing a front view InSAR interference from the front view SLC image pair and a rear view InSAR interference from the rear view SLC image pair, and producing a MAI interference from the front view InSAR interference and the rear view InSAR interference. It is preferable.

그리고, 상기 한쌍의 전방 관측 SLC 영상 중 하나를 주영상 및 나머지를 부영상으로 결정하고, 상기 한 쌍의 후방 관측 SLC 영상 중 하나를 주영상 및 나머지를 부영상으로 설정하며, 상기 부영상의 컬러복소수를 상기 부영상에 대응되는 주영상에 곱함으로써 상기 전방관측 InSAR 간섭도 및 상기 후방관측 InSAR 간섭도를 제작하며, 상기 후방관측 InSAR 간섭도의 컬레복소수를 상기 전방관측 InSAR 간섭도에 곱함으로써 상기 MAI 간섭도를 제작하는 것이 바람직하다.One of the pair of forward-viewing SLC images is determined as the main image and the rest as sub-pictures, and one of the pair of backward-viewing SLC images is set as the main-image and the rest as sub-pictures, and the color complex number of the sub-images is determined. Multiplying the main image corresponding to the sub-image to produce the front-view InSAR interference and the back-view InSAR interference, and multiply the COLE complex number of the back-view InSAR interference by the front-view InSAR interference. It is desirable to produce an interference degree.

또한, 상기 MAI 간섭도로부터 MAI 간섭도의 위상을 결정하는 단계; 및 상기 MAI 간섭도에 대한 위상, 두 레이더 센서 탑재체가 동일 지역을 촬영하는 시간간 격, 비행방향으로의 레이더 안테나 길이 및 정규화된 스퀸스 조정 변수(normalized squint adjustment parameter)를 이용하여 상기 비행방향의 속도를 결정하는 단계;를 포함하는 것이 바람직하다.In addition, determining the phase of the MAI interference from the MAI interference; And the phase of the MAI interference, the time interval for the two radar sensor payloads to photograph the same area, the radar antenna length in the flight direction, and the normalized squint adjustment parameter. Determining a speed; preferably.

그리고, 상기 MAI 간섭도의 위상을 결정하는 단계는, 상기 이동물체가 지표관측 방향으로 이동하고 있는 경우, 상기 지표관측방향으로의 속도를 이용하여 상기 MAI 간섭도에 대한 위상을 보정하기 위한 보정값을 결정하는 단계; 및 상기 MAI 간섭도에서 상기 보정값의 차를 최종 MAI 간섭도로 결정하는 단계;를 포함하는 것이 바람직하다.The determining of the phase of the MAI interference may include: a correction value for correcting the phase of the MAI interference using the speed in the surface observation direction when the moving object moves in the surface observation direction. Determining; And determining the difference between the correction values in the MAI interference degree as a final MAI interference degree.

또한, 상기 보정값은, 비행방향으로의 레이더 안테나 길이, 상기 지표관측 방향으로의 속도, 정규화된 스퀸스 조정 변수(normalized squint adjustment parameter), 기선의 수평면에 투영된 길이 및 레이더 센서 탑재체의 속도를 이용하여 결정되는 것이 바람직하다.Further, the correction value may include the radar antenna length in the flight direction, the speed in the surface observation direction, the normalized squint adjustment parameter, the length projected on the horizontal plane of the base line, and the speed of the radar sensor payload. It is preferable to determine using.

한편, 상기 목적을 달성하기 위한 본 발명에 따른, 합성 개구 레이더 영상을 이용한 물체의 속도 측정 장치는, 합성 개구 레이더 영상 자료를 이용하여 ATI(along-track SAR interferometry)간섭도를 제작하는 제1 간섭도 제작부(110); 상기 ATI 간섭도로부터 이동물체에 대한 지표관측방향의 속도를 관측하는 제1 속도 관측부; 상기 합성 개구 레이더 영상 자료를 이용하여 MAI(multiple aperture SAR interferometry)간섭도를 제작하는 제2 간섭도 제작부; 상기 MAI 간섭도로부터 상기 이동물체에 대한 비행방향의 속도를 관측하는 제2 속도 관측부; 및 상기 지표관측방향의 속도 및 상기 비행방향의 속도를 결합하여 상기 이동물체의 2차원 속도를 결정하는 속도 결정부;를 포함한다.On the other hand, according to the present invention for achieving the above object, the apparatus for measuring the velocity of an object using a synthetic aperture radar image, the first interference to produce the ATI (along-track SAR interferometry) interference using the synthetic aperture radar image data Road production unit 110; A first speed observing unit for observing the speed in the surface observation direction with respect to a moving object from the ATI interference degree; A second interference producing unit configured to manufacture multiple aperture SAR interferometry (MAI) interference using the synthesized aperture radar image data; A second speed observing unit for observing the speed of the flying direction with respect to the moving object from the MAI interference degree; And a speed determiner configured to determine the two-dimensional speed of the moving object by combining the speed in the surface observation direction and the speed in the flight direction.

그리고, 상기 제1 간섭도 제작부는, 레이더 센서를 앞 및 뒤에 구비한 한 쌍의 센서탑재체에서 한 쌍의 전방관측 SLC 영상과 한 쌍의 후방관측 SLC 영상을 이용하여 ATI 간섭도를 제작하는 것이 바람직하다.In addition, the first interference manufacturing unit, it is preferable to produce an ATI interference degree using a pair of front observation SLC image and a pair of rear observation SLC image in a pair of sensor mounting body having the radar sensor front and back. Do.

또한, 상기 ATI 간섭도는 상기 전방관측 SLC 영상과 상기 후방관측 SLC 영상의 컬레복소수를 곱함으로써 결정되는 것이 바람직하다.In addition, the ATI interference degree is preferably determined by multiplying the complex complex number of the front observation SLC image and the rear observation SLC image.

그리고, 상기 제1 속도 관측부는, 상기 ATI 간섭도로부터 ATI 간섭도의 위상을 결정하고, 상기 ATI 간섭도의 위상, 레이더의 파장, 상기 센서탑재체의 파장 및 상기 앞의 레이더 센서와 상기 뒤의 레이더 센서간의 거리를 이용하여 상기 이동물체에 대한 시선방향의 속도를 결정하며, 상기 시선 방향의 속도 및 상기 레이더의 입사각을 이용하여 상기 이동물체에 대한 지표관측 방향의 속도를 관측하는 것이 바람직하다.The first velocity observation unit determines the phase of the ATI interference degree from the ATI interference degree, and determines the phase of the ATI interference degree, the wavelength of the radar, the wavelength of the sensor mounting body, the front radar sensor and the rear radar. It is preferable to determine the speed in the eye direction with respect to the moving object using the distance between the sensors, and to observe the speed in the surface observation direction with respect to the moving object using the speed in the eye direction and the incident angle of the radar.

또한, 상기 제2 간섭도 제작부는, 레이더 센서를 앞 및 뒤에 구비한 한 쌍의 센서탑재체에서 한 쌍의 전방관측(forward-looking) SLC영상과 한 쌍의 후방관측(backward-looking) SLC 영상을 제작하고, 상기 전방관측 SLC 영상 쌍으로부터 전방관측 InSAR 간섭도와 상기 후방관측 SLC 영상 쌍으로부터 후방관측 InSAR 간섭도를 제작하며, 상기 전방관측 InSAR 간섭도 및 상기 후방관측 InSAR 간섭도로부터 MAI 간섭도를 제작하는 것이 바람직하다.In addition, the second interference degree production unit, a pair of forward-looking SLC image and a pair of backward-looking SLC image in a pair of sensor mounting body having the radar sensor front and back Fabricate and produce a front view InSAR interference from the front view SLC image pair and a rear view InSAR interference from the rear view SLC image pair, and produce a MAI interference from the front view InSAR interference and the rear view InSAR interference. It is desirable to.

그리고, 상기 제2 간섭도 제작부는, 상기 한쌍의 전방 관측 SLC 영상 중 하나를 주영상 및 나머지를 부영상으로 결정하고, 상기 한 쌍의 후방 관측 SLC 영상 중 하나를 주영상 및 나머지를 부영상으로 설정하며, 상기 부영상의 컬러복소수를 상기 부영상에 대응되는 주영상에 곱함으로써 상기 전방관측 InSAR 간섭도 및 상기 후방관측 InSAR 간섭도를 제작하며, 상기 후방관측 InSAR 간섭도의 컬레복소수를 상기 전방관측 InSAR 간섭도에 곱함으로써 상기 MAI 간섭도를 제작하는 것이 바람직하다.The second interference degree producing unit may determine one of the pair of forward-viewing SLC images as a main image and the other as a sub-image, and one of the pair of backward-viewing SLC images as a main image and the rest as a sub-image. Set the color complex number of the sub-image to the main image corresponding to the sub-image to produce the front observation InSAR interference degree and the rear observation InSAR interference degree, and the front complex observation of the rear complex InSAR interference degree. It is preferable to produce the MAI interference level by multiplying the InSAR interference degree.

또한, 상기 제2 속도 관측부는, 상기 MAI 간섭도로부터 MAI 간섭도의 위상을 결정하며, 상기 MAI 간섭도의 위상, 두 레이더 센서 탑재체가 동일 지역을 촬영하는 시간간격, 비행방향으로의 레이더 안테나 길이 및 정규화된 스퀸스 조정 변수(normalized squint adjustment parameter)를 이용하여 상기 비행방향의 속도를 관측하는 것이 바람직하다.In addition, the second speed observation unit, the phase of the MAI interference from the MAI interference degree, the phase of the MAI interference degree, the time interval between the two radar sensor payload to shoot the same area, the radar antenna length in the flight direction And using the normalized squint adjustment parameter to observe the velocity of the flight direction.

그리고, 상기 제2 속도 관측부는, 상기 이동물체가 지표관측 방향으로 이동하고 있는 경우, 상기 지표관측방향으로의 속도를 이용하여 상기 MAI 간섭도에 대한 위상을 보정하기 위한 보정값을 결정하고, 상기 MAI 간섭도에서 상기 보정값의 차를 최종 MAI 간섭도로 결정하며, 상기 최종 MAI 간섭도로부터 상기 이동물체에 대한 비행방향의 속도를 관측하는 것이 바람직하다.And, when the moving object is moving in the surface observation direction, the second speed observation unit determines a correction value for correcting the phase for the MAI interference degree using the speed in the surface observation direction, It is preferable to determine the difference of the correction value in the MAI interference degree as the final MAI interference degree, and to observe the speed of the flight direction with respect to the moving object from the final MAI interference degree.

또한, 상기 보정값은, 비행방향으로의 레이더 안테나 길이, 상기 지표관측 방향으로의 속도, 정규화된 스퀸스 조정 변수(normalized squint adjustment parameter), 기선의 수평면에 투영된 길이 및 레이더 센서 탑재체의 속도를 이용하여 결정되는 것이 바람직하다. Further, the correction value may include the radar antenna length in the flight direction, the speed in the surface observation direction, the normalized squint adjustment parameter, the length projected on the horizontal plane of the base line, and the speed of the radar sensor payload. It is preferable to determine using.

본 발명에 의하면, ATI간섭도와 MAI간섭도를 이용하여 이동물체의 지표관측방향 속도와 비행방향 속도를 관측하여 2차원 속도를 결정하는 것이 가능하다. According to the present invention, it is possible to determine the two-dimensional speed by observing the surface observation speed and the flight direction speed of the moving object using the ATI interference and the MAI interference.

이는 교통량 관측, 자동차 속도 관측, 빙하 이동, 해수면 이동 등에 활용될 수 있다.This can be used for traffic observations, car speed observations, glacier movements and sea level movements.

본 출원에서 사용한 용어는 단지 특정한 실시예를 설명하기 위해 사용된 것으로, 본 발명을 한정하려는 의도가 아니다. 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한 복수의 표현을 포함한다. 본 출원에서, "포함하다" 또는 "가지다" 등의 용어는 명세서상에 기재된 특징, 숫자, 단계, 동작, 구성요소, 부품 또는 이들의 조합한 것이 존재함을 지정하려는 것이지, 하나 또는 그 이상의 다른 특징들이나 숫자, 단계, 동작, 구성요소, 부품 또는 이들을 조합한 것들의 존재 또는 부가 가능성을 미리 배제하지 않는 것으로 이해되어야 한다.The terminology used herein is for the purpose of describing particular example embodiments only and is not intended to be limiting of the present invention. A singular expression includes a plural expression unless the context clearly indicates otherwise. In this application, the terms "comprise" or "having" are intended to indicate that there is a feature, number, step, operation, component, part, or combination thereof described in the specification, and one or more other It is to be understood that the present invention does not exclude the possibility of the presence or the addition of features, numbers, steps, operations, components, parts, or a combination thereof.

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

도 1은 본 발명의 일 실시예에 따른 합성 개구 레이더 영상을 이용하여 물체 속도를 측정하는 장치의 블록도이고, 도 2는 본 발명의 일 실시예에 따른 합성 개구 레이더 영상을 이용한 물체 속도 측정 방법을 수행하는 과정을 나타내는 순서도이다.1 is a block diagram of an apparatus for measuring object velocity using a synthetic aperture radar image according to an embodiment of the present invention, and FIG. 2 is a method of measuring object velocity using a composite aperture radar image according to an embodiment of the present invention. This is a flow chart showing the process of performing.

도 1 및 도 2에 도시된 바와 같이, 제1 간섭도 제적부는 합성 개구 레이더 영상 자료를 이용하여 ATI(along-track SAR interferometry)간섭도를 제적하고(S210), 제2 간섭도 제작부(120)는 합성 개구 레이더 영상 자료를 이용하여 MAI(multiple aperture SAR interferometry)간섭도를 제작한다(S220). 그리고, 제1 속도 관측부(130)는 ATI 간섭도로부터 지표관측방향(ground range direction)의 이동물체 속도를 관측하고(S230), 제2 속도 관측부(140)는 MAI간섭도로부터 비행방향(along-track direction)의 이동물체 속도를 관측한다(S240). 마지막으로, 속도 결정부(150)는 지표관측방향의 속도 및 비행방향의 속도들을 결합하여 이동물체의 2차원 속도를 결정한다(S250). As shown in FIG. 1 and FIG. 2, the first interference canceling unit removes the ATI (along-track SAR interferometry) interference using the synthetic aperture radar image data (S210), and the second interference manufacturing unit 120. In operation S220, multiple aperture SAR interferometry (MAI) coherence is produced using the synthetic aperture radar image data. In addition, the first speed observer 130 observes the moving object speed in the ground range direction from the ATI interference degree (S230), and the second speed observer 140 measures the flight direction from the MAI interference degree (S230). The moving object velocity in the along-track direction is observed (S240). Finally, the speed determination unit 150 determines the two-dimensional speed of the moving object by combining the speed in the surface observation direction and the speed in the flight direction (S250).

이하에서는 도면을 참조하여 각 단계에 대하여 보다 상세하게 설명한다.Hereinafter, each step will be described in more detail with reference to the accompanying drawings.

먼저 이동물체에 대한 지표관측방향의 속도를 결정하는 과정을 설명한다. First, the process of determining the velocity in the surface observation direction with respect to the moving object will be described.

도 3은 2차원 이동물체의 속도를 관측하기 위한 센서탑재체의 기하를 나타내는 도면이다. 각 앞 및 뒤 센서는 일정한 거리()를 두고 센서탑채체에 탑재되며, 센서탑재체1과 센서탑재체2는 일정한 시간간격()을 두고 동일지역의 레이더 영상을 얻는다. 3 is a view showing the geometry of the sensor mounting body for observing the speed of the two-dimensional moving object. Each of the front and rear sensors is mounted on the sensor tower body at a certain distance (), and the sensor mounting body 1 and the sensor mounting body 2 obtain radar images of the same area with a certain time interval ().

이때 지표면의 이동물체에 대한 지표관측방향의 속도는 각 센서탑재체의 앞 및 뒤 센서로부터 각각 SLC 영상을 제작한 후, 이들로부터 만들어진 ATI 영상으로부터 계산된다. 반면에, 비행방향의 속도는 센서탑재체1의 앞 및 뒤 센서로부터 얻어진 SLC 영상 중 하나를 주영상으로 하고, 센서탑재체2의 앞 및 뒤 센서로부터 얻어진 SLC 영상 중 하나를 부영상으로 하여 얻어진 MAI 영상으로부터 계산된다. 상기한 앞 및 뒤 센서의 간격은 1.2m이며, 센서탑재체1과 센서탑재체2의 시간간격은 최대 약 0.08초인 것이 바람직하다.At this time, the velocity of the surface observation direction with respect to the moving object on the ground surface is calculated from the ATI images made from the SLC images after the SLC images are produced from the front and rear sensors of each sensor mounting body. On the other hand, the velocity of the flight direction is a MAI image obtained by using one of the SLC images obtained from the front and rear sensors of the sensor mounting body 1 as a main image, and one of the SLC images obtained from the front and rear sensors of the sensor mounting body 2 as a sub image. Is calculated from The distance between the front and rear sensors is 1.2m, and the time interval between the sensor mounting body 1 and the sensor mounting body 2 is preferably at most about 0.08 seconds.

구체적으로, 제1 간섭도 제작부(110)는 센서탑재체1 및 센서탑재체2 각각에 마련된 앞 및 뒤 센서로부터 얻어진 SLC 영상쌍으로부터 ATI 간섭도를 제작하는데, ATI 간섭도(

Figure 112009073910722-PAT00001
)는 앞의 레이더 센서에 의해 촬영된 전방관측 SLC 영상 자료(
Figure 112009073910722-PAT00002
)와 뒤의 합성 개구 레이더 영상 센서에 의해 촬영된 후방 관측 SLC 영상 자료(
Figure 112009073910722-PAT00003
)에 의하여 하기 수학식 1과 같이 결정된다.Specifically, the first interference manufacturing unit 110 produces the ATI interference degree from the SLC image pairs obtained from the front and rear sensors provided on the sensor mounting body 1 and the sensor mounting body 2, respectively.
Figure 112009073910722-PAT00001
) Is the front view SLC image data (
Figure 112009073910722-PAT00002
And rearview SLC image data captured by the composite aperture radar image sensor behind
Figure 112009073910722-PAT00003
Is determined by Equation 1 below.

Figure 112009073910722-PAT00004
Figure 112009073910722-PAT00004

여기서, *는 복소켤레(complex conjugate)이다. Where * is a complex conjugate.

또한 제1 속도 관측부(130)는 ATI 간섭도의 위상(phase)을 하기 수학식2와 같이 결정한다. In addition, the first speed observation unit 130 determines the phase of the ATI interference degree as shown in Equation 2 below.

Figure 112009073910722-PAT00005
Figure 112009073910722-PAT00005

여기서, 는 복소수의 각을 의미한다. Where is the angle of the complex number.

한편, 두 개의 레이더 안테나 중 하나는 송수신기 기능만, 다른 하나는 수신기 기능만을 지니는 레이더 시스템의 경우, 제1 속도 관측부(130)는 하기 수학식3과 같이 ATI 위상으로부터 시선방향으로의 속도(

Figure 112009073910722-PAT00006
)를 결정한다.On the other hand, in the case of a radar system having one of the two radar antennas only the transceiver function, the other only the receiver function, the first speed observer 130 is the speed from the ATI phase to the visual direction (Equation 3 below)
Figure 112009073910722-PAT00006
Is determined.

Figure 112009073910722-PAT00007
Figure 112009073910722-PAT00007

여기서,

Figure 112009073910722-PAT00008
는 레이더 파장(radar's wavelength)이고,
Figure 112009073910722-PAT00009
는 레이더 센서 탑재체(예를 들어, 위성, 비행기 등)의 속도이며,
Figure 112009073910722-PAT00010
는 앞 및 뒤 센서간의 거리이다. here,
Figure 112009073910722-PAT00008
Is the radar's wavelength,
Figure 112009073910722-PAT00009
Is the speed of the radar sensor payload (e.g., satellite, airplane, etc.)
Figure 112009073910722-PAT00010
Is the distance between the front and rear sensors.

또한, 제1 속도 관측부(130)는 수학식 3의 시선방향 속도를 하기 수학식 4와 같이 지표관측방향 속도(

Figure 112009073910722-PAT00011
)로 변환한다.In addition, the first speed observer 130 may measure the visual direction speed of Equation 3 as shown in Equation 4 below.
Figure 112009073910722-PAT00011
To.

Figure 112009073910722-PAT00012
Figure 112009073910722-PAT00012

여기서,

Figure 112009073910722-PAT00013
는 레이더파의 입사각을 의미한다.here,
Figure 112009073910722-PAT00013
Denotes the incident angle of the radar wave.

도 4는 본 발명의 일 실시에에 따른 ATI 간섭도를 이용하여 이동물체에 대한 지표관측방향의 속도를 결정한 예를 도시한 도면이다. 도 4와 같이 이동물체에 대한 지표관측방향의 속도를 결정하기 위해 지표면에 존재하는 네 이동물체를 촬영한 합성 개구 레이더 영상 자료를 모의(simulation)한 후 ATI 간섭도를 제작하고, 지표관측방향 속도를 계산한 것이다. 4 is a diagram illustrating an example of determining the speed of the surface observation direction with respect to a moving object using the ATI interference degree according to one embodiment of the present invention. In order to determine the speed of the surface observation direction with respect to the moving object as shown in FIG. 4, after simulating the synthetic aperture radar image data photographing the four moving objects on the surface, an ATI interference degree was produced, and the surface observation speed was measured. Is calculated.

도4의 (a)는 시뮬레이션에 이용된 네 이동물체에 대한 속도와 방향벡터를 나타낸 것으로서, 지표관측방향으로의 속도를 나타낸 도면이다. 도4의 (b)는 수학식 2를 이용하여 각 이동물체에 대한 ATI 간섭도의 위상값을 나타낸 도면이다. 그리고, 도 4의 (c)는 수학식 3을 이용하여 ATI 위상을 시선방향속도(

Figure 112009073910722-PAT00014
)로 변환한 값을 나타낸 도면이고, 도 4의 (d)는 수학식 4를 이용하여 시선방향속도를 지표관 측방향 속도(
Figure 112009073910722-PAT00015
)로 변환한 값을 나타낸 도면이다. FIG. 4 (a) shows the velocity and direction vectors for four moving objects used in the simulation, and shows the velocity in the surface observation direction. Figure 4 (b) is a diagram showing the phase value of the ATI interference degree for each moving object using the equation (2). 4 (c) shows the ATI phase by using the equation 3,
Figure 112009073910722-PAT00014
Figure 4 is a diagram showing the value converted to, Figure 4 (d) is a surface direction lateral velocity (
Figure 112009073910722-PAT00015
The figure which converted the value into ().

네 이동물체의 실제 지표관측방향 속도는 각각 -5, 14.14, 0, 10m/s였으며, ATI간섭도로부터 측정된 지표관측방향속도는 -4.94, 14.21, 0.05, 10.09m/s였다. 그러므로 모의된 네 이동물체의 오차는 0.06, 0.07, 0.05, 0.09m/s로 평균 약 0.068m/s(0.25km/h)의 오차로 ATI 간섭도로부터 지표면의 이동물체를 관측할 수 있다는 것을 확인할 수 있다. 그러나 이 값은 모의영상을 이용한 것으로 실제영상을 이용할 경우, 좀 더 큰 오차를 지닐 수 있음은 물론이다. 도면에서 도 4의 (a)에 있는 이동물체 위치와 도 4의 (b (c) 및 (d)에 있는 이동물체 위치가 다른 것은 이동물체의 도플러 효과에 기인한 것이다. The actual surface observation velocities of the four moving objects were -5, 14.14, 0, and 10 m / s, respectively, and the surface observation velocities measured from the ATI coherence were -4.94, 14.21, 0.05, and 10.09 m / s, respectively. Therefore, the simulated four moving objects are 0.06, 0.07, 0.05, 0.09 m / s, and the average of about 0.068 m / s (0.25 km / h) can be used to confirm that the moving objects on the surface can be observed from the ATI interference degree. Can be. However, this value uses a simulated image. Of course, when using an actual image, it may have a larger error. In FIG. 4, the position of the moving object in FIG. 4A and the position of the moving object in FIG. 4B (c) and (d) is due to the Doppler effect of the moving object.

다음은 이동 물체에 대한 비행방향의 속도를 결정하는 방법을 설명한다. The following describes how to determine the velocity of the flight direction for the moving object.

한편, 제2 간섭도 제작부(120)는 수직기선(perpendicular baseline)이 짧은 한 쌍의 합성 개구 레이더 영상 자료로부터 한 쌍의 전방관측(forward-looking) SLC영상과 한 쌍의 후방관측(backward-looking) SLC 영상을 제작하고, 상기 전방관측 SLC 영상 쌍으로부터 전방관측 InSAR 간섭도와 상기 후방관측 SLC 영상 쌍으로부터 후방관측 InSAR 간섭도를 제작하며, 상기한 두 InSAR(interferometric SAR)간섭도로부터 MAI 간섭도를 제작한다.On the other hand, the second interference producing unit 120 is a pair of forward-looking SLC image and a pair of backward-looking from a pair of composite aperture radar image data having a short vertical baseline (perpendicular baseline) A SLC image is produced, and a front-side InSAR interference from the front-side SLC image pair and a rear-side InSAR interference from the rear-view SLC image pair are produced, and MAI interference is obtained from the two inter-SAR (interferometric SAR) interference. To make.

수직기선(baseline)을 지니는 합성 개구 레이더 영상 쌍을 이용할 경우, MAI 간섭도는 시선방향으로 움직이는 물체에 영향을 받게 되므로, MAI 간섭도로부터 비행방향의 이동물체 속도를 결정할 수 있다. 한편, 합성 개구 레이더 영상 쌍은 실 제 서로 다른 위치에 놓이기 때문에 비행방향의 이동물체 속도를 지표관측방향 속도로부터 보정할 필요가 있다.When using a composite aperture radar image pair having a baseline, the MAI interference is influenced by an object moving in the visual direction, and thus the moving object velocity in the flight direction can be determined from the MAI interference. On the other hand, since the composite aperture radar image pairs are actually located at different positions, it is necessary to correct the moving object velocity in the flight direction from the surface observation direction velocity.

이하 도 5를 참조하여 보다 상세하게 설명한다. 도 5는 비행방향으로의 이동물체에 대한 MAI 위상 기하를 나타내는 도면이다.Hereinafter, a detailed description will be given with reference to FIG. 5. 5 is a diagram showing the MAI phase geometry for the moving object in the flight direction.

제2 간섭도 제작부(120)는 도 5에 도시된 센서탑재체1의 앞 및 뒤 센서로부터 얻어진 SLC 영상 중 하나를 주영상(master image)으로, 센서탑재체2의 앞 및 뒤 센서로부터 얻어진 SLC 영상 중 하나를 부영상(slave image)으로 하여 한 쌍의 전방관측 SLC영상과 한 쌍의 후방관측 SLC 영상을 제작한다. 여기서 동일한 지역에 대해 먼저 촬영된 영상을 주영상으로 하고, 이후에 촬영된 영상을 부영상으로 한다. 그리고, 상기 전방관측 SLC 영상 쌍으로부터 전방관측 InSAR 간섭도와 상기 후방관측 SLC 영상 쌍으로부터 후방관측 InSAR 간섭도를 제작하고, 이 두 InSAR간섭도로부터 MAI간섭도를 제작한다. 여기서 전 및 후방관측 SLC영상은 전 및 후방관측 도플러 중심주파수(

Figure 112009073910722-PAT00016
), 평균 도플러 중심주파수(
Figure 112009073910722-PAT00017
), 서브어퍼쳐 도플러 대역폭(Subaperture Doppler bandwidth;
Figure 112009073910722-PAT00018
)을 이용하여 제작되며, 일반적인 영상압축과정을 이용하여 제작되지만, 하기 두 가지 점에서 일반적인 영상압축과정과 다르다.The second interference degree producing unit 120 is a master image of one of the SLC images obtained from the front and rear sensors of the sensor mounting body 1 illustrated in FIG. 5, and among the SLC images obtained from the front and rear sensors of the sensor mounting body 2. One pair of front-view SLC images and one pair of rear-view SLC images are produced as slave images. In this case, the first image captured for the same region is used as the main image, and the image captured afterwards is the sub image. Then, the front observation InSAR interference from the front observation SLC image pair and the rear observation InSAR interference from the rear observation SLC image pair are produced, and the MAI interference is produced from the two InSAR interferences. Here, the front and rear observation SLC images are the front and rear observation Doppler center frequencies (
Figure 112009073910722-PAT00016
), Average Doppler center frequency (
Figure 112009073910722-PAT00017
), Subaperture Doppler bandwidth;
Figure 112009073910722-PAT00018
It is produced by using) and is produced by using a general image compression process, but is different from the general image compression process in the following two points.

첫째, 일반적으로 영상압축처리는 레이더 센서가 지니는 도플러 중심주파수와 도플러 대역폭을 이용하여 수행된다. 그러나 전방관측 SLC자료를 위한 영상압축처리는 전방관측 도플러 중심주파수(

Figure 112009073910722-PAT00019
)와 서브어퍼쳐 도플러 대역폭(
Figure 112009073910722-PAT00020
)을 이용하여 처리되며, 후방관측 SLC자료를 위한 영상압축처리는 후방관측 도플러 중심주파수(
Figure 112009073910722-PAT00021
)와 서브어퍼쳐 도플러 대역폭(
Figure 112009073910722-PAT00022
)를 이용하여 처리된다.First, image compression is generally performed using the Doppler center frequency and the Doppler bandwidth of the radar sensor. However, the image compression processing for forward observation SLC data does
Figure 112009073910722-PAT00019
) And subaperture Doppler bandwidth (
Figure 112009073910722-PAT00020
Image compression processing for the posterior observation SLC data.
Figure 112009073910722-PAT00021
) And subaperture Doppler bandwidth (
Figure 112009073910722-PAT00022
Is processed using).

둘째, 합성 개구 레이더 영상은 한 물체를 수천 번 촬영하여 이를 압축하는 과정을 통하여 SLC영상을 얻는다. 그러나 합성 개구 레이더 영상 자료에서 한 물체에 대한 정보는 센서와 물체사이의 거리의 변화에 의해서 직선의 형태로 나타나지 않고 곡선의 형태를 지니게 되는데 이를 레인지 셀 마이그레이션(range cell migration)이라 한다. 또한, 이러한 현상을 직선으로 보정해주는 과정을 레인지 셀 마이그레이션 보정(RCMC, range cell migration correction)이라 한다. 일반적으로 RCMC는 도플러 중심주파수를 이용하여 직선화한다. 그러나 전 및 후방관측 SLC자료제작시 서로 다른 도플러 중심주파수를 이용하기 때문에 영상의 가로방향(range방향)으로 위치차이를 지니게 되므로, 전 및 후방관측 SLC자료에서 동일 물체를 동일 위치로 만들기 위하여 평균 도플러 중심주파수(

Figure 112009073910722-PAT00023
)를 이용하여 RCMC를 적용한다.Second, the synthetic aperture radar image is obtained through the process of compressing the object thousands of times to obtain the SLC image. However, the information about an object in the synthetic aperture radar image data does not appear in the form of a straight line due to the change of the distance between the sensor and the object. It is called range cell migration. In addition, the process of straightening this phenomenon is called range cell migration correction (RCMC). In general, RCMC linearizes using the Doppler center frequency. However, since the different Doppler center frequencies are used in the production of front and rear observation SLC data, they have a positional difference in the horizontal direction of the image. Therefore, the average Doppler is used to make the same object in the front and rear observation SLC data. Center frequency
Figure 112009073910722-PAT00023
Apply RCMC using).

제2 간섭도 제작부(120)는 전방관측 SLC 주영상(master image,

Figure 112009073910722-PAT00024
)과 후방관측 SLC 부영상(slave image,
Figure 112009073910722-PAT00025
)을 이용하여 전방관측 InSAR간섭도(
Figure 112009073910722-PAT00026
)를 제작한다. 특히, 전방관측 InSAR간섭도는 하기 수학식 5와 같이 제작되는 것이 바람직하다. The second interference producing unit 120 is a front view SLC main image (master image,
Figure 112009073910722-PAT00024
) And posterior observation SLC sub-images
Figure 112009073910722-PAT00025
InSAR coherence using
Figure 112009073910722-PAT00026
). In particular, the front observation InSAR interference is preferably produced as shown in Equation 5.

Figure 112009073910722-PAT00027
Figure 112009073910722-PAT00027

그리고, 제2 간섭도 제작부(120)는 후방관측 SLC 주영상(

Figure 112009073910722-PAT00028
)과 후방관측 SLC 부영상(
Figure 112009073910722-PAT00029
)를 이용하여 하기 수학식 6와 같은 후방관측 InSAR간섭도(
Figure 112009073910722-PAT00030
)를 제작한다. And, the second interference degree production unit 120 is the rear view SLC main image (
Figure 112009073910722-PAT00028
) And posterior observation SLC sub-image (
Figure 112009073910722-PAT00029
BackSurveillance InSAR interference (
Figure 112009073910722-PAT00030
).

Figure 112009073910722-PAT00031
Figure 112009073910722-PAT00031

또한, 제2 간섭도 제작부(120)는 전방관측 InSAR간섭도와 후방관측 InSAR 간섭도를 이용하여 하기 수학식 7과 같은 MAI 간섭도(

Figure 112009073910722-PAT00032
)를 제작한다.In addition, the second interference degree manufacturing unit 120 uses the MAI interference (Equation 7) using the front observation InSAR interference and the rear observation InSAR interference degree (
Figure 112009073910722-PAT00032
).

Figure 112009073910722-PAT00033
Figure 112009073910722-PAT00033

제2 속도 관측부(140)는 MAI 간섭도의 위상(phase)을 MAI 간섭도로부터 하기 수학식8과 같이 결정한다.The second speed observing unit 140 determines the phase of the MAI interference degree from the MAI interference degree as shown in Equation 8 below.

Figure 112009073910722-PAT00034
Figure 112009073910722-PAT00034

한편, 제2 속도 관측부(140)는 MAI간섭도의 위상을 기초로 비행방향속도를 관측하는데, 이하 보다 구체적으로 설명한다.On the other hand, the second speed observation unit 140 observes the flight direction speed based on the phase of the MAI interference, it will be described in more detail below.

도 6은 본 발명의 일 실시예에 따른 비행방향으로의 이동물체에 대한 MAI위상기하를 나타내는 도면이다. 6 is a view showing the MAI phase up and down for the moving object in the flight direction according to an embodiment of the present invention.

도 6에서

Figure 112009073910722-PAT00035
은 비행방향으로의 레이더 안테나 길이(Effective antenna length in azimuth direction)이고,
Figure 112009073910722-PAT00036
는 레이다 파장(radar wavelength)이며,
Figure 112009073910722-PAT00037
은 정규화된 스퀸트 조정 변수(normalized squint adjustment parameter)이다. 그리고,
Figure 112009073910722-PAT00038
는 각각 레이더 센서 탑재체와 비행방향으로의 이동물체속도를 나타낸다. 여기서, 정규화된 스퀸트 조정변수는 안테나 어퍼쳐(aperture)의 크기를 전기적으로 조정하는 변수이다. In Figure 6
Figure 112009073910722-PAT00035
Is the effective antenna length in azimuth direction,
Figure 112009073910722-PAT00036
Is the radar wavelength,
Figure 112009073910722-PAT00037
Is a normalized squint adjustment parameter. And,
Figure 112009073910722-PAT00038
Are the radar sensor payloads and the moving object velocity in the flight direction, respectively. Here, the normalized squint adjustment variable is a variable for electrically adjusting the size of the antenna aperture.

일반적으로 합성 개구 레이더 영상 자료는 지표면의 한 물체를 수천 번 촬영하여 압축하여 고해상의 영상을 제작한다. 이때 전방관측 합성 개구 레이더 영상 자료는 안테나 어퍼쳐(aperture)의 크기 중 전방으로 관측된 어퍼쳐의 크기에 해당하는 자료를 이용하여 영상을 제작한 것이고, 후방관측 합성 개구 레이더 영상 자료는 후방으로 관측된 자료를 이용하여 제작한 것이다. 일반적으로 정규화된 스퀸드 조정변수는 0.5 (정확히 안테나 어퍼쳐를 반으로 나눔)를 이용한다. In general, synthetic aperture radar image data is compressed by photographing an object on the ground thousands of times to produce a high resolution image. In this case, the front observation composite aperture radar image data is produced by using the data corresponding to the size of the aperture observed forward among the size of the antenna aperture, and the rear observation composite aperture radar image data is observed backward. It was produced using the data. Normally, the normalized squind adjustment variable uses 0.5 (exactly dividing the antenna aperture in half).

또한,

Figure 112009073910722-PAT00039
는 각각 전방관측과 후방관측 주영상에 대한 경사거리(slant range distance)를 의미하며,
Figure 112009073910722-PAT00040
는 각각 전방관측과 후방관측 부영상에 대한 경사거리를 의미하고,
Figure 112009073910722-PAT00041
는 센서탑재체와 지표면간의 최단경사거리를 의미한다.
Figure 112009073910722-PAT00042
는 각각 전방관측과 후방관측영상을 촬영할 때의 센서탐재체의 위치를 나타낸다. 또한,
Figure 112009073910722-PAT00043
은 주영상의 전방관측시 이동물체의 위치,
Figure 112009073910722-PAT00044
는 주영상의 후방관측시 이동물체의 위치,
Figure 112009073910722-PAT00045
는 부영상의 전방관측시 이동물체의 위치 및
Figure 112009073910722-PAT00046
는 부영상의 후방관측시 이동물체의 위치를 보인다.Also,
Figure 112009073910722-PAT00039
Denotes the slant range distance for the front and rear main images, respectively.
Figure 112009073910722-PAT00040
Are the inclination distances for the anterior and posterior subimages, respectively.
Figure 112009073910722-PAT00041
Means the shortest inclination distance between the sensor mounting body and the ground surface.
Figure 112009073910722-PAT00042
Indicates the positions of the sensor probes when the front and rear observation images are taken, respectively. Also,
Figure 112009073910722-PAT00043
Is the position of the moving object in front view of the main image,
Figure 112009073910722-PAT00044
Is the position of the moving object in the rear view of the main image,
Figure 112009073910722-PAT00045
Is the position of the moving object in front view of the sub-image and
Figure 112009073910722-PAT00046
Shows the position of the moving object in the rear view of the sub-image.

도 6에서 전방관측 InSAR 간섭도의 경사거리(

Figure 112009073910722-PAT00047
)는 전방관측 주영상의 경사거리(
Figure 112009073910722-PAT00048
)와 전방관측 부영상의 경사거리(
Figure 112009073910722-PAT00049
)의 차에 의해서 하기 수학식 9와 같이 결정된다. In Fig. 6 the inclination distance of the front observation InSAR interference degree (
Figure 112009073910722-PAT00047
) Is the inclination distance (
Figure 112009073910722-PAT00048
) And the inclination distance of the anterior subview
Figure 112009073910722-PAT00049
Is determined as shown in Equation 9 below.

Figure 112009073910722-PAT00050
Figure 112009073910722-PAT00050

여기서 전방관측 주영상의 경사거리(

Figure 112009073910722-PAT00051
)와 전방관측 부영상의 경사거리(
Figure 112009073910722-PAT00052
)는 각각 하기 수학식 10 및 수학식 11과 같이 결정된다. Where the inclination distance of the main observation image (
Figure 112009073910722-PAT00051
) And the inclination distance of the anterior subview
Figure 112009073910722-PAT00052
) Are determined as in Equations 10 and 11, respectively.

Figure 112009073910722-PAT00053
Figure 112009073910722-PAT00053

Figure 112009073910722-PAT00054
Figure 112009073910722-PAT00054

한편, 후방관측 InSAR 간섭도의 경사거리(

Figure 112009073910722-PAT00055
)도 후방관측 주영상의 경사거리(
Figure 112009073910722-PAT00056
)와 후방관측 부영상의 경사거리(
Figure 112009073910722-PAT00057
)의 차에 의해서 하기 수학식 12와 같이 결정된다.On the other hand, the inclination distance of the rearview InSAR interference degree (
Figure 112009073910722-PAT00055
Degrees of inclination of the posterior main image
Figure 112009073910722-PAT00056
) And the oblique distance of the posterior subimage
Figure 112009073910722-PAT00057
Is determined as shown in Equation 12 below.

Figure 112009073910722-PAT00058
Figure 112009073910722-PAT00058

여기서 후방관측 주영상의 경사거리(

Figure 112009073910722-PAT00059
)와 후방관측 부영상의 경사거리(
Figure 112009073910722-PAT00060
)는 각각 하기 수학식 13 및 수학식 14와 같이 결정된다.Where the oblique distance of the posterior main image (
Figure 112009073910722-PAT00059
) And the oblique distance of the posterior subimage
Figure 112009073910722-PAT00060
) Are determined as in Equations 13 and 14, respectively.

Figure 112009073910722-PAT00061
Figure 112009073910722-PAT00061

Figure 112009073910722-PAT00062
Figure 112009073910722-PAT00062

결국 MAI간섭도의 경사거리(

Figure 112009073910722-PAT00063
)는 수학식 9와 수학식 12를 이용하여 하기 수학식 15와 같이 결정된다.The slope distance of MAI interference
Figure 112009073910722-PAT00063
) Is determined as in Equation 15 using Equations 9 and 12.

Figure 112009073910722-PAT00064
Figure 112009073910722-PAT00064

여기서,

Figure 112009073910722-PAT00065
인 것과 탑재체에 비해 이동물체의 속도가 일반적으로 매우 작기 때문에
Figure 112009073910722-PAT00066
이므로 하기 수학식 16과 같이 단순화 할 수 있다.here,
Figure 112009073910722-PAT00065
Because the speed of the moving object is generally very small compared to
Figure 112009073910722-PAT00066
Therefore, it can be simplified as shown in Equation 16 below.

Figure 112009073910722-PAT00067
Figure 112009073910722-PAT00067

제2 속도 관측부(140)는 MAI 간섭도의 위상을 수학식 16을 이용하여 하기 수학식 17과 같이 결정한다.The second speed observation unit 140 determines the phase of the MAI interference degree using Equation 16 as shown in Equation 17 below.

Figure 112009073910722-PAT00068
Figure 112009073910722-PAT00068

결국 제2 속도 관측부(140)는 비행방향으로 이동하는 물체의 속도(

Figure 112009073910722-PAT00069
)를수학식 17에 정의된 MAI간섭도 위상(
Figure 112009073910722-PAT00070
)를 이용하여 하기 수학식 18와 같이 결 정한다. As a result, the second speed observer 140 measures the speed of the object moving in the flight direction (
Figure 112009073910722-PAT00069
MAI coherence phase (
Figure 112009073910722-PAT00070
) Is determined as in Equation 18 below.

Figure 112009073910722-PAT00071
Figure 112009073910722-PAT00071

그러므로, 비행방향 속도는 상기 제작된 MAI 간섭도에서 수학식 8을 이용하여 MAI간섭도의 위상을 계산하고, 이 위상을 수학식 18을 이용하여 변환하는 과정을 통하여 얻어진다. Therefore, the flight direction velocity is obtained by calculating the phase of MAI interference using Equation 8 and converting the phase by using Equation 18 from the prepared MAI interference.

한편, MAI간섭도의 위상은 비행방향으로 움직이는 이동물체의 속도 뿐 아니라 지표관측방향(즉, 비행방향과 교차방향)으로 움직이는 이동물체의 속도에 의해서도 변하기 때문에 이를 보정해주어야 한다.On the other hand, the phase of the MAI interference should be corrected because it changes not only by the speed of the moving object moving in the flight direction but also by the speed of the moving object moving in the surface observation direction (that is, the cross direction with the flying direction).

비행방향 속도보정은 제2 간섭도 제작부(120)에서 제작된 MAI간섭도의 위상으로부터 지표관측방향으로 움직이는 이동물체의 속도에 의한 위상을 보정하는 단계를 포함한다.Flight direction speed correction includes a step of correcting the phase by the speed of the moving object moving in the surface observation direction from the phase of the MAI interference produced by the second interference producing unit 120.

도 7은 본 발명의 일 실시예에 따른 지표관측방향으로의 이동물체에 대한 MAI위상기하를 나타내는 도면이다.7 is a view showing the MAI phase up and down for the moving object in the surface observation direction according to an embodiment of the present invention.

도 7에서 각 기호에 대한 설명은 도 5와 동일하며, 추가적으로

Figure 112009073910722-PAT00072
는 지표관측방향으로의 이동물체속도를 나타내고,
Figure 112009073910722-PAT00073
는 각각 주영상과 부영상에서 센서탑재체와 지표면간의 최단경사거리를 의미한다.
Figure 112009073910722-PAT00074
Figure 112009073910722-PAT00075
는 각각 주영 상과 부영상에서
Figure 112009073910722-PAT00076
의 지표면 상의 거리를 나타내며,
Figure 112009073910722-PAT00077
는 각각 주영상과 부영상의 위성고도이다.
Figure 112009073910722-PAT00078
는 주영상에 대한 전방관측과 후방관측영상을 촬영할 때의 센서탐재체의 위치이며,
Figure 112009073910722-PAT00079
는 부영상에 대한 전방관측과 후방관측영상을 촬영할 때의 센서탑재체의 위치를 나타낸다.Description of each symbol in FIG. 7 is the same as FIG. 5, and additionally
Figure 112009073910722-PAT00072
Represents the moving object velocity in the surface observation direction,
Figure 112009073910722-PAT00073
Denotes the shortest inclination distance between the sensor mounting body and the ground surface in the main and sub-images, respectively.
Figure 112009073910722-PAT00074
and
Figure 112009073910722-PAT00075
In the main and sub-images
Figure 112009073910722-PAT00076
Represents the distance on Earth's surface,
Figure 112009073910722-PAT00077
Are the altitudes of the main and sub images respectively.
Figure 112009073910722-PAT00078
Is the position of the sensor probe when the front and rear observation images of the main image are taken.
Figure 112009073910722-PAT00079
and Indicates the position of the sensor mounting body when photographing the front and rear observation images for the sub-image.

도 7에서 전방관측 InSAR 간섭도의 경사거리(

Figure 112009073910722-PAT00081
)는 전방관측 주영상의 경사거리(
Figure 112009073910722-PAT00082
)와 전방관측 부영상의 경사거리(
Figure 112009073910722-PAT00083
)의 차에 의해서 결정된다.The inclination distance of the forward observation InSAR interference degree in FIG.
Figure 112009073910722-PAT00081
) Is the inclination distance (
Figure 112009073910722-PAT00082
) And the inclination distance of the anterior subview
Figure 112009073910722-PAT00083
) Is determined by the difference.

Figure 112009073910722-PAT00084
Figure 112009073910722-PAT00084

여기서 전방관측 주영상의 경사거리(

Figure 112009073910722-PAT00085
)와 전방관측 부영상의 경사거리(
Figure 112009073910722-PAT00086
)는 각각 하기 수학식으로 정의될 수 있다.Where the inclination distance of the main observation image (
Figure 112009073910722-PAT00085
) And the inclination distance of the anterior subview
Figure 112009073910722-PAT00086
) May be defined by the following equations, respectively.

Figure 112009073910722-PAT00087
Figure 112009073910722-PAT00087

Figure 112009073910722-PAT00088
Figure 112009073910722-PAT00088

또한, 후방관측 InSAR 간섭도의 경사거리(

Figure 112009073910722-PAT00089
)는 후방관측 주영상의 경사거리(
Figure 112009073910722-PAT00090
)와 후방관측 부영상의 경사거리(
Figure 112009073910722-PAT00091
)의 차에 의해서 결정된다.Also, the inclination distance of the rearview InSAR interference degree (
Figure 112009073910722-PAT00089
) Is the inclination distance (
Figure 112009073910722-PAT00090
) And the oblique distance of the posterior subimage
Figure 112009073910722-PAT00091
) Is determined by the difference.

Figure 112009073910722-PAT00092
Figure 112009073910722-PAT00092

여기서

Figure 112009073910722-PAT00093
는 각각 하기 수학식으로 정의될 수 있다.here
Figure 112009073910722-PAT00093
Each can be defined by the following equation.

Figure 112009073910722-PAT00094
Figure 112009073910722-PAT00094

Figure 112009073910722-PAT00095
Figure 112009073910722-PAT00095

결국 MAI간섭도의 경사거리(

Figure 112009073910722-PAT00096
)는 수학식 19와 수학식 22를 이용하여 하기 수학식으로 정의할 수 있다.The slope distance of MAI interference
Figure 112009073910722-PAT00096
) Can be defined by the following equation (19) and (22).

Figure 112009073910722-PAT00097
Figure 112009073910722-PAT00097

여기서,

Figure 112009073910722-PAT00098
인 것을 이용하면 하기 수학식과 같이 단순화할 수 있다.here,
Figure 112009073910722-PAT00098
Using can be simplified as in the following equation.

Figure 112009073910722-PAT00099
Figure 112009073910722-PAT00099

또한, 도 7로부터 기선(baseline)

Figure 112009073910722-PAT00100
와 수평면과 기선이 이루는 각
Figure 112009073910722-PAT00101
를 이용하면
Figure 112009073910722-PAT00102
으로 대체할 수 있고, 상기 값에 비하여
Figure 112009073910722-PAT00103
는 매우 작은 값이므로 무시할 수 있다. 이와 같은 단순화를 통하여 수학식 26은 하기와 같이 표현된다.Also, baseline from FIG. 7.
Figure 112009073910722-PAT00100
Angle between the horizontal plane and the base line
Figure 112009073910722-PAT00101
If you use
Figure 112009073910722-PAT00102
Can be replaced with
Figure 112009073910722-PAT00103
Is a very small value and can be ignored. Through this simplification, Equation 26 is expressed as follows.

Figure 112009073910722-PAT00104
Figure 112009073910722-PAT00104

MAI 간섭도의 위상은 수학식 27을 이용하여 하기 수학식과 같이 정의된다.The phase of the MAI interference is defined by the following equation (27).

Figure 112009073910722-PAT00105
Figure 112009073910722-PAT00105

결국 지표관측방향으로 이동하는 물체의 속도(

Figure 112009073910722-PAT00106
)에 의하여 MAI 간섭도의 위상은 수학식 28과 같이 변환된다. Eventually, the velocity of the object moving in the surface observation direction (
Figure 112009073910722-PAT00106
), The phase of the MAI interference is converted as shown in Equation 28.

즉, MAI간섭도의 위상은 비행방향의 속도를 계산하기 위하여 이용되는 것으로 지표관측방향의 이동속도에 따른 MAI 간섭도 위상은 보정되어야 한다.That is, the phase of MAI interference is used to calculate the velocity in the flight direction. The phase of MAI interference according to the movement speed in the surface observation direction should be corrected.

이와 같은 보정과정은 지표관측방향 속도관측 단계(S1)에 의하여 추출된 지표관측방향으로 이동하는 물체의 속도(

Figure 112009073910722-PAT00107
)로부터 변형된 MAI 위상값을 수학식 28로부터 계산한다. 그리고, 수학식 8에 의하여 제작된 MAI 간섭도에서 변형된 MAI 위상값을 뺀 값을 최종 MAI 간섭도의 위상으로 결정하고, 이를 수학식 18을 적용하여 비행방향의 물체 속도(
Figure 112009073910722-PAT00108
)를 계산한다.This correction process is performed by the speed of the object moving in the surface observation direction extracted by the surface observation direction observation step (S1).
Figure 112009073910722-PAT00107
The MAI phase value modified from) is calculated from Equation 28. Then, the value obtained by subtracting the modified MAI phase value from the MAI interference degree produced by Equation 8 is determined as the phase of the final MAI interference degree, and the equation 18 is applied to the object velocity in the flight direction (
Figure 112009073910722-PAT00108
Calculate

도 8은 본 발명의 일 실시예에 따른 지표면에 존재하는 네 이동물체에 대한 앞 및 뒤의 합성 개구 레이더 영상 자료를 모의(simulation)한 후 MAI 간섭도를 제작하고 비행방향 속도를 계산한 결과를 도시한 도면이다.FIG. 8 illustrates simulation results of the front and rear composite aperture radar image data of four moving objects existing on the ground according to an embodiment of the present invention, and then a MAI interference calculation and flight direction velocity calculations. Figure is shown.

도9의 (a)는 시뮬레이션에 이용된 네 이동물체에 대한 속도와 방향벡터를 보이며, 비행방향으로의 속도를 나타내고, 도9의 (b)는 각 이동물체의 MAI 간섭도의 위상값을 나타낸다. 도9의 (c)는 수학식 18을 이용하여 MAI 위상을 비행방향속도(

Figure 112009073910722-PAT00109
)로 변환한 값을 나타낸다. 도9의 (d)는 수학식 28를 이용하여 지표관측방향 속도(
Figure 112009073910722-PAT00110
)에 의하여 변형된 MAI 위상을 보정한 비행방향속도(
Figure 112009073910722-PAT00111
)를 나타낸다. FIG. 9 (a) shows the velocity and direction vector for the four moving objects used in the simulation, shows the speed in the flight direction, and FIG. 9 (b) shows the phase value of the MAI interference degree of each moving object. . 9 (c) shows the MAI phase in flight direction velocity (Equation 18).
Figure 112009073910722-PAT00109
) Is converted to). 9 (d) shows the surface observation direction velocity (Equation 28).
Figure 112009073910722-PAT00110
Flight direction velocity corrected MAI phase modified by
Figure 112009073910722-PAT00111
).

네 이동물체의 실제 비행방향 속도는 각각 0, 14.14, 10, 0m/s였으며, MAI간섭도로부터 측정된 최종 비행방향속도는 0.01, 14.51, 9.90, -0.18m/s였다. 그러므로 모의된 네 이동물체의 오차는 0.01, 0.37, 0.1, -0.18m/s로 평균 약 0.11m/s(0.39km/h)의 오차로 지표면의 이동물체의 비행방향 속도를 관측할 수 있다는 것을 알 수 있다. 그러나 이 값은 모의영상을 이용한 것으로 실제영상을 이용할 경우, 좀 더 큰 오차를 지닐 수 있음은 물론이다. 도 8의 (a)의 이동물체 위치와 도 8의 (b), (c) 및 (d)에서의 이동물체 위치가 다른 것은 이동물체의 도플러 효과에 의하여 나타난다.The actual flight speeds of the four moving objects were 0, 14.14, 10 and 0m / s, respectively, and the final flight speeds measured from MAI interference were 0.01, 14.51, 9.90 and -0.18m / s. Therefore, the simulated error of the four moving objects is 0.01, 0.37, 0.1, -0.18m / s, and the average velocity of the moving object on the ground can be observed with an error of about 0.11m / s (0.39km / h). Able to know. However, this value uses a simulated image. Of course, when using an actual image, it may have a larger error. The position of the moving object in FIG. 8 (a) and the position of the moving object in FIGS. 8 (b), (c) and (d) is indicated by the Doppler effect of the moving object.

또한, 이상에서는 본 발명의 바람직한 실시예에 대하여 도시하고 설명하였지만, 본 발명은 상술한 특정의 실시예에 한정되지 아니하며, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 당해 발명이 속하는 기술분야에서 통상의 지식을 가진자에 의해 다양한 변형실시가 가능한 것은 물론이고, 이러한 변형실시들은 본 발명의 기술적 사상이나 전망으로부터 개별적으로 이해되어져서는 안될 것이다.In addition, although the preferred embodiment of the present invention has been shown and described above, the present invention is not limited to the specific embodiments described above, but the technical field to which the invention belongs without departing from the spirit of the invention claimed in the claims. Of course, various modifications can be made by those skilled in the art, and these modifications should not be individually understood from the technical spirit or the prospect of the present invention.

도 1은 본 발명의 일 실시예에 따른 합성 개구 레이더 영상을 이용하여 물체 속도를 측정하는 장치의 블록도,1 is a block diagram of an apparatus for measuring an object velocity using a composite aperture radar image according to an embodiment of the present invention;

도 2는 본 발명의 일 실시예에 따른 합성 개구 레이더 영상을 이용한 물체 속도 측정 방법을 수행하는 과정을 나타내는 순서도,2 is a flowchart illustrating a process of performing an object velocity measuring method using a composite aperture radar image according to an embodiment of the present invention;

도 3은 2차원 이동물체의 속도를 관측하기 위한 센서탑재체의 기하를 나타내는 도면,3 is a view showing the geometry of the sensor mounting body for observing the speed of the two-dimensional moving object,

도 4는 본 발명의 일 실시에에 따른 ATI 간섭도를 이용하여 이동물체에 대한 지표관측방향의 속도를 결정한 예를 도시한 도면,4 is a view showing an example of determining the speed of the surface observation direction for a moving object using the ATI interference degree according to an embodiment of the present invention,

도 5는 비행방향으로의 이동물체에 대한 MAI 위상 기하를 나타내는 도면,5 is a diagram showing the MAI phase geometry for the moving object in the flight direction,

도 6은 본 발명의 일 실시예에 따른 비행방향으로의 이동물체에 대한 MAI위상기하를 나타내는 도면,6 is a view showing the MAI phase upper and lower for the moving object in the flight direction according to an embodiment of the present invention,

도 7은 본 발명의 일 실시예에 따른 지표관측방향으로의 이동물체에 대한 MAI위상기하를 나타내는 도면, 그리고,7 is a view showing the MAI phase up and down for the moving object in the surface observation direction according to an embodiment of the present invention, and

도 8은 본 발명의 일 실시예에 따른 지표면에 존재하는 네 이동물체에 대한 앞 및 뒤의 합성 개구 레이더 영상 자료를 모의(simulation)한 후 MAI 간섭도를 제작하고 비행방향 속도를 계산한 결과를 도시한 도면이다.FIG. 8 illustrates simulation results of the front and rear composite aperture radar image data of four moving objects existing on the ground according to an embodiment of the present invention, and then a MAI interference calculation and flight direction velocity calculations. Figure is shown.

* 도면의 주요 부분에 대한 부호의 설명 *Explanation of symbols on the main parts of the drawings

110: 제1 간섭도 제작부 120: 제2 간섭도 제작부 110: first interference degree production unit 120: second interference degree production unit

130: 제1 속도 관측부 140: 제2 속도 관측부 130: first speed observing unit 140: second speed observing unit

150: 속도 결정부 150: speed determination unit

Claims (18)

합성 개구 레이더 영상을 이용하여 ATI(along-track SAR interferometry)간섭도를 제작하는 제1 제작 단계;A first manufacturing step of producing an ATI (along-track SAR interferometry) interference using a synthetic aperture radar image; 상기 합성 개구 레이더 영상 자료를 이용하여 MAI(multiple aperture SAR interferometry)간섭도를 제작하는 제2 제작 단계;A second manufacturing step of manufacturing multiple aperture SAR interferometry (MAI) coherence using the synthesized aperture radar image data; 상기 ATI 간섭도로부터 이동물체에 대한 지표관측방향의 속도를 관측하는 제1 관측 단계;A first observation step of observing a velocity in a surface observation direction with respect to a moving object from the ATI interference degree; 상기 MAI 간섭도로부터 상기 이동물체에 대한 비행방향의 속도를 관측하는 제2 관측 단계; 및A second observation step of observing a velocity of a flying direction with respect to the moving object from the MAI interference degree; And 상기 지표관측방향의 속도 및 상기 비행방향의 속도를 결합하여 상기 이동물체의 2차원 속도를 결정하는 단계;를 포함하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법.And determining a two-dimensional velocity of the moving object by combining the velocity in the surface observation direction and the velocity in the flight direction. 제 1항에 있어서,The method of claim 1, 상기 제1 제작 단계는,The first production step, 레이더 센서를 앞 및 뒤에 구비한 한 쌍의 센서탑재체에서 한 쌍의 전방관측 SLC 영상과 한 쌍의 후방관측 SLC 영상을 이용하여 ATI 간섭도를 제작하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법. An ATI interference degree is produced by using a pair of front-view SLC images and a pair of rear-view SLC images in a pair of sensor mounting bodies equipped with the front and rear of the radar sensor. How to measure speed. 제 2항에 있어서,3. The method of claim 2, 상기 ATI 간섭도는 The ATI interference degree is 상기 전방관측 SLC 영상과 상기 후방관측 SLC 영상의 컬레복소수를 곱함으로써 결정되는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법. The object velocity measuring method using a composite aperture radar image, characterized in that determined by multiplying the Cull complex number of the front observation SLC image and the rear observation SLC image. 제 1항에 있어서, The method of claim 1, 상기 제1 관측 단계는,The first observation step, 상기 ATI 간섭도로부터 ATI 간섭도의 위상을 결정하는 단계; Determining a phase of an ATI interference degree from the ATI interference degree; 상기 ATI 간섭도의 위상, 레이더의 파장, 상기 센서탑재체의 파장 및 상기 앞의 레이더 센서와 상기 뒤의 레이더 센서간의 거리를 이용하여 상기 이동물체에 대한 시선방향의 속도를 결정하는 단계; 및 Determining a speed in the visual direction with respect to the moving object using the phase of the ATI interference degree, the wavelength of the radar, the wavelength of the sensor mounting body, and the distance between the front radar sensor and the rear radar sensor; And 상기 시선 방향의 속도 및 상기 레이더의 입사각을 이용하여 상기 이동물체에 대한 지표관측 방향의 속도를 관측하는 단계;를 포함하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법. And observing the speed in the surface observation direction with respect to the moving object by using the speed in the eye direction and the incident angle of the radar. 제 1항에 있어서,The method of claim 1, 상기 제2 제작 단계는,The second production step, 레이더 센서를 앞 및 뒤에 구비한 한 쌍의 센서탑재체에서 한 쌍의 전방관측(forward-looking) SLC영상과 한 쌍의 후방관측(backward-looking) SLC 영상을 제작하고, 상기 전방관측 SLC 영상 쌍으로부터 전방관측 InSAR 간섭도와 상기 후방관측 SLC 영상 쌍으로부터 후방관측 InSAR 간섭도를 제작하며, 상기 전방관측 InSAR 간섭도 및 상기 후방관측 InSAR 간섭도로부터 MAI 간섭도를 제작하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법.A pair of forward-looking SLC images and a pair of backward-looking SLC images are produced from a pair of sensor mounting bodies equipped with radar sensors at the front and the rear, and from the pair of front-view SLC images A composite aperture radar image is generated from the front observation InSAR interference and the rear observation SLC image pair. Method of measuring the speed of the object used. 제 5항에 있어서,The method of claim 5, 상기 한쌍의 전방 관측 SLC 영상 중 하나를 주영상 및 나머지를 부영상으로 결정하고, 상기 한 쌍의 후방 관측 SLC 영상 중 하나를 주영상 및 나머지를 부영상으로 설정하며,One of the pair of forward-viewing SLC images is determined as a main image and the rest as sub-pictures, and one of the pair of backward-viewing SLC images is set as a main-image and the rest as sub-pictures, 상기 부영상의 컬러복소수를 상기 부영상에 대응되는 주영상에 곱함으로써 상기 전방관측 InSAR 간섭도 및 상기 후방관측 InSAR 간섭도를 제작하며,Producing the front observation InSAR interference and the rear observation InSAR interference by multiplying the color complex number of the sub-image by the main image corresponding to the sub-image, 상기 후방관측 InSAR 간섭도의 컬레복소수를 상기 전방관측 InSAR 간섭도에 곱함으로써 상기 MAI 간섭도를 제작하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법. And producing the MAI interference by multiplying the front complex InSAR interference degree by the COLE complex number of the rear observation InSAR interference degree. 제 1항에 있어서,The method of claim 1, 상기 MAI 간섭도로부터 MAI 간섭도의 위상을 결정하는 단계; 및 Determining a phase of the MAI interference from the MAI interference; And 상기 MAI 간섭도에 대한 위상, 두 레이더 센서 탑재체가 동일 지역을 촬영하는 시간간격, 비행방향으로의 레이더 안테나 길이 및 정규화된 스퀸스 조정 변수(normalized squint adjustment parameter)를 이용하여 상기 비행방향의 속도를 결정하는 단계;를 포함한느 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법. The velocity of the flight direction is determined using the phase of the MAI interference, the time interval between the two radar sensor payloads photographing the same area, the length of the radar antenna in the direction of flight, and a normalized squint adjustment parameter. Determining the velocity of the object using the composite aperture radar image. 제 7항에 있어서,The method of claim 7, wherein 상기 MAI 간섭도의 위상을 결정하는 단계는, Determining the phase of the MAI interference degree, 상기 이동물체가 지표관측 방향으로 이동하고 있는 경우, 상기 지표관측방향으로의 속도를 이용하여 상기 MAI 간섭도에 대한 위상을 보정하기 위한 보정값을 결정하는 단계; 및 Determining a correction value for correcting a phase with respect to the MAI interference degree using the speed in the surface observation direction when the moving object is moving in the surface observation direction; And 상기 MAI 간섭도에서 상기 보정값의 차를 최종 MAI 간섭도로 결정하는 단계;를 포함하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법. And determining a difference between the correction values in the MAI interference degree as a final MAI interference degree. 제 7항에 있어서,The method of claim 7, wherein 상기 보정값은,The correction value is, 비행방향으로의 레이더 안테나 길이, 상기 지표관측 방향으로의 속도, 정규화된 스퀸스 조정 변수(normalized squint adjustment parameter), 기선의 수평면에 투영된 길이 및 레이더 센서 탑재체의 속도를 이용하여 결정되는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법. Radar antenna length in the direction of flight, speed in the surface observation direction, normalized squint adjustment parameter, length projected on the horizontal plane of the baseline, and speed of the radar sensor payload. An object velocity measurement method using a composite aperture radar image. 합성 개구 레이더 영상 자료를 이용하여 ATI(along-track SAR interferometry)간섭도를 제작하는 제1 간섭도 제작부;A first interference producing unit for producing an ATI (along-track SAR interferometry) interference using synthetic aperture radar image data; 상기 합성 개구 레이더 영상 자료를 이용하여 MAI(multiple aperture SAR interferometry)간섭도를 제작하는 제2 간섭도 제작부;A second interference producing unit configured to manufacture multiple aperture SAR interferometry (MAI) interference using the synthesized aperture radar image data; 상기 ATI 간섭도로부터 이동물체에 대한 지표관측방향의 속도를 관측하는 제1 속도 관측부;A first speed observing unit for observing the speed in the surface observation direction with respect to a moving object from the ATI interference degree; 상기 MAI 간섭도로부터 상기 이동물체에 대한 비행방향의 속도를 관측하는 제2 속도 관측부; 및 A second speed observing unit for observing the speed of the flying direction with respect to the moving object from the MAI interference degree; And 상기 지표관측방향의 속도 및 상기 비행방향의 속도를 결합하여 상기 이동물체의 2차원 속도를 결정하는 속도 결정부;를 포함하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 장치.And a speed determiner configured to determine the two-dimensional speed of the moving object by combining the speed in the surface observation direction and the speed in the flight direction. 제 10항에 있어서,The method of claim 10, 상기 제1 간섭도 제작부는, The first interference degree production unit, 레이더 센서를 앞 및 뒤에 구비한 한 쌍의 센서탑재체에서 한 쌍의 전방관측 SLC 영상과 한 쌍의 후방관측 SLC 영상을 이용하여 ATI 간섭도를 제작하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 장치. An ATI interference degree is produced by using a pair of front-view SLC images and a pair of rear-view SLC images in a pair of sensor mounting bodies equipped with the front and rear of the radar sensor. Speed measuring device. 제 11항에 있어서,The method of claim 11, 상기 ATI 간섭도는 The ATI interference degree is 상기 전방관측 SLC 영상과 상기 후방관측 SLC 영상의 컬레복소수를 곱함으로 써 결정되는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 장치. The apparatus for measuring the velocity of an object using a composite aperture radar image, which is determined by multiplying the COLE complex number of the front observation SLC image and the rear observation SLC image. 제 11항에 있어서, The method of claim 11, 상기 제1 속도 관측부는, The first speed observation unit, 상기 ATI 간섭도로부터 ATI 간섭도의 위상을 결정하고, 상기 ATI 간섭도의 위상, 레이더의 파장, 상기 센서탑재체의 파장 및 상기 앞의 레이더 센서와 상기 뒤의 레이더 센서간의 거리를 이용하여 상기 이동물체에 대한 시선방향의 속도를 결정하며, 상기 시선 방향의 속도 및 상기 레이더의 입사각을 이용하여 상기 이동물체에 대한 지표관측 방향의 속도를 관측하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 장치. The phase of the ATI interference degree is determined from the ATI interference degree, and the moving object is determined using the phase of the ATI interference degree, the wavelength of the radar, the wavelength of the sensor mounting body, and the distance between the front radar sensor and the rear radar sensor. The speed of the object is determined using a composite aperture radar image, the speed of which is determined in the visual direction, and the speed of the surface observation direction with respect to the moving object is observed using the speed in the eye direction and the incident angle of the radar. Device. 제 10항에 있어서,The method of claim 10, 상기 제2 간섭도 제작부는, The second interference degree production unit, 레이더 센서를 앞 및 뒤에 구비한 한 쌍의 센서탑재체에서 한 쌍의 전방관측(forward-looking) SLC영상과 한 쌍의 후방관측(backward-looking) SLC 영상을 제작하고, 상기 전방관측 SLC 영상 쌍으로부터 전방관측 InSAR 간섭도와 상기 후방관측 SLC 영상 쌍으로부터 후방관측 InSAR 간섭도를 제작하며, 상기 전방관측 InSAR 간섭도 및 상기 후방관측 InSAR 간섭도로부터 MAI 간섭도를 제작하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 장치.A pair of forward-looking SLC images and a pair of backward-looking SLC images are produced from a pair of sensor mounting bodies equipped with radar sensors at the front and the rear, and from the pair of front-view SLC images A composite aperture radar image is produced from a front observation InSAR interference map and a rear view InSAR interference map from the rear view SLC image pair, and a MAI interference from the front view InSAR interference map and the rear view InSAR interference map. Speed measuring device of the used object. 제 14항에 있어서,15. The method of claim 14, 상기 제2 간섭도 제작부는, The second interference degree production unit, 상기 한쌍의 전방 관측 SLC 영상 중 하나를 주영상 및 나머지를 부영상으로 결정하고, 상기 한 쌍의 후방 관측 SLC 영상 중 하나를 주영상 및 나머지를 부영상으로 설정하며,One of the pair of forward-viewing SLC images is determined as a main image and the rest as sub-pictures, and one of the pair of backward-viewing SLC images is set as a main-image and the rest as sub-pictures, 상기 부영상의 컬러복소수를 상기 부영상에 대응되는 주영상에 곱함으로써 상기 전방관측 InSAR 간섭도 및 상기 후방관측 InSAR 간섭도를 제작하며,Producing the front observation InSAR interference and the rear observation InSAR interference by multiplying the color complex number of the sub-image by the main image corresponding to the sub-image, 상기 후방관측 InSAR 간섭도의 컬레복소수를 상기 전방관측 InSAR 간섭도에 곱함으로써 상기 MAI 간섭도를 제작하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 장치. And the MAI interference level is produced by multiplying the front complex InSAR interference degree by the COLE complex number of the rearview InSAR interference degree. 제 10항에 있어서,The method of claim 10, 상기 제2 속도 관측부는, The second speed observation unit, 상기 MAI 간섭도로부터 MAI 간섭도의 위상을 결정하며, 상기 MAI 간섭도의 위상, 두 레이더 센서 탑재체가 동일 지역을 촬영하는 시간간격, 비행방향으로의 레이더 안테나 길이 및 정규화된 스퀸스 조정 변수(normalized squint adjustment parameter)를 이용하여 상기 비행방향의 속도를 관측하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 방법. The phase of the MAI interference is determined from the MAI interference degree, the phase of the MAI interference degree, the time interval at which the two radar sensor payloads photograph the same region, the radar antenna length in the flight direction, and the normalized sequence adjustment parameter (normalized) and a squint adjustment parameter to observe the velocity of the flight direction. 제 10항에 있어서,The method of claim 10, 상기 제2 속도 관측부는,, The second speed observation unit, 상기 이동물체가 지표관측 방향으로 이동하고 있는 경우, 상기 지표관측방향으로의 속도를 이용하여 상기 MAI 간섭도에 대한 위상을 보정하기 위한 보정값을 결정하고, When the moving object is moving in the surface observation direction, a correction value for correcting a phase with respect to the MAI interference degree is determined using the speed in the surface observation direction, 상기 MAI 간섭도에서 상기 보정값의 차를 최종 MAI 간섭도로 결정하며, 상기 최종 MAI 간섭도로부터 상기 이동물체에 대한 비행방향의 속도를 관측하는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 장치. The speed difference of the object using the composite aperture radar image is determined by determining the difference between the correction values in the MAI interference degree as the final MAI interference degree and observing the velocity of the flight direction with respect to the moving object from the final MAI interference degree. Device. 제 17항에 있어서,The method of claim 17, 상기 보정값은,The correction value is, 비행방향으로의 레이더 안테나 길이, 상기 지표관측 방향으로의 속도, 정규화된 스퀸스 조정 변수(normalized squint adjustment parameter), 기선의 수평면에 투영된 길이 및 레이더 센서 탑재체의 속도를 이용하여 결정되는 것을 특징으로 하는 합성 개구 레이더 영상을 이용한 물체의 속도 측정 장치. Radar antenna length in the direction of flight, speed in the surface observation direction, normalized squint adjustment parameter, length projected on the horizontal plane of the baseline, and speed of the radar sensor payload. An apparatus for measuring the velocity of an object using a synthetic aperture radar image.
KR1020090117256A 2009-11-30 2009-11-30 The method for measurign object's velocity using synthetic aperture radar image and the apparatus thereof KR101135070B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020090117256A KR101135070B1 (en) 2009-11-30 2009-11-30 The method for measurign object's velocity using synthetic aperture radar image and the apparatus thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020090117256A KR101135070B1 (en) 2009-11-30 2009-11-30 The method for measurign object's velocity using synthetic aperture radar image and the apparatus thereof

Publications (2)

Publication Number Publication Date
KR20110060626A true KR20110060626A (en) 2011-06-08
KR101135070B1 KR101135070B1 (en) 2012-04-13

Family

ID=44395359

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020090117256A KR101135070B1 (en) 2009-11-30 2009-11-30 The method for measurign object's velocity using synthetic aperture radar image and the apparatus thereof

Country Status (1)

Country Link
KR (1) KR101135070B1 (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102955158A (en) * 2012-06-04 2013-03-06 北京航空航天大学 Multiple baseline design method for improving indication performance of ground moving target
CN103728621A (en) * 2014-01-24 2014-04-16 中国科学院电子学研究所 Onboard SAL imaging method adopting along-track interference to inhibit platform vibration
WO2014104550A1 (en) * 2012-12-28 2014-07-03 서울시립대학교 산학협력단 Method for correcting ion distortions of satellite radar interferogram and apparatus therefor
US9753133B2 (en) 2014-05-23 2017-09-05 Agency For Defense Development Apparatus and method for converting multi-channel tracking information for integrated processing of flight data
CN107505614A (en) * 2017-07-20 2017-12-22 上海无线电设备研究所 A kind of object detection method based on the fusion of ATI SAR image width phase information
CN108346292A (en) * 2018-04-17 2018-07-31 吉林大学 City expressway real-time traffic index calculation method based on bayonet data
CN108646238A (en) * 2018-03-06 2018-10-12 中国船舶重工集团公司第七二四研究所 A kind of interference source tracking based on sidelobe cancellation coefficient mapping
US10459061B2 (en) 2014-05-08 2019-10-29 Src, Inc. FFT-based displaced phase center array/along-track interferometry architecture
KR20200122890A (en) * 2019-04-19 2020-10-28 한국전자통신연구원 Apparatus and method for generating a front-view radar image of a vehicle based on a time-domain correlation
KR102211580B1 (en) * 2019-10-01 2021-02-03 엘아이지넥스원 주식회사 Method and system for estimating moving target speed using SAR mounted on multi unmanned aerial vehicles
CN113687352A (en) * 2021-08-05 2021-11-23 南京航空航天大学 Inversion method for down-track interferometric synthetic aperture radar sea surface flow field
CN113805146A (en) * 2021-07-23 2021-12-17 北京空间飞行器总体设计部 High-orbit SAR image space-variant phase error estimation method

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101456185B1 (en) * 2013-05-16 2014-10-31 포항공과대학교 산학협력단 Method and apparatus for yielding radar imaging
KR101605450B1 (en) 2014-08-04 2016-03-22 서울시립대학교산학협력단 Method of stacking multi-temporal MAI interferogram and Apparatus Thereof
KR101873732B1 (en) 2016-03-14 2018-07-03 서울시립대학교 산학협력단 Apparatus and Method for Correction of Distortion on Displacement Map generated from Synthetic Aperture Radar Offset Tracking Method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100529401B1 (en) 2005-07-20 2005-11-17 중앙항업(주) Apparatus and method of dem generation using synthetic aperture radar(sar) data
KR100966289B1 (en) * 2008-02-25 2010-06-28 삼성중공업 주식회사 Target detection method applied in marine radar system, and record carrier, arpa board and radar using the same

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102955158A (en) * 2012-06-04 2013-03-06 北京航空航天大学 Multiple baseline design method for improving indication performance of ground moving target
WO2014104550A1 (en) * 2012-12-28 2014-07-03 서울시립대학교 산학협력단 Method for correcting ion distortions of satellite radar interferogram and apparatus therefor
KR101490981B1 (en) * 2012-12-28 2015-02-09 서울시립대학교 산학협력단 Method for Correction of Ionospheric Distortion of Synthetic Aperture Radar Interferogram and Apparatus Thereof
US9971030B2 (en) 2012-12-28 2018-05-15 University Of Seoul Industry Cooperation Foundation Method and apparatus for correcting ionic distortion of satellite radar interferogram
CN103728621A (en) * 2014-01-24 2014-04-16 中国科学院电子学研究所 Onboard SAL imaging method adopting along-track interference to inhibit platform vibration
US10459061B2 (en) 2014-05-08 2019-10-29 Src, Inc. FFT-based displaced phase center array/along-track interferometry architecture
US9753133B2 (en) 2014-05-23 2017-09-05 Agency For Defense Development Apparatus and method for converting multi-channel tracking information for integrated processing of flight data
CN107505614A (en) * 2017-07-20 2017-12-22 上海无线电设备研究所 A kind of object detection method based on the fusion of ATI SAR image width phase information
CN107505614B (en) * 2017-07-20 2020-07-31 上海无线电设备研究所 Target detection method based on ATI-SAR image amplitude and phase information fusion
CN108646238A (en) * 2018-03-06 2018-10-12 中国船舶重工集团公司第七二四研究所 A kind of interference source tracking based on sidelobe cancellation coefficient mapping
CN108346292A (en) * 2018-04-17 2018-07-31 吉林大学 City expressway real-time traffic index calculation method based on bayonet data
KR20200122890A (en) * 2019-04-19 2020-10-28 한국전자통신연구원 Apparatus and method for generating a front-view radar image of a vehicle based on a time-domain correlation
KR102211580B1 (en) * 2019-10-01 2021-02-03 엘아이지넥스원 주식회사 Method and system for estimating moving target speed using SAR mounted on multi unmanned aerial vehicles
CN113805146A (en) * 2021-07-23 2021-12-17 北京空间飞行器总体设计部 High-orbit SAR image space-variant phase error estimation method
CN113805146B (en) * 2021-07-23 2024-04-02 北京空间飞行器总体设计部 Space-variant phase error estimation method for high-orbit SAR image
CN113687352A (en) * 2021-08-05 2021-11-23 南京航空航天大学 Inversion method for down-track interferometric synthetic aperture radar sea surface flow field

Also Published As

Publication number Publication date
KR101135070B1 (en) 2012-04-13

Similar Documents

Publication Publication Date Title
KR101135070B1 (en) The method for measurign object's velocity using synthetic aperture radar image and the apparatus thereof
Zebker et al. Topographic mapping from interferometric synthetic aperture radar observations
US9417323B2 (en) SAR point cloud generation system
KR101111689B1 (en) The method for three-dimensional deformation measurement and the apparatus thereof
CN102967859B (en) Forward-looking scanning radar imaging method
CN108776342B (en) High-speed platform SAR slow moving target detection and speed estimation method
CN110007302B (en) Spaceborne double-antenna squint in-orbit interference SAR ocean current velocity measurement method
US8159387B1 (en) Multi-transmitter interferometry
EP2413158A1 (en) A method for monitoring terrain and man-made feature displacements using ground-based synthetic aperture radar (GBSAR) data
JP2003500658A (en) Procedures for radar measurements of urban area and landslide motion.
CN105929399A (en) Interference SAR data imaging and elevation estimation method
CN103630905B (en) The overlapping sub-aperture imaging method of array antenna SAR polar coordinates
CN103412310A (en) Bistatic forward-looking synthetic aperture radar ground moving target detecting method and imaging method
CN103744081B (en) A kind of airborne circular track synthetic aperture radar high-precision three-dimensional imaging compensating method
CN113238225B (en) Radar for moving target detection
CN110068817B (en) Terrain mapping method, instrument and system based on laser ranging and InSAR
CN108226925A (en) A kind of overlapped sub-aperture algorithm suitable for missile-borne strabismus time-varying parameter SAR before big
CN110554377B (en) Single-channel SAR two-dimensional flow field inversion method and system based on Doppler center offset
CN103454633A (en) Interference SAR movable base line processing method based on back-projection algorithm
KR102151362B1 (en) Image decoding apparatus based on airborn using polar coordinates transformation and method of decoding image using the same
RU2597148C1 (en) Method of measuring vector field speed of ocean and river streams in space sar
CN104122552B (en) A kind of slidingtype dual station circumferential synthetic aperture radar imaging method
KR101856826B1 (en) A terrain-aided navigation apparatus using a multi-look angle radar altimeter
US8723718B2 (en) Method and device for determining aspect angle progression
CN107728144B (en) Interference SAR imaging method based on forward-looking double-basis mode

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: 20160302

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20170324

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20180323

Year of fee payment: 7

FPAY Annual fee payment

Payment date: 20190329

Year of fee payment: 8