KR100252345B1 - Distortion complement method of radar system - Google Patents

Distortion complement method of radar system Download PDF

Info

Publication number
KR100252345B1
KR100252345B1 KR1019970038540A KR19970038540A KR100252345B1 KR 100252345 B1 KR100252345 B1 KR 100252345B1 KR 1019970038540 A KR1019970038540 A KR 1019970038540A KR 19970038540 A KR19970038540 A KR 19970038540A KR 100252345 B1 KR100252345 B1 KR 100252345B1
Authority
KR
South Korea
Prior art keywords
ground
point
control point
distance
ground control
Prior art date
Application number
KR1019970038540A
Other languages
Korean (ko)
Other versions
KR19990016091A (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 KR1019970038540A priority Critical patent/KR100252345B1/en
Publication of KR19990016091A publication Critical patent/KR19990016091A/en
Application granted granted Critical
Publication of KR100252345B1 publication Critical patent/KR100252345B1/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
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9027Pattern recognition for feature extraction
    • 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
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric techniques
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

PURPOSE: A method of correcting topographical distortion in a synthetic aperture radar system is provided which corrects the topographical distortion due to foreshortening with regard to reflectivity of a corresponding point when a video signal is reconstructed in the SAR system. CONSTITUTION: The position of a ground range image is calculated using the start point and end point of the shortest range image with respect to a synthetic aperture radar system when a signal obtained by correcting RF measurement distortion for a target is applied(301). The height of a predetermined ground control point within the ground range image is calculated with reference to a digital elevation map of a corresponding area(302). The shortest range corresponding to the ground control point is calculated using the distance between the actual position of the ground control point in the ground range and the synthetic aperture radar system(303). The reflectivity of the ground control point is calculated for a pixel placed in the shortest range(304). A value obtained by multiplying the reflectivity by the pixel value located in the shortest range is mapped into the pixel value of the ground control point(305). The step of calculating the height information to the mapping step are repeated until the mapping of the ground range image for the obtained shortest range image is completed(306).

Description

합성개구 레이다시스템에 있어서 지형적인 왜곡 보정방법Terrain Distortion Correction Method in Synthetic Opening Radar System

본 발명은 합성개구 레이다(Synthetic Aperture Radar, 이하 SAR이라 약함)시스템에 있어서 지형적인 왜곡을 보정하는 방법에 관한 것으로, 특히 포쇼트닝(Foreshortning)에 의한 SAR영상의 지형적인 왜곡을 보정하기 위한 지형적인 왜곡보정방법에 관한 것이다.The present invention relates to a method of correcting topographical distortion in a Synthetic Aperture Radar (hereinafter referred to as SAR) system. It relates to a distortion correction method.

SAR시스템은 거리에 관계없이 일정한 방위각 해상도를 유지할 수 있으며, 초고주파를 사용하여 영상을 얻어내기 때문에 가시광의 유무와 날씨, 구름 등에 영향을 받지 않고 소기의 지역에 대한 정보를 얻어낼 수 있는 특징이 있다. 특히 도플러 효과를 이용하기 때문에 기존의 SLAR(Side Looking Aperture Radar)시스템보다 방위각방향으로 훨씬 더 높은 고해상도를 실현한 시스템이다.SAR system can maintain constant azimuth resolution regardless of distance, and because it acquires image using ultra high frequency, it can get information about desired area without being affected by visible light, weather and cloud. . In particular, the Doppler effect is used to achieve much higher resolution in the azimuth direction than the existing SLAR (Side Looking Aperture Radar) system.

이러한 SAR시스템은 목표물에 대해 도 1에 도시된 바와 같이 최단 거리(Slant Range)에 근거하여 영상을 획득하기 때문에 SAR시스템의 높이 H에 의해 최단거리 영상은 지상거리(Ground Range)영상과 다른 값을 가지게 된다. 특히, 도 1에 도시된 바와같이 SAR시스템 방향으로 α와 같은 경사면이 존재하는 경우에, α경사면에 대한 지상거리상에서의 영상의 길이 a보다 최단 거리 영상의 길이 a'가 작게 표현된다. 이러한 현상을 포쇼트닝이라 한다.Since the SAR system acquires the image based on the shortest range (Slant Range) for the target as shown in FIG. 1, the shortest distance image is different from the ground range image by the height H of the SAR system. Have. In particular, when there is an inclined surface such as α in the SAR system direction as shown in FIG. 1, the length a 'of the shortest distance image is smaller than the length a of the image on the ground distance to the α inclined surface. This phenomenon is called four shortening.

이러한 포쇼트닝에 따른 지형적인 왜곡을 보정하기 위하여 기존의 SAR시스템에서는 사전에 취득된 그 지역에 대한 표본화된 높이 정보로 이루어진 DEM(Digital Elevation Map, 이하 DEM이라 약함)을 이용하여 최단거리 영상을 지상거리 영상으로 매핑한다. 예를 들어 지상거리 영상에서의 b지점에 대해 DEM을 참조하여 해당 높이 h를 얻어낸다. 그리고 얻어낸 높이 h를 이용하여 최단거리 b'를 계산해 낸다. 그리고 계산된 b'지점의 최단 거리 영상의 화소값을 지상거리 영상의 화소값 b에 대입하도록 한다.In order to correct the topographical distortion due to the shortening, the existing SAR system uses the digital elevation map (DEM), which is composed of sampled height information of the area previously acquired, to ground the shortest distance image. Map to distance image. For example, for the point b in the ground distance image, the height h is obtained by referring to the DEM. The shortest distance b 'is calculated using the height h obtained. The calculated pixel value of the shortest distance image at point b 'is substituted into the pixel value b of the ground distance image.

그러나 도 1에 도시된 α와 같은 경사면의 경우 상술한 바와 같이 최단 거리상에서의 거리 a'보다 지상 거리상에서의 거리 a가 크기 때문에 지상거리상의 몇 개의 점들이 최단거리 영상의 한점에 모여 최단 거리상의 해당 화소값을 크게 할 가능성이 높다. 즉, 지상거리상의 b지점외에 경사면에 존재하는 인접한 지점에 대해 상술한 DEM을 참조하여 해당 최단거리를 계산한 결과, b'지점에 해당되는 경우에, 상술한 b지점을 비롯한 인접한 지점의 반사된 신호에 의해 최단거리상의 b' 지점의 화소값이 결정되어 b'지점의 화소값은 큰 값을 갖게 되고, 이러한 b'지점의 큰 화소값이 상술한 매핑과정에 의해 지상거리상의 b지점을 비롯한 인접한 지점의 화소값으로 매핑되어, 지상거리 영상의 b지점을 비롯한 인접한 점들의 화소값은 원래의 화소값보다 더 큰 값으로 표현되는 지형적인 왜곡이 발생된다.However, in the case of the inclined surface such as α shown in FIG. 1, as described above, since the distance a on the ground distance is larger than the distance a 'on the shortest distance, several points on the ground distance are gathered at one point of the shortest distance image and thus on the shortest distance. It is highly possible to increase the pixel value. That is, when the shortest distance is calculated with reference to the above-described DEM for the adjacent points existing on the inclined plane in addition to the b point on the ground distance, and when the b 'point corresponds to the b' point, the reflected point of the adjacent point including the above-described b point is reflected. The pixel value at point b 'on the shortest distance is determined by the signal, and the pixel value at point b' has a large value, and the large pixel value at point b 'includes the point b on the ground distance by the above-described mapping process. Mapping to pixel values of adjacent points, topographical distortion is generated in which the pixel values of adjacent points, including point b of the ground-distance image, are represented by larger values than the original pixel values.

본 발명은 상술한 결점을 개선하기 위하여 안출한 것으로서, SAR시스템에 있어서 영상신호 재구성시 포쇼트닝에 의한 지형적인 왜곡을 해당 지점의 반사기여도를 고려하여 보정하기 위한 지형적인 왜곡 보정방법을 제공하는데 그 목적이 있다.The present invention has been made to solve the above-described drawbacks, and provides a topographic distortion correction method for correcting topographical distortion due to shortening when reconstructing an image signal in a SAR system in consideration of reflectivity of a corresponding point. There is a purpose.

상기 목적을 달성하기 위하여 본 발명에 따른 지형적인 왜곡보정방법은, 소정의 목표물에 대해 획득된 최단거리영상을 지상거리영상으로 매핑하여 영상을 획득하는 합성개구 레이다시스템에서 지형적인 왜곡을 보정하는 방법에 있어서, 목표물에 대한 고주파측정 왜곡을 보정한 신호가 인가되면, 합성개구 레이다시스템에 대한 최단거리영상의 시작점과 끝점을 이용하여 지상거리영상의 위치를 계산하는 단계; 지상거리영상내의 소정의 지상제어점의 높이를 해당 지역의 표본화된 높이정보로 이루어진 맵(DEM)을 참조하여 계산하는 단계; 높이계산단계에서 획득된 소정 지상제어점의 높이를 이용하여 소정 지상제어점의 지상거리상에서의 실질적인 위치와 합성개구 레이다시스템간의 거리를 이용하여 소정의 지상제어점에 대응되는 최단거리를 계산하는 단계; 대응되는 최단거리에 위치한 화소에 대해 소정 지상제어점의 반사기여도를 계산하는 단계; 반사기여도와 대응되는 최단거리에 위치한 화소값을 곱하여 얻어진 값을 소정 지상제어점의 화소값으로 매핑하는 단계; 획득된 최단거리영상에 대한 지상거리영상의 매핑이 모두 완료될 때까지 높이정보를 계산하는 단계 내지 매핑하는 단계를 반복 수행하는 단계를 포함하는 포함하여 수행되는 것을 특징으로 한다.In order to achieve the above object, the topographical distortion correction method according to the present invention is a method for correcting topographical distortion in a synthetic opening radar system for obtaining an image by mapping the shortest distance image acquired for a predetermined target to a ground distance image. The method of claim 1, further comprising: calculating a position of the ground distance image using a start point and an end point of the shortest distance image of the synthesized opening radar system when a signal corrected for the high frequency measurement distortion of the target is applied; Calculating a height of a predetermined ground control point in the ground distance image by referring to a map (DEM) of sampled height information of a corresponding region; Calculating a shortest distance corresponding to a predetermined ground control point by using the height of the predetermined ground control point obtained in the height calculation step using a distance between the actual position on the ground distance of the predetermined ground control point and the synthetic opening radar system; Calculating reflectivity of a predetermined ground control point with respect to the pixel located at the shortest distance; Mapping a value obtained by multiplying a reflectance degree by a pixel value corresponding to the shortest distance to a pixel value of a predetermined ground control point; Comprising the step of repeating the step of calculating the height information until the mapping of the ground distance image to the obtained shortest distance image is completed, characterized in that it is performed.

도 1은 기존의 합성개구 레이다시스템을 통해 얻어지는 영상의 기하학적인 모델 예시도,1 is a view illustrating a geometric model of an image obtained through a conventional synthetic opening radar system,

도 2는 일반적인 합성개구 레이다시스템의 개략적인 블럭도,2 is a schematic block diagram of a general synthetic opening radar system;

도 3은 합성개구 레이다 시스템에 적용되는 본 발명에 따른 지형적인 왜곡보정방법의 동작 흐름도,3 is an operation flowchart of a topographic distortion correction method according to the present invention applied to a synthetic opening radar system;

도 4는 도 3의 동작 흐름도에 따른 동작을 설명하기 위한 합성개구 레이다시스템을 통해 얻어지는 영상의 기하학적인 모델 예시도.4 is an exemplary view of a geometric model of an image obtained through a synthetic opening radar system for explaining the operation according to the operation flowchart of FIG.

<도면의 주요부분에 대한 부호의 설명><Description of the symbols for the main parts of the drawings>

200:원(Raw)신호 제공부 210:압축신호 처리부200: raw signal providing unit 210: compressed signal processing unit

220:고주파 측정(Radiometric) 왜곡 보정부 230:지형적 왜곡 보정부220: high-frequency measurement distortion correction unit 230: terrain distortion correction unit

본 발명의 상술한 목적 및 기타 목적과 여러 가지 장점은 이 기술 분야에 숙련된 사람들에 의해 첨부된 도면을 참조하여 하기에 기술되는 발명의 바람직한 실시예로부터 더욱 명확하게 될 것이다.The above and other objects and various advantages of the present invention will become more apparent from the preferred embodiments of the invention described below with reference to the accompanying drawings by those skilled in the art.

이하, 첨부된 도면을 참조하여 본 발명에 따른 바람직한 실시예를 상세하게 설명하기로 한다.Hereinafter, exemplary embodiments of the present invention will be described in detail with reference to the accompanying drawings.

도 2는 일반적인 합성개구 레이다(이하 SAR이라 약함)시스템의 블럭도로서, 원하는 지형으로 발사된 초고주파의 반사되는 신호를 수신하여 소정의 원(Raw)신호를 생성하여 제공하는 원신호 제공부(200), 원신호 제공부(200)로부터 제공되는 원신호에 대한 거리방향 압축, 방위각방향 압축처리를하여 해상도를 향상시키는 압축 신호처리부(210), 압축 신호처리부(210)에서 압축된 원신호에 대한 출력변화, 수신단 이득변환, 안테나 패턴 및 거리에 따른 손실과 같은 고주파측정 왜곡에 대한 보정처리를하는 고주파 측정(Radiometric) 왜곡 보정부(220), 고주파 측정 왜곡 보정부(220)에서 보정되어 출력되는 원신호의 지형적인 왜곡을 본 발명에 따라 보정하기 위한 지형적 왜곡 보정부(230)로 구성된다.FIG. 2 is a block diagram of a general synthetic opening radar system (hereinafter, referred to as SAR) system, and includes: an original signal providing unit 200 for generating and providing a predetermined raw signal by receiving a reflected signal of a very high frequency emitted to a desired terrain; ), The compression signal processing unit 210 to improve the resolution by performing distance-direction compression and azimuth compression processing on the original signal provided from the original signal providing unit 200, and to the original signal compressed by the compression signal processing unit 210. The high frequency measurement distortion correction unit 220 and the high frequency measurement distortion correction unit 220 correcting and outputting the high frequency measurement distortion such as the output change, the receiver gain conversion, the antenna pattern and the loss according to the distance It is composed of a topographic distortion correction unit 230 for correcting the topographical distortion of the original signal according to the present invention.

도 3은 본 발명에 따른 지형적인 왜곡 보정방법에 대한 동작 흐름도이고, 도 4는 본 발명에 따른 지형적인 왜곡보정을 설명하기 위한 SAR영상의 기하학적인 모델 예시도이다.3 is a flowchart illustrating a topographic distortion correction method according to the present invention, and FIG. 4 is a diagram illustrating a geometric model of a SAR image for explaining the topographic distortion correction according to the present invention.

그러면 도 2 및 도 4를 참조하여 도 3에 도시된 본 발명에 따른 실시예를 상세하게 설명하기로 한다.Next, the embodiment according to the present invention shown in FIG. 3 will be described in detail with reference to FIGS. 2 and 4.

우선, 원신호 제공부(200)와 압축신호 처리부(210) 및 고주파 측정 왜곡 보정부(220)는 기존의 SAR시스템에서 이루어진 바와 동일하므로 여기서 이에 대한 구체적인 설명은 생략하고, 지형적인 왜곡 보정부(230)의 동작을 중심으로 설명하기로 한다.First, since the original signal providing unit 200, the compressed signal processing unit 210, and the high frequency measurement distortion correction unit 220 are the same as those of the conventional SAR system, a detailed description thereof will be omitted, and the topographic distortion correction unit ( The operation of 230 will be described below.

즉, 지형적 왜곡 보정부(230)는 고주파측정 왜곡 보정부(220)로부터 고주파 왜곡이 보정된 영상이 전송되면, 제 301 단계로 진행되어 지상거리의 영상의 위치를 계산한다. 지상거리의 영상의 위치는 최단거리의 영상의 시작점과 끝점을 이용하여 계산한다. 즉, SAR시스템으로부터 최단거리상의 영상의 시작점간의 거리와 대응되는 SAR시스템과 지상거리상의 거리를 갖는 지상거리상의 지점을 지상거리 영상의 시작점으로 하고 최단거리상의 영상의 끝점과 일치하는 지상거리상의 지점을 지상거리 영상의 끝지점으로 한다.That is, when the image distortion-corrected image is transmitted from the high frequency measurement distortion correction unit 220, the topographic distortion correction unit 230 proceeds to step 301 to calculate the position of the image of the ground distance. The position of the image of the ground distance is calculated using the start point and the end point of the image of the shortest distance. That is, the point on the ground distance that corresponds to the distance between the start point of the image on the shortest distance from the SAR system and the point on the ground distance having the distance on the ground distance and the corresponding SAR system as the start point of the ground distance image, and the point on the ground distance that matches the end point of the image on the shortest distance. Is the end point of the ground distance image.

이와 같이 지상거리 영상에 대한 위치가 계산되면, 지형적 왜곡 보정부(230)는 제 302 단계로 진행되어 현 지상제어점(Ground Control Point)의 높이를 계산한다. 여기서 현 지상제어점은 지상거리영상의 시작점부터 끝지점중 임의의 한 지점으로 순차적으로 선택된다. 예를 들어 현 지상제어점이 도 4에 도시된 바와 같이 c지점인 경우에 지형적 왜곡 보정부(230)는 h'에 해당되는 높이를 계산하게 된다. 높이는 상술한 바와 같이 사전에 획득된 해당 지역에 대한 표본화된 높이정보로 구성된 DEM을 참조한다.When the position of the ground distance image is calculated as described above, the terrain distortion correction unit 230 proceeds to step 302 to calculate the height of the current ground control point. The current ground control point is sequentially selected to any one of the starting point and the ending point of the ground distance image. For example, when the current ground control point is point c as shown in FIG. 4, the topographic distortion correction unit 230 calculates a height corresponding to h ′. As described above, the height refers to a DEM composed of sampled height information about a corresponding region acquired in advance.

즉, DEM은 2차원의 맵의 형태로 해당 지형에 대한 표본화된 높이정보를 보유하고 있으므로 지상거리 영상상에서 현 지상제어점의 2차원 좌표계를 이용하여 해당 DEM정보를 액세스하게 된다. 이 때, 해당 좌표지점에 대한 DEM정보가 존재하지 않는 경우에는 인접한 높이정보들의 인터폴레이션(Interpolation)에 의해 얻어진 높이정보를 현 지상제어점에서의 높이 h'로 획득한다. 해당 높이 h'가 획득되면, 지형적 왜곡 보정부(230)는 제 303 단계로 진행된다.That is, since the DEM holds sampled height information on the terrain in the form of a two-dimensional map, the DEM information is accessed using the two-dimensional coordinate system of the current ground control point on the ground distance image. At this time, if there is no DEM information for the corresponding coordinate point, the height information obtained by interpolation of adjacent height information is obtained as the height h 'at the current ground control point. When the height h 'is obtained, the topographical distortion correction unit 230 proceeds to step 303.

제 303 단계에서 지형적 왜곡 보정부(230)는 해당되는 최단거리 RA'를 계산한다. 이 최단거리 계산은 c지점의 구해진 높이 h'를 이용하여 얻어진 SAR시스템과 실질적인 지상제어점간의 거리(RA)계산에 의해 얻어진다. 즉, 거리(RA)는 최단거리 RA'와 동일하기 때문에 지상거리상에서 높이 h'인 지점과 SAR시스템간의 거리가 구해지면 SAR시스템와 최단거리간의 거리로 대입하면 된다. 이러한 거리계산은 최단거리상의 해당 화소의 값( Pc′ )을 획득하기 위함이다. 그러나 이 때, 해당 지점에 화소값이 존재하지 않는 경우에는 인접한 위치에 높인 화소들의 인터폴레이션에 의해 해당 화소값을 결정한다.In step 303, the topographic distortion correction unit 230 calculates a corresponding shortest distance R A ′. This shortest distance calculation is obtained by calculating the distance R A between the SAR system and the actual ground control point obtained using the height h 'obtained at point c. In other words, since the distance R A is the same as the shortest distance R A ', the distance between the SAR system and the shortest distance is obtained when the distance between the point h' on the ground distance and the SAR system is obtained. This distance calculation is based on the value of the corresponding pixel on the shortest distance ( P c ′ To obtain). However, in this case, when the pixel value does not exist at the corresponding point, the pixel value is determined by interpolation of the pixels raised at the adjacent position.

그 다음 제 304 단계로 진행되어 최단거리 RA'를 갖는 c'화소값에 대한 c지점의 반사기여도를 계산한다. 이 때, 반사기여도는 최단거리가 동일한 것으로 계산된 지상거리상의 인접한 점들의 반사도를 고려하여 계산된다. 동일한 최단거리상의 화소값을 형성하는 지점 검출은 지상거리 영상의 시작점에서 끝점까지에 존재하는 모든 점들에 대해 제 303단계와 같이 최단거리르 계산하여 지상제어점과 같은 최단거리를 가지는 지점을 찾아내는 방법으로 이루어진다. 즉 도 4의 경우에는 지상거리에서 c, d, e지점이 최단거리상의 c'지점으로 반사되므로 c지점의 c'지점에 대한 반사기여도는 수학식 1과 같이 계산된다.The flow proceeds to step 304 to calculate the reflectivity of the point c with respect to the pixel value c 'having the shortest distance R A '. At this time, the reflectivity is calculated in consideration of the reflectivity of adjacent points on the ground distance calculated as the shortest distance. The point detection that forms the same pixel value on the shortest distance is to find the point with the shortest distance like the ground control point by calculating the shortest distance for all points existing from the start point to the end point of the ground distance image as in step 303. Is done. That is, in the case of FIG. 4, since c, d, and e points are reflected from the ground distance to the c 'point on the shortest distance, the reflectivity of the c' c point of the c point is calculated as in Equation 1.

Figure 1019970038540_B1_M0001
Figure 1019970038540_B1_M0001

이와 같은 반사기여도는 d지점에 대해 계산될 때에는 수학식 1의 분자가 γc 에서 γd 로 변경되어 계산되고, e지점에 대해 계산될 때에는 수학식 1의 분자가 가 γc 에서 γe 로 변경되어 계산된다.Such reflectivity may be calculated by the numerator of γ c in γ d When calculated for point e, the numerator of Equation 1 is γ c in γ e It is changed to and calculated.

c지점에 대한 반사기여도가 계산되면, 제 305 단계로 진행되어 제 303 단계에서 계산된 최단거리상에 위치한 c'화소의 값 Pc′ 와 제 304 단계에서 계산된 반사기여도를 수학식 2와 같이 결합하여 지상거리의 c지점에 매핑할 화소값을 결정한다.When the reflectivity of the point c is calculated, the flow proceeds to step 305 where the value of c'pixel located on the shortest distance calculated in step 303 is obtained. P c ′ And the reflectivity of the reflector calculated in step 304 are combined as shown in Equation 2 to determine a pixel value to be mapped to point c of the ground distance.

Figure 1019970038540_B1_M0002
Figure 1019970038540_B1_M0002

이에 따라 지상거리상의 c지점에 대한 최단거리 영상간의 매핑처리가 완료된다.This completes the mapping process between the shortest-range images of point c on the ground distance.

그러나 전 지상거리상의 모든 픽셀에 대한 최단거리 영상과의 매핑처리를 위해 지형적 왜곡보정부는 제 306 단계로 진행되어 모든 지상거리상의 화소값이 결정되었는 지를 체크하고, 결정되지 않은 경우에는 제 302 단계로 진행되어 상술한 과정을 반복 수행한다. 이 때, 그 다음의 지상제어점은 스캐닝순서에 의해 결정된다. 그리고 제 306 단계의 체크결과, 지상거리영상의 모든 화소값이 결정되면 지형적인 왜곡보정작업을 종료한다.However, in order to process the mapping with the shortest-range image of all the pixels on all ground distances, the terrain distortion correction proceeds to step 306 to check whether the pixel values on all ground distances are determined, and if not, to step 302. It proceeds and repeats the above process. At this time, the next ground control point is determined by the scanning order. As a result of the check in step 306, when all the pixel values of the terrestrial distance image are determined, the terrain distortion correction operation is terminated.

이상에서 설명한 바와 같이 본 발명은 합성개구 레이다시스템에 있어서 포쇼트닝에 의한 지형적인 왜곡을 해당 최단거리의 화소값에 대한 해당 지상거리 지점의 반사기여도를 이용하여 보정함으로써, 거의 원래의 값에 근사한 값을 갖는 SAR영상을 획득할 수 있는 이점이 있다.As described above, the present invention is a value close to the original value by correcting the topographical distortion due to the shortening in the composite opening radar system using the reflectivity of the ground distance point with respect to the pixel value of the shortest distance. There is an advantage that can obtain a SAR image having a.

본 발명은 상술한 실시예로서 설명되었으나 당업자라면 본 발명의 기술사상을 일탈하지 아니하는 범위에서 다양한 변경 및 수정이 가능함을 알 수 있을 것이다. 예를 들어 상술한 실시예에서는 3개의 지상거리상의 지점이 동일한 최단거리를 갖는 경우를 예시하였으나 동일한 최단거리를 갖는 지점이 3개 미만이거나 3개 이상이 존재하는 경우에는 동일한 방식으로 반사기여도를 계산할 수 있다. 따라서 본 발명의 기술적 범위는 명세서의 상세한 설명에 기재된 내용으로 한정되는 것이 아니라 특허청구의 범위에 의하여 정하여져야 한다.Although the present invention has been described as the above-described embodiment, those skilled in the art will appreciate that various changes and modifications can be made without departing from the spirit of the present invention. For example, the above-described embodiment illustrates the case where three points on the ground distance have the same shortest distance, but when there are less than three points or three or more points having the same shortest distance, the reflectivity is calculated in the same manner. Can be. Therefore, the technical scope of the present invention should not be limited to the contents described in the detailed description of the specification but should be defined by the claims.

Claims (2)

소정의 목표물에 대해 획득된 최단거리영상을 지상거리영상으로 매핑하여 영상을 획득하는 합성개구 레이다시스템에서 지형적인 왜곡을 보정하는 방법에 있어서,In the method for correcting the topographical distortion in a synthetic opening radar system for obtaining an image by mapping the shortest distance image obtained for a predetermined target to a ground distance image, 상기 목표물에 대한 고주파측정 왜곡을 보정한 신호가 인가되면, 상기 합성개구 레이다시스템에 대한 상기 최단거리영상의 시작점과 끝점을 이용하여 지상거리영상의 위치를 계산하는 단계;Calculating a position of the ground distance image by using a start point and an end point of the shortest distance image with respect to the synthesis opening radar system when a signal corrected for the high frequency measurement distortion of the target is applied; 상기 지상거리영상내의 소정의 지상제어점의 높이를 해당 지역의 표본화된 높이정보로 이루어진 맵(DEM)을 참조하여 계산하는 단계;Calculating a height of a predetermined ground control point in the ground distance image with reference to a map (DEM) of sampled height information of a corresponding region; 상기 높이계산단계에서 획득된 상기 소정 지상제어점의 높이를 이용하여 상기 소정 지상제어점의 지상거리상에서의 실질적인 위치와 상기 합성개구 레이다시스템간의 거리를 이용하여 상기 소정의 지상제어점에 대응되는 최단거리를 계산하는 단계;Using the height of the predetermined ground control point obtained in the height calculation step, the shortest distance corresponding to the predetermined ground control point is calculated using the actual position on the ground distance of the predetermined ground control point and the distance between the synthetic opening radar system. Doing; 상기 대응되는 최단거리에 위치한 화소에 대해 상기 소정 지상제어점의 반사기여도를 계산하는 단계;Calculating a reflectivity of the predetermined ground control point for the pixel located at the corresponding shortest distance; 상기 반사기여도와 상기 대응되는 최단거리에 위치한 화소값을 곱하여 얻어진 값을 상기 소정 지상제어점의 화소값으로 매핑하는 단계;Mapping a value obtained by multiplying the reflectivity by the pixel value located at the corresponding shortest distance to the pixel value of the predetermined ground control point; 상기 획득된 최단거리영상에 대한 상기 지상거리영상의 매핑이 모두 완료될 때까지 상기 높이정보를 계산하는 단계 내지 상기 매핑하는 단계를 반복 수행하는 단계를 포함하여 수행되는 것을 특징으로 하는 합성개구 레이다시스템에 있어서 지형적인 왜곡 보정방법.Comprising the step of repeating the step of calculating the height information until the mapping of the ground distance image to the obtained shortest distance image is completed, characterized in that it is performed Terrain Distortion Correction Method 제 1 항에 있어서, 상기 반사기여도를 계산하는 단계는 상기 지상거리상에서 상기 소정 지상제어점과 동일한 최단거리를 갖는 지점의 반사도( γn )와 상기 소정 지상제어점의 반사도( γc )를 하기 수학식과 같이 고려하여 계산하는 것을 특징으로 하는 합성개구 레이다시스템에 있어서 지형적인 왜곡 보정방법.The method of claim 1, wherein the calculating of the reflectivity of the reflector comprises: reflectance of the point having the shortest distance equal to the predetermined ground control point on the ground distance; γ n ) And the reflectivity of the predetermined ground control point ( γ c Topographic distortion correction method for a synthetic opening radar system, characterized in that it is calculated by considering the following equation.
Figure 1019970038540_B1_M0003
Figure 1019970038540_B1_M0003
KR1019970038540A 1997-08-13 1997-08-13 Distortion complement method of radar system KR100252345B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1019970038540A KR100252345B1 (en) 1997-08-13 1997-08-13 Distortion complement method of radar system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1019970038540A KR100252345B1 (en) 1997-08-13 1997-08-13 Distortion complement method of radar system

Publications (2)

Publication Number Publication Date
KR19990016091A KR19990016091A (en) 1999-03-05
KR100252345B1 true KR100252345B1 (en) 2000-04-15

Family

ID=19517375

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1019970038540A KR100252345B1 (en) 1997-08-13 1997-08-13 Distortion complement method of radar system

Country Status (1)

Country Link
KR (1) KR100252345B1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102449672B1 (en) * 2021-03-25 2022-09-30 한국항공우주연구원 Method and system for correcting sar images
WO2023004628A1 (en) * 2021-07-28 2023-02-02 深圳市速腾聚创科技有限公司 Reflectivity correction method and apparatus, computer readable storage medium, and terminal device

Also Published As

Publication number Publication date
KR19990016091A (en) 1999-03-05

Similar Documents

Publication Publication Date Title
Filipponi Sentinel-1 GRD preprocessing workflow
CN107024683B (en) DEM-based SAR image direct orthorectification method and system
US7551119B1 (en) Flight path-driven mitigation of wavefront curvature effects in SAR images
US7532150B1 (en) Restoration of signal to noise and spatial aperture in squint angles range migration algorithm for SAR
KR101757885B1 (en) Apparatus and method for compensating SAR image
CN109597074B (en) SAR image geometric positioning parameter correction method and system
CN112882030B (en) InSAR imaging interference integrated processing method
US11818361B2 (en) Data compression apparatus, data compression method, and program
CN111260542A (en) SAR image splicing method based on sub-block registration
CN114137519A (en) High-resolution SAR imaging parameter calculation method
CN110703252B (en) Digital elevation model correction method for interferometric synthetic aperture radar shadow area
CN115951350A (en) Permanent scatterer point extraction method, device, equipment and medium
KR100252345B1 (en) Distortion complement method of radar system
Roth et al. Geocoding of TerraSAR-X data
US9035820B2 (en) Measurement device, measurement system, measurement method, and program
Ferretti et al. Multibaseline phase unwrapping for InSAR topography estimation
Frey et al. Processing SAR data of rugged terrain by time-domain back-projection
CN111220981B (en) Medium-orbit satellite-borne SAR imaging method based on non-orthogonal non-linear coordinate system output
Parcharidis et al. Landsat TM and ERS data fusion: a statistical approach evaluation for four different methods
KR100252344B1 (en) High frequency measurement distortion complement method of radar system
Sefercik et al. Advanced analysis of differences between C and X bands using SRTM data for mountainous topography
US5485162A (en) Method of performing beam compression process on antenna pattern of radar
KR102666369B1 (en) Method for calculating power level reflectivity of object on ground using SAR image
Bezvesilniy et al. Multi-look radiometric correction of SAR images
US11425324B2 (en) Time-of-flight down-up sampling using a compressed guide

Legal Events

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

Payment date: 20110103

Year of fee payment: 12

LAPS Lapse due to unpaid annual fee