CN111325682B - Phase recovery improvement method and device for self-correlation signal with interference - Google Patents
Phase recovery improvement method and device for self-correlation signal with interference Download PDFInfo
- Publication number
- CN111325682B CN111325682B CN202010071872.2A CN202010071872A CN111325682B CN 111325682 B CN111325682 B CN 111325682B CN 202010071872 A CN202010071872 A CN 202010071872A CN 111325682 B CN111325682 B CN 111325682B
- Authority
- CN
- China
- Prior art keywords
- reconstruction result
- fourier
- phase recovery
- reconstruction
- target object
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000011084 recovery Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000006872 improvement Effects 0.000 title abstract description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Discrete Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
A phase recovery improvement method and device for a band interference autocorrelation signal, the method includes: extracting an initial Fourier amplitude of a target object from a known target object autocorrelation signal, and taking a random Fourier phase as an input to obtain a phase recovery reconstruction result; performing space domain denoising treatment on the reconstruction result, assigning a part smaller than the maximum energy setting percentage of the reconstruction result to zero to obtain a denoised reconstruction result, and calculating the mean square error between the Fourier amplitude of the denoised reconstruction result and the initial Fourier amplitude of the target object; repeating the steps by taking a plurality of different random Fourier phases as inputs to obtain a plurality of de-noised reconstruction results and a plurality of mean square errors; and taking the mean square error as a judging index of the quality of the reconstruction result, selecting the optimal reconstruction result from the plurality of denoised reconstruction results as a final output, and realizing the phase recovery of the target object. The invention can realize stable phase recovery for the self-correlation signal with interference.
Description
Technical Field
The invention relates to the field of computer vision and digital image processing, in particular to an improved phase recovery method and device for an auto-correlation signal with interference.
Background
The traditional phase recovery algorithm can be used for recovering and reconstructing the object to be observed from the autocorrelation signal of the object, namely, firstly, the Fourier transformation is utilized to extract the Fourier amplitude of the object from the autocorrelation signal; under the condition that the Fourier amplitude of the target object is known and the phase is unknown, taking the random phase as an initial guess of the Fourier phase of the target object, and iteratively recovering the real Fourier phase of the target object and reconstructing the spatial distribution of the target object. However, due to randomness of initial phase guess, the spatial distribution of the target object reconstructed by the traditional phase recovery algorithm and the actual fourier phase of the recovered target object have great instability, and particularly when partial interference exists in the autocorrelation signal of the target object, the reconstruction instability is larger, so that the phase recovery result is disordered, the accurate and effective spatial information of the target object is difficult to obtain from the single phase recovery result, and the applicability of the traditional phase recovery algorithm to the autocorrelation signal with interference is limited.
Disclosure of Invention
The main object of the present invention is to overcome at least one of the above technical drawbacks and to provide an improved method and device for phase recovery for an auto-correlation signal with interference.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
a phase recovery method for a band-disturbance-oriented autocorrelation signal, said method comprising the steps of:
a1: extracting an initial Fourier amplitude of a target object from a known target object autocorrelation signal, and taking a random Fourier phase as an input to obtain a phase recovery reconstruction result;
a2: performing space domain denoising treatment on the reconstruction result, assigning a part smaller than the maximum energy setting percentage of the reconstruction result to zero to obtain a denoised reconstruction result, and calculating the mean square error between the Fourier amplitude of the denoised reconstruction result and the initial Fourier amplitude of the target object;
a3: repeating the steps A1 and A2 by taking a plurality of different random Fourier phases as inputs to obtain a corresponding plurality of de-noised reconstruction results and a corresponding plurality of mean square errors;
a4: and taking the mean square error as a judging index of the quality of the reconstruction result, selecting an optimal reconstruction result from the plurality of denoised reconstruction results as a final output, and realizing the phase recovery of the target object.
Further:
in step A2, the part smaller than 20% of the maximum energy of the reconstruction result is assigned zero.
In step A1, for the initial target object autocorrelation signal with interference, the process of extracting the initial fourier amplitude of the target object by fourier transform and evolution operation is as follows:
wherein ,representing an initial target autocorrelation signal with interference, c representing the dryDisturbing signal, O represents the object to be reconstructed, < ->Representing an autocorrelation operation. />
In step A1, the initial Fourier amplitude |F (O) | of the object is known, and the random Fourier phase theta is used m As an initial guess of the phase, iteratively reconstructing the spatial distribution O of the target object m And obtaining a phase recovery reconstruction result.
In step A2, for the reconstruction result O m Performing space domain denoising processing, and assigning a part of which the maximum energy is less than 20% of the reconstruction result to be zero to obtain a denoised reconstruction result P m According to the following formula:
wherein (x, y) represents the two-dimensional coordinate distribution of the target airspace.
In step A2, for the denoised reconstruction result P m Calculating the Fourier amplitude |F (P) m ) Mean square error between i and initial fourier magnitude of the object, F (O)According to the following formula:
wherein, (s, t) represents the two-dimensional coordinate distribution of the Fourier domain of the target object; H. v represents the maximum range of the two-dimensional coordinates of the fourier domain of the object, respectively.
In step A3, a plurality of different random Fourier phases θ m Repeating steps A1 and A2 with m=1, 2, …, M as input to obtain a plurality of corresponding denoised reconstruction results P m M=1, 2, …, M and corresponding multiple mean square errorsIn step A4, the mean square error is used as a criterion of the quality of the reconstruction result, and the optimal reconstruction result is selected from the plurality of denoised reconstruction results according to the following formula:
where m=1, 2, …, M represents the number of random fourier phases input.
In step A1, the initial target autocorrelation signal a with interference is taken as the known target autocorrelation signal.
A phase recovery device for a signal with interference autocorrelation, comprising at least one memory and at least one processor;
the memory stores at least one executable program;
the executable program, when executed by the processor, implements the method.
The invention has the following beneficial effects:
the invention provides a phase recovery method and a phase recovery device for an autocorrelation signal with interference, which can realize stable and effective reconstruction of spatial distribution of a target object and Fourier phase of the target object under the condition that the autocorrelation signal is known to have certain interference. The method comprises the steps of carrying out spatial domain denoising on a reconstruction result of a traditional phase recovery algorithm according to a certain percentage proportion of maximum energy, particularly preferably 20%, taking constraint as a constraint, calculating a Mean Square Error (MSE) between a Fourier amplitude of the denoised reconstruction result and a Fourier amplitude of an initial target object, taking the MSE as a judging standard of the denoised reconstruction result, and selecting an optimal reconstruction result from a plurality of denoising reconstruction results with different random phases as inputs so as to realize stable phase recovery and target reconstruction. By using the phase recovery method and the phase recovery device, stable phase recovery can be realized for the autocorrelation signal with interference.
Drawings
Fig. 1 is a flowchart of a phase recovery method for an auto-correlation signal with interference according to an embodiment of the present invention.
Detailed Description
The following describes embodiments of the present invention in detail. It should be emphasized that the following description is merely exemplary in nature and is in no way intended to limit the scope of the invention or its applications.
Referring to fig. 1, an embodiment of the present invention provides a phase recovery method for an auto-correlation signal with interference, which includes the following steps:
a1: extracting an initial Fourier amplitude of a target object from a known target object autocorrelation signal, and taking a random Fourier phase as an input to obtain a phase recovery reconstruction result;
a2: performing space domain denoising treatment on the reconstruction result, assigning a part of which the maximum energy is less than 20% of the reconstruction result to be zero, obtaining a denoised reconstruction result, and calculating the mean square error between the Fourier amplitude of the denoised reconstruction result and the initial Fourier amplitude of the target object;
a3: repeating the steps A1 and A2 by taking a plurality of different random Fourier phases as inputs to obtain a corresponding plurality of de-noised reconstruction results and a corresponding plurality of mean square errors;
a4: and taking the mean square error as a judging index of the quality of the reconstruction result, selecting the optimal reconstruction result from a plurality of reconstruction results as a final output, and realizing the phase recovery of the target object.
In the step A1, the phase recovery method may be implemented using the initial target autocorrelation signal a with interference as a known target autocorrelation signal.
In a specific preferred embodiment, in step A1, for an initial target autocorrelation signal with interference, the process of extracting the initial fourier amplitude of the target using fourier transform and squaring operation is as follows:
wherein ,an object autocorrelation signal representing an initial disturbance, c representing a disturbance signal, O representing an object to be reconstructed,/->Representing an autocorrelation operation.
Further preferably, the initial Fourier amplitude |F (O) | of the object is known as the random Fourier phase θ m As an initial guess of the phase, iteratively reconstructing the spatial distribution O of the target object m And obtaining a phase recovery reconstruction result.
In a particularly preferred embodiment, in step A2, the fraction smaller than 20% of the maximum energy of the reconstruction result is assigned zero.
In a particularly preferred embodiment, in step A2, for the reconstruction result O m Performing space domain denoising processing, and assigning a part of which the maximum energy is less than 20% of the reconstruction result to be zero to obtain a denoised reconstruction result P m :
Wherein (x, y) represents the two-dimensional coordinate distribution of the target airspace.
Further preferably, for the denoised reconstruction result P m Calculating the Fourier amplitude |F (P) m ) Mean square error between i and initial fourier magnitude of the object, F (O)/>
Wherein, (s, t) represents the two-dimensional coordinate distribution of the Fourier domain of the target object; H. v represents the maximum range of the two-dimensional coordinates of the fourier domain of the object, respectively.
In a particularly preferred embodiment, in step A3, a plurality of different random Fourier phases θ m Repeating steps A1 and A2 with m=1, 2, …, M as input to obtain a plurality of corresponding denoised reconstruction results P m M=1, 2, …, M and corresponding multiple mean square errors
In step A4, the mean square error is used as a judging index of the quality of the reconstruction result, the optimal reconstruction result is selected from a plurality of reconstruction results, and the following formula is adopted:
where m=1, 2, …, M represents the number of random fourier phases input.
The phase recovery method for the self-correlation signal with interference provided by the embodiment of the invention can realize stable and effective reconstruction of the spatial distribution of the target object and the Fourier phase of the target object under the condition that the self-correlation signal with certain interference is known. The method comprises the steps of carrying out spatial domain denoising on a reconstruction result of a traditional phase recovery algorithm according to a certain percentage proportion of maximum energy, particularly preferably 20%, taking constraint as a constraint, calculating a Mean Square Error (MSE) between a Fourier amplitude of the denoised reconstruction result and a Fourier amplitude of an initial target object, taking the MSE as a judging standard of the denoised reconstruction result, and selecting an optimal reconstruction result from a plurality of denoising reconstruction results with different random phases as inputs so as to realize stable phase recovery and target reconstruction. By using the phase recovery method and the phase recovery device provided by the embodiment of the invention, stable phase recovery can be realized for the self-correlation signal with interference.
It should be noted that the specific preferred methods (such as phase recovery algorithm, generation of random phase, etc.) in the above embodiments are only illustrative, and the scope of the present invention is not limited to the methods.
The background section of the present invention may contain background information about the problems or environments of the present invention and is not necessarily descriptive of the prior art. Accordingly, inclusion in the background section is not an admission of prior art by the applicant.
The foregoing is a further detailed description of the invention in connection with specific/preferred embodiments, and it is not intended that the invention be limited to such description. It will be apparent to those skilled in the art that several alternatives or modifications can be made to the described embodiments without departing from the spirit of the invention, and these alternatives or modifications should be considered to be within the scope of the invention. In the description of the present specification, reference to the terms "one embodiment," "some embodiments," "preferred embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Those skilled in the art may combine and combine the features of the different embodiments or examples described in this specification and of the different embodiments or examples without contradiction. Although embodiments of the present invention and their advantages have been described in detail, it should be understood that various changes, substitutions and alterations can be made herein without departing from the scope of the invention as defined by the appended claims.
Claims (6)
1. A method for phase recovery of an auto-correlation signal with interference, the method comprising the steps of:
a1: extracting an initial Fourier amplitude of a target object from an initial target object autocorrelation signal with interference, and taking a random Fourier phase as input to obtain a phase recovery reconstruction result;
a2: performing space domain denoising treatment on the reconstruction result, assigning a part which is smaller than 20% of the maximum energy of the reconstruction result to be zero, obtaining a denoised reconstruction result, and calculating the mean square error between the Fourier amplitude of the denoised reconstruction result and the initial Fourier amplitude of the target object;
in step A2, for the reconstruction result O m Performing space domain denoising processing, and assigning a part of which the maximum energy is less than 20% of the reconstruction result to be zero to obtain a denoised reconstruction result P m According to the following formula:
wherein, (x, y) represents two-dimensional coordinate distribution of the airspace of the target object;
a3: repeating the steps A1 and A2 by taking a plurality of different random Fourier phases as inputs to obtain a corresponding plurality of de-noised reconstruction results and a corresponding plurality of mean square errors;
a4: and taking the mean square error as a judging index of the quality of the reconstruction result, selecting an optimal reconstruction result from the plurality of denoised reconstruction results as a final output, and realizing the phase recovery of the target object.
2. The phase recovery method for an auto-correlation signal with interference according to claim 1, wherein in step A1, for an initial target auto-correlation signal with interference, the process of extracting an initial fourier amplitude of the target by using fourier transform and squaring operation is as follows:
3. The phase recovery method for an auto-correlation signal with interference according to claim 2, wherein in step A1, the initial fourier amplitude |f (O) | of the target object is known, and the random fourier phase θ is used m As an initial guess of the phase, iteratively reconstructing the spatial distribution O of the target object m And obtaining a phase recovery reconstruction result.
4. The phase recovery method for an interfering autocorrelation signal of claim 1 wherein in step A2, for the denoised reconstruction result P m Calculating the Fourier amplitude |F (P) m ) Mean square error between i and initial fourier magnitude of the object, F (O)According to the following formula:
wherein, (s, t) represents the two-dimensional coordinate distribution of the Fourier domain of the target object; H. v represents the maximum range of the two-dimensional coordinates of the fourier domain of the object, respectively.
5. The phase recovery method for an interfering autocorrelation signal of claim 4 wherein in step A3, a plurality of different random fourier phases θ are used m Repeating steps A1 and A2 with m=1, 2, … and m as input to obtain a plurality of corresponding denoised reconstruction results P m M=1, 2, …, M and corresponding multiple mean square errorsIn step A4, toThe mean square error is used as a judging index of the quality of the reconstruction result, the optimal reconstruction result is selected from the plurality of denoised reconstruction results, and the following formula is adopted:
where m=1, 2, …, M represents the number of random fourier phases input.
6. A phase recovery device for a signal with interference autocorrelation, comprising at least one memory and at least one processor;
the memory stores at least one executable program;
the executable program, when executed by the processor, implements the method of any one of claims 1 to 5.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010071872.2A CN111325682B (en) | 2020-01-21 | 2020-01-21 | Phase recovery improvement method and device for self-correlation signal with interference |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010071872.2A CN111325682B (en) | 2020-01-21 | 2020-01-21 | Phase recovery improvement method and device for self-correlation signal with interference |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111325682A CN111325682A (en) | 2020-06-23 |
CN111325682B true CN111325682B (en) | 2023-05-23 |
Family
ID=71167225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010071872.2A Active CN111325682B (en) | 2020-01-21 | 2020-01-21 | Phase recovery improvement method and device for self-correlation signal with interference |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111325682B (en) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106204466A (en) * | 2016-06-24 | 2016-12-07 | 南京理工大学 | A kind of self-adaptive solution method for Fourier lamination micro-imaging technique |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10393865B2 (en) * | 2013-12-11 | 2019-08-27 | Airbus Sas | Phase retrieval algorithm for generation of constant time envelope with prescribed fourier transform magnitude signal |
CN107004259A (en) * | 2014-12-10 | 2017-08-01 | 皇家飞利浦有限公司 | Statistical weight regularization in many contrast imagings |
CN108680953B (en) * | 2018-05-16 | 2019-07-26 | 中国海洋石油集团有限公司 | A kind of seismic data based on inverse proportion model interpolation and denoising method simultaneously |
CN110047048B (en) * | 2019-04-17 | 2021-03-02 | 清华大学深圳研究生院 | Phase recovery improved algorithm based on MSE (mean square error) optimization |
-
2020
- 2020-01-21 CN CN202010071872.2A patent/CN111325682B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106204466A (en) * | 2016-06-24 | 2016-12-07 | 南京理工大学 | A kind of self-adaptive solution method for Fourier lamination micro-imaging technique |
Also Published As
Publication number | Publication date |
---|---|
CN111325682A (en) | 2020-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7163504B2 (en) | Image processing method and its apparatus, computer program and electronic equipment | |
CN103971346B (en) | SAR (Synthetic Aperture Radar) image spot-inhibiting method based on spare domain noise distribution constraint | |
CN110018438B (en) | Direction-of-arrival estimation method and device | |
CN101141560A (en) | Synthetic aperture radar image noise-eliminating method based on independent component analysis based image | |
CN111010356A (en) | Underwater acoustic communication signal modulation mode identification method based on support vector machine | |
JPWO2020091891A5 (en) | ||
Yang et al. | One-dimensional deep attention convolution network (ODACN) for signals classification | |
CN102184451A (en) | Multi-wavelet fusion feature based neuron action feature extraction method | |
Wang et al. | Low probability of intercept radar waveform recognition based on dictionary leaming | |
CN117692074B (en) | Low-frequency aliasing noise suppression method suitable for unsteady-state underwater sound target signal | |
CN111325682B (en) | Phase recovery improvement method and device for self-correlation signal with interference | |
CN110047048B (en) | Phase recovery improved algorithm based on MSE (mean square error) optimization | |
CN114565772B (en) | Method and device for extracting set features, electronic equipment and storage medium | |
CN114358050B (en) | Intelligent radar radiation source identification method based on bicubic interpolation and WVD feature square matrix | |
CN114047481A (en) | Robust adaptive beam forming method based on subspace orthogonality | |
CN110332987B (en) | Voiceprint signal imaging method and microphone array signal imaging method | |
CN110473210B (en) | Image segmentation method and device based on confidence propagation | |
CN108256581B (en) | Gabor wavelet domain copula model image classification method | |
CN102509268A (en) | Immune-clonal-selection-based nonsubsampled contourlet domain image denoising method | |
CN113553901A (en) | Improved synchronous extraction time-frequency analysis method based on self-adaptive window length | |
Salehi et al. | A robust hybrid filter based on evolutionary intelligence and fuzzy evaluation | |
US20240153118A1 (en) | Method and device for estimating a depth map associated with a digital hologram representing a scene and computer program associated | |
CN117706514B (en) | Clutter elimination method, system and equipment based on generation countermeasure network | |
CN117872379B (en) | Underwater target ranging method, medium and system under shallow sea internal wave condition | |
CN115932732B (en) | Combined underwater sound multi-target blind separation method based on improved NMF and FastICA |
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 |