Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular system configurations, techniques, etc. in order to provide a thorough understanding of the disclosed embodiments. However, it will be apparent to one skilled in the art that the present disclosure may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present disclosure with unnecessary detail.
For the purposes of promoting an understanding of the principles and advantages of the disclosure, reference will now be made to the embodiments illustrated in the drawings.
Referring to fig. 1, fig. 1 is a flow chart of a three-dimensional deformation resolving method of a ground-based radar according to an embodiment of the disclosure, where the method includes:
s101, multi-angle observation is carried out on deformation information of the same observation target based on a plurality of radars, and multi-angle one-dimensional vision line deformation data of the same observation target are obtained.
In this embodiment, deformation information of a certain object is observed, and if the observed object is a landslide, one-dimensional line of sight deformation data of the landslide is calculated.
The deformation information is a description of the relevant conditions of the object (namely, the observation object) which is reflected by the change of the shape, the size and the like of the object under the action of various internal and external forces.
The one-dimensional line-of-sight deformation data is data which only reflects the deformation condition of the target object in a single dimension direction. In practical applications, the line of sight is often used as a single dimension to acquire data, i.e. only the displacement or deformation amount information of the target object in a specific straight line direction is concerned, for example, the straight line direction determined by the radar wave transmitting direction and the receiving reflected wave direction when the radar is observed.
Therefore, the radar acquires the deformation information of the target object in the sight line direction, and one-dimensional sight line deformation data of the target object can be obtained. In addition, the embodiment can acquire multi-angle one-dimensional visual line deformation data from multiple directions by utilizing multiple radars.
S102, resolving the multi-angle one-dimensional visual line deformation data based on a parameterized least square method to obtain target three-dimensional deformation data of the same observation target, wherein a resolving weight matrix of the parameterized least square method is determined by a variance component estimation method.
In the present embodiment, parameterized least squares is a method of estimating model parameters by minimizing the sum of squares of residuals between observations and model predictions.
The solution is a process of in turn deriving unknown parameters or physical quantities in the model from known observation data and mathematical models. According to the calculated one-dimensional visual line deformation data, the three-dimensional deformation data of the object to be observed are calculated through a parameterized least square method.
The solution weight matrix is a matrix used for measuring the importance of different observation data or model parameters in the parameterized least square solution process. The larger the weight of each element in the solution weight matrix corresponding to the corresponding observed data or parameter, the greater the influence of the data or parameter in the fitting process is, and the more important the contribution to the final result is. Therefore, the embodiment can make the parameterized least square method pay more attention to certain key data or parameters by reasonably setting the resolving weight matrix, thereby improving the resolving accuracy and reliability.
The variance component estimation method is a statistical method for estimating variance components of different observed data sources or different error sources. Under the condition that multiple kinds of observation data participate, different types of data can have different precision and reliability, and the variance component estimation method can analyze and determine the variance of each group of data according to the actual condition of the observation data, so that weights are reasonably distributed according to the variances, and the overall data processing and analyzing precision is improved.
Specifically, the steps of the present embodiment may be as follows:
firstly, taking the calculated multi-angle one-dimensional line of sight as input to deform data, and adopting a parameterized least square method to calculate;
secondly, calculating initial three-dimensional deformation data;
Then, in the parameterized least square method solution process, one key step is to determine a solution weight matrix;
Finally, the determination of the solution weight matrix can be obtained through a variance component estimation method, and the determination of the solution weight matrix is beneficial to the subsequent obtaining of target three-dimensional deformation data which are more accurate than the initial three-dimensional deformation data.
In this embodiment, a solution weight matrix is determined according to the variance component estimation method, and the parameterized least square method is updated through the solution weight matrix, so that target three-dimensional deformation data of the observation target can be obtained. Therefore, the initial three-dimensional deformation data is also updated. Compared with the initial three-dimensional deformation data, the accuracy and the reliability of the target three-dimensional deformation data and the degree of fit between the target three-dimensional deformation data and the actual observed object three-dimensional deformation data are improved.
From the above, the method and the device adopt the parameterized least square method and determine the resolving weight matrix by means of the variance component estimation method to carry out three-dimensional deformation resolving, can comprehensively consider the influence of various factors on different data, reasonably allocate weights, enable the resolving process to be more accurate, improve the reliability of target three-dimensional deformation data, can be accurately used in the fields of infrastructure health detection monitoring, geological disaster early warning and the like, and provide powerful and accurate data support for related decisions. Therefore, the three-dimensional deformation resolving precision can be improved, and the monitoring precision and the early warning capability of geological disasters are improved.
In one embodiment of the present disclosure, performing multi-angle observation on deformation information of a same observation target based on a plurality of radars, obtaining multi-angle one-dimensional line of sight deformation data of the same observation target, includes:
Measuring the deformation of the real deformation of the observation target point projected to the sight line direction, and obtaining multidirectional sight line deformation data, wherein the multidirectional sight line deformation data are multi-angle one-dimensional sight line deformation data;
wherein, the radar has at least three parts.
In this embodiment, the radar can emit electromagnetic wave signals and receive signals reflected by the target object, so as to monitor deformation of the target object, for example, in a landslide deformation monitoring scene, the radar can be used to observe displacement of a target point on a landslide.
The line of sight deformation data is calculated based on the observation principle of the radar, and the displacement, deformation and other deformation conditions of the target object are generated by the target point in the straight line direction (namely, the sight direction) determined by the radar wave transmitting direction and the receiving reflected wave direction. Meanwhile, the sight line deformation data are also one-dimensional sight line deformation data.
The radar according to this embodiment has three parts, and can obtain deformation data of one-dimensional vision lines in three directions, see fig. 2, where fig. 2 (a) is a radar 1 vision line direction observation result, fig. 2 (b) is a radar 2 vision line direction observation result, and fig. 2 (c) is a radar 3 vision line direction observation result.
According to the method, multiple radars can be accurately monitored from a specific direction, basic data support is provided for subsequent three-dimensional deformation calculation, and rapid development of deformation analysis of an observation target is facilitated.
In one embodiment of the present disclosure, resolving multi-angle one-dimensional gaze direction deformation data based on a parameterized least square method to obtain target three-dimensional deformation data of the same observation target, including:
determining a least square method observation equation;
determining fusion space parameters;
Determining a parameterized least square method observation equation based on the least square method observation equation and the fusion space parameter, wherein the parameterized least square method observation equation is a parameterized least square method;
And inputting the multi-angle one-dimensional view line into the deformation data to determine the target three-dimensional deformation data of the observation target by using a parameterized least square method.
In the present embodiment, the least square observation equation is:
Wherein, Representation radarThe one-dimensional line-of-sight deformation data of the target object observed at each point in time,,The matrix of the unit direction is represented,Representing the three-dimensional deformation data that is sought,The residual is represented by a representation of the residual,Representing the deformation data of the one-dimensional line of sight observed by the radar 1,Representing the deformation data of the one-dimensional line of sight observed by the radar 2,Representing the deformation data of the one-dimensional line of sight observed by the radar 3.
The fusion space parameter is a parameter related to a space characteristic, and the space characteristic relates to different directions (such as east-west, north-south, vertical directions and the like) and relations among different positions and the like under a scene of three-dimensional deformation of a target object.
Specifically, determining fusion space parameters) Comprising:
Defining a row vector As a fused spatial parameterWhere k represents the number of target points, k= {1,2,..,Represents the coordinate value of the target point, G is defined as follows,
Fusion of spatial parametersA design matrix is also provided.
The parameterized least squares observation equation is:
Wherein, Representation radarFor one-dimensional line of sight deformation data of the observation target point,The matrix of coefficients is represented and,Representing the design matrix.
Expressed as basis functionsCan be expressed as a linear superposition of:
Wherein, Representing a three-dimensional deformation vector, i.e. three-dimensional deformation data,The design matrix is represented by a representation of the design matrix,The parameter vector is represented by a vector of parameters,The values of the parameters are represented by the values,Representing the degree of the introduced polynomial,Representation ofIs the first of (2)Column vectors.
Is a block matrix composed ofThe composition of the composite material comprises the components,Representing target scatterers to each radarIs used for the matrix of unit vectors of (a),AndThe following is shown:
The purpose of the parameter least square method is to obtain To the point ofA kind of electronic device:
The multi-angle one-dimensional visual line direction deformation data is used as the input of the parameterized least square method, and the target three-dimensional deformation data can be determined.
The target three-dimensional deformation data is a solution result of the present disclosure, see fig. 3, where fig. 3 (a) is a solution result of a parameterized least square method X direction based on variance component estimation, fig. 3 (b) is a solution result of a parameterized least square method Y direction based on variance component estimation, and fig. 3 (c) is a solution result of a parameterized least square method Z direction based on variance component estimation.
According to the method, the fusion space parameters are tightly combined with the three-dimensional space characteristics, so that the constructed parameterized least square observation equation and the corresponding method can effectively process the calculation from deformation data of a plurality of radar one-dimensional vision lines to three-dimensional deformation data of a target, and reliable initial data support is provided for follow-up accurate analysis of the overall deformation trend of the observed target, evaluation of stability, prediction of geological disaster risk and the like.
In one embodiment of the present disclosure, a method for three-dimensional deformation resolution of a ground-based radar further includes:
Determining an initial variance based on a parameterized least squares method;
an initial weight matrix is determined based on the initial variance.
In this embodiment, a weighting matrix is introduced,Is thatIs used for the matrix of the matrix,Is a diagonal matrix.
According to the known conditions, byCan calculate the parameter vector and then calculateCan calculate three-dimensional deformation result and weighting matrixThe initial weight matrix is represented, and is composed of initial variances.
From the above, the present embodiment determines the initial variance and quantifiable data discrete degree, and knows the fluctuation characteristics of the data. And an initial weight matrix is constructed based on the initial variance, so that different weights can be given according to the stability and reliability differences of the data. Therefore, in the resolving process, more stable and reliable data are strengthened, the influence of unstable data is weakened, and the resolving accuracy and reliability are effectively improved.
In one embodiment of the present disclosure, the solution weight matrix of the parameterized least squares method is determined by a variance component estimation method, comprising:
iterating the initial variance based on the variance component estimation method;
if the iteration is carried out until the difference value between the two adjacent variances is greater than or equal to a first threshold value, continuing the iteration;
if the iteration is carried out until the difference value between the two adjacent variances is smaller than a first threshold value, stopping the iteration;
After iteration is completed, obtaining a target variance;
a solution weight matrix is determined based on the target variance.
In the present embodiment of the present invention, in the present embodiment,Representing a covariance matrix, defined as follows:
Wherein, Representing an initial standard deviation, obtained by a parameterized least squares method,The initial weight matrix representing the radar observation of section b, typically consists of an identity matrix,Representing an identity matrix of size,Indicating the observation time value.
Residual vectorThe calculation of (a) depends on the residual projection matrix,
The residual vector is defined as the sum of the values,
Construction matrixSum vector:
Wherein, each element is defined as follows:
Wherein, The quadratic vector representing the observed value correction is updated for variance by:
Iterating the process until The difference is less than a thresholdThe selection of the threshold value can be determined according to the empirical variance level of various observed values. With the algorithm iterated, the updated variance components are used to recalculate the weight matrix and again to calculate the regression coefficients to ensure that the regression model is gradually optimized. Obtaining a new weight matrix by using the updated variance valueI.e. a solution weight matrix for the subsequent three-dimensional deformation solution.
From the above, the present embodiment updates the weight matrix by the variance component estimation method, so as to solve the problem of unreasonable weight matrix setting. And then, the determined resolving weight matrix can reasonably distribute data weights in the parameterized least square method, so that the accuracy of the subsequent three-dimensional deformation resolving is improved, and the geological research and disaster early warning are facilitated.
In one embodiment of the present disclosure, a method for three-dimensional deformation resolution of a ground-based radar further includes:
Obtaining a target parameter vector based on the solution weight matrix;
and updating the initial three-dimensional deformation data based on the target parameter vector, and determining the target three-dimensional deformation data.
In this embodiment, first, each weight value in the solution weight matrix affects different data elements or parameters to different extents, so as to guide the generation of the target parameter vector.
And then, updating the initial three-dimensional deformation data by using the obtained target parameter vector. And correcting and perfecting deformation values (such as displacement amounts, strain values and the like in different directions) of each dimension in the initial three-dimensional deformation data. Because the target parameter vector is generated by integrating more accurate data weight relation reflected by the target weight matrix, the method can provide more reasonable and more practical basis for updating the initial three-dimensional deformation data.
And finally, determining target three-dimensional deformation data which is more accurate and can more truly reflect the actual three-dimensional deformation condition of the observed target object. The target three-dimensional deformation data can be used for subsequent important application scenes such as infrastructure health detection, geological disaster early warning and the like, and reliable data support is provided for related decisions.
In one embodiment of the present disclosure, a method for three-dimensional deformation resolution of a ground-based radar further includes:
Determining real three-dimensional deformation data of an observation target;
and calculating based on the real three-dimensional deformation data and the target three-dimensional deformation data, and evaluating the accuracy of the calculation result.
In this embodiment, the real three-dimensional deformation data is used as standard or reference three-dimensional deformation data for evaluating the resolution accuracy of the target three-dimensional deformation data.
The real three-dimensional deformation data of the target observation point is shown in fig. 4, wherein fig. 4 (a) is the real deformation in the X direction, fig. 4 (b) is the real deformation in the Y direction, and fig. 4 (c) is the real deformation in the Z direction.
The three-dimensional deformation data of the target can be obtained by resolving in steps S101-S102.
The present embodiment is used to evaluate the accuracy of the solution by calculating the root mean square error (Root Mean Squared Error, RMSE).
Wherein, Representing the actual three-dimensional deformation data,Representing three-dimensional deformation data of the object,Indicating the observation time value.
It can be obtained from the above that, in the embodiment, the real three-dimensional deformation data and the target three-dimensional deformation data obtained by resolving are compared and calculated, so that the resolving precision can be evaluated, the accuracy of the resolving result can be intuitively checked, a reliable basis can be provided for geological disaster early warning, the relevant departments can accurately decide, and the disaster loss is reduced.
Referring to fig. 6, fig. 6 is a schematic diagram of a specific flow of a three-dimensional deformation resolving method of a ground-based radar.
Specifically, the step of parameterizing the least squares method is as follows:
firstly, calculating one-dimensional deformation data of multiple radar vision lines;
Secondly, establishing a least square method model;
Third, spatial parameters are introduced by popularization Sum parameter vector;
Fourth, three-dimensional displacement vectorFrom the following componentsLinear superposition;
fifth, design matrix Defining the number k of observation points and polynomial degree,Is a unit vector matrix;
sixth, calculate the initial variance and initial weight matrix
Seventh, calculate the parameter vector;
Eighth, solve to obtain three-dimensional displacement vector。
Specifically, the variance component estimation method includes the steps of:
First, constructing a covariance matrix: ;
Second, calculate the residual vector MatrixVector quantity;
Third, update variance component
Fourth, if it does not meetContinuing to execute the third step;
Fifth, if it meets The variance and weight matrix are updated.
In this embodiment, the initial variance and the initial weight matrix are updated, that is, the updated variance (target variance) and the weight matrix (solution weight matrix) may be obtained, and a new parameter vector may be obtained, so as to finally obtain the target three-dimensional deformation data.
According to the method, the target parameter vector is obtained based on the resolving weight matrix, the weight can be reasonably set, the target parameter vector can be more accurate, the final three-dimensional deformation resolving precision can be improved, and the method is beneficial to geological disaster prediction and prevention.
Fig. 6 is a block diagram of a three-dimensional deformation resolving system of a ground-based radar according to an embodiment of the present disclosure. For ease of illustration, only portions relevant to embodiments of the present disclosure are shown. Referring to fig. 6, the three-dimensional deformation resolving system 20 of the ground-based radar includes a data acquisition module 21 and a resolving module 22.
The data acquisition module 21 performs multi-angle observation on deformation information of the same observation target based on multiple radars to obtain multi-angle one-dimensional vision line deformation data of the same observation target;
The resolving module 22 performs comprehensive resolving on the multi-angle one-dimensional visual line deformation data based on a parameterized least square method to obtain target three-dimensional deformation data of the same observation target, wherein a resolving weight matrix of the parameterized least square method is determined by a variance component estimation method.
In one embodiment of the present disclosure, the data acquisition module 21 is specifically further configured to:
Measuring the deformation of the real deformation of the observation target point projected to the sight line direction, and obtaining multidirectional sight line deformation data, wherein the multidirectional sight line deformation data are multi-angle one-dimensional sight line deformation data;
wherein, the radar has at least three parts.
In one embodiment of the present disclosure, the resolving module 22 is specifically further configured to:
determining a least square method observation equation;
determining fusion space parameters;
Determining a parameterized least square method observation equation based on the least square method observation equation and the fusion space parameter, wherein the parameterized least square method observation equation is a parameterized least square method;
And inputting the multi-angle one-dimensional view line into the deformation data to determine target three-dimensional deformation data of the same observation target by using a parameterized least square method.
In one embodiment of the present disclosure, a ground-based radar three-dimensional deformation solution system 20 further comprises an initial weight determination module;
An initial weight determination module for determining an initial variance based on a parameterized least squares method;
an initial weight matrix is determined based on the initial variance.
In one embodiment of the present disclosure, the resolving module 22 is specifically further configured to:
iterating the initial variance based on the variance component estimation method;
if the iteration is carried out until the difference value between the two adjacent variances is greater than or equal to a first threshold value, continuing the iteration;
if the iteration is carried out until the difference value between the two adjacent variances is smaller than a first threshold value, stopping the iteration;
After iteration is completed, obtaining a target variance;
a solution weight matrix is determined based on the target variance.
In one embodiment of the present disclosure, the resolving module 23 is specifically further configured to:
Obtaining a target parameter vector based on the solution weight matrix;
and updating the initial three-dimensional deformation data based on the target parameter vector, and determining the target three-dimensional deformation data.
In one embodiment of the present disclosure, a ground-based radar three-dimensional deformation solution system 20 further comprises a solution result evaluation module;
the resolving result evaluation module is used for determining the real three-dimensional deformation data of the observation target;
and calculating based on the real three-dimensional deformation data and the target three-dimensional deformation data, and evaluating the accuracy of the calculation result.
Referring to fig. 7, fig. 7 is a schematic block diagram of an electronic device according to an embodiment of the disclosure. The electronic device 300 in this embodiment as shown in fig. 7 may include one or more processors 301, one or more input devices 302, one or more output devices 303, and one or more memories 304. The processor 301, the input device 302, the output device 303, and the memory 304 communicate with each other via a communication bus 305. The memory 304 is used to store a computer program comprising program instructions. The processor 301 is configured to execute program instructions stored in the memory 304. Wherein the processor 301 is configured to invoke program instructions to perform the functions of the modules/units of the system embodiments described above, such as the functions of the modules 21 to 22 shown in fig. 6.
It should be appreciated that in the disclosed embodiments, the Processor 301 may be a central processing unit (Central Processing Unit, CPU), which may also be other general purpose processors, digital signal processors (DIGITAL SIGNAL processors, DSPs), application SPECIFIC INTEGRATED Circuits (ASICs), field-Programmable gate arrays (Field-Programmable GATE ARRAY, FPGA) or other Programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The input device 302 may include a touch pad, a fingerprint collection sensor (for collecting fingerprint information of a user and direction information of a fingerprint), a microphone, etc., and the output device 303 may include a display (LCD, etc.), a speaker, etc.
The memory 304 may include read only memory and random access memory and provides instructions and data to the processor 301. A portion of memory 304 may also include non-volatile random access memory. For example, the memory 304 may also store information of device type.
In a specific implementation, the processor 301, the input device 302, and the output device 303 described in the embodiments of the present disclosure may perform the implementation manners described in the first embodiment and the second embodiment of the three-dimensional deformation resolving method for a ground-based radar provided in the embodiments of the present disclosure, and may also perform the implementation manners of the electronic device described in the embodiments of the present disclosure, which are not described herein again.
In another embodiment of the disclosure, a computer readable storage medium is provided, where the computer readable storage medium stores a computer program, where the computer program includes program instructions, where the program instructions, when executed by a processor, implement all or part of the procedures in the method embodiments described above, or may be implemented by instructing related hardware by the computer program, where the computer program may be stored in a computer readable storage medium, where the computer program, when executed by the processor, implements the steps of each of the method embodiments described above. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, executable files or in some intermediate form, etc. The computer readable medium may include any entity or device capable of carrying computer program code, recording medium, USB flash disk, removable hard disk, magnetic disk, optical disk, computer Memory, read-Only Memory (ROM), random access Memory (RAM, random Access Memory), electrical carrier signals, telecommunications signals, and software distribution media, among others.
The computer readable storage medium may be an internal storage unit of the electronic device of any of the foregoing embodiments, such as a hard disk or a memory of the electronic device. The computer readable storage medium may also be an external storage device of the electronic device, such as a plug-in hard disk provided on the electronic device, a smart memory card (SMART MEDIA CARD, SMC), a Secure Digital (SD) card, a flash memory card (FLASH CARD), or the like. Further, the computer-readable storage medium may also include both internal storage units and external storage devices of the electronic device. The computer-readable storage medium is used to store a computer program and other programs and data required for the electronic device. The computer-readable storage medium may also be used to temporarily store data that has been output or is to be output.
Those of ordinary skill in the art will appreciate that the elements and algorithm steps described in connection with the embodiments disclosed herein may be embodied in electronic hardware, in computer software, or in a combination of the two, and that the elements and steps of the examples have been generally described in terms of function in the foregoing description to clearly illustrate the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present disclosure.
It will be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working process of the electronic device and unit described above may refer to the corresponding process in the foregoing method embodiment, which is not repeated herein.
In the several embodiments provided in the present application, it should be understood that the disclosed electronic device and method may be implemented in other manners. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of elements is merely a logical functional division, and there may be additional divisions of actual implementation, e.g., multiple elements or components may be combined or integrated into another system, or some features may be omitted, or not performed. In addition, the coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection via some interfaces or units, or may be an electrical, mechanical, or other form of connection.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed over a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purposes of the embodiments of the present disclosure.
In addition, each functional unit in each embodiment of the present disclosure may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The foregoing is merely a specific embodiment of the present disclosure, but the protection scope of the present disclosure is not limited thereto, and any equivalent modifications or substitutions will be apparent to those skilled in the art within the scope of the present disclosure, and these modifications or substitutions should be covered in the scope of the present disclosure. Therefore, the protection scope of the present disclosure shall be subject to the protection scope of the claims.