CN111522005A - Deformation monitoring and terrain reconstruction method - Google Patents
Deformation monitoring and terrain reconstruction method Download PDFInfo
- Publication number
- CN111522005A CN111522005A CN202010504347.5A CN202010504347A CN111522005A CN 111522005 A CN111522005 A CN 111522005A CN 202010504347 A CN202010504347 A CN 202010504347A CN 111522005 A CN111522005 A CN 111522005A
- Authority
- CN
- China
- Prior art keywords
- monitoring
- radar
- deformation
- receiving antenna
- terrain reconstruction
- 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
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
- 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
-
- 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/86—Combinations of radar systems with non-radar systems, e.g. sonar, direction finder
-
- 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/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
-
- 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/411—Identification of targets based on measurements of radar reflectivity
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention provides a deformation monitoring and terrain reconstruction method. The method comprises the steps of obtaining a radar image through a radar monitoring device provided with two or more receiving antennas; carrying out time sequence differential interference processing on the radar images received by each receiving antenna to obtain deformation information of a monitored scene; erecting a plurality of corner reflectors in a monitoring scene range, carrying out interferometric synthetic aperture radar imaging on radar images acquired by different receiving antennas, and carrying out terrain reconstruction on a monitoring scene; and measuring position information of the centers of the corner reflectors through the GNSS to obtain error correction parameters of the elevation measured by the radar monitoring device, and finally outputting corrected terrain reconstruction information. The method and the device integrate the deformation information received by each receiving antenna to obtain the deformation information of the side slope, and simultaneously perform terrain reconstruction and correction on radar images acquired by different receiving antennas to complete the monitoring of the side slope.
Description
Technical Field
The invention relates to the technical field of radar monitoring, in particular to a deformation monitoring and terrain reconstruction method.
Background
Slope geological disaster early warning is one of the key tasks of applications such as mine monitoring, geological survey, emergency rescue always, and for conventional monitoring technology means such as displacement meter, big dipper difference monitoring facilities, surveyor's level, the mode that adopts ground radar to carry out remote sensing monitoring does not receive weather influences such as rain, snow, fog and the like and illumination influences such as highlight, night, measurement accuracy reaches the submillimeter level, can carry out comprehensive and quick measurement to scene on a large scale, and the comprehensive use cost is low, is one kind promising side slope deformation monitoring technology means.
The existing terrain reconstruction algorithm is mainly as follows: after obtaining the data of the two receiving channels, two multi-view radar images can be formed, and one of the multi-view radar images can be generally defined as a main image and the other as an auxiliary image. Through the main image and the auxiliary image of the radar, the elevation can be estimated by utilizing the phase interference technology, as shown in fig. 4, but when the baseline distance between the receiving antennas does not meet the condition that the distance is far less than the slant distance, a large error exists when the elevation is calculated by using the method, and the calculation precision requirement cannot be met.
In view of the above, there is a need for a deformation monitoring and terrain reconstruction method to solve the problems in the prior art.
Disclosure of Invention
The invention aims to provide a deformation monitoring and terrain reconstruction method to solve the problems of deformation monitoring and terrain reconstruction in slope monitoring.
In order to achieve the above object, the present invention provides a deformation monitoring and terrain reconstruction method, comprising the following steps:
step A: acquiring a radar image through a radar monitoring device, wherein the radar monitoring device is provided with two or more than two receiving antennas;
and B: carrying out time sequence differential interference processing on radar images received in the monitoring process of each receiving antenna, and then integrating real deformation values obtained by each receiving antenna to form deformation information of a monitoring scene;
and C: erecting a plurality of corner reflectors in a monitoring scene range, performing terrain reconstruction in the monitoring range through a radar monitoring device, performing interferometric synthetic aperture radar imaging on radar images acquired by different receiving antennas to form elevation information, and performing terrain reconstruction on a monitoring scene;
step D: and C, measuring position information of the centers of the corner reflectors through the GNSS to obtain error correction parameters of the elevation measured by the radar monitoring device, correcting the error of the elevation information in the step C, and finally outputting corrected terrain reconstruction information.
Preferably, in the step B, the number of the receiving antennas is n, and n is more than or equal to 2; deformation value delta D corresponding to ith receiving antennaiThe radar image phase interference method is obtained by phase interference calculation of two adjacent radar images, and adopts an expression 1) to calculate:
where c is the speed of light, fcFor the radar operating frequency, phitAnd phit-1Phase values measured at the current moment and the previous moment respectively; radar carrier wavelength is λ, when | Δ DiWhen | ≦ λ, Δ DiThe true deformation value corresponding to the ith receiving antenna is equal to delta DiFor Δ D, | > λiAnd after phase unwrapping, obtaining a true deformation value corresponding to the ith receiving antenna.
Preferably, in the step B, the deformation information D of the monitoring scene is calculated by adopting an expression 2):
preferably, in said step CAdopting an expression 7) to carry out terrain reconstruction: selecting two receiving antennas S1And S2Performing terrain reconstruction on the monitoring point T to obtain the elevation h of the monitoring point TT:
Wherein H1Is a receiving antenna S1Height of center of (R)1And R2Are respectively a receiving antenna S1And S2Distance to monitoring point T, B being S1To S2Length of the base line of (e), theta being the receiving antenna S1Center to ground connection and receiving antenna S1The angle between the two connecting lines to the monitoring point T is β, which is the attitude angle between the base line and the horizontal plane.
Preferably, in the step D, the number of the corner reflectors is m, and m is more than or equal to 3; obtaining the elevation error delta h of the jth corner reflector by the expression 8)j:
Δhj=hGj-hJj8);
Wherein h isGjFor the elevation of the jth corner reflector measured by GNSS, hJjThe elevation of the jth angular reflector obtained by topographic reconstruction of expression 7);
j th corner reflector to receiving antenna S1And S2Has a mean distance of LjThen the elevation correction parameter delta h of the terrain reconstructionTUsing expression 9) to calculate:
altitude value H after error correctionTUsing expression 10) calculate:
HT=hT+ΔhT10)。
preferably, in the step C, the receiving antenna S1And S2Distance R to monitoring point T1And R2Is obtained by a synthetic aperture imaging algorithm.
Preferably, in the step D, coordinates of the center of the corner reflector in the WGS84 coordinate system are measured through GNSS, and coordinates of the center of the corner reflector in the radar monitoring coordinate system are measured through a radar monitoring device, so as to obtain affine transformation parameters of coordinates in the WGS84 coordinate system and the radar monitoring coordinate system, and coordinates of all radar monitoring points in the radar monitoring coordinate system are converted into longitude and latitude coordinates in the WGS84 coordinate system.
The technical scheme of the invention has the following beneficial effects:
(1) according to the method, a plurality of receiving antennas are arranged, time sequence differential interference processing is carried out on radar images received in the monitoring process of each receiving antenna, deformation information received by each receiving antenna is integrated to obtain deformation information of a side slope, meanwhile, interference synthetic aperture radar imaging is carried out on radar images obtained by different receiving antennas to form elevation information and carry out topographic reconstruction on a scene, and monitoring of the side slope is completed.
(2) In the invention, the slope deformation information monitored by a single receiving antenna can be quickly obtained by performing phase interference calculation on two frames of radar images of the single receiving antenna.
(3) In the invention, the final deformation information D of the monitoring scene is obtained by comprehensively solving the real deformation values of all the receiving antennas, so that the reliability of the deformation information can be improved, and the monitoring error can be reduced.
(4) In the invention, two receiving antennas S are selected1And S2The elevation value of the monitoring point is solved through the cosine law, the process is simple and convenient, the constraint of the baseline distance and the slope distance is avoided, and the application range of the terrain reconstruction algorithm is expanded.
(5) According to the invention, by arranging the corner reflector and measuring the elevation of the corner reflector by using the GNSS, the elevation correction parameter of terrain reconstruction can be obtained, the elevation value in the terrain reconstruction can be corrected and then output, and the accuracy of the terrain reconstruction is improved.
(6) In the invention, after the synthetic aperture imaging algorithm is used for calculation, the phase unwrapping and the error correction are carried out to obtain the slant distance from the receiving antenna to the monitoring point, so that the precision of the slant distance measurement can be improved.
(7) In the invention, the coordinates of the corner reflector during GNSS measurement and the coordinates of the corner reflector during radar monitoring device measurement are compared to obtain affine transformation parameters of the coordinates, and the coordinates of all radar monitoring points under a radar monitoring coordinate system can be converted into longitude and latitude coordinates under a WGS84 coordinate system.
In addition to the objects, features and advantages described above, other objects, features and advantages of the present invention are also provided. The present invention will be described in further detail below with reference to the drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the invention and, together with the description, serve to explain the invention and not to limit the invention. In the drawings:
FIG. 1 is a schematic structural diagram of a radar monitoring device in an embodiment of the present application;
FIG. 2 is a side view of a radar monitoring device in an embodiment of the present application;
FIG. 3 is a schematic diagram of a terrain reconstruction algorithm in an embodiment of the present application;
FIG. 4 is a flow chart of prior art terrain reconstruction;
the antenna comprises a base 1, a base 2, a rotating shaft 3, a rotating arm 4, an antenna frame 5, a transmitting antenna 6 and a receiving antenna.
Detailed Description
Embodiments of the invention will be described in detail below with reference to the drawings, but the invention can be implemented in many different ways, which are defined and covered by the claims.
Example (b):
referring to fig. 1 to 4, a deformation monitoring and terrain reconstruction method is applied to slope monitoring in the present embodiment.
A deformation monitoring and terrain reconstruction method adopts a radar monitoring device as shown in figure 1, the radar monitoring device comprises a rotating shaft 2 arranged on a base 1 and a rotating arm 3 connected with the rotating shaft 2, and the rotating arm 3 can rotate 360 degrees around the rotating shaft 2 and can also rotate in a designated sector area according to measurement requirements; with pivoted arms 3 remote from axis 2One end of the antenna frame 4 is hinged with the antenna frame 4, the pitching angle of the antenna frame 4 can be adjusted relative to the rotating arm 3, and the adjusting range is +/-45 degrees; the antenna frame 4 is provided with a transmitting antenna 5 and a receiving antenna 6, and in this embodiment, the radar monitoring device is provided with one transmitting antenna 5 and two receiving antennas 6 (receiving antennas S respectively)1And a receiving antenna S2) As shown in fig. 2; in this embodiment, the transmitting antenna 5 and the two receiving antennas 6 both adopt square horn antennas, wherein the two receiving antennas 6 are arranged in a vertical column, and the distance between the centers of the two receiving antennas 6 is at least twice the side length of the square horn antenna; the center of the transmitting antenna 5 is arranged on the perpendicular bisector of the central lines of the two receiving antennas 6, and the distance between the central point of the transmitting antenna 5 and the connecting line of the central points of the receiving antennas 6 is not less than 3 times of the side length of the antennas, so that interference between signals received by the receiving antennas 6 can not be formed, a base line between the receiving antennas 6 is formed in a measuring range, the terrain reconstruction can be effectively carried out, and the accuracy requirement is met. The method for monitoring the slope deformation and reconstructing the terrain by adopting the radar monitoring device comprises the following steps:
step A: acquiring a radar image through a radar monitoring device, wherein the radar monitoring device is provided with two or more than two receiving antennas;
and B: carrying out time sequence differential interference processing on radar images received in the monitoring process of each receiving antenna, and synthesizing real deformation values obtained by each receiving antenna to form deformation information of a monitoring scene;
and C: erecting a plurality of corner reflectors in a monitoring scene range, performing terrain reconstruction in the monitoring range through a radar monitoring device, performing interferometric synthetic aperture radar imaging on radar images acquired by different receiving antennas to form elevation information, and performing terrain reconstruction on a monitoring scene;
step D: and C, measuring position information of the centers of the corner reflectors through the GNSS to obtain error correction parameters of the elevation measured by the radar monitoring device, correcting the errors of the elevation information in the step C, and finally outputting terrain reconstruction information.
Wherein, the number of the receiving antennas in the step B is n, and n is more than or equal to 2; shape corresponding to ith receiving antennaVariable value Delta DiThe radar image is obtained by phase interference calculation from two adjacent frames of radar images, and the receiving antenna S is calculated by adopting an expression 1)1And a receiving antenna S2Measured deformation values:
where c is the speed of light, fcFor the radar operating frequency, phitAnd phit-1Phase values measured at the current moment and the previous moment respectively; radar carrier wavelength is λ, when | Δ DiWhen | ≦ λ, Δ DiThe true deformation value corresponding to the ith receiving antenna is equal to delta DiFor Δ D, | > λiAnd after phase unwrapping, obtaining a true deformation value corresponding to the ith receiving antenna.
The true deformation value corresponding to the ith antenna (i ═ 1, 2.... n) of the radar monitoring device is Δ DiAnd averaging the real deformation values of the n receiving antennas, wherein the output deformation information D of the monitoring scene is as follows:
respectively to receiving antennas S1And a receiving antenna S2Carrying out time series differential interference processing on the received radar image, and obtaining a receiving antenna S through expressions 1) to 2)1The obtained true deformation value Delta D1Obtaining a receiving antenna S2The obtained true deformation value Delta D2And finally outputting the deformation information D of the monitoring scene (side slope) as follows:
step C, selecting a geological stable point in a radar monitoring range, erecting more than 3 corner reflectors, and measuring the longitude and latitude and the elevation under a WGS84 coordinate system (World geodetic System-1984 CoordinateSystem) at the center of each corner reflector through GNSS or a total station and the like;
when the radar monitoring range is determined, the angular resolution and the distance resolution of the radar monitoring device are unchanged, so that the position of each radar monitoring point is also fixed. Each frame of radar image is firstly subjected to PS (permanent scattering) point screening on a monitoring area through a plurality of thresholds of comprehensive coherence coefficients, signal-to-noise ratios and phase stability, and the position of the corner reflector in a radar monitoring coordinate system is determined. The radar monitoring coordinate system is a polar coordinate system taking a radar rotating shaft as a center;
calculating coordinate affine transformation parameters converted from a radar monitoring coordinate system to a WGS84 coordinate system through coordinates of each corner reflector in the radar monitoring coordinate system and coordinates of the WGS84 coordinate system; converting coordinates of monitoring points in all radar monitoring scene ranges in a radar monitoring coordinate system into longitude and latitude coordinates in a WGS84 coordinate system through coordinate affine transformation parameters, and outputting the longitude and latitude coordinates as longitude and latitude coordinate results of the monitoring points;
by the terrain reconstruction method, the elevation values (including the elevation values of the corner reflectors) of all monitoring points in the radar monitoring scene range are calculated: as shown in fig. 3, the positive direction of the x-axis is directed out of the paper, the positive direction of the y-axis is horizontally to the right, and the positive direction of the z-axis is vertically upward; receiving antenna S1And a receiving antenna S2The slope distances to the monitoring points T are respectively R1And R2The distance value is obtained by phase unwrapping and error correction after the distance value is obtained through a synthetic aperture imaging algorithm; h1Is a receiving antenna S1The center height of (d); b is S1To S2Length of the base line of (e), theta being the receiving antenna S1Center to ground connection and receiving antenna S1An angle between the two connecting lines to the monitoring point T, β is an attitude angle formed by a base line and a horizontal plane, and gamma is a receiving antenna S1For the depression angle of the monitoring point T, the elevation h of the monitoring point TTComprises the following steps:
hT=H1-R1cosθ 4);
θ=90°-γ 5);
simultaneous expressions 4) to 6) can be obtained:
in the step D), the elevation of each monitoring point in the radar monitoring scene range is obtained through an expression 7), wherein the number of the corner reflectors is m, and m is more than or equal to 3; obtaining the elevation error delta h of the jth corner reflector by the expression 8)j:
Δhj=hGj-hJj8);
Wherein h isGjFor the elevation of the jth corner reflector measured by GNSS, hJjThe elevation of the jth angular reflector obtained by topographic reconstruction of expression 7); because the elevation value measured by the GNSS has higher precision in hGjIs true value;
j th corner reflector to receiving antenna S1And S2Has a mean distance of L j1,2, the elevation correction parameter Δ h for terrain reconstructionTComprises the following steps:
altitude value H after error correctionTComprises the following steps:
HT=hT+ΔhT10)。
the final output information of the radar monitoring device comprises deformation information and terrain reconstruction information of a monitoring scene, and comprises the following steps: the deformation D (the direction far away from the sight line of the radar is positive, and the direction close to the sight line of the radar is negative) of the monitoring point in the radar monitoring range and the elevation H of the monitoring pointTAnd the longitude and latitude coordinates of the monitoring point under the WGS84 coordinate system.
The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention, and various modifications and changes may be made by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (7)
1. A deformation monitoring and terrain reconstruction method is characterized by comprising the following steps:
step A: acquiring a radar image through a radar monitoring device, wherein the radar monitoring device is provided with two or more than two receiving antennas;
and B: carrying out time sequence differential interference processing on radar images received in the monitoring process of each receiving antenna, and then integrating real deformation values obtained by each receiving antenna to form deformation information of a monitoring scene;
and C: erecting a plurality of corner reflectors in a monitoring scene range, performing terrain reconstruction in the monitoring range through a radar monitoring device, performing interferometric synthetic aperture radar imaging on radar images acquired by different receiving antennas to form elevation information, and performing terrain reconstruction on a monitoring scene;
step D: and C, measuring position information of the centers of the corner reflectors through the GNSS to obtain error correction parameters of the elevation information measured by the radar monitoring device, correcting the errors of the elevation information in the step C, and finally outputting corrected terrain reconstruction information.
2. A deformation monitoring and terrain reconstructing method as claimed in claim 1, wherein in the step B, the number of receiving antennas is n, n is greater than or equal to 2; deformation value delta D corresponding to ith receiving antennaiThe radar image phase interference method is obtained by phase interference calculation of two adjacent radar images, and adopts an expression 1) to calculate:
where c is the speed of light, fcFor the radar operating frequency, phitAnd phit-1Phase values measured at the current moment and the previous moment respectively; radar carrier wavelength is λ, when | Δ DiWhen | ≦ λ, Δ DiThe true deformation value corresponding to the ith receiving antenna is equal to | Δ DiFor Δ D, | > λiAnd after phase unwrapping, obtaining a true deformation value corresponding to the ith receiving antenna.
4. a deformation monitoring and terrain reconstruction method according to claim 1, wherein in the step C, the terrain reconstruction is performed by using expression 7): selecting two receiving antennas S1And S2Performing terrain reconstruction on the monitoring point T to obtain the elevation h of the monitoring point TT:
Wherein H1Is a receiving antenna S1Height of center of (R)1And R2Are respectively a receiving antenna S1And S2Distance to monitoring point T, B being S1To S2Length of the base line of (e), theta being the receiving antenna S1Center to ground connection and receiving antenna S1The angle between the two connecting lines to the monitoring point T is β, which is the attitude angle between the base line and the horizontal plane.
5. A deformation monitoring and terrain reconstruction method as set forth in claim 4, wherein in the step D, the number of corner reflectors is m, and m is greater than or equal to 3; obtaining the elevation error delta h of the jth corner reflector by the expression 8)j:
Δhj=hGj-hJj8);
Wherein h isGjFor the elevation of the jth corner reflector measured by GNSS, hJjIs the jth terrain reconstruction obtained by expression 7)The elevation of the corner reflector;
j th corner reflector to receiving antenna S1And S2Has a mean distance of LjThen the elevation correction parameter delta h of the terrain reconstructionTUsing expression 9) to calculate:
altitude value H after error correctionTUsing expression 10) calculate:
HT=hT+ΔhT10)。
6. a deformation monitoring and terrain reconstruction method as set forth in claim 4, wherein in the step C, the receiving antenna S1And S2Distance R to monitoring point T1And R2Is obtained by a synthetic aperture imaging algorithm.
7. A deformation monitoring and terrain reconstruction method as set forth in claim 1, wherein in the step D, coordinates of the center of the corner reflector in the WGS84 coordinate system are measured by GNSS, and coordinates of the center of the corner reflector in the radar monitoring coordinate system are measured by a radar monitoring device, affine transformation parameters of the coordinates in the WGS84 coordinate system and the radar monitoring coordinate system are obtained, and coordinates of all radar monitoring points in the radar monitoring coordinate system are converted into longitude and latitude coordinates in the WGS84 coordinate system.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010504347.5A CN111522005A (en) | 2020-06-05 | 2020-06-05 | Deformation monitoring and terrain reconstruction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010504347.5A CN111522005A (en) | 2020-06-05 | 2020-06-05 | Deformation monitoring and terrain reconstruction method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111522005A true CN111522005A (en) | 2020-08-11 |
Family
ID=71909488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010504347.5A Pending CN111522005A (en) | 2020-06-05 | 2020-06-05 | Deformation monitoring and terrain reconstruction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111522005A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112883000A (en) * | 2021-03-17 | 2021-06-01 | 中国有色金属长沙勘察设计研究院有限公司 | Deformation monitoring radar data file storage method |
CN113534081A (en) * | 2021-08-17 | 2021-10-22 | 中国有色金属长沙勘察设计研究院有限公司 | Detection method and device for deformation monitoring radar precision |
CN114428247A (en) * | 2022-03-11 | 2022-05-03 | 深圳航天科技创新研究院 | Single antenna ultra-wideband radar system for imaging applications |
CN114545348A (en) * | 2022-02-25 | 2022-05-27 | 中电科技扬州宝军电子有限公司 | SVD-based radar system error calibration method |
CN114966601A (en) * | 2022-08-01 | 2022-08-30 | 南京隼眼电子科技有限公司 | Mountain landslide prediction method based on millimeter wave radar and electronic equipment |
CN115079166A (en) * | 2022-07-27 | 2022-09-20 | 南京隼眼电子科技有限公司 | Millimeter wave radar disaster monitoring method and system and electronic equipment |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182339A (en) * | 2015-09-25 | 2015-12-23 | 昆明理工大学 | Method for correcting environmental influences at slope deformation monitoring on the basis of corner reflector |
CN106093938A (en) * | 2016-05-17 | 2016-11-09 | 长安大学 | A kind of mining area based on manual corner reflector side-play amount deformation monitoring method |
CN110618425A (en) * | 2019-08-19 | 2019-12-27 | 中国电力科学研究院有限公司 | Baseline phase determination method and system for ground-based interferometric radar in deformation monitoring |
-
2020
- 2020-06-05 CN CN202010504347.5A patent/CN111522005A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182339A (en) * | 2015-09-25 | 2015-12-23 | 昆明理工大学 | Method for correcting environmental influences at slope deformation monitoring on the basis of corner reflector |
CN106093938A (en) * | 2016-05-17 | 2016-11-09 | 长安大学 | A kind of mining area based on manual corner reflector side-play amount deformation monitoring method |
CN110618425A (en) * | 2019-08-19 | 2019-12-27 | 中国电力科学研究院有限公司 | Baseline phase determination method and system for ground-based interferometric radar in deformation monitoring |
Non-Patent Citations (3)
Title |
---|
李志林等: "合成孔径雷达的干涉测量法", 《数字高程模型 第3版》 * |
邢学敏: "CRInSAR与PSInSAR联合监测矿区时序地表形变研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅰ辑》 * |
邢学敏等: "CRInSAR与PSInSAR融合探测地表线性变形速率", 《大地测量与地球动力学》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112883000A (en) * | 2021-03-17 | 2021-06-01 | 中国有色金属长沙勘察设计研究院有限公司 | Deformation monitoring radar data file storage method |
CN112883000B (en) * | 2021-03-17 | 2022-04-15 | 中国有色金属长沙勘察设计研究院有限公司 | Deformation monitoring radar data file storage method |
CN113534081A (en) * | 2021-08-17 | 2021-10-22 | 中国有色金属长沙勘察设计研究院有限公司 | Detection method and device for deformation monitoring radar precision |
WO2023108544A1 (en) * | 2021-12-15 | 2023-06-22 | 深圳航天科技创新研究院 | Single-antenna ultra-wideband radar system for imaging application |
CN114545348A (en) * | 2022-02-25 | 2022-05-27 | 中电科技扬州宝军电子有限公司 | SVD-based radar system error calibration method |
CN114428247A (en) * | 2022-03-11 | 2022-05-03 | 深圳航天科技创新研究院 | Single antenna ultra-wideband radar system for imaging applications |
CN114428247B (en) * | 2022-03-11 | 2022-09-27 | 深圳航天科技创新研究院 | Single antenna ultra-wideband radar system for imaging applications |
CN115079166A (en) * | 2022-07-27 | 2022-09-20 | 南京隼眼电子科技有限公司 | Millimeter wave radar disaster monitoring method and system and electronic equipment |
CN115079166B (en) * | 2022-07-27 | 2022-11-01 | 南京隼眼电子科技有限公司 | Millimeter wave radar disaster monitoring method and system and electronic equipment |
CN114966601A (en) * | 2022-08-01 | 2022-08-30 | 南京隼眼电子科技有限公司 | Mountain landslide prediction method based on millimeter wave radar and electronic equipment |
CN114966601B (en) * | 2022-08-01 | 2022-10-21 | 南京隼眼电子科技有限公司 | Mountain landslide prediction method based on millimeter wave radar and electronic equipment |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111522005A (en) | Deformation monitoring and terrain reconstruction method | |
GREJNER‐BRZEZINSKA | Direct exterior orientation of airborne imagery with GPS/INS system: Performance analysis | |
CN108868268B (en) | Unmanned parking space posture estimation method based on point-to-surface distance and cross-correlation entropy registration | |
WO2022214114A2 (en) | Bridge deformation monitoring method fusing gnss data and insar technology | |
US9927513B2 (en) | Method for determining the geographic coordinates of pixels in SAR images | |
CN110108984B (en) | Spatial relationship synchronization method for multiple sensors of power line patrol laser radar system | |
CN103323855B (en) | A kind of precision acquisition methods of baseline dynamic measurement system | |
CN113359097B (en) | Millimeter wave radar and camera combined calibration method | |
Mostafa et al. | A multi-sensor system for airborne image capture and georeferencing | |
JP2590689B2 (en) | Interferometric synthetic aperture radar system and terrain change observation method | |
WO2020151213A1 (en) | Air and ground combined intertidal zone integrated mapping method | |
CN110068817B (en) | Terrain mapping method, instrument and system based on laser ranging and InSAR | |
CN112882030B (en) | InSAR imaging interference integrated processing method | |
JP2596364B2 (en) | Topographic map generator using three-dimensional information obtained from interferometric synthetic aperture radar | |
CN110986888A (en) | Aerial photography integrated method | |
CN112130151B (en) | Arc synthetic aperture ground radar coordinate projection rapid calculation method | |
Behan et al. | Steps towards quality improvement of airborne laser scanner data | |
Holecz et al. | Height model generation, automatic geocoding and a mosaicing using airborne AeS-1 InSAR data | |
CN113483739B (en) | Offshore target position measuring method | |
CN115951369A (en) | Multi-sensor fusion positioning method for complex port environment | |
CN116202410A (en) | Geological disaster monitoring method and device, electronic equipment and storage medium | |
CN106371096B (en) | Airborne double-antenna InSAR three-dimensional configuration model construction method | |
Zhang | Photogrammetric processing of low altitude image sequences by unmanned airship | |
Kartal et al. | Comperative analysis of different geometric correction methods for very high resolution pleiades images | |
CN111856464B (en) | DEM extraction method of vehicle-mounted SAR (synthetic aperture radar) based on single control point information |
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: 20200811 |
|
RJ01 | Rejection of invention patent application after publication |