CN102680956B - Energy inhibiting system for azimuth first fuzzy region echo signal of space-borne synthetic aperture radar (SAR) - Google Patents

Energy inhibiting system for azimuth first fuzzy region echo signal of space-borne synthetic aperture radar (SAR) Download PDF

Info

Publication number
CN102680956B
CN102680956B CN 201210149529 CN201210149529A CN102680956B CN 102680956 B CN102680956 B CN 102680956B CN 201210149529 CN201210149529 CN 201210149529 CN 201210149529 A CN201210149529 A CN 201210149529A CN 102680956 B CN102680956 B CN 102680956B
Authority
CN
China
Prior art keywords
orientation
distance
grid
region
confusion region
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
Application number
CN 201210149529
Other languages
Chinese (zh)
Other versions
CN102680956A (en
Inventor
于泽
刘敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN 201210149529 priority Critical patent/CN102680956B/en
Publication of CN102680956A publication Critical patent/CN102680956A/en
Application granted granted Critical
Publication of CN102680956B publication Critical patent/CN102680956B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses an energy inhibiting system for an azimuth first fuzzy region echo signal of a space-borne synthetic aperture radar (SAR). The energy inhibiting system runs in an SAR signal processor; and the raw echo of the SAR is processed in the SAR signal processor, the energy of the azimuth first fuzzy region echo signal is eliminated, and a residue signal is imaged, so that a non-fuzzy imaging result is obtained. In the system, the relevance between a first fuzzy region echo signal and raw echo data is utilized, the scattering coefficient of a first azimuth fuzzy region is resolved on the basis of a least square method, and an echo signal rebuilt by using the first fuzzy region scattering coefficient is eliminated from an acquired echo signal, so that a non-fuzzy main echo signal is obtained. As proved by simulation and imaging of a point target and a face target of a typical space-borne SAR parameter, the energy inhibiting system designed in the invention can be suitable for any space-borne SAR system.

Description

The satellite-borne synthetic aperture radar orientation suppresses system to the first confusion region echoed signal energy
Technical field
The present invention relates to a kind of disposal system of the echoed signal to synthetic-aperture radar, more particularly say, refer to that a kind of azimuth ambiguity for the satellite-borne synthetic aperture radar imaging processing suppresses system.
Background technology
Synthetic-aperture radar (Synthetic Aperture Radar, SAR) system is a kind of active remote sensing equipment, its is initiatively to target emission electromagnetic wave, utilizes to receive echoed signal under the target scattering characteristics effect and carry out imaging.Satellite-borne synthetic aperture radar is a very complicated system, and systematic parameter is numerous and interrelated, mutual restriction.To, radar emission linear FM signal, through echoed signal is carried out pulse compression, thereby obtain higher distance to resolution in distance; In the orientation to, Texas tower via satellite in the variation of process wave beam and target relative position of advancing, form the equivalent linear array antenna, thereby improved azimuthal resolution.
One of the key factor that influences the picture quality of satellite-borne SAR is exactly the azimuth ambiguity problem.This pulsed mode with satellite-borne SAR is relevant, and (Pulse Repetitive Frequency, PRF) selection and azimuth ambiguity are closely related its pulse repetition rate.Because Doppler effect, the echo that SAR receives along the orientation to can equivalence being that (sampling frequency is PRF to a linear frequency modulation for linear frequency modulation, LFM) signal, and this LFM signal disperses.The orientation to the LFM signal by the orientation to the antenna radiation pattern weighting, but antenna radiation pattern all exists secondary lobe and is on the higher frequency band.
" radar imagery technology " that in October, 1999 publishing house of Harbin Institute of Technology publication, Liu Yongtan write pointed out, the Doppler frequency spectrum of synthetic-aperture radar is non-band limit, the target echo spectrum repeats with pulse repetition rate, therefore the echoed signal outside the primary spectrum will be folded into the primary spectrum district, thereby cause azimuth ambiguity.Synthetic-aperture radar (SAR, Synthetic Aperture Radar) is installed on the motion platform, and repetition frequency emission, received pulse according to certain form echoed signal.The structured flowchart of SAR system as shown in Figure 1, the SAR system includes radar system on the star, satellite platform and data down transmission system and ground system three parts, and the synthetic aperture Radar Imaging Processing is finished in ground system.The echoed signal that ground system issues by ground receiving station reception satellite platform and data down transmission system, this echoed signal is carried out imaging processing through the SAR signal processor, obtains the SAR image; Described SAR image is stored in and is equipped with in the shelves operating system.Provided typical satellite-borne SAR azimuth ambiguity imaging results in Gabban A., Greidanus H., Smith A.J.E., Anitori L., Thoorens F.-X., Malorqui J. deliver on the 3rd TerraSAR-X Science Team Meeting in 2008 " Ship surveillance with terrasar-x scansar ".In theory, the azimuth ambiguity district of satellite-borne SAR has countless a plurality of, " the On suppressing azimuth ambiguities of Synthetic Aperture Radar by three filters " that Zhimin Z., Zhensong W. deliver on the calendar year 2001 CIE International Conference on Radar point out blurred signal concentration of energy more than 85% in the orientation to first confusion region.Wherein, the orientation comprises the orientation to+1 confusion region and two zones ,-1 confusion region to first confusion region.
Suppress about azimuth ambiguity, mainly contain two kinds of methods at present.1993 the 31st volume " IEEE Transactions on Geoscience and Remote Sensing " Moreira A. " the Suppressing the Azimuth Ambiguities in Synthetic Aperture Radar Images " that deliver the method that homophase offsets has been proposed, it must be point target that this method requires the source of generation azimuth ambiguity, and will accurately estimate its position.2005 the 43rd volume " IEEE Transactions on Geoscience and Remote Sensing " Guarnieri A.M. " the Adaptive removal of azimuth ambiguities in SAR images " that deliver adopt adaptive wiener filter from imaging results, to estimate not have fuzzy imaging results, this method requires observation scene to obey Gaussian process again.
Summary of the invention
The characteristics that the limitation that exists at said method and azimuth ambiguity energy are mainly derived from first confusion region, the present invention handles in the SAR signal processor the SAR original echo, at first reject the orientation to the energy of the first confusion region echoed signal, then residue signal is carried out imaging processing, obtain not having fuzzy imaging results.The present invention proposes a kind of to the satellite-borne synthetic aperture radar orientation to system that the energy of the first confusion region echoed signal suppresses to handle, this system has utilized the correlativity of the first confusion region echoed signal and original echo data, parse the scattering coefficient of first orientation confusion region based on least square method, from the echoed signal of obtaining, reject the echoed signal of utilizing the first confusion region scattering coefficient to rebuild then, thereby obtain not having fuzzy primary area echoed signal.By point target and appearance mark emulation and the imaging to typical satellite-borne SAR parameter, proved that the energy inhibition system of the present invention's design can be applicable to Spaceborne SAR System arbitrarily.
A kind of satellite-borne synthetic aperture radar of the present invention orientation suppresses system to the first confusion region echoed signal energy, this energy inhibition system operates in the SAR signal processor, and described SAR signal processor includes echoed signal and arranges module (1), main imaging area space orientation module (2), the orientation is to+1 confusion region space orientation module (3), the orientation is to+1 confusion region spatial discretization module (4), the orientation is to+1 confusion region inverting and eliminate obfuscation module (5), the orientation is to-1 confusion region space orientation module (6), the orientation is to-1 confusion region spatial discretization module (7), the orientation is to-1 confusion region inverting and eliminate obfuscation module (8), main imaging area spatial discretization module (9) and orientation are to main imaging area processing module (10);
(A) arrange the spaceborne echo data E to receiving in the module (1) in echoed signal The two dimension echoBe arranged as one-dimensional data E The one dimension echo
(B) in main imaging area space orientation module (2) according to the systematic parameter of satellite-borne synthetic aperture radar, determine the locus of the main imaging region of echo data correspondence;
(C) in the orientation in+1 confusion region space orientation module (3) according to the position of main imaging region, determine the orientation to the position of+1 confusion region;
In+1 confusion region spatial discretization module (4), spatial discretization is carried out to+1 confusion region in the orientation in the orientation, make up and obtain the orientation to the measurement matrix Φ of the first confusion region correspondence + 1
In the orientation to+1 confusion region inverting with eliminate in the obfuscation module (5) according to least square method from measuring matrix Φ + 1With one dimension echo data E The one dimension echoIn obtain the orientation to the scattering coefficient of+1 confusion region
Figure GDA00003427071100021
Then from one dimension echo data E The one dimension echoIn deduct the orientation and obtain residue signal E to the echoed signal of+1 confusion region correspondence Residue signal 1
(D) in the orientation in-1 confusion region space orientation module (6) according to the position of main imaging region, determine the orientation to the position of-1 confusion region;
In-1 confusion region spatial discretization module (7), spatial discretization is carried out to-1 confusion region in the orientation in the orientation, make up and obtain the orientation to the measurement matrix Φ of the first confusion region correspondence -1
In the orientation to-1 confusion region inverting with eliminate in the obfuscation module (8) according to least square method from measuring matrix Φ -1With residue signal E Residue signal 1In obtain the orientation to the scattering coefficient of-1 confusion region
Figure GDA00003427071100022
, then from residue signal E Residue signal 1In deduct the orientation and obtain residue signal E to the echoed signal of-1 confusion region correspondence Residue signal 2
(E) in main imaging area spatial discretization module (9), main imaging region is carried out spatial discretization, make up the measurement matrix φ of main imaging region correspondence;
(F) in the orientation in the main imaging area processing module (10) according to measurement matrix φ and the residue signal E of least square method from main imaging region correspondence Residue signal 2Obtain final echo-wave imaging I as a result Abs
In the satellite-borne synthetic aperture radar azimuth ambiguity suppresses to handle, use method of the present invention and have following advantage:
1. echoed signal is rearranged for one dimension by two dimension, make up the observing matrix of+1 azimuth ambiguity district ,-1 azimuth ambiguity district, main imaging region correspondence respectively, the echoed signal of from echo data, separating+1 azimuth ambiguity district and-1 azimuth ambiguity district successively, then residue signal is handled, suppressed the azimuth ambiguity energy in the final process result significantly.
2. compare with the azimuth ambiguity inhibition method that homophase offsets, the present invention does not require that the azimuth ambiguity source is necessary for point target, the target that can have any shape, and do not need its positional information accurately.
3. compare with the azimuth ambiguity inhibition technology that adopts adaptive wiener filter, the present invention does not have any specific (special) requirements to the statistical property of scene, is applicable to various types of imaging regions.
Description of drawings
Fig. 1 is the SAR system chart.
Fig. 2 is the structured flowchart of each module in the SAR signal processor of the present invention.
Fig. 3 A is that the orientation is to the regional schematic diagram of+1 confusion region.
Fig. 3 B is the regional schematic diagram of main imaging area.
Fig. 3 C is that the orientation is to the regional schematic diagram of-1 confusion region.
Fig. 3 D is the azimuth ambiguity synoptic diagram that has the angle of squint.
Fig. 4 is the imaging results 3-D view that has the azimuth ambiguity energy;
Fig. 5 is the imaging results 3-D view after employing the inventive method.
Embodiment
The present invention is described in further detail below in conjunction with accompanying drawing and simulation example.
Referring to shown in Figure 2, a kind of satellite-borne synthetic aperture radar orientation that the present invention proposes suppresses system to the first confusion region echoed signal energy, this energy inhibition system operates in the SAR signal processor, and described SAR signal processor includes echoed signal and arranges module 1, main imaging area space orientation module 2, the orientation is to+1 confusion region space orientation module 3, the orientation is to+1 confusion region spatial discretization module 4, the orientation is to+1 confusion region inverting and eliminate obfuscation module 5, the orientation is to-1 confusion region space orientation module 6, the orientation is to-1 confusion region spatial discretization module 7, the orientation is to-1 confusion region inverting and eliminate obfuscation module 8, main imaging area spatial discretization module 9 and orientation are to main imaging area processing module 10.The present invention is rearranged for one dimension with echoed signal by two dimension, make up the observing matrix of+1 azimuth ambiguity district ,-1 azimuth ambiguity district, main imaging region correspondence respectively, the echoed signal of from echo data, separating+1 azimuth ambiguity district and-1 azimuth ambiguity district successively, then residue signal is handled, suppressed the azimuth ambiguity energy in the final process result (echo-wave imaging) significantly.
The function that each module realizes in the energy inhibition system of the present invention's design is:
(1) echoed signal is arranged module 1
In the present invention, echoed signal is arranged the spaceborne echo data E that 1 pair of module receives The two dimension echoBe arranged as the spaceborne data E of one dimension The one dimension echo
In the present invention, with the spaceborne echo data E of synthetic-aperture radar The two dimension echoThe employing matrix form is expressed as
Figure GDA00003427071100031
N The orientationExpression along the orientation to sampling number; N DistanceThe expression along the distance to sampling number; a 1-1The 1st orientation that expression upwards collects along orientation distance constantly to the 1st sampled point; a 1-2The 1st orientation that expression upwards collects along orientation distance constantly to the 2nd sampled point;
Figure GDA00003427071100032
The 1st orientation that expression upwards collects along orientation distance constantly to N DistanceIndividual sampled point; a 2-1The 2nd orientation that expression upwards collects along orientation distance constantly to the 1st sampled point; a 2-2The 2nd orientation that expression upwards collects along orientation distance constantly to the 2nd sampled point; The 2nd orientation that expression upwards collects along orientation distance constantly to N DistanceIndividual sampled point;
Figure GDA00003427071100034
The N that expression upwards collects along the orientation The orientationIndividual orientation distance constantly to the 1st sampled point;
Figure GDA00003427071100045
The N that expression upwards collects along the orientation The orientationIndividual orientation distance constantly to the 2nd sampled point;
Figure GDA00003427071100046
The N that expression upwards collects along the orientation The orientationIndividual orientation distance constantly to N DistanceIndividual sampled point.
In the present invention, will Be arranged as one dimension echo data E The one dimension echoMatrix form is expressed as:
Figure GDA00003427071100042
In the present invention, with spaceborne echo data E The two dimension echoBe arranged as the spaceborne data E of one dimension The one dimension echoBe conducive to import the calculating of carrying out scattering coefficient in the least square method.
(2) main imaging area space orientation module 2
According to the spaceborne synthetic aperture radar (SAR) system parameter, calculate spaceborne echo data E The two dimension echoThe locus of corresponding main imaging region (shown in Fig. 3 B), the central point of main imaging region is designated as O Main, adopt coordinate to be expressed as on four summits of main imaging region:
The upper left corner coordinate of main imaging region is designated as AA (P A left side, P On),
The upper right corner coordinate of main imaging region is designated as BB (P Right, P On),
The lower left corner coordinate of main imaging region is designated as CC (P A left side, P Down),
The lower right corner coordinate of main imaging region is designated as DD (P Right, P Down).Wherein:
P OnThe peak coordinate of representing the locus of main imaging region, and
P DownThe minimum point coordinate of representing the locus of main imaging region, and
Figure GDA00003427071100044
P A left sideThe leftmost side point coordinate of representing the locus of main imaging region, and
P RightThe rightmost side point coordinate of representing the locus of main imaging region, and
Figure GDA00003427071100052
λ represents the radar system wavelength; L aThe expression orientation is to antenna size; R MinThe biography antenna is to the bee-line of main imaging region under the expression; R MaxThe biography antenna is to the longest distance of main imaging region under the expression; Velocity equivalent when V represents satellite platform flight; PRT indicating impulse repetition interval; τ pRepresent exomonental time width; F sThe expression distance is to the impulse sampling rate; C represents the light velocity.
In the present invention, utilize coordinate to represent the zone, the institute in this zone a bit can satisfy the irradiation of satellite-borne synthetic aperture radar full aperture.
(3) orientation is to+1 confusion region space orientation module 3
According to the position (shown in Fig. 3 B) of main imaging region, determine the orientation to+1 confusion region (shown in Fig. 3 A) corresponding main imaging region central point O MainThe position, the orientation is designated as O to the central point of+1 confusion region Side+1, adopt coordinate to be expressed as to four summits of+1 confusion region in the orientation:
The orientation is designated as AA to the upper left corner coordinate of+1 confusion region + 1(P + 1 left side, P On+1),
The orientation is designated as BB to the upper right corner coordinate of+1 confusion region + 1(P + 1 right side, P On+1),
The orientation is designated as CC to the lower left corner coordinate of+1 confusion region + 1(P + 1 left side, P + 1 time),
The orientation is designated as DD to the lower right corner coordinate of+1 confusion region + 1(P + 1 right side, P + 1 time).
Wherein:
Figure GDA00003427071100053
Figure GDA00003427071100054
Figure GDA00003427071100055
Figure GDA00003427071100056
The skew (shown in Fig. 3 D) that the expression orientation makes progress in the orientation to the main imaging region center of+1 confusion region correspondence, and
Figure GDA00003427071100059
The skew (shown in Fig. 3 D) that the expression orientation makes progress in distance to the main imaging region center of+1 confusion region correspondence, and
θ +Represent the orientation to the angle of squint of the main imaging region center correspondence of+1 confusion region correspondence, and
θ + = arcsin [ sin θ + λ × f p 2 × V ] ;
R MinThe biography antenna is to the bee-line at main imaging region center under the expression; f pThe indicating impulse repetition frequency; θ represents the angle of squint at main imaging region center.
(4) orientation is to+1 confusion region spatial discretization module 4
In the present invention, (shown in Fig. 3 A) carries out spatial discretization to+1 confusion region to the orientation to adopt mesh generation.
Step 401: adopt mesh generation to carry out the orientation and divide processing to the territory of+1 confusion region, described orientation is divided into to+1 confusion region
Figure GDA00003427071100061
Individual grid;
Step 402: the orientation to the orientation of+1 confusion region to net point
Figure GDA00003427071100062
With the distance to net point
Figure GDA00003427071100063
The grid element center coordinate be designated as
Figure GDA00003427071100064
Figure GDA00003427071100065
The expression orientation to the grid number, and The grid number rounds up;
Figure GDA00003427071100067
Expression distance to the grid number, and
Figure GDA00003427071100068
The grid number rounds up;
Step 403: at any time t of satellite platform flight kDown to the grid element center coordinate
Figure GDA00003427071100069
Carry out oblique distance and calculate, this oblique distance is
Figure GDA000034270711000610
Step 404: at t kIn time, inscribe, the grid element center coordinate Corresponding observation vector is
Figure GDA000034270711000612
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure GDA000034270711000614
Distance to the 1st sampled point observed reading;
Figure GDA000034270711000615
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure GDA000034270711000616
Distance to the 2nd sampled point observed reading;
Figure GDA000034270711000617
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure GDA000034270711000618
Distance to n sampled point observed reading;
Described in the orientation n sampled point observed reading in+1 confusion region
The b frequency modulation rate of representing to transmit; I represents imaginary unit; N represents the grid element center coordinate
Figure GDA000034270711000620
Distance to sampled point; Represent distance to starting point, and
Figure GDA000034270711000622
Round up;
Figure GDA000034270711000623
Represent distance to terminal point, and
Figure GDA000034270711000624
Round downwards;
Figure GDA000034270711000625
Be illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to N DistanceIndividual sampled point observed reading;
Step 405: be structured in total movement time interior orientation to
Figure GDA000034270711000626
Individual, distance is to the
Figure GDA000034270711000627
The observation vector of individual grid element center correspondence, as follows:
Wherein,
Figure GDA00003427071100072
Expression
Figure GDA00003427071100073
Transposition;
Step 406: make up the orientation and measure matrix to+1 confusion region, as follows:
Figure GDA00003427071100074
Figure GDA00003427071100075
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region the observation vector to the 1st grid, distance to the 1st grid element center point;
Figure GDA00003427071100076
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region the observation vector to the 1st grid, distance to the 2nd grid element center point;
Figure GDA00003427071100077
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to the 1st grid, distance to the
Figure GDA00003427071100078
The observation vector of individual grid element center point;
Figure GDA00003427071100079
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region the observation vector to the 2nd grid, distance to the 1st grid element center point;
Figure GDA000034270711000710
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region the observation vector to the 2nd grid, distance to the 2nd grid element center point;
Figure GDA000034270711000711
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to the 2nd grid, distance to the
Figure GDA000034270711000712
The observation vector of individual grid element center point; Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to
Figure GDA000034270711000714
Individual grid, distance are to the observation vector of the 1st grid element center point;
Figure GDA000034270711000715
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to
Figure GDA000034270711000716
Individual grid, distance are to the observation vector of the 2nd grid element center point;
Figure GDA000034270711000717
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to
Figure GDA000034270711000718
Individual grid, distance are to the
Figure GDA000034270711000719
The observation vector of individual grid element center point.
(5) orientation is to+1 confusion region inverting and elimination obfuscation module 5
According to least square method
Figure GDA000034270711000720
Calculate the orientation to the scattering coefficient of+1 confusion region
Figure GDA000034270711000721
As follows:
Figure GDA00003427071100081
α 1-1The expression orientation is the scattering coefficient to the 1st grid, distance to the 1st grid element center point to the orientation of+1 confusion region; α 1-2The expression orientation is the scattering coefficient to the 1st grid, distance to the 2nd grid element center point to the orientation of+1 confusion region
Figure GDA00003427071100082
The expression orientation to the orientation of+1 confusion region to the 1st grid, distance to the
Figure GDA00003427071100083
The scattering coefficient of individual grid element center point; α 2-1The expression orientation is the scattering coefficient to the 2nd grid, distance to the 1st grid element center point to the orientation of+1 confusion region; α 2-2The expression orientation is the scattering coefficient to the 2nd grid, distance to the 2nd grid element center point to the orientation of+1 confusion region;
Figure GDA00003427071100084
The expression orientation to the orientation of+1 confusion region to the 2nd grid, distance to the
Figure GDA00003427071100085
The scattering coefficient of individual grid element center point;
Figure GDA00003427071100086
The expression orientation to the orientation of+1 confusion region to the Individual grid, distance are to the scattering coefficient of the 1st grid element center point;
Figure GDA00003427071100088
The expression orientation to the orientation of+1 confusion region to the
Figure GDA00003427071100089
Individual grid, distance are to the scattering coefficient of the 2nd grid element center point;
Figure GDA000034270711000810
The expression orientation to the orientation of+1 confusion region to the Individual grid, distance are to the
Figure GDA000034270711000812
The scattering coefficient of individual grid element center point; From echo data E The one dimension echoThe middle orientation of separating obtains residue signal E to the echoed signal of+1 confusion region correspondence Residue signal 1For
Figure GDA000034270711000813
(6) orientation is to-1 confusion region space orientation module 6
According to the position (shown in Fig. 3 B) of main imaging region, determine the orientation to-1 confusion region (shown in Fig. 3 C) corresponding main imaging region central point O MainThe position, the orientation is designated as O to the central point of-1 confusion region Side-1, adopt coordinate to be expressed as to four summits of-1 confusion region in the orientation:
The orientation is designated as AA to the upper left corner coordinate of-1 confusion region -1(P -1 left side, P On-1),
The orientation is designated as BB to the upper right corner coordinate of-1 confusion region -1(P -1 right side, P On-1),
The orientation is designated as CC to the lower left corner coordinate of-1 confusion region -1(P -1 left side, P -1 time),
The orientation is designated as DD to the lower right corner coordinate of-1 confusion region -1(P -1 right side, P -1 time).
Wherein,
Figure GDA000034270711000814
Figure GDA000034270711000816
Figure GDA000034270711000817
Figure GDA00003427071100091
The skew (shown in Fig. 3 D) that the expression orientation makes progress in the orientation to the main imaging region center of-1 confusion region correspondence, and
Figure GDA00003427071100092
The skew (shown in Fig. 3 D) that the expression orientation makes progress in distance to the main imaging region center of-1 confusion region correspondence, and
Figure GDA00003427071100094
θ -Represent the orientation to the angle of squint of the main imaging region center correspondence of-1 confusion region correspondence, and
θ - = arcsin [ sin θ - λ × f p 2 × V ] ;
(7) orientation is to-1 confusion region spatial discretization module 7
In the present invention, (shown in Fig. 3 C) carries out spatial discretization to-1 confusion region to the orientation to adopt mesh generation.
Step 701: adopt mesh generation to carry out the orientation and divide processing to the territory of-1 confusion region, described orientation is divided into to-1 confusion region
Figure GDA00003427071100096
Individual grid;
Step 702: the orientation to the orientation of-1 confusion region to net point
Figure GDA00003427071100097
With the distance to net point
Figure GDA00003427071100098
The grid element center coordinate be designated as
Figure GDA00003427071100099
Figure GDA000034270711000910
The expression orientation to the grid number, and
Figure GDA000034270711000911
The grid number rounds up;
Figure GDA000034270711000912
Expression distance to the grid number, and The grid number rounds up;
Step 703: at any time t of satellite platform flight kDown to the grid element center coordinate
Figure GDA000034270711000914
Carry out oblique distance and calculate, this oblique distance is
Step 704: at t kIn time, inscribe, the grid element center coordinate
Figure GDA000034270711000916
Corresponding observation vector is
Figure GDA000034270711000917
Figure GDA000034270711000918
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure GDA000034270711000919
Distance to the 1st sampled point observed reading;
Figure GDA000034270711000920
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure GDA000034270711000921
Distance to the 2nd sampled point observed reading;
Be illustrated in t kIn time, inscribe, the grid element center coordinate Distance to n sampled point observed reading;
Described in the orientation n sampled point observed reading in-1 confusion region
Figure GDA00003427071100101
N represents the grid element center coordinate
Figure GDA00003427071100102
Distance to sampled point;
Figure GDA00003427071100103
Represent distance to starting point, and
Figure GDA00003427071100104
Round up;
Figure GDA00003427071100105
Represent distance to terminal point, and
Figure GDA00003427071100106
Round downwards; Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure GDA00003427071100108
Distance to N DistanceIndividual sampled point observed reading;
Step 705: be structured in total movement time interior orientation to
Figure GDA00003427071100109
Individual, distance is to the
Figure GDA000034270711001010
The observation vector of individual grid element center correspondence, as follows:
Wherein,
Figure GDA000034270711001012
Expression
Figure GDA000034270711001013
Transposition;
Step 706: make up the orientation and measure matrix to-1 confusion region, as follows:
Figure GDA000034270711001014
Figure GDA000034270711001015
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region the observation vector to the 1st grid, distance to the 1st grid element center point;
Figure GDA000034270711001016
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region the observation vector to the 1st grid, distance to the 2nd grid element center point;
Figure GDA000034270711001017
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region to the 1st grid, distance to the
Figure GDA000034270711001018
The observation vector of individual grid element center point; Be illustrated in total movement time interior orientation to the orientation of-1 confusion region the observation vector to the 2nd grid, distance to the 1st grid element center point;
Figure GDA000034270711001020
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region the observation vector to the 2nd grid, distance to the 2nd grid element center point;
Figure GDA000034270711001021
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region to the 2nd grid, distance to the
Figure GDA000034270711001022
The observation vector of individual grid element center point;
Figure GDA000034270711001023
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region to
Figure GDA000034270711001024
Individual grid, distance are to the observation vector of the 1st grid element center point;
Figure GDA000034270711001025
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to
Figure GDA000034270711001026
Individual grid, distance are to the observation vector of the 2nd grid element center point;
Figure GDA000034270711001027
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region to
Figure GDA000034270711001028
Individual grid, distance are to the
Figure GDA00003427071100111
The observation vector of individual grid element center point.
(8) orientation is to-1 confusion region inverting and elimination obfuscation module 8
According to least square method
Figure GDA00003427071100112
Calculate the orientation to the scattering coefficient of-1 confusion region
Figure GDA00003427071100113
As follows:
Figure GDA00003427071100114
β 1-1The expression orientation is the scattering coefficient to the 1st grid, distance to the 1st grid element center point to the orientation of-1 confusion region; β 1-2The expression orientation is the scattering coefficient to the 1st grid, distance to the 2nd grid element center point to the orientation of-1 confusion region;
Figure GDA00003427071100115
The expression orientation to the orientation of-1 confusion region to the 1st grid, distance to the
Figure GDA00003427071100116
The scattering coefficient of individual grid element center point; β 2-1The expression orientation is the scattering coefficient to the 2nd grid, distance to the 1st grid element center point to the orientation of-1 confusion region; β 2-2The expression orientation is the scattering coefficient to the 2nd grid, distance to the 2nd grid element center point to the orientation of-1 confusion region;
Figure GDA00003427071100117
The expression orientation to the orientation of-1 confusion region to the 2nd grid, distance to the
Figure GDA00003427071100118
The scattering coefficient of individual grid element center point;
Figure GDA00003427071100119
The expression orientation to the orientation of-1 confusion region to the
Figure GDA000034270711001110
Individual grid, distance are to the scattering coefficient of the 1st grid element center point;
Figure GDA000034270711001111
The expression orientation to the orientation of-1 confusion region to the
Figure GDA000034270711001112
Individual grid, distance are to the scattering coefficient of the 2nd grid element center point;
Figure GDA000034270711001113
The expression orientation to the orientation of-1 confusion region to the
Figure GDA000034270711001114
Individual grid, distance are to the
Figure GDA000034270711001115
The scattering coefficient of individual grid element center point.
From residue signal E Residue signal 1The middle orientation of separating obtains residue signal E to the echoed signal of-1 confusion region correspondence Residue signal 2For
Figure GDA000034270711001116
In the present invention, the orientation all utilizes coordinate to represent the zone to+1 confusion region and orientation to-1 confusion region, and this zone is being carried out discrete, the mesh generation of regional space, thereby structure obtains measuring matrix.Last echoed signal according to least square method, separation confusion region correspondence, and then obtain residue signal.
(9) main imaging area spatial discretization module 9
In the present invention, adopt mesh generation that main imaging region (shown in Fig. 3 B) is carried out spatial discretization.
Step 901: the territory of adopting mesh generation to carry out main imaging area is divided and is handled, and described main imaging region is divided into D The orientation* D DistanceIndividual grid;
Step 902: the orientation of main imaging region is to net point I aWith the distance to net point I rThe grid element center coordinate be designated as
Figure GDA00003427071100121
D The orientationThe expression orientation to the grid number, and
Figure GDA00003427071100122
The grid number rounds up; D DistanceExpression distance to the grid number, and
Figure GDA00003427071100123
The grid number rounds up;
Step 903: at any time t of satellite platform flight kDown to grid element center coordinate W GridCarry out oblique distance and calculate, this oblique distance is
Figure GDA00003427071100124
Step 904: at t kIn time, inscribe, grid element center coordinate W GridCorresponding observation vector is
Figure GDA00003427071100125
e K-1Be illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to the 1st sampled point observed reading; e K-2Be illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to the 2nd sampled point observed reading; e K-nBe illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to n sampled point observed reading;
Described n sampled point observed reading in main imaging area
Figure GDA00003427071100126
N represents grid element center coordinate W GridDistance to sampled point; N BeginRepresent distance to starting point, and
Figure GDA00003427071100127
Round up; N EndRepresent distance to terminal point, and
Figure GDA00003427071100128
Round downwards;
Figure GDA00003427071100129
Be illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to N DistanceIndividual sampled point observed reading;
Step 905: be structured in total movement time interior orientation to I aIndividual, the distance to I rThe observation vector of individual grid element center correspondence, as follows:
Wherein, Expression
Figure GDA00003427071100133
Transposition;
Step 906: make up main imaging region and measure matrix, as follows:
Figure GDA00003427071100134
(10) orientation is to main imaging area processing module 10
According to least square method E The scattering of one dimension primary area=(φ H* φ) -1* φ H* E Residue signal 2Calculate main imaging region scattering coefficient one-dimensional representation E The scattering of one dimension primary area, it is as follows namely to obtain final imaging results:
Figure GDA00003427071100135
σ 1-1Represent that the orientation of main imaging area is to the 1st grid, apart from the scattering coefficient to the 1st grid element center point; σ 1-2Represent that the orientation of main imaging area is to the 1st grid, apart from the scattering coefficient to the 2nd grid element center point;
Figure GDA00003427071100136
The orientation of representing main imaging area to the 1st grid, distance to D DistanceThe scattering coefficient of individual grid element center point; σ 2-1Represent that the orientation of main imaging area is to the 2nd grid, apart from the scattering coefficient to the 1st grid element center point; σ 2-2Represent that the orientation of main imaging area is to the 2nd grid, apart from the scattering coefficient to the 2nd grid element center point;
Figure GDA00003427071100137
The orientation of representing main imaging area to the 2nd grid, distance to D DistanceThe scattering coefficient of individual grid element center point; Represent that the orientation of main imaging area is to D The orientationIndividual grid, distance are to the scattering coefficient of the 1st grid element center point;
Figure GDA00003427071100139
Represent that the orientation of main imaging area is to D The orientationIndividual grid, distance are to the scattering coefficient of the 2nd grid element center point;
Figure GDA00003427071100141
Represent that the orientation of main imaging area is to D The orientationIndividual grid, distance are to D DistanceThe scattering coefficient of individual grid element center point.
With E The scattering of one dimension primary areaBe converted to two-dimensional matrix E by one-dimensional vector Two dimension primary area scattering, can be expressed as:
Figure GDA00003427071100142
To matrix E Two dimension primary area scatteringAsk absolute value, obtain the imaging results of imaging region, be expressed as:
Figure GDA00003427071100143
In the present invention, utilize least square method to carry out imaging processing, accurately the scattering coefficient of target in the estimated image.
To describe each module below in detail and be the processing procedure of how spaceborne echoed signal being carried out:
(A) arrange the spaceborne echo data E to receiving in the module 1 in echoed signal The two dimension echoBe arranged as one-dimensional data E The one dimension echo
(B) in main imaging area space orientation module 2 according to the systematic parameter of satellite-borne synthetic aperture radar, determine the locus of the main imaging region of echo data correspondence;
(C) in the orientation in+1 confusion region space orientation module 3 according to the position of main imaging region, determine the orientation to the position of+1 confusion region;
In+1 confusion region spatial discretization module 4, spatial discretization is carried out to+1 confusion region in the orientation in the orientation, make up and obtain the orientation to the measurement matrix Φ of the first confusion region correspondence + 1
In the orientation to+1 confusion region inverting with eliminate in the obfuscation module 5 according to least square method from measuring matrix Φ + 1With one dimension echo data E The one dimension echoIn obtain the orientation to the scattering coefficient of+1 confusion region
Figure GDA00003427071100144
Then from one dimension echo data E The one dimension echoIn deduct the orientation and obtain residue signal E to the echoed signal of+1 confusion region correspondence Residue signal 1
(D) in the orientation in-1 confusion region space orientation module 6 according to the position of main imaging region, determine the orientation to the position of-1 confusion region;
In-1 confusion region spatial discretization module 7, spatial discretization is carried out to-1 confusion region in the orientation in the orientation, make up and obtain the orientation to the measurement matrix Φ of the first confusion region correspondence -1
In the orientation to-1 confusion region inverting with eliminate in the obfuscation module 8 according to least square method from measuring matrix Φ -1With residue signal E Residue signal 1In obtain the orientation to the scattering coefficient of-1 confusion region
Figure GDA00003427071100145
Then from residue signal E Residue signal 1In deduct the orientation and obtain residue signal E to the echoed signal of-1 confusion region correspondence Residue signal 2
(E) in main imaging area spatial discretization module 9, main imaging region is carried out spatial discretization, make up the measurement matrix φ of main imaging region correspondence;
(F) in the orientation in the main imaging area processing module 10 according to measurement matrix φ and the residue signal E of least square method from main imaging region correspondence Residue signal 2Obtain final imaging results I Abs
The simulating, verifying example
Adopt the satellite-borne synthetic aperture radar orientation of the present invention's design to the first confusion region echoed signal energy inhibition system spaceborne echo data to be handled.
Spaceborne synthetic aperture radar (SAR) system parameter such as table 1:
The setting of table 1 imaging parameters
The light velocity 3×10 8m/s
Wavelength 0.03m
Orbit altitude 500km
The angle of squint 0°
Radar is to the bee-line of main imaging region 610452.540074783m
Radar is to the longest distance of main imaging region 610325.294381233m
Radar is to the bee-line distance at main imaging region center 610387.294380728m
Velocity equivalent 7500m/s
Signal bandwidth 90MHz
Sampling rate 99MHz
Pulse width 1μs
Pulse repetition rate 3000Hz
Pulse-recurrence time 1/3000s
The orientation is to antenna length 4.8m
The orientation is to sampled point 1616
Distance is to sampled point 167
The imaging results of synthetic-aperture radar has been reacted the scattering situation that is imaged the zone.Fig. 4 is the result that adopts BP formation method (seeing " the BP imaging algorithm of ultra broadband LFM signal " article that Zhu Guofu calendar year 2001 delivered in the 17 volume fifth phase of information processing).Among this result, the target scattering information of main imaging region is submerged in the orientation fully in the signal of first confusion region, can't judge whether there is target in the main imaging region.Fig. 5 is the result who uses after the inventive method is handled same original echo data.This result shows that azimuth ambiguity inhibition of the present invention is very obvious, can't see significantly fuzzy energy and has existed, and 3 * 3 dot matrix target has obtained recovery in the main imaging region.

Claims (5)

1. a satellite-borne synthetic aperture radar orientation suppresses system to the first confusion region echoed signal energy, this energy inhibition system operates in the SAR signal processor, it is characterized in that: described SAR signal processor includes echoed signal and arranges module (1), main imaging area space orientation module (2), the orientation is to+1 confusion region space orientation module (3), the orientation is to+1 confusion region spatial discretization module (4), the orientation is to+1 confusion region inverting and eliminate obfuscation module (5), the orientation is to-1 confusion region space orientation module (6), the orientation is to-1 confusion region spatial discretization module (7), the orientation is to-1 confusion region inverting and eliminate obfuscation module (8), main imaging area spatial discretization module (9) and orientation are to main imaging area processing module (10);
(A) arrange the spaceborne echo data E to receiving in the module (1) in echoed signal The two dimension echoBe arranged as one-dimensional data E The one dimension echo
Spaceborne echo data with synthetic-aperture radar
Figure FDA00003427071000011
Be arranged as one-dimensional data
N The orientationExpression along the orientation to sampling number; N DistanceThe expression along the distance to sampling number; a 1-1The 1st orientation that expression upwards collects along orientation distance constantly to the 1st sampled point; a 1-2The 1st orientation that expression upwards collects along orientation distance constantly to the 2nd sampled point;
Figure FDA00003427071000013
The 1st orientation that expression upwards collects along orientation distance constantly to N DistanceIndividual sampled point; a 2-1The 2nd orientation that expression upwards collects along orientation distance constantly to the 1st sampled point; a 2-2The 2nd orientation that expression upwards collects along orientation distance constantly to the 2nd sampled point;
Figure FDA00003427071000014
The 2nd orientation that expression upwards collects along orientation distance constantly to N DistanceIndividual sampled point;
Figure FDA00003427071000015
The N that expression upwards collects along the orientation The orientationIndividual orientation distance constantly to the 1st sampled point;
Figure FDA00003427071000016
The N that expression upwards collects along the orientation The orientationIndividual orientation distance constantly to the 2nd sampled point; The N that expression upwards collects along the orientation The orientationIndividual orientation distance constantly to N DistanceIndividual sampled point;
(B) in main imaging area space orientation module (2) according to the systematic parameter of satellite-borne synthetic aperture radar, determine the locus of the main imaging region of echo data correspondence;
Described spaceborne echo data E The two dimension echoThe locus of corresponding main imaging region, adopt coordinate to be expressed as: the upper left corner coordinate of main imaging region is designated as AA (P A left side, P On), the upper right corner coordinate of main imaging region is designated as BB (P Right, P On), the lower left corner coordinate of main imaging region is designated as CC (P A left side, P Down), the lower right corner coordinate of main imaging region is designated as DD (P Right, P Down); Wherein: P OnThe peak coordinate of representing the locus of main imaging region, and
Figure FDA00003427071000022
P DownThe minimum point coordinate of representing the locus of main imaging region, and P A left sideThe leftmost side point coordinate of representing the locus of main imaging region, and
Figure FDA00003427071000024
P RightThe rightmost side point coordinate of representing the locus of main imaging region, and
Figure FDA00003427071000025
λ represents the radar system wavelength; L aThe expression orientation is to antenna size; R MinThe biography antenna is to the bee-line of main imaging region under the expression; R MaxThe biography antenna is to the longest distance of main imaging region under the expression; Velocity equivalent when V represents satellite platform flight; PRT indicating impulse repetition interval; τ pRepresent exomonental time width; F sThe expression distance is to the impulse sampling rate; C represents the light velocity;
(C) in the orientation in+1 confusion region space orientation module (3) according to the position of main imaging region, determine the orientation to the position of+1 confusion region;
In+1 confusion region spatial discretization module (4), spatial discretization is carried out to+1 confusion region in the orientation in the orientation, make up and obtain the orientation to the measurement matrix Φ of the first confusion region correspondence + 1
In the orientation to+1 confusion region inverting with eliminate in the obfuscation module (5) according to least square method from measuring matrix Φ + 1With one dimension echo data E The one dimension echoIn obtain the orientation to the scattering coefficient of+1 confusion region
Figure FDA00003427071000026
Then from one dimension echo data E The one dimension echoIn deduct the orientation and obtain residue signal E to the echoed signal of+1 confusion region correspondence Residue signal 1
According to the position of main imaging region, determine that the orientation is to the main imaging region central point O of+1 confusion region correspondence MainThe position, the orientation is designated as O to the central point of+1 confusion region Side+1, adopt coordinate to be expressed as to four summits of+1 confusion region in the orientation: the orientation is designated as AA to the upper left corner coordinate of+1 confusion region + 1(P + 1 left side, P On+1), the orientation is designated as BB to the upper right corner coordinate of+1 confusion region + 1(P + 1 right side, P On+1), the orientation is designated as CC to the lower left corner coordinate of+1 confusion region + 1(P + 1 left side, P + 1 time), the orientation is designated as DD to the lower right corner coordinate of+1 confusion region + 1(P + 1 right side, P + 1 time); Wherein:
Figure FDA00003427071000027
Figure FDA00003427071000029
The skew that the expression orientation makes progress in the orientation to the main imaging region center of+1 confusion region correspondence, and
Figure FDA00003427071000031
Figure FDA00003427071000032
The skew that the expression orientation makes progress in distance to the main imaging region center of+1 confusion region correspondence, and
Figure FDA00003427071000033
θ +Represent the orientation to the angle of squint of the main imaging region center correspondence of+1 confusion region correspondence, and
Figure FDA00003427071000034
R MinThe biography antenna is to the bee-line at main imaging region center under the expression; f pThe indicating impulse repetition frequency; θ represents the angle of squint at main imaging region center;
(D) in the orientation in-1 confusion region space orientation module (6) according to the position of main imaging region, determine the orientation to the position of-1 confusion region;
In-1 confusion region spatial discretization module (7), spatial discretization is carried out to-1 confusion region in the orientation in the orientation, make up and obtain the orientation to the measurement matrix Φ of the first confusion region correspondence -1
In the orientation to-1 confusion region inverting with eliminate in the obfuscation module (8) according to least square method from measuring matrix Φ -1With residue signal E Residue signal 1In obtain the orientation to the scattering coefficient of-1 confusion region
Figure FDA00003427071000035
Then from residue signal E Residue signal 1In deduct the orientation and obtain residue signal E to the echoed signal of-1 confusion region correspondence Residue signal 2
According to the position of main imaging region, determine the orientation to-1 confusion region (shown in Fig. 3 C) corresponding main imaging region central point O MainThe position, the orientation is designated as O to the central point of-1 confusion region Side-1, adopt coordinate to be expressed as to four summits of-1 confusion region in the orientation: the orientation is designated as AA to the upper left corner coordinate of-1 confusion region -1(P -1 left side, P On-1), the orientation is designated as BB to the upper right corner coordinate of-1 confusion region -1(P -1 right side, P On-1), the orientation is designated as CC to the lower left corner coordinate of-1 confusion region -1(P -1 left side, P -1 time), the orientation is designated as DD to the lower right corner coordinate of-1 confusion region -1(P -1 right side, P -1 time); Wherein,
Figure FDA000034270710000313
Figure FDA000034270710000314
Figure FDA000034270710000316
The skew that the expression orientation makes progress in the orientation to the main imaging region center of-1 confusion region correspondence, and
Figure FDA00003427071000036
Figure FDA00003427071000037
The skew that the expression orientation makes progress in distance to the main imaging region center of-1 confusion region correspondence, and
Figure FDA00003427071000038
θ -Represent the orientation to the angle of squint of the main imaging region center correspondence of-1 confusion region correspondence, and θ - = arcsin [ sin θ - λ + f p 2 × V ] ;
(E) in main imaging area spatial discretization module (9), main imaging region is carried out spatial discretization, make up the measurement matrix φ of main imaging region correspondence;
(F) in the orientation in the main imaging area processing module (10) according to measurement matrix φ and the residue signal E of least square method from main imaging region correspondence Residue signal 2Obtain final echo-wave imaging I as a result Abs
2. satellite-borne synthetic aperture radar according to claim 1 orientation suppresses system to the first confusion region echoed signal energy, it is characterized in that the step that adopts mesh generation that spatial discretization is carried out to+1 confusion region in the orientation has:
Step 401: adopt mesh generation to carry out the orientation and divide processing to the territory of+1 confusion region, described orientation is divided into to+1 confusion region
Figure FDA000034270710000310
Individual grid;
Step 402: the orientation to the orientation of+1 confusion region to net point
Figure FDA000034270710000311
With the distance to net point
Figure FDA000034270710000312
The grid element center coordinate be designated as
Figure FDA00003427071000041
Figure FDA00003427071000042
The expression orientation to the grid number, and The grid number rounds up; Expression distance to the grid number, and The grid number rounds up;
Step 403: at any time t of satellite platform flight kDown to the grid element center coordinate
Figure FDA00003427071000046
Carry out oblique distance and calculate, this oblique distance is
Figure FDA00003427071000047
Step 404: at t kIn time, inscribe, the grid element center coordinate
Figure FDA00003427071000048
Corresponding observation vector is
Figure FDA00003427071000049
Figure FDA000034270710000410
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure FDA000034270710000411
Distance to the 1st sampled point observed reading;
Figure FDA000034270710000412
Be illustrated in t kIn time, inscribe, the grid element center coordinate Distance to the 2nd sampled point observed reading;
Figure FDA000034270710000414
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure FDA000034270710000415
Distance to n sampled point observed reading;
Described in the orientation n sampled point observed reading in+1 confusion region The b frequency modulation rate of representing to transmit; I represents imaginary unit; N represents the grid element center coordinate
Figure FDA000034270710000417
Distance to sampled point; Represent distance to starting point, and
Figure FDA000034270710000419
Round up;
Figure FDA000034270710000420
Represent distance to terminal point, and
Figure FDA000034270710000421
Round downwards;
Figure FDA000034270710000422
Be illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to N DistanceIndividual sampled point observed reading;
Step 405: be structured in total movement time interior orientation to
Figure FDA000034270710000423
Individual, distance is to the
Figure FDA000034270710000424
The observation vector of individual grid element center correspondence, as follows:
Figure FDA00003427071000051
Wherein,
Figure FDA00003427071000052
Expression Transposition;
Step 406: make up the orientation and measure matrix to+1 confusion region, as follows:
Figure FDA00003427071000054
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region the observation vector to the 1st grid, distance to the 1st grid element center point; Be illustrated in total movement time interior orientation to the orientation of+1 confusion region the observation vector to the 1st grid, distance to the 2nd grid element center point;
Figure FDA00003427071000057
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to the 1st grid, distance to the
Figure FDA00003427071000058
The observation vector of individual grid element center point;
Figure FDA00003427071000059
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region the observation vector to the 2nd grid, distance to the 1st grid element center point;
Figure FDA000034270710000510
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region the observation vector to the 2nd grid, distance to the 2nd grid element center point;
Figure FDA000034270710000511
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to the 2nd grid, distance to the
Figure FDA000034270710000512
The observation vector of individual grid element center point;
Figure FDA000034270710000513
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to Individual grid, distance are to the observation vector of the 1st grid element center point;
Figure FDA000034270710000515
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to Individual grid, distance are to the observation vector of the 2nd grid element center point;
Figure FDA000034270710000517
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to
Figure FDA000034270710000521
Individual grid, distance are to the
Figure FDA000034270710000518
The observation vector of individual grid element center point;
According to least square method
Figure FDA000034270710000519
Calculate the orientation to the scattering coefficient of+1 confusion region
Figure FDA000034270710000520
As follows:
Figure FDA00003427071000061
α 1-1The expression orientation is the scattering coefficient to the 1st grid, distance to the 1st grid element center point to the orientation of+1 confusion region; α 1-2The expression orientation is the scattering coefficient to the 1st grid, distance to the 2nd grid element center point to the orientation of+1 confusion region
Figure FDA00003427071000062
The expression orientation to the orientation of+1 confusion region to the 1st grid, distance to the
Figure FDA00003427071000063
The scattering coefficient of individual grid element center point; α 2-1The expression orientation is the scattering coefficient to the 2nd grid, distance to the 1st grid element center point to the orientation of+1 confusion region; α 2-2The expression orientation is the scattering coefficient to the 2nd grid, distance to the 2nd grid element center point to the orientation of+1 confusion region;
Figure FDA00003427071000064
The expression orientation to the orientation of+1 confusion region to the 2nd grid, distance to the
Figure FDA00003427071000065
The scattering coefficient of individual grid element center point;
Figure FDA00003427071000066
The expression orientation to the orientation of+1 confusion region to the
Figure FDA00003427071000067
Individual grid, distance are to the scattering coefficient of the 1st grid element center point;
Figure FDA00003427071000068
The expression orientation to the orientation of+1 confusion region to the
Figure FDA00003427071000069
Individual grid, distance are to the scattering coefficient of the 2nd grid element center point;
Figure FDA000034270710000610
The expression orientation to the orientation of+1 confusion region to the
Figure FDA000034270710000611
Individual grid, distance are to the
Figure FDA000034270710000612
The scattering coefficient of individual grid element center point; From echo data E The one dimension echoThe middle orientation of separating obtains residue signal E to the echoed signal of+1 confusion region correspondence Residue signal 1For
Figure FDA000034270710000613
3. satellite-borne synthetic aperture radar according to claim 1 orientation suppresses system to the first confusion region echoed signal energy, and it is characterized in that adopting mesh generation that the spatial discretization step is carried out to-1 confusion region in the orientation has:
Step 701: adopt mesh generation to carry out the orientation and divide processing to the territory of-1 confusion region, described orientation is divided into to-1 confusion region Individual grid;
Step 702: the orientation to the orientation of-1 confusion region to net point
Figure FDA000034270710000615
With the distance to net point
Figure FDA000034270710000616
The grid element center coordinate be designated as
Figure FDA000034270710000617
Figure FDA000034270710000618
The expression orientation to the grid number, and
Figure FDA000034270710000619
The grid number rounds up;
Figure FDA00003427071000071
Expression distance to the grid number, and The grid number rounds up;
Step 703: at any time t of satellite platform flight kDown to the grid element center coordinate
Figure FDA00003427071000073
Carry out oblique distance and calculate, this oblique distance is
Figure FDA00003427071000074
Step 704: at t kIn time, inscribe, the grid element center coordinate
Figure FDA00003427071000075
Corresponding observation vector is
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure FDA00003427071000078
Distance to the 1st sampled point observed reading;
Figure FDA00003427071000079
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure FDA000034270710000710
Distance to the 2nd sampled point observed reading;
Figure FDA000034270710000711
Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure FDA000034270710000712
Distance to n sampled point observed reading;
Described in the orientation n sampled point observed reading in-1 confusion region
Figure FDA000034270710000713
N represents the grid element center coordinate
Figure FDA000034270710000714
Distance to sampled point;
Figure FDA000034270710000715
Represent distance to starting point, and
Figure FDA000034270710000716
Round up;
Figure FDA000034270710000717
Represent distance to terminal point, and
Figure FDA000034270710000718
Round downwards; Be illustrated in t kIn time, inscribe, the grid element center coordinate
Figure FDA000034270710000720
Distance to N DistanceIndividual sampled point observed reading;
Step 705: be structured in total movement time interior orientation to
Figure FDA000034270710000721
Individual, distance is to the The observation vector of individual grid element center correspondence, as follows:
Figure FDA000034270710000723
Wherein,
Figure FDA000034270710000724
Expression Transposition;
Step 706: make up the orientation and measure matrix to-1 confusion region, as follows:
Figure FDA00003427071000081
Figure FDA00003427071000082
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region the observation vector to the 1st grid, distance to the 1st grid element center point;
Figure FDA00003427071000083
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region the observation vector to the 1st grid, distance to the 2nd grid element center point; Be illustrated in total movement time interior orientation to the orientation of-1 confusion region to the 1st grid, distance to the The observation vector of individual grid element center point;
Figure FDA00003427071000086
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region the observation vector to the 2nd grid, distance to the 1st grid element center point;
Figure FDA00003427071000087
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region the observation vector to the 2nd grid, distance to the 2nd grid element center point;
Figure FDA00003427071000088
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region to the 2nd grid, distance to the The observation vector of individual grid element center point;
Figure FDA000034270710000810
Be illustrated in total movement time interior orientation to the orientation of-1 confusion region to
Figure FDA000034270710000811
Individual grid, distance are to the observation vector of the 1st grid element center point;
Figure FDA000034270710000812
Be illustrated in total movement time interior orientation to the orientation of+1 confusion region to
Figure FDA000034270710000813
Individual grid, distance are to the observation vector of the 2nd grid element center point; Be illustrated in total movement time interior orientation to the orientation of-1 confusion region to
Figure FDA000034270710000815
Individual grid, distance are to the
Figure FDA000034270710000816
The observation vector of individual grid element center point;
According to least square method
Figure FDA000034270710000817
Calculate the orientation to the scattering coefficient of-1 confusion region As follows:
Figure FDA000034270710000819
β 1-1The expression orientation is the scattering coefficient to the 1st grid, distance to the 1st grid element center point to the orientation of-1 confusion region; β 1-2The expression orientation is the scattering coefficient to the 1st grid, distance to the 2nd grid element center point to the orientation of-1 confusion region;
Figure FDA00003427071000091
The expression orientation to the orientation of-1 confusion region to the 1st grid, distance to the
Figure FDA00003427071000092
The scattering coefficient of individual grid element center point; β 2-1The expression orientation is the scattering coefficient to the 2nd grid, distance to the 1st grid element center point to the orientation of-1 confusion region; β 2-2The expression orientation is the scattering coefficient to the 2nd grid, distance to the 2nd grid element center point to the orientation of-1 confusion region;
Figure FDA00003427071000093
The expression orientation to the orientation of-1 confusion region to the 2nd grid, distance to the
Figure FDA00003427071000094
The scattering coefficient of individual grid element center point;
Figure FDA00003427071000095
The expression orientation to the orientation of-1 confusion region to the Individual grid, distance are to the scattering coefficient of the 1st grid element center point; The expression orientation to the orientation of-1 confusion region to the
Figure FDA00003427071000098
Individual grid, distance are to the scattering coefficient of the 2nd grid element center point; The expression orientation to the orientation of-1 confusion region to the
Figure FDA000034270710000910
Individual grid, distance are to the
Figure FDA000034270710000911
The scattering coefficient of individual grid element center point;
From residue signal E Residue signal 1The middle orientation of separating obtains residue signal E to the echoed signal of-1 confusion region correspondence Residue signal 2For
4. satellite-borne synthetic aperture radar according to claim 1 orientation suppresses system to the first confusion region echoed signal energy, and it is characterized in that adopting mesh generation that main imaging region is carried out the spatial discretization step has:
Step 901: the zone of adopting mesh generation to carry out main imaging area is divided and is handled, and described main imaging region is divided into D The orientation* D DistanceIndividual grid;
Step 902: the orientation of main imaging region is to net point I aWith the distance to net point I rThe grid element center coordinate be designated as
Figure FDA000034270710000913
D The orientationThe expression orientation to the grid number, and
Figure FDA000034270710000914
The grid number rounds up; D DistanceExpression distance to the grid number, and
Figure FDA000034270710000915
The grid number rounds up;
Step 903: at any time t of satellite platform flight kDown to grid element center coordinate W GridCarry out oblique distance and calculate, this oblique distance is
Figure FDA000034270710000916
Step 904: at t kIn time, inscribe, grid element center coordinate W GridCorresponding observation vector is e K-1Be illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to the 1st sampled point observed reading; e K-2Be illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to the 2nd sampled point observed reading; e K-nBe illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to n sampled point observed reading;
Described n sampled point observed reading in main imaging area
Figure FDA00003427071000101
N represents grid element center coordinate W GridDistance to sampled point; N BeginRepresent distance to starting point, and
Figure FDA00003427071000102
Round up; N EndRepresent distance to terminal point, and
Figure FDA00003427071000103
Round downwards; Be illustrated in t kIn time, inscribe, grid element center coordinate W GridDistance to N DistanceIndividual sampled point observed reading;
Step 905: be structured in total movement time interior orientation to I aIndividual, the distance to I rThe observation vector of individual grid element center correspondence, as follows:
Figure FDA00003427071000105
Wherein,
Figure FDA00003427071000106
Expression
Figure FDA00003427071000107
Transposition;
Step 906: make up main imaging region and measure matrix, as follows:
Figure FDA00003427071000108
5. satellite-borne synthetic aperture radar according to claim 4 orientation suppresses system to the first confusion region echoed signal energy, it is characterized in that: according to least square method E The scattering of one dimension primary area=(φ H* φ) -1* φ H* E Residue signal 2Calculate main imaging region scattering coefficient one-dimensional representation E The scattering of one dimension primary area
Described main imaging region scattering coefficient one dimension
Figure FDA00003427071000111
σ 1-1Represent that the orientation of main imaging area is to the 1st grid, apart from the scattering coefficient to the 1st grid element center point; σ 1-2Represent that the orientation of main imaging area is to the 1st grid, apart from the scattering coefficient to the 2nd grid element center point;
Figure FDA00003427071000112
The orientation of representing main imaging area to the 1st grid, distance to D DistanceThe scattering coefficient of individual grid element center point; σ 2-1Represent that the orientation of main imaging area is to the 2nd grid, apart from the scattering coefficient to the 1st grid element center point; σ 2-2Represent that the orientation of main imaging area is to the 2nd grid, apart from the scattering coefficient to the 2nd grid element center point;
Figure FDA00003427071000113
The orientation of representing main imaging area to the 2nd grid, distance to D DistanceThe scattering coefficient of individual grid element center point;
Figure FDA00003427071000114
Represent that the orientation of main imaging area is to D The orientationIndividual grid, distance are to the scattering coefficient of the 1st grid element center point;
Figure FDA00003427071000115
Represent that the orientation of main imaging area is to D The orientationIndividual grid, distance are to the scattering coefficient of the 2nd grid element center point;
Figure FDA00003427071000116
Represent that the orientation of main imaging area is to D The orientationIndividual grid, distance are to D DistanceThe scattering coefficient of individual grid element center point;
With E The scattering of one dimension primary areaBe converted to two-dimensional matrix E by one-dimensional vector Two dimension primary area scattering, can be expressed as:
Figure FDA00003427071000117
To matrix E Two dimension primary area scatteringAsk absolute value, obtain the imaging results of imaging region, be expressed as:
Figure FDA00003427071000118
CN 201210149529 2012-05-14 2012-05-14 Energy inhibiting system for azimuth first fuzzy region echo signal of space-borne synthetic aperture radar (SAR) Expired - Fee Related CN102680956B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210149529 CN102680956B (en) 2012-05-14 2012-05-14 Energy inhibiting system for azimuth first fuzzy region echo signal of space-borne synthetic aperture radar (SAR)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210149529 CN102680956B (en) 2012-05-14 2012-05-14 Energy inhibiting system for azimuth first fuzzy region echo signal of space-borne synthetic aperture radar (SAR)

Publications (2)

Publication Number Publication Date
CN102680956A CN102680956A (en) 2012-09-19
CN102680956B true CN102680956B (en) 2013-08-14

Family

ID=46813171

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210149529 Expired - Fee Related CN102680956B (en) 2012-05-14 2012-05-14 Energy inhibiting system for azimuth first fuzzy region echo signal of space-borne synthetic aperture radar (SAR)

Country Status (1)

Country Link
CN (1) CN102680956B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103605119B (en) * 2013-09-18 2015-12-02 北京航空航天大学 A kind of satellite-borne synthetic aperture radar band pattern azimuth ambiguity suppression method
CN106093945B (en) * 2016-08-22 2018-04-06 西安电子科技大学 A kind of SAR radar imaging methods for removing orientation ghost image
CN108399607B (en) * 2017-11-21 2021-08-17 北京航空航天大学 Method, device and equipment for inhibiting image orientation blur and computer readable storage medium
CN112198510B (en) * 2020-08-28 2023-06-06 中国科学院空天信息创新研究院 Method for inhibiting azimuth ambiguity of dual-receiving-channel hybrid full-polarization radar system
CN118169688A (en) * 2024-05-16 2024-06-11 中国电子科技集团公司第十四研究所 Channel error correction method and device for satellite-borne multichannel SAR

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4853699A (en) * 1987-11-13 1989-08-01 Hughes Aircraft Company Method for cancelling azimuth ambiguity in a SAR receiver
CN101776755A (en) * 2010-02-05 2010-07-14 西安电子科技大学 Method for preprocessing high speed platform ultra-high resolution stop and reserves (SAR) data
CN101957449A (en) * 2010-04-06 2011-01-26 北京航空航天大学 Optimization method for azimuth ambiguity in space-borne TOPSAR mode

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008261720A (en) * 2007-04-12 2008-10-30 Mitsubishi Electric Corp Ambiguity processing device

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4853699A (en) * 1987-11-13 1989-08-01 Hughes Aircraft Company Method for cancelling azimuth ambiguity in a SAR receiver
CN101776755A (en) * 2010-02-05 2010-07-14 西安电子科技大学 Method for preprocessing high speed platform ultra-high resolution stop and reserves (SAR) data
CN101957449A (en) * 2010-04-06 2011-01-26 北京航空航天大学 Optimization method for azimuth ambiguity in space-borne TOPSAR mode

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Albert0 Moreira.Suppressing the Azimuth Ambiguities in Synthetic Aperture Radar Images.《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》.1993,第31卷(第4期),885-895.
Andrea Monti Guarnieri.Adaptive Removal of Azimuth Ambiguities in SAR Images.《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》.2005,第43卷(第3期), *
JP特开2008-261720A 2008.10.30
Suppressing the Azimuth Ambiguities in Synthetic Aperture Radar Images;Albert0 Moreira;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;19930731;第31卷(第4期);885-895 *
于泽等.基于LCMV算法的星载SAR距离模糊抑制方法.《北京航空航天大学学报》.2006,第32卷(第9期),
基于LCMV算法的星载SAR距离模糊抑制方法;于泽等;《北京航空航天大学学报》;20060930;第32卷(第9期);1035-1038 *

Also Published As

Publication number Publication date
CN102680956A (en) 2012-09-19

Similar Documents

Publication Publication Date Title
CN100470255C (en) Single-channel synthetic aperture radar moving-target detection method based on multi-apparent subimage paire
CN105785327B (en) Frequency diversity array synthetic aperture radar high-resolution wide swath imaging method
US8861588B2 (en) Apparatus and method for sampling and reconstruction of wide bandwidth signals below Nyquist rate
CN102680956B (en) Energy inhibiting system for azimuth first fuzzy region echo signal of space-borne synthetic aperture radar (SAR)
CN102183762A (en) Method for acquiring and imaging data of compressive sensing synthetic aperture radar
CN107037429B (en) Linear array SAR three-dimensional imaging method based on threshold gradient tracking algorithm
CN103149561B (en) Microwave imaging method based on scenario block sparsity
CN102323583B (en) Super-resolution linear array three-dimensional synthetic aperture radar imaging method
Gogineni et al. Bed topography of Jakobshavn Isbræ, Greenland, and Byrd Glacier, Antarctica
CN106970386A (en) A kind of optimization method of RADOP beam sharpening
CN108226927B (en) SAR imaging method based on weighted iteration minimum sparse Bayesian reconstruction algorithm
CN106405548A (en) Inverse synthetic aperture radar imaging method based on multi-task Bayesian compression perception
CN102221697A (en) Airborne multi-antenna SAR chromatography three dimensional imaging system and imaging method thereof
CN102914773B (en) Multi-pass circumference SAR three-dimensional imaging method
CN104950305A (en) Real beam scanning radar angle super-resolution imaging method based on sparse constraint
CN105676190B (en) A kind of method and apparatus of correction synthetic aperture radar echo data
CN107192992B (en) Radar calibration method, calibrator, radar device and system
CN102721965B (en) Omni-directional side lobe suppression method for synthetic aperture radar imaging processing
CN102998673A (en) Compressive sensing imaging method for synthetic aperture radar
CN103630903B (en) The method of flow field, sea radial velocity is measured based on straight rail interference SAR
CN104678386A (en) Method for detecting target by utilizing relevant power of GNSS(global navigation satellite system) sea surface reflection signal
Cheng et al. Multipath scattering of typical structures in urban areas
Jezek et al. Remote sensing of sea ice thickness and salinity with 0.5–2 GHz microwave radiometry
KR102151362B1 (en) Image decoding apparatus based on airborn using polar coordinates transformation and method of decoding image using the same
CN102890277A (en) Range migration imaging method of shift invariant bi-static synthetic aperture radar

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: 20130814