CN113640793A - MRF-based real aperture scanning radar super-resolution imaging method - Google Patents

MRF-based real aperture scanning radar super-resolution imaging method Download PDF

Info

Publication number
CN113640793A
CN113640793A CN202110765492.3A CN202110765492A CN113640793A CN 113640793 A CN113640793 A CN 113640793A CN 202110765492 A CN202110765492 A CN 202110765492A CN 113640793 A CN113640793 A CN 113640793A
Authority
CN
China
Prior art keywords
super
resolution imaging
objective function
target
mrf
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
Application number
CN202110765492.3A
Other languages
Chinese (zh)
Other versions
CN113640793B (en
Inventor
谭珂
杨建超
陆星宇
张婷婷
苏卫民
顾红
戴峥
虞文超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN202110765492.3A priority Critical patent/CN113640793B/en
Publication of CN113640793A publication Critical patent/CN113640793A/en
Application granted granted Critical
Publication of CN113640793B publication Critical patent/CN113640793B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging

Abstract

The invention discloses a MRF-based real aperture scanning radar super-resolution imaging method, which comprises the following steps: introducing structural prior information represented by a hidden Markov random field into an objective function through a regularization frame; and solving the optimization problem of the objective function by using a fast iterative threshold contraction method, and finally obtaining the super-resolution imaging result of the target. According to the method, the two-dimensional spatial correlation of the pixels can be better described by introducing the Markov random field prior model, so that the shape of a scene can be better recovered.

Description

MRF-based real aperture scanning radar super-resolution imaging method
Technical Field
The invention belongs to the technical field of radars, and particularly relates to a real aperture scanning radar super-resolution imaging method based on MRF.
Background
Forward looking radar imaging has a wide range of application scenarios such as aircraft navigation, landing and airport surveillance. However, the sharply decreased doppler bandwidth and left/right doppler ambiguity in the forward looking direction make the conventional single station synthetic aperture radar and doppler beam sharpening techniques difficult to apply in forward looking imaging.
The convolution inversion technology carries out super-resolution processing on RASR from the signal processing angle, and can break through the limitation of radar real aperture on the position resolution. The method has the advantages of simple imaging mode, high compatibility, no increase of physical aperture, capability of distinguishing a plurality of targets in the same beam and the like, and gradually becomes a great research hotspot of radar forward-looking imaging.
However, many deconvolution methods are sensitive to noise due to the ill-conditioned nature of the deconvolution problem. The regularization method can effectively relieve the ill-posed property of the convolution inversion problem, constructs different regularization punishment items according to the known prior information to apply constraint on a target solution, and can solve the problem that the convolution inversion solution is sensitive to noise and discontinuously depends on observation data. For example, the document "Iterative non-systematic angular super resolution" (Richards M. radius Conference, 1988, Proceedings of the 1988 IEEE national. IEEE, 1988:100-1The norm constraint regularization algorithm carries out super-resolution imaging on the target, because l1The sparse representation capability of the norm on the solution can effectively distinguish sparse targets in a scene, but the background contour information is damaged when the number of optimization iterations is too many; the document "Bayesian deconstruction for Angular Super-Resolution in Forward-Looking Scanning Raar" (ZHa Y, Huang Y, Sun Z, et al. Sensors,2015,15(3):6924-6946) utilizes a Super-Resolution algorithm of total variation plus sparse constraint, which can better retain the contour of a scene while distinguishing sparse objects. However, due to the limitation of the sequential processing mode according to the distance units, the methods only use the prior information in the one-dimensional distance units. In the imaging task, each pixel is considered as a result determined by surrounding pixels, so the scene usually has a two-dimensional structure characteristic, and therefore, the two-dimensional structure prior information of the scene is still to be mined and utilized.
Disclosure of Invention
The invention aims to provide an RASR super-resolution imaging method based on MRF (Markov random field) aiming at the defects in the background technology.
The technical solution for realizing the purpose of the invention is as follows: a MRF-based real aperture scanning radar super-resolution imaging method comprises the following steps:
introducing structural prior information represented by a hidden Markov random field into an objective function through a regularization frame;
and solving the optimization problem of the objective function by using a fast iterative threshold contraction method, and finally obtaining the super-resolution imaging result of the target.
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the above-mentioned super-resolution imaging method when executing the computer program.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the above-mentioned super-resolution imaging method.
Compared with the prior art, the invention has the following remarkable advantages: the invention introduces the structural prior information expressed by the hidden Markov random field into the target function through the regularization frame, then solves the optimization problem of the target function by using the FIST method, and finally obtains the super-resolution imaging result of the target. The method has the advantages that the two-dimensional spatial correlation of the pixels can be better described by introducing the Markov random field prior model, so that the shape of the scene can be better recovered.
Drawings
FIG. 1 is a block flow diagram of a method provided by the present invention.
Fig. 2 is a diagram of the spatial relationship of motion scanning of a real aperture scanning radar used in the embodiment of the present invention.
FIG. 3 is a schematic diagram of a simulated imaging scene used in accordance with an embodiment of the present invention.
Fig. 4 is an echo signal superimposed with white gaussian noise according to an embodiment of the present invention.
Fig. 5 shows the signal-to-noise ratio of the signal after correction for range-wise compression and migration, which is 18 dB.
Fig. 6 is an enlarged view of the target area of fig. 5.
FIG. 7 is a graphical representation of the imaging results after processing by the method of the present invention and by a comparison algorithm.
Detailed Description
The solution of the invention is to apply a Markov random field prior model to the real aperture scanning radar super-resolution imaging. Introducing structural prior information expressed by a hidden Markov random field into a target function through a regularization frame, then solving the optimization problem of the target function by using a fast iterative threshold shrinkage (FIST) method, and finally obtaining a super-resolution imaging result of the target.
For the convenience of describing the contents of the present invention, the following terms are first explained:
the term 1: real aperture scanning radar
The Real Aperture Scanning Radar (RASR) is a radar system relative to a Synthetic Aperture Radar (SAR), and is different from the SAR in that a synthetic aperture mode is adopted in the azimuth direction to improve the resolution, a beam system is adopted in the azimuth direction to obtain a two-dimensional echo, no signal processing is carried out in the azimuth direction, and only real beams are utilized to distinguish a target.
The term 2: markov random field
Markov (MRF) means that under certain conditions dictated by "current", the result produced by the "last" will not affect the "future" result. In the one-dimensional case, the state of a random chain is only related to the state of its neighboring pixel points, that is, only affected by the state of the neighboring pixel points. For a two-dimensional image, the image may be defined in a two-dimensional space while being considered as a random field and a two-dimensional markov random field with markov properties.
The invention provides a RASR super-resolution imaging method, as shown in figure 1, comprising:
the method comprises the following steps: imaging system parameter initialization
The platform moving speed is v, the antenna wave beam scanning speed is omega, and the radar lower visual angle is
Figure BDA0003150886440000031
The coordinates of the target point are denoted as P (r)00) T is azimuth time, target point distance history rP(x, y) is
Figure BDA0003150886440000032
Due to the fast scanning speed and small imaging sector, the slope distance obtained by the approximation of Taylor expansion and first order term is
Figure BDA0003150886440000033
The expression of the point target echo is recorded as sP(τ,t),
Figure BDA0003150886440000034
Where τ is the distance-wise time variable, σ0For point target scattering intensity, h (t) for the dual-pass beam pattern of the antenna, t0Rect [. for the time when the beam center hits point P]Representing the distance time window, c is the speed of light, TrRepresenting the distance-time pulse width, f0Is the carrier frequency, KrIs the time chirp rate of the transmitted signal.
Step two: radial pulse compression
For a scattering function of
Figure BDA0003150886440000041
For a surface target, the echo is the integral of the echoes of all point targets, wherein
Figure BDA0003150886440000042
Is the initial slope distance of the target,
Figure BDA0003150886440000043
is the azimuth coordinate of the target. Because t is (theta-theta)a) ω and τ 2r/c, where θaFor the initial scan angle, the surface target echo can be obtained according to equations (2) and (3):
Figure BDA0003150886440000044
aiming at the acquired two-dimensional echo data s (r, theta), a matched filter is constructed by using radar emission signal parameters, the range direction is subjected to pulse compression, and the data s after range compression is obtainedout(r,θ)。
Figure BDA0003150886440000045
Step three: linear migration correction
After compressing in the range direction, for the target plane
Figure BDA0003150886440000046
Up located on arc line
Figure BDA0003150886440000047
The center of the echo position of all point targets is overlapped on a straight line
Figure BDA0003150886440000048
The above. For the subsequent super-resolution processing, the echo needs to be corrected by linear migration, and the order is
Figure BDA0003150886440000049
The variable r is scaled. Eliminating range migration corrected echo as
Figure BDA00031508864400000410
Step four: constructing an azimuthal convolution signal model
Only the bearing signal amplitude is studied, the bearing signal can be written as,
s(θ)=σ(θ)*h(θ) (6)
where is the convolution operation. Considering the influence of noise, the azimuth signal expression is,
s(θ)=σ(θ)*h(θ)+n(θ) (7)
and (3) constructing an azimuth convolution signal vector model aiming at the model:
s=Hσ+n (8)
wherein s ═ s1,…,sN]TAnd σ ═ σ [ σ ]1,...,σN]TRespectively representing the echo data and the target scattering coefficient in a range gate, n ═ n1,...,nN]TN is the number of vectors, and H is the convolution matrix constructed by the antenna directional diagram H (theta).
Step five: constructing an imaging objective function
For the problem of noise sensitivity in deconvolution operation, a regularization method is adopted to solve the problem, and a standard regularization equation for solving the problem is as follows:
Figure BDA0003150886440000051
wherein the content of the first and second substances,
Figure BDA0003150886440000052
for data fidelity terms, | Γ (σ) | purpleqFor the regularization term, Γ (σ) is the specific operation of the regularization term, q is the norm of the regularization term, and λ is used to balance the fidelity of the measured data with the effect of the constraints of the regularization term, which can be determined by the L-curve method.
In this patent, Huber-Markov is chosen as the regularization term
Figure BDA0003150886440000053
Where ρ isT(. cndot.) is a function of Huber,
Figure BDA0003150886440000056
is a coefficient vector of the cluster c, tau is a temperature parameter, Z is a normalized constant, superscript (·)tDenoted as a transpose operation. The roughness of the image is obtained from the second derivative:
Figure BDA0003150886440000054
the Huber function is defined as:
Figure BDA0003150886440000055
the threshold T penalizes the gray scale variation of the image. The square penalty term smoothes out small-scale noise, e.g., σ ≦ T, for regions with smooth gray-scale changes, while the linear penalty term is used to smooth out boundary regions with strong gray-scale changes in the image, e.g., σ > T.
Substituting equations (10), (11) and (12) into equation (9), the objective function of RASR super-resolution imaging can become:
Figure BDA0003150886440000061
step six: determining the gradient of an objective function
Because the obtained objective function is a convex function, the optimization problem of the formula (13) is solved by adopting a fast iterative shrinkage/threshold (FIST) algorithm to obtain the gradient of the objective function
Figure BDA0003150886440000062
Wherein, (.)mThe mth element of the vector in parentheses, diag (. cndot.) is the diagonal matrix, ε is a solutionSmall constants introduced due to the non-microminiaturization caused by the non-smoothness of the objective function. The one-time iterative threshold shrinkage operation process is as follows:
Figure BDA0003150886440000063
wherein the content of the first and second substances,
Figure BDA0003150886440000064
is a shrink operator, defined as:
Figure BDA0003150886440000065
step seven: fast iterative solution of objective function
The specific solving process is as follows:
step k (k > 1):
σk=ψ(yk) (17)
Figure BDA0003150886440000066
Figure BDA0003150886440000071
wherein σkIs the result of the k-th cycle, ykIs a prediction vector, t, for accelerating the iterative processkIs a parameter that controls the rate of convergence.
Compared with the existing angular super-resolution method, the method can better recover the scene outline and the details by utilizing the two-dimensional space prior information. By selecting a proper neighborhood system, the MRF can model the structural features of the image; then, utilizing a regularization frame to absorb MRF prior information to obtain a target function; and finally, solving the objective function through a fast iterative threshold shrinkage algorithm, and realizing the super-resolution imaging of the RASR.
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
In one embodiment, the method is mainly verified by a simulation experiment method, and all the steps and conclusions are verified to be correct on Matlab 2018. The present invention will be described in further detail with reference to specific embodiments.
The method comprises the following steps: fig. 2 is a schematic view of the scanning radar imaging according to the embodiment, in which the radar transmits a chirp signal while scanning, and parameters required for simulation are shown in table 1.
TABLE 1 simulation parameter table for real aperture scanning radar imaging
Parameter(s) Numerical value
Platform velocity 100m/s
Pulse repetition frequency 2000Hz
Main lobe beamwidth 2.5°
Scanning speed of antenna 60°/s
Antenna scanning range ±10°
Minimum pitch 3km
Target scene as shown in fig. 3, the simulation scene is two X-shaped planar targets composed of 84 point targets, the scattering coefficients of which are both 1, the azimuth distribution is-3.7 ° to-0.7 ° and 3.7 ° to 0.7 °, and the distance distribution is 2980m to 3020 m. At a sampling rate frThe range Echo is sampled at 100MHz to obtain an Echo matrix, denoted Echo (τ, t), with a size Mr×NaAs shown in fig. 4.
Step two: distance direction pulse compression is carried out on the echo waves, and a distance direction pulse compression reference function refer (tau) ═ exp (i pi k) is constructed according to system parametersrτ2) Wherein k isr-1. mu.s.ltoreq.τ.ltoreq.1. mu.s at 20 MHz/. mu.s. Fourier transform is carried out on the echo in the step one along the distance direction, and the echo is subjected to Fourier transform with a conjugate refer of the distance direction reference function*(f) Multiplying, then inversely transforming the multiplied result to the time domain to complete the range-wise pulse compression of the simulated echo, and recording the data after range compression as sout(τ, t), the size of which remains unchanged.
Step three: performing migration correction on the data obtained in the step two, and constructing a distance frequency vector f ═ fr/2,fr/Mr,…,fr/2-fr/Mr]T,MrIs the distance vector number and time vector t ═ 0, Tscan/Na,2×Tscan/Na,…,Tscan]Tscan is the scan time; and then constructing a migration correction function correct (f, t) ═ exp (-i2f × dr (t)/c) according to the system parameters, wherein dr (t) ═ vt. Fourier transform is carried out on the echo along the distance direction, and the transform is recorded as sout(f, t) and multiplying with a migration correction function correct (f, t), then reversely converting the multiplied result into a time domain to finish the range-direction pulse compression of the simulation echo, and recording the data after the migration correction as sc(τ, t) as shown in fig. 5 and 6. The third step is to theoretically describe the migration correction in the time domain, and the third step is to perform the correction in the distance frequency domain because the third step is performed in the time domainCan be achieved by phase multiplication in the frequency domain.
Step four: discretizing the scanning antenna pattern according to the antenna pattern width theta, the pulse repetition frequency prf and the scanning speed omega
Figure BDA0003150886440000081
Obtaining the point number m of the discretized antenna directional diagram, then h (theta) is h0 h1...hm-1]From which a convolution matrix H is constructed
Figure BDA0003150886440000082
Step five: to scompressAnd (tau, t) taking the amplitude of the jth data (j is 1) which is the original echo s to be processed, constructing an MRF imaging target function according to an expression (13), wherein s is a vector obtained by taking a specific distance unit from the data after the migration correction in the third step, and H is a convolution matrix constructed in the fourth step.
Step six: the gradient of the objective function is obtained according to the expression (14), and epsilon is 10-8The regularization parameter λ is determined using an L-curve method. Obtaining a basic iterative formula of FIST solution by a fixed point solution strategy:
Figure BDA0003150886440000091
wherein sigmakAnd σk+1For the results of two adjacent iterations, γ is taken as 1 and δ is taken as 0.001.
Step seven: iteratively solving an objective function, and assigning the target echo s to be processed as an iteration initial value sigma0. Obtaining an iteration result sigma by using the formula (20)1(ii) a Then using the obtained sigma0And σ1Performing one-step prediction according to the formulas (18) and (19), and substituting the prediction into the formula (17) to obtain the next iteration result sigma2,tkThe initial assignment is 1. When the results of two iterations satisfy the following formula
||σ21||2<κ (21)
The result converges, i.e. a target solution is found, where κ is set to 0.025 based on empirical values. If not, then sigma is added2、σ1、tk+1Respectively assign to sigma1、σ0And tkAnd continuing to perform the iteration step seven until the iteration termination condition is met to obtain the target solution.
The MRF-based real aperture scanning radar super-resolution imaging processing is completed, and the imaging result is shown in fig. 7. As can be seen from the figure, the method provided by the invention can well realize the super-resolution imaging of the real aperture radar, and the recovery of the X shape is far better than the result of the comparison algorithm.
The foregoing illustrates and describes the principles, general features, and advantages of the present invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.

Claims (10)

1. A MRF-based real aperture scanning radar super-resolution imaging method is characterized by comprising the following steps:
introducing structural prior information represented by a hidden Markov random field into an objective function through a regularization frame;
and solving the optimization problem of the objective function by using a fast iterative threshold contraction method, and finally obtaining the super-resolution imaging result of the target.
2. The MRF-based real aperture scanning radar super-resolution imaging method according to claim 1, wherein the structural prior information represented by the hidden markov random field is introduced into the objective function through the regularization framework, and firstly, parameters of the imaging system are initialized, specifically:
let the platform moving speed be v, the antenna beam scanning speed omega, and the radar down-view angle be
Figure FDA0003150886430000017
The coordinates of the target point are denoted as P (r)0,θ0) T is azimuth time, target point distance history rP(x, y) is
Figure FDA0003150886430000011
The slope distance obtained by approximation of Taylor's expansion and first order term is
Figure FDA0003150886430000012
The expression of the point target echo is recorded as sP(τ,t):
Figure FDA0003150886430000013
Where τ is the distance-wise time variable, σ0For point target scattering intensity, h (t) for the dual-pass beam pattern of the antenna, t0Rect [. for the time when the beam center hits point P]Representing the distance time window, c is the speed of light, TrRepresenting the distance-time pulse width, f0Is the carrier frequency, KrIs the time chirp rate of the transmitted signal.
3. The MRF-based real aperture scanning radar super-resolution imaging method according to claim 2, wherein after the initialization of the imaging system parameters, the range-wise pulse compression is performed on the echo, specifically as follows:
for a scattering function of
Figure FDA0003150886430000014
For a surface target, the echo is the integral of the echoes of all point targets, wherein
Figure FDA0003150886430000015
Is the initial slope distance of the target,
Figure FDA0003150886430000016
is the azimuth coordinate of the target; because t is (theta-theta)a) ω and τ 2r/c, where θaFor the initial scan angle, the surface target echo can be obtained according to equations (2) and (3):
Figure FDA0003150886430000021
aiming at the acquired two-dimensional echo data s (r, theta), a matched filter is constructed by using radar emission signal parameters, the range direction is subjected to pulse compression, and the data s after range compression is obtainedout(r,θ);
Figure FDA0003150886430000022
4. The MRF-based real aperture scanning radar super-resolution imaging method according to claim 3, wherein linear migration correction is performed after range-wise pulse compression, specifically:
after compressing in the range direction, for the target plane
Figure FDA0003150886430000023
Up located on arc line
Figure FDA0003150886430000024
The center of the echo position of all point targets is overlapped on a straight line
Figure FDA0003150886430000025
The above step (1); the variable r can be subjected to linear scale transformation to eliminate distance walking generated by platform motion
Figure FDA0003150886430000026
5. The MRF-based real aperture scanning radar super-resolution imaging method according to claim 4, wherein after linear migration correction, an orientation convolution signal model is constructed, specifically:
writing an orientation signal as
s(θ)=σ(θ)*h(θ) (6)
Wherein is a convolution operation; considering the influence of noise, the azimuth signal expression is:
s(θ)=σ(θ)*h(θ)+n(θ) (7)
and (3) constructing an azimuth convolution signal vector model aiming at the model:
s=Hσ+n (8)
wherein s ═ s1,…,sN]TAnd σ ═ σ [ σ ]1,...,σN]TRespectively representing the echo data and the target scattering coefficient in a range gate, n ═ n1,...,nN]TN is the number of vectors, and H is the convolution matrix constructed by the antenna directional diagram H (theta).
6. The MRF-based real aperture scanning radar super-resolution imaging method according to claim 5, wherein after constructing the azimuth convolution signal model, constructing an imaging objective function, specifically:
for the problem of noise sensitivity in deconvolution operation, a regularization method is adopted to solve the problem, and a standard regularization equation for solving the problem is as follows:
Figure FDA0003150886430000031
wherein the content of the first and second substances,
Figure FDA0003150886430000032
for data fidelity terms, | Γ (σ) | purpleqThe method is characterized in that the method is a regular term, gamma (sigma) is a specific operation of the regular term, q is a norm of the regular term, and lambda is used for balancing the fidelity of measured data and the influence of the constraint of the regular term and can be determined by an L-curve method;
selecting Huber-Markov as regularization term
Figure FDA0003150886430000033
Where ρ isT(. cndot.) is a function of Huber,
Figure FDA0003150886430000034
is a coefficient vector of the cluster c, tau is a temperature parameter, Z is a normalized constant, superscript (·)tDenoted as a transpose operation; the roughness of the image is obtained from the second derivative:
Figure FDA0003150886430000035
the Huber function is defined as:
Figure FDA0003150886430000036
the threshold value T can punish the gray level change of the image;
substituting equations (10), (11), (12) into equation (9), the objective function of RASR super-resolution imaging can become:
Figure FDA0003150886430000041
7. the MRF-based real aperture scanning radar super-resolution imaging method according to claim 6, wherein a complete image objective function is constructed to obtain an objective function gradient, specifically as follows:
because the obtained objective function is a convex function, the optimization problem of the formula (13) is solved by adopting a fast iterative shrinkage/threshold algorithm, and the gradient of the objective function is obtained
Figure FDA0003150886430000042
Wherein, (.)mThe mth element of the vector in brackets, diag (·) is a diagonal matrix, ε is a constant introduced to account for the irreducibility due to the non-smoothness of the objective function; the one-time iterative threshold shrinkage operation process is as follows:
Figure FDA0003150886430000043
wherein the content of the first and second substances,
Figure FDA0003150886430000044
is a shrink operator, defined as:
Figure FDA0003150886430000045
8. the MRF-based real aperture scanning radar super-resolution imaging method according to claim 7, wherein after the gradient of the objective function is solved, the objective function is solved by fast iteration, and the specific solving process is as follows:
step k:
σk=ψ(yk) (17)
Figure FDA0003150886430000046
Figure FDA0003150886430000051
wherein k is>1,σkIs the result of the k-th cycle, ykIs a prediction vector, t, for accelerating the iterative processkIs a parameter that controls the rate of convergence.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the computer program implements the steps of the super-resolution imaging method according to any one of claims 1 to 8.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the super-resolution imaging method according to any one of claims 1 to 8.
CN202110765492.3A 2021-07-06 2021-07-06 MRF-based real aperture scanning radar super-resolution imaging method Active CN113640793B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110765492.3A CN113640793B (en) 2021-07-06 2021-07-06 MRF-based real aperture scanning radar super-resolution imaging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110765492.3A CN113640793B (en) 2021-07-06 2021-07-06 MRF-based real aperture scanning radar super-resolution imaging method

Publications (2)

Publication Number Publication Date
CN113640793A true CN113640793A (en) 2021-11-12
CN113640793B CN113640793B (en) 2023-12-26

Family

ID=78416780

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110765492.3A Active CN113640793B (en) 2021-07-06 2021-07-06 MRF-based real aperture scanning radar super-resolution imaging method

Country Status (1)

Country Link
CN (1) CN113640793B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117420553A (en) * 2023-12-13 2024-01-19 南京理工大学 Super-resolution imaging method for sea surface target scanning radar

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103487802A (en) * 2013-09-18 2014-01-01 电子科技大学 Scanning radar angle super-resolution imaging method
CN104732480A (en) * 2015-03-03 2015-06-24 中国空间技术研究院 Remote-sensed image super-resolution method based on non-local regularization model
US20160005399A1 (en) * 2012-07-23 2016-01-07 University Of Southern California Noise speed-ups in hidden markov models with applications to speech recognition
CN105550993A (en) * 2016-01-18 2016-05-04 中国空间技术研究院 Multiple transform domain based super-resolution reconstruction method
CN105699969A (en) * 2016-01-29 2016-06-22 电子科技大学 A maximum posterior estimated angle super-resolution imaging method based on generalized Gaussian constraints
CN110780274A (en) * 2019-11-04 2020-02-11 电子科技大学 Improved L for scanning radar 1Regularization azimuth super-resolution imaging method
CN112037136A (en) * 2020-09-18 2020-12-04 中国科学院国家天文台南京天文光学技术研究所 Super-resolution imaging method based on aperture modulation

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160005399A1 (en) * 2012-07-23 2016-01-07 University Of Southern California Noise speed-ups in hidden markov models with applications to speech recognition
CN103487802A (en) * 2013-09-18 2014-01-01 电子科技大学 Scanning radar angle super-resolution imaging method
CN104732480A (en) * 2015-03-03 2015-06-24 中国空间技术研究院 Remote-sensed image super-resolution method based on non-local regularization model
CN105550993A (en) * 2016-01-18 2016-05-04 中国空间技术研究院 Multiple transform domain based super-resolution reconstruction method
CN105699969A (en) * 2016-01-29 2016-06-22 电子科技大学 A maximum posterior estimated angle super-resolution imaging method based on generalized Gaussian constraints
CN110780274A (en) * 2019-11-04 2020-02-11 电子科技大学 Improved L for scanning radar 1Regularization azimuth super-resolution imaging method
CN112037136A (en) * 2020-09-18 2020-12-04 中国科学院国家天文台南京天文光学技术研究所 Super-resolution imaging method based on aperture modulation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
焦淑红等: "迭代收缩阈值雷达前视成像方法", 《航空学报》, vol. 35, no. 12, pages 3384 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117420553A (en) * 2023-12-13 2024-01-19 南京理工大学 Super-resolution imaging method for sea surface target scanning radar
CN117420553B (en) * 2023-12-13 2024-03-12 南京理工大学 Super-resolution imaging method for sea surface target scanning radar

Also Published As

Publication number Publication date
CN113640793B (en) 2023-12-26

Similar Documents

Publication Publication Date Title
CN102967859B (en) Forward-looking scanning radar imaging method
CN103487803B (en) Airborne scanning radar imaging method in iteration compression mode
CN110726992B (en) SA-ISAR self-focusing method based on structure sparsity and entropy joint constraint
CN109613532B (en) Airborne radar real-time Doppler beam sharpening super-resolution imaging method
CN107621635B (en) Forward-looking sea surface target angle super-resolution method
CN105137425B (en) The preceding visual angle ultra-resolution method of scanning radar based on Deconvolution principle
CN111551928B (en) Wall-penetrating radar imaging method based on low-rank sparse constraint of wall
CN107607945B (en) Scanning radar foresight imaging method based on spatial embedding mapping
Yang et al. Doppler beam sharpening imaging based on fast iterative adaptive approach
CN112147608A (en) Rapid Gaussian gridding non-uniform FFT through-wall imaging radar BP method
Mao et al. Target fast reconstruction of real aperture radar using data extrapolation-based parallel iterative adaptive approach
CN112198506A (en) Method, device and system for learning and imaging ultra-wideband through-wall radar and readable storage medium
Mao et al. Angular superresolution of real aperture radar using online detect-before-reconstruct framework
CN109188436B (en) Efficient bistatic SAR echo generation method suitable for any platform track
CN113640793B (en) MRF-based real aperture scanning radar super-resolution imaging method
CN105759264A (en) Micro-motion target defect echo high-resolution imaging method based on time-frequency dictionary
CN104182942A (en) Method for synthetic aperture radar (SAR) image orientation blurring inhibition
CN115902886A (en) ITR-DTV radar forward-looking super-resolution imaging method based on Renyi entropy
CN115453523A (en) Scanning radar sparse target batch processing super-resolution method
CN115480245A (en) Airborne radar forward-looking three-dimensional imaging technology based on two-dimensional super-resolution algorithm
CN115453536A (en) Forward-looking azimuth pitching two-dimensional super-resolution imaging method for motion platform
CN113030964B (en) Bistatic ISAR (inverse synthetic aperture radar) thin-aperture high-resolution imaging method based on complex Laplace prior
CN112859074B (en) Multi-band multi-view ISAR fusion imaging method
Zha et al. An iterative shrinkage deconvolution for angular superresolution imaging in forward-looking scanning radar
Wu et al. Forward-looking angular super-resolution for moving radar platform with complex deconvolution

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