CN103777188B - A kind of InISAR formation method based on the process of double frequency conjugation - Google Patents
A kind of InISAR formation method based on the process of double frequency conjugation Download PDFInfo
- Publication number
- CN103777188B CN103777188B CN201410026269.7A CN201410026269A CN103777188B CN 103777188 B CN103777188 B CN 103777188B CN 201410026269 A CN201410026269 A CN 201410026269A CN 103777188 B CN103777188 B CN 103777188B
- Authority
- CN
- China
- Prior art keywords
- target
- image
- imaging
- signal
- interferometric phase
- 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.)
- Expired - Fee Related
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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/583—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
-
- 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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
-
- 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
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9064—Inverse SAR [ISAR]
Abstract
Based on an InISAR formation method for double frequency conjugation process, describe one and double frequency conjugation treatment technology and interference technique are combined, be applied in ISAR, realize the method for moving target location and imaging, the present invention includes step: subband is split; Sub-band burst compression and the process of double frequency conjugation; Composite signal is utilized to interfere bigness scale angle to target; Original signal imaging and interference treatment; The result of interference of original signal image realizes target scattering point Superprecision Angle Measurement in conjunction with bigness scale angle result, and then realize target is located and imaging.The present invention mainly solves the problem that target Doppler frequency ambiguity and target interfere angle measurement fuzzy under InISAR pattern, can realize target location and imaging simultaneously.
Description
Technical field
The present invention relates to a kind of InISAR(Interferemetric InverseSynthetic Aperture Radar based on the process of double frequency conjugation) formation method, for obtaining true geometric size and the true spatial location of target, belong to radar imagery and Radar Signal Processing Technology field.
Background technology
Inverse synthetic aperture radar (ISAR) (Inverse Synthetic Aperture Radar, ISAR) range resolution is obtained by launching broadband signal, the relative motion between radar and target is utilized to form synthetic aperture, obtain Doppler by signal processing technology and tie up resolution, thus obtain target image.ISAR is the main radar system to motive target imaging, and from it is born, namely people have recognized that ISAR is imaged on the aspects such as astronomical sight, tactical weapon, strategic defensive and has significant application value.
Along with the continuous growth of application demand, the defect of conventional single-antenna ISAR imaging technique also embodies gradually.On the one hand, ISAR obtains target image mainly in distance-Doppler territory, can not obtain the space size information of target in transverse direction; On the other hand, single antenna ISAR image can not provide the drift angle information of target in wave beam, namely can not realize target orientation to location.Need the InISAR system considering to use multiple antennas to solve relevant issues for this reason.
On the basis of radar interference angle measurement technique, along with InSAR(interference synthetic aperture aperture radar) development of technology, the research work of InISAR technology have also been obtained to be paid attention to widely.1996, SoumekhM. propose a kind of for the InISAR method of ISAR image by tiny inclination angle during interference technique detection aircraft landing in Automatic aircraft landing using interferometric inverse synthetic apertureradar imaging mono-literary composition first, the various research about InISAR technology subsequently worldwide launches.Before and after 2000, the people such as the Wang Genyuan of Delaware university of the U.S. and the Victor C.Chen of Naval Research Labratory have studied and utilize three slave antennas forming vertical parallax to carry out the method for three-dimensional imaging (see Wang G to target, Xia X-G, Chen V C.Three-dimensionalISAR imaging of maneuvering targets using three receivers [J] .Image Processing, IEEE Transactions on, 2001).Along with the development of InISAR 3 Dimension Image Technique, the research for practical application request such as the estimation of moving target three-dimensional motion parameter, high precision imagings is also progressively goed deep into, and reports relevant progress in a series of periodical meeting.See with Publication about Document:
[1]Zhang Q,Yeo T S,Du G,et al.Estimation of three-dimensional motionparameters in interferometric ISAR imaging[J].Geoscience and RemoteSensing,IEEE Transactions on,2004,42(2):292-300.
[2]Zhang Q,Yeo T S.Three-dimensional SAR imaging of a ground movingtarget using the InISAR technique[J].Geoscience and Remote Sensing,IEEE Transactions on,2004,42(9):1818-1828.
[3]Given J A,Schmidt W R.Generalized ISAR-part II:interferometrictechniques for three-dimensional location of scatterers[J].ImageProcessing,IEEE Transactions on,2005
[4] Yin Jianfeng, Li Daojing, Wu Yirong. based on Space Object Detection and the imaging [J] of spaceborne millimetre-wave radar. aerospace journal, 2007,28 (6): 6.
In order to obtain the horizontal space dimensional information of target, the target rotation angle that No. CN101498788, Tsing-Hua University patent discloses a kind of inverse synthetic aperture radar (ISAR) is estimated and horizontal calibrating method, the method utilizes segmentation imaging results to be optimized combinational estimation target relative to the velocity of rotation of radar and relevant corner to scattering center, thus determines the breadth wise dimension of image.The method advantage does not adopt principle of interference, and the few system of receiving cable is simple.Shortcoming needs, to target imaging under short time data, therefore to require higher to echo signal to noise ratio (S/N ratio); Do not solve the problem of target localization, can not determine the absolute position of target in radar beam.
No. CN101000374, Xian Electronics Science and Technology University patent discloses a kind of InISAR formation method based on many special aobvious points, the method is carried out spy's aobvious some phase place main value unwrapping by screening aobvious of multiple spies of not glimmering and then calibrates the aobvious point of spy, simulate the relation of Doppler frequency and lateral separation, utilize this relation view picture ISAR to be obtained to the InISAR image of target true geometric size by pixel process.The method is calibrated based on the picture point after ISAR imaging, and require lower to echo signal to noise ratio (S/N ratio), in patent specification, the accompanying drawing of simulation result shows that the method can be applicable to the situation of low signal-to-noise ratio, and robustness is better.But, the method also has its significant defect: when to interferometric phase unwrapping, require according to Doppler's increasing or decreasing path integral unwrapping, but, different integration starting points may cause uncoiled result different, therefore, the actual value of interferometric phase that the method obtains is relatively accurate.That is, utilize this value to judge the distance between target scattering point, but the absolute position of each scattering point can not be obtained, namely can not to target localization.
The paper that double frequency conjugation treatment technology first time is seen in report one section " moving air target of satellite-borne SAR detects and imaging " by name that to be Yin Jianfeng deliver in " electronic letters, vol ".This article describes the method that subband double frequency conjugation ought to be used for moving target solution doppler ambiguity.The process of subband double frequency conjugation is a kind of difference frequency disposal route in essence, and the composite signal carrier frequency generated after process is less than the bandwidth of former linear FM signal, and compare former linear FM signal, carrier frequency reduces at least one the order of magnitude, and then reduces doppler centroid.But it is fuzzy fuzzy with doppler bandwidth that solution doppler ambiguity needs to solve doppler centroid, and this article does not have multianalysis simultaneously; Double frequency conjugation treatment technology can increase carrier wavelength, and this is highly beneficial to interferometric phase ambiguity solution, and this article does not deeply excavate.
Summary of the invention
Technology of the present invention is dealt with problems and is: overcome the deficiencies in the prior art, a kind of InISAR formation method based on the process of double frequency conjugation is provided, this method solves existing motive target imaging method can not to the problem of moving target location, optimize existing moving target size estimation method, and realize target is determined appearance and has been expanded the application of double frequency conjugation treatment technology in interference on the basis of location simultaneously.
Technical solution of the present invention is:
Based on an InISAR formation method for double frequency conjugation process, comprise the imaging processing of the Moving Target Return signal of two antenna channels, step is as follows:
Step S1: the Moving Target Return signal of an antenna channels is done Fourier transform to the fast time, m-fast frequency field when full time-domain signal is transformed into slow, and do in fast frequency field and divide equally, be divided into two subband linear FM signals;
Step S2: two subband linear FM signals in step S1 are made distance respectively to pulse compression, then does the process of double frequency conjugation and obtains composite signal, and described double frequency conjugation process is about to two distances to the subband signal conjugate multiplication after pulse compression;
Step S3: utilize ISAR imaging method to do imaging processing to the signal after synthesis, obtain the image of moving target in distance-Doppler territory, and enter step S5;
Step S4: the Moving Target Return signal of another one antenna channels is repeated to the signal transacting in step S1-S3, and enter step S5;
Step S5: utilize two ranges that obtain in step S3 and S4 to carry out radial velocity without blur estimation from the image of-Doppler domain to moving target, radial velocity is averaged without the result of blur estimation and is obtained target radial speed; The image in the distance-Doppler territory of obtain in step S3 and S4 two antenna channels done conjugate multiplication simultaneously and then extract scattering point interferometric phase, utilizing scattering point interferometric phase realize target drift angle bigness scale amount;
Step S6: utilize the original motion target echo signal of ISAR imaging method to two antenna channels to do imaging processing, then does conjugate multiplication to image after process, extracts scattering point interferometric phase;
Step S7: the result of the target drift angle bigness scale amount that the target image interferometric phase integrating step S5 utilizing step S6 to obtain obtains, carries out target drift angle accurate measurement amount, and then to target scattering point location, realizes InISAR imaging simultaneously.
In described step S5, the method for realize target drift angle bigness scale amount is as follows:
A, complete image registration: image after two antenna channels imagings is done interpolation processing respectively, translation is done in picture position after one of them antenna channels interpolation, image after another one antenna channels interpolation remains unchanged, calculate cross-correlation simultaneously, when cross-correlation is maximum, translational movement is translational movement needed for registration, then the interferometric phase image of down-sampled acquisition imaging region;
B, utilize interpolation after the imaging results of the composite signal of antenna channels that remains unchanged of image do threshold value threshold window, this window is utilized to do filtering to imaging region interferometric phase image, obtain target interferometric phase image, and utilize this target interferometric phase image realize target drift angle bigness scale amount:
Wherein, φ
nrepresent the interferometric phase of the n-th scattering point in target, λ
drepresent effective wavelength, θ is antenna beam scanning angle, d
12it is the distance of two antenna channels phase centers.
In described step S7, target scattering point location and target imaging complete simultaneously.
The present invention's beneficial effect is compared with prior art:
(1) the present invention requires lower to Moving Target Return signal to noise ratio (S/N ratio), can realize the bigness scale amount to target drift angle under Low SNR.
(2) the present invention can solve moving target doppler centroid and the fuzzy problem of doppler bandwidth simultaneously, can realize to target radial speed without blur estimation.
(3) the present invention adopts the method for thick smart two-stage angle measurement, the result of target drift angle bigness scale amount is utilized to provide true interferometric phase information initial value when accurately measuring, avoid the uncoiled work of conventional I nISAR interferometric phase, can realize the orientation of target to location.
Accompanying drawing explanation
Fig. 1 is formation method process flow diagram of the present invention;
Fig. 2 is imaging geometry model schematic of the present invention;
Fig. 3 is the subband signal generating principle figure of double frequency conjugation of the present invention process;
Fig. 4 is the target ISAR imaging results that the composite signal utilizing the process of double frequency conjugation to obtain obtains;
Fig. 5 is the target image inverting positioning result that obtains of InISAR formation method that the present invention introduces and simulation objectives actual position comparison diagram.
Embodiment
Specific embodiment below in conjunction with accompanying drawing and emulation experiment is further described in detail the specific embodiment of the present invention.
As shown in Figure 1, the concrete steps of a kind of InISAR formation method based on the process of double frequency conjugation of the present invention are as follows:
Step S1: the Moving Target Return signal of an antenna channels is done Fourier transform to the fast time, m-fast frequency field when full time-domain signal is transformed into slow, and do in fast frequency field and divide equally, be divided into two subband linear FM signals;
Step S2: two subband linear FM signals in step S1 are made distance respectively to pulse compression, then does the process of double frequency conjugation and obtains composite signal, and described double frequency conjugation process is about to two distances to the subband signal conjugate multiplication after pulse compression;
Step S3: utilize ISAR imaging method to do imaging processing to the signal after synthesis, obtain the image of moving target in distance-Doppler territory, and enter step S5;
Step S4: the Moving Target Return signal of another one antenna channels is repeated to the signal transacting in step S1-S3, and enter step S5;
Step S5: utilize two ranges that obtain in step S3 and S4 to carry out radial velocity without blur estimation from the image of-Doppler domain to moving target, radial velocity is averaged without the result of blur estimation and is obtained target radial speed; The image in the distance-Doppler territory of obtain in step S3 and S4 two antenna channels done conjugate multiplication simultaneously and then extract scattering point interferometric phase, utilizing scattering point interferometric phase realize target drift angle bigness scale amount;
Step S6: utilize the original motion target echo signal of ISAR imaging method to two antenna channels to do imaging processing, then does conjugate multiplication to image after process, extracts scattering point interferometric phase;
Step S7: the result of the target drift angle bigness scale amount that the target image interferometric phase integrating step S5 utilizing step S6 to obtain obtains, carries out target drift angle accurate measurement amount, and then to target scattering point location, realizes InISAR imaging simultaneously.
The course of work of the present invention is further illustrated below in conjunction with a specific embodiment:
Simulate signal is carrier frequency f
c=35GHz, bandwidth B
r=400MHz, time wide 2 μ s linear FM signal, system PRF is 4kHz, imaging geometry model as shown in Figure 2, θ=30 °, antenna beam scanning angle, the distance d of two antenna channels phase centers
12=3m, the synthetic aperture time is 0.25s, moving target radial velocity-61.6m/s, transverse velocity 78.8m/s, target forms " ten " font by 5 scattering points and distributes, lateral dimension 20m, radial dimension 10m, at imaging moment (intermediate time), target is 10km to the distance of antenna, drift angle in the wave beam of target's center's point
In described step S1, target echo is coupled fast time and slow time, and slow time-sampling point obtains the echo-pulse that a fast time series represents; Moving Target Return signal is done Fourier transform to the fast time and namely Fourier transform is done to each echo-pulse, m-fast frequency field when full time-domain signal will be transformed into slow.Each pulse is done in fast frequency field and divides equally, obtain height frequency two segment signal, and in relevant position respectively zero padding do inverse Fourier transform, obtain two subband linear FM signals, as shown in Figure 3.
In described step S2, the matched filter that the two subband linear FM signals adopting the method for matched filtering to build respectively to obtain with step S1 adapt, two subband linear FM signals are done distance respectively to pulse compression, and two signals are done the process of double frequency conjugation, by the subband signal conjugate multiplication after two pulse pressures, obtain the signal after difference frequency:
Wherein
represent low frequency sub-band,
represent high-frequency sub-band, " * " represents conjugation,
and t
mrepresent the speed time respectively.
This signal by each scattering point signal in target self and the polynomial expression that forms of cross term sum.Because the time difference of two subband linear FM signals is usually only at musec order, the displacement of moving target is far smaller than a Range resolution unit during this period, so can think that the Range Profile of moving target to aim in this two width image.Meanwhile, consider the character of sinc function, as long as the oblique distance difference of two scattering points is greater than a Range resolution unit, the amplitude of its cross term be just far smaller than self, cross term just can be ignored, especially when resolution is higher.
Based on above analysis, after the process of double frequency conjugation, composite signal can be expressed as
In formula, C represents the light velocity, and N represents the number of target scattering point, R
n(t
m) represent that the n-th scattering point is at t
mmoment to the distance of antenna, ξ
n1and ξ
n2represent first phase, the centre frequency of this signal is the difference of two subband linear FM signal centre frequencies, i.e. B
r/ 2, f
c2represent the carrier frequency of high-frequency sub-band linear FM signal, f
c1represent the carrier frequency of low frequency sub-band linear FM signal.The target Doppler centre frequency that then composite signal is corresponding is
The variation range of target Doppler bandwidth and target Doppler frequency within the synthetic aperture time.Target radial speed is not contributed doppler bandwidth, and the Doppler frequency difference that target lateral speed produces at imaging start/stop time is target Doppler bandwidth.Consider maximum Doppler bandwidth, namely known target maximum lateral speed is v
amax, the synthetic aperture time is T
s, then doppler bandwidth is
Wherein R represents that antenna arrives the reference distance of target, and λ represents the wavelength that carrier wave is corresponding.After the process of double frequency conjugation, the doppler bandwidth of composite signal is
visible, the doppler bandwidth of composite signal reduces greatly.
According to above-mentioned analysis, the parameter in conjunction with specific embodiments, target Doppler centre frequency is reduced to 82.1Hz by 14.4kHz; Target Doppler bandwidth is reduced to 0.2Hz by 37.1Hz.Visible, the process of double frequency conjugation makes the effective wavelength of composite signal greatly increase, and solves the doppler ambiguity problem of general motion target.
In described step S3, composite signal carrier frequency reduces, and doppler bandwidth reduces, and corresponding target image also significantly reduces in the scattered band of Doppler domain, and target ISAR image outline is lost, as shown in Figure 4.According to parameter in embodiment, target is reduced to 0.4Hz in the scattered band of Doppler domain by 73.5Hz, and therefore, target in the diagram its orientation can not identification to the position relationship between scattering point, and objective contour is lost.Distance is still preserved to scattering point information.
In described step S4, because composite signal effective wavelength is comparatively large, target Doppler frequency without fuzzy, Doppler frequency and corresponding target radial speed
Target Doppler frequency is 83Hz, and the target velocity that correspondence obtains is-62.25m/s, and actual value is relatively, achieve to target radial speed without blur estimation.
In described step S5, the method for realize target drift angle bigness scale amount is as follows:
A, complete image registration: image after two antenna channels imagings is done interpolation processing respectively, translation is done in picture position after one of them antenna channels interpolation, image after another one antenna channels interpolation remains unchanged, calculate cross-correlation simultaneously, when cross-correlation is maximum, translational movement is translational movement needed for registration, then the interferometric phase image of down-sampled acquisition imaging region;
B, utilize interpolation after the imaging results of the composite signal of antenna channels that remains unchanged of image do threshold value threshold window, this window is utilized to do filtering to imaging region interferometric phase image, obtain target interferometric phase image, and utilize this target interferometric phase image realize target drift angle bigness scale amount:
Wherein, φ
nrepresent the interferometric phase of the n-th scattering point in target, λ
drepresent effective wavelength, θ is antenna beam scanning angle, d
12it is the distance of two antenna channels phase centers.
In conjunction with specific embodiments, the process of double frequency conjugation makes wavelength greatly increase, consider θ=30 °, beam scanning angle, original signal angle measurement not fuzzy ranges expands to 14.4775 ° ~ 48.5904 ° from 29.9055 ° ~ 30.0946 °, solve and interfere angle measurement fuzzy problem, the result of target drift angle bigness scale amount is the average of all target scatterings point drift angle
In the present embodiment, threshold value thresholding is taken as-8dB, and corresponding n is 5.
In described step S6, the method for imaging and extraction interferometric phase can refer step S3 and step S5.
In described step S7, target localization and InISAR imaging need the result utilizing target drift angle bigness scale amount
true interferometric phase information initial value φ is provided,
corresponding true interferometric phase is known, and initial value φ is:
(k
1for positive integer)
Wherein, φ
0represent the main value of φ; In like manner, the n-th scattering point drift angle in target
corresponding interferometric phase actual value can be expressed as
(k
2for positive integer)
φ in formula
n0represent φ
nmain value.
Due to target scattering point drift angle
?
near, the difference not fuzzy ranges of the angle measurement more than both can thinking time remote.Therefore, can think that namely the difference of interferometric phase main value represents the difference of interferometric phase, in target, the interferometric phase actual value of the n-th scattering point can be asked
And then the smart measurement result of target scattering point drift angle can be tried to achieve:
Obtain the drift angle of each scattering point, in conjunction with the range information R that echo time delay characterizes
ntarget scattering point position can be finally inversed by and realize location:
The positioning instant of each scattering point of realize target achieves target imaging (namely the position being aware of each point defines target image) simultaneously, as shown in Figure 5.Give the contrast of inverting position and actual position in figure, location, distance 10km place mean square deviation is less than 2m, and illustrate that this method imaging effect is better, target location accuracy is higher.
Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; be not limited to the present invention; within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
The content be not described in detail in instructions of the present invention belongs to the known technology of those skilled in the art.
Claims (2)
1., based on an InISAR formation method for double frequency conjugation process, comprise the imaging processing of the Moving Target Return signal of two antenna channels, it is characterized in that step is as follows:
Step S1: the Moving Target Return signal of an antenna channels is done Fourier transform to the fast time, m-fast frequency field when full time-domain signal is transformed into slow, and do in fast frequency field and divide equally, be divided into two subband linear FM signals;
Step S2: two subband linear FM signals in step S1 are made distance respectively to pulse compression, then does the process of double frequency conjugation and obtains composite signal, and described double frequency conjugation process is about to two distances to the subband signal conjugate multiplication after pulse compression;
Step S3: utilize ISAR imaging method to do imaging processing to the signal after synthesis, obtain the image of moving target in distance-Doppler territory, and enter step S5;
Step S4: the Moving Target Return signal of another one antenna channels is repeated to the signal transacting in step S1-S3, and enter step S5;
Step S5: utilize two ranges that obtain in step S3 and S4 to carry out radial velocity without blur estimation to moving target respectively from the image of-Doppler domain, radial velocity is averaged without the result of blur estimation and is obtained target radial speed; The image in the distance-Doppler territory of obtain in step S3 and S4 two antenna channels done conjugate multiplication simultaneously and then extract scattering point interferometric phase, utilizing scattering point interferometric phase realize target drift angle bigness scale amount;
Step S6: utilize the original motion target echo signal of ISAR imaging method to two antenna channels to do imaging processing, then does conjugate multiplication to image after process, extracts scattering point interferometric phase;
Step S7: the result of the target drift angle bigness scale amount that the target image interferometric phase integrating step S5 utilizing step S6 to obtain obtains, carries out target drift angle accurate measurement amount, and then to target scattering point location, realizes InISAR imaging simultaneously;
In described step S5, the method for realize target drift angle bigness scale amount is as follows:
A, complete image registration: image after two antenna channels imagings is done interpolation processing respectively, translation is done in picture position after one of them antenna channels interpolation, image after another one antenna channels interpolation remains unchanged, calculate cross-correlation simultaneously, when cross-correlation is maximum, translational movement is translational movement needed for registration, then the interferometric phase image of down-sampled acquisition imaging region;
B, utilize interpolation after the imaging results of the composite signal of antenna channels that remains unchanged of image do threshold value threshold window, this window is utilized to do filtering to imaging region interferometric phase image, obtain target interferometric phase image, and utilize this target interferometric phase image realize target drift angle bigness scale amount:
Wherein, φ
nrepresent the interferometric phase of the n-th scattering point in target, λ
drepresent effective wavelength, θ is antenna beam scanning angle, d
12it is the distance of two antenna channels phase centers.
2. a kind of InISAR formation method based on the process of double frequency conjugation according to claim 1, is characterized in that: in described step S7, target scattering point location and target imaging complete simultaneously.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410026269.7A CN103777188B (en) | 2014-01-21 | 2014-01-21 | A kind of InISAR formation method based on the process of double frequency conjugation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410026269.7A CN103777188B (en) | 2014-01-21 | 2014-01-21 | A kind of InISAR formation method based on the process of double frequency conjugation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103777188A CN103777188A (en) | 2014-05-07 |
| CN103777188B true CN103777188B (en) | 2015-11-04 |
Family
ID=50569654
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410026269.7A Expired - Fee Related CN103777188B (en) | 2014-01-21 | 2014-01-21 | A kind of InISAR formation method based on the process of double frequency conjugation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103777188B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109031295A (en) * | 2018-07-17 | 2018-12-18 | 中国人民解放军国防科技大学 | ISAR image registration method based on wave path difference compensation |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104267378A (en) * | 2014-09-27 | 2015-01-07 | 励盼攀 | Interferometric phase simulation method |
| DE102018100632A1 (en) * | 2017-10-11 | 2019-04-11 | Symeo Gmbh | Radar method and system for determining the angular position, the location and / or the, in particular vectorial, speed of a target |
| CN107966688B (en) * | 2017-11-09 | 2021-04-20 | 东南大学 | Broadband radar target speed ambiguity resolving method based on phase interference technology |
| CN108802706B (en) * | 2018-06-19 | 2022-02-01 | 中国人民解放军63889部队 | Frequency modulation stepping radar signal target extraction method based on position calibration |
| CN111638516B (en) * | 2019-03-01 | 2022-02-18 | 中国科学院电子学研究所 | Terahertz frequency band SAR motion compensation algorithm based on double-frequency conjugate processing technology |
| CN111103585B (en) * | 2019-11-29 | 2023-05-26 | 西安电子科技大学 | Synthetic aperture broadband signal source reconnaissance imaging method based on double-channel combined processing |
| CN112558055B (en) * | 2020-11-13 | 2023-11-17 | 中国电子科技集团公司电子科学研究院 | Target positioning method, target positioning device, GMTI system and readable storage medium |
| CN113835089A (en) * | 2021-08-27 | 2021-12-24 | 中国空间技术研究院 | Slowly-varying elevation inversion method based on difference frequency InSAR |
| CN117289277B (en) * | 2023-11-27 | 2024-01-30 | 中山大学 | Multi-frequency radar three-dimensional imaging method and system based on subband segmentation synthesis |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4855747A (en) * | 1987-08-17 | 1989-08-08 | Trustees Of The University Of Pennsylvania | Method of target imaging and identification |
| US6023235A (en) * | 1997-06-05 | 2000-02-08 | Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. | Method for generating microwave-resolution images of moving objects by inverse synthetic aperture radar |
| CN101000374A (en) * | 2006-12-27 | 2007-07-18 | 西安电子科技大学 | Interference reverse synthetic aperture radarimaging method based on multiple-spot |
| CN103091682A (en) * | 2011-11-04 | 2013-05-08 | 中国科学院电子学研究所 | Interferometric inverse synthetic aperture radar (InISAR) hyperactivity target-oriented imaging and motion trail reconstruction method based on time frequency analysis |
-
2014
- 2014-01-21 CN CN201410026269.7A patent/CN103777188B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4855747A (en) * | 1987-08-17 | 1989-08-08 | Trustees Of The University Of Pennsylvania | Method of target imaging and identification |
| US6023235A (en) * | 1997-06-05 | 2000-02-08 | Deutsches Zentrum Fur Luft-Und Raumfahrt E.V. | Method for generating microwave-resolution images of moving objects by inverse synthetic aperture radar |
| CN101000374A (en) * | 2006-12-27 | 2007-07-18 | 西安电子科技大学 | Interference reverse synthetic aperture radarimaging method based on multiple-spot |
| CN103091682A (en) * | 2011-11-04 | 2013-05-08 | 中国科学院电子学研究所 | Interferometric inverse synthetic aperture radar (InISAR) hyperactivity target-oriented imaging and motion trail reconstruction method based on time frequency analysis |
Non-Patent Citations (2)
| Title |
|---|
| 一种高速运动目标长时间积累方法;张涛,吕彦伟,张琳,王晓东;《现代电子技术》;20130701;第36卷(第13期);28-30 * |
| 基于毫米波InISAR成像的运动目标探测与定位;刘波； 潘舟浩； 李道京； 乔明;《红外与毫米波学报》;20120630;第31卷(第3期);258-264 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109031295A (en) * | 2018-07-17 | 2018-12-18 | 中国人民解放军国防科技大学 | ISAR image registration method based on wave path difference compensation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103777188A (en) | 2014-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103777188B (en) | A kind of InISAR formation method based on the process of double frequency conjugation | |
| Martorella et al. | 3D interferometric ISAR imaging of noncooperative targets | |
| Xu et al. | Three-dimensional interferometric ISAR imaging for target scattering diagnosis and modeling | |
| CN103487809B (en) | A kind of based on BP algorithm and time become the airborne InSAR data disposal route of baseline | |
| CN107085212B (en) | Spin target time-varying three-dimensional imaging method based on linear frequency modulation stepping signal | |
| CN104808204A (en) | Moving-target detecting method and imaging method of stationary transmitter bistatic forward-looking synthetic aperture radar (SAR) | |
| Guo et al. | Geometry reconstruction based on attributes of scattering centers by using time-frequency representations | |
| CN103675816B (en) | A kind of interference inverse synthetic aperture radar imaging method based on compressive sensing theory | |
| KR102142674B1 (en) | Method and Apparatus for Synthetic Aperture Radar Phase Unwrapping based on SAR Offset Tracking Displacement Model | |
| Li et al. | A coarse-to-fine autofocus approach for very high-resolution airborne stripmap SAR imagery | |
| KR20110060626A (en) | The method for measurign object's velocity using synthetic aperture radar image and the apparatus thereof | |
| CN103630903B (en) | The method of flow field, sea radial velocity is measured based on straight rail interference SAR | |
| CN103792535A (en) | Method for measuring TEC value of ionized layer by utilizing SAR satellite | |
| CN103091682B (en) | Interferometric inverse synthetic aperture radar (InISAR) hyperactivity target-oriented imaging and motion trail reconstruction method based on time frequency analysis | |
| CN103809180A (en) | Azimuth pre-filtering processing method for Interferometric Synthetic Aperture Radar (InSAR) topographic survey | |
| CN103885062A (en) | Double-base foresight SAR moving target imaging method and moving target speed estimation method | |
| Felguera-Martin et al. | Interferometric ISAR imaging on maritime target applications: Simulation of realistic targets and dynamics | |
| CN104459634A (en) | Interferometric phase truth value computing method for UWB InSAR | |
| Xu et al. | A variable PRF imaging method for high squint diving SAR | |
| Li et al. | An advanced DSS-SAR InSAR terrain height estimation approach based on baseline decoupling | |
| Noviello et al. | ISAR motion compensation based on a new Doppler parameters estimation procedure | |
| Shi et al. | Study on the signature of ground moving target for airborne squint SAR imaging | |
| Qi et al. | Application of microwave imaging in regional deformation monitoring using ground based SAR | |
| CN105842695A (en) | Squint InSAR moving-static target separation method based on mixed baselines | |
| Magnard et al. | High resolution MEMPHIS SAR data processing and applications |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151104 Termination date: 20170121 |