CN113866768A - Time sequence interference radar phase optimization estimation method - Google Patents
Time sequence interference radar phase optimization estimation method Download PDFInfo
- Publication number
- CN113866768A CN113866768A CN202111455165.4A CN202111455165A CN113866768A CN 113866768 A CN113866768 A CN 113866768A CN 202111455165 A CN202111455165 A CN 202111455165A CN 113866768 A CN113866768 A CN 113866768A
- Authority
- CN
- China
- Prior art keywords
- phase
- likelihood function
- estimation
- value
- observation vector
- 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.)
- Granted
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
- 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
- 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
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 discloses a time sequence interference radar phase optimization estimation method, which comprises the following steps: acquiring an observation vector transmitted by a time sequence interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar; obtaining a phase likelihood function corresponding to the observation vector according to the observation vector; and carrying out maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimation phase. According to the embodiment of the invention, the phase quality of the scatterer in the time sequence interference radar can be improved by carrying out maximum likelihood estimation and optimization processing on the phase of the observation vector transmitted by the time sequence interference radar.
Description
Technical Field
The invention relates to the technical field of interference radars, in particular to a time sequence interference radar phase optimization estimation method.
Background
Due to the rapid development of satellite radar technology in recent years, the synthetic aperture radar interferometry (InSAR) technology has been moved from the research and development stage to practical engineering applications. The technique uses multiple synthetic aperture radar images to generate a map of the earth's surface deformation from the phase differences of echoes received by a satellite or aircraft. The technology has the characteristics of all weather and all day long, can penetrate through cloud layer dense fog smoke dust and can obtain ground surface deformation information in a large area, and is particularly suitable for areas with difficult imaging of the traditional optical sensor, and the monitoring precision can reach millimeter level. Based on these advantages, the InSAR technology has played a great role in the fields of emergency and daily monitoring.
In sequential synthetic aperture radar interferometric InSAR processing, sequential interferometric radars typically involve two types of point targets. The prior art has poor effect of optimizing and estimating the phase of a part of point targets in the time sequence interference radar.
Thus, there is still a need for improvement and development of the prior art.
Disclosure of Invention
The technical problem to be solved by the present invention is to provide a method for optimizing and estimating a phase of a time-series interference radar, aiming at solving the problem of poor effect of optimizing and estimating the phase of a part of point targets in the time-series interference radar in the prior art.
The technical scheme adopted by the invention for solving the problems is as follows:
in a first aspect, an embodiment of the present invention provides a method for estimating phase optimization of a time-series interference radar, where the method includes:
acquiring an observation vector transmitted by a time sequence interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar;
obtaining a phase likelihood function corresponding to the observation vector according to the observation vector;
and carrying out maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimation phase.
In one implementation, the obtaining, according to the observation vector, a phase likelihood function corresponding to the observation vector includes:
acquiring a real number covariance matrix; the real covariance matrix is used for representing a scattering point real covariance value transmitted by the time sequence interference radar;
performing conditional probability operation on the observation vector and the real number covariance matrix to obtain a conditional probability function;
and carrying out logarithmic operation on the conditional probability function to obtain a phase likelihood function corresponding to the observation vector.
In one implementation, the performing maximum likelihood estimation and optimization on the phase likelihood function to obtain a target estimated phase includes:
carrying out derivation operation on the phase likelihood function to obtain a first derivative function;
solving the first derivative function after zero setting to obtain a maximum likelihood estimation value corresponding to the real covariance matrix;
and optimizing the phase likelihood function based on the maximum likelihood estimation value to obtain a target estimation phase.
In one implementation, the optimizing the phase likelihood function based on the maximum likelihood estimation value to obtain the target estimated phase includes:
substituting the maximum likelihood estimation value into the phase likelihood function to obtain an updated phase likelihood function;
and optimizing the updated phase likelihood function to obtain the target estimated phase.
In one implementation, the updated phase likelihood function includes a first sub-term, a second sub-term, and a third sub-term; the optimizing the updated phase likelihood function to obtain the target estimated phase includes:
calculating the minimum value of the second sub-item of the updated phase likelihood function based on a preset optimization algorithm;
and substituting the minimum value into the updated phase likelihood function to obtain the target estimated phase.
In one implementation, the optimizing the updated phase likelihood function to obtain the target estimated phase further includes:
acquiring a preset sample complex covariance matrix;
solving the absolute value of the sample complex covariance matrix to obtain the modulus of the sample complex covariance matrix;
substituting the modulus of the sample complex covariance matrix as a real covariance matrix into a preset first formula for calculation to obtain a first approximate phase;
and carrying out iterative computation on the first approximate phase based on the updated phase likelihood function and a preset first formula to obtain a target estimated phase.
In one implementation, the iteratively calculating the first approximate phase based on the updated phase likelihood function and a preset first formula to obtain a target estimated phase includes:
substituting the first approximate phase into a preset second formula to obtain an approximate maximum likelihood estimation value;
substituting the approximate maximum likelihood estimation value serving as a real covariance matrix into a preset first formula for calculation to obtain a second approximate phase;
substituting the second approximate phase into the second sub-term to obtain a second sub-term value;
and repeatedly executing the step of substituting the first approximate phase into a preset second formula to obtain an approximate maximum likelihood estimation value, stopping iteration when the second subentry value is smaller than a preset threshold value, and substituting the second subentry value when the iteration is stopped into an updated phase likelihood function to obtain a target estimated phase.
In one implementation, the optimizing the updated phase likelihood function to obtain the target estimated phase further includes:
carrying out derivation operation on the updated phase likelihood function to obtain a second derivative function;
and solving the second derivative function after the zero setting to estimate the phase by the target.
In a second aspect, an embodiment of the present invention further provides a time-series interferometric radar phase optimization estimation apparatus, where the apparatus includes:
the observation vector acquisition module is used for acquiring an observation vector transmitted by the time sequence interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar;
the phase likelihood function determining module is used for obtaining a phase likelihood function corresponding to the observation vector according to the observation vector;
and the target estimation phase determining module is used for carrying out maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimation phase.
In a third aspect, an embodiment of the present invention further provides an intelligent terminal, including a memory, and one or more programs, where the one or more programs are stored in the memory, and configured to be executed by one or more processors includes a processor configured to execute the method for estimating phase optimization of time-series interferometric radar according to any one of the above items.
In a fourth aspect, embodiments of the present invention further provide a non-transitory computer-readable storage medium, where instructions of the storage medium, when executed by a processor of an electronic device, enable the electronic device to perform the time-series interference radar phase optimization estimation method according to any one of the above.
The invention has the beneficial effects that: firstly, acquiring an observation vector transmitted by a time sequence interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar; then according to the observation vector, obtaining a phase likelihood function corresponding to the observation vector; finally, carrying out maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimation phase; therefore, the phase quality of the point target in the time sequence interference radar can be improved by carrying out maximum likelihood estimation and optimization processing on the phase of the observation vector transmitted by the time sequence interference radar in the embodiment of the invention.
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, it is obvious that the drawings in the following description are only some embodiments described in the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for phase optimization estimation of a time-series interference radar according to an embodiment of the present invention.
Fig. 2 is a schematic block diagram of a time-series interferometric radar phase optimization estimation apparatus according to an embodiment of the present invention.
Fig. 3 is a schematic block diagram of an internal structure of an intelligent terminal according to an embodiment of the present invention.
Detailed Description
The invention discloses a time sequence interference radar phase optimization estimation method, which is further described in detail below by referring to the attached drawings and embodiments in order to make the purpose, technical scheme and effect of the invention clearer and clearer. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
As used herein, the singular forms "a", "an", "the" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. It will be understood that when an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may also be present. Further, "connected" or "coupled" as used herein may include wirelessly connected or wirelessly coupled. As used herein, the term "and/or" includes all or any element and all combinations of one or more of the associated listed items.
It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
Since in the prior art, two types of point targets are typically involved in the time series InSAR process. One is a Permanent Scatterer (PS) that can have a stable phase, and the other is a Distributed Scatterer (DS) that has a relatively poor phase quality. The PS is easy to process due to good phase quality, and the time sequence InSAR algorithm for processing the PS point is abundant. The DS point has poor phase quality and great processing challenge, and cannot be processed together with the PS point.
In order to solve the problems in the prior art, the embodiment provides a phase optimization estimation method for a time-series interference radar, which can improve the phase quality of scatterers in the time-series interference radar by performing maximum likelihood estimation and optimization processing on the phase of an observation vector transmitted by the time-series interference radar. In specific implementation, firstly, acquiring an observation vector transmitted by a time sequence interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar; then according to the observation vector, obtaining a phase likelihood function corresponding to the observation vector; and finally, carrying out maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimation phase. Exemplary method
The embodiment provides a phase optimization estimation method for a time-series interferometric radar, which can be applied to an intelligent terminal of the interferometric radar. As shown in fig. 1 in detail, the method includes:
s100, acquiring an observation vector transmitted by a time sequence interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar;
in particular, observation vectorsIs oneThe vector value of the dimension, wherein,the data are transmitted by the time-series interference radar and mainly comprise data of scattering points in a time dimension. Phase information can be extracted from the observation vector, so that the obtained observation vector is prepared for obtaining the target estimated phase subsequently.
After the observation vector is obtained, the following steps can be performed as shown in fig. 1: s200, obtaining a phase likelihood function corresponding to the observation vector according to the observation vector;
correspondingly, the step of obtaining the phase likelihood function corresponding to the observation vector according to the observation vector includes the following steps:
s201, acquiring a real number covariance matrix; the real covariance matrix is used for representing a scattering point real covariance value transmitted by the time sequence interference radar;
s202, performing conditional probability operation on the observation vector and the real covariance matrix to obtain a conditional probability function;
s203, carrying out logarithm operation on the conditional probability function to obtain a phase likelihood function corresponding to the observation vector.
In particular, a real covariance matrixThe method is characterized in that a real covariance matrix of scattering points transmitted by a time sequence interference radar gives phases of the scattering points under the assumption of complex circle Gaussian distribution(in this embodiment, a closed phase) and a real covariance matrix G, performing conditional probability operation on the observation vector and the real covariance matrix to obtain a conditional probability function; the formula of the conditional probability function is as follows:
wherein the content of the first and second substances,in the form of a matrix of a closed vector of the scattering points in time dimension,(ii) a Complex covariance matrixThe relationship to the closed-phase and real covariance matrices is as follows:
wherein the content of the first and second substances,is a scattering point complex covariance matrix.
The method for calculating the target estimated phase in the prior art is as follows: suppose a real covariance matrixKnown as closed phaseThe maximum likelihood estimation of (c) can be achieved by maximizing the following expression:
wherein the content of the first and second substances,is a time-dimensional closed phase vector of the scattering points,an estimate of a random vector is characterized.Is a sample complex covariance matrix given by the following equation:
wherein the content of the first and second substances,is a set of adjacent homogenous points and is,number of homogenous spots.
Due to real covariance matrixNot in fact known, existing methods use the modulus of a sample complex covariance matrixTo replace the real covariance matrixThe phase optimization is performed according to equation (3). The significant drawback of the prior art method is that if the sample complex covariance matrix deviates too far from the true value, the estimation effect is greatly affected and the optimization effect is poor.
The invention adopts a new phase estimation method to avoid using a sample complex covariance matrixTo replace the real covariance matrixAnd realizing accurate maximum likelihood estimation. Therefore, after obtaining the conditional probability function, performing a logarithm operation on the conditional probability function to obtain a phase likelihood function corresponding to the observation vector, where the phase likelihood function is as follows:
after obtaining the phase likelihood function, the following steps can be performed as shown in fig. 1: s300, carrying out maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimation phase. Correspondingly, the maximum likelihood estimation and optimization processing of the phase likelihood function to obtain the target estimated phase includes the following steps:
s301, carrying out derivation operation on the phase likelihood function to obtain a first derivative function;
s302, solving the first derivative function after zero setting to obtain a maximum likelihood estimation value corresponding to the real number covariance matrix;
and S303, optimizing the phase likelihood function based on the maximum likelihood estimation value to obtain a target estimation phase.
converting the likelihood function to the following formula:
wherein tr is a matrix trace, a derivation operation is performed on the phase likelihood function of the formula (7) to obtain a first derivative function, and then the first derivative function after being zeroed is solved to obtain a maximum likelihood estimation value corresponding to the real number covariance matrix; in the present embodiment, the pair of likelihood functionsTaking the derivative and let it be 0, the following equation is obtained:
solving the above equation can obtain the maximum likelihood estimation value corresponding to the real covariance matrix:
And then, optimizing the phase likelihood function based on the maximum likelihood estimation value to obtain a target estimation phase. Correspondingly, the optimizing the phase likelihood function based on the maximum likelihood estimation value to obtain the target estimation phase includes the following steps: substituting the maximum likelihood estimation value into the phase likelihood function to obtain an updated phase likelihood function; and optimizing the updated phase likelihood function to obtain the target estimated phase.
Specifically, the maximum likelihood estimation value is substituted into the phase likelihood function, and the order is givenObtaining an updated phase likelihood function:
wherein p is a p-dimensional identity matrix. The formula (10) is an innovative key point of the invention, because a phase optimization estimation mathematical formula of the distributed scattering points is established, an accurate maximum likelihood estimation value is obtained, and the assumed conditions in the existing method are removed, thereby overcoming the influence on the phase estimation result caused by too large deviation between the assumption and the actual in the existing method. And after the phase likelihood function is obtained, optimizing the updated phase likelihood function to obtain the target estimation phase.
In one implementation, the updated phase likelihood function includes a first sub-term, a second sub-term, and a third sub-term; the optimizing the updated phase likelihood function to obtain the target estimated phase comprises the following steps: calculating the minimum value of the second sub-item of the updated phase likelihood function based on a preset optimization algorithm; and substituting the minimum value into the updated phase likelihood function to obtain the target estimated phase.
Specifically, after the updated phase likelihood function is maximized, the target estimated phase (i.e. the maximum likelihood estimation value of the closed phase) can be obtained, and the first sub-term of the phase likelihood function isSecond sonThe item isThe third sub-item isTo maximize the updated phase likelihood function, a second sub-term is soughtIn this embodiment, the preset optimization algorithm is a Newton-Raphson method, a simulated annealing method, and a search for the minimum value ofIs a minimum value ofSubstituting the minimum value of the target phase into the updated phase likelihood function to obtain the target estimated phase. In practice, because the minimum value falls into a local optimal solution, the existing method and other phase optimization methods can be used to obtain a plurality of minimum values, the minimum values are iterated, the minimum value during convergence is taken as a final minimum value, and the final minimum value is substituted into the updated phase likelihood function to obtain the target estimated phase.
In one implementation, the optimizing the updated phase likelihood function to obtain the target estimated phase further includes: acquiring a preset sample complex covariance matrix; solving the absolute value of the sample complex covariance matrix to obtain the modulus of the sample complex covariance matrix; substituting the modulus of the sample complex covariance matrix as a real covariance matrix into a preset first formula for calculation to obtain a first approximate phase; and carrying out iterative computation on the first approximate phase based on the updated phase likelihood function and a preset first formula to obtain a target estimated phase.
Specifically, a preset sample complex covariance matrix is obtainedSolving the absolute value of the sample complex covariance matrix to obtain the modulus of the sample complex covariance matrix, wherein the absolute value is obtained by a formula (4); then substituting the modulus of the sample complex covariance matrix as a real covariance matrix into a preset first formula for calculation to obtain a first approximate phase; in this embodiment, the first formula is formula (3), based on existing methods, using the norm of the sample complex covariance matrixTo replace the real covariance matrixThe initial phase value, i.e., the first approximate phase, is estimated by equation (3). Substituting the first approximate phase into a preset second formula to obtain an approximate maximum likelihood estimation value; wherein the second formula is formula (9), and the approximate maximum likelihood estimation value is used as a real covariance matrixSubstituting the approximate maximum likelihood estimation value into a preset first formula for calculation, namely using the approximate maximum likelihood estimation value as a real covariance matrixSubstituting the formula (3) to obtain a second approximate phase; substituting the second approximate phase into the second sub-term to obtain a second sub-term value; and repeatedly executing the step of substituting the first approximate phase into a preset second formula to obtain an approximate maximum likelihood estimation value, stopping iteration when the iteration times exceed a certain number (such as 5000 times) or the second subentry value is smaller than a preset threshold (such as 0.001), and substituting the second subentry value when the iteration is stopped into an updated phase likelihood function to obtain a target estimation phase.
In one implementation, the optimizing the updated phase likelihood function to obtain the target estimated phase further includes: carrying out derivation operation on the updated phase likelihood function to obtain a second derivative function; and solving the second derivative function after the zero setting to estimate the phase by the target.
Specifically, a derivative operation is performed on the updated phase likelihood function to obtain a second derivative function; after the second derivative function is set to zero, the following expression is obtained:
wherein the content of the first and second substances,is the real part of the complex number,the imaginary part of the complex number. Solving the above expression can obtain the value of W, and,from this, can obtainTo do soSo that a closed phase can be obtainedI.e. the target estimated phase.
The three methods are used for optimizing the updated phase likelihood function and reliably and accurately solving the phase of the DS point, so that the obtained target estimation phase has high calculation efficiency and calculation accuracy.
In another implementation, the three methods for optimizing the updated phase likelihood function to obtain the target estimated phase may be used in a cross-combination manner, that is, the target estimated phase obtained by the first method may be used as an input of the second method or an input of the third method, an output of the second method may be used as an input of the first method or an input of the third method, and an output of the third method may be used as an input of the first method or an input of the second method. The method of the invention can be extended to other types of radar scatterers.
Exemplary device
As shown in fig. 1, the embodiment of the present invention provides a time-series interference radar phase optimization estimation apparatus, which includes an observation vector acquisition module 401, a phase likelihood function determination module 402, and a target estimation phase determination module 403, wherein,
an observation vector obtaining module 401, configured to obtain an observation vector transmitted by the time-series interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar;
a phase likelihood function determining module 402, configured to obtain, according to the observation vector, a phase likelihood function corresponding to the observation vector;
and a target estimated phase determining module 403, configured to perform maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimated phase.
Based on the above embodiment, the present invention further provides an intelligent terminal, and a schematic block diagram thereof may be as shown in fig. 2. The intelligent terminal comprises a processor, a memory, a network interface, a display screen and a temperature sensor which are connected through a system bus. Wherein, the processor of the intelligent terminal is used for providing calculation and control capability. The memory of the intelligent terminal comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The network interface of the intelligent terminal is used for being connected and communicated with an external terminal through a network. The computer program is executed by a processor to implement a time-series interferometric radar phase optimization estimation method. The display screen of the intelligent terminal can be a liquid crystal display screen or an electronic ink display screen, and the temperature sensor of the intelligent terminal is arranged inside the intelligent terminal in advance and used for detecting the operating temperature of internal equipment.
It will be understood by those skilled in the art that the schematic diagram in fig. 3 is only a block diagram of a part of the structure related to the solution of the present invention, and does not constitute a limitation to the intelligent terminal to which the solution of the present invention is applied, and a specific intelligent terminal may include more or less components than those shown in the figure, or combine some components, or have different arrangements of components.
In one embodiment, an intelligent terminal is provided that includes a memory, and one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for: acquiring an observation vector transmitted by a time sequence interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar;
obtaining a phase likelihood function corresponding to the observation vector according to the observation vector;
and carrying out maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimation phase.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, databases, or other media used in embodiments provided herein may include non-volatile and/or volatile memory. Non-volatile memory can include read-only memory (ROM), Programmable ROM (PROM), Electrically Programmable ROM (EPROM), Electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
In summary, the present invention discloses a method for estimating phase optimization of a time-series interference radar, the method includes: acquiring an observation vector transmitted by a time sequence interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar; obtaining a phase likelihood function corresponding to the observation vector according to the observation vector; and carrying out maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimation phase. According to the embodiment of the invention, the phase quality of the scatterer in the time sequence interference radar can be improved by carrying out maximum likelihood estimation and optimization processing on the phase of the observation vector transmitted by the time sequence interference radar.
Based on the above embodiments, the present invention discloses a method for estimating phase optimization of a time-series interference radar, and it should be understood that the application of the present invention is not limited to the above examples, and it is obvious to those skilled in the art that modifications and variations can be made in the light of the above description, and all such modifications and variations are intended to fall within the scope of the appended claims.
Claims (10)
1. A method for phase-optimized estimation of a time-series interferometric radar, the method comprising:
acquiring an observation vector transmitted by a time sequence interference radar; the observation vector is used for representing a scattering point time dimension observation value transmitted by the time sequence interference radar;
obtaining a phase likelihood function corresponding to the observation vector according to the observation vector;
and carrying out maximum likelihood estimation and optimization processing on the phase likelihood function to obtain a target estimation phase.
2. The method according to claim 1, wherein the obtaining a phase likelihood function corresponding to the observation vector according to the observation vector comprises:
acquiring a real number covariance matrix; the real covariance matrix is used for representing a scattering point real covariance value transmitted by the time sequence interference radar;
performing conditional probability operation on the observation vector and the real number covariance matrix to obtain a conditional probability function;
and carrying out logarithmic operation on the conditional probability function to obtain a phase likelihood function corresponding to the observation vector.
3. The time-series interference radar phase optimization estimation method according to claim 1, wherein the performing maximum likelihood estimation and optimization processing on the phase likelihood function to obtain the target estimation phase comprises:
carrying out derivation operation on the phase likelihood function to obtain a first derivative function;
solving the first derivative function after zero setting to obtain a maximum likelihood estimation value corresponding to a real number covariance matrix;
and optimizing the phase likelihood function based on the maximum likelihood estimation value to obtain a target estimation phase.
4. The time-series interference radar phase optimization estimation method according to claim 3, wherein the optimizing the phase likelihood function based on the maximum likelihood estimation value to obtain the target estimation phase comprises:
substituting the maximum likelihood estimation value into the phase likelihood function to obtain an updated phase likelihood function;
and optimizing the updated phase likelihood function to obtain the target estimated phase.
5. The time-series interference radar phase optimization estimation method of claim 4, wherein the updated phase likelihood function includes a first sub-term, a second sub-term, and a third sub-term; the optimizing the updated phase likelihood function to obtain the target estimated phase includes:
calculating the minimum value of the second sub-item of the updated phase likelihood function based on a preset optimization algorithm;
and substituting the minimum value into the updated phase likelihood function to obtain the target estimated phase.
6. The method according to claim 5, wherein the optimizing the updated phase likelihood function to obtain the target estimated phase further comprises:
acquiring a preset sample complex covariance matrix;
solving the absolute value of the sample complex covariance matrix to obtain the modulus of the sample complex covariance matrix;
substituting the modulus of the sample complex covariance matrix as a real covariance matrix into a preset first formula for calculation to obtain a first approximate phase;
and carrying out iterative computation on the first approximate phase based on the updated phase likelihood function and a preset first formula to obtain a target estimated phase.
7. The time-series interference radar phase optimization estimation method according to claim 6, wherein the iteratively calculating the first approximate phase based on the updated phase likelihood function and a preset first formula to obtain a target estimation phase comprises:
substituting the first approximate phase into a preset second formula to obtain an approximate maximum likelihood estimation value;
substituting the approximate maximum likelihood estimation value serving as a real covariance matrix into a preset first formula for calculation to obtain a second approximate phase;
substituting the second approximate phase into the second sub-term to obtain a second sub-term value;
and repeatedly executing the step of substituting the first approximate phase into a preset second formula to obtain an approximate maximum likelihood estimation value, stopping iteration when the second subentry value is smaller than a preset threshold value, and substituting the second subentry value when the iteration is stopped into an updated phase likelihood function to obtain a target estimated phase.
8. The method according to claim 4, wherein the optimizing the updated phase likelihood function to obtain the target estimated phase further comprises:
carrying out derivation operation on the updated phase likelihood function to obtain a second derivative function;
and solving the second derivative function after the zero setting to estimate the phase by the target.
9. An intelligent terminal comprising a memory, and one or more programs, wherein the one or more programs are stored in the memory, and wherein the one or more programs being configured to be executed by the one or more processors comprises instructions for performing the method of any of claims 1-8.
10. A non-transitory computer-readable storage medium, wherein instructions in the storage medium, when executed by a processor of an electronic device, enable the electronic device to perform the method of any of claims 1-8.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111455165.4A CN113866768B (en) | 2021-12-02 | 2021-12-02 | Time sequence interference radar phase optimization estimation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111455165.4A CN113866768B (en) | 2021-12-02 | 2021-12-02 | Time sequence interference radar phase optimization estimation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113866768A true CN113866768A (en) | 2021-12-31 |
CN113866768B CN113866768B (en) | 2022-04-15 |
Family
ID=78985514
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111455165.4A Active CN113866768B (en) | 2021-12-02 | 2021-12-02 | Time sequence interference radar phase optimization estimation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113866768B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115267772A (en) * | 2022-05-30 | 2022-11-01 | 杭州电子科技大学 | Self-adaptive multi-temporal interferometry method and system based on complex covariance matrix |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2213025T3 (en) * | 1999-05-25 | 2004-08-16 | Politecnico Di Milano | PROCEDURE FOR MEASURING THE MOVEMENT OF URBAN AREAS AND LAND RUNNING AREAS. |
CN103616686A (en) * | 2013-12-05 | 2014-03-05 | 中国测绘科学研究院 | Optimal phase-position estimating method for complete polarization interferometric synthetic aperture radar image based on mixed mode |
CN103713287A (en) * | 2013-12-26 | 2014-04-09 | 中国科学院电子学研究所 | Elevation reestablishing method and device based on coprime of multiple base lines |
CN104808203A (en) * | 2015-03-03 | 2015-07-29 | 电子科技大学 | Multi-baseline InSAR phase unwrapping method by iterating maximum likelihood estimation |
CN106950556A (en) * | 2017-05-03 | 2017-07-14 | 三亚中科遥感研究所 | Heritage area deformation monitoring method based on distributed diffusion body sequential interference SAR technology |
CN108051810A (en) * | 2017-12-01 | 2018-05-18 | 南京市测绘勘察研究院股份有限公司 | A kind of InSAR distributed diffusions body phase optimization method |
CN108535698A (en) * | 2018-04-04 | 2018-09-14 | 西安电子科技大学 | The low elevation estimate method of metre wave radar based on beam space |
CN109932698A (en) * | 2019-03-10 | 2019-06-25 | 西安电子科技大学 | The low elevation estimate method of metre wave radar based on terrain information |
CN110763187A (en) * | 2019-09-30 | 2020-02-07 | 中国科学院测量与地球物理研究所 | Stable ground settlement monitoring method based on radar distributed target |
CN111856459A (en) * | 2020-06-18 | 2020-10-30 | 同济大学 | Improved DEM maximum likelihood constraint multi-baseline InSAR phase unwrapping method |
CN112034457A (en) * | 2020-07-21 | 2020-12-04 | 西安电子科技大学 | Multi-baseline elevation interference phase estimation method based on interference fringe direction |
-
2021
- 2021-12-02 CN CN202111455165.4A patent/CN113866768B/en active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2213025T3 (en) * | 1999-05-25 | 2004-08-16 | Politecnico Di Milano | PROCEDURE FOR MEASURING THE MOVEMENT OF URBAN AREAS AND LAND RUNNING AREAS. |
CN103616686A (en) * | 2013-12-05 | 2014-03-05 | 中国测绘科学研究院 | Optimal phase-position estimating method for complete polarization interferometric synthetic aperture radar image based on mixed mode |
CN103713287A (en) * | 2013-12-26 | 2014-04-09 | 中国科学院电子学研究所 | Elevation reestablishing method and device based on coprime of multiple base lines |
CN104808203A (en) * | 2015-03-03 | 2015-07-29 | 电子科技大学 | Multi-baseline InSAR phase unwrapping method by iterating maximum likelihood estimation |
CN106950556A (en) * | 2017-05-03 | 2017-07-14 | 三亚中科遥感研究所 | Heritage area deformation monitoring method based on distributed diffusion body sequential interference SAR technology |
CN108051810A (en) * | 2017-12-01 | 2018-05-18 | 南京市测绘勘察研究院股份有限公司 | A kind of InSAR distributed diffusions body phase optimization method |
CN108535698A (en) * | 2018-04-04 | 2018-09-14 | 西安电子科技大学 | The low elevation estimate method of metre wave radar based on beam space |
CN109932698A (en) * | 2019-03-10 | 2019-06-25 | 西安电子科技大学 | The low elevation estimate method of metre wave radar based on terrain information |
CN110763187A (en) * | 2019-09-30 | 2020-02-07 | 中国科学院测量与地球物理研究所 | Stable ground settlement monitoring method based on radar distributed target |
CN111856459A (en) * | 2020-06-18 | 2020-10-30 | 同济大学 | Improved DEM maximum likelihood constraint multi-baseline InSAR phase unwrapping method |
CN112034457A (en) * | 2020-07-21 | 2020-12-04 | 西安电子科技大学 | Multi-baseline elevation interference phase estimation method based on interference fringe direction |
Non-Patent Citations (6)
Title |
---|
CHRISTOPHE MAGNARD等: "Analysis of a Maximum Likelihood Phase Estimation Method for Airborne Multibaseline SAR", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING》 * |
M.S.SEYMOUR等: "Maximum Likelihood Estimation For SAR Interferometry", 《PROCEEDINGS OF IGARSS" 94-1994 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM》 * |
MI JIANG等: "Distributed Scatterer Interferometry With the Refinement of Spatiotemporal Coherence", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
张秋玲等: "利用多基线数据融合提高分布式卫星InSAR系统的干涉相位精度", 《电子与信息学报》 * |
李毅: "融合分布式目标的矿区长时序InSAR地表形变监测", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑(月刊)》 * |
杨康等: "基于相位估计的InSAR信号优化处理", 《指挥信息系统与技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115267772A (en) * | 2022-05-30 | 2022-11-01 | 杭州电子科技大学 | Self-adaptive multi-temporal interferometry method and system based on complex covariance matrix |
CN115267772B (en) * | 2022-05-30 | 2023-08-29 | 杭州电子科技大学 | Self-adaptive multi-time interference measurement method and system based on complex covariance matrix |
Also Published As
Publication number | Publication date |
---|---|
CN113866768B (en) | 2022-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cao et al. | Infrared dim target detection via mode-k1k2 extension tensor tubal rank under complex ocean environment | |
CN113866768B (en) | Time sequence interference radar phase optimization estimation method | |
CN116433707B (en) | Accurate extraction method and system for optical center sub-pixels of line structure under complex background | |
Peng et al. | Infrared small-target detection based on multi-directional multi-scale high-boost response | |
CN114494371A (en) | Optical image and SAR image registration method based on multi-scale phase consistency | |
Xu et al. | A new shadow tracking method to locate the moving target in SAR imagery based on KCF | |
Havangi | Intelligent adaptive unscented particle filter with application in target tracking | |
Sirish Kumar et al. | Implementation of new navigation algorithm based on cross-correntropy for precise positioning in low latitude regions of South India | |
CN116381686A (en) | Terahertz video synthetic aperture radar moving target imaging method based on time-frequency analysis | |
CN116047459A (en) | Array radar echo signal recovery method and related equipment in pulse interference environment | |
Zhou et al. | High precision cross-range scaling and 3D geometry reconstruction of ISAR targets based on geometrical analysis | |
CN114743150A (en) | Target tracking method and device, electronic equipment and storage medium | |
Zhan et al. | SAR image super-resolution reconstruction based on an optimize iterative method for regularization | |
Musso et al. | A Laplace-based particle filter for track-before-detect | |
WO2024113425A1 (en) | Time series sar data filtering method suitable for farmland fragmentation region | |
Ivanova et al. | Restoration of orientation distribution function using texture components with radial normal distributions | |
Denisenko et al. | Reconstruction of the height profiles of the electron concentration based on vertical sounding data with the IRI model | |
Wang et al. | Three‐dimensional point cloud reconstruction of inverse synthetic aperture radar image sequences based on back projection and iterative closest point fusion | |
Gao et al. | Exploitation of SRTM DEM in InSAR data processing and its application to phase unwrapping | |
Li et al. | LSD and skeleton extraction combined with farmland ridge detection | |
Li et al. | Improved Doppler parameter estimation of squint SAR based on slope detection | |
CN112130144B (en) | Microwave correlation imaging method and imaging system based on dynamic grid | |
Gao et al. | Parallel processing of sliding spotlight mode SAR imaging based on GPU | |
Sun et al. | Extended Target Tracking Using Non-linear Observations | |
Chen et al. | Adaptive clutter nulling approach for heterogeneous environments |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |