CN112034463A - High-precision settlement monitoring method and device for SAR data - Google Patents
High-precision settlement monitoring method and device for SAR data Download PDFInfo
- Publication number
- CN112034463A CN112034463A CN202011226759.3A CN202011226759A CN112034463A CN 112034463 A CN112034463 A CN 112034463A CN 202011226759 A CN202011226759 A CN 202011226759A CN 112034463 A CN112034463 A CN 112034463A
- Authority
- CN
- China
- Prior art keywords
- sar data
- target
- processed
- processing
- phase
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000012544 monitoring process Methods 0.000 title claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 76
- 238000012937 correction Methods 0.000 claims abstract description 25
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 24
- 239000006185 dispersion Substances 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims description 24
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000001228 spectrum Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 7
- 238000012952 Resampling Methods 0.000 claims description 6
- 238000012806 monitoring device Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 230000001427 coherent effect Effects 0.000 claims description 3
- 238000004891 communication Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
-
- 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/9094—Theoretical aspects
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/418—Theoretical aspects
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a high-precision settlement monitoring method and device of SAR data, which relate to the technical field of data processing and comprise the following steps: acquiring SAR data to be processed input by a user, wherein the SAR data to be processed is SAR data in a TOPS mode; carrying out high-precision registration processing on SAR data to be processed to obtain a multi-scene image; determining preselected points of scatterers of the multi-scene images by using an amplitude dispersion threshold algorithm; carrying out phase correction processing on the preselected points to obtain target preset points, and determining a target scatterer by using the target preset points; the phase unwrapping processing and the atmosphere and orbit correction processing are sequentially carried out on the target scatterer, so that the linear deformation information of the region corresponding to the SAR data to be processed is obtained, and the technical problem that the existing settlement monitoring method for the SAR data is low in monitoring accuracy of the linear deformation information of the monitored region is solved.
Description
Technical Field
The invention relates to the technical field of data processing, in particular to a high-precision settlement monitoring method and device for SAR data.
Background
The prior art discloses a double-sample-based estimation method for the direction offset of enhanced spectral diversity (application number: CN 201910597877.6). firstly, a sentinel No. 1 TOPS wide image stack is geometrically registered on the basis of precise orbit data and a digital elevation model to obtain a single-view complex image stack, and a prior statistical characteristic of gamma distribution is utilized to extract a homogeneous pixel from a beam overlapping area in a TOPS imaging mode; carrying out accurate coherence estimation on the overlapped wave beams by using the homogeneous pixel extracted by the double samples; obtaining an interference pattern by utilizing an overlapping area of adjacent beams of the main image and the auxiliary image, and continuously carrying out complex conjugation on the two generated interference patterns to obtain the azimuth offset of the enhanced spectrum diversity; and performing coherence weighting on the enhanced spectrum diversity by using the calculated accurate coherence to obtain accurate azimuth offset of the enhanced spectrum diversity, and then extracting azimuth offset of geometric registration residue.
The TOPS mode data registration method is mainly used for data registration, and is considered that certain error interference still exists in the differential phase in the interference difference of the registered images, and the requirement on ground monitoring precision is high in production practice, so that the differential interference phase obtained by registering SAR data is poor in precision, the automation degree is low, and the production requirement cannot be met.
No effective solution has been proposed to the above problems.
Disclosure of Invention
In view of this, the present invention aims to provide a high-precision settlement monitoring method and device for SAR data, so as to alleviate the technical problem that the existing settlement monitoring method for SAR data has low monitoring precision on linear deformation information of a monitored area.
In a first aspect, an embodiment of the present invention provides a high-precision settlement monitoring method for SAR data, including: acquiring SAR data to be processed input by a user, wherein the SAR data to be processed is SAR data in a TOPS mode; carrying out high-precision registration processing on the SAR data to be processed to obtain a multi-scene image; determining preselected points of scatterers of the multi-scene images by using an amplitude dispersion threshold algorithm; carrying out phase correction processing on the preselected points to obtain target preset points, and determining a target scatterer by using the target preset points; and sequentially carrying out phase unwrapping processing and atmospheric and orbit correction processing on the target scatterer to obtain linear deformation information of an area corresponding to the SAR data to be processed.
Further, the SAR data to be processed comprises: the method comprises the following steps of carrying out high-precision registration processing on SAR data to be processed by main image data and auxiliary image data to obtain a multi-scene image, wherein the high-precision registration processing comprises the following steps: performing geometric transformation on the SAR data to be processed by using a digital elevation model and satellite orbit data to obtain an initial registration offset coefficient of the SAR data to be processed, wherein the initial registration offset coefficient is used for representing the coordinate difference of the pixels of the same-name points between the main image and the auxiliary image; resampling the auxiliary image by using the initial registration offset coefficient to obtain a target auxiliary image; and registering the target auxiliary image by using an enhanced spectrum grading registration algorithm to obtain the multi-scene image.
Further, determining preselected points of scatterers of the multi-view image by using an amplitude dispersion threshold algorithm, comprising: extracting amplitude information of the multi-scene image; calculating the average amplitude value and the amplitude standard deviation of each pixel in the multi-scene image by using the amplitude information; calculating the ratio of the average amplitude value to the amplitude standard deviation of each pixel; and determining the pixel of which the ratio is smaller than a first preset threshold value as the preselected point.
Further, performing phase correction processing on the preselected point to obtain a target preselected point, including: calculating the spatial correlation component of the preset point by using a spatial adaptive filtering algorithm; calculating a coherent value and a spatial uncorrelated error of the preset point by using the spatial correlation component; and determining a target preset point by using the coherence value of the preset point and a second preset threshold, wherein the second preset threshold is the proportion of the pseudo scatterers in the acceptable scatterer set.
Further, sequentially performing phase unwrapping processing and atmospheric and orbit correction processing on the target scatterer to obtain linear deformation information of a region corresponding to the SAR data to be processed, including: performing the phase unwrapping processing on the target scatterer to obtain an unwrapped phase; calculating a target matrix by using a least square algorithm and the unwrapping phase, wherein the first column of the target matrix is an atmospheric and orbital error phase of the main image, and the third column is a linear deformation rate; removing the atmospheric and orbit error phase in the unwrapping phase to obtain a target unwrapping phase; and carrying out low-pass filtering processing on the target unwrapping phase to obtain linear deformation information of the region corresponding to the SAR data to be processed.
Further, the low-pass filtering processing is performed on the target unwrapping phase to obtain linear deformation information of a region corresponding to the SAR data to be processed, and the processing includes: performing low-pass filtering processing in a time domain on the target unwrapping phase to obtain a first filter, wherein the first filter includes: high-frequency filtering and low-frequency filtering; determining the low-frequency filtering as a deformation phase of a region corresponding to the SAR data to be processed; and carrying out low-pass filtering processing and denoising processing in a spatial domain on the high-frequency filtering to obtain the atmospheric and orbital errors of the auxiliary image.
In a second aspect, an embodiment of the present invention further provides a high-precision settlement monitoring device for SAR data, where the device includes: the SAR processing device comprises an acquisition unit, a registration unit, a first determination unit, a second determination unit and a processing unit, wherein the acquisition unit is used for acquiring SAR data to be processed input by a user, and the SAR data to be processed is SAR data in a TOPS mode; the registration unit is used for carrying out high-precision registration processing on the SAR data to be processed to obtain a multi-scene image; the first determining unit is used for determining preselected points of the scatterers of the multi-scene images by using an amplitude dispersion threshold algorithm; the second determining unit is used for performing phase correction processing on the preselected point to obtain a target preset point and determining a target scatterer by using the target preset point; and the processing unit is used for sequentially carrying out phase unwrapping processing and atmospheric and orbit correction processing on the target scatterer to obtain linear deformation information of an area corresponding to the SAR data to be processed.
Further, the SAR data to be processed comprises: the registration unit is configured to: performing geometric transformation on the SAR data to be processed by using a digital elevation model and satellite orbit data to obtain an initial registration offset coefficient of the SAR data to be processed, wherein the initial registration offset coefficient is used for representing the coordinate difference of the pixels of the same-name points between the main image and the auxiliary image; resampling the auxiliary image by using the initial registration offset coefficient to obtain a target auxiliary image; and registering the target auxiliary image by using an enhanced spectrum grading registration algorithm to obtain the multi-scene image.
In a third aspect, an embodiment of the present invention further provides an electronic device, including a memory and a processor, where the memory is used to store a program that supports the processor to execute the method in the first aspect, and the processor is configured to execute the program stored in the memory.
In a fourth aspect, the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to perform the steps of the method in the first aspect.
In the embodiment of the invention, SAR data to be processed input by a user is acquired; carrying out high-precision registration processing on SAR data to be processed to obtain a multi-scene image; determining preselected points of scatterers of the multi-scene images by using an amplitude dispersion threshold algorithm; carrying out phase correction processing on the preselected points to obtain target preset points, and determining a target scatterer by using the target preset points; the phase unwrapping processing and the atmosphere and orbit correction processing are sequentially carried out on the target scatterer, so that the linear deformation information of the region corresponding to the SAR data to be processed is obtained, the purpose of monitoring the linear deformation information of the monitored region is achieved, the technical problem that the monitoring precision of the existing settlement monitoring method for the SAR data on the linear deformation information of the monitored region is low is solved, and the technical effect of carrying out high-precision monitoring on the linear deformation information of the monitored region is achieved.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a flowchart of a high-precision settlement monitoring method for SAR data according to an embodiment of the present invention;
fig. 2 is a detailed flowchart of a high-precision settlement monitoring method for SAR data according to an embodiment of the present invention;
fig. 3 is a schematic diagram of a high-precision settlement monitoring device for SAR data according to an embodiment of the present invention;
fig. 4 is a schematic diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The first embodiment is as follows:
according to an embodiment of the present invention, there is provided an embodiment of a method for high-precision settlement monitoring of SAR data, it is noted that the steps illustrated in the flowchart of the figure may be performed in a computer system such as a set of computer executable instructions, and that although a logical order is illustrated in the flowchart, in some cases the steps illustrated or described may be performed in an order different than here.
Fig. 1 is a flowchart of a method for monitoring high-precision settlement of SAR data according to an embodiment of the present invention, as shown in fig. 1, the method includes the following steps:
step S102, obtaining SAR data to be processed input by a user, wherein the SAR data to be processed is SAR data in a TOPS mode;
step S104, carrying out high-precision registration processing on the SAR data to be processed to obtain a multi-scene image;
step S106, determining preselected points of scatterers of the multi-scene images by using an amplitude dispersion threshold algorithm;
step S108, carrying out phase correction processing on the preselected points to obtain target preset points, and determining a target scatterer by using the target preset points;
and step S110, sequentially carrying out phase unwrapping processing and atmospheric and orbit correction processing on the target scatterer to obtain linear deformation information of an area corresponding to the SAR data to be processed.
In the embodiment of the invention, SAR data to be processed input by a user is acquired; carrying out high-precision registration processing on SAR data to be processed to obtain a multi-scene image; determining preselected points of scatterers of the multi-scene images by using an amplitude dispersion threshold algorithm; carrying out phase correction processing on the preselected points to obtain target preset points, and determining a target scatterer by using the target preset points; the phase unwrapping processing and the atmosphere and orbit correction processing are sequentially carried out on the target scatterer, so that the linear deformation information of the region corresponding to the SAR data to be processed is obtained, the purpose of monitoring the linear deformation information of the monitored region is achieved, the technical problem that the monitoring precision of the existing settlement monitoring method for the SAR data on the linear deformation information of the monitored region is low is solved, and the technical effect of carrying out high-precision monitoring on the linear deformation information of the monitored region is achieved.
The above steps S102 to S108 are described in detail with reference to fig. 1 and fig. 2, and fig. 2 is a detailed flowchart of S102 to S108.
In an embodiment of the present invention, the to-be-processed SAR data includes: the main image data and the auxiliary image data, step S104 includes the following steps:
step S11, performing geometric transformation on the SAR data to be processed by using a digital elevation model and satellite orbit data to obtain an initial registration offset coefficient of the SAR data to be processed, wherein the initial registration offset coefficient is used for representing the coordinate difference of the pixels of the same-name points between the main image and the auxiliary image;
step S12, resampling the auxiliary image by using the initial registration offset coefficient to obtain a target auxiliary image;
and step S13, registering the target auxiliary image by using an enhanced spectrum grading registration algorithm to obtain the multi-scene image.
In the embodiment of the invention, the registration process of the OPS mode SAR data is divided into two parts of geometric registration and enhanced spectral diversity technology registration, firstly, the DEM is used for assisting to perform coarse registration between main and auxiliary images in a mode of searching for a same-name point, namely, the geometric registration, and on the basis of the geometric registration, the enhanced spectral diversity technology is used for performing fine registration on TOPS mode data to realize the registration precision of one thousandth of pixels in the azimuth direction.
Firstly, registering auxiliary DEM (Digital Elevation Model) and SAR data to be processed, realizing geocoding conversion of the DEM to a map coordinate system through geocoding under the assistance of satellite orbit data, realizing geometric transformation, selecting a series of points in uniform distribution of a main image, solving the coordinates of the points by using a distance equation and a Doppler equation and an earth Model, reversely calculating the pixel coordinates of the points in an image coordinate system of an auxiliary image, and calculating the offset by calculating the pixel coordinate difference of the same-name points between the main image and the auxiliary image.
Considering the efficiency problem, the offset of all pixels is obtained by establishing a Delaunay triangulation network interpolation method, and an initial registration offset coefficient is obtained. On the basis, resampling is carried out on the auxiliary image to obtain the target auxiliary image.
The coarse registration precision can only reach the pixel level, the precision required by interference needs 1/10 pixels, if the TOPS mode data needs to realize that the interference phase does not jump, the registration precision in the azimuth direction needs one thousandth of pixel precision, the common method cannot achieve the registration, and the high-precision registration is realized by designing and utilizing the spectrum diversity enhancement technology.
Considering that each adjacent Burst has a certain overlapping area when TOPS mode data acquisition is carried out, the relation between the registration error delta y and the interference phase difference of each pixel i overlapping with the Burst can be expressed as follows:
wherein,
expressed as the doppler center frequency variation in the overlap region,
for the sampling frequency of the azimuth direction, when the beam attitude of the TOPS mode data changes, the Doppler center can reach 5KHz in the edge region, the azimuth frequency is 486Hz, the interference phase obtained by the registered image does not have obvious phase jump between adjacent bursts, the interference phase difference needs to reach within 3 degrees, and therefore the registration precision needs to be controlled within one thousandth of pixel. The method comprises the following steps of firstly calculating to obtain delta y by utilizing an interference phase difference technology, calculating by utilizing an ESD algorithm through an average value of phase differences and an average value of Doppler center variation, and expressing a formula as follows:
wherein,
indicating that the mean is being taken.
The resulting registration accuracy for enhanced spectral diversity registration can be expressed as:
wherein
The coherence is expressed, in the practical production practice, the requirement of the registration accuracy cannot be realized once by considering the spectrum diversity enhancement technology, the design adopts a multi-iteration mode for calculation, through a large number of test experiments, the registration success can be ensured on the premise of ensuring the calculation efficiency by adopting the maximum 10 times of iteration spectrum diversity enhancement registration, and finally the obtained offset parameters ensure that the registration accuracy in the azimuth direction reaches one thousandth.
In the embodiment of the present invention, step S106 includes the following steps:
step S21, extracting the amplitude information of the multi-scene image;
step S22, calculating the average amplitude value and the amplitude standard deviation of each pixel in the multi-scene image by using the amplitude information;
step S23, calculating the ratio of the average amplitude value to the amplitude standard deviation of each pixel;
in step S24, the pixel with the ratio smaller than the first preset threshold is determined as the preselected point.
In the embodiment of the present invention, firstly, the amplitude information of the multi-view image is extracted, and the average amplitude value and the amplitude standard deviation of each pixel in the multi-view image are calculated by using the amplitude information, and the calculation formula is as follows:
wherein,
expressed as a standard deviation of the amplitude information,
is the amplitude mean.
Expressed as amplitude deviation values. With a number of practical analyses, it is preferable to set the first preset threshold value to 0.36. And when the pixel is smaller than the first preset threshold, determining the pixel to be changed as a scatterer preselected point.
In the embodiment of the present invention, step S108 includes the following steps:
step S31, performing phase correction processing on the pre-selected point to obtain a target pre-selected point, including:
step S32, calculating the space correlation component of the preset point by using a space adaptive filtering algorithm;
step S33, calculating the coherent value and the spatial uncorrelated error of the preset point by using the spatial correlation component;
step S34, determining the target preset point by using the coherence value of the preset point and a second preset threshold, where the second preset threshold is a ratio of the pseudo scatterers in the acceptable scatterer set.
In the embodiment of the invention, firstly, a grid is established, a multi-scene image is segmented, and the space correlation components of all scatterer preselected points in the grid are considered
Consistently, the grid size was chosen to be 100 x 100m2 by a number of practical analyses, taking into account image resolution.
Weighting calculation is carried out on the scatterer preselected points in each grid, and the calculation formula is as follows:
wherein
In order to be the weight, the weight is,
the amplitude deviation threshold value is used for the first time to define, namely 1/DA.
Calculating to obtain phase values in all grids, taking N-N grids as input data, and designing a space adaptive filtering method to obtain space related components in the region
。
Then, calculating the pre-selected point coherence value of the scatterer
And spatially uncorrelated view angle errors
。
Firstly, phase recombination is carried out on the differential interference phase to obtain a differential phase with the space correlation component removed.
Establishing a coherence factor
The calculation formula can be expressed as:
using a space search method to record when a certain one is
So that
At maximum, record the time
And
the value is obtained.
And finally, setting a second preset threshold according to the distribution of the coherence coefficients, and screening out the target scatterer.
Constructing a probability density function from the coherence values can be expressed as:
wherein,
is the weight coefficient of the weight of the image,
the representation is a scatterer or a volume of interest,
expressed as the probability of scatterers.
Designing and randomly generating 50 ten thousand noise points to obtain the probability density of the noise points
And calculating to obtain weighting coefficients
. Given on the basis of
The value (i.e. the second preset threshold), expressed as the proportion of false scatterers in the set of acceptable scatterers, can be expressed as:
on the basis of which the coherence coefficient can be determined
And (5) threshold value, and screening the scatterers to obtain the target scatterers.
In the embodiment of the present invention, step S110 includes the following steps:
step S41, performing the phase unwrapping process on the target scatterer to obtain an unwrapped phase;
step S42, calculating a target matrix by using a least square algorithm and the unwrapping phase, wherein the first row of the target matrix is the atmospheric and orbital error phase of the main image, and the third row is the linear deformation rate;
step S43, removing the atmospheric and orbit error phase in the unwrapping phase to obtain a target unwrapping phase;
and step S44, performing low-pass filtering processing on the target unwrapping phase to obtain linear deformation information of a region corresponding to the SAR data to be processed.
In an embodiment of the present invention, the expression of scatterer phase unwrapping is as follows:
then is provided with
An interference pair
A diffuser. Constructing a matrix of interfering relative vertical baseline and time baseline squares
. Unwrapping phase
Covariance matrix
And based on the covariance matrix, solving by using a least square method:
solving to obtain an object matrix
Wherein the object matrix
The first column of (a) represents the atmospheric and orbital error phase of the main image, and the third column the linear deformation rate.
And finally, carrying out low-pass filtering in a time domain on the unwrapping phase of the atmosphere and orbit errors with the spatial correlation visual angle errors and the main images removed, and extracting a low-frequency part to obtain a deformation phase.
And performing low-pass filtering and denoising on the high-frequency part obtained by the low-pass filtering in a spatial domain to obtain the atmospheric and orbital error of the auxiliary image.
Example two:
the embodiment of the invention also provides a high-precision settlement monitoring device for the SAR data, which is used for executing the high-precision settlement monitoring method for the SAR data provided by the embodiment of the invention, and the following is a specific introduction of the high-precision settlement monitoring device for the SAR data provided by the embodiment of the invention.
As shown in fig. 3, fig. 3 is a schematic diagram of the high-precision settlement monitoring device for SAR data, where the high-precision settlement monitoring for SAR data includes: an acquisition unit 10, a registration unit 20, a first determination unit 30, a second determination unit 40 and a processing unit 50.
The acquiring unit 10 is configured to acquire to-be-processed SAR data input by a user, where the to-be-processed SAR data is SAR data in a TOPS mode;
the registration unit 20 is configured to perform high-precision registration processing on the SAR data to be processed to obtain a multi-view image;
the first determining unit 30 is configured to determine a preselected point of a scatterer of the multi-view image by using an amplitude dispersion threshold algorithm;
the second determining unit 40 is configured to perform phase correction processing on the preselected point to obtain a target preset point, and determine a target scatterer by using the target preset point;
the processing unit 50 is configured to sequentially perform phase unwrapping processing and atmospheric and orbit correction processing on the target scatterer to obtain linear deformation information of a region corresponding to the to-be-processed SAR data.
In the embodiment of the invention, SAR data to be processed input by a user is acquired; carrying out high-precision registration processing on SAR data to be processed to obtain a multi-scene image; determining preselected points of scatterers of the multi-scene images by using an amplitude dispersion threshold algorithm; carrying out phase correction processing on the preselected points to obtain target preset points, and determining a target scatterer by using the target preset points; the phase unwrapping processing and the atmosphere and orbit correction processing are sequentially carried out on the target scatterer, so that the linear deformation information of the region corresponding to the SAR data to be processed is obtained, the purpose of monitoring the linear deformation information of the monitored region is achieved, the technical problem that the monitoring precision of the existing settlement monitoring method for the SAR data on the linear deformation information of the monitored region is low is solved, and the technical effect of carrying out high-precision monitoring on the linear deformation information of the monitored region is achieved.
Example three:
an embodiment of the present invention further provides an electronic device, including a memory and a processor, where the memory is used to store a program that supports the processor to execute the method described in the first embodiment, and the processor is configured to execute the program stored in the memory.
Referring to fig. 4, an embodiment of the present invention further provides an electronic device 100, including: a processor 60, a memory 61, a bus 62 and a communication interface 63, wherein the processor 60, the communication interface 63 and the memory 61 are connected through the bus 62; the processor 60 is arranged to execute executable modules, such as computer programs, stored in the memory 61.
The Memory 61 may include a high-speed Random Access Memory (RAM) and may also include a non-volatile Memory (non-volatile Memory), such as at least one disk Memory. The communication connection between the network element of the system and at least one other network element is realized through at least one communication interface 63 (which may be wired or wireless), and the internet, a wide area network, a local network, a metropolitan area network, and the like can be used.
The bus 62 may be an ISA bus, PCI bus, EISA bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one double-headed arrow is shown in FIG. 4, but that does not indicate only one bus or one type of bus.
The memory 61 is used for storing a program, the processor 60 executes the program after receiving an execution instruction, and the method executed by the apparatus defined by the flow process disclosed in any of the foregoing embodiments of the present invention may be applied to the processor 60, or implemented by the processor 60.
The processor 60 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware or instructions in the form of software in the processor 60. The Processor 60 may be a general-purpose Processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; the device can also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field-Programmable Gate Array (FPGA), or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components. The various methods, steps and logic blocks disclosed in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The storage medium is located in a memory 61, and the processor 60 reads the information in the memory 61 and, in combination with its hardware, performs the steps of the above method.
Example four:
the embodiment of the present invention further provides a computer-readable storage medium, where a computer program is stored on the computer-readable storage medium, and when the computer program is executed by a processor, the computer program performs the steps of the method in the first embodiment.
In addition, in the description of the embodiments of the present invention, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, and are only for convenience of description and simplicity of description, but do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units is only one logical division, and there may be other divisions when actually implemented, and for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection of devices or units through some communication interfaces, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (10)
1. A high-precision settlement monitoring method for SAR data is characterized by comprising the following steps:
acquiring SAR data to be processed input by a user, wherein the SAR data to be processed is SAR data in a TOPS mode;
carrying out high-precision registration processing on the SAR data to be processed to obtain a multi-scene image;
determining preselected points of scatterers of the multi-scene images by using an amplitude dispersion threshold algorithm;
carrying out phase correction processing on the preselected points to obtain target preset points, and determining a target scatterer by using the target preset points;
and sequentially carrying out phase unwrapping processing and atmospheric and orbit correction processing on the target scatterer to obtain linear deformation information of an area corresponding to the SAR data to be processed.
2. The method of claim 1, wherein the SAR data to be processed comprises: the method comprises the following steps of carrying out high-precision registration processing on SAR data to be processed by main image data and auxiliary image data to obtain a multi-scene image, wherein the high-precision registration processing comprises the following steps:
performing geometric transformation on the SAR data to be processed by using a digital elevation model and satellite orbit data to obtain an initial registration offset coefficient of the SAR data to be processed, wherein the initial registration offset coefficient is used for representing the coordinate difference of the pixels of the same-name points between the main image and the auxiliary image;
resampling the auxiliary image by using the initial registration offset coefficient to obtain a target auxiliary image;
and registering the target auxiliary image by using an enhanced spectrum grading registration algorithm to obtain the multi-scene image.
3. The method of claim 1, wherein determining preselected points of scatterers of the multi-view image using an amplitude dispersion threshold algorithm comprises:
extracting amplitude information of the multi-scene image;
calculating the average amplitude value and the amplitude standard deviation of each pixel in the multi-scene image by using the amplitude information;
calculating the ratio of the average amplitude value to the amplitude standard deviation of each pixel;
and determining the pixel of which the ratio is smaller than a first preset threshold value as the preselected point.
4. The method of claim 1, wherein performing a phase correction process on the preselected point to obtain a target preset point comprises:
calculating the spatial correlation component of the preset point by using a spatial adaptive filtering algorithm;
calculating a coherent value and a spatial uncorrelated error of the preset point by using the spatial correlation component;
and determining a target preset point by using the coherence value of the preset point and a second preset threshold, wherein the second preset threshold is the proportion of the pseudo scatterers in the acceptable scatterer set.
5. The method according to claim 2, wherein the step of sequentially performing phase unwrapping processing and atmospheric and orbit correction processing on the target scatterer to obtain linear deformation information of a region corresponding to the SAR data to be processed includes:
performing the phase unwrapping processing on the target scatterer to obtain an unwrapped phase;
calculating a target matrix by using a least square algorithm and the unwrapping phase, wherein the first column of the target matrix is an atmospheric and orbital error phase of the main image, and the third column is a linear deformation rate;
removing the atmospheric and orbit error phase in the unwrapping phase to obtain a target unwrapping phase;
and carrying out low-pass filtering processing on the target unwrapping phase to obtain linear deformation information of the region corresponding to the SAR data to be processed.
6. The method according to claim 5, wherein the low-pass filtering processing is performed on the target unwrapping phase to obtain linear deformation information of a region corresponding to the SAR data to be processed, and the method comprises the following steps:
performing low-pass filtering processing in a time domain on the target unwrapping phase to obtain a first filter, wherein the first filter includes: high-frequency filtering and low-frequency filtering;
determining the low-frequency filtering as a deformation phase of a region corresponding to the SAR data to be processed;
and carrying out low-pass filtering processing and denoising processing in a spatial domain on the high-frequency filtering to obtain the atmospheric and orbital errors of the auxiliary image.
7. A high accuracy settlement monitoring device of SAR data, the device comprising: an acquisition unit, a registration unit, a first determination unit, a second determination unit and a processing unit, wherein,
the acquisition unit is used for acquiring SAR data to be processed input by a user, wherein the SAR data to be processed is SAR data in a TOPS mode;
the registration unit is used for carrying out high-precision registration processing on the SAR data to be processed to obtain a multi-scene image;
the first determining unit is used for determining preselected points of the scatterers of the multi-scene images by using an amplitude dispersion threshold algorithm;
the second determining unit is used for performing phase correction processing on the preselected point to obtain a target preset point and determining a target scatterer by using the target preset point;
and the processing unit is used for sequentially carrying out phase unwrapping processing and atmospheric and orbit correction processing on the target scatterer to obtain linear deformation information of an area corresponding to the SAR data to be processed.
8. The apparatus of claim 7, wherein the to-be-processed SAR data comprises: the registration unit is configured to:
performing geometric transformation on the SAR data to be processed by using a digital elevation model and satellite orbit data to obtain an initial registration offset coefficient of the SAR data to be processed, wherein the initial registration offset coefficient is used for representing the coordinate difference of the pixels of the same-name points between the main image and the auxiliary image;
resampling the auxiliary image by using the initial registration offset coefficient to obtain a target auxiliary image;
and registering the target auxiliary image by using an enhanced spectrum grading registration algorithm to obtain the multi-scene image.
9. An electronic device comprising a memory for storing a program that enables a processor to perform the method of any of claims 1 to 6 and a processor configured to execute the program stored in the memory.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the method of any one of the preceding claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011226759.3A CN112034463A (en) | 2020-11-06 | 2020-11-06 | High-precision settlement monitoring method and device for SAR data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011226759.3A CN112034463A (en) | 2020-11-06 | 2020-11-06 | High-precision settlement monitoring method and device for SAR data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112034463A true CN112034463A (en) | 2020-12-04 |
Family
ID=73572845
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011226759.3A Pending CN112034463A (en) | 2020-11-06 | 2020-11-06 | High-precision settlement monitoring method and device for SAR data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112034463A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112986949A (en) * | 2021-04-21 | 2021-06-18 | 北京东方至远科技股份有限公司 | SAR high-precision time sequence deformation monitoring method and device for diagonal reflector |
| CN113884054A (en) * | 2021-09-29 | 2022-01-04 | 广东电网有限责任公司 | Tower displacement settlement monitoring method and system |
| CN117632264A (en) * | 2024-01-24 | 2024-03-01 | 中国科学院空天信息创新研究院 | SAR data rapid parallel processing method |
-
2020
- 2020-11-06 CN CN202011226759.3A patent/CN112034463A/en active Pending
Non-Patent Citations (6)
| Title |
|---|
| 卢学辉: "PS-InSAR技术在北京典型地区地面沉降监测中的应用研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 康琪: "TOPS模式数据InSAR形变监测配准方法", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 武松超: "基于PS-InSAR技术的地面沉降监测研究及应用", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 裴媛媛 等: "利用时序DInSAR监测填海造陆地区地表沉降", 《武汉大学学报 信息科学版》 * |
| 高二涛: "基于高分辨率SAR影像监测高速铁路沿线形变研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
| 黄其欢 等: "Sentinel-1 TOPS影像干涉处理与形变监测", 《大地测量与地球动力学》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112986949A (en) * | 2021-04-21 | 2021-06-18 | 北京东方至远科技股份有限公司 | SAR high-precision time sequence deformation monitoring method and device for diagonal reflector |
| CN113884054A (en) * | 2021-09-29 | 2022-01-04 | 广东电网有限责任公司 | Tower displacement settlement monitoring method and system |
| CN113884054B (en) * | 2021-09-29 | 2023-09-26 | 广东电网有限责任公司 | Tower displacement settlement monitoring method and system |
| CN117632264A (en) * | 2024-01-24 | 2024-03-01 | 中国科学院空天信息创新研究院 | SAR data rapid parallel processing method |
| CN117632264B (en) * | 2024-01-24 | 2024-05-07 | 中国科学院空天信息创新研究院 | SAR data rapid parallel processing method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112034463A (en) | High-precision settlement monitoring method and device for SAR data | |
| Xu et al. | A refined strategy for removing composite errors of SAR interferogram | |
| US8587471B2 (en) | Process for identifying statistically homogeneous pixels in SAR images acquired on the same area | |
| CN109752715B (en) | SAR data total-dispersion body detection method and device | |
| CN106772380A (en) | A kind of circumferential synthetic aperture radar imaging method | |
| CN111239736B (en) | Single-baseline-based surface elevation correction method, device, equipment and storage medium | |
| CN115629384B (en) | Correction method of time sequence InSAR error and related equipment | |
| CN111709876B (en) | Image splicing method, device, equipment and storage medium | |
| CN108919350A (en) | Diffracted wave separation method and device | |
| CN110297243B (en) | Phase error compensation method and device in synthetic aperture radar tomography three-dimensional imaging | |
| CN109886910B (en) | DEM (digital elevation model) correction method and device for external digital elevation model | |
| CN112363165B (en) | Method, device, equipment and medium for forest subsurface shape inversion | |
| CN111965645B (en) | Multi-baseline vegetation height inversion method and device considering geometric constraints | |
| CN110657742A (en) | Aquifer deformation signal separation method, aquifer deformation signal separation device, aquifer deformation signal separation equipment and readable storage medium | |
| CN110297242A (en) | Compressed sensing based synthetic aperture radar chromatography three-D imaging method and device | |
| CN111222544A (en) | Ground simulation test system for influence of satellite flutter on camera imaging | |
| CN112419205B (en) | Goldstein pyramid construction method for SAR interferogram processing | |
| CN113554011B (en) | Method, device and equipment for detecting change of remote sensing image and storage medium | |
| WO2022055608A2 (en) | Geometric tools and methods to measure closure phase for robust feature recognition in interferometric images | |
| CN111738135A (en) | SAR image feature extraction method considering slant range and geographic coordinate mapping deviation | |
| CN111681299A (en) | Method and device for generating digital surface model based on InSAR unwrapping phase | |
| CN116342438B (en) | Spatial spectrum self-adaptive filtering method and device | |
| CN116028480B (en) | Filling method for improving space-time coverage of remote sensing soil moisture product | |
| CN115391721B (en) | Magnetic moment calculation method and device based on magnetic total field vertical gradient and related equipment | |
| Goodrick | Image reconstruction in radio astronomy with non-coplanar synthesis arrays |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201204 |
|
| RJ01 | Rejection of invention patent application after publication |