CN103792536B - A kind of satellite-borne synthetic aperture radar slip beam bunching mode orientation is to parameter acquiring method - Google Patents

A kind of satellite-borne synthetic aperture radar slip beam bunching mode orientation is to parameter acquiring method Download PDF

Info

Publication number
CN103792536B
CN103792536B CN201410041182.7A CN201410041182A CN103792536B CN 103792536 B CN103792536 B CN 103792536B CN 201410041182 A CN201410041182 A CN 201410041182A CN 103792536 B CN103792536 B CN 103792536B
Authority
CN
China
Prior art keywords
orientation
point target
centerdot
oblique distance
data
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.)
Active
Application number
CN201410041182.7A
Other languages
Chinese (zh)
Other versions
CN103792536A (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.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
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 Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201410041182.7A priority Critical patent/CN103792536B/en
Publication of CN103792536A publication Critical patent/CN103792536A/en
Application granted granted Critical
Publication of CN103792536B publication Critical patent/CN103792536B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • G01S13/9052Spotlight mode
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system
    • G01S7/4017Means for monitoring or calibrating of parts of a radar system of HF systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes
    • G01S13/9047Doppler beam sharpening mode

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The present invention provides a kind of satellite-borne synthetic aperture radar slip beam bunching mode orientation to parameter acquiring method, it is to utilize orientation to initial imaging belt length, antenna beamwidth, initial angle of strabismus, terminate angle of strabismus and build, with target oblique distance, the reference oblique distance model that wave beam center of rotation is corresponding, obtain and according to the effective irradiation time model corresponding with reference to one point target of oblique distance data construct, obtain and build effective Doppler frequency corresponding to this point target and doppler bandwidth model according to the most initial irradiation time of this point target with effectively terminating irradiation time, obtain and according to initial effective Doppler frequency corresponding to this point target, terminate effective Doppler frequency and doppler bandwidth, resolution descent coefficient builds the azimuth resolution model that this point target is corresponding, obtain the azimuth resolution that this point target is corresponding;Maximum orientation is then exported to imaging belt length data when azimuth resolution is more than resolution index request.

Description

A kind of satellite-borne synthetic aperture radar slip beam bunching mode orientation is to parameter acquiring method
Technical field
The invention belongs to satellite-borne synthetic aperture radar imaging field, relate to satellite-borne synthetic aperture radar (SAR) acquiring technology of systematic parameter.
Background technology
Satellite-borne synthetic aperture radar slip beam bunching mode is to study and realized high-resolution width in recent years to cover The Main Patterns of lid satellite-borne synthetic aperture radar imaging, and defended at Germany TerraSAR-X in 2007 Obtain imaging applications on star first, obtain a large amount of high-quality satellite-borne synthetic aperture radar image.But It is that TerraSAR-X is limited to scan capability due to antenna bearingt, its resolution and imaging belt length Index Design the highest (1m × 5km), by increasing scan capability and the power supply of satellite of antenna Ability, the performance of satellite-borne synthetic aperture radar slip beam bunching mode can be greatly improved that (resolution is excellent In 0.3m, imaging belt length is better than 10km).But performance raising can cause parameter designing, echo Signal simulation is greatly increased with imaging processing difficulty and complexity.
Satellite-borne synthetic aperture radar slip beam bunching mode distance is to the basic phase of parameter designing and band pattern With, but orientation is more complex to parameter designing.Satellite-borne synthetic aperture radar slip beam bunching mode side Position to parameter designing mainly to resolution, imaging belt length, orientation to beam angle, orientation to Initial angle of strabismus and orientation trade off analysis, to obtain work to terminating restriction relation between angle of strabismus The high-resolution width that can realize in journey cover diameter radar image azimuth resolution and orientation to Re-imaging length index.The parameter designing result the most accurately ground integrated echo-signal of star to be passed through Emulation, imaging processing and point target Performance Evaluation are verified, are better than 0.3m high-resolution spaceborne Synthetic aperture radar slip beam bunching mode echo simulation implements the most multiple with image-processing algorithms Miscellaneous, operand is big, the longest.Parameter designing process then can be when tentative programme designs accurately First avoid echo simulation and imaging process, and can ensure that parameter designing reliable results, carry The efficiency of high conceptual design.
The satellite-borne synthetic aperture radar slip beam bunching mode orientation occurred in the outer document of Present Domestic is to ginseng Number design processes are excessively simple, have following 2 deficiencies for engineer applied:
(1) calculate and time and the change relevant with Doppler frequency with single star ground relative velocity Amount, this approximation causes the calculated azimuth resolution of simple parameter and star ground integrated emulation knot Fruit differs greatly, and makes parameter designing result reliability reduce.
(2) orientation to antenna radiation pattern, imaging processing weighting, Doppler's parameter estimate error and becomes The azimuth resolution decline caused as algorithm approximation four factors does not has in relevant computing formula Embodying, discounting for the impact of these four factor, the resolution that simple parameter designing obtains exists Practical Project possibly cannot realize.
In satellite-borne SAR slip beam bunching mode engineering design process azimuth resolution with orientation to becoming Image-tape length is index request, user propose;Antenna bearingt to beam angle by radar operating wave Long and antenna bearingt determines to size, it is contemplated that the antenna length of satellite-borne synthetic aperture radar satellite is past Toward being determined by band pattern azimuth resolution, at this it is believed that orientation is known to antenna beamwidth; Antenna bearingt is determined to scan capability by antenna bearingt to initial and termination scan angle.In engineering, Making demands slip beam bunching mode orientation to antenna scanning ability is topmost design requirement.
Summary of the invention
(1) to solve the technical problem that
Complicated in order to solve prior art computing, length computationally intensive, time-consuming, inefficient defect, The purpose of the present invention is to propose to a kind of spaceborne slip beam bunching mode orientation of engineer applied that is applicable to one-tenth Image-tape Design of length flow process, can first avoid echo simulation and imaging processing when tentative programme designs Process, and can ensure that the reliable results that gets parms, improve efficiency getparms.
(2) technical scheme
In order to realize the purpose of the present invention, it is poly-that the present invention provides a kind of satellite-borne synthetic aperture radar to slide Bundle mode orientation is to parameter acquiring method, and it is as follows that the method comprising the steps of:
Step S1: set orientation to initial imaging belt length data, orientation to antenna beamwidth, Orientation to initial angle of strabismus data, orientation to terminating angle of strabismus data and target place oblique distance;
Step S2: utilize orientation to initial imaging belt length data, orientation to antenna beamwidth Data, orientation to initial angle of strabismus data, orientation to terminating angle of strabismus data and target oblique distance structure point The reference oblique distance model that wave beam center of rotation during target oblique distance is corresponding, it is thus achieved that during point target oblique distance The reference oblique distance data that wave beam center of rotation is corresponding;
Step S3: according to the reference oblique distance data that wave beam center of rotation during point target oblique distance is corresponding Building oblique distance is effective irradiation time model corresponding to the point target of R, it is thus achieved that this point target Effectively initial radiation time data and this point target effectively terminate radiation time data;
Step S4: building oblique distance according to the most initial described irradiation time and termination irradiation time is Effective Doppler Frequency Model corresponding to a point target and doppler bandwidth model, obtain this Initial effective doppler frequency data corresponding to point target, terminate effective doppler frequency data and many General Le band data;
Step S5: according to initial effective doppler frequency data and the effective Doppler frequency number of termination According to, the resolution descent coefficient that causes of doppler bandwidth data, antenna radiation pattern, imaging processing add Resolution descent coefficient that the resolution descent coefficient that causes of power, Doppler's parameter estimate error cause, It is corresponding that the resolution descent coefficient that imaging algorithm approximation causes builds the point target that oblique distance is R Azimuth resolution model, it is thus achieved that the azimuth resolution that this point target is corresponding;
Step S6: judge whether azimuth resolution meets less than resolution index request, if Show that orientation set in advance is less than normal to imaging belt length data, then perform step S1 and reset Orientation, to imaging belt length data value, repeats step S1~S5 until azimuth resolution is more than dividing Resolution index request, then the maximum orientation of output is to imaging belt length data.
(3) beneficial effect
The present invention provides satellite-borne synthetic aperture radar slip beam bunching mode certain to scan angle in orientation Under conditions of obtain the method for imaging belt length meeting resolution requirement, profit in this way can be straight Connect the orientation obtained to initial angle of strabismus, terminate angle of strabismus, beam angle, resolution and imaging belt Restriction relation between length, it is to avoid complicated echo simulation and imaging processing proof procedure, simplifies Satellite-borne synthetic aperture radar slip beam bunching mode conceptual design flow process, improves engineering design efficiency, tool Higher engineer applied is had to be worth.
Accompanying drawing explanation
Fig. 1 is satellite-borne synthetic aperture radar slip beam bunching mode operation principle rough schematic view.
Fig. 2 is that satellite-borne synthetic aperture radar slip beam bunching mode orientation obtains to imaging belt length Flow chart.
Detailed description of the invention
Below in conjunction with appended drawings 1 and Fig. 2, the present invention is described in detail.
The ultimate principle of satellite-borne SAR slip beam bunching mode is as it is shown in figure 1, at one section of imaging time In, satellite points to beam center with certain speed along orbital flight, control azimuth, makes orientation To beam center with certain initial angle of strabismus ψstartStarting to rotate around rotary middle point, angle of strabismus reaches To terminating angle of strabismus ψendTime terminate.In engineering, rotary middle point is be under earth's surface Virtual point, the oblique distance that center of rotation is corresponding is more than the oblique distance at impact point place, on the ground, wave beam Center is moved from left to right, slips over all targets of imaging region successively, and in figure, middle heavy black is Imaging region.In Fig. 1, in imaging region, solid large circle point is any point target, and this point target is permissible It is in imaging region any position.Definition oblique distance is that any point target of R is in 3dB beam angle The interior time period is effective irradiation time, is in the moment of 3dB wave beam left side edge for initial effective Irradiate the moment, be in the moment of 3dB wave beam right side edge for terminating effectively irradiating the moment.To being in The point target echo data of different irradiation time sections carries out imaging processing, can obtain meeting certain resolution The point target image that rate requires.
In Fig. 1 and Fig. 2, each variable-definition is as follows:
ψstartFor orientation to initial angle of strabismus, it is defined as orientation of Polaroid period in home beam The heart points to.
ψendFor orientation to terminating angle of strabismus, it is defined as orientation of Polaroid period to terminating in wave beam The heart points to.
θaFor orientation to beam angle, it is defined as antenna bearingt to transmitting-receiving equivalent directions figure half-power width Degree.
WsceneFor orientation to imaging belt length, it is defined as orientation to the one-tenth meeting resolution index request As region.
R is point target oblique distance, and during being defined as imaging, radar is to the nearest oblique distance of target.
RrotFor center of rotation oblique distance, during being defined as imaging, radar is to the nearest oblique distance of center of rotation.
XpFor any point target imaging region location variable (relative to imaging region center, And to define imaging region center location variable be 0).
taFor the time variable of Polaroid period satellite position, definition distance by radar center of rotation is Moment time near is 0.
tA, start, 0For the most initial irradiation time t that any point target is correspondingA, start, 0, it is defined as at target In moment of beam angle left hand edge (relative to ta=0)。
tA, end, 0For effectively termination irradiation time t that any point target is correspondingA, end, 0, it is defined as at target In moment of beam angle right hand edge (relative to ta=0)。
ψStart, 0It is defined as any point target effective and initiates orientation corresponding to irradiation time to angle of strabismus.
ψEnd, 0It is defined as any point target effective and terminates orientation corresponding to irradiation time to angle of strabismus.
fA, start, 0Initiate effective Doppler frequency for any point target, be defined as any point target and initiate The Doppler frequency that effectively irradiation time is corresponding.
fA, end, 0For the effective Doppler frequency of any point target termination, it is defined as any point target termination The Doppler frequency that effectively irradiation time is corresponding.
fdFor the doppler bandwidth that any point target is corresponding, it is defined as fd=fA, start, 0-fA, end, 0
ρaFor azimuth resolution.Resolution is defined as the orientation 3dB to point target impulse response Width.
k1The resolution descent coefficient caused is weighted for antenna radiation pattern.
k2The resolution descent coefficient caused for imaging processing weighting (for forcing down secondary lobe).
k3The resolution descent coefficient caused for Doppler's parameter estimate error.
k4The resolution descent coefficient caused is approximated for imaging algorithm.
k1、k2、k3And k4Being the parameter relevant with imaging performance, its value thinks known at this.
VgFor wave beam ground speed (for scalar).
VsFor satellite velocities (for scalar).
λ is the wavelength that synthetic aperture radar launches signal.
The echo samples initial time worked according to satellite orbit parameter and satellite-borne synthetic aperture radar, Satellite velocities V during the most a certain location point of satellite transit can be calculateds, star ground relatively Speed VgWith oblique distance R of target to radar, to Vs、VgUsed with the calculating process of R, recognized at this For known.
In orientation to angle of strabismus ψstart, angle of strabismus ψendWith 3dB beam angle θaUnder conditions of Yi Ding, Obtain orientation as follows to the calculation procedure of the maximum imaging belt length meeting resolution requirement:
Step S1: set orientation to initial imaging belt length data Wscene, orientation is to antenna beam Width θa, orientation is to initial angle of strabismus data ψstart, orientation is to terminating angle of strabismus data ψendAnd mesh Mark place oblique distance R;
Step S2: according to geometrical relationship between target and radar during satellite flight, utilizes Orientation is to initial imaging belt length data Wscene, orientation is to antenna 3dB beam angle data θa, side Position is to initial angle of strabismus data ψstart, orientation is to terminating angle of strabismus data and target oblique distance R ψend, Build reference oblique distance R that when point target oblique distance is R, wave beam center of rotation is correspondingrot, computing formula is such as Under:
R rot = R · [ tan ( ψ end + θ a / 2 ) - tan ( ψ start - θ a / 2 ) ] - W scene [ tan ( ψ end ) - tan ( ψ start ) ] - - - ( 1 )
Step S3: according to the reference oblique distance data that wave beam center of rotation during point target oblique distance R is corresponding Building oblique distance is effective irradiation time model corresponding to the point target of R, it is thus achieved that this point target Effectively initial radiation time data tA, start, 0With this point target effectively terminate radiation time data tA, end, 0;Calculating oblique distance is the most initial irradiation time t corresponding to the point target of RA, start, 0And end Only irradiation time tA, end, 0, computing formula is as follows:
V g · t a , start , 0 R rot = tan ( ψ start , 0 ) V g · t a , start , 0 - X p R rot = tan ( ψ start , 0 - θ a / 2 ) ⇒ tan ( ψ start , 0 ) = X p 2 · R rot + R / R rot - 1 2 · tan ( θ a / 2 ) + [ X p 2 · R rot + R / R rot - 1 2 · tan ( θ a / 2 ) ] 2 - R R rot + X p R rot · tan ( θ a / 2 ) ⇒ t a , start , 0 = R rot · tan ( ψ start , 0 ) V g - - - ( 2 )
V g · t a , end , 0 R rot = tan ( ψ end , 0 ) V g · t a , end , 0 - X p R rot = tan ( ψ end , 0 + θ a / 2 ) ⇒ tan ( ψ end , 0 ) = X p 2 · R rot + 1 - R / R rot 2 · tan ( θ a / 2 ) + [ X p 2 · R rot + 1 - R / R rot 2 · tan ( θ a / 2 ) ] 2 - R R rot + X p R rot · tan ( θ a / 2 ) ⇒ t a , end , 0 = R rot · tan ( ψ end , 0 ) V g - - - ( 3 )
Wherein: this point target is in location variable X of imaging regionpIt is expressed asSolve equation (2) and obtain tA, start, 0Data, solve equation (3) and obtain tA, end, 0 Data.
In the present invention, point target effectively initiates irradiation time and uses star in the calculating terminating irradiation time Ground ground velocity relatively, makes the calculating about the time the most accurate.General document uses star ground relatively The speed approximate calculation time, introduce bigger deviation.
Step S4: according to the most initial described irradiation time tA, start, 0With termination irradiation time tA, end, 0Structure Build effective Doppler Frequency Model corresponding to a point target that oblique distance is R and doppler bandwidth mould Type, obtains initial effective doppler frequency data f that this point target is correspondingA, start, 0, terminate effectively Doppler frequency data fA, end, 0With doppler bandwidth data fd;Calculate corresponding the rising of any point target Begin effective Doppler frequency fA, start, 0With effective Doppler frequency f of terminationA, end, 0With doppler bandwidth fd
f a , start , 0 = - 2 · V s λ · sin ( ψ start , 0 - θ a / 2 ) = - 2 · V s λ · V g · t a , start , 0 - X p R 2 + ( V g · t a , start , 0 - X p ) 2 - - - ( 4 )
f a , end , 0 = - 2 · V s λ · sin ( ψ start , 0 + θ a / 2 ) = - 2 · V s λ · V g · t a , end , 0 - X p R 2 + ( V g · t a , end , 0 - X p ) 2 - - - ( 5 )
f d = f a , start , 0 - f a , end , 0 = - 2 · V s λ [ V g · t a , start , 0 - X p R 2 + ( V g · t a , start , 0 - X p ) 2 - V g · t a , end , 0 - X p R 2 + ( V g · t a , end , 0 - X p ) 2 ] - - - ( 6 )
Initial effective Doppler frequency and the effective Doppler frequency of termination that point target is corresponding calculate and use The speed of satellite, makes the calculating about Doppler frequency more accurate.General document uses star ground Relative velocity approximate calculation Doppler frequency, introduces bigger deviation.
Step S5: according to initial effective doppler frequency data fA, start, 0With termination effective Doppler frequency Rate data fA, end, 0, doppler bandwidth data fd, the resolution descent coefficient k that causes of antenna radiation pattern1、 Imaging processing weights the resolution descent coefficient k caused2, Doppler's parameter estimate error cause point Resolution descent coefficient k3, imaging algorithm approximate the resolution descent coefficient k that causes4Structure oblique distance is R Azimuth resolution model corresponding to a point target, it is thus achieved that orientation corresponding to this point target to Resolution;Simplify geometric model according to satellite-borne synthetic aperture radar slip beam bunching mode, calculate acquisition side Position oblique distance in imaging belt is the azimuth resolution that any point target of R is corresponding, is simultaneously introduced sky The weighting of line directional diagram, imaging processing, Doppler's parameter estimate error and imaging algorithm approximation cause Azimuth resolution broadening.Described azimuth resolution model ρaIt is expressed as follows:
ρ a = k 1 · k 2 · k 3 · k 4 · V g f d = - k 1 · k 2 · k 3 · k 4 · λ · V g 2 · V s [ V g · t a , start , 0 - X p R 2 + ( V g · t a , start , 0 - X p ) 2 - V g · t a , start , 0 - X p R 2 + ( V g · t a , end , 0 - X p ) 2 ] - 1 - - - ( 7 )
It is to set satellite in a period of time that described satellite-borne SAR slip beam bunching mode simplifies geometric model Flight path is approximately straight line, and point target is static target;During imaging, control azimuth refers to wave beam To, make beam center point to a virtual point below away from imaging region, described virtual point all the time Being wave beam center of rotation, virtual point is in below ground;Wave beam the most slowly slips over imaging area Territory WXscene, in imaging region, oblique distance is that point target X of R goes through 3dB beam angle, from direction Figure right side edge enters into left side edge and leaves, initiateing during point target approach axis figure right side edge Angle of strabismus is ψStart, 0, termination angle of strabismus during point target departure direction figure left side edge is ψEnd, 0.Institute State orientation to initial angle of strabismus data ψstartWith orientation to terminating angle of strabismus data ψendSize and just Negative not restriction, therefore, it is adaptable to positive side-looking satellite-borne SAR slip beam bunching mode, is also applied for Stravismus satellite-borne SAR slip beam bunching mode.
In the present invention, azimuth resolution corresponding to point target calculates and introduces what antenna radiation pattern caused Azimuth resolution descent coefficient k1, imaging processing weight the azimuth resolution descent coefficient that causes k2, the azimuth resolution descent coefficient k that causes of Doppler's parameter estimate error3, and imaging algorithm The azimuth resolution descent coefficient k that approximation causes4.In engineering, it is better than 0.5 meter in resolution Spaceborne slip beam bunching mode in, it is must that point target azimuth resolution that these four factor causes declines The design content that need consider, otherwise, the resolution of point target image, peak sidelobe ratio and integration Secondary lobe does not reaches index request than design load.General document is not the most mentioned these four factor pair Orientation is to the combined influence of imaging performance.But these four factor is respectively to satellite-borne SAR image performance Analysis process and the influence degree size of impact are not belonging to present disclosure.
Step S6: judge azimuth resolution ρaWhether meet less than index request, if orientation To resolution ρaMeet less than resolution index request, show that orientation set in advance is to imaging belt length Degree WsceneLess than normal, then perform step S1 and reset orientation to imaging belt length data WsceneValue, Then step S2~S5 are repeated, until ρaMore than resolution index request, then orientation is to imaging belt Length reaches maximum.Loop iteration method is utilized to meet orientation to resolution when determining target oblique distance for R The maximum imaging belt length that rate requires.
Need to provide following 2 explanations at this:
(1) oblique distance be the point target of R in orientation when the difference of position resolution be slightly changed, ginseng Number design result to ensure that worst resolution meets index request.
(2) the calculation method of parameters derivation of the present invention there is a precondition, it is simply that false During being located at imaging, satellite flight track is approximately straight line.If imaging time is long, satellite flies Row track can not be approximately straight line, then the parameter utilizing the method for the invention to obtain has certain journey The deviation of degree.The method of the invention can ensure to set in the case of meeting three below condition at the same time The accuracy of meter parameter: first, orientation is not less than 0.4 ° to beam angle;Second, orientation is to one-tenth Image-tape length is not more than 15 kilometers;Azimuth resolution is more than or equal to 0.3 meter.The most at present and not In coming 10 years from the point of view of the Development Trends of satellite-borne synthetic aperture slip beam bunching mode, in engineering Attainable spaceborne synthetic aperture radar (SAR) system meets three above condition the most simultaneously, therefore, and this Bright content has higher engineer applied and is worth.

Claims (10)

1. satellite-borne synthetic aperture radar slip beam bunching mode orientation is to a parameter acquiring method, and it is special Levy and be to comprise the following steps that
Step S1: set orientation to initial imaging belt length data, orientation to antenna beamwidth, Orientation to initial angle of strabismus data, orientation to terminating angle of strabismus data and target place oblique distance;
Step S2: utilize orientation to initial imaging belt length data, orientation to antenna beamwidth Data, orientation build to termination angle of strabismus data and target oblique distance to initial angle of strabismus data, orientation The reference oblique distance model that wave beam center of rotation during point target oblique distance is corresponding, it is thus achieved that during point target oblique distance Reference oblique distance data corresponding to wave beam center of rotation;
Step S3: according to the reference oblique distance data that wave beam center of rotation during point target oblique distance R is corresponding Build effective irradiation time model that a point target of oblique distance R is corresponding, it is thus achieved that this point target Effectively initial radiation time data and this point target effectively terminate radiation time data;
Step S4: building oblique distance according to the most initial described irradiation time and termination irradiation time is R Effective Doppler Frequency Model corresponding to a point target and doppler bandwidth model, obtain this Initial effective doppler frequency data corresponding to point target, terminate effective doppler frequency data and many General Le band data;
Step S5: according to initial effective doppler frequency data and the effective Doppler frequency number of termination According to, the resolution descent coefficient that causes of doppler bandwidth data, antenna radiation pattern, imaging processing add Resolution descent coefficient that the resolution descent coefficient that causes of power, Doppler's parameter estimate error cause, It is corresponding that the resolution descent coefficient that imaging algorithm approximation causes builds the point target that oblique distance is R Azimuth resolution model, it is thus achieved that the azimuth resolution that this point target is corresponding;
Step S6: judge whether azimuth resolution meets less than resolution index request, if Show that orientation set in advance is less than normal to imaging belt length data, then perform step S1 and reset Orientation, to imaging belt length data value, repeats step S1~S5 until azimuth resolution is more than dividing Resolution index request, then the maximum orientation of output is to imaging belt length data.
2. satellite-borne synthetic aperture radar slip beam bunching mode orientation obtains to parameter as claimed in claim 1 Access method, it is characterised in that described point target oblique distance is the ginseng that wave beam center of rotation during R is corresponding Examine oblique distance model RrotIt is expressed as follows:
R r o t = R · [ t a n ( ψ e n d + θ a / 2 ) - t a n ( ψ s t a r t - θ a / 2 ) ] - W s c e n e [ t a n ( ψ e n d ) - t a n ( ψ s t a r t ) ] - - - ( 1 )
Wherein: ψendFor orientation to terminating angle of strabismus data, θaFor orientation to antenna beam width degrees of data, ψstartFor orientation to initial angle of strabismus data, WsceneFor orientation to initial imaging belt length data.
3. satellite-borne synthetic aperture radar slip beam bunching mode orientation obtains to parameter as claimed in claim 1 Access method, it is characterised in that described oblique distance is effective irradiation time that the point target of R is corresponding Model includes effectively initiateing irradiation time model tA, start, 0With effectively termination irradiation time model tA, end, 0, Solve tA, start, 0And tA, end, 0Equation group be expressed as follows respectively:
V g · t a , s t a r t , 0 R r o t = t a n ( ψ s t a r t , 0 ) V g · t a , s t a r t , 0 - X p R r o t = t a n ( ψ s t a r t , 0 - θ a / 2 ) - - - ( 2 )
V g · t a , e n d , 0 R r o t = t a n ( ψ e n d , 0 ) V g · t a , e n d , 0 - X p R r o t = t a n ( ψ e n d , 0 + θ a / 2 ) - - - ( 3 )
Wherein: VgFor wave beam ground speed, ψStart, 0Irradiation time is initiateed corresponding for this point target effective Orientation to angle of strabismus, XpFor this point target in the location variable of imaging region,θaFor orientation to beam angle, RrotFor center of rotation oblique distance, ψEnd, 0Fixed Justice effectively terminates orientation corresponding to irradiation time to angle of strabismus for this point target;Solve equation (2) to obtain Obtain tA, start, 0Data, solve equation (3) and obtain tA, end, 0Data.
4. satellite-borne synthetic aperture radar slip beam bunching mode orientation obtains to parameter as claimed in claim 1 Access method, it is characterised in that described oblique distance is effective Doppler frequency that the point target of R is corresponding Rate model includes initial effective Doppler frequency that this point target is corresponding with doppler bandwidth model Model fA, start, 0, terminate effective Doppler frequency fA, end, 0Model, wherein:
f a , s t a r t , 0 = - 2 · V s λ · s i n ( ψ s t a r t , 0 - θ a / 2 ) = - 2 · V s λ · V g · t a , s t a r t , 0 - X p R 2 + ( V g · t a , s t a r t , 0 - X p ) 2 - - - ( 4 )
f a , e n d , 0 = - 2 · V s λ · s i n ( ψ e n d , 0 + θ a / 2 ) = - 2 · V s λ · V g · t a , e n d , 0 - X p R 2 + ( V g · t a , e n d , 0 - X p ) 2 - - - ( 5 )
Wherein: VsIt is that synthetic aperture radar launches the wavelength of signal, ψ for satellite velocities, λStart, 0For oblique distance For the point target of R effectively initiate orientation corresponding to irradiation time to angle of strabismus, θaFor orientation to Beam angle, VgFor wave beam ground speed, tA, start, 0For the most initial irradiation that this point target is corresponding Time, XpFor this point target in the location variable of imaging region,R is Point target oblique distance, ψEnd, 0It is defined as this point target effective and terminates orientation corresponding to irradiation time to tiltedly Visual angle, tA, end, 0For this point target corresponding effectively terminate irradiation time.
5. satellite-borne synthetic aperture radar slip beam bunching mode orientation obtains to parameter as claimed in claim 4 Access method, it is characterised in that described oblique distance is the doppler bandwidth mould that the point target of R is corresponding Type is expressed as follows:
f d = f a , s t a r t , 0 - f a , e n d , 0 = - 2 · V s λ [ V g · t a , s t a r t , 0 - X p R 2 + ( V g · t a , s t a r t , 0 - X p ) 2 - V g · t a , e n d , 0 - X p R 2 + ( V g · t a , e n d , 0 - X p ) 2 ] - - - ( 6 )
Wherein: VsIt is that synthetic aperture radar launches the wavelength of signal, V for satellite velocities, λgFor wave beam ground Face velocity, tA, start, 0For the most initial irradiation time, X that the point target that oblique distance is R is correspondingpFor This point target in the location variable of imaging region,R is this point mesh Mark oblique distance, tA, end, 0For this point target corresponding effectively terminate irradiation time.
6. satellite-borne synthetic aperture radar slip beam bunching mode orientation obtains to parameter as claimed in claim 1 Access method, it is characterised in that described azimuth resolution model ρaIt is expressed as follows:
ρ a = k 1 · k 2 · k 3 · k 4 · V g f d - - - ( 7 )
Wherein: k1The resolution descent coefficient that causes for antenna radiation pattern, k2Cause for imaging processing weighting Resolution descent coefficient, k3The resolution descent coefficient that causes for Doppler's parameter estimate error, k4 Resolution descent coefficient, the f caused is approximated for imaging algorithmdFor the point target pair that oblique distance is R The doppler bandwidth answered, VgFor wave beam ground speed.
7. satellite-borne synthetic aperture radar slip beam bunching mode orientation obtains to parameter as claimed in claim 1 Access method, it is characterised in that simplify geometry mould according to satellite-borne synthetic aperture radar slip beam bunching mode Type obtains the azimuth resolution of any point target that orientation oblique distance in imaging belt is R.
8. satellite-borne synthetic aperture radar slip beam bunching mode orientation obtains to parameter as claimed in claim 1 Access method, it is characterised in that utilize loop iteration method to meet orientation when determining target oblique distance for R To the maximum imaging belt length of resolution requirement.
9. satellite-borne synthetic aperture radar slip beam bunching mode orientation as claimed in claim 1 is to parameter Acquisition methods, it is characterised in that in setting a period of time, satellite flight track is approximately straight line, point Target is static target;During imaging, control azimuth is to beam position, makes beam center point to all the time Away from a virtual point below imaging region, described virtual point is wave beam center of rotation, virtual point It is in below ground;Wave beam the most slowly slips over imaging region, and in imaging region, oblique distance is R Point target go through 3dB beam angle, enter into left side edge from directional diagram right side edge and leave, Initial angle of strabismus during point target approach axis figure right side edge, point target departure direction figure left side Termination angle of strabismus during edge.
10. satellite-borne synthetic aperture radar slip beam bunching mode orientation as claimed in claim 1 is to parameter Acquisition methods, it is characterised in that described orientation is looked side ways to termination to initial angle of strabismus data and orientation The size of angular data and positive and negative not restriction, therefore, it is adaptable to positive side-looking satellite-borne SAR slides poly- Bundle pattern, is also applied for looking side ways satellite-borne SAR slip beam bunching mode.
CN201410041182.7A 2014-01-28 2014-01-28 A kind of satellite-borne synthetic aperture radar slip beam bunching mode orientation is to parameter acquiring method Active CN103792536B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410041182.7A CN103792536B (en) 2014-01-28 2014-01-28 A kind of satellite-borne synthetic aperture radar slip beam bunching mode orientation is to parameter acquiring method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410041182.7A CN103792536B (en) 2014-01-28 2014-01-28 A kind of satellite-borne synthetic aperture radar slip beam bunching mode orientation is to parameter acquiring method

Publications (2)

Publication Number Publication Date
CN103792536A CN103792536A (en) 2014-05-14
CN103792536B true CN103792536B (en) 2016-08-17

Family

ID=50668396

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410041182.7A Active CN103792536B (en) 2014-01-28 2014-01-28 A kind of satellite-borne synthetic aperture radar slip beam bunching mode orientation is to parameter acquiring method

Country Status (1)

Country Link
CN (1) CN103792536B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104880699B (en) * 2015-06-17 2017-10-13 西安电子科技大学 A kind of determination method of the pulse recurrence frequency lower limit of synthetic aperture radar
CN106226768B (en) * 2016-08-09 2018-08-21 北京空间飞行器总体设计部 Ultrahigh resolution agility SAR satellites slide beam bunching mode System Parameter Design method
CN113702971B (en) * 2021-07-14 2024-02-20 西安电子科技大学 Radar beam width design method for synthetic aperture passive positioning system
CN114114266B (en) * 2022-01-24 2022-05-17 北京宏锐星通科技有限公司 Detection method and device for synthetic aperture radar
CN115792907B (en) * 2023-02-07 2023-04-14 中国科学院空天信息创新研究院 Method for designing azimuth imaging parameters of spaceborne SAR squint sliding bunching mode
CN116699613B (en) * 2023-08-04 2023-10-13 中国科学院空天信息创新研究院 On-orbit accurate calculation method for scanning mode imaging starting time

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102176016A (en) * 2011-01-25 2011-09-07 北京航空航天大学 Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001141821A (en) * 1999-11-11 2001-05-25 Mitsubishi Electric Corp Radar signal processor

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102176016A (en) * 2011-01-25 2011-09-07 北京航空航天大学 Large squint sliding spotlight SAR (synthetic aperture radar) imaging processing method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
星载滑动聚束SAR成像模型误差校正方法研究;韩冰 等;《电子与信息学报》;20110731;第33卷(第7期);第1694-1699页 *
条带模式、聚束模式和滑动聚束模式的比较;刘寒艳 等;《中国科学院研究生学报》;20110531;第28卷(第3期);第410-417页 *

Also Published As

Publication number Publication date
CN103792536A (en) 2014-05-14

Similar Documents

Publication Publication Date Title
CN103792536B (en) A kind of satellite-borne synthetic aperture radar slip beam bunching mode orientation is to parameter acquiring method
CN101858976B (en) Moving target detection method based on multiple sub-apertures of single-channel SAR
CN104076353B (en) A kind of Area Objects echo beam center speed measurement method
CN102346249B (en) Implementation method for wide swath earth observation step scanning mode of synthetic aperture radar
CN104035095B (en) Based on the low level wind shear velocity estimation method of optimal processor during sky
CN102508243B (en) Beam position design method of inclined geosynchronous orbit synthetic aperture radar
CN104076343B (en) Satellite-borne three-channel SAR-GMTI self-adaptive clutter suppression method
CN102288964A (en) Imaging processing method for spaceborne high-resolution synthetic aperture radar
CN103364782B (en) Geosynchronous orbit synthetic aperture radar velocity spatial variability compensating method
CN100561250C (en) The clutter base-band analog signal production method of spaceborne bistatic radar
CN103675775B (en) Background ionosphere is to the analytical approach of GEO SAR Imaging
CN103076607B (en) Method for realizing sliding spotlight mode based on SAR (Synthetic Aperture Radar) satellite attitude control
CN107607936B (en) High-frequency sky-ground wave radar ocean surface flow inversion method
CN103235309B (en) Near space low-speed platform SAR (Synthetic Aperture Radar) imaging method
CN109444882B (en) Double-station SAR imaging method based on variable squint elliptical beam synchronous model
CN106556822A (en) Spaceborne Sliding spotlight SAR pointing accuracy Orbital detection method
CN102169173A (en) Method for analyzing ambiguity of inclined geo-synchronization orbit synthetic aperture radar
CN106093932A (en) A kind of high-resolution radar scatterometer of scanning beam
CN104076338A (en) Airborne radar clutter simulation method based on digital elevation and digital ground coverage
CN107976674A (en) The PRF systems of selection of dualbeam SAR obscured based on top optimization direction a kind of and system
CN105652255A (en) Spatial aligning method for radar networking system
CN102736073B (en) Method for computing range ambiguity of satellite-borne synthetic aperture radar (SAR) in universal mode
CN110018474A (en) Three-D imaging method based on geostationary orbit synthetic aperture radar chromatographic technique
CN108267736A (en) A kind of GEO SAR staring imagings mode orientation fuzziness determines method
CN105738887A (en) Airborne radar clutter power spectrum optimization method based on Doppler channel division

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