The content of the invention
Therefore, the present invention provides a kind of radar calibration scheme, effectively solve it is above-mentioned at least one problem.
According to an aspect of the present invention there is provided a kind of radar calibration device, including:Reception antenna, suitable for receiving radar
Pulse signal Srx;Demodulating unit, suitable for extracting intermediate-freuqncy signal Sdc from Srx;Signal processing unit, suitable for based on intermediate-freuqncy signal
Sdc generates frequency shift signal and delay time signal;Transmitter unit, suitable for behaviour is mixed and amplified to frequency shift signal and delay time signal
Make, to generate transmitting time delay to be launched and frequency shift signal;Transmitting antenna, suitable for transmitting time delay and frequency shift signal, so as to radar
Receive the signal.
Alternatively, in the calibrator according to the present invention, demodulating unit is suitable to filter carrier wave from Srx according to following manner
Signal, radar and the intermediate-freuqncy signal Sdc of calibrator distance are included to extract:
Alternatively, in the calibrator according to the present invention, the signal processing unit includes:Frequency shift unit, suitable for generation
The frequency shift signal of the bright spot of orientation in the radar image of correspondence radar generation;Time delay elements, suitable for generation correspondence radar
The delay time signal of the bright spot of the radar image middle-range descriscent of generation.
According to a further aspect of the invention there is provided a kind of radar installations, including:Reception antenna, thunder is come from suitable for receiving
Up to the time delay and frequency shift signal of calibrator;Wave filter, suitable for eliminating time delay and frequency based on Wigner-Ville distribution and Hankel
Time frequency signal interferes with each other in shifting signal;Geometry correction unit, suitable for generating radar image according to time delay and frequency shift signal, and
According to the coordinate pair of radar image lieutenant colonel on schedule, the radar image carries out geometric correction;Radiant correction unit, suitable for according to time delay
With the energy balane radiant correction coefficient of frequency shift signal, and in radar image pixel carry out radiant correction.
Alternatively, in the radar installations according to the present invention, the time delay and frequency shift signal received is suitable to pass through following sides
Formula is represented:
Wherein, σ:The radar cross section (RCS) of object
λ:Radar wavelength
η:Traveling time of the radar along orientation
T is rapid time, that is, signal propagation time.Oblique distance R is one when being moved on radar installations under speed u
The slow time function.
Alternatively, in the radar installations according to the present invention, geometry correction unit is suitable to according to following manner to radar map
As carrying out geometric correction:The default spacing conduct in radar image bright spot is produced according to the frequency shift unit and the time delay elements
Precision reference control point carries out geometric correction, wherein default spacing is by setting the frequency shift amount and delay volume of radar calibration device come really
It is fixed.
Alternatively, according to the present invention radar installations in, radiant correction unit be suitable to according to following manner come according to when
Prolong the energy balane radiant correction COEFFICIENT K with frequency shift signal:
σI, j=KDNI, j
Wherein σI, j:Pixel ij backscattering system, Dij:Pixel ij GTG angle value,
K:The backscattering conversion coefficient (radiant correction coefficient) of Synthetic Aperture Radar images;
Radiant correction COEFFICIENT K is calculated according to following formula,
Wherein, PrTThe reflection power of bright spot target, ρ are produced for calibratorrg、ρazRespectively the oblique distance and orientation of image are divided
Resolution, σTFor the radar cross section of object, α is radar incidence angle.
Alternatively, in the radar installations according to the present invention, wave filter is suitable to according to following manner come based on Wigner-
Ville is distributed and Hankel eliminates interfering for time frequency signal:
Using displacement of the corresponding Doppler frequency of maximum relative to transfer point of the WVD correlation functions exported,
Produce the high-resolution simultaneously in time delay and frequency displacement;
Multiple signals are disassembled using Hankel conversions, each single signal is other separated, then input to WVD
In, to eliminate interfering for time frequency signal, and the linear FM characteristic in stick signal.
According to a further aspect of the invention there is provided a kind of radar system, including according to the present invention radar calibration device and
Radar installations.
According to a further aspect of the invention there is provided a kind of radar calibration method, suitable for being performed in radar calibration device.Should
Method includes:Receive the pulse signal Srx of radar;Intermediate-freuqncy signal Sdc is extracted from Srx;Based on intermediate-freuqncy signal Sdc, generation frequency
Shifting signal and delay time signal;Frequency shift signal and delay time signal are mixed and amplifieroperation, with generate time delay to be launched and
Frequency shift signal;Launch time delay and frequency shift signal, so that radar receives the time delay and frequency shift signal.
According to a further aspect of the invention there is provided a kind of method that radar installations is performed, including:Receive and come from radar school
The time delay and frequency shift signal of quasi- device;Time frequency signal in time delay and frequency shift signal is eliminated based on Wigner-Ville distribution and Hankel
Interfere with each other;According to time delay and frequency shift signal generation radar image, and according to the coordinate pair of the radar image lieutenant colonel on schedule thunder
Geometric correction is carried out up to image;Radiant correction unit, suitable for the energy balane radiant correction coefficient according to time delay and frequency shift signal,
And radiant correction is carried out to pixel in radar image.
, can be along range direction to received radar by radar calibration device according to the radar calibration technology of the present invention
Pulse signal carries out time delay operations, and carries out frequency displacement operation to radar pulse signal along azimuth direction.So, radar installations root
The original image generated according to time delay and frequency shift signal can include two-dimension netted calibration point.Each calibration point is to be based on frequency displacement
Operation and time delay operations and the result that occurs.On this basis, radar installations can improve radiation according to netted calibration point
Correction and the precision of geometric correction.
Embodiment
The exemplary embodiment of the disclosure is more fully described below with reference to accompanying drawings.Although showing the disclosure in accompanying drawing
Exemplary embodiment, it being understood, however, that may be realized in various forms the disclosure without should be by embodiments set forth here
Limited.On the contrary, these embodiments are provided to facilitate a more thoroughly understanding of the present invention, and can be by the scope of the present disclosure
Complete conveys to those skilled in the art.
Fig. 1 shows the application scenarios schematic diagram of radar system 100 according to some embodiments of the invention.
As shown in figure 1, radar installations 110 (rader) is deployed in such as satellite, aircraft on aircraft (not shown).Allusion quotation
Type, radar installations 110 is, for example, synthetic aperture radar (SAR).Based on the three-dimensional system of coordinate (x, y, z) shown in Fig. 1, radar
Device is h relative to the flying height on ground.Radar installations is moved with speed u along x-axis forward direction, and with incident angle α and along distance
Launch pulse signal in direction.Here pulse signal is, for example, chirp pulse signal.Chirp (Linear
Frequency Modulation, LFM) it is the signal that a frequency changes over time, it is referred to as chirping of birds communication signal
(chirp signal).The frequency of signal can rise (up chirp) over time during modulation or decline (down
Chirp change) and linearly.In one embodiment, the mode of signal frequency change is as follows:
ft=f0+aTt
f0Represent initial frequency, aTRepresent the rate of change (chirp rate) of frequency.
According to one embodiment of the invention, the mathematic(al) representation of LFM pulse signals is:
S (t)=exp (j π art2)
The scanning area (swath) of radar installations 110 includes location point P.P position point is suitable to thunder of the arrangement according to the present invention
Up to calibrator 120 (TF-IRC).Here, the position coordinates (longitude and latitude) of radar calibration device 120 is determined.Radar calibration device 120 can
To carry out time delay operations to received radar pulse signal along range direction, and along azimuth direction to radar pulse signal
Carry out frequency displacement operation.So, the original image that radar installations 110 is generated according to time delay and frequency shift signal can include two-dimensional mesh
The calibration point of shape.Each calibration point is the result occurred based on frequency displacement operation and time delay operations.On this basis, radar installations
110 can improve the precision of radiant correction and geometric correction according to netted calibration point.With reference to Fig. 2 to according to the present invention
Radar calibration device it is illustrative.
Fig. 2 shows the schematic diagram of radar calibration device 200 according to some embodiments of the invention.
As shown in Fig. 2 radar calibration device 200 includes reception antenna 210, demodulating unit 220, signal processing unit 230, hair
Penetrate unit 240 and transmitting antenna 250.Here, reception antenna 210 both can be discrete or multiplexing, and the present invention is right
This does not do excessive limitation.After reception antenna 210 receives radar pulse signal Srx, demodulating unit 220 is suitable to carry from Srx
Take intermediate-freuqncy signal Sdc.Signal processing unit 230 can generate two paths of signals based on Sdc.Wherein, all the way to pass through time delay operations
Signal St, another road be by frequency displacement operate signal Sp.According to one embodiment of the invention, the quilt of signal processing unit 230
It is configured to include frequency shift unit and time delay elements.Frequency shift unit is suitable to orientation in the radar image of generation correspondence radar generation
Bright spot frequency shift signal.Time delay elements are suitable to the time delay of the bright spot of the radar image middle-range descriscent of generation correspondence radar generation
Signal.Here, signal processing unit 230 can be application specific integrated circuit (ASIC), can also be realized based on fpga chip, this
In repeat no more.For example, by the FPGA of programmable, radar system 100 is able to produce the school of different frequency displacement and time shift distances
On schedule.Frequency displacement and time shift distance due to generation can be adjusted voluntarily, are very helpful for the geometric correction on image.Simultaneously
Radar calibration device 200 is difficult the characteristic of decay based on digital data transmission so that the power output of frequency displacement and time shift point is more steady
It is fixed.
Transmitter unit 240 is suitable to St and Sp being mixed and amplifieroperation, is then launched by transmitting antenna 250 and returned
To the signal of radar installations.
Fig. 3 shows the schematic diagram of demodulating unit 220 according to an embodiment of the invention.Demodulating unit 220 is following
Formula carries out signal transacting.
Fig. 4 shows the schematic diagram of radar installations 400 according to some embodiments of the invention.As shown in figure 4, radar installations
400 include reception antenna 410, wave filter 420, geometry correction unit 430 and radiant correction unit 440.Reception antenna 410, is fitted
In receiving time delay from radar calibration device and frequency shift signal.Wave filter 420 be suitable to based on Wigner-Ville distribution and
Time frequency signal interferes with each other in Hankel elimination time delays and frequency shift signal.Geometry correction unit 430 is suitable to according to time delay and frequency
Shifting signal generates radar image, and the radar image carries out geometric correction according to the coordinate pair of radar image lieutenant colonel on schedule.Radiation
Unit 440 is corrected suitable for the energy balane radiant correction coefficient according to time delay and frequency shift signal, and to pixel in radar image
Carry out radiant correction.
Fig. 5 shows radar image schematic diagram according to an embodiment of the invention.As shown in figure 5, figure include it is multiple
Bright spot.Each bright spot represents a calibration point (calibrator for example as shown in Figure 2 carries out frequency displacement and time delay operations and produced), i.e.,
Determine the radar cross section (Radar cross section, RCS) in place.So, although Fig. 5 radar image has not been
Whole earth's surface image (in other words, radar can not quick obtaining to whole radar image coordinate), radar installations can pass through school
(its gray value and coordinate are determined) to carry out radiant correction and geometric correction to whole radar image on schedule.Wherein, geometric correction
Operation be mainly based upon the coordinate value of calibration point, according to the signal delay time of calibration point, longitude and latitude is done to whole radar image
Registration process is spent to realize.
It is illustrative to geometry and radiant correction operation with reference to Fig. 1,4 and 5.
Firstly the need of explanation, the time delay received and frequency shift signal of radar installations are as follows:
Wherein, σ:The radar cross section (RCS) of object
λ:Radar wavelength
η:Traveling time of the radar along orientation
T is rapid time, that is, signal propagation time.Oblique distance R is one and moved down on radar installations 100 in speed u
The function of slow time when dynamic.
The phase of wherein orientation is exactlyDistance to phase be The side of being
Position chirp rate.
In an embodiment in accordance with the invention, geometry correction unit 430 is entered with specific reference to following manner to radar image
Row geometric correction:The default spacing in radar image (Fig. 5) bright spot is produced as precision ginseng according to frequency shift unit and time delay elements
Examine control point and carry out geometric correction.Wherein, spacing is preset to determine by setting the frequency shift amount and delay volume of radar calibration device.
In addition, can be derived from the backward energy mathematical formulae that radar installations 400 is received by above-mentioned radar equation:
Wherein,
Pr:The echo power that SAR system is received
Pt:The power of SAR system transmitting
G:The gain of antenna
λ:The wavelength of radar wave
G:The distance of antenna and object
Ls:The loss of air in transmitting procedure
σ:The radar cross section (RCS) of object
For the region (Homogeneous region) of a piece of homogeneous atural object, RCS can be write as:
α:Incidence angle of the incidence wave to object
ρrg:Oblique distance directional resolution
ρaz:Azimuth direction resolution ratio
σ0:Backscattering coefficient
Thus, it is possible to which backward energy equation is as follows:
The method for calculating the backscattering coefficient of homogeneous area is as follows.Known radar section is put in the region of a homogeneous
Product is σTObject, and the image energy of object be PrT, and homogeneous area is AHHomogeneous region, received energy
For PrH.Then the backscattering coefficient of this homogeneous area can be obtained so
For homogeneous area AH
K:The backscattering conversion coefficient (radiant correction coefficient) of Synthetic Aperture Radar images.Homogeneous region and known radar
Sectional area σTThe distance of object to radar antenna have to close or equal, can just obtain close, correct correction coefficient
K, if oblique distance gap is excessive on image, may set up multiple known RCS values along oblique distance direction to calculate object, finally
Obtain the correction coefficient of multiple different oblique distances.
By above-mentioned formula, we can obtain radar image backscattering conversion coefficient:
PrT:The reflection power of object
ρrg:Oblique distance directional resolution
ρaz:Azimuth direction resolution ratio
σT:The radar cross section of object
α:Radar incidence angle of the incidence wave to object
Calculate K value modes as follows:The image of Digital Synthetic Aperture Radar includes the imaging contexts of homogeneous area and object.
One piece of area is first cut out from image for AHHomogenous area.N is contained in the homogenous areaHIndividual pixel, and ensure this homogeneous
Without other objects in region.Now, the reflection general power expression of the homogenous area is:
αij:Quantized value (digital number) on image
Power density expression is:
Then cut out from image and contain N in another piece of region containing object and background, this regionpIndividual pixel,
And the reflection general power expression for containing object and background is:
It can finally calculate and comprise only the reflection power of object and be:
The general power that can calculate object by above formula is PrT, add ρrg、ρazIt can be learnt in image, and object
Radar cross section σTAnd incident angle α is, it is known that therefore correction coefficient K can be calculated.
In an embodiment in accordance with the invention, radiant correction unit 440 be suitable to according to following manner come according to time delay and
The energy balane radiant correction COEFFICIENT K of frequency shift signal:
σI, j=KDNI, j
Wherein σI, j:Pixel ij backscattering system, Dij:Pixel ij GTG angle value,
K:The backscattering conversion coefficient (radiant correction coefficient) of Synthetic Aperture Radar images;
Radiant correction COEFFICIENT K is calculated according to following formula,
Wherein, PrTThe reflection power of bright spot target, ρ are produced for calibratorrg、ρazRespectively the oblique distance and orientation of image are divided
Resolution, σTFor the radar cross section of object, α is radar incidence angle.
Wigner-Ville (WVD) distributions can all reach high-resolution over time and frequency.But for synthesis hole
For the radar of footpath, time delay and frequency shift signal be in the range of radar beam from return relevant of each individual goal and.This will
Interference is formed between multiple signals.Fortunately, Fourier's Bezier expansion can be decomposed preferably compared to WVD analyses
These multiple signals.
Radar installations 400 also includes being used for eliminating based on the Novel Filter that Wigner-Ville distribution and Hankel are converted
The interference of time frequency signal.
The operation principle of wave filter is described as follows:
The power spectrum P (ω) for describing a time domain communication signal s (t) is the auto-correlation function of this communication signal in itself
Fourier transform, its expression is as follows:
Wherein R (τ) is time domain communication signal s (t) and the auto-correlation function of time independence, is expressed as follows:
However, this power spectrum only shows the composition of communication signal frequency in this time range, and without Faxian
The composition for showing frequency how is distributed in this time range, that is, frequency is how to change over time.Therefore
Redefine its expression with auto-correlation function R (t, the τ) of time correlation as follows:
The R (τ) of (1) formula is replaced with R (t, τ) formula of (3) formula, following formula is can obtain
After the replacement of auto-correlation function, power spectrum has become the function with time correlation, then by arranging, can
Define WVD as follows:
Then the mathematic(al) representation of LFM communication signals is:
S (t)=exp (j π art2) (6)
Then WVD is:
And the time delay and frequency shift signal for adding Doppler frequency can be write as following formula:
Wherein fDopplerFor Doppler frequency
(9) formula of utilization makees correlation output with (10) formula and derives such as following formula:
Wherein TaFor generated time
By (11), formula is known that Scorr(toffset) maximum must be in toffsetShift time (offset time) isWhen have maximum, therefore it is recognised that Doppler frequency is relative to the displacement of transfer point
Formula is shown below:
In an embodiment in accordance with the invention, when wave filter 420 is based on Wigner-Ville distribution and Hankel eliminations
The course of work interfered of frequency signal is as follows.Utilize the corresponding Doppler frequency of maximum of the WVD correlation functions exported
Relative to the displacement of transfer point, the high-resolution simultaneously in time delay and frequency displacement is produced.Disassembled using Hankel conversions
Multiple signals, each single signal is other separated, then input into WVD, to eliminate interfering for time frequency signal,
And the linear FM characteristic in stick signal.
Using displacement of the corresponding Doppler frequency of maximum relative to transfer point of the WVD correlation functions exported,
Very high resolution ratio can be produced in time delay and frequency displacement simultaneously.Because calibrator signal is multiple signal, i.e., multiple time delays
Point and many frequency displacement points, can be interfered in WVD times and frequency.The interference of time frequency signal can be disassembled using Hankel conversions
Multiple signals, it is each single signal is other separated, then input into WVD, you can the problem of to solve distracter, and
Linear FM characteristic in signal can also completely retain.
In addition, the radar cross section (RCS) for providing maximum according to the amplifying power of the radar calibration device of the present invention can be full
Sufficient radar system can receive the RCS of maximum.It is therefore assumed that the dynamic range that radar system can detect maximum is backscattering system
Number ρ0 max, radar system resolution ratio is A, then the maximum radar cross section that radar system can be detected is:
σDR_max=ρ0 max×A (13)
The radar cross section of active radar calibration device is represented by:
When the sectional area of active radar calibration device is equal to the maximum radar cross section that may detect that, we can assess
The gain amplifier G of active radar calibration deviceampFor:
Attenuation range:
The attenuation range of the program-controlled attenuator of design system, is according to the dynamic range that can be detected in radar system.
Peak power output:
When designing active radar calibration device, it is necessary to its peak power output is considered, to avoid intraware
Beyond working range, so first understanding the communication signal intensity that active radar calibration device receives maximum first, relation is as follows:
P in above formulatIt is the transmission power of SAR system, G is the gain of SAR system transmitting antenna, RnearBe it is airborne to ground
The distance of table low coverage point.Have receive communication signal size, we combine 1.25 and 1.26 formulas can learn we maximum
Power output is:
Therefore, the power that the component of radar calibration device of the invention can bear, will be more than peak power output, this
Sample just can ensure that component operational is normal.
Time delay
In SAR width Wide, temporal delay also has limitation.If delay exceedes greatly very much covering scope, that is in shadow
Do not show not come out as upper, so time difference τ of the scope of our time delays no more than low coverage point to remote point yet1, such as following formula
It is shown:
Radar system echo communication signal emulator:
Radar system echo communication signal emulator is with the maximum difference of active radar calibration device before, emulator
When doing echo simulation to radar system indoors, not via antenna but be directly connected to inside loop because radar system
Communication signal out is very high for rate, so communication signal must first pass through high power attenuator, decaying to can be with emulator work
The scope of work.Pad value is as follows:
Attenuation range be also the dynamic range according to radar system as foundation, the part of time delay because being echo
Emulator, so being the echo communication signal between emulating communication signal from low coverage point to remote point, so time delay τ2Such as
Under:
In the specification that this place is provided, numerous specific details are set forth.It is to be appreciated, however, that the implementation of the present invention
Example can be put into practice in the case of these no details.In some instances, known method, knot is not been shown in detail
Structure and technology, so as not to obscure the understanding of this description.
Similarly, it will be appreciated that in order to simplify the disclosure and help to understand one or more of each inventive aspect, exist
Above in the description of the exemplary embodiment of the present invention, each feature of the invention is grouped together into single implementation sometimes
In example, figure or descriptions thereof.However, the method for the disclosure should be construed to reflect following intention:It is i.e. required to protect
The application claims of shield are than the feature more features that is expressly recited in each claim.More precisely, as following
As claims reflect, inventive aspect is all features less than single embodiment disclosed above.Therefore, abide by
Thus the claims for following embodiment are expressly incorporated in the embodiment, wherein each claim is in itself
It is used as the separate embodiments of the present invention.
Those skilled in the art should be understood the module or unit or group of the equipment in example disclosed herein
Part can be arranged in equipment as depicted in this embodiment, or alternatively can be positioned at and the equipment in the example
In different one or more equipment.Module in aforementioned exemplary can be combined as a module or be segmented into addition multiple
Submodule.
Those skilled in the art, which are appreciated that, to be carried out adaptively to the module in the equipment in embodiment
Change and they are arranged in one or more equipment different from the embodiment.Can be the module or list in embodiment
Member or component be combined into a module or unit or component, and can be divided into addition multiple submodule or subelement or
Sub-component.In addition at least some in such feature and/or process or unit exclude each other, it can use any
Combination is disclosed to all features disclosed in this specification (including adjoint claim, summary and accompanying drawing) and so to appoint
Where all processes or unit of method or equipment are combined.Unless expressly stated otherwise, this specification (including adjoint power
Profit is required, summary and accompanying drawing) disclosed in each feature can or similar purpose identical, equivalent by offer alternative features come generation
Replace.
Although in addition, it will be appreciated by those of skill in the art that some embodiments described herein include other embodiments
In included some features rather than further feature, but the combination of the feature of be the same as Example does not mean in of the invention
Within the scope of and form different embodiments.For example, in the following claims, times of embodiment claimed
One of meaning mode can be used in any combination.
In addition, be described as herein can be by the processor of computer system or by performing for some in the embodiment
Method or the combination of method element that other devices of the function are implemented.Therefore, with for implementing methods described or method
The processor of the necessary instruction of element forms the device for implementing this method or method element.In addition, device embodiment
Element described in this is the example of following device:The device is used to implement as in order to performed by implementing the element of the purpose of the invention
Function.
As used in this, unless specifically stated so, come using ordinal number " first ", " second ", " the 3rd " etc.
Description plain objects are merely representative of the different instances for being related to similar object, and are not intended to imply that the object being so described must
Must have the time it is upper, spatially, in terms of sequence or given order in any other manner.
Although describing the present invention according to the embodiment of limited quantity, above description, the art are benefited from
It is interior it is clear for the skilled person that in the scope of the present invention thus described, it can be envisaged that other embodiments.Additionally, it should be noted that
The language that is used in this specification primarily to readable and teaching purpose and select, rather than in order to explain or limit
Determine subject of the present invention and select.Therefore, in the case of without departing from the scope and spirit of the appended claims, for this
Many modifications and changes will be apparent from for the those of ordinary skill of technical field.For the scope of the present invention, to this
The done disclosure of invention is illustrative and not restrictive, and it is intended that the scope of the present invention be defined by the claims appended hereto.