CN112986994B - SAR chromatographic reference network rapid generation method - Google Patents

SAR chromatographic reference network rapid generation method Download PDF

Info

Publication number
CN112986994B
CN112986994B CN202110175560.0A CN202110175560A CN112986994B CN 112986994 B CN112986994 B CN 112986994B CN 202110175560 A CN202110175560 A CN 202110175560A CN 112986994 B CN112986994 B CN 112986994B
Authority
CN
China
Prior art keywords
sub
block
sar
reference network
generating
Prior art date
Legal status (The legal status 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 status listed.)
Active
Application number
CN202110175560.0A
Other languages
Chinese (zh)
Other versions
CN112986994A (en
Inventor
余安喜
朱小祥
董臻
张永胜
吴曼青
何峰
何志华
李德鑫
金光虎
孙造宇
王先涛
王莘景
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 filed Critical National University of Defense Technology
Priority to CN202110175560.0A priority Critical patent/CN112986994B/en
Publication of CN112986994A publication Critical patent/CN112986994A/en
Application granted granted Critical
Publication of CN112986994B publication Critical patent/CN112986994B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

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)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

The invention provides a method for quickly generating a reference network based on a segmented search block networking. The technical scheme is that the SAR image is uniformly divided into a plurality of sub-blocks; generating a sub-block communication network in each sub-block by adopting a searching mode based on segmented ant colony; establishing a connecting edge between adjacent sub-blocks; and generating the SAR chromatographic reference network by using the information. The method has high calculation efficiency, and the generated SAR chromatographic reference network has strong full-scene coverage capability.

Description

SAR chromatographic reference network rapid generation method
Technical Field
The invention belongs to the technical field of remote sensing imaging, and relates to a reference model which is needed to be utilized for eliminating atmospheric influence when a plurality of satellite-borne SAR (Synthetic Aperture Radar ) images are utilized for high-dimensional imaging.
Background
The space-borne SAR three-dimensional imaging is also commonly called space-borne SAR tomography, namely, by utilizing multi-angle coherent observation on the same target, a synthetic aperture is formed in the height direction, the space-borne SAR three-dimensional imaging has high-directional resolution, and true three-dimensional imaging on the observed target can be realized. The main data source of SAR chromatography comes from a satellite-borne multi-navigation system, and APS (Atmospheric Phase Screen, atmospheric phase error) of each navigation SAR image obtained by the system changes randomly along with the observation time, so that the SAR chromatography can be seriously influenced. In order to eliminate the atmospheric influence, the SAR chromatographic processing mainly uses a reference network technology to realize the full-scene atmospheric phase correction, see document one: robust Detection of Single and Double Persistent Scatterers in Urban Built Environments, p.f.ma et al. The basic idea of the technology is as follows: based on the spatial ramp characteristic of the APS, the APS approximation of the adjacent pixels in the SAR image can be considered, and the APS error can be largely counteracted through the phase difference of the adjacent pixels. By utilizing the basic idea, the reference network technology is to connect single PS (Persistent Scatterer, permanent scattering points) covering the whole scene into a network (i.e. a reference network) covering the whole scene, and the rest pixels take the single PS point nearest to the reference network as a reference to perform phase difference so as to realize the APS correction and PS point detection of the whole scene. For reference network technology, the quality of the reference network will directly determine whether full scene APS valid correction can be achieved. Theoretically, it is desirable that the generated reference network can uniformly cover the full scene.
The reference network provided in reference one is the maximum connected network after the Delaunay triangulation network formed based on Delaunay triangulation is optimized. However, due to low network redundancy of the Delaunay triangle network, the problems of PS isolated points and isolated islands are easy to form after network optimization, so that the reference network does not have full scene coverage capability, and full scene APS correction cannot be realized based on the reference network. In addition, the maximum connected network search process is extremely time consuming, which can result in inefficient reference network generation.
Disclosure of Invention
The invention provides a rapid generation method of a reference network based on segmented search and block networking, the generated reference network has full scene coverage capability, and the generation efficiency of the reference network is very high.
The technical scheme of the invention is as follows: a rapid generation method of a reference network based on segmented search and block networking is characterized in that SAR images are uniformly divided into a plurality of sub-blocks; generating a sub-block communication network in each sub-block; establishing a connecting edge between adjacent sub-blocks; and generating the SAR chromatographic reference network by using the information.
Further, a sub-block communication network is generated in each sub-block in a mode based on segmented ant colony search.
The invention has the technical effects that: the rapid generation method of the reference network does not need to carry out the maximum connected network search, reduces the calculation time waste caused by directly carrying out single PS search in the whole search space through the segmentation search, can rapidly connect all the sub-blocks to form a large connected network through a block networking mode, and can greatly accelerate the generation efficiency of the large scene reference network in two modes, namely, the calculation efficiency of the method is very high. In addition, an ant colony network uniformly covering each sub-block is generated in each sub-block in an ant colony searching mode, and then the sub-block communication networks can be communicated into a maximum reference network uniformly covering the whole scene by establishing high-quality edges among the sub-block communication networks, namely, the SAR chromatographic reference network obtained based on the method has the full scene coverage capability.
Drawings
FIG. 1 is a schematic diagram of a SAR chromatographic reference network rapid generation method based on segmented search and block networking provided by the invention;
FIG. 2 is a flowchart of a method for rapidly generating SAR chromatographic reference network based on segmented search and block networking provided by the invention;
fig. 3 to 7 show results of actual measurement data processing according to the embodiment of the present invention.
Detailed Description
The SAR chromatography reference network generation method based on ant colony search provided by the invention is described in detail below with reference to the accompanying drawings.
The whole process is divided into three steps. Step1, region segmentation; step2, quickly generating a sub-block communication network; step3: the sub-block communication net is communicated. The specific description is as follows:
step1: region partitioning
And uniformly dividing the whole scene, namely the whole SAR image, into a plurality of sub-blocks. Suppose the SAR image size is: x is X a ×X r (Unit: pixel), the sub-block size set in the partitioning process is x a ×x r . The size of the sub-block is determined according to the actual situation. Through laboratory experiments, it was found that the selection is generallyPreferably, wherein: />Representing a rounding operator.
Step2: fast generation of sub-block communication network
A sub-block communication network is generated within each sub-block. Typically, the generation of the sub-block communication network may be formed using existing Delaunay triangulation-based techniques. The invention recommends to use a segmented ant colony search mode to generate a subblock communication network, and the process is as follows:
step (1): threshold distance D thres Divided into L subsections, and search space D corresponding to the ith subsection (i is more than or equal to 1 and less than or equal to L) i The definition is as follows: with a distance d from the picture element as the center i In the region of all picture elements in whichDistance threshold D thres The value of (2) is determined according to the actual situation, and is usually in the order of hundreds of meters. Experiments prove that the optimal selection range of the segmentation number is more than or equal to 5 and less than or equal to 20.
Step (2): let i=1, choose any pixel in the sub-block as the initial ant, in the search space corresponding to the 1 st sub-segment, adopt patent number: 202011621544.1 the method provides for the generation and screening of satisfactory P around the starting ant i A high-quality single PS side is formed, and the end point of the screened single PS side is used as a new starting ant;
step (3): let i=i+1, in the search space corresponding to the ith (1. Ltoreq.i.ltoreq.L) sub-segment, patent number: 202011621544.1 the method comprises the steps of generating and screening out high-quality single PS edges meeting the requirements around each initial ant, and taking the end point of the screened single PS edge as a new initial ant;
step (4): if the total number of the high-quality single PS sides screened in the first i sections of search space meets the cycle ending condition: i.e. p=p 1 +…+P i If not less than F (F represents a network redundancy threshold), ending the cycle; otherwiseAnd (3) returning to the step (3) to continue to search for the high-quality single PS side in the search space corresponding to the next sub-segment, and taking the end point of the screened single PS side as a new starting ant until the cycle end condition is met or all the sub-segments (i.e. i=L) are searched.
By searching in a sectional expansion range in the constraint space around each ant, the network construction time required by the sub-block connected network can be greatly reduced on the premise of not influencing the quality of the reference network edge construction. In the processing process of the step, all the sub-blocks are mutually independent, and the establishing process of the sub-block communication network is not affected, so that the processing can be performed in a parallel mode.
Step3: sub-block communication net communication
And establishing a connecting edge between the sub-block communication networks to form an SAR chromatographic reference network which uniformly covers the whole scene. I.e. between adjacent sub-block communication networks, using reference two: the method provided by New Approaches For Robust and Efficient Detection of Persistent Scatterers in SAR Tomography, X.X.Zhu et al establishes the connecting edges.
Fig. 3 to 7 show the results of actual measurement data processing.
FIG. 3 is a terraSAR-X average SAR intensity map. The horizontal direction represents the azimuth direction, the vertical direction represents the distance direction, and the resolutions of the azimuth direction and the distance direction are respectively 0.24m and 0.6m. The image is obtained by averaging the amplitude images of the SAR images of the 26-scene terra SAR-X satellite staring spotlight mode.
Fig. 4 is a reference mesh generated using the prior art method provided in reference one, which can only cover a portion of the viewing area. As can be seen from the figure, neither area 1 nor area 2 is effectively covered, and the line-shaped road object in area 3 is not covered.
Fig. 5 is a reference network established using the proposed method of the present invention. Distance-azimuth size X of SAR image a ×X r Size x of subblock 1370×3826 a ×x r For 688×1913, the whole SAR image is divided into 2×2 blocks, distance threshold D thres The number of segments is set to l=10, =300 meters. Compared with the reference network generated by the prior method, the reference network generated based on the method provided by the invention can cover the whole scene and has the whole sceneScene coverage capability.
Fig. 6 is a reference grid PS-point elevation map generated based on the prior art method provided in reference one. The horizontal direction represents the azimuth direction, the vertical direction represents the distance direction, the gray scale represents the height, and the hexagon in the figure represents the gray scale value of the zero elevation reference point. Since the reference network generated by the existing method cannot cover the current road targets in the area 1, the area 2 and the area 3, the PS point heights of the areas cannot be effectively inverted.
Fig. 7 is a reference network PS point elevation map generated based on the proposed method of the present invention. The horizontal direction represents the azimuth direction, the vertical direction represents the distance direction, the gray scale represents the height, and the hexagon in the figure represents the gray scale value of the zero elevation reference point. The reference network generated based on the method provided by the invention can cover the full scene, and the PS point heights of the full scene can be effectively inverted, so that the effectiveness of the method provided by the invention is proved.
The calculation time of fig. 4 is 832.8 seconds by using the existing reference network generation method provided by the reference document one, the calculation time of fig. 5 is 131.6 seconds by using the reference network generation method provided by the invention, the calculation efficiency of the reference network generation method provided by the invention is improved by about 6.3 times, and the high efficiency of the method provided by the invention is proved.

Claims (5)

1. A rapid generation method of a reference network based on segmented search and block networking is characterized in that SAR images are uniformly divided into a plurality of sub-blocks; generating a sub-block communication network in each sub-block; establishing a connecting edge between adjacent sub-blocks; generating a chromatographic reference network of the SAR image by utilizing the information; wherein SAR refers to synthetic aperture radar;
the process of generating the sub-block communication network in each sub-block is as follows:
step (1): threshold distance D thres Divided into L subsections, search space D corresponding to the ith subsection i The definition is as follows: with a distance d from the picture element as the center i In the region of all picture elements in which
Step (2): let i=1, select any one of the pixels in the sub-block as the initial ant, generate and screen out the P meeting the requirement around the initial ant in the search space corresponding to the 1 st sub-segment i The end point of the single PS side is used as a new starting ant, and PS refers to a permanent scattering point;
step (3): generating and screening out high-quality single PS edges meeting the requirements around each initial ant in a search space corresponding to the ith subsection, and taking the end point of the screened single PS edge as a new initial ant;
step (4): if the total number of the high-quality single PS sides screened in the first i sections of search space meets the cycle ending condition: i.e. p=p 1 +…+P i F represents a network redundancy threshold when the value is not less than F, and the cycle is ended; otherwise, returning to the step (3), continuing to search the high-quality single PS side in the search space corresponding to the next sub-section, and taking the end point of the screened single PS side as a new starting ant until the cycle ending condition is met or all the sub-sections are searched.
2. The method for quickly generating the reference network based on the segmented search and block networking according to claim 1, wherein the following process is adopted to uniformly divide the image into a plurality of sub-blocks:
suppose the SAR image size is: x is X a ×X r ,X a And X r The sub-block size is x a ×x r Representing the number of pixels in two directions of the sub-block,wherein: />Representing a rounding operator.
3. The method for quickly generating the reference network based on the segmented search and block networking according to claim 2, wherein the selection range of the number of segments of the sub-segments is 5-20.
4. A method for high-dimensional imaging using a plurality of spaceborne SAR images, wherein a chromatographic reference network of SAR images is established using the method of any of claims 1 to 3.
5. A system for high-dimensional imaging using a plurality of satellite-borne SAR images, wherein a chromatographic reference network of SAR images is established using the method of any of claims 1 to 3.
CN202110175560.0A 2021-02-06 2021-02-06 SAR chromatographic reference network rapid generation method Active CN112986994B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110175560.0A CN112986994B (en) 2021-02-06 2021-02-06 SAR chromatographic reference network rapid generation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110175560.0A CN112986994B (en) 2021-02-06 2021-02-06 SAR chromatographic reference network rapid generation method

Publications (2)

Publication Number Publication Date
CN112986994A CN112986994A (en) 2021-06-18
CN112986994B true CN112986994B (en) 2023-11-17

Family

ID=76392490

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110175560.0A Active CN112986994B (en) 2021-02-06 2021-02-06 SAR chromatographic reference network rapid generation method

Country Status (1)

Country Link
CN (1) CN112986994B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014126638A1 (en) * 2013-02-15 2014-08-21 Raytheon Company Sar image formation
CN105954747A (en) * 2016-06-20 2016-09-21 中国电力工程顾问集团中南电力设计院有限公司 Tower foundation stability analyzing method based on three-dimensional deformation monitoring of unfavorable geologic body of power grid
DE102016107065B3 (en) * 2016-04-15 2017-07-27 Deutsches Zentrum für Luft- und Raumfahrt e.V. Method and system for determining movements of long-term stable radar spreaders on a surface of a planet by means of SAR interferometry
CN108387899A (en) * 2018-04-17 2018-08-10 南京师范大学 Ground control point automatically selecting method in synthetic aperture radar interferometry
CN109212528A (en) * 2018-10-30 2019-01-15 中国科学院遥感与数字地球研究所 A kind of multidimensional chromatography historic site SAR ruins deformation monitoring method
CN110221298A (en) * 2019-06-20 2019-09-10 中国人民解放军国防科技大学 Self-focusing method for low-frequency-band spaceborne SAR ionosphere scintillation effect
CN110412574A (en) * 2019-09-05 2019-11-05 河海大学 A kind of distributed object InSAR timing sequence process method and apparatus of temporal and spatial coherence enhancing
CN111812647A (en) * 2020-07-11 2020-10-23 桂林电子科技大学 Phase unwrapping method for interferometric synthetic aperture radar

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6670909B2 (en) * 2001-01-16 2003-12-30 Time Domain Corporation Ultra-wideband smart sensor interface network and method
CN104765026A (en) * 2015-04-29 2015-07-08 天津市测绘院 Method for extracting ground attribute data in interferometry synthetic aperture radar data
US10359514B2 (en) * 2015-12-10 2019-07-23 The Climate Corporation Generating estimates of uncertainty for radar based precipitation estimates
CN108627833B (en) * 2018-05-15 2021-08-24 电子科技大学 Atmospheric phase compensation method based on GB-InSAR
US11231715B2 (en) * 2018-05-22 2022-01-25 King Fahd University Of Petroleum And Minerals Method and system for controlling a vehicle
CN110261394B (en) * 2019-06-24 2022-09-16 内蒙古工业大学 Online real-time diagnosis system and method for damage of fan blade

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014126638A1 (en) * 2013-02-15 2014-08-21 Raytheon Company Sar image formation
DE102016107065B3 (en) * 2016-04-15 2017-07-27 Deutsches Zentrum für Luft- und Raumfahrt e.V. Method and system for determining movements of long-term stable radar spreaders on a surface of a planet by means of SAR interferometry
CN105954747A (en) * 2016-06-20 2016-09-21 中国电力工程顾问集团中南电力设计院有限公司 Tower foundation stability analyzing method based on three-dimensional deformation monitoring of unfavorable geologic body of power grid
CN108387899A (en) * 2018-04-17 2018-08-10 南京师范大学 Ground control point automatically selecting method in synthetic aperture radar interferometry
CN109212528A (en) * 2018-10-30 2019-01-15 中国科学院遥感与数字地球研究所 A kind of multidimensional chromatography historic site SAR ruins deformation monitoring method
CN110221298A (en) * 2019-06-20 2019-09-10 中国人民解放军国防科技大学 Self-focusing method for low-frequency-band spaceborne SAR ionosphere scintillation effect
CN110412574A (en) * 2019-09-05 2019-11-05 河海大学 A kind of distributed object InSAR timing sequence process method and apparatus of temporal and spatial coherence enhancing
CN111812647A (en) * 2020-07-11 2020-10-23 桂林电子科技大学 Phase unwrapping method for interferometric synthetic aperture radar

Non-Patent Citations (12)

* Cited by examiner, † Cited by third party
Title
A Tutorial on Speckle Reduction in Synthetic Aperture Radar Images;Fabrizio Argenti 等;《IEEE Geoscience and Remote Sensing Magazine》;第6-35页 *
Advanced differential interferometry synthetic aperture radar techniques for deformation monitoring: a review on sensors and recent research development;O.Idrees Mohammed 等;《Geocarto International》;第536-553页 *
Phase unwrapping for SAR interferometry based on an ant colony optimization algorithm;Z.-Q. Wei 等;《International Journal of Remote Sensing》;第711-725页 *
Super-Resolution Power and Robustness of Compressive Sensing for Spectral Estimation With Application to Spaceborne Tomographic SAR;Xiao Xiang Zhu 等;《IEEE Transactions on Geoscience and Remote Sensing》;第247-258页 *
一种新的SAR图像分块匹配算法;张东兴 等;《电讯技术》;第726-729页 *
一种用于SAR图像统计模型参数估计的蚁群算法;曹兰英 等;《数据采集与处理》;第78-81页 *
基于自适应蚁群的FCM聚类优化算法研究;谈玲珑 等;《测控技术》;第46-50页 *
基于虚拟特显点的高精度InSAR性能评估方法;余安喜 等;《地球物理学报》;第865-870页 *
多方位角多基线星载SAR三维成像方法研究;匡辉 等;《雷达学报》;第685-695页 *
时序合成孔径雷达干涉测量的多级化解算;胡凤鸣 等;《同济大学学报(自然科学版)》;第1817-1824页 *
永久散射体差分干涉测量技术中SAR影像精配准的一种新方法;陶秋香 等;《测绘学报》;第69-73页 *
高低轨协同多基SAR关键技术分析;于海锋 等;《第三届高分辨率对地观测学术年会(天基对地观测技术分会)优秀论文集》;第1-7页 *

Also Published As

Publication number Publication date
CN112986994A (en) 2021-06-18

Similar Documents

Publication Publication Date Title
CN112946628A (en) Road running state detection method and system based on radar and video fusion
CN110244302B (en) Three-dimensional transformation method for image pixel coordinates of foundation synthetic aperture radar
JPH09189762A (en) Method for observing variation of ground surface using radar and synthetic aperture radar and transponder therefor
CN110568433A (en) High-altitude parabolic detection method based on millimeter wave radar
CN114067223B (en) Night city excessive illumination monitoring method based on noctilucent remote sensing data
CN112558076B (en) Volume scanning mode calculation method based on networking weather radar coverage area and application
LU504084B1 (en) A non-rigid registration method of vehicle-borne LiDAR point clouds by fusing multiple road features
CN105180852B (en) GB SAR deformation monitoring methods based on triple steppings
CN115267815A (en) Road side laser radar group optimization layout method based on point cloud modeling
JP2001328600A (en) Landing point searching device, flying object using therewith and landing point evaluating device
CN112986994B (en) SAR chromatographic reference network rapid generation method
CN114463521B (en) Building target point cloud rapid generation method for air-ground image data fusion
CN112857312B (en) Fusion method for measuring ground settlement by different time sequence differential interference according to settlement rate
Rieke-Zapp et al. Digital photogrammetry for measuring soil surface roughness
CN106910178B (en) Multi-angle SAR image fusion method based on tone statistical characteristic classification
KR102185307B1 (en) Method and system for high resolving object response of sar images
CN113008202B (en) Ground settlement monitoring method integrating different synthetic aperture radar interferometry
Einfalt et al. Comparison of radar and raingauge measurements during heavy rainfall
CN110969576A (en) Highway pavement image splicing method based on roadside PTZ camera
CN115908294A (en) High-voltage isolating switch opening and closing state distinguishing method based on 3D laser radar
FI129476B (en) Antenna radiation pattern extraction using sparse field measurements
CN115760933A (en) SAR image registration method, device, equipment and medium
CN114355347A (en) Image characteristic enhancement method and system based on multi-azimuth SAR observation
CN112799065B (en) SAR (synthetic Aperture Radar) chromatographic reference network generation method based on ant colony search
Brown et al. UK planning model for digital terrestrial television coverage

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant