CN102221696B - Sparse microwave imaging method - Google Patents

Abstract

The invention discloses a sparse microwave imaging method, relating to the information acquisition and processing technology. By utilizing sparsity of microwave imaging observation, through introducing sparse signal processing theory into microwave imaging technology, after signal processing and information extraction, the geometrical and physical characteristics of an observation object such as space position, scattering characteristic and motion characteristic and the like are obtained, wherein the sparse signal processing theory is through searching sparse representation field of the observation object, obtaining sparse microwave signal of the observation object in space, time, frequency spectrum or polarizing field sparse sampling. According to the sparse microwave imaging method in the invention, the bottleneck problems existed in the prior art such as difficult system realization, complex imaging processing method, redundant information and difficult feature extraction and the like based on Nyquist sampling theorem and classic digital signal processing theory are solved, thus microwave imaging system structure and imaging complexity is reduced.

Description

Sparse microwave imaging method

Technical field

The information of the present invention relates to is obtained and processing technology field, it is a kind of sparse microwave imaging method, utilize the sparse property of microwave imaging observation scene data set, when signal is launched in conjunction with observation scene sparse property, and utilize antenna to be lower than the data sampling of Nyquist's theorem when data are obtained, it is integrated to be embodied as picture and information extraction when signal is handled.

Background technology

Microwave Imaging Technique be electromagnetic wave with the microwave spectral coverage as detection means, utilize the microwave imaging sensor to obtain the technology of object being observed scattering signatures and relevant information.Compare with optical image technology, microwave imaging is not subjected to the restriction of sunshine and weather condition, has round-the-clock, round-the-clock ability to target or scene observation, has developed into the important means of remote sensing applications such as resources survey, environmental monitoring and disaster assessment.In over half a century in the past, microwave imaging system resolution develops into a centimetre magnitude by tens meters magnitudes; The imaging system develops into bunching type, scan-type, slip pack, two/multistation and the three-dimensional imaging that faces the future and moving-target imaging pattern by initial band imaging; Polarization mode is developed to complete polarization from single polarization; Application mode develops into DEM measurement, atural object parameter measurement, ocean wave parameter measurement and parameters of target motion measurement etc. by the qualitative decipher of single image.

Existing microwave imaging system is to be based upon radar to differentiate on theoretical and the nyquist sampling theorem basis.Radar is differentiated theoretical resolution=light velocity/(2 * signal bandwidth), and the sample frequency during recovery that nyquist sampling theorem is undistorted is greater than 2 times bandwidth.Radar is differentiated theoretical and nyquist sampling theorem is the theoretical foundation of the universality of microwave imaging, can not break through.It has determined scale and the complexity of microwave imaging system.Resolution is more high, the mapping bandwidth is more big, more multisystem is just more complicated for port number.

The core that existing microwave imaging is handled is to make up distance-Doppler's two dimension according to the microwave system observing matrix to be decoupled into the picture operator, thereby obtains the microwave imagery of object of observation.Obtaining distance to aspect the high resolving power, generally all adopt the pulse compression system; Obtaining the orientation to aspect the high resolving power, the main employing reduced the orientation to the method for antenna length.

It is one of most active branch of signal process field in this century that sparse signal is handled.The goal in research of this branch is to extract the least possible observation data from original signal, keeps contained information in the original signal simultaneously to greatest extent, and original signal is effectively approached and recovered.

1986, Santosa and Symes clearly proposed the concept of sparse signal the earliest, and the basic mathematic model that sparse signal is handled is linear regression model (LRM).Suitably select sparse observing matrix Ф, then the sparse signal transaction module can be described with following equation,

Y＝ФX+N

Wherein Y is the observed samples data, and N is noise, and X is vector to be estimated.In general, the non-zero entry number of X seldom that is to say that X has sparse characteristic.Simultaneously, the dimension of Y is much smaller than the dimension of X.Therefore, this is a underdetermined equation, and infinite a plurality of solution is arranged generally speaking.The effective ways of finding the solution above-mentioned " sparse linear recurrence " problem are to use low-order mode that common second order error is carried out regularization, namely to objective function || and y-Φ x|| ₂+ λ || x|| _lBe optimized, this optimization problem can be converted into following problem and find the solution:

min||x|| _l s.t.||y-Фx|| ₂≤ε，l＝0

Here || || _lThe l rank norm of expression variable, min () expression minimizes operation, s.t. expression makes and satisfies, ε is the threshold value that arranges, and the zero norm of x is represented the non-zero entry number of x, can be used for describing the sparse characteristic for the treatment of estimate vector, under existing theory of foundations of mathematics condition, finding the solution of zero norm is very difficult, is similar to l ∈ [0,2] rank norm usually.

2006, Donoho has systematically discussed the association between l 1 (being l=1) optimization and the sparse property, also used the concept of Compressed Sensing (being called for short CS) first, at treating characteristic that estimated signal can sparse sign in certain space, adopt specific dimensionality reduction compression sampling, utilizing optimization method to realize signal reconstruction, signals sampling, recovery and information extraction directly are based upon on the signal sparse characteristic Foundation of Representation.2006, Candes and Tao proved that the RIP condition is l 1 and l 0 optimization problem adequate condition of equal value, has disclosed the association between l 1 optimization and the sparse signal reconstruction.On their working foundation, the scholars that signal is handled the boundary have carried out research deeply and widely in the relevant issues of in recent years sparse signal being handled.At present, the numerical solution of l 1 optimization problem mainly contains two big class methods: a class is based on the method for iteration optimizing technology, as protruding optimization and linear programming etc.; The another kind of method that is based on greedy algorithm (Greedy Algorithm).

The achievement in research that the sparse signal treatment theory is recent and the breakthrough that obtains make it obtain comparatively successful application at optical imagery (as single pixel optics camera of Duarte, Davenport, Takhar etc.).In these were used, object of observation had tangible sparse characteristic, can obtain comparatively ideal reconstructed results by sparse base commonly used, adopted sample relation between sparse property and the sparse property of object being observed of statistical method approximate description.

Existing microwave imaging system is not considered following factor:

At first, existing microwave imaging system does not consider that microwave imaging has sparse characteristic.The sparse property of microwave imaging can be embodied directly in the transform domain of microwave imagery or image, and this is because the scene itself that is observed often has stronger correlativity.Microwave imagery is under the specific microwave observation condition, the scene echoes data are through the relevant synthetic sign of handling back electromagnetic scattering characteristic, as the different observation data collection that are observed scene (the microwave imaging raw data, through the part imaging processing data, image or at its transform domain) have a possibility of rarefaction sign.

Secondly, existing algorithm data to obtain with imaging processing, imaging processing and information extraction be separate between any two.Because the observation complicacy of scene and the diversity of microwave imaging system, existing microwave imaging theory are not still supported to conduct a research from the obtaining of microwave imaging data, processing and the integrated angle of information extraction systemicly.

The 3rd, the application of sparse signal treatment theory in microwave imaging at present all is from existing radar system and working system, be mainly used in addressing the problem: in real work, because situations such as condition is limited or be interfered, when the data of obtaining can not satisfy sampling constraint or segmental defect, by interpolation and the extrapolation to the default signal, recover to satisfy the default sampling of Nyquist's theorem; Recycling classical signals disposal route is embodied as picture and information extraction.Handle at the sparse signal that this meaning is introduced, still can't consider when the microwave imaging data are obtained that application demand carries out on purpose sparse sampling, can't solve fundamentally that application demand and microwave imaging system are realized difficulty, magnanimity information is redundant, signal is handled and characteristic information extracts contradiction between the difficulty.

List of references

[1]F.Santosa，and W.W.Symes，Linear inversion of band-limited reflectionseismograms，SIAM J.Sci.Statist Comput 7(1986)，1307-1330.

[2]D.L.Donoho，Compressed Sensing，IEEE Trans.Inform.Theory，52，no.4，(2006)，pp.1289-1306.

[3]D.L.Donoho，For most large underdetermined systems of linear equationsthe minimal l 1-norm solution is also the sparsest solution，Comm.PureAppl.Math.59，no.6(2006)，pp.797-829.

[4]E.J.Candès and T.Tao，Decoding by linear programming，IEEE Trans.Inform.Theory 51，4203-4215，2006.

Summary of the invention

The objective of the invention is to disclose a kind of sparse microwave imaging method, utilize the sparse property of microwave imaging, by reducing data transfer rate, make algorithm quick, easy, reduced the complexity that system realizes, system cost reduced, and be embodied as picture and feature extraction integrated.

For achieving the above object, technical solution of the present invention is:

A kind of sparse microwave imaging method, it is handled sparse signal and introduces microwave imaging, and organically combines the microwave imaging that forms;

By seeking the sparse sign territory of object being observed, in the space, time, frequency spectrum or polarizing field sparse sampling, obtain the sparse microwave signal of object being observed, through sparse microwave signal is handled and information extraction, obtain locus, scattering signatures and the kinetic characteristic of object being observed, by optimizing algorithm, restoration scenario target information;

According to the sparse microwave imaging method of the sparse property structure of microwave imaging observation data collection, the data volume of collection is lacked than the data volume based on the systematic sampling of nyquist sampling theorem, thereby reduces the complexity of microwave imaging system.

Described sparse microwave imaging method, its described time is sparse, is that the microwave imaging middle distance is sparse to the sampling time two dimension to sampling time, orientation.

Described sparse microwave imaging method, its described space is sparse, is and position of platform, the array element relevant sparse characteristic that distributes, and by the overall treatment to each passage echoed signal, the information that realizes being observed the zone is recovered and the extraction of target signature information.

Described sparse microwave imaging method, its described frequency spectrum is sparse, be that sparse spectrum signal is adopted in radar emission and reception, handle the recovery broadband signal by sparse signal, and super-resolution image and the multiple spectra characteristic information of object being observed are obtained in the sparse microwave imaging of the sparse microwave imaging of binding time or space.

Described sparse microwave imaging method, its described polarization is sparse, is to utilize polarization, the mixed polarization technology of condensing, and realizes that by sparse POLARIZATION CHANNEL the full polarimetric SAR data of target is obtained, to reduce the requirement to radar system transceiver channel number and pulse repetition rate.

Described sparse microwave imaging method, its described joint sparse is to utilize wherein both or both above sparse characteristic of space, time, frequency spectrum or polarization, realizes the joint sparse microwave imaging.

Described sparse microwave imaging method, it comprises step:

Step S1: the signal of emission is linear FM signal, or the random series signal, or bidimensional coded signal when empty, and the signal of emission is:

In the formula, t is the time, and a (t) is the amplitude that transmits, and f is the carrier frequency that transmits,

It is the phase place that transmits;

Step S2: signal is launched by the antenna of sparse microwave imaging radar system after amplifying by power amplifier;

Step S3: set up echo model according to the geometric relationship between the form that transmits and Texas tower and the target;

Step S4: the sample mode of sparse microwave imaging radar system is even low speed sampling, stochastic sampling or process pre-service sampling afterwards;

Step S5: set up sparse microwave observation equation, to observing the recovery of scene:

S51: for the scene X with obvious sparse characteristic, the observation equation of sparse microwave is:

Y＝ФX +N

Wherein, Y is echo data, and Ф is observing matrix, and N is noise; According to sparse signal treatment theory, l ₁Optimization can well provide K element of scene X intermediate value maximum

$\hat{X}=\arg \min {\left|\right|X\left|\right|}_{1}s.t.{\left|\right|Y-\mathrm{\ΦX}\left|\right|}_2\≤\mathrm{\ϵ}$

Wherein, || || _lThe l rank norm of expression variable, min () expression minimizes operation, and ε is to be to optimize the threshold value that convergence is set when existing for noise;

S52: for there not being obviously sparse scene X, must comprise space, time, frequency spectrum, polarization or various dimensions joint sparse characteristic according to its sparse base, the sparse map table of scene X is shown:

X＝Ψθ

Wherein Ψ is sparse transformation matrix, and θ is the expression of X under the Ψ transformation matrix, and therefore the observation equation of sparse microwave imaging radar is expressed as:

Y＝ФΨθ+N

According to l ₁Optimum theory, K maximum sparse base system number of value provided by following formula:

$\widehat{\mathrm{\&theta;}}=\arg \min {\left|\right|\mathrm{\&theta;}\left|\right|}_{1}s.t.{\left|\right|Y-\mathrm{\&Phi;\&Psi;\&theta;}\left|\right|}_2<\mathrm{\&epsiv;}$

Here, || || _lThe l rank norm of expression variable, min () expression minimizes operation, and s.t. represents to make satisfied; Wherein, ε is when existing for noise, for optimizing the thresholding that convergence is set, recovers after the sparse coefficient, and scene is expressed as:

$X=\mathrm{\&Psi;}\widehat{\mathrm{\&theta;}};$

S53: obtain among the scene information X positions of elements corresponding to the locus of object being observed according to above-mentioned steps, among the X value of element corresponding to scattering properties and the kinetic characteristic of object being observed, under the transform domain among the θ value of element corresponding to characteristic information.

Described sparse microwave imaging method, the sparse microwave imaging among its described step S5, it is relevant with the sparse property in the observation data set transformation territory of the sparse property of the observation data collection of observing scene or observation scene that observation scene sparse table is levied; Utilize observation data collection space, time, frequency, polarization or the associating various dimensions of observation scene sparse; Its transform domain relation is unit transformation, discrete cosine transform, wavelet transformation or Walsh transform.

Described sparse microwave imaging method, the observation data collection of its described observation scene comprises the transform domain data of microwave imaging raw data, the data through the part imaging processing, view data or above-mentioned data.

Described sparse microwave imaging method, the waveform that transmits among its described step S1 are to depend on the observation scene characteristics, this waveform be linear FM signal, random series signal, nonlinear frequency modulation signal or when empty the bidimensional coded signal one of them.

Described sparse microwave imaging method, data capture method among its described step S4 is sampled for being lower than Nyquist rate ADC, and there not to be the fuzzy target information of recovering, its sampling rate is to determine that according to the complexity of scene more strong its ADC sampling rate of the sparse property of scene is more low; This data capture method adopts stochastic sampling, uniform sampling or the pretreated sampling of process.

Described sparse microwave imaging method, observing matrix Ф among its described step S51 is to be determined by the waveform that transmits, antenna parameter and position, data capture method, operation wavelength, pulse repetition rate, operating distance, work visual angle, position of platform, platform kinetic characteristic.

Described sparse microwave imaging method, sparse microwave signal among its described step S5 is handled, be to adopt protruding optimization and linear programming, or based on the method for greedy algorithm, by the coefficient in transform domain of this signal processing with radar scattering characteristic or the scene objects radar scattering characteristic of restoration scenario target, thereby realization imaging processing and feature extraction are integrated.

Described sparse microwave imaging method, it utilizes the sparse property of observation scene target signature, sparse microwave imaging method observing matrix and described sparse microwave signal to handle for observing the extraction of scene target information and imaging integrated.

Described sparse microwave imaging method, it is used for carrying out moving object detection, utilizes the sparse property of moving target, sparse microwave imaging method observing matrix and described sparse microwave signal to handle.

Described sparse microwave imaging method, it is used for wide swath ocean target imaging, utilizes the sparse property of target on the ocean, sparse microwave imaging method observing matrix and described sparse microwave signal to handle.

Described sparse microwave imaging method, it is used for carrying out synthetic aperture radar image-forming, utilizes the sparse property of the observation data set transformation domain coefficient of observation scene, sparse microwave imaging method observing matrix and described sparse microwave signal to handle.

Described sparse microwave imaging method, it is used for carrying out inverse synthetic aperture radar imaging, utilizes the sparse property of aerial target, sparse microwave imaging method observing matrix and described sparse microwave signal to handle.

Described sparse microwave imaging method, it is used for carrying out circumferential synthetic aperture radar imaging (CSAR), utilizes the sparse property of the observation data set transformation domain coefficient of observation scene, sparse microwave imaging method observing matrix and described sparse microwave signal to handle.

Described sparse microwave imaging method, it is used for the hyperchannel radar imagery, comprise that interference synthetic aperture radar imaging, three-dimensional imaging, MIMO (Multiple-Input Multiple-Out-put) imaging, station-keeping radar satellite imagery, multistatic radar imaging, digital beam form, space-time adaptive is handled, utilize the sparse property of each passage observation data set transformation domain coefficient of observation scene, the correlativity of observing each passage observation data of scene, sparse microwave imaging method observing matrix and described sparse microwave signal to handle.

Sparse microwave imaging method of the present invention has reduced data transfer rate, makes algorithm quick, easy, has reduced system cost, and be embodied as the picture and feature extraction carry out simultaneously.

Description of drawings

Fig. 1 is a kind of sparse microwave imaging method synoptic diagram of the present invention;

Fig. 2 is a kind of sparse microwave imaging schematic flow sheet of the present invention.

Embodiment

Sparse microwave imaging refers to the sparse signal treatment theory is introduced microwave imaging and organically combined microwave imaging new theory, New System and the new method that forms, namely by seeking the sparse sign territory of object being observed, in the space, time, frequency spectrum or polarizing field sparse sampling obtain the sparse microwave signal of object being observed, handle and information extraction through signal, obtain how much of locus, scattering signatures and kinetic characteristics etc. and the physical features of object being observed.

Specific embodiment one: based on sparse sparse microwave imaging method (seeing also Fig. 1) of time

In sparse microwave imaging, exist to have obviously sparse scene X, that is to say among the X nonzero element seldom in other words most element all much larger than zero, for N element and obvious sparse scene are arranged:

||X|| ₀＜＜N

Wherein, || X|| _p, the p norm of expression X.

In sparse microwave imaging, though there is information redundancy in the not obvious sparse very strong correlativity that often has of the scene that is observed, therefore also has sparse property.We will observe scene be expressed as X, and there is sparse territory in it, that is to say that there is sparse transformation matrix Ψ in it, make

X＝Ψ·α

||α|| ₀＜＜N

It is the sparse sign of X that sparse coefficient vector α just is called, and it is nonzero element vector seldom.The conversion of this sparse transformation matrices correspondence can be unit transformation, discrete cosine transform, wavelet transformation or Walsh transform.Show and the microwave power image of decipher being generally used for, can directly adopt wavelet field or discrete cosine transform domain etc. as sparse territory.But because power diagram similarly is that complex pattern has been carried out result after the Nonlinear Processing (mould square), complex pattern is sought suitable Ψ, make complex pattern available sparse coefficient vector α under Ψ come sparse sign.In the data acquisition of microwave imaging, the different observation data collection of object of observation comprise microwave imaging raw data, the X as a result after the data of part imaging processing (to matched filtering handle as distance after data), image or their transform domain (can be unit transformation, discrete cosine transform, wavelet transformation or Walsh transform) data can be thought X carried out conversion process _f=f (X) just can be considered as these data sets to observe the function of scene, so also contain the intrinsic sparse characteristic of object of observation, but to the X in the different application _f, seek different sparse transformation matrix Ψ _f, to satisfy

X _f＝f(X)＝Ψ _f·α

||α|| ₀＜＜N

In a word, in actual treatment, can carry out conversion to the microwave power image, also can carry out conversion to complex pattern, also can carry out conversion to the microwave data that has passed through section processes.

Step S1: the signal of emission is bidimensional coded signal in the time of can being linear FM signal or random series signal or sky, and the signal of emission is

In the formula, a (t) is the amplitude that transmits, and f is the carrier frequency that transmits,

It is the phase place that transmits.

Step S2: signal is launched by the antenna of sparse microwave imaging radar system after amplifying by power amplifier.

Step S3: set up echo model according to the geometric relationship between the form that transmits and Texas tower and the target.Be the target of R for being positioned at distance by radar, the available echo shaping of system is

In the formula, C is the light velocity; ρ and φ are respectively amplitude and the phase place of target; 2R/C is that to transmit through target range be the two-way time of R; N (t) is observation noise.

So for the scene that has I target, after the echoed signal demodulation be

Step S4: the sample mode of sparse microwave imaging radar system can be the sampling after even low speed sampling, stochastic sampling and the process pre-service.Mode with even low rate sampling is the observation equation that example makes up sparse microwave imaging radar below.

For sparse microwave imaging radar, the signals sampling rate is determined according to the sparse property of signal.Degree of rarefication is defined as nonzero element shared proportion in the element of observation scene, for not obvious sparse scene, is defined as nonzero element and the ratio of observing situation elements in the sparse sign vector of scene.

p＝K/N

The sampling rate of writing sufficient nyquist sampling theorem all over is F _N, and the sampling rate of sparse microwave imaging radar is determined for following formula can be arranged:

F _C＝O(plog(1/p))·F _N

Its sampling interval is Δ T=1/F _C, then the point target sampling equation for unit amplitude can be expressed as

The form of being write as column vector is,

Because the radar return of diverse location can be regarded the time delay that transmits as and take advantage of amplitude and phase place in target, therefore adopt uniform sampling each column vector of observing matrix Ф be the time delay after the discretize of transmitting, if only fetch the limited echo of target scene region in the wave datum, adopt the observation equation of the sparse microwave imaging radar of uniform sampling to be

Wherein, the column vector number of Ф is the number behind the target scene uniform sampling.Target scene note X is

X＝[x ₁ x ₂ … x _M] ^T

Wherein, M is the length of X.x _i＝ρ _iexp(jφ _i)。But the scattering properties of representing the i scene objects.Therefore, the signal acquisition based on sparse microwave imaging radar can be expressed as

$y\left(n\right)=S\&CircleTimes;X=\underset{m=1}{\overset{M}{\mathrm{\&Sigma;}}}s(n-m)x\left(m\right)=\underset{m=1}{\overset{M}{\mathrm{\&Sigma;}}}s(n-m){\mathrm{\&rho;}}_M\exp \left(j{\mathrm{\&phi;}}_M\right)$

Wherein

The expression convolution, the form that is expressed as matrix is

Therefore, for scene X sparse on transform domain Ψ, its observation equation can be described as

Y＝ФX＝ФΨα

Step S5: to the recovery of observation scene

S51: for the scene X with obvious sparse characteristic, the observation equation of sparse microwave is

Y＝ФX+N

According to sparse signal treatment theory, l ₁Optimization can well provide K element of scene X intermediate value maximum

$\hat{X}=\arg \min {\left|\right|X\left|\right|}_1\mathrm{st}.{\left|\right|Y-\mathrm{\&Phi;X}\left|\right|}_2\&le;\mathrm{\&epsiv;}$

Wherein ε is when existing for noise, is to optimize the thresholding that convergence is set.Recover Can well represent scene.

S52: for there not being obviously sparse scene X, must comprise space, time, frequency spectrum, polarization or various dimensions joint sparse characteristic according to its sparse base, the sparse conversion of scene X can be expressed as

X＝Ψθ

The observation equation of sparse microwave imaging radar can be expressed as

Y＝ФΨθ+N

According to l ₁Optimum theory, K maximum sparse base system number of value can be provided by following formula:

$\widehat{\mathrm{\&theta;}}=\arg \min {\left|\right|\mathrm{\&theta;}\left|\right|}_{1}s.t.{\left|\right|Y-\mathrm{\&Phi;\&Psi;\&theta;}\left|\right|}_2\&le;\mathrm{\&epsiv;}$

Wherein ε is when existing for noise, and s.t. represents to make and satisfies, and is to optimize the thresholding that convergence is set.Recover after the sparse coefficient, scene can be expressed as

$X=\mathrm{\&Psi;}\widehat{\mathrm{\&theta;}}$

Specific embodiment two: based on the sparse microwave imaging of time and space joint sparse

Step 1-step 4: adopt the echo data step S5 that obtains a plurality of passages with embodiment one identical method, can obtain the echo samples formula according to step 1-step 4

Y＝ФX+N

Following formula is described is the observed quantity obtain manner of time during sparse sampling, below sparse and sparse signals collecting mode and the disposal route in space of analysis joint time.Consider that sparse microwave imaging system has I sampled aperture (space).The data of each sampled aperture are obtained equation

T _i＝Ф _iX+N _i i＝1，…，I

Wherein, Y _i, Ф _iAnd N _iBe respectively observed quantity, observing matrix and the observation noise in i aperture.All sampled points can be expressed as a new vector

$Y^{\&prime;}=\left[\begin{array}{c}{Y}_1\\ Y_2\\ .\\ .\\ .\\ T_I\end{array}\right]=\left[\begin{array}{c}{\mathrm{\&Phi;}}_1\\ {\mathrm{\&Phi;}}_2\\ .\\ .\\ .\\ {\mathrm{\&Phi;}}_I\end{array}\right]X+\left[\begin{array}{c}{N}_1\\ N_2\\ .\\ .\\ .\\ N_I\end{array}\right]$

The space is sparse to be referred to for being some apertures of sparse distribution on the space to the scene observed samples.Sampled data after the sparse sampling of space can be expressed as

Wherein, Δ is space sparse sampling matrix, is used for describing space sparse sampling problem.Δ _iBe and Ф _iEqual-sized matrix, it is unit matrix or null matrix, it is unit matrix when i channel sample, is null matrix when not participating in sampling.If order

Then the data of room and time joint sparse are obtained the step that can be expressed as among the S5.

Y＝ФX+N

Step S6: to the recovery of observation scene.The implementation method of this step is identical with step S5 among the embodiment one.

Claims (14)

1. a sparse microwave imaging method is characterized in that, is sparse signal to be handled introduce microwave imaging and organically combine the microwave imaging that forms;

Sparse sign territory by object being observed, in the space, time, frequency spectrum or polarizing field sparse sampling, obtain the sparse microwave signal of object being observed, through sparse microwave signal is handled and information extraction, obtain locus, scattering signatures and the kinetic characteristic of object being observed, by optimizing algorithm, restoration scenario target information;

According to the described sparse microwave imaging method of the sparse property structure of microwave imaging observation data, the data volume of collection is lacked than the data volume based on the systematic sampling of nyquist sampling theorem, thereby reduces the complexity of microwave imaging system;

Described sparse microwave imaging method comprises step:

Step S1: the signal of emission is linear FM signal, or the random series signal, or bidimensional coded signal when empty;

Step S2: signal is launched by the antenna of sparse microwave imaging radar system after amplifying by power amplifier;

Step S3: set up echo model according to the geometric relationship between the form that transmits and Texas tower and the target;

Step S4: the sample mode of sparse microwave imaging radar system is even low speed sampling, stochastic sampling or process pre-service sampling afterwards;

Step S5: set up sparse microwave observation equation, realize the recovery to the observation scene:

S51: for the scene X with obvious sparse characteristic, the observation equation of sparse microwave imaging radar system is:

Y＝ΦX+N

Wherein, Y is echo data, and Φ is the observing matrix that scene has sparse microwave imaging radar system under the condition of obvious sparse characteristic, and N is noise; According to sparse signal treatment theory, l ₁Optimization can well provide K element of scene X intermediate value maximum:

$\hat{X}=\arg \min {\left|\right|X\left|\right|}_{1}s.t.{\left|\right|Y-\mathrm{\&Phi;X}\left|\right|}_2\&le;\mathrm{\&epsiv;}$

Wherein, || || _lThe l rank norm of expression variable, min () expression minimizes operation, and s.t. represents to make and satisfies, and ε is to be to optimize the threshold value that convergence is set when existing for noise;

S52: for the scene X that does not have obvious sparse characteristic, must comprise space, time, frequency spectrum, polarization or various dimensions joint sparse characteristic according to its sparse base, the sparse map table that does not have the scene X of obvious sparse characteristic is shown:

X＝Ψθ

Wherein Ψ is sparse transformation matrix, scene X the expression under sparse transformation matrix Ψ of θ for not having obvious sparse characteristic, and therefore the observation equation of sparse microwave imaging radar system is expressed as:

Y＝ΦΨθ+N

According to l ₁Optimum theory, K maximum sparse base system number of value provided by following formula:

$\widehat{\mathrm{\&theta;}}=\arg \min {\left|\right|\mathrm{\&theta;}\left|\right|}_{1}s.t.{\left|\right|Y-\mathrm{\&Phi;\&Psi;\&theta;}\left|\right|}_2<\mathrm{\&epsiv;}$

Here, || || _lThe l rank norm of expression variable, min () expression minimizes operation, and s.t. represents to make satisfied; Wherein, ε is when existing for noise, for optimizing the thresholding that convergence is set, recovers after the sparse base system number, and the scene that does not have obvious sparse characteristic is expressed as:

$X=\mathrm{\&Psi;}\widehat{\mathrm{\&theta;}};$

S53: positions of elements is corresponding to the locus of object being observed among the scene X that obtains according to above-mentioned steps S51 or step S52, the value of element is corresponding to scattering properties and the kinetic characteristic of object being observed among the scene X, under the transform domain among the θ value of element corresponding to characteristic information.

2. sparse microwave imaging method as claimed in claim 1, it is characterized in that, the sparse microwave observation equation set up described in the described step S5 realizes in the recovery to the observation scene, and it is relevant with the sparse property in the observation data set transformation territory of the sparse property of the observation data collection of observing scene or observation scene that observation scene sparse table is levied; Utilize the sparse property of observation data collection space, time, frequency, polarization or the associating various dimensions of observation scene; The transform domain relation of observation data collection is unit transformation, discrete cosine transform, wavelet transformation or Walsh transform.

3. sparse microwave imaging method as claimed in claim 2, it is characterized in that, the observation data collection of described observation scene comprises microwave imaging raw data, the data through the part imaging processing, view data or above-mentioned microwave imaging raw data, through the data of part imaging processing, the transform domain data of view data.

4. sparse microwave imaging method according to claim 1, it is characterized in that, the waveform that transmits among the described step S1 depends on the observation scene characteristics, this waveform be linear FM signal, random series signal, nonlinear frequency modulation signal or when empty the bidimensional coded signal one of them.

5. sparse microwave imaging method according to claim 1, it is characterized in that, the method of sampling among the described step S4, for being lower than Nyquist rate ADC sampling, and there is not fuzzy recovery target information, its sampling rate is to determine that according to the complexity of scene more strong its ADC sampling rate of the sparse property of scene is more low; This method of sampling adopts stochastic sampling, uniform sampling or the pretreated sampling of process.

6. sparse microwave imaging method according to claim 1, it is characterized in that, observing matrix Φ among the described step S51 is to be determined by the waveform that transmits, antenna parameter and position, data capture method, operation wavelength, pulse repetition rate, operating distance, work visual angle, Texas tower position, Texas tower kinetic characteristic.

7. sparse microwave imaging method as claimed in claim 1, it is characterized in that, the sparse microwave observation equation of setting up among the described step S5 is realized observing the recovery of scene, be to adopt protruding optimization and linear programming, or based on the method for greedy algorithm, handle to recover the coefficient in transform domain of the scattering properties of the scattering properties of object being observed or object being observed by this signal, thereby realize that imaging processing and feature extraction are integrated.

8. sparse microwave imaging method according to claim 1, it is characterized in that, be used for observing the extraction of scene target information and imaging integrated, utilize the observing matrix of the sparse property of observation scene target signature, sparse microwave imaging radar system and described sparse microwave signal to handle.

9. sparse microwave imaging method according to claim 1 is characterized in that, is used for carrying out moving object detection, utilizes the observing matrix of the sparse property of moving target, sparse microwave imaging radar system and described sparse microwave signal to handle.

10. sparse microwave imaging method according to claim 1 is characterized in that, is used for wide swath ocean target imaging, utilizes the observing matrix of the sparse property of observed object on the ocean, sparse microwave imaging radar system and described sparse microwave signal to handle.

11. sparse microwave imaging method according to claim 1, it is characterized in that, be used for carrying out synthetic aperture radar image-forming, utilize the observing matrix of the sparse property of the observation data set transformation domain coefficient of observation scene, sparse microwave imaging radar system and described sparse microwave signal to handle.

12. sparse microwave imaging method according to claim 1 is characterized in that, is used for carrying out inverse synthetic aperture radar imaging, utilizes the observing matrix of the sparse property of aerial target, sparse microwave imaging radar system and described sparse microwave signal to handle.

13. sparse microwave imaging method according to claim 1, it is characterized in that, be used for carrying out the circumferential synthetic aperture radar imaging, utilize the observing matrix of the sparse property of the observation data set transformation domain coefficient of observation scene, sparse microwave imaging radar system and described sparse microwave signal to handle.

14. sparse microwave imaging method according to claim 1, it is characterized in that, be used for the hyperchannel radar imagery, comprise that interference synthetic aperture radar imaging, three-dimensional imaging, MIMO (Multiple-Input Multiple-Out-put) imaging, station-keeping radar satellite imagery, multistatic radar imaging, digital beam form, space-time adaptive is handled, utilize the observing matrix of the sparse property of each passage observation data set transformation domain coefficient of observation scene, the correlativity of observing each passage observation data of scene, sparse microwave imaging radar system and described sparse microwave signal to handle.

Images