WO2005098473A1 - A real-imaginary correlation method for generating a speckless interferogram of sar - Google Patents

A real-imaginary correlation method for generating a speckless interferogram of sar Download PDF

Info

Publication number
WO2005098473A1
WO2005098473A1 PCT/CN2005/000346 CN2005000346W WO2005098473A1 WO 2005098473 A1 WO2005098473 A1 WO 2005098473A1 CN 2005000346 W CN2005000346 W CN 2005000346W WO 2005098473 A1 WO2005098473 A1 WO 2005098473A1
Authority
WO
WIPO (PCT)
Prior art keywords
window
real
complex image
complex
image data
Prior art date
Application number
PCT/CN2005/000346
Other languages
French (fr)
Chinese (zh)
Inventor
Qifeng Yu
Sihua Fu
Original Assignee
National University Of Defense Technology Of People's Liberation Army
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 National University Of Defense Technology Of People's Liberation Army filed Critical National University Of Defense Technology Of People's Liberation Army
Publication of WO2005098473A1 publication Critical patent/WO2005098473A1/en

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • G01S13/9023—SAR image post-processing techniques combined with interferometric techniques

Definitions

  • the present invention relates to the fields of image processing technology, remote sensing measurement and the like, and further refers to a method for correlating real and imaginary parts of a fringe contour window and a rectangle and other windows for generating an interference synthesis aperture radar interference phase map.
  • Interferometric Synthetic Aperture Radar (Abbreviation: InSAR or IFSAR) three-dimensional imaging technology is a newly developed space observation technology. It is a traditional Synthetic Aperture Radar (abbreviation: SM) remote sensing technology. Combining with radio astronomy interference technology. Sensing Synthetic Aperture Radar technology uses two antennas to observe at the same time or one antenna to perform two parallel observations to receive echo signals from ground targets. After imaging processing, two coherent synthetic aperture radar images (including (Intensity information and phase information), and then the interference phase images of the two complex images are obtained after interference processing, and the three-dimensional information of the observation target is obtained by transforming according to a certain geometric relationship. At present, the generation of interference phase maps mainly adopts the method of conjugate multiplication. Let ⁇ ( ⁇ , ⁇ ) and V 2 (r, x) be the complex values of the two complex images after exact registration, respectively, iJ:
  • V x (?% X) - ⁇ ⁇ ⁇ ⁇ A cos ⁇ + iA x sin ⁇ ⁇ (1)
  • a and E are the echo phases of the two antennas respectively
  • Middle * indicates complex conjugate. From equation (3), we get-
  • Figure 1 shows the interference phase map obtained by this method. Because of interferometric synthetic aperture radar systems The influence of speckle noise, spatial miscorrelation, temporal miscorrelation, and data processing noise makes the signal to noise of the interference phase map relatively low. The existence of such high noise has seriously affected the phase unwrapping and the acquisition of high-precision digital elevation maps, which has become a major bottleneck in the interference synthesis? L-path radar data processing.
  • the technical problem to be solved by the present invention is that, in view of the defects existing in the prior art, in the process of the interference synthetic aperture radar data processing, the method of correlating the real and imaginary parts of the window of the fringe contour and the method of correlating the real and imaginary parts of the window such as a rectangle is used instead of the traditional
  • the conjugate multiplication method generates an interference phase map that is completely free of coherent speckle noise from the complex image after accurate registration of the two frames. This greatly improves the reliability of phase unwrapping and the accuracy of generating digital elevation maps, and solves a major difficulty in the processing of interferometric synthetic aperture radar data.
  • Correlation operation is performed on the real part of the first complex image data and the imaginary part of the second complex image data by the same method as above to find the correlation coefficient C 2 , or the imaginary part of the first complex image data Perform a correlation operation with the real part of the second complex image data to find the number of mesh relations-C 2 , and then obtain C 2;
  • the window may be a stripe contour window, or may be a rectangular or other window.
  • the correlation formula used can be a variety of mathematical expressions such as covariance (mean normalized) correlation, standardized covariance correlation, and so on.
  • covariance mean normalized correlation
  • standardized covariance correlation and so on.
  • the following uses the standardized covariance Mumuguan formula as an example to derive.
  • the standardized covariance correlation formula is as follows:
  • ⁇ ->, IM1 operator means to average a certain variable in the range of mxn pixels.
  • Equation (4) sxct ⁇ C 2 lC x ) (19)
  • Equation (4) sxct ⁇ C 2 lC x ) (19)
  • Equation (19) By comparing Equation (4) and Equation (19), it can be seen that the phase distributions of the phase diagrams obtained by the two are completely consistent.
  • the correlation operation between the imaginary part of the first picture and the imaginary part of the second picture can also be obtained (, and the correlation operation between the imaginary part of the first picture and the real part of the second picture can be obtained -c 2 To obtain C 2 .
  • the optimal window for correlation operations is the fringe contour window, which can be used to obtain the interference phase map with the best quality and accuracy. You can also choose rectangular and other forms of windows, and still get better approximate results.
  • the fringe contour window is a curve window. It is a window formed by taking a certain length and width along the fringe direction of the interference phase map. It can be obtained from the interference phase map generated by the conjugate multiplication method, or from a rectangle. The interference phase map generated by the window correlation method is obtained.
  • the stripe contour window can be obtained by using different algorithms, which is not the content of the present invention and will not be repeated here.
  • the present invention is a method for generating a real-imaginary part correlation method of an interference-free spot synthetic aperture radar interference phase map.
  • the method of speckle interference synthetic aperture radar interferometry; in the method of the present invention, in the process of the data processing of the interference synthetic aperture radar, the method of correlating the real and imaginary parts of the fringe contour window and the method of correlating the real and imaginary parts of the rectangle and other windows is used instead of the conventional common.
  • the yoke multiplication method generates interference phase maps that are completely free of coherent speckle noise from two accurately registered complex images. This greatly improves the reliability of phase unwrapping and the accuracy of generating digital elevation maps, and solves a major difficulty in data processing of interferometric synthetic aperture radar.
  • Figure 1 is an interference phase map generated by conjugate multiplication of two complex images after accurate registration.
  • the image size is 2048 X 2048 pixels.
  • Figure 2 is the interference phase map generated by the rectangular window correlation method for the two complex images after accurate registration, and the size of the rectangular window used is 19 ⁇ 19.
  • Figure 3 is the interference phase generated by the stripe contour window correlation method for the two complex images after accurate registration.
  • the stripe contour window size used is 41 X 5.

Landscapes

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

Abstract

A Real-Imaginary correlation method for generating a speckless interferogram of Synthetic Aperture Radar (SAR) comprises the following steps: applying a window on each of two complex images during the Interferometric Synthetic Aperture Radar (INSAR) data processing, wherein the window is a mxn (m>0, n>0) isoline-window or any other suitable window and is centered on each corresponding pixel in the two complex images; correlating the real parts or the imaginary parts of the data of the two images in the windows with each other to obtain their correlation coefficient C1; obtaining the correlation coefficient C1 of the real part of the first complex image with the imaginary part of the second complex image in the same way, or computing the correlation coefficient -C2 of the imaginary part of the first complex image with the real part of the second complex image, thereby obtaining C2; finally, obtaining the interferogram from the arc tangent of the ratio of C2 to C1.

Description

生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法  Real-imaginary part correlation method for generating interference-free spot synthetic aperture radar interference phase map
技术领域 Technical field
本发明涉及图像处理技术、遥感测量等领域， 进一步是指用于生成干涉合 成孔径雷达干涉相位图的条紋等值线窗口实虚部相关方法、矩形等窗口实虚部 相关方法。  The present invention relates to the fields of image processing technology, remote sensing measurement and the like, and further refers to a method for correlating real and imaginary parts of a fringe contour window and a rectangle and other windows for generating an interference synthesis aperture radar interference phase map.
背景技术 Background technique
千涉合成孑 L径雷达 ( Interf erometric Synthetic Aperture Radar, 缩写： InSAR或 IFSAR)三维成像技术是新近发展起来的空间观测技术， 它是传统的 合成孔径雷达 （Synthetic Aperture Radar, 缩写： SM)遥感技术与射电天 文干涉技术结合的产物。千涉合成孔径雷达技术通过两幅天线同时观测或一副 天线两次平行观测，接收地面目标的回波信号，经成像处理后得到同一观测区 域两幅具有相干性的合成孔径雷达复数图像 （包括强度信息和相位信息）， 再 经过干涉处理后得到这两幅复图像的干涉相位图像，按照一定的几何关系进行 变换， 进而得到观测目标的三维信息。 目前干涉相位图的生成主要采用共轭相乘的方法。 记 ^(Γ,Χ)和 V2(r，x)分 别为精确配准后的两幅复图像的复值， 贝 iJ: Interferometric Synthetic Aperture Radar (Abbreviation: InSAR or IFSAR) three-dimensional imaging technology is a newly developed space observation technology. It is a traditional Synthetic Aperture Radar (abbreviation: SM) remote sensing technology. Combining with radio astronomy interference technology. Sensing Synthetic Aperture Radar technology uses two antennas to observe at the same time or one antenna to perform two parallel observations to receive echo signals from ground targets. After imaging processing, two coherent synthetic aperture radar images (including (Intensity information and phase information), and then the interference phase images of the two complex images are obtained after interference processing, and the three-dimensional information of the observation target is obtained by transforming according to a certain geometric relationship. At present, the generation of interference phase maps mainly adopts the method of conjugate multiplication. Let ^ (Γ, ×) and V 2 (r, x) be the complex values of the two complex images after exact registration, respectively, iJ:
Vx (?%x) - Αχείφ = A cos φ + iAx sin φ{ (1) V2 (r, x) = A e' - = Α = A2 cos(^ + Δ^) + iA2 sin ( + Δ^) (2) 式中 Α和 分别为两根天线的回波相位， Δ^=ίί2 为回波相位差，即干涉相位。 将两幅复图像对应像素进行复共轭相乘得： ,*(r,x)'-『/2(r,x) = Α,Α,ε'^'^ = Α,Α (3) 式中 *表示复共轭。 由式 （3) 可得-V x (?% X) -Α χ ε ίφ = A cos φ + iA x sin φ { (1) V 2 (r, x) = A e '-= Α = A 2 cos (^ + Δ ^) + iA 2 sin (+ Δ ^) (2) where A and E are the echo phases of the two antennas respectively, and Δ ^ = ίί 2 is the echo phase difference, that is, the interference phase. Multiplying the corresponding pixels of two complex images by complex conjugate multiplication:, * (r, x) '-『/ 2 (r, x) = Α, Α, ε' ^ '^ = Α, Α (3) Middle * indicates complex conjugate. From equation (3), we get-
Λ . r lm( * (/·, x) - (ί% χ)) Ί Λ . R lm (* (/ ·, x)-(ί% χ)) Ί
Δ二 arctan[ ~ ^ ~~ - ~ -、 . 1 (4) e(F, (r,x)-F2(r,x)) Δ 二 arctan [~ ^ ~~-~-,. 1 (4) e (F, (r, x) -F 2 (r, x))
图 1所示为用这种方法得到的干涉相位图。 由于干涉合成孔径雷达系统中 斑点噪声、空间失相关、 时间失相关、数据处理噪声等因素的影响， 使得干涉 相位图的信噪比较低。这种高噪声的存在严重地影响了相位解缠的进行和高精 度数字高程图的获取，已经成为干涉合成? L径雷达数据处理中一个最主要的瓶 颈。 Figure 1 shows the interference phase map obtained by this method. Because of interferometric synthetic aperture radar systems The influence of speckle noise, spatial miscorrelation, temporal miscorrelation, and data processing noise makes the signal to noise of the interference phase map relatively low. The existence of such high noise has seriously affected the phase unwrapping and the acquisition of high-precision digital elevation maps, which has become a major bottleneck in the interference synthesis? L-path radar data processing.
发明内容 Summary of the invention
本发明要解决的技术问题是，针对现有技术存在的缺陷， 在干涉合成孔径 雷达数据处理过程中， 用条紋等值线窗口实虚部相关方法、 矩形等窗口实虚部 相关方法代替传统的共轭相乘方法， 由两巾畐精确配准后的复图像生成完全免除 相干斑噪声的干涉相位图。从而极大地提高相位解缠的可靠性及生成数字高程 图的精度， 解决了干涉合成孔径雷达数据处理中一个最主要的困难。  The technical problem to be solved by the present invention is that, in view of the defects existing in the prior art, in the process of the interference synthetic aperture radar data processing, the method of correlating the real and imaginary parts of the window of the fringe contour and the method of correlating the real and imaginary parts of the window such as a rectangle is used instead of the traditional The conjugate multiplication method generates an interference phase map that is completely free of coherent speckle noise from the complex image after accurate registration of the two frames. This greatly improves the reliability of phase unwrapping and the accuracy of generating digital elevation maps, and solves a major difficulty in the processing of interferometric synthetic aperture radar data.
本发明的技术方案是， 所述生成无干涉斑合成孔径雷达干涉相位图的实虚 部相关方法是：  The technical solution of the present invention is that the method for generating the real and imaginary parts of the interference-free spot synthetic aperture radar interference phase map is:
( 1 )对于两幅复数图像中对应的每一像素点， 以当前点为中心， 取一大小 为讓 (m n为大于 0的数） 的等值线窗口或其它形式窗口， 在窗口内对两 幅复图像数据的实部或虚部进行相关运算， 求出其相关系数 C1 ; (1) For each corresponding pixel point in the two complex images, take the current point as the center, and take a contour window or other window with a size of (mn is a number greater than 0). Correlate the real or imaginary part of the complex image data to find its correlation coefficient C 1;
(2)用上述相同方法对第一幅复图像数据的实部和第二幅复图像数据的虚 部进行相关运算，求出其相关系数 C2，或对第一幅复图像数据的虚部和第二幅 复图像数据的实部进行相关运算， 求出其 目关系数- C2， 进而得到 C2 ; (2) Correlation operation is performed on the real part of the first complex image data and the imaginary part of the second complex image data by the same method as above to find the correlation coefficient C 2 , or the imaginary part of the first complex image data Perform a correlation operation with the real part of the second complex image data to find the number of mesh relations-C 2 , and then obtain C 2;
( 3 ) 对 C2与 的比值求反正切即可获得干涉相位图像。 (3) Obtain the interference phase image by inverse tangent of the ratio of C 2 and.
所述窗口可以是条纹等值线窗口， 也可以是矩形及其它形式窗口。  The window may be a stripe contour window, or may be a rectangular or other window.
以下对本发明做出进一步说明。  The invention is further explained below.
关于本发明原理的公式推导。  Derivation of formulas on the principles of the invention.
使用的相关公式可以是协方差（均值归一化）相关、标准化协方差相关等 多种数学表达形式。 下面以标准化协方差木目关公式为例来加以推导。  The correlation formula used can be a variety of mathematical expressions such as covariance (mean normalized) correlation, standardized covariance correlation, and so on. The following uses the standardized covariance Mumuguan formula as an example to derive.
标准化协方差相关公式如下：
Figure imgf000004_0001
The standardized covariance correlation formula is as follows:
Figure imgf000004_0001
< - >,IM1算符代表在 m x n像素范围内对某一变量求均值。 对于两幅复图像 V r, Χ)和 V2 (r, Χ)分别取其实部数据， 记为：
Figure imgf000005_0001
由于图像精确配准， 可以认为 4=4= ， Α和 为随机分布变量， 由散斑统计 理论， 可以认为在满足一定尺度 mx" 的窗口上： <->, IM1 operator means to average a certain variable in the range of mxn pixels. For the two complex images V r, X ) and V 2 (r, X ), the real data is taken separately and recorded as:
Figure imgf000005_0001
Due to the accurate registration of the image, 4 = 4 = can be considered, and Α sum is a randomly distributed variable. From the speckle statistical theory, it can be considered that the window satisfies a certain scale mx ":
< cos φ, > cos +Α ) > =0 (8 并假设相位变化量 在 MX «个像素窗口区域内保持不变， 由式（6)和式（7) 可得- =0 (10) <cos φ,> cos + Α)> = 0 (8 and assuming that the amount of phase change remains unchanged in the MX «pixel window area, from equations (6) and (7) we can get-= 0 (10)
<( i-</i > )(/2一 =< 1/2 »>„ cos <(I - </ i> ) (/ 2 a = <1/2 »>" cos
<( -< fl f fx >mxl, = < COS2 Φχ<(-<fl f fx> mxl , = <COS 2 Φχ
Figure imgf000005_0002
Figure imgf000005_0003
将（12) (13) (14)代入（5) 中， 得
Figure imgf000005_0002
Figure imgf000005_0003
Substituting (12) (13) (14) into (5), we get
4
Figure imgf000005_0004
4
Figure imgf000005_0004
同理， 分别取 V,(r，X^G V2(r，x)的实部与虚部， 记为: Similarly, take the real and imaginary parts of V, (r, X ^ GV 2 (r, x) respectively, and record them as:
j = Ax cos (16)
Figure imgf000005_0005
j = A x cos (16)
Figure imgf000005_0005
用上面同样的方法进行推导可以得到 -
Figure imgf000006_0001
Using the same method as above, we can get −
Figure imgf000006_0001
由式 （15 ) 与式 （18)， 可以得到:  From (15) and (18), we can get:
^ = sxct {C2 lCx) (19) 对比式（4 )与式 （19 )可知二者得到 涉相位图的相位分布完全一致。 另外， 将第一幅图的虚部与第二幅图的虚部进行相关运算也可以得到 ( ， 将第一幅图的虚部与第二幅图的实部进行相关运算可以得到 -c2，进而得到 C2。 ^ = sxct {C 2 lC x ) (19) By comparing Equation (4) and Equation (19), it can be seen that the phase distributions of the phase diagrams obtained by the two are completely consistent. In addition, the correlation operation between the imaginary part of the first picture and the imaginary part of the second picture can also be obtained (, and the correlation operation between the imaginary part of the first picture and the real part of the second picture can be obtained -c 2 To obtain C 2 .
关于本发明的相关窗口选取。  Relevant window selection related to the present invention.
由上述推导可知， 上述结论是在条紋等值线窗口上才严格成立。 因此， 相 关运算的最优窗口是条纹等值线窗口，用它可以得到质量和精度最好的干涉相 位图。也可以选取矩形及其它形式窗口， 仍可得到比较好的近似结果。条紋等 值线窗口是一种曲线窗口，它是沿干涉相位图的条纹走向取一定长度和宽度所 形成的窗口，可以从共轭相乘方法生成的干涉相位图得到，也可以从矩形等窗 口相关方法生成的干涉相位图得到。可以用不同算法获取条紋等值线窗口，它 不属本发明的内容， 此处不赘述。  It can be known from the above derivation that the above conclusion is strictly established only on the fringe contour window. Therefore, the optimal window for correlation operations is the fringe contour window, which can be used to obtain the interference phase map with the best quality and accuracy. You can also choose rectangular and other forms of windows, and still get better approximate results. The fringe contour window is a curve window. It is a window formed by taking a certain length and width along the fringe direction of the interference phase map. It can be obtained from the interference phase map generated by the conjugate multiplication method, or from a rectangle. The interference phase map generated by the window correlation method is obtained. The stripe contour window can be obtained by using different algorithms, which is not the content of the present invention and will not be repeated here.
由以上可知，本发明为生成无干涉斑合成孔径雷达干涉相位图的实虛部相 关方法，它是一种用基于条纹等值线窗口、矩形等窗口的复图像实虛部相关方 法生成无干涉斑干涉合成孔径雷达干涉 目位图的方法；本发明方法在干涉合成 孔径雷达数据处理过程中，用条紋等值线窗口实虚部相关方法、矩形等窗口实 虚部相关方法代替传统的共轭相乘方法，由两幅精确配准后的复图像生成完全 免除相干斑噪声的干涉相位图。从而极大地提高相位解缠的可靠性及生成数字 高程图的精度， 解决了干涉合成孔径雷达数据处理中一个最主要的困难。 附图说明  It can be known from the above that the present invention is a method for generating a real-imaginary part correlation method of an interference-free spot synthetic aperture radar interference phase map. The method of speckle interference synthetic aperture radar interferometry; in the method of the present invention, in the process of the data processing of the interference synthetic aperture radar, the method of correlating the real and imaginary parts of the fringe contour window and the method of correlating the real and imaginary parts of the rectangle and other windows is used instead of the conventional common The yoke multiplication method generates interference phase maps that are completely free of coherent speckle noise from two accurately registered complex images. This greatly improves the reliability of phase unwrapping and the accuracy of generating digital elevation maps, and solves a major difficulty in data processing of interferometric synthetic aperture radar. BRIEF DESCRIPTION OF THE DRAWINGS
图 1 是精确配准后两幅复图像用共轭相乘方法生成的干涉相位图， 图像 大小为 2048 X 2048像素。  Figure 1 is an interference phase map generated by conjugate multiplication of two complex images after accurate registration. The image size is 2048 X 2048 pixels.
图 2 是精确配准后两幅复图像用矩形窗口相关方法生成的干涉相位图， 所用矩形窗口大小为 19 X 19。  Figure 2 is the interference phase map generated by the rectangular window correlation method for the two complex images after accurate registration, and the size of the rectangular window used is 19 × 19.
图 3 是精确配准后两幅复图像用条紋等值线窗口相关方法生成的干涉相 位图， 所用条纹等值线窗口大小为 41 X 5。 Figure 3 is the interference phase generated by the stripe contour window correlation method for the two complex images after accurate registration. For bitmaps, the stripe contour window size used is 41 X 5.
具体实施方式 detailed description
1. 获取具有相干性的两幅 SAR复图像数据， 并将其精确配准；  1. Obtain two SAR complex image data with coherence and register them accurately;
2.对于两幅复数图像中对应的每一像素点， 以当前点为中心， 取一大小为 的窗口 （可以是条紋等值线窗口， 也可以是矩形及其它形式窗口）， 在窗 口内对两幅复图像数据的实部（或虚部）进行相关运算， 求出其相关系数 C1 ; 2. For each pixel point in the two complex images, take the current point as the center and take a window of size (can be a stripe contour window, or a rectangle and other windows), in the window Perform a correlation operation on the real (or imaginary) part of the two complex image data to find the correlation coefficient C 1;
3.用相同的方法对第一幅复图像数据的实部和第二幅复图像数据的虚部 进行相关运算，求出其相关系数 C2，或是对第一幅复图像数据的虚部和第二幅 复图像数据的实部进行相关运算， 求出其相关系数- C2， 进而得到 C2; 3. Correlate the real part of the first complex image data with the imaginary part of the second complex image data by the same method, and find the correlation coefficient C 2 , or the imaginary part of the first complex image data. Perform a correlation operation with the real part of the second complex image data to find its correlation coefficient-C 2 , and then obtain C 2;
4.对 (：2与 的比值求反正切即可获得干涉相位图像。 4. Find the interference phase image by inverse tangent of the ratio of (: 2 and.

Claims

权 利 要 求 Rights request
1 . 生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法， 其特征是， 该方法为： 1. A method for generating real and imaginary parts of a non-interfering spot synthetic aperture radar interference phase map, characterized in that the method is:
- ( 1 )对于两幅复数图像中对应的每一像素点， 以当前点为中心， 取一大小 为 mxn (m、 n为大于 0的数） 的窗口， 在窗口内对两幅复图像数据的实部或 虛部进行相关运算， 求出其相关系数 C -(1) For each pixel point corresponding to the two complex images, take the current point as the center, take a window of size mxn ( m , n is a number greater than 0), and compare the two complex image data in the window. Perform the correlation operation on the real or imaginary part to find the correlation coefficient C
(2)用上述相同方法对第一幅复图像数据的实部和第二幅复图像数据的虚 部进行相关运算，求出其相关系数 C2，或对第一幅复图像数据的虚部和第二幅 复图像数据的实部进行相关运算， 求出其相关系数 -c2， 进而得到 C2; (2) Correlation operation is performed on the real part of the first complex image data and the imaginary part of the second complex image data by the same method as above to find the correlation coefficient C 2 , or the imaginary part of the first complex image data Perform a correlation operation with the real part of the second complex image data to find its correlation coefficient -c 2 , and then obtain C 2;
(3 ) 对 C2与 的比值求反正切即可获得干涉相位图像。 (3) Obtain the interference phase image by inverse tangent of the ratio of C 2 and.
2.根据权利要求 1所述的生成无干涉斑合成孔径雷达干涉相位图的实虚部 相关方法， 其特征是， 所述窗口是条纹等值线窗口或矩形及其它形式窗口。  The method according to claim 1, wherein the window is a fringe contour window or a rectangular or other window.
3.根据权利要求 2所述的生成无干涉斑合成孔径雷达干涉相位图的实虚部 相关方法， 其特征是， 所述条紋等值线窗口为曲线窗口， 它是沿干涉相位图的 条紋走向取一定长度和宽度所形成的窗口。  The method according to claim 2, wherein the fringe contour window is a curve window, and it is a bar along the interference phase map. The pattern is a window formed by a certain length and width.
4.根据权利要求 3所述的生成无干涉斑合成孔径雷达干涉相位图的实虚部 相关方法，其特征是，所述条紋等值线窗口由传统复图像共轭相乘方法生成的 干涉相位图得到。  The method according to claim 3, wherein the fringe contour window is an interference generated by a conventional complex image conjugate multiplication method The phase map is obtained.
5.根据权利要求 3所述的生成无干涉斑合成孔径雷达干涉相位图的实虚部 相关方法，其特征是，所述条纹等值线窗口由矩形及其它形式窗口实虚部相关 方法生成的干涉相位图得到。  The method according to claim 3, wherein the fringe contour window is generated by a method related to the real and imaginary parts of the rectangle and other forms of windows Interference phase map is obtained.
PCT/CN2005/000346 2004-04-09 2005-03-21 A real-imaginary correlation method for generating a speckless interferogram of sar WO2005098473A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN200410023098.9 2004-04-09
CNB2004100230989A CN1307428C (en) 2004-04-09 2004-04-09 Real-virtual relating method for generating interference phase pattern of synthetic aperture radar without interference spots

Publications (1)

Publication Number Publication Date
WO2005098473A1 true WO2005098473A1 (en) 2005-10-20

Family

ID=35067481

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2005/000346 WO2005098473A1 (en) 2004-04-09 2005-03-21 A real-imaginary correlation method for generating a speckless interferogram of sar

Country Status (2)

Country Link
CN (1) CN1307428C (en)
WO (1) WO2005098473A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108663678A (en) * 2018-01-29 2018-10-16 西北农林科技大学 More baseline InSAR phase unwrapping algorithms based on mixed integer optimization model
CN111461999A (en) * 2020-03-13 2020-07-28 西安工程大学 SAR image speckle suppression method based on super-pixel similarity measurement
CN114910909A (en) * 2022-06-17 2022-08-16 苏州大学 Synthetic aperture radar interferogram phase recovery method and system under data driving
US11933883B2 (en) 2021-09-24 2024-03-19 Aloft Sensing, Inc. System and method for self-contained high-precision navigation

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007071113A1 (en) * 2005-12-22 2007-06-28 Qifeng Yu A three-images isoline correlation interferometric method for generating an interferogram of radar
CN102073992B (en) * 2010-12-09 2012-10-03 国网电力科学研究院 High-resolution SAR satellite image speckle de-noising method
CN102323570B (en) * 2011-05-24 2013-03-13 中国人民解放军国防科学技术大学 Method for estimating magnitude-phase characteristics of radar target echo signal simulator
CN103308914B (en) * 2012-12-06 2014-11-05 北京理工大学 One-station fixed bistatic interference synthetic aperture radar (SAR) processing method
CN103871030B (en) * 2014-02-17 2017-06-30 中国科学院电子学研究所 The filtering method and equipment of a kind of interference image
CN104463874A (en) * 2014-12-12 2015-03-25 中国人民解放军装备学院 Interferometric phase image noise reduction method adopting self-adaptive complex dictionary and sparse coding
CN109708584B (en) * 2018-12-30 2020-07-28 南京航空航天大学 Oil film interference method fringe spacing image identification method
CN112098997B (en) * 2020-09-18 2021-10-15 欧必翼太赫兹科技(北京)有限公司 Three-dimensional holographic imaging security inspection radar image foreign matter detection method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5677693A (en) * 1994-05-26 1997-10-14 Hughes Aircraft Company Multi-pass and multi-channel interferometric synthetic aperture radars
CN1393827A (en) * 2001-06-26 2003-01-29 中国科学院遥感应用研究所 Method for using computer to recreate image containing spot noise
US6532073B2 (en) * 2000-09-13 2003-03-11 Fuji Photo Optical Co., Ltd. Fringe analysis error detection method and fringe analysis error correction method
US6583751B1 (en) * 1999-05-25 2003-06-24 Politecnico Di Milano Process for radar measurements of the movement of city areas and landsliding zones

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2877106B2 (en) * 1996-11-18 1999-03-31 日本電気株式会社 Along track interferometry SAR
CN1148067C (en) * 2001-06-01 2004-04-28 中国科学院计算技术研究所 Data compressing method for complex image of synthetic apertre radar
GB0219645D0 (en) * 2002-08-23 2002-10-02 Secr Defence Radar imaging

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5677693A (en) * 1994-05-26 1997-10-14 Hughes Aircraft Company Multi-pass and multi-channel interferometric synthetic aperture radars
US6583751B1 (en) * 1999-05-25 2003-06-24 Politecnico Di Milano Process for radar measurements of the movement of city areas and landsliding zones
US6532073B2 (en) * 2000-09-13 2003-03-11 Fuji Photo Optical Co., Ltd. Fringe analysis error detection method and fringe analysis error correction method
CN1393827A (en) * 2001-06-26 2003-01-29 中国科学院遥感应用研究所 Method for using computer to recreate image containing spot noise

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108663678A (en) * 2018-01-29 2018-10-16 西北农林科技大学 More baseline InSAR phase unwrapping algorithms based on mixed integer optimization model
CN111461999A (en) * 2020-03-13 2020-07-28 西安工程大学 SAR image speckle suppression method based on super-pixel similarity measurement
CN111461999B (en) * 2020-03-13 2023-02-14 西安工程大学 SAR image speckle suppression method based on super-pixel similarity measurement
US11933883B2 (en) 2021-09-24 2024-03-19 Aloft Sensing, Inc. System and method for self-contained high-precision navigation
CN114910909A (en) * 2022-06-17 2022-08-16 苏州大学 Synthetic aperture radar interferogram phase recovery method and system under data driving

Also Published As

Publication number Publication date
CN1307428C (en) 2007-03-28
CN1680826A (en) 2005-10-12

Similar Documents

Publication Publication Date Title
WO2005098473A1 (en) A real-imaginary correlation method for generating a speckless interferogram of sar
US9417323B2 (en) SAR point cloud generation system
EP3229038B1 (en) Wavelet domain insar interferometric phase filtering method in combination with local frequency estimation
Plyer et al. A new coregistration algorithm for recent applications on urban SAR images
Bhattacharya et al. Stereo-particle image velocimetry uncertainty quantification
EP3022582A1 (en) Method for filtering of interferometric data acquired by synthetic aperture radar (sar)
CN100545676C (en) Method for interfering synthetic aperture radar interferometric phase estimation based on related weighing
CN104730505A (en) Multi-channel SAR ground target detecting and positioning method and device
CN103439708A (en) Polarized InSAR interferogram estimation method based on generalized scattering vector
Duman et al. Moving target artifacts in bistatic synthetic aperture radar images
CN109270527B (en) Circular SAR sub-aperture image sequence combined correlation DEM extraction method
Serafino SAR image coregistration based on isolated point scatterers
WO2007071113A1 (en) A three-images isoline correlation interferometric method for generating an interferogram of radar
Feng et al. A phase calibration method based on phase gradient autofocus for airborne holographic SAR imaging
CN112986949A (en) SAR high-precision time sequence deformation monitoring method and device for diagonal reflector
Torgrimsson et al. SAR processing without a motion measurement system
CN115685200A (en) High-precision large-front-squint SAR imaging motion compensation and geometric correction method
CN103792534B (en) SAR two-dimension autofocus method based on prior phase structure knowledge
CN103630898A (en) Method for estimating multi-baseline interferometry SAR phase bias
CN108008382B (en) A kind of method of more base spaceborne interferometric SAR systematic survey mountain terrains
CN110554377A (en) Single-channel SAR two-dimensional flow field inversion method and system based on Doppler center offset
CN112505647A (en) Moving target azimuth speed estimation method based on sequential sub-image sequence
CN104833971A (en) Bistatic/multistatic radar image PS (permanent scatterer) point associating method based on sliding scattering center
Peterson et al. Registration of multi-frequency SAR imagery using phase correlation methods
KR100425283B1 (en) Topography Measurement Method using InSAR

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SM SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

DPEN Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed from 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
NENP Non-entry into the national phase

Ref country code: DE

WWW Wipo information: withdrawn in national office

Country of ref document: DE

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 69(1) EPC

122 Ep: pct application non-entry in european phase