CN107102326A - Quick relevance imaging method and device with data fusion is split based on band - Google Patents
Quick relevance imaging method and device with data fusion is split based on band Download PDFInfo
- Publication number
- CN107102326A CN107102326A CN201710385506.2A CN201710385506A CN107102326A CN 107102326 A CN107102326 A CN 107102326A CN 201710385506 A CN201710385506 A CN 201710385506A CN 107102326 A CN107102326 A CN 107102326A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mover
- band
- echo
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
The quick relevance imaging method split based on band with data fusion that the present invention is provided, orthogonal, the independent Stochastic narrow pulse signal of each antenna radiation unit synchronized transmissions space-time of radiating curtain, multiple bands are divided into by the imaging region of the tiltedly lower apparent direction of aerial radiation array；The space-time random radiation formed in each band interacts to form the pulse echo of each band with corresponding observed object；Each pulse echo is divided into multiple echo range gates and echo subphylum, is sampled in the specific echo subphylum in each echo range gate, obtains multigroup echo samples value；Using the space-time random radiation of each band and corresponding echo samples value, the relevance imaging equation of each band is constructed and solved, the imaging results of each banded zone are obtained；And select image interfusion method to seamlessly transit adjacent ribbons, obtain the overall inversion result of imaging region.The relevance imaging equation scale of present invention single band per treatment is smaller, realizes fast imaging.
Description
Technical field
The present invention relates to Radar Technology field, more particularly to it is a kind of split based on band and data fusion it is quick
Relevance imaging method and device.
Background technology
Radar imagery (Radar Imaging) technology is a leap in radar development history, and it is initial that it has expanded radar
Detection (Detection) and ranging (Ranging) function so that radar can using obtain electromagnetic scattering information, obtain
To the panoramic radar image of scene.
Synthetic aperture radar (SAR) is although with higher azimuth resolution, and its revisiting period is long, when needing long
Between it is continuous stare under the scene of observation and imaging, there is inevitable congenital inferior position.Traditional real aperture radar staring imaging,
Its angular resolution is limited by actual antennas array aperture, limits its application in practice.
Microwave stares relevance imaging because with surmounting the excellent of real aperture radar imaging resolution limit and fast imaging
Point, faster development is achieved in recent years.The core of the imaging method is that construction one has space-time bidimensional stochastic behaviour
Radiation field irradiation target area, make the echo of the time-varying of target scattering independence positioned at diverse location, finally to echo and pre-
The radiation field put is associated processing and obtains inversion chart picture.
When microwave is stared into relevance imaging applied to large scene imaging, the area of imaging region is kilometer or tens kilometers
Magnitude, resolution requirement is generally meter level or decimeter grade, therefore, it is necessary to carry out fine and closely woven during to imaging region progress sliding-model control
Mesh generation.Existing relevance imaging method for reconstructing is the radiation field at each mesh point of two-dimensional scene and target back scattering system
Number is all concatenated into one-dimensional vector, and now observing matrix and target information size to be reconstructed are very big, rebuild large-sized mesh
Mark takes very much, in addition, large scale observing matrix needs to occupy larger memory space, also has higher to calculator memory capacity
It is required that.
Due to needing longer reconstruction time and larger memory headroom, so high-quality using existing relevance imaging method
It is extremely difficult that amount ground rebuilds larger image scene.
The content of the invention
In consideration of it, the invention provides a kind of quick relevance imaging method and dress split based on band with data fusion
Put, solve microwave and stare the problem of longer reconstruction time and larger content space are needed in relevance imaging technology.
In order to realize foregoing invention purpose, the embodiments of the invention provide following technical scheme：
A kind of quick relevance imaging method split based on band with data fusion, including：
Control the Stochastic narrow pulse letter that each antenna radiation unit synchronized transmissions space-time of aerial radiation array is orthogonal, independent
Number；
By the imaging region of the tiltedly lower apparent direction of the aerial radiation array in distance to being divided into multiple bands, the day
The space-time random radiation that beta radiation array is formed in each band interacts with the observed object in respective strap
Form the pulse echo of each band；
Each pulse echo is divided into multiple echo range gates, each echo range gate is divided into multiple echoes
Sampled in subphylum, the specific echo subphylum in each echo range gate, obtain corresponding with band each described
Multigroup echo samples value；
Using the space-time random radiation of each band and corresponding multigroup echo samples value, construct and solve each
The relevance imaging equation of the band, obtains the imaging results of each band；
Imaging results to band each described are spliced, and obtain the overall imaging results of the imaging region.
It is preferred that, each antenna radiation unit is uniformly distributed but is not limited to be uniformly distributed in the aerial radiation array, institute
The position for stating antenna radiation unit isN=1 ..., the summation that N, N are the antenna radiation unit quantity, single-geophone receiver machine
Position beEach described antenna radiation unit synchronized transmissions random frequency hopping but the narrow pulse signal for being not limited to random frequency hopping, the
The transmission signal of the n antenna radiation units is expressed as：
Wherein, fnlThe tranmitting frequency of l-th of pulse of antenna radiation unit described in n-th, l=1 ..., L, L are total
Pulse number, τ is burst pulse duration, TpFor pulse spacing, Tp＞ Tmin+ τ+Δ τ, wherein, Δ τ is the exhibition of burst pulse echo
Wide time, Δ τ=2 (Rmax-Rmin)/c=Tmax-Tmin≈ 2Wsin θ/c, Rmax、RminRespectively described aerial radiation array is to institute
State the distance of imaging region proximally and distally, c=3 × 108M/s, W be imaging region distance to width, θ is the aerial radiation
Oblique downwards angle of visibility of the array relative to the imaging region.
It is preferred that, the imaging region of the tiltedly lower apparent direction by the aerial radiation array is more apart from being divided into upwards
Individual band, including：
According to the range gate feature of short-pulse radar, the imaging region of the tiltedly lower apparent direction of the aerial radiation array is existed
Distance is divided into multiple different bands upwards, and the strip width of duration τ pulse covering is w=c τ/2sin θ, band
Number is K=[W/w]=[2Wsin θ/c τ], and [x] is round up function, c=3 × 108m/s；
Spatial discretization processing is carried out to the imaging region, each band M grid, M=P × Q is divided into, P is
The imaging region orientation resolution cell number, Q is the band distance to resolution cell number, ρa、ρrRespectively described imaging region
Orientation and distance are to the grid spacing of discretization, and the position of j-th of grid dot center described in k-th in band is designated asIt is described
The target backscattering coefficient of position isK-th of band of observation area is expressed as
It is preferred that, space-time random radiation and respective strap that the aerial radiation array is formed in each band
Interior observed object interacts to form the pulse echo of each band, is specially：
Each institute of the time-varying random field for the random signal formation that each antenna radiation unit is produced in wave cover
State in band and be overlapped, it is any in kth (k=1,2 ..., K) individual described band during synthesis, empty bidimensional random radiation
A bitThe in-field at place is：
The in-field interacts the scattered field to be formed by the second pass in free space with the observed object
It is multicast to up to receiver, is for locusThe pulse echo of k-th of band that receives of receiver be：
It is defined into the amendment radiation field up to the receiver, it is any one in kth (k=1,2 ..., K) individual described band
PointThe amendment radiation field at place can be expressed as：
The scatter echo of k-th of band has following relation with amendment radiation field：
It is preferred that, it is described that each pulse echo is divided into multiple echo range gates, each echo range gate is drawn
It is divided into multiple echo subphylums, is sampled, obtained and each institute in the specific echo subphylum in each echo range gate
The corresponding multigroup echo samples value of band is stated, including：
By echo arrival point (l-1) T of l (l=2,3 ..., L) the individual Stochastic narrow pulse signalP+(τk,j)min
(l-1) T is moved to left on time shaftP, so as to which each burst pulse echo is alignd, wherein,k
=1,2 ..., K, j=1,2 ..., M；
Each burst pulse echo after alignment is divided into multiple echo range gates, each echo range gate is divided
For multiple echo subphylums, each band distance to discrete resolution cell number to be differentiated be Q, i.e. each band
It is ρ containing Q widthrSub- band, the corresponding echo subphylum width of each sub- band is the ρ of δ=2rSin θ/c, Mei Gesuo
Echo range gate width is stated for τ；
Is sampled in the specific echo subphylum in each echo range gate after alignment, obtain with described in each
The corresponding multigroup echo samples value of band.
It is preferred that, sampled, obtained in the specific echo subphylum in each echo range gate after alignment
To the multigroup echo samples value corresponding with band each described, any one in following two sampling policies is met：
Sampling policy one：When being not present overlapping between adjacent ribbons, each burst pulse echo after alignment it is every
Sampled in last echo subphylum in individual echo range gate, obtain the multigroup echo corresponding with band each described
Sampling instant t in sampled value, k-th of range gate of burst pulse echo described in l-thl,kMeet τ * k- δ ＜ tl,k＜ τ * k；
Sampling policy two：When adjacent ribbons distance is to overlapping q (q=1,2 ..., Q-1) individual sub- band, described in l-th
K-th of sampling instant t of burst pulse echol,kMeet tl,k=tl,1+ Δ t × (k-1), wherein, tl,1The burst pulse described in l-th
Sampling instant in first range gate of echo, τ-δ ＜ tl,1<τ, Δ t are sampling interval, Δ t=τ-δ × q.
It is preferred that, it is described to utilize the space-time random radiation of each band and corresponding multigroup echo samples value, structure
The relevance imaging equation of each band is made and solved, the imaging results of each band are obtained, including：
When being carried out to the reception echo of each band, after empty discretization, the echo samples of each band are obtained
Value matrix
When being carried out to the space-time random radiation of each band, after empty discretization, obtain the sight of each band
Survey matrix
According to the observing matrix of each bandAnd corresponding echo samples value matrixBuild each described
The relevance imaging equation of band：I.e.
Relevance imaging equation based on default relevance imaging algorithm and each band, is calculated in each band
The inverting value of target scattering coefficientAnd by the inverting value of the target scattering coefficientAs the imaging results of respective strap,Wherein,For the operator of the default relevance imaging algorithm.
Optionally, it is described to each when methods described, which is met, is not present overlapping between the sampling policy one, adjacent ribbons
The imaging results of the band are spliced, and obtain the overall imaging results of the imaging region, are specially：
Imaging results directly to band each described are spliced, and obtain the overall imaging results of the imaging region.
Optionally, when methods described meets the sampling policy two, adjacent ribbons distance is to overlapping1
It is individual) sub- band when, the imaging results to band each described are spliced, and obtain the overall imaging knot of the imaging region
Really, including：
Obtain the overlapping region between every group of adjacent ribbons in the imaging region；
It is without any processing to the imaging results outside overlapping region in the imaging region directly to retain, to every group of adjacent strip
Imaging results in two adjacent ribbons in band in overlapping region are weighted summation process, obtain in corresponding overlapping region
Inversion result, and using the inversion result as the imaging results in the overlapping region, be specially：
Wherein, S1And S2Two respectively adjacent bands,WithRespectively two neighboring strips S1And S2
PositionThe object inversion result at place,For in overlapping regionThe object inversion result at place, w1And w2It is respectively two neighboring
Strips S1And S2The corresponding weight coefficient in overlapping region in inversion chart picture, and meet w1+w2=1, w1,w2∈(0,1)；
Imaging results outside overlapping region in the imaging region and the imaging results each described in overlapping region are entered
Row splicing, obtains the overall imaging results of the imaging region.
A kind of quick relevance imaging device split based on band with data fusion, including：
Control unit, for controlling each antenna radiation unit synchronized transmissions space-time of aerial radiation array orthogonal, independent
Stochastic narrow pulse signal；
Cutting unit, for the imaging region of the tiltedly lower apparent direction of the aerial radiation array is more to being divided into distance
Individual band, space-time random radiation and the observation in respective strap that the aerial radiation array is formed in each band
Objectives interation forms the pulse echo of each band；
Sampling unit, for each pulse echo to be divided into multiple echo range gates, by each echo range gate
Multiple echo subphylums are divided into, is sampled, obtained and each in the specific echo subphylum in each echo range gate
The corresponding multigroup echo samples value of the band；
Structural unit, space-time random radiation and corresponding multigroup echo samples value for utilizing each band,
The relevance imaging equation of each band is constructed and solved, the imaging results of each band are obtained；
Concatenation unit, splices for the imaging results to band each described, obtains the entirety of the imaging region
Imaging results.
Relative to prior art, beneficial effects of the present invention are as follows：
The quick relevance imaging method and device based on band segmentation and data fusion that the present invention is provided, radiating curtain
Orthogonal, the independent Stochastic narrow pulse signal of each antenna synchronized transmissions space-time, by the imaging of the tiltedly lower apparent direction of aerial radiation array
Region is divided into the band that multiple width are range gate upwards in distance；What aerial radiation array was formed in each band
Space-time random radiation interacts to form the pulse echo of each band with the observed object in respective strap；By each pulse
Echo is divided in multiple range gates and echo subphylum, the specific echo subphylum in each range gate and sampled, and is obtained and each institute
State the corresponding multigroup echo samples value of band；Utilize the space-time random radiation of each band and corresponding multigroup echo samples
Value, constructs and solves the relevance imaging equation of each band, obtain the imaging results of each banded zone；In strips mosaic, selection
Image interfusion method seamlessly transits adjacent ribbons, obtains the overall inversion result of imaging region.The present invention provide based on bar
Relevance imaging method with segmentation and data fusion, because the relevance imaging equation scale of single band per treatment is smaller, because
This can realize fast imaging.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the accompanying drawing used required in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
The embodiment of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can also basis
The accompanying drawing of offer obtains other accompanying drawings.
Fig. 1 is a kind of quick relevance imaging method stream split based on band with data fusion disclosed in the embodiment of the present invention
Cheng Tu；
Fig. 2 is that the alignment of burst pulse echo and range gate divide schematic diagram disclosed in the embodiment of the present invention；
Fig. 3 is a kind of quick relevance imaging method split based on band with data fusion disclosed in the embodiment of the present invention
Flow chart；
Fig. 4 is strips mosaic method schematic diagram disclosed in the embodiment of the present invention；
Fig. 5 be disclosed in the embodiment of the present invention it is a kind of based on band split with the quick relevance imaging method of data fusion into
As schematic diagram of a scenario；
Fig. 6 is simulation objectives model schematic disclosed in the embodiment of the present invention；
Fig. 7 is imaging time disclosed in the embodiment of the present invention with the variation relation schematic diagram for dividing band number；
Fig. 8 is that object module recovers image schematic diagram when not dividing band disclosed in the embodiment of the present invention；
Fig. 9 recovers image schematic diagram when being bar interband no overlap when the disclosed division band number of the embodiment of the present invention is 8；
Figure 10 is that pulse width τ=70ns, bar interband disclosed in the embodiment of the present invention recover image schematic diagram when having overlapping；
Figure 11 is a kind of quick relevance imaging device split based on band with data fusion disclosed in the embodiment of the present invention
Structural representation.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, the technical scheme in the embodiment of the present invention is carried out clear, complete
Site preparation is described, it is clear that described embodiment is only a part of embodiment of the invention, rather than whole embodiments.It is based on
Embodiment in the present invention, it is every other that those of ordinary skill in the art are obtained under the premise of creative work is not made
Embodiment, belongs to the scope of protection of the invention.
Referring to Fig. 1, present embodiment discloses a kind of quick relevance imaging method split based on band with data fusion,
Specifically include following steps：
S101：Control orthogonal, independent random narrow of each antenna radiation unit synchronized transmissions space-time of aerial radiation array
Pulse signal；
It should be noted that each antenna radiation unit is uniformly distributed but is not limited to uniform divide in the aerial radiation array
Cloth, the position of the antenna radiation unit isN=1 ..., the summation that N, N are the antenna radiation unit quantity, single-point connect
The position of receipts machine isEach described antenna radiation unit synchronized transmissions random frequency hopping but the burst pulse letter for being not limited to random frequency hopping
Number, the transmission signal of antenna radiation unit described in n-th is expressed as：
Wherein, fnlThe tranmitting frequency of l-th of pulse of antenna radiation unit described in n-th, l=1 ..., L, L are total
Pulse number, τ is burst pulse duration, TpFor pulse spacing, Tp＞ Tmin+ τ+Δ τ, wherein, Δ τ is the exhibition of burst pulse echo
Wide time, Δ τ=2 (Rmax-Rmin)/c=Tmax-Tmin≈ 2Wsin θ/c, Rmax、RminRespectively described aerial radiation array is to institute
State the distance of imaging region proximally and distally, c=3 × 108M/s, W be imaging region distance to width, θ is the aerial radiation
Oblique downwards angle of visibility of the array relative to the imaging region.
S102：The imaging region of the tiltedly lower apparent direction of the aerial radiation array is divided into multiple upwards in distance
Band, space-time random radiation and the observed object in respective strap that the aerial radiation array is formed in each band
Interaction forms the pulse echo of each band；
Specifically, according to the range gate feature of short-pulse radar, by the tiltedly lower apparent direction of the aerial radiation array into
As region in distance to multiple different bands are divided into, the strip width of duration τ pulse covering is w=c τ/2sin
θ, band number is K=[W/w]=[2Wsin θ/c τ], and [x] is round up function, c=3 × 108m/s；
Spatial discretization processing is carried out to the imaging region, each band M grid, M=P × Q is divided into, P is
The imaging region orientation resolution cell number, Q is the band distance to resolution cell number, ρa、ρrRespectively described imaging area
Domain orientation and distance are to the grid spacing of discretization, and the position of j-th of grid dot center described in k-th in band is designated as
The target backscattering coefficient of the position isK-th of band of observation area is expressed as
When mesh generation enough to it is small when, can be with the target scattering coefficient of grid element center position come approximate whole grid.
Each band of the time-varying random field that the random signal that N number of antenna radiation unit is produced is formed in far field in wave cover
Inside it is overlapped, during synthesis, empty bidimensional random radiation, any point in kth (k=1,2 ..., K) individual bandPlace
In-field is：
Wherein,For the scalar Green's function of free space, it is assumed that transmitting antenna is preferable point source, propagate empty
Between be vacuum,It can be expressed as：
Wherein, δ is that Dick draws function.
The second pass that the scattered field of in-field and objectives interation formation passes through in free space is multicast to up to receiver,
It is for locusReceiver, the echo of its k-th of band received is：
Radiation field in view of the aerial radiation array emitter by Radiative antenna elements to target arrives receiver again
After round trip free-space propagation, it is defined into the amendment radiation field up to the receiver, kth (k=1,2 ..., K) individual band
Any pointThe amendment radiation field at place can be expressed as：
Revised space-time random radiation in quick relevance imaging method with data fusion, each band is split based on band
Field and objectives interation, form the scatter echo of each band, and scatter echo and the amendment radiation field of k-th of band have such as
Lower relation：
S103：Each pulse echo is divided into multiple echo range gates, each echo range gate is divided into many
Sampled, obtained and band each described in individual echo subphylum, the specific echo subphylum in each echo range gate
Corresponding multigroup echo samples value；
Different from traditional range gate be distance to minimum resolution cell, in each echo range gate only comprising distance to
The information of one resolution cell, splits the quick relevance imaging method with data fusion disclosed in the present embodiment based on band
In, each band is Q to discrete resolution cell number to be differentiated in distance, i.e., each band contains Q width for ρrSub- bar
Band, the corresponding echo subphylum width of sub- band is the ρ of δ=2rSin θ/c, in this case, the information that each echo samples value is included is
The coupling of resolution cell information in Q sub- bands.As can be seen from Figure 2 during Q=3 in each echo range gate distance to echo information
In each echo range gate comprising 3 width it is the ρ of δ=2 in situation, figurerSin θ/c echo subphylum, second echo range gate
Last echo subphylum in echo samples value be the sub- echo information of band 4,5,6 coupling result, this little band position
In in second band of imaging region, for the information for the sub- band 4,5,6 being coupled, pass through the association process in later stage
The resolution of single resolution cell can be realized.
S103 specific implementation procedure is as follows：
By echo arrival point (l-1) T of l (l=2,3 ..., L) the individual Stochastic narrow pulse signalP+(τk,j)min
(l-1) T is moved to left on time shaftP, so as to which each burst pulse echo is alignd, wherein,k
=1,2 ..., K, j=1,2 ..., M；
Each burst pulse echo after alignment is divided into multiple echo range gates, each echo range gate is divided
For multiple echo subphylums, each band distance to discrete resolution cell number to be differentiated be Q, i.e. each band
It is ρ containing Q widthrSub- band, the corresponding echo subphylum width of each sub- band is the ρ of δ=2rSin θ/c, Mei Gesuo
Echo range gate width is stated for τ,；
Is sampled in the specific echo subphylum in each echo range gate after alignment, obtain with described in each
The corresponding multigroup echo samples value of band.
It should be noted that present embodiment discloses two kinds of sampling policies, each echo after alignment away from
From being sampled in the specific echo subphylum in door, the multigroup echo samples value corresponding with band each described is obtained, is met
Any one in following two sampling policies：
Sampling policy one：When being not present overlapping between adjacent ribbons, each burst pulse echo after alignment it is every
Sampled in last echo subphylum in individual echo range gate, obtain the multigroup echo corresponding with band each described
Sampling instant t in sampled value, k-th of range gate of burst pulse echo described in l-thl,kMeet τ * k- δ ＜ tl,k＜ τ * k；
Sampling policy two：When adjacent ribbons distance is to overlapping q (q=1,2 ..., Q-1) individual sub- band, described in l-th
K-th of sampling instant t of burst pulse echol,kMeet tl,k=tl,1+ Δ t × (k-1), wherein, tl,1The burst pulse described in l-th
Sampling instant in first range gate of echo, τ-δ ＜ tl,1<τ, Δ t are sampling interval, Δ t=τ-δ × q.
S104：Using the space-time random radiation of each band and corresponding multigroup echo samples value, construct and ask
The relevance imaging equation of each band of solution, obtains the imaging results of each band；
Referring to Fig. 3, S104 implementation procedure is as follows：
S201：When being carried out to the reception echo of each band, empty discretization, the echo for obtaining each band adopts
Sample matrices
Wherein, L is transmitting burst pulse number, tk,lFor k-th of band in l-th of pulse echo corresponding echo samples
Moment,Represent the echo samples value of k-th of band in l-th of pulse echo.
S202：When being carried out to the space-time random radiation of each band, after empty discretization, obtain each band
Observing matrix
Wherein, the l rows of the observing matrix represent the M in band described in k-th that l-th of Stochastic narrow pulse is acted on
M-th of mesh point in amendment radiation field vector at individual mesh point, m row k-th of band of expression is in l-th sampling instant
Correct radiation field vector.
S203：According to the observing matrix of each bandAnd corresponding echo samples value matrixBuild each
The relevance imaging equation of the band：
The relevance imaging equationSpecially：
S204：Relevance imaging equation based on default relevance imaging algorithm and each band, calculates each described
The inverting value of target scattering coefficient in bandAnd by the inverting value of the target scattering coefficientIt is used as the imaging of respective strap
As a result,Wherein,For the operator of the default relevance imaging algorithm.
Calculated it should be noted that the default relevance imaging algorithm can be the association of single order field strength or the association of high-order field strength
Method, meets the optimization algorithm under the conditions of target sparse and has base to follow the trail of (BP) algorithm, orthogonal matching pursuit (OMP) algorithm, sparse
Bayesian learning (SBL) etc., therefore not to repeat here for specific algorithm content.
S105：Imaging results to band each described are spliced, and obtain the overall imaging results of the imaging region.
It should be noted that the disclosed quick relevance imaging method split based on band with data fusion of the present embodiment,
When overlapping region is not present between adjacent ribbons, sampled according to sampling policy one：Each burst pulse after alignment
Sampled in last echo subphylum in each echo range gate of echo, obtain corresponding with band each described
Sampling instant t in multigroup echo samples value, k-th of range gate of burst pulse echo described in l-thl,kMeet τ * k- δ ＜ tl,k
＜ τ * k.In strips mosaic, the imaging results directly to band each described are spliced, you can obtain the imaging region
Overall imaging results, splicing schematic diagram such as Fig. 4 (a) is shown.
In order that band is seamlessly transitted in splicing, it can also allow between the adjacent ribbons of division and exist overlapping, work as adjacent strip
When band distance is to overlapping q (q=1,2 ..., Q-1) individual sub- band, sampled according to above-mentioned sampling policy two：Described in l-th
K-th of sampling instant t of burst pulse echol,kMeet tl,k=tl,1+ Δ t × (k-1), wherein, tl,1The burst pulse described in l-th
Sampling instant in first range gate of echo, τ-δ ＜ tl,1<τ, Δ t be the sampling interval, Δ t=τ-δ × q, now, each
May more than one sampled point in echo range gate.
In strips mosaic, it can use but be not limited to average weighted image interfusion method and handled.Using weighting
When average image interfusion method is handled, the imaging results to band each described are spliced, obtain it is described into
As the overall imaging results in region, including：
Obtain the overlapping region between every group of adjacent ribbons in the imaging region；
It is without any processing to the imaging results outside overlapping region in the imaging region directly to retain, to every group of adjacent strip
Imaging results in two adjacent ribbons in band in overlapping region are weighted summation process, obtain in corresponding overlapping region
Inversion result, and using the inversion result as the imaging results in the overlapping region, be specially：
Wherein, S1And S2Two respectively adjacent bands,WithRespectively two neighboring strips S1And S2In place
PutThe object inversion result at place,For in overlapping regionThe object inversion result at place, w1And w2Respectively two neighboring bar
Band S1And S2The corresponding weight coefficient in overlapping region in inversion chart picture, and meet w1+w2=1, w1,w2∈(0,1)；
Imaging results outside overlapping region in the imaging region and the imaging results each described in overlapping region are entered
Row splicing, obtains the overall imaging results of the imaging region.Splice shown in schematic diagram such as Fig. 4 (b).
It should be noted that mean value method is the special case that above-mentioned weighted sum is handled, i.e. w1=w2=0.5, average value
Method keeps the inversion result outside overlapping region constant, the inversion result averaging in overlapping region is new as overlapping region
Value,
Relative to above-mentioned weighted sum processing, the advantage of mean value method is to calculate simple, and arithmetic speed is fast.
The quick relevance imaging method based on band segmentation and data fusion that the present embodiment is provided, each day of radiating curtain
Orthogonal, the independent Stochastic narrow pulse signal of line locking transmitting space-time, by the imaging region of the tiltedly lower apparent direction of aerial radiation array
The band that multiple width are range gate is divided into upwards in distance；The space-time that aerial radiation array is formed in each band
Random radiation interacts to form the pulse echo of each band with the observed object in respective strap；By each pulse echo
Multiple range gates and multiple echo subphylums are divided, samples, obtains and each institute in the specific echo subphylum in each range gate
State the corresponding multigroup echo samples value of band；Utilize the space-time random radiation of each band and corresponding multigroup echo samples
Value, constructs and solves the relevance imaging equation of each band, obtain the imaging results of each banded zone；In strips mosaic, selection
Image interfusion method seamlessly transits adjacent ribbons, obtains the overall inversion result of imaging region.Sum should be split based on band
According to the relevance imaging method of fusion, because the relevance imaging equation scale of single band per treatment is smaller, therefore can be quick
Imaging.
With reference to specific image scene, to being associated to disclosed in embodiment based on band segmentation and the quick of data fusion
The principle and effect of image space method are described in detail, still, and the practical application and protection domain of the present invention should not be limited with this.
Image scene is as shown in figure 5, systematic parameter is as shown in table 1, and object module is as shown in Figure 6.
The simulated conditions parameter setting of table 1
| Systematic parameter | Parameter setting |
| Aerial radiation array caliber size D | 1.5m×1.5m |
| Antenna radiation unit number N | 25 |
| Imaging plane S and the vertical range of aerial radiation array plane | H=350m |
| Oblique downwards angle of visibility of the aerial radiation array plane to imaging region | θ=45 ° |
| Two-dimensional imaging region area W × W | 120m×120m |
| Imaging region mesh generation number | 40 × 40=1600 |
| Grid spacing ρa=ρr | 3m |
| Transmission signal form | Random frequency hopping |
| The transmission signal pulse spacing | Tp=2us |
| Signal frequency hopping bandwidth | B=500MHz |
| Transmission signal carrier frequency f0 | 10GHz |
In order to illustrate the validity of the quick relevance imaging method based on band segmentation and data fusion disclosed in embodiment,
Recover image quality evaluation standard to recover the signal to noise ratio of image, computational methods are as follows：
Wherein, P, Q are respectively that imaging region orientation and distance are counted to discrete grid block, σ (ri,j),It is respectively former
Beginning target information and inverting target information are in point ri,jThe value at place.
Emulation content
Under above-mentioned simulated conditions, emulated as follows：
Emulation 1, relevance imaging is emulated during bar interband no overlap：
As shown in table 2, the band echo samples moment is located at last echo in each echo range gate to simulation parameter
In subphylum, the bar interband no overlap that now imaging region is divided, only one of which sampled point, overall inverting image in each range gate
Obtained by each band inversion result direct splicing.
Simulation parameter is set during 2 interband no overlaps of table
| Pulse width τ/ns | Divide band number K | Resolution cell number M in band | Umber of pulse L |
| 560 | 1 | 1600 | 2000 |
| 280 | 2 | 800 | 1000 |
| 140 | 4 | 400 | 500 |
| 70 | 8 | 200 | 250 |
| 56 | 10 | 160 | 200 |
By emulation, Fig. 7 gives curve map of the imaging time with bin number change is divided, and Fig. 8 and Fig. 9 are not respectively
The target for dividing band and dividing when band number is 8 recovers image.From simulation result as can be seen that with the increasing for dividing band number
Plus, imaging time is quickly reduced, and compared to band is not divided, is divided after band, image quality is also improved, the result verification
The validity for the quick relevance imaging method split based on band.
Emulation 2, the overlapping data fusion relevance imaging emulation of band：
As shown in table 3, when dividing band, adjacent ribbons distance dimension overlaps two sub- bands to simulation parameter, overlapping part
Discrete resolution cell number is 40 × 2=80, now, and echo samples interval takes the ρ of Δ t=τ-δ × 2=τ -28ns, wherein δ=2rSin θ/c ≈ 14ns, is spliced using mean value method to band.
Simulation parameter is set when the band of table 3 is overlapping
By emulation, the object module that Figure 10 gives when having overlapping between pulse width τ=70ns and adjacent ribbons recovers
Image.From simulation result as can be seen that when bar interband no overlap compared with relevance imaging simulation result, the overlapping data of band are melted
To close relevance imaging can improve image quality, what the result verification was carried herein split based on band and data fusion it is quick
The validity of relevance imaging method.
Figure 11 is referred to, image space is associated to the quick of data fusion based on splitting disclosed in above-described embodiment based on band
Method, the present embodiment correspondence discloses a kind of quick relevance imaging device split based on band with data fusion, including：
Control unit 101, for controlling each antenna radiation unit synchronized transmissions space-time of aerial radiation array orthogonal, only
Vertical Stochastic narrow pulse signal；
Cutting unit 102, for the imaging region of the tiltedly lower apparent direction of the aerial radiation array to be divided upwards in distance
It is cut into multiple bands, the space-time random radiation and respective strap that the aerial radiation array is formed in each band
Observed object interact to form the pulse echo of each band；
Sampling unit 103, for each pulse echo to be divided into multiple echo range gates and multiple echo subphylums,
Sampled in specific echo subphylum in each echo range gate, obtain corresponding with band each described multigroup
Echo samples value；
Structural unit 104, for utilizing the space-time random radiation of each band and corresponding multigroup echo samples
Value, constructs and solves the relevance imaging equation of each band, obtain the imaging results of each band；
Concatenation unit 105, splices for the imaging results to band each described, obtains the whole of the imaging region
Body imaging results.
The quick relevance imaging device based on band segmentation and data fusion that the present embodiment is provided, each day of radiating curtain
Orthogonal, the independent Stochastic narrow pulse signal of line locking transmitting space-time, by the imaging region of the tiltedly lower apparent direction of aerial radiation array
The band that multiple width are range gate is divided into upwards in distance；The space-time that aerial radiation array is formed in each band
Random radiation interacts to form the pulse echo of each band with the observed object in respective strap；By each pulse echo
Multiple echo range gates are divided, while each echo range gate is divided into multiple echo subphylums, in each echo range gate
Specific echo subphylum in sample, obtain the multigroup echo samples value corresponding with band each described；Utilize each band
Space-time random radiation and corresponding multigroup echo samples value, construct and solve the relevance imaging equation of each band, obtain each bar
The imaging results of region；In strips mosaic, selection image interfusion method seamlessly transits adjacent ribbons, obtains imaging region
Overall inversion result.The disclosed relevance imaging device split based on band with data fusion of the present embodiment, because per treatment
Single band relevance imaging equation scale it is smaller, therefore fast imaging can be realized.
The foregoing description of the disclosed embodiments, enables professional and technical personnel in the field to realize or using the present invention.
A variety of modifications to these embodiments will be apparent for those skilled in the art, as defined herein
General Principle can be realized in other embodiments without departing from the spirit or scope of the present invention.Therefore, it is of the invention
The embodiments shown herein is not intended to be limited to, and is to fit to and principles disclosed herein and features of novelty phase one
The most wide scope caused.
Claims (10)
1. a kind of quick relevance imaging method split based on band with data fusion, it is characterised in that including：
The Stochastic narrow pulse signal for controlling each antenna radiation unit synchronized transmissions space-time of aerial radiation array orthogonal, independent；
By the imaging region of the tiltedly lower apparent direction of the aerial radiation array in distance to being divided into multiple bands, the antenna spoke
The space-time random radiation that array formed in each band is penetrated to interact to be formed with the observed object in respective strap
The pulse echo of each band；
Each pulse echo is divided into multiple echo range gates, it is sub- that each echo range gate is divided into multiple echoes
Sampled in door, the specific echo subphylum in each echo range gate, obtain corresponding with band each described
Multigroup echo samples value；
Using the space-time random radiation of each band and corresponding multigroup echo samples value, construct and solve each described
The relevance imaging equation of band, obtains the imaging results of each band；
Imaging results to band each described are spliced, and obtain the overall imaging results of the imaging region.
2. according to the method described in claim 1, it is characterised in that each antenna radiation unit is equal in the aerial radiation array
Even to be distributed but be not limited to be uniformly distributed, the position of the antenna radiation unit isN=1 ..., N, N are the aerial radiation list
The summation of first quantity, the position of single-geophone receiver machine isEach described antenna radiation unit synchronized transmissions random frequency hopping but it is not limited to
The narrow pulse signal of random frequency hopping, the transmission signal of antenna radiation unit described in n-th is expressed as：
<mrow>
<msub>
<mi>f</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>L</mi>
</munderover>
<mi>r</mi>
<mi>e</mi>
<mi>c</mi>
<mi>t</mi>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
</mrow>
<mi>&tau;</mi>
</mfrac>
<mo>&rsqb;</mo>
<mi>exp</mi>
<mo>{</mo>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mrow>
<mi>n</mi>
<mi>l</mi>
</mrow>
</msub>
<mo>&lsqb;</mo>
<mi>t</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>T</mi>
<mi>p</mi>
</msub>
<mo>&rsqb;</mo>
<mo>}</mo>
<mo>,</mo>
</mrow>
Wherein, fnlThe tranmitting frequency of l-th of pulse of antenna radiation unit described in n-th, l=1 ..., L, L is overall pulse
Number, τ is burst pulse duration, TpFor pulse spacing, Tp＞ Tmin+ τ+Δ τ, wherein, when Δ τ is the broadening of burst pulse echo
Between, Δ τ=2 (Rmax-Rmin)/c=Tmax-Tmin≈ 2Wsin θ/c, Rmax、RminRespectively described aerial radiation array to it is described into
As the distance of region proximally and distally, c=3 × 108M/s, W be imaging region distance to width, θ is the aerial radiation array
Relative to the oblique downwards angle of visibility of the imaging region.
3. according to the method described in claim 1, it is characterised in that the tiltedly lower apparent direction by the aerial radiation array
Imaging region is divided into multiple bands upwards in distance, including：
According to the range gate feature of short-pulse radar, by the imaging region of the tiltedly lower apparent direction of the aerial radiation array in distance
Multiple different bands are divided into upwards, and the strip width of duration τ pulse covering is w=c τ/2sin θ, and band number is K
=[W/w]=[2Wsin θ/c τ], [x] is round up function, c=3 × 108m/s；
Spatial discretization processing is carried out to the imaging region, each band M grid, M=P × Q is divided into, P is institute
State imaging region orientation resolution cell number, Q is the band distance to resolution cell number, ρa、ρrRespectively described imaging region
Orientation and distance are to the grid spacing of discretization, and the position of j-th of grid dot center described in k-th in band is designated asInstitute
The target backscattering coefficient that rheme is put isK-th of band of observation area is expressed as
4. according to the method described in claim 1, it is characterised in that the aerial radiation array is formed in each band
Space-time random radiation interact to form the pulse echo of each band with the observed object in respective strap, specially：
Each the described bar of the time-varying random field for the random signal formation that each antenna radiation unit is produced in wave cover
Band is interior to be overlapped, during synthesis, empty bidimensional random radiation field, any point in kth (k=1,2 ..., K) individual described bandThe in-field at place is：
<mrow>
<msup>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mfrac>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>/</mo>
<mi>c</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
1
The in-field and the observed object scattered field to be formed that interacts are multicast to by the second pass in free space
Up to receiver, it is for locusThe pulse echo of k-th of band that receives of receiver be：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mfrac>
<mrow>
<msup>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mi>c</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>|</mo>
<mo>/</mo>
<mi>c</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mfrac>
<mrow>
<mi>f</mi>
<mo>&lsqb;</mo>
<mi>t</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>c</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
It is defined into the amendment radiation field up to the receiver, any point in kth (k=1,2 ..., K) individual described band
The amendment radiation field at place can be expressed as：
<mrow>
<msup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>d</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mfrac>
<mrow>
<mi>f</mi>
<mo>&lsqb;</mo>
<mi>t</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
</mrow>
<mo>|</mo>
<mo>+</mo>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>|</mo>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>c</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>n</mi>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mo>|</mo>
</mrow>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
The scatter echo of k-th of band has following relation with amendment radiation field：
<mrow>
<msup>
<mi>E</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>d</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
5. according to the method described in claim 1, it is characterised in that described that each pulse echo is divided into multiple echoes
Range gate, multiple echo subphylums are divided into by each echo range gate, and the specific echo in each echo range gate is sub-
Sampled in door, obtain the multigroup echo samples value corresponding with band each described, including：
By echo arrival point (l-1) T of l (l=2,3 ..., L) the individual Stochastic narrow pulse signalP+(τk,j)minIn time shaft
On move to left (l-1) TP, so as to which each burst pulse echo is alignd, wherein,K=1,
2 ..., K, j=1,2 ..., M；
Each burst pulse echo after alignment is divided into multiple echo range gates, each echo range gate is divided into many
Individual echo subphylum, each band distance to discrete resolution cell number to be differentiated be Q, i.e. each band contains Q
Individual width is ρrSub- band, the corresponding echo subphylum width of each sub- band is the ρ of δ=2rSin θ/c, each echo
Range gate width is τ；
Sampled, obtained and band each described in the specific echo subphylum in each echo range gate after alignment
Corresponding multigroup echo samples value.
6. method according to claim 5, it is characterised in that in each echo range gate after alignment
Sampled in specific echo subphylum, obtain the multigroup echo samples value corresponding with band each described, meet following two
Any one in sampling policy：
Sampling policy one：When being not present overlapping between adjacent ribbons, each of each burst pulse echo after alignment returns
Sampled in last echo subphylum in ripple range gate, obtain the multigroup echo samples corresponding with band each described
Sampling instant t in value, k-th of range gate of burst pulse echo described in l-thl,kMeet τ * k- δ ＜ tl,k＜ τ * k；
Sampling policy two：When adjacent ribbons distance is to overlapping q (q=1,2 ..., Q-1) individual sub- band, narrow arteries and veins described in l-th
K-th of sampling instant t of refunds ripplel,kMeet tl,k=tl,1+ Δ t × (k-1), wherein, tl,1The burst pulse echo described in l-th
First range gate in sampling instant, τ-δ ＜ tl,1<τ, Δ t are sampling interval, Δ t=τ-δ × q.
7. according to the method described in claim 1, it is characterised in that the space-time random radiation using each band
And corresponding multigroup echo samples value, the relevance imaging equation of each band is constructed and solved, each band is obtained
Imaging results, including：
When being carried out to the reception echo of each band, after empty discretization, the echo samples value square of each band is obtained
Battle array
<mrow>
<msubsup>
<mi>E</mi>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>E</mi>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>E</mi>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>E</mi>
<mi>k</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>a</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>L</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>;</mo>
</mrow>
When being carried out to the space-time random radiation of each band, after empty discretization, obtain the observation square of each band
Battle array
According to the observing matrix of each bandAnd corresponding echo samples value matrixBuild each band
Relevance imaging equation：I.e.
Relevance imaging equation based on default relevance imaging algorithm and each band, calculates target in each band
The inverting value of scattering coefficientAnd by the inverting value of the target scattering coefficientAs the imaging results of respective strap,Wherein,For the operator of the default relevance imaging algorithm.
8. method according to claim 6, it is characterised in that when methods described meets the sampling policy one, adjacent strip
When interband is not present overlapping, the imaging results to band each described are spliced, and obtain the entirety of the imaging region
Imaging results, be specially：
Imaging results directly to band each described are spliced, and obtain the overall imaging results of the imaging region.
9. method according to claim 6, it is characterised in that when methods described meets the sampling policy two, adjacent strip
When band distance is to overlapping q (q=1,2 ..., Q-1) individual sub- band, the imaging results to band each described are spliced,
The overall imaging results of the imaging region are obtained, including：
Obtain the overlapping region between every group of adjacent ribbons in the imaging region；
It is without any processing to the imaging results outside overlapping region in the imaging region directly to retain, in every group of adjacent ribbons
Two adjacent ribbons in imaging results in overlapping region be weighted summation process, obtain the inverting in corresponding overlapping region
As a result, and using the inversion result as the imaging results in the overlapping region, it is specially：
<mrow>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>w</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>w</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>&Element;</mo>
<mo>(</mo>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
<mo>&cap;</mo>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>&Element;</mo>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Wherein, S1And S2Two respectively adjacent bands,WithRespectively two neighboring strips S1And S2In position
The object inversion result at place,For in overlapping regionThe object inversion result at place, w1And w2Respectively two neighboring strips S1
And S2The corresponding weight coefficient in overlapping region in inversion chart picture, and meet w1+w2=1, w1,w2∈(0,1)；
Imaging results outside overlapping region in the imaging region and the imaging results each described in overlapping region are spelled
Connect, obtain the overall imaging results of the imaging region.
10. a kind of quick relevance imaging device split based on band with data fusion, it is characterised in that including：
Control unit, for control each antenna radiation unit synchronized transmissions space-time of aerial radiation array it is orthogonal, it is independent with
Machine narrow pulse signal；
Cutting unit, for by the imaging region of the tiltedly lower apparent direction of the aerial radiation array in distance to being divided into multiple
Band, space-time random radiation and the observed object in respective strap that the aerial radiation array is formed in each band
Interaction forms the pulse echo of each band；
Sampling unit, for each pulse echo to be divided into multiple echo range gates, each echo range gate is divided
For multiple echo subphylums, sampled in the specific echo subphylum in each echo range gate, obtain with described in each
The corresponding multigroup echo samples value of band；
Structural unit, for utilizing the space-time random radiation of each band and corresponding multigroup echo samples value, construction
And the relevance imaging equation of each band is solved, obtain the imaging results of each band；
Concatenation unit, splices for the imaging results to band each described, obtains the overall imaging of the imaging region
As a result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710385506.2A CN107102326B (en) | 2017-05-26 | 2017-05-26 | Rapid correlation imaging method and device based on strip segmentation and data fusion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710385506.2A CN107102326B (en) | 2017-05-26 | 2017-05-26 | Rapid correlation imaging method and device based on strip segmentation and data fusion |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107102326A true CN107102326A (en) | 2017-08-29 |
| CN107102326B CN107102326B (en) | 2020-03-31 |
Family
ID=59668998
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710385506.2A Active CN107102326B (en) | 2017-05-26 | 2017-05-26 | Rapid correlation imaging method and device based on strip segmentation and data fusion |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107102326B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107976673A (en) * | 2017-11-17 | 2018-05-01 | 中国科学技术大学 | Improve the MIMO radar imaging method of large scene target imaging quality |
| CN108828593B (en) * | 2018-04-26 | 2020-06-02 | 电子科技大学 | Random radiation radar correlation imaging method |
| CN111273228A (en) * | 2019-05-09 | 2020-06-12 | 哈尔滨工程大学 | Multi-coherent source positioning method based on traversal search strategy |
| CN112630796A (en) * | 2021-01-07 | 2021-04-09 | 苏州蛟视智能科技有限公司 | Large-view-field associated imaging device and imaging method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103837873A (en) * | 2014-03-14 | 2014-06-04 | 中国科学技术大学 | Microwave and stare correlated imaging system and method based on floating platform and intensive array antennae |
| US20160195613A1 (en) * | 2015-01-05 | 2016-07-07 | Robert M. Knox | Dual Mode Undercarriage Vehicle Inspection System |
| CN106597434A (en) * | 2016-11-28 | 2017-04-26 | 中国人民解放军国防科学技术大学 | Agile satellite target decomposition method and system based on push-scan trajectories |
-
2017
- 2017-05-26 CN CN201710385506.2A patent/CN107102326B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103837873A (en) * | 2014-03-14 | 2014-06-04 | 中国科学技术大学 | Microwave and stare correlated imaging system and method based on floating platform and intensive array antennae |
| US20160195613A1 (en) * | 2015-01-05 | 2016-07-07 | Robert M. Knox | Dual Mode Undercarriage Vehicle Inspection System |
| CN106597434A (en) * | 2016-11-28 | 2017-04-26 | 中国人民解放军国防科学技术大学 | Agile satellite target decomposition method and system based on push-scan trajectories |
Non-Patent Citations (2)
| Title |
|---|
| YUANYUE GUO等: "Research on Sensing Matrix Characteristics in Microwave Staring Correlated Imaging Based on Compressed Sensing", 《2014 IEEE INTERNATIONAL CONFERENCE ON IMAGING SYSTEMS AND TECHNIQUES (IST) PROCEEDINGS》 * |
| 何学智: "微波凝视关联成像的信息处理方法与仿真", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107976673A (en) * | 2017-11-17 | 2018-05-01 | 中国科学技术大学 | Improve the MIMO radar imaging method of large scene target imaging quality |
| CN107976673B (en) * | 2017-11-17 | 2020-02-07 | 中国科学技术大学 | MIMO radar imaging method for improving imaging quality of large-scene target |
| CN108828593B (en) * | 2018-04-26 | 2020-06-02 | 电子科技大学 | Random radiation radar correlation imaging method |
| CN111273228A (en) * | 2019-05-09 | 2020-06-12 | 哈尔滨工程大学 | Multi-coherent source positioning method based on traversal search strategy |
| CN112630796A (en) * | 2021-01-07 | 2021-04-09 | 苏州蛟视智能科技有限公司 | Large-view-field associated imaging device and imaging method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107102326B (en) | 2020-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107102326A (en) | Quick relevance imaging method and device with data fusion is split based on band | |
| CN107976673A (en) | Improve the MIMO radar imaging method of large scene target imaging quality | |
| CN103633452B (en) | A kind of antenna and wireless signal sending, receiving method | |
| CN104597444B (en) | Microwave gaze high resolution imaging method based on intensity association | |
| DE112018004451T5 (en) | CUSTOMIZE TRANSMISSIONS BASED ON DIRECTLY SCANING THE IONOSPHERE | |
| CN107678028A (en) | Microwave under Low SNR stares relevance imaging method | |
| CN108008380A (en) | A kind of microwave based on satellites formation stares relevance imaging method and system | |
| CN105841705A (en) | Time-attitude based method for decomposition and synthesis of imaging satellite observation task | |
| CN103323845B (en) | Image inversion method of non-uniform sampling comprehensive bore diameter radiometer | |
| CN105589070A (en) | Frequency diversity array-based radar target imaging method | |
| CN105137430A (en) | Forward-looking array SAR echo sparse acquisition and three-dimensional imaging method | |
| Malanowski et al. | Real-time high-resolution SAR processor using CUDA technology | |
| CN102879768A (en) | Satellite-borne synthetic aperture radar (SAR) high-fidelity echo simulation method based on steady-state radar cross section (RCS) | |
| CN105044714A (en) | Satellite formation-based imaging method | |
| CN104199020A (en) | Multi-frame information fusion based meter wave array radar target elevation measuring method | |
| CN107102328B (en) | FPGA-based real-time imaging signal processing method and FPGA | |
| CN103439888B (en) | Multi-channel space-borne SAR antenna control error compensation method | |
| Yue et al. | Preliminary experimental results by the prototype of Sanya Incoherent Scatter Radar | |
| Janches et al. | An initial meteoroid stream survey in the southern hemisphere using the Southern Argentina Agile Meteor Radar (SAAMER) | |
| CN102998655B (en) | Optimal wave position selecting method of synthetic aperture radar under scanning working mode | |
| CN108335263A (en) | A kind of microwave remote sensing image super-resolution rebuilding method based on VDSR | |
| Kravchenko et al. | The jet of S5 0716+ 71 at μas scales with RadioAstron | |
| CN105320787A (en) | Method for constructing dynamic evaluation model of network link reliability | |
| Alam et al. | Refractivity variations and propagation at ultra high frequency | |
| CN106908793A (en) | Ka wave band along-trade INSAR systems and its method of work |
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 |