CN103675760B - A kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method - Google Patents
A kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method Download PDFInfo
- Publication number
- CN103675760B CN103675760B CN201310651459.3A CN201310651459A CN103675760B CN 103675760 B CN103675760 B CN 103675760B CN 201310651459 A CN201310651459 A CN 201310651459A CN 103675760 B CN103675760 B CN 103675760B
- Authority
- CN
- China
- Prior art keywords
- angle
- omega
- guiding
- attitude
- resolution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9058—Bistatic or multistatic SAR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
Abstract
The invention discloses a kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method.Use the present invention that less attitude guiding angle can be utilized to realize optimum ground range resolution.First the present invention differentiates the expression formula of ellipse area according to the condition acquisition satellite velocities directions such as known orbit parameter, downwards angle of visibility and distance, when then utilizing distance resolution ellipse area minimum, ground range resolution is optimum, calculate the top optimization direction angle of optimum ground range resolution, direction, place, beam position top optimization direction angle is made finally by pitching roll guiding or roll pitching guiding, complete the guiding of optimal resolution attitude, attitude steering angle is little.
Description
Technical field
The present invention relates to Synthetic Aperture Radar Technique field, be specifically related to a kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method.
Background technology
Synthetic aperture radar (SAR) is a kind of round-the-clock, round-the-clock high-resolution microwave remote sensing imaging radar, can be arranged on the flying platforms such as aircraft, satellite, spaceship.Application in environmental monitoring, oceanographic observation, resource exploration, Crop Estimation, mapping and military affairs etc. has unique advantage, the effect that other remote sensings are difficult to play can be played.Apply so widely for SAR, its application quality quality generally depends on the image resolution ratio of SAR.
Spaceborne geostationary orbit (GEO) synthetic-aperture radar refer to be positioned at geostationary orbit, take satellite as the synthetic-aperture radar of carrier.Spaceborne geostationary orbit synthetic-aperture radar (GEOSAR) has the advantages such as overlay area is wide, revisit time is short.In the orbital period of a GEOSAR, the direction variation range that earth center rotates (ECR) coordinate system speed is large.When GEOSAR is positioned at equator, ECR velocity reversal is close to range direction.Now, ambiguity function is projected in the two-dimentional secondary lobe on ground will be no longer orthogonal.This causes orientation to ground range resolution much larger than orientation to slant range resolution, ground range resolution severe exacerbation.
Synthetic-aperture radar can improve range resolution by improving transmitted bandwidth, can improve azimuth resolution by extending the synthetic aperture time.But because distance is no longer orthogonal to bidimensional to, orientation, even may overlap, thus only by raising transmitted bandwidth with extend the synthetic aperture time and be difficult to effectively improve resolution.
Synthetic-aperture radar can also change beam position by attitude guiding, and then improves ground range resolution.Attitude guiding refers to the adjustment for reaching the attitude of satellite that certain ideal behavior is carried out.In earth low orbit (LEO) synthetic-aperture radar, the object of attitude guiding normally reduces Doppler frequency.The attitude guidance method that current LEOSAR is conventional has one dimension Yaw steering and complete zero Doppler to guide, and the former can reduce Doppler frequency, and the latter can make Doppler frequency be zero.In GEOSAR, Yaw steering effectively can reduce the deterioration of ground range resolution, but required guiding angle is large, high to the requirement of platform stance control system.Complete zero Doppler's guiding not only needs larger guiding angle, and can not reduce the deterioration of ground range resolution.
Summary of the invention
In view of this, the invention provides a kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method, less attitude guiding angle can be utilized to realize optimum ground range resolution.
Spaceborne geostationary orbit synthetic-aperture radar attitude guidance method of the present invention, comprises the following steps:
Step 1, definitely apart from differentiating ellipse area:
Wherein, S is the area of-3dB ground range resolution ellipse, ρ
afor oblique distance azimuth resolution; ρ
rfor oblique distance range resolution; α is the angle of pitch of radar; β is the downwards angle of visibility of radar;
for the position angle of radar;
Wherein,
Wherein, ω
sfor orbit angular velocity, ω
efor rotational-angular velocity of the earth, R
sfor the height in satellite distance the earth's core,
for R
sto the derivative of time, i is orbit inclination, and u is latitude argument, and μ is Gravitational coefficient of the Earth, and a is semi-major axis of orbit, and e is orbital eccentricity, and ω is argument of perigee;
Step 2, base area is apart from resolution ellipse area determination top optimization direction angle:
Solve
draw top optimization direction angle
for
Wherein, positive sign represents that right side is looked, and negative sign represents that left side is looked;
Step 3, determines attitude guiding angle according to top optimization direction angle, carries out attitude guiding:
According to the guidance mode of pitching after first roll, then roll guiding angle θ
r, pitching guiding angle θ
pbe respectively
According to the guidance mode of roll after first pitching, then roll guiding angle θ
r, pitching guiding angle θ
pbe respectively
Wherein,
Beneficial effect:
When the present invention utilizes distance resolution ellipse area minimum, ground range resolution is optimum, calculate the top optimization direction angle of optimum ground range resolution, then guided by pitching roll or roll pitching guiding make direction, place, beam position top optimization direction angle, complete the guiding of optimal resolution attitude, attitude steering angle is little.
Accompanying drawing explanation
Fig. 1 is process flow diagram of the present invention.
Fig. 2 is coordinate angle schematic diagram of the present invention.
Fig. 3 is that the attitude guiding roll angle of pitching guidance mode after adopting first roll changes.
Fig. 4 is that the attitude guiding angle of pitch of pitching guidance mode after adopting first roll changes.
Fig. 5 is distance ambiguity function after attitude guiding.
Embodiment
To develop simultaneously embodiment below in conjunction with accompanying drawing, describe the present invention.
The invention provides a kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method, first differentiate the expression formula of ellipse area according to the condition acquisition satellite velocities directions such as known orbit parameter, downwards angle of visibility and distance.Then base area is apart from the expression formula differentiating ellipse area, solves the top optimization direction angle realized needed for minimum resolution ellipse area.Finally, guided by pitching roll or roll pitching guiding make direction, place, beam position above-mentioned top optimization direction angle, as shown in Figure 1.Concrete steps are as follows.
Step 1, differentiates ellipse area according to orbit parameter and downwards angle of visibility determination satellite velocities direction and distance.
Step 1.1, determines satellite velocities direction.
The velocity reversal of satellite in earth fixed coordinate system is determined according to semi-major axis of orbit a, orbit inclination i, orbital eccentricity e, argument of perigee ω, latitude argument u.Wherein, earth fixed coordinate system initial point is positioned at earth centroid, and x-axis points to Greenwich meridian in earth equatorial plane, and z-axis is along the sensing arctic, earth rotation direction, and y-axis and x-axis, z-axis form right hand rectangular coordinate system.
Under earth fixed coordinate system, the direction of satellite velocities can use crab angle φ and angle of pitch α to represent, crab angle φ and angle of pitch α computing formula as follows:
Wherein, ω
efor rotational-angular velocity of the earth, ω
sfor orbit angular velocity, R
sfor the height in satellite distance the earth's core,
for R
sto the derivative of time, can be expressed as
Wherein, μ is Gravitational coefficient of the Earth, and numerical value is 398600.5km
3/ s
2; Sign () is sign function, can be expressed as
Step 1.2, definitely apart from differentiating ellipse area
Ground range resolution can use the cartographic represenation of area of-3dB resolution ellipse, and its formula is
Wherein, S is the area of-3dB resolution ellipse, ρ
afor oblique distance azimuth resolution, ρ
rfor oblique distance range resolution.Angle of squint γ is the angle of beam position and ECR velocity reversal, position angle
for beam position is at the projection of surface level and the angle of ECR speed between the projection of surface level, β is downwards angle of visibility, position angle
determine according to the beam position of reality with downwards angle of visibility β.
As shown in Figure 2, wherein V is ECR speed to above-mentioned angle, and R is beam position, and O point is the earth's core.
Angle of squint γ and position angle
between there is following relation:
Thus, ground range resolution area can be expressed as
Step 2, base area is apart from resolution ellipse area determination top optimization direction angle.
If ground range resolution area is minimum, demand fulfillment equation
Solving equation (7), can obtain the top optimization direction angle corresponding to optimum ground range resolution:
Wherein,
be top optimization direction angle, positive sign represents that right side is looked, and negative sign represents that left side is looked.
Step 3, determines attitude guiding angle according to top optimization direction angle.
In initial time beam position substar, position, beam position top optimization direction angle after attitude guiding.Attitude guiding strategy is the guiding of roll pitching two dimension.Beam position before attitude guiding can represent
u'=[0,0,-1]
T(9)
Beam position after attitude guiding can be expressed as
Attitude guided procedure can be expressed as
u'=Au(11)
Wherein, A is the rotation matrix of attitude guiding correspondence.
According to the strategy of pitching after first roll, rotation matrix can be expressed as
Wherein, θ
rfor roll guiding angle, θ
pfor pitching guiding angle.Rotation matrix is substituted into equation (11), solving equation (11) can obtain attitude steering angle and be
According to the strategy of roll after first pitching, rotation matrix can be expressed as
Rotation matrix substitutes into equation (11), and solving equation (11) can obtain attitude steering angle and be
Since then, a kind of satellite-borne synthetic aperture radar attitude guidance method is just achieved.
Satellite transit on elliptical orbit, shown in each parameter of track is specific as follows: earth radius is 6371.004km; Semi-major axis of orbit is 42100km; Orbit inclination is 50 °; Orbital eccentricity is 0.1; Perigee of orbit argument is 90 °; Downwards angle of visibility is 7 °; Slant range resolution is 1m; Rotational-angular velocity of the earth is 7.292115 × 10
-5rads.The mode of pitching after the first roll of attitude guiding employing.As shown in Figure 3, the angle of pitch as shown in Figure 4 for roll angle needed for attitude guiding.Ground range resolution ambiguity function after attitude guides as shown in Figure 5.From Fig. 3, Fig. 4, the attitude steering angle needed for geostationary orbit synthetic-aperture radar optimal resolution attitude guidance method is less.As shown in Figure 5, after carrying out the guiding of geostationary orbit synthetic-aperture radar optimal resolution attitude, it is vertical that distance differentiates two secondary lobes, and ground range resolution is optimum.
In sum, these are only preferred embodiment of the present invention, be not intended to limit protection scope of the present invention.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (1)
1. a spaceborne geostationary orbit synthetic-aperture radar attitude guidance method, is characterized in that, comprise the following steps:
Step 1, definitely apart from differentiating ellipse area:
Wherein, S is the area of-3dB ground range resolution ellipse, ρ
afor oblique distance azimuth resolution; ρ
rfor oblique distance range resolution; α is the angle of pitch of radar; β is the downwards angle of visibility of radar;
for the position angle of radar;
Wherein, sign () is sign function, can be expressed as
Wherein, ω
sfor orbit angular velocity, ω
efor rotational-angular velocity of the earth, R
sfor the height in satellite distance the earth's core,
for R
sto the derivative of time, i is orbit inclination, and u is latitude argument, and μ is Gravitational coefficient of the Earth, and a is semi-major axis of orbit, and e is orbital eccentricity, and ω is argument of perigee;
Step 2, base area is apart from resolution ellipse area determination top optimization direction angle:
Solve
draw top optimization direction angle
for
Wherein, positive sign represents that right side is looked, and negative sign represents that left side is looked;
Step 3, determines attitude guiding angle according to top optimization direction angle, carries out attitude guiding:
According to the guidance mode of pitching after first roll, then roll guiding angle θ
r, pitching guiding angle θ
pbe respectively
According to the guidance mode of roll after first pitching, then roll guiding angle θ
r, pitching guiding angle θ
pbe respectively
Wherein,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310651459.3A CN103675760B (en) | 2013-12-03 | 2013-12-03 | A kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310651459.3A CN103675760B (en) | 2013-12-03 | 2013-12-03 | A kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103675760A CN103675760A (en) | 2014-03-26 |
CN103675760B true CN103675760B (en) | 2015-12-02 |
Family
ID=50313903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310651459.3A Expired - Fee Related CN103675760B (en) | 2013-12-03 | 2013-12-03 | A kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103675760B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001762A (en) * | 2018-06-08 | 2018-12-14 | 上海微小卫星工程中心 | A kind of the gas light suppressing method and system of synchronous belt objective optics observation satellite |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104597446B (en) * | 2014-10-10 | 2017-02-15 | 北京理工大学 | Space-borne synthetic aperture radar ground range resolution representation and parameter design method |
CN104375511A (en) * | 2014-10-28 | 2015-02-25 | 上海卫星工程研究所 | Geosynchronous orbit SAR satellite off-course guide method based on wave beam cooperative control |
CN104765023B (en) * | 2015-03-20 | 2017-03-29 | 北京理工大学 | A kind of satellite-borne synthetic aperture radar height resolution computational methods |
CN105403888A (en) * | 2015-12-30 | 2016-03-16 | 中国科学院电子学研究所 | Geosynchronous orbit circular track SAR imaging method using beam pointing control |
CN108958272A (en) * | 2018-06-15 | 2018-12-07 | 上海卫星工程研究所 | Yaw steering method is imaged in fixed statellite |
CN111796277B (en) * | 2020-02-28 | 2022-05-20 | 北京理工大学 | Through-wall radar rapid imaging method based on unmanned aerial vehicle platform |
CN112255606A (en) * | 2020-09-29 | 2021-01-22 | 北京空间飞行器总体设计部 | Method for calculating front side-view imaging attitude angle of Geo-SAR (synthetic aperture radar) satellite based on single reflector antenna |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3114600C2 (en) * | 1981-04-10 | 1990-06-21 | Telefunken Systemtechnik Gmbh, 7900 Ulm, De | |
EP1965223A1 (en) * | 2007-03-02 | 2008-09-03 | Saab Ab | Subsurface Imaging radar |
CN101464511B (en) * | 2009-01-16 | 2011-07-27 | 北京航空航天大学 | Working wave position decision method for satellite-loaded synthetic aperture radar |
CN102608576A (en) * | 2012-03-20 | 2012-07-25 | 北京理工大学 | Geometric correction method for large rake forward synthetic aperture radar return image |
CN103364782A (en) * | 2013-07-04 | 2013-10-23 | 北京理工大学 | Geosynchronous orbit synthetic aperture radar velocity spatial variability compensating method |
-
2013
- 2013-12-03 CN CN201310651459.3A patent/CN103675760B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3114600C2 (en) * | 1981-04-10 | 1990-06-21 | Telefunken Systemtechnik Gmbh, 7900 Ulm, De | |
EP1965223A1 (en) * | 2007-03-02 | 2008-09-03 | Saab Ab | Subsurface Imaging radar |
CN101256235A (en) * | 2007-03-02 | 2008-09-03 | Saab公司 | Subsurface imaging radar |
CN101464511B (en) * | 2009-01-16 | 2011-07-27 | 北京航空航天大学 | Working wave position decision method for satellite-loaded synthetic aperture radar |
CN102608576A (en) * | 2012-03-20 | 2012-07-25 | 北京理工大学 | Geometric correction method for large rake forward synthetic aperture radar return image |
CN103364782A (en) * | 2013-07-04 | 2013-10-23 | 北京理工大学 | Geosynchronous orbit synthetic aperture radar velocity spatial variability compensating method |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001762A (en) * | 2018-06-08 | 2018-12-14 | 上海微小卫星工程中心 | A kind of the gas light suppressing method and system of synchronous belt objective optics observation satellite |
Also Published As
Publication number | Publication date |
---|---|
CN103675760A (en) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103675760B (en) | A kind of spaceborne geostationary orbit synthetic-aperture radar attitude guidance method | |
CN113439220B (en) | Remote sensing system and constellation system based on satellite formation | |
ES2951990T3 (en) | Device and method of navigation | |
CN101893440B (en) | Celestial autonomous navigation method based on star sensors | |
CN107607947B (en) | On-line estimation method for imaging parameters of satellite-borne radar based on Kalman filtering | |
CN102508243B (en) | Beam position design method of inclined geosynchronous orbit synthetic aperture radar | |
CN103217987B (en) | Agile satellite dynamic imaging posture adjustment method | |
CN101430379B (en) | Synthetic aperture radar three-dimensional microwave imaging method for circular track of earth synchronization orbit | |
CN102169173B (en) | Method for analyzing ambiguity of inclined geo-synchronization orbit synthetic aperture radar | |
CN106871932A (en) | The in-orbit sensing calibration method of satellite borne laser based on Pyramidal search terrain match | |
CN101750067B (en) | Imaging type method for correcting earth oblateness by earth sensor | |
CN102323571B (en) | Distribution method of satellite-borne dual-antenna SAR (Synthetic Aperture Radar) interferometric calibrator with comprehensive overall parameter | |
CN108710379A (en) | Fixed statellite is imaged Yaw steering angle computational methods | |
CN102621994B (en) | Control method of geosynchronous earth orbit (GEO) synthetic aperture radar (SAR) for covering all over China | |
CN105004354A (en) | Unmanned aerial vehicle visible light and infrared image target positioning method under large squint angle | |
CN105242285A (en) | Method based on satellite communication for identifying deception jamming acting on navigation data of unmanned plane | |
CN105403888A (en) | Geosynchronous orbit circular track SAR imaging method using beam pointing control | |
CN101513939A (en) | Two dimentional attitude control system of synthetic aperture radar satellite | |
CN103344958B (en) | Based on the satellite-borne SAR high-order Doppler parameter evaluation method of almanac data | |
CN112255606A (en) | Method for calculating front side-view imaging attitude angle of Geo-SAR (synthetic aperture radar) satellite based on single reflector antenna | |
CN103487032A (en) | Low earth orbit space camera free-pointing image motion vector calculation method | |
CN103940429A (en) | Real-time measuring method of carrier attitude in transverse coordinate system of inertial navigation system | |
US20030006927A1 (en) | Satellite configuration for interferometric and/or tomographic remote sensing by means of synthetic aperture radar (SAR) | |
CN104597446B (en) | Space-borne synthetic aperture radar ground range resolution representation and parameter design method | |
CN103245948B (en) | Image match navigation method for double-area image formation synthetic aperture radars |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |
Granted publication date: 20151202 Termination date: 20161203 |
|
CF01 | Termination of patent right due to non-payment of annual fee |