CN105044719A - Terahertz high-precision vertical curved surface imaging method based on circumference SAR - Google Patents
Terahertz high-precision vertical curved surface imaging method based on circumference SAR Download PDFInfo
- Publication number
- CN105044719A CN105044719A CN201510347185.8A CN201510347185A CN105044719A CN 105044719 A CN105044719 A CN 105044719A CN 201510347185 A CN201510347185 A CN 201510347185A CN 105044719 A CN105044719 A CN 105044719A
- Authority
- CN
- China
- Prior art keywords
- theta
- signal
- radar
- max
- frequency domain
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 claims abstract description 48
- 238000006467 substitution reaction Methods 0.000 claims abstract description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 230000002123 temporal effect Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 2
- 230000006835 compression Effects 0.000 abstract 1
- 238000007906 compression Methods 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 2
- 208000004350 Strabismus Diseases 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/904—SAR modes
- G01S13/9088—Circular SAR [CSAR, C-SAR]
-
- 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
-
- 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/9004—SAR image acquisition techniques
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)
Abstract
The invention belongs to the synthetic aperture radar (SAR) imaging processing technology field and provides a circumference SAR curved surface imaging method based on a terahertz wave band. According to a characteristic of a wave number domain method, a circumference SAR is combined so that imaging of a point target on a curved surface is realized. Through analyzing a circumference SAR model, a distance formula r (ralTn) after approximation is acquired. An echo signal s (theta, k) after distance compression is directly established in a fast time wave number domain. Then Fourier transform about an azimuth variable theta is performed on the echo signal so that a two-dimension wave number domain signal s (k theta, k) is acquired. Then variable substitution is performed and a complex radical item is eliminated so that inverse Fourier transform can be performed in form. Meanwhile, interpolation on a ky dimension is performed on a signal s (kx, ky) after substitution so that rectangularization is performed on a data lattice of the echo signal. And the signal can smoothly carry out two-dimensional Fourier inverse transform so that an original point target can be recovered.
Description
Technical field
The invention belongs to synthetic-aperture radar (SAR) Imaging processing techniques field, namely imaging technique Analysis of Radar echo model is used, therefrom recover the process of point target, be specifically related to the vertical curved surface formation method of the synthetic-aperture radar (SAR) under nearly terahertz wave band.
Background technology
Terahertz Technology is the focus of academic circles at present, all receives showing great attention to of countries in the world research institution in various field.Terahertz Technology is applied to quality and the degree of accuracy that detection in target detection, will be improved greatly.The first, the carrier frequency of Terahertz radar is high, can the high signal of transmitted bandwidth, the range resolution of radar is increased, thus can identify more tiny object.Secondly, Terahertz frequency range can provide extremely narrow antenna beam, thus substantially increases the azimuth resolution of radar, also can detect for the indiscoverable Stealth weapon of normal radar.3rd, the penetrability of THz wave is strong, also has good effect to the detection of space high-speed moving object.
Synthetic-aperture radar (SAR) under circumference mode pore diameter, relative to the synthetic-aperture radar under general mode, it can carry out the observation of 360 degree to target, therefore can carry out the observation of longer time to scene, thus obtain higher resolution.And the more important thing is, be positioned at differing heights face for synchronization and identical two targets of oblique distance, along with the movement of radar, the oblique distance of these two targets can become no longer identical.Therefore, circumference SAR can utilize such change obtain object height to information, thus there is the ability of three-dimensional imaging.
The wavenumber domain method of synthetic-aperture radar (SAR) and traditional temporal frequency (t, f) method such as range Doppler (R-D) method in territory is different, the echoed signal of temporal frequency domain is transformed to distance wave number (x, k by wavenumber domain method
x) territory, both have certain duality relation.But time quantum is one dimension, represents with scalar, and amount of space is multidimensional, uses vector representation.The more important thing is, in slow time domain and spatial domain, the echoed signal received by a certain moment (place) is from all directions in wave beam; Doppler domain and wavenumber domain then different, the echo corresponding to a certain frequency not from synchronization, but from same direction.Therefore, regardless of the angle of squint of radar and the size of scene, method under wavenumber domain can both not add the focusing effect of other approximate conditions realization without geometric deformation to whole region based on scatter times, no matter be below the far field of microwave region, or under the near field pattern of terahertz wave band.In principle, it is the optimum realization of SAR imaging.
Circumference SAR vertical curved surface formation method before this, circumference SAR for terahertz wave band does not almost relate to: 2009, FerraraM, JacksonJA, propose with AustinC circumference SAR pattern of navigating more, namely carrier of radar platform at various height plane Circular test observation is carried out to same observation area.This pattern is carried out under microwave region, this pattern by be increased in height to sampling improve imaging system in the resolution perpendicular to direction of visual lines, strengthen the three-dimensional imaging ability of circumference SAR system.In theory, this method can increase radar in image quality highly upwards.But in the model of Terahertz radar curved planar reformation, because sensor can only carry out DATA REASONING round a fixing height track, can not plane take multiple measurements at various height, this is the shortcoming of this method.For traditional rear orientation projection's formation method, its resolution have received great restriction, if use nearly Terahertz frequency range as carrier frequency, so the operand of method and arithmetic speed will be slowly, therefore in the ordinary course of things, we can not use rear orientation projection's method under nearly Terahertz frequency range; If transfer to make in this way under microwave frequency band, resolution situation high not enough can be faced again.Above-mentioned content all can cause the difficulty of this classic method on application.
Summary of the invention
The object of the invention is, for overcoming existing circumference SAR for the detection of target on vertical curved surface face and imaging not accurately and clearly shortcoming, specially to provide a kind of circumference SAR curved planar reformation method based on terahertz wave band.The basic ideas of the method are the wavenumber domain interpolation methods utilizing synthetic-aperture radar, the circumferential synthetic aperture data of wide bandwidth are processed, and by imaging results reconstruct on curved surface.
Detailed technology scheme of the present invention is as follows:
Based on a Terahertz high precision vertical curved surface formation method of circumference SAR, comprise the following steps:
Step 1: the circumference SAR high precision vertical curved surface formation method geometric model of structure terahertz wave band:
The flying height of radar is H, take radius as R
gcircular arc make uniform circular motion, radar position position angle be θ, so θ ∈ (-θ
max, θ
max) represent the angular range of synthetic aperture focusing target under the slow time;
It is r that set-point target is all positioned at radius
0cambered surface on, and the position angle of the n-th target is
be highly z
n;
Under this model, use the conversion formula of polar coordinates and rectangular coordinate, the positional representation of radar is: r
a=(R
gcos θ, R
gsin θ, H), the impact point on reconstructed surface region is expressed as:
Step 2: calculate the distance of radar to impact point:
Abbreviation is:
Wherein, m=(r
0-R
g)
2+ (H-z
n)
2, n=R
gr
0;
Step 3: the range formula calculated according to step 2, derivation point target is spatial frequency domain signal under polar coordinates:
Setting radar transmit as linear FM signal is:
p(t)=rect(t/T)exp(-jπk
rt
2);
Wherein,
t is the fast time, and T is the pulse persistance cycle, k
rit is frequency modulation rate;
When transmitting through distance r (r
a| T
n) and after returning, the echo complex baseband signal that radar receives is:
s(t,θ)=rect([t-2(r(r
a|T
n)-r
ref)/c]/T)exp(-jπk
r[t-2(r(r
a|T
n)-r
ref)/c]
2)·exp(-j4πfr(r
a|T
n)/c);
Wherein, f is fast temporal frequency, and c is the light velocity, r
reffor reference distance;
The method of matched filtering is used to carry out Range compress to signal above, the Received signal strength after obtaining compressing:
s(r
a,k)=δ(t-2(r(r
a|T
n)-r
ref)/c)exp(-j2kr(r
a|T
n));
Wherein, k=2 π f/c is called as fast time wave number, and scope is k ∈ [2 π f
min/ c, 2 π f
max/ c]; f
minwith f
maxfor minimum frequency and the maximum frequency of linear FM signal;
Brought into by range formula, the Received signal strength obtaining scattering point is:
Above formula is also referred to as the spatial frequency domain signal of target under polar coordinates;
Step 4: work about azimuth angle theta is done about the Fourier transform of azimuth angle theta to step 3 gained spatial frequency domain signal:
The site in phasing of above-mentioned azimuth Fourier transform integration is obtained according to theorem in phase bit:
Substituted in expression formula in phase bit, obtained the result that azimuthal θ asks Fourier transform:
Wherein, k
θbe called as position angle frequency domain wave number, above formula is called as the two-dimensional frequency signal of position angle frequency domain-fast time frequency domain;
Step 5: computer azimuth angle frequency domain k
θexpression formula and scope;
In the two-dimensional frequency signal expression that step 4 obtains, the scope of k is k ∈ [2 π f
min/ c, 2 π f
max/ c];
First the expression formula of computer azimuth angle frequency domain:
Calculate k again
θscope:
Wherein, k
max=2 π f
max/ c is called as maximum wave number, and θ is the position angle of radar position;
Step 6: substitution of variable is done to two-dimensional frequency signal as follows:
The k that step 4 is obtained
θthe wavenumber spectrum of-k transforms to k
x-k
yplane:
Wherein, k
x, k
ybe respectively the wave number of imaging region transverse axis, the longitudinal axis;
Step 7: to converting the k obtained
x-k
ythe net point of plane carries out k
ydimension interpolation;
Work as k
θ=k
xwhen=0, k
y=2k is that interpolation exports ordinate;
Draw the size of interpolation amount:
Step 8: carry out two-dimentional inverse Fourier transform to the signal of gained after interpolation, obtains the original position distribution of final target, exports final image.
So core of the present invention is proposing the vertical curved surface formation method of a kind of circumference SAR under nearly terahertz wave band (85GHz-105GHz), because Terahertz wave frequency is high, so some the traditional formation methods under microwave region just cannot use at terahertz wave band, the present invention, according to the characteristic of wavenumber domain method, achieves the imaging for point target on curved surface in conjunction with circumference SAR.By the analysis of circumference SAR model is obtained approximate after range formula r (r
a| T
n), and the echoed signal s (θ, k) after fast Time Wave number field directly sets up Range compress, then this echoed signal is done about the Fourier transform of orientation to variable θ, obtain two-dimentional wavenumber domain signal s (k
θ, k); Next substitution of variable is carried out to it, eliminate complicated radical item, thus allow it carry out inverse Fourier transform in form; Meanwhile, to the signal s (k after replacement
x, k
y) carry out k
ydimension interpolation above, makes the data dot matrix of echoed signal squaring, enables signal carry out two-dimentional inverse Fourier transform smoothly thus recover original point target.The maximum feature of the present invention is that the vertical curved surface formation method based on circumference SAR proposed, under working in nearly terahertz wave band, can significantly improve the degree of accuracy of target detection and the sharpness of target imaging.
In addition, method provided by the present invention is by the regional restructuring of target imaging on curved surface, has better imaging effect for the target on vertical curved surface; Because the structure of human body is similar to its structure, the present invention is also for the safety check at airport, station provides basic theoretical foundation.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of a kind of Terahertz high precision vertical curved surface formation method based on circumference SAR of the present invention.
Fig. 2 is the geometric model of a kind of Terahertz high precision vertical curved surface formation method based on circumference SAR of the present invention.
Fig. 3-1 is echoed signal schematic diagram in the Support of two-dimentional wavenumber domain before substitution of variable.
Fig. 3-2 is after substitution of variable, and echoed signal is at k
x-k
yplanar interpolation obtains the schematic diagram of rectangular node point data.
Embodiment
Below in conjunction with accompanying drawing and embodiment, the present invention is described further.
Based on a Terahertz high precision vertical curved surface formation method process flow diagram of circumference SAR, as shown in Figure 1, comprise the following steps:
Step 1: the circumference SAR high precision vertical curved surface formation method geometric model of structure terahertz wave band:
Based on a Terahertz high precision vertical curved surface formation method geometric model of circumference SAR, as shown in Figure 2:
The flying height of radar is H, take radius as R
gcircular arc make uniform circular motion; Radar position position angle be θ, so θ ∈ (-θ
max, θ
max) represent the angular range of synthetic aperture focusing target under the slow time;
It is r that postulated point target is all positioned at radius
0cambered surface on, and the position angle of the n-th target is
be highly z
n;
Under this model, use the conversion formula of polar coordinates and rectangular coordinate, the position of radar can be expressed as r
a=(R
gcos θ, R
gsin θ, H), for the target on reconstructed surface region,
Step 2: calculate the distance of radar to impact point:
Taylor's formula form according to cosx:
Range formula can abbreviation be:
Wherein, m=(r
0-R
g)
2+ (H-z
n)
2, n=R
gr
0;
Step 3: the approximate distance formula calculated according to step 2, derive point target spatial frequency domain signal under polar coordinates:
That supposes radar transmits as linear FM signal:
p(t)=rect(t/T)exp(-jπk
rt
2);
Wherein,
t is the fast time, and T is the pulse persistance cycle, k
rit is frequency modulation rate;
When transmitting through distance r (r
a| T
n) and after returning, the echo complex baseband signal that radar receives can be obtained:
s(t,θ)
=rect([t-2(r(r
a|T
n)-r
ref)/c]/T)exp(-jπk
r[t-2(r(r
a|T
n)-r
ref)/c]
2)·exp(-j4πfr(r
a|T
n)/c);
Wherein, f is fast temporal frequency, and c is the light velocity, r
reffor reference distance, refer to that radar arrives the distance of scene central point here;
The method of matched filtering is used to carry out Range compress to signal above, the Received signal strength after obtaining compressing:
s(r
a,k)=δ(t-2(r(r
a|T
n)-r
ref)/c)exp(-j2kr(r
a|T
n));
Wherein, k=2 π f/c is called as fast time wave number, and its scope is k ∈ [2 π f
min/ c, 2 π f
max/ c]; The Support of k determined by system bandwidth B, B=f
max-f
min, f
minwith f
maxbe called as minimum frequency and the maximum frequency of linear FM signal;
Brought into by range formula after abbreviation, the Received signal strength that can obtain scattering point is:
Above formula is also referred to as the spatial frequency domain signal of target under polar coordinates;
Step 4: the position angle-fast time frequency-region signal next drawn step 3 makes the Fourier transform about azimuth angle theta:
Common conclusions according to theorem in phase bit:
X in above formula
*∈ x is called as site in phasing, by g ' (x
*site in phasing can be obtained in)=0;
The site in phasing of above-mentioned azimuth Fourier transform integration can be obtained thus, as follows:
Substituted in the expression formula of staying phase bit, the approximation that azimuthal asks Fourier transform can be obtained:
Wherein, k
θbe called as position angle frequency domain wave number; Above formula is called as the two-dimensional frequency signal of position angle frequency domain-fast time frequency domain;
Step 5: calculate position angle frequency domain k
θexpression formula, estimate its scope, and draw the relation of itself and fast time frequency domain k, in two-dimentional wave number plane coordinate system, draw k
θthe relation of-k;
In the two-dimensional frequency expression formula that step 4 obtains, the scope of k is k ∈ [2 π f
min/ c, 2 π f
max/ c];
Before application interpolation, need the Support calculating angle frequency-region signal, first need the expression formula of computer azimuth angle frequency domain:
K
θthe scope of Support is determined by the azimuth coverage of reconstruction region and the length in aperture, its scope of the estimation that can be similar to thus:
Wherein, k
max=2 π f
max/ c is called as maximum wave number, and θ is the position angle of radar position;
K can be known by above analysis
θbe proportional with k, this relation can be drawn in two-dimentional wave number plane coordinate system, as shown in figure 3-1;
Step 6: substitution of variable is done to two-dimensional frequency signal as follows:
According to the character of wavenumber domain, in analytical model before and after substitution of variable in wavenumber domain, the relation between sampled point, reaches a conclusion and draws net point distribution plan;
The k that step 4 is obtained
θthe wavenumber spectrum of-k transforms to k
x-k
yplane:
Wherein, k
x, k
ybe respectively the wave number of imaging region transverse axis, the longitudinal axis;
Can see, transform to k from k
yterritory is nonlinear transformation, and the data point pointwise in Fig. 3-1 is transformed to k
x-k
yplane, now net point distribution is the solid section in Fig. 3-2; Data point is now elliptic systems;
Step 7: to converting the k obtained
x-k
ythe net point of plane carries out k
ydimension interpolation;
Due to the k that step 6 obtains
x-k
yplane net lattice point is elliptic systems, if we want the original position of recovery point objectives, smooth imaging, must do two-dimentional inverse Fourier transform to obtained signal, and such signal must meet the distribution of rectangular node point data; Therefore, we need to pass through k
ythe interpolation of dimension obtains the rectangle net grid point distribution as Fig. 3-2 hollow parts;
Work as k
θ=k
xwhen=0, k
y=2k is that required interpolation exports ordinate position;
According to Fig. 3-2, we can draw the size of interpolation amount:
The signal obtained after only having interpolation can reflect the time-frequency characteristic of target echo accurately;
Step 8: carry out two-dimentional inverse Fourier transform to the signal of gained after interpolation, obtains the original position distribution of final target; Export final image.
Point target being placed on radius in emulation experiment is r
0on the vertical curved surface of=20cm, test effect of the present invention, suppose that carrier frequency is about 0.1THz, the position angle of reconstruction region, from-50 ° to 50 °, is highly 0 to 0.5m.Radar is 0.5m at radius, highly for uniform circular motion made by the track of 2m.To echo according to step tectonic model of the present invention, find in time domain, signal is the distribution of several sinusoidal curves, shows the correctness of model.At two-dimensional frequency domain, the distribution of echo signal data point is in trapezoidal, consistent with theoretical analysis.
Adopt method of the present invention to make interpolation to two-dimensional frequency signal, the data point of the rectangular distribution after obtaining substitution of variable, more two-dimentional inverse Fourier transform is done to the signal after process, just obtain the point target of reconstruction.Can see, the point target of reconstruction is distributed on curved surface, and coordinate distributes consistent with the coordinate of hypothesis.This illustrates that the present invention has good imaging effect to target.
Compared with traditional rear orientation projection method, this method more has superiority in arithmetic speed.This is because rear orientation projection's method needs to carry out imaging scene the division of grid, and the cumulative realization that will be concerned with to each orientation data is upwards to the imaging of scene objects, and all to need from echo data interpolation to estimate the data of this net point for each net point.This will inevitably cause increase significantly its operation time.And on imaging effect, method used in the present invention is also better than rear orientation projection's method, more superior from the effect focused on.
The above, be only the specific embodiment of the present invention, arbitrary feature disclosed in this specification, unless specifically stated otherwise, all can be replaced by other equivalences or the alternative features with similar object; Step in disclosed all features or all methods or process, except mutually exclusive feature and/or step, all can be combined in any way.
Claims (1)
1., based on a Terahertz high precision vertical curved surface formation method of circumference SAR, comprise the following steps:
Step 1: the circumference SAR high precision vertical curved surface formation method geometric model of structure terahertz wave band:
The flying height of radar is H, take radius as R
gcircular arc make uniform circular motion, radar position position angle be θ, so θ ∈ (-θ
max, θ
max) represent the angular range of synthetic aperture focusing target under the slow time;
It is r that set-point target is all positioned at radius
0cambered surface on, and the position angle of the n-th target is
be highly z
n;
Under this model, use the conversion formula of polar coordinates and rectangular coordinate, the positional representation of radar is: r
a=(R
gcos θ, R
gsin θ, H), the impact point on reconstructed surface region is expressed as:
Step 2: calculate the distance of radar to impact point:
Abbreviation is:
Wherein, m=(r
0-R
g)
2+ (H-z
n)
2, n=R
gr
0;
Step 3: the range formula calculated according to step 2, derivation point target is spatial frequency domain signal under polar coordinates:
Setting radar transmit as linear FM signal is:
p(t)=rect(t/T)exp(-jπk
rt
2);
Wherein,
t is the fast time, and T is the pulse persistance cycle, k
rit is frequency modulation rate;
When transmitting through distance r (r
a| T
n) and after returning, the echo complex baseband signal that radar receives is:
s(t,θ)=rect([t-2(r(r
a|T
n)-r
ref)/c]/T)exp(-jπk
r[t-2(r(r
a|T
n)-r
ref)/c]
2)·exp(-j4πfr(r
a|T
n)/c);
Wherein, f is fast temporal frequency, and c is the light velocity, r
reffor reference distance;
The method of matched filtering is used to carry out Range compress to signal above, the Received signal strength after obtaining compressing:
s(r
a,k)=δ(t-2(r(r
a|T
n)-r
ref)/c)exp(-j2kr(r
a|T
n));
Wherein, k=2 π f/c is called as fast time wave number, and scope is k ∈ [2 π f
min/ c, 2 π f
max/ c]; f
minwith f
maxfor minimum frequency and the maximum frequency of linear FM signal;
Brought into by range formula, the Received signal strength obtaining scattering point is:
Above formula is also referred to as the spatial frequency domain signal of target under polar coordinates;
Step 4: work about azimuth angle theta is done about the Fourier transform of azimuth angle theta to step 3 gained spatial frequency domain signal:
The site in phasing of above-mentioned azimuth Fourier transform integration is obtained according to theorem in phase bit:
Substituted in expression formula in phase bit, obtained the result that azimuthal θ asks Fourier transform:
Wherein, k
θbe called as position angle frequency domain wave number, above formula is called as the two-dimensional frequency signal of position angle frequency domain-fast time frequency domain;
Step 5: computer azimuth angle frequency domain k
θexpression formula and scope;
In the two-dimensional frequency signal expression that step 4 obtains, the scope of k is k ∈ [2 π f
min/ c, 2 π f
max/ c];
First the expression formula of computer azimuth angle frequency domain:
Calculate k again
θscope:
Wherein, k
max=2 π f
max/ c is called as maximum wave number, and θ is the position angle of radar position;
Step 6: substitution of variable is done to two-dimensional frequency signal as follows:
The k that step 4 is obtained
θthe wavenumber spectrum of-k transforms to k
x-k
yplane:
Wherein, k
x, k
ybe respectively the wave number of imaging region transverse axis, the longitudinal axis;
Step 7: to converting the k obtained
x-k
ythe net point of plane carries out k
ydimension interpolation;
Work as k
θ=k
xwhen=0, k
y=2k is that interpolation exports ordinate;
Draw the size of interpolation amount:
Step 8: carry out two-dimentional inverse Fourier transform to the signal of gained after interpolation, obtains the original position distribution of final target, exports final image.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510347185.8A CN105044719B (en) | 2015-06-23 | 2015-06-23 | A kind of high-precision vertical curved surface imaging method of the Terahertz based on circumference SAR |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510347185.8A CN105044719B (en) | 2015-06-23 | 2015-06-23 | A kind of high-precision vertical curved surface imaging method of the Terahertz based on circumference SAR |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105044719A true CN105044719A (en) | 2015-11-11 |
CN105044719B CN105044719B (en) | 2017-10-17 |
Family
ID=54451394
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510347185.8A Expired - Fee Related CN105044719B (en) | 2015-06-23 | 2015-06-23 | A kind of high-precision vertical curved surface imaging method of the Terahertz based on circumference SAR |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105044719B (en) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105676223A (en) * | 2016-01-25 | 2016-06-15 | 电子科技大学 | SAR (Synthetic Aperture Radar) imaging method of non-uniform-speed terahertz radar platform |
CN105785368A (en) * | 2016-04-29 | 2016-07-20 | 电子科技大学 | Hidden hazardous material detecting method based on video SAR |
CN106772380A (en) * | 2017-03-31 | 2017-05-31 | 电子科技大学 | A kind of circumferential synthetic aperture radar imaging method |
CN106908781A (en) * | 2017-02-27 | 2017-06-30 | 中国科学院电子学研究所 | The velocity acquisition methods of the linear motion moving-target based on single channel circular track SAR |
CN106950555A (en) * | 2017-05-03 | 2017-07-14 | 中国人民解放军国防科学技术大学 | A kind of Area Objects imaging method based on Terahertz aperture coded imaging system |
CN107255815A (en) * | 2017-05-19 | 2017-10-17 | 北京空间飞行器总体设计部 | A kind of target surface reconstructing method based on bistatic scattering center time-frequency characteristics |
WO2017223386A1 (en) * | 2016-06-22 | 2017-12-28 | Duke University | Multiple-input-multiple-output (mimo) imaging systems and methods for performing massively parallel computation |
CN107730463A (en) * | 2017-09-30 | 2018-02-23 | 西安电子科技大学 | SAR image Enhancement Method based on the empty sunken copped wave of mutative scale and background filtering |
CN108872985A (en) * | 2018-04-10 | 2018-11-23 | 西北工业大学 | A kind of near field circumference SAR rapid three dimensional imaging process |
CN109143239A (en) * | 2018-10-12 | 2019-01-04 | 黑龙江八农垦大学 | A kind of imaging method of the circumferential synthetic aperture radar based on one-dimensional range profile |
CN109557541A (en) * | 2018-12-17 | 2019-04-02 | 中国人民解放军国防科技大学 | Holographic penetration imaging radar polar coordinate data processing method |
CN109633636A (en) * | 2018-10-29 | 2019-04-16 | 上海无线电设备研究所 | A kind of kinematic error is to Terahertz circular track SAR image quality impact analysis method |
CN110554387A (en) * | 2019-09-04 | 2019-12-10 | 中国科学院电子学研究所 | synthetic aperture interference near-field active source imaging method and device |
CN111220979A (en) * | 2020-01-16 | 2020-06-02 | 电子科技大学 | Imaging method for curved synthetic aperture radar |
CN113239773A (en) * | 2021-05-08 | 2021-08-10 | 重庆邮电大学 | Human body posture non-contact sensing method based on CSAR three-dimensional imaging |
CN115982571A (en) * | 2023-01-12 | 2023-04-18 | 中国人民解放军军事科学院系统工程研究院 | Geometric representation method and device for space terahertz communication signals |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070114418A1 (en) * | 2005-09-20 | 2007-05-24 | Mueller Eric R | Security portal with THz trans-receiver |
CN102313885A (en) * | 2010-06-30 | 2012-01-11 | 中国科学院电子学研究所 | Multi-dimensional microwave imaging system and method |
-
2015
- 2015-06-23 CN CN201510347185.8A patent/CN105044719B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070114418A1 (en) * | 2005-09-20 | 2007-05-24 | Mueller Eric R | Security portal with THz trans-receiver |
CN102313885A (en) * | 2010-06-30 | 2012-01-11 | 中国科学院电子学研究所 | Multi-dimensional microwave imaging system and method |
Non-Patent Citations (4)
Title |
---|
张彪 等: "采用格林函数分解的太赫兹逆合成孔径雷达近场成像算法", 《信号处理》 * |
李晋 等: "太赫兹频段目标微多普勒信号特征分析", 《电子测量与仪器学报》 * |
闵锐: "机载SAR三维成像理论及关键技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
魏明贵 等: "太赫兹时域雷达成像研究", 《雷达学报》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105676223A (en) * | 2016-01-25 | 2016-06-15 | 电子科技大学 | SAR (Synthetic Aperture Radar) imaging method of non-uniform-speed terahertz radar platform |
CN105785368A (en) * | 2016-04-29 | 2016-07-20 | 电子科技大学 | Hidden hazardous material detecting method based on video SAR |
WO2017223386A1 (en) * | 2016-06-22 | 2017-12-28 | Duke University | Multiple-input-multiple-output (mimo) imaging systems and methods for performing massively parallel computation |
CN106908781A (en) * | 2017-02-27 | 2017-06-30 | 中国科学院电子学研究所 | The velocity acquisition methods of the linear motion moving-target based on single channel circular track SAR |
CN106908781B (en) * | 2017-02-27 | 2019-10-18 | 中国科学院电子学研究所 | The velocity vector acquisition methods of linear motion moving-target based on single channel circular track SAR |
CN106772380A (en) * | 2017-03-31 | 2017-05-31 | 电子科技大学 | A kind of circumferential synthetic aperture radar imaging method |
CN106950555B (en) * | 2017-05-03 | 2019-05-17 | 中国人民解放军国防科学技术大学 | A kind of Area Objects imaging method based on Terahertz aperture coded imaging system |
CN106950555A (en) * | 2017-05-03 | 2017-07-14 | 中国人民解放军国防科学技术大学 | A kind of Area Objects imaging method based on Terahertz aperture coded imaging system |
CN107255815A (en) * | 2017-05-19 | 2017-10-17 | 北京空间飞行器总体设计部 | A kind of target surface reconstructing method based on bistatic scattering center time-frequency characteristics |
CN107730463A (en) * | 2017-09-30 | 2018-02-23 | 西安电子科技大学 | SAR image Enhancement Method based on the empty sunken copped wave of mutative scale and background filtering |
CN107730463B (en) * | 2017-09-30 | 2021-02-12 | 西安电子科技大学 | SAR image enhancement method based on variable-scale space-trapped chopping and background filtering |
CN108872985A (en) * | 2018-04-10 | 2018-11-23 | 西北工业大学 | A kind of near field circumference SAR rapid three dimensional imaging process |
CN108872985B (en) * | 2018-04-10 | 2022-04-05 | 西北工业大学 | Near-field circumference SAR rapid three-dimensional imaging method |
CN109143239A (en) * | 2018-10-12 | 2019-01-04 | 黑龙江八农垦大学 | A kind of imaging method of the circumferential synthetic aperture radar based on one-dimensional range profile |
CN109143239B (en) * | 2018-10-12 | 2022-05-20 | 黑龙江八一农垦大学 | Imaging method of circumferential synthetic aperture radar based on one-dimensional range profile |
CN109633636A (en) * | 2018-10-29 | 2019-04-16 | 上海无线电设备研究所 | A kind of kinematic error is to Terahertz circular track SAR image quality impact analysis method |
CN109557541A (en) * | 2018-12-17 | 2019-04-02 | 中国人民解放军国防科技大学 | Holographic penetration imaging radar polar coordinate data processing method |
CN110554387A (en) * | 2019-09-04 | 2019-12-10 | 中国科学院电子学研究所 | synthetic aperture interference near-field active source imaging method and device |
CN111220979A (en) * | 2020-01-16 | 2020-06-02 | 电子科技大学 | Imaging method for curved synthetic aperture radar |
CN113239773A (en) * | 2021-05-08 | 2021-08-10 | 重庆邮电大学 | Human body posture non-contact sensing method based on CSAR three-dimensional imaging |
CN113239773B (en) * | 2021-05-08 | 2022-07-22 | 重庆邮电大学 | Human body perception method based on CSAR three-dimensional imaging |
CN115982571A (en) * | 2023-01-12 | 2023-04-18 | 中国人民解放军军事科学院系统工程研究院 | Geometric representation method and device for space terahertz communication signals |
Also Published As
Publication number | Publication date |
---|---|
CN105044719B (en) | 2017-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105044719A (en) | Terahertz high-precision vertical curved surface imaging method based on circumference SAR | |
CN102393518B (en) | Airborne SAR (synthetic aperture radar) imaging method suitable for large squint angle | |
CN102778681B (en) | Method for imaging stationary transmitter bistatic foresight synthetic aperture radar (ST-BFSAR) | |
CN103454638B (en) | Circular synthetic aperture radar three-dimension layer tomographic imaging method | |
CN103760558B (en) | Terahertz radar ISAR imaging method | |
CN102749621B (en) | Bistatic synthetic aperture radar (BSAR) frequency domain imaging method | |
CN103869311A (en) | Real beam scanning radar super-resolution imaging method | |
CN105487074B (en) | A kind of double-base synthetic aperture radar numerical distance Doppler imaging method | |
CN102914773B (en) | Multi-pass circumference SAR three-dimensional imaging method | |
CN102004250B (en) | Frequency domain expansion based spaceborne/airborne hybrid bistatic synthetic aperture radar imaging method | |
CN103424741A (en) | Smooth procession cone parameter estimation method based on high-resolution ISAR imaging | |
CN104749570A (en) | Shift invariant airborne bistatic synthetic aperture radar target positioning method | |
CN104535978A (en) | Three-dimensional InISAR image registration and fusion method based on mutual information | |
CN103308913A (en) | Foresight SAR ambiguity resolving algorithm of double antennas carried by high-speed aircraft | |
CN103969628A (en) | Compressive-sensing-based synthetic aperture radar PFA imaging algorithm | |
CN104166134A (en) | Real beam foresight scanning radar target two-dimension locating method | |
CN105334540A (en) | Active microwave holographic security check instrument system | |
CN104166129A (en) | Real beam radar iteration minimum mean square error angle super-resolution method | |
CN103675816A (en) | Compressed sensing theory based interferometric inverse synthetic aperture radar imaging method | |
CN103630899B (en) | Method for high-resolution radar compressed sensing imaging of moving object on ground | |
CN102928841A (en) | Series inversion-based airborne circular scanning SAR (Synthetic Aperture Radar) imaging method | |
CN104407349A (en) | Frequency domain imaging method for one-stationary bistatic low-frequency ultra-wide band SAR | |
CN104515989A (en) | Three-dimensional holographic imaging method and system for close-range millimeter waves | |
CN106054152A (en) | Non-ideal locus SAR echo obtaining method based on inverse extended Omega-K algorithm | |
CN112285707B (en) | Passive multi-base high-resolution imaging method based on GPS navigation signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171017 |