CN110596704A - Satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation - Google Patents

Satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation Download PDF

Info

Publication number
CN110596704A
CN110596704A CN201910765450.2A CN201910765450A CN110596704A CN 110596704 A CN110596704 A CN 110596704A CN 201910765450 A CN201910765450 A CN 201910765450A CN 110596704 A CN110596704 A CN 110596704A
Authority
CN
China
Prior art keywords
imaging
satellite
center
determining
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910765450.2A
Other languages
Chinese (zh)
Other versions
CN110596704B (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.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
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 Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201910765450.2A priority Critical patent/CN110596704B/en
Publication of CN110596704A publication Critical patent/CN110596704A/en
Application granted granted Critical
Publication of CN110596704B publication Critical patent/CN110596704B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q1/00—Details of, or arrangements associated with, antennas
    • H01Q1/27—Adaptation for use in or on movable bodies
    • H01Q1/28—Adaptation for use in or on aircraft, missiles, satellites, or balloons
    • H01Q1/288—Satellite antennas
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q3/00—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system
    • H01Q3/02—Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system using mechanical movement of antenna or antenna system as a whole

Landscapes

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

Abstract

The invention relates to a satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation, comprising the steps of s1, determining a satellite platform attitude maneuver mode of primary imaging, namely determining the position of a satellite platform at the initial imaging moment, the satellite attitude and a virtual rotation center in the primary imaging process; s2, performing sliding bunching imaging, and determining the satellite position at the next imaging starting moment; s3, determining satellite attitudes in the maneuvering process between two adjacent imaging and a virtual rotation center in the next imaging process according to the satellite position of the next imaging starting moment and the distance between the center of the beam after maneuvering and the scene, and adjusting the beam direction to the azimuth direction starting position of the scene in the attitude maneuvering adjustment process; s4, repeating the steps s2 and s3 until the imaging work is finished; in the sliding beamforming imaging process, the beam center of the satellite platform always points to the virtual rotation center determined in the imaging process.

Description

Satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation
Technical Field
The invention belongs to the technical field of space microwave remote sensing, and particularly relates to a satellite platform attitude maneuver method.
Background
In order to improve the reconnaissance efficiency of the satellite-borne SAR, the multi-azimuth repeated observation capability of the satellite-borne SAR on the military sensitive area in a single transit process is required. If the phased array antenna system is used for realizing the purpose, the phased array antenna system realizes beam scanning and repeated observation by controlling the phase shifting codes of the receiving and transmitting channels, and the method has the advantages that the satellite platform does not need to be maneuvered, but has the defect that the gain of antenna beams is reduced along with the increase of the rotation angle. This will result in a loss of the echo signal-to-noise ratio and a drastic degradation of the SAR image quality. If the method is realized by using the reflector antenna, the method needs to be realized by maneuvering of the satellite platform, and the method has the advantages that the peak antenna gain can be always pointed to an observation area and cannot be changed along with the change of the rotation angle, but the attitude maneuvering mode of the satellite platform needs to be designed to ensure that the SAR load stably works in orbit. There is no suitable solution to the above problems.
Disclosure of Invention
The technical problem solved by the invention is as follows: the method overcomes the defects of the prior art and provides a satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation.
The technical scheme of the invention is as follows: a satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation comprises the following steps:
s1, determining the attitude maneuver mode of the satellite platform for the first imaging, namely determining the position of the satellite platform at the initial imaging moment, the satellite attitude and the virtual rotation center in the first imaging process;
s2, performing sliding bunching imaging, and determining the satellite position at the next imaging starting moment;
s3, determining satellite attitudes in the maneuvering process between two adjacent imaging and a virtual rotation center in the next imaging process according to the satellite position of the next imaging starting moment and the distance between the center of the beam after maneuvering and the scene, and adjusting the beam direction to the azimuth direction starting position of the scene in the attitude maneuvering adjustment process;
s4, repeating the steps s2 and s3 until the imaging work is finished;
in the sliding beamforming imaging process, the beam center of the satellite platform always points to the virtual rotation center determined in the imaging process.
Preferably, the virtual center of rotation during the first imaging in step s1 is determined by:
initial squint angle theta by sliding bunch imaging1sSatellite platform flight speed VsAnd beam footprint velocity V in stripe modegAnd the intersection point of the antenna beam center and the ground at the initial imaging moment of the satelliteDetermining satellite imaging start time position
Starting time position by satellite imagingIntersection of antenna beam center and groundSatellite Z-axis pointing to determine starting imaging time
From satellitesImaging start time positionSatellite Z-axis pointing at start timeAnd determining a virtual rotation center position by using the sliding bunching resolution improvement factor A
Preferably, the satellite imaging start time positionThe determination is made by:
determining the satellite imaging initial time, the intersection point of the antenna beam center and the groundAnd scene centerThe distance in the azimuth direction is further determined, and the intersection point of the antenna beam center and the ground is further determinedA location;
according to the intersection point of the antenna beam center and the groundPosition, determining the intersection pointSatellite position at zero doppler time
Determining connectionsAndvector and join of two pointsAndangle theta between vectors of two pointsmax；
According to the angle thetamaxDetermining the imaging start time position
Above, LaIs the azimuth length of the scene, RcPositive side-view skew distance, theta, at zero Doppler time from the center of the sceneazThe azimuth beamwidth of the antenna.
Preferably, the included angle θmaxAccording to the relationAnd (4) determining.
Preferably, the virtual center of rotation positionThe following relationship is satisfied:
to be started from a satellitePointing to first sliding spotlight imaging virtual center of rotationThe vector of (2).
Preferably, step s2 is implemented by:
determining the satellite position of the imaging ending time according to the position of the imaging starting time of the satellite, the imaging time length and the flight orbit;
determining the satellite attitude at the imaging termination time and the intersection point of the satellite attitude and the ground according to the determined satellite position at the imaging termination time and the virtual rotation center of the imaging process;
and determining the satellite position at the moment of starting the next imaging after maneuvering according to the determined satellite position at the moment of ending the imaging and given attitude maneuvering time between two times of imaging.
Preferably, step s3 is implemented by:
from the satellite position at the end of this imagingSatellite position at the start of next imagingAnd the intersection point of the beam center and the ground at the imaging termination momentDetermining the plane determined by the three points and the distance D between the center of the beam and the scene after maneuvering to determine the intersection point between the center of the beam and the ground after the attitude maneuvering is adjusted
Satellite position adjusted by attitude maneuverAnd intersection of beam center and groundObtaining the Z-axis direction of the satellite at the maneuvering ending moment;
by vectorAndto determine the virtual center of rotation of the satellite platform during the maneuvering between two adjacent images
By virtual centre of rotation of a motorised processDetermining the three-axis attitude of the satellite in the whole maneuvering process;
byAnd determining the virtual rotation center of the next sliding bunching imaging by using the sliding bunching resolution improvement factor A
Preferably, the distance D is calculated as follows:
wherein theta isendThe oblique angle of the beam center at the end of the sliding beam bunching imaging is RcPositive side-view skew distance, theta, at zero Doppler time from the center of the sceneazThe azimuth beamwidth of the antenna.
Preferably, the beam center is located at the groundPoint of intersection of planesIs determined by the following steps:
firstly, a vector parallel to the azimuth edge of the scene is made at a position which is at a distance D from the azimuth front edge of the scene
Then, the vector is determinedAnd plane surfaceThe intersection point of the antenna beam center pointing to the ground at the moment of termination of the maneuver
As described aboveThe plane is the satellite position of the current imaging end timeSatellite position at the starting time of next imagingAnd the intersection point of the beam center and the ground at the current imaging end momentThe plane determined by the three points.
Preferably, the beam rotation angular velocity during the attitude maneuver adjustment is greater than the beam rotation angular velocity during the sliding beamforming imaging.
Compared with the prior art, the invention has the beneficial effects that:
in order to improve the reconnaissance efficiency of the satellite-borne SAR, the multi-azimuth repeated observation capability of the satellite-borne SAR on the military sensitive area in a single transit process is required. Aiming at the aim, the invention provides a method realized by attitude maneuver of a satellite platform. Through attitude maneuver of the satellite platform, the peak value of an antenna directional diagram can always point to an observation scene, so that the highest echo signal-to-noise ratio and the optimal image quality are obtained. The method is designed based on the system parameters of the load, the longitude and latitude position information of an observation target, and finally realized relevant indexes such as resolution, breadth, the number of times of repeated observation and the like: 1) a satellite attitude stable maneuvering mode in each observation process; 2) and the satellite attitude is quickly adjusted in a maneuvering way between two observation processes so as to ensure that the satellite can stably work in an imaging way.
The method has the following two advantages:
1) in the attitude maneuver process of satellite imaging, the designed three-axis pointing of the satellite ensures that all range gates have the same Doppler frequency in the echo received each time, thereby simplifying the difficulty in later ground imaging processing.
2) In the satellite platform attitude rapid maneuvering process between two adjacent imaging, the Z axis with the fastest maneuvering is always positioned in a two-dimensional plane determined by the satellite position at the current imaging termination time, the satellite position at the next imaging starting time and the current Z axis direction. Therefore, compared with maneuvering in three-dimensional space, the method can minimize the maneuvering amplitude and speed of the satellite, so that two adjacent imaging can be transited more smoothly.
Drawings
FIG. 1 is a schematic diagram of the multi-azimuth observation mode of the present invention;
FIG. 2 is a schematic view of the whole process of two sliding bunching observations and the rapid attitude maneuver adjustment beam pointing between them of the present invention;
FIG. 3 is a schematic diagram of the start time and position of a satellite according to the present invention;
FIG. 4 is a flow chart of the present invention for determining a satellite platform position and attitude at an initial imaging time;
FIG. 5 is a flowchart illustrating the method for determining the virtual point position of the beam pointing direction during the first imaging process and the position of the beam pointing direction at the imaging end time according to the present invention;
FIG. 6 is a schematic diagram of a process of the satellite platform rapidly maneuvering to achieve antenna beam pointing from a scene terminal to a starting end;
FIG. 7 is a diagram of the intersection point of the center of the antenna beam and the ground at two momentsAndis a schematic diagram of the position relation;
FIG. 8 is a flow chart of the design of the fast maneuvering mode of the satellite platform between two adjacent sliding beamforming imaging according to the present invention.
Detailed Description
The invention is further illustrated by the following examples.
The new multi-azimuth observation mode requires the satellite platform to maneuver within a large rotation angle range of the azimuth direction so as to realize the imaging of the multiple sliding bunching mode in the military area. Between two adjacent imaging, the satellite platform needs to be quickly maneuvered, the antenna beam direction is adjusted to the initial position of the scene again, and the next sliding beam-bunching imaging process is started again.
Figure 1 shows the working principle of the multi-azimuth observation mode. In the stage 1 noted in the figure, the satellite-borne SAR starts to perform azimuth super-large squint angle sliding bunching imaging on the scene; after the work of the stage is completed, the satellite platform starts to rapidly maneuver, the beam direction is adjusted to the initial position of the scene again, and the sliding beamforming imaging of the 2 nd stage is started. And analogizing in sequence until the whole satellite finishes all n times of sliding bunching imaging, thereby realizing multi-azimuth repeated observation of the SAR load on the area.
The following will give an example of the first and second imaging observation processes of the mode and the maneuvering process between them, how to design the attitude maneuvering mode of the satellite in the whole working process according to multiple conditions such as imaging resolution, imaging width, antenna beam width, and maneuvering time constraint.
In fig. 2, the whole process of two sliding beamforming observations and the fast attitude maneuver between them to adjust the beam pointing is shown. It can be seen from the figure that the two adjacent sliding bunching imaging processes are attitude rotation around different virtual rotation centers (the first time is O)1The second time is O2) (ii) a In addition, the rotating angular speed of the wave beam in the process of quick maneuvering is obviously greater than that of the wave beam in the process of sliding beam bunching imaging. In attached table 1, the working parameters and constraint conditions that need to be given in advance in the design process of the attitude maneuver are given:
TABLE 1 relevant parameters for Multi-azimuthal repeat Observation
Serial number System/imaging parameters Symbol Remarks for note
1. Scene center side view slope distance Rc System parameter
2. Satellite flight velocity Vs System parameter
3. Ground speed of travel of a beam(strip mode) Vg System parameter
4. System carrier frequency wavelength λ System parameter
5. Antenna azimuth beam width θa System parameter
6. Azimuthal resolution requirement ρa Constraint conditions
7. Azimuth length of scene La Constraint conditions
8. Azimuth beam center squint angle at initial sliding convergence moment θ1s Initial conditions
9. Angular range of azimuthal strabismus observation ±θsq_max Constraint conditions
10. Platform maneuver time between two sliding gather imaging Tm Constraint conditions
Under the condition that the system parameters and the constraint conditions are given, the attitude change mode of the satellite platform in the sliding bunching imaging process and the maneuvering process between the sliding bunching imaging process and the maneuvering process is designed to ensure that continuous multiple times of observation can be realized.
(1) Step 1: satellite platform attitude maneuver mode for determining primary imaging
During the first step of analysis, the following three parameters need to be determined according to the system parameters and constraints listed in attached table 1:
the platform flight distance in the first sliding bunching imaging process is as follows: l iss1；
Beam center squint angle at the first sliding beamforming imaging termination time: theta1e；
Beam rotation angular velocity at the first sliding beamforming termination time: omega1e；
When the load is in the strip mode, its azimuthal resolution may be ρa_stripCan be expressed as
Wherein VsIs the satellite platform flight speed, VgFor space-borne SAR loadingBeam footprint travel speed when operating in stripe mode.
Based on the formula (1) to give the azimuth resolution of the strip mode and the azimuth resolution required by the sliding bunching mode, the improvement factor A of the sliding bunching mode can be obtained as
Then combines with the center slant distance R of the front side view fieldcThe slant distance R of the sliding bunching rotation center can be obtainedtotIs composed of
When the azimuth oblique angle of the antenna beam center is theta1sThe length L of the area covered by the beam in azimuth directionbIs composed of
Therefore, the travel distance L of the beam on the ground in the first sliding beam bunching imaging stage1Is the azimuth length L of the sceneaAnd beam azimuth coverage length LbSum of
L1＝La+Lb (5)
Combining with the sliding bunching resolution improvement factor A, the flight distance L of the satellite platform can be obtaineds1Is composed of
Distance of flight L through satellite platforms1Beam center starting squint angle theta of antenna1sAnd center of rotation slope distance R of sliding bunchingtotThe beam center squint angle theta at the sliding bunching termination moment can be obtained1eIs composed of
At the end time of the first sliding beam bunching imaging, the corresponding beam center rotation angular velocity omega1eIs composed of
After the first imaging is finished, the satellite platform starts to rapidly maneuver, and the beam direction is adjusted to the azimuth initial position of the scene.
In the specific implementation process of step 1, the scene center position based on the scene is firstly determinedScene size La×LrAnd maximum squint angle operation capability of antenna beam [ -theta [ ]1s,θ1s]To determine: 1) position of satellite platform at initial imaging moment2) Attitude of satellite at that time3) Beam center pointing virtual point in primary imaging process
When the satellite is at the starting positionWhen the antenna beam has a maximum squint angle theta1sAnd the beam front edge just meets the scene edge.
Based on the above analysis, the intersection point of the antenna beam center and the ground at the satellite initial position can be obtainedIn azimuth of the scene centerThe distance in the direction is:
wherein L isaIs the azimuth length of the scene, RcPositive side-view skew distance, theta, at zero Doppler time from the center of the sceneazThe azimuth beamwidth of the antenna. In addition, the first and second substrates are,andare positioned in the same distance door.
At the determined intersection pointAfter the position of (2), its corresponding zero doppler time can be determined, as well as the position of the satellite at that timeThe starting position of the satellite can then be determined on the orbit Satisfies the relationship: connection ofAndvector and join of two pointsAndclamp between vectors of two pointsAngle thetamax(as noted in FIG. 3), θmaxIs thetasq_max(i.e. theta)1s) The relationship between is
Wherein VsIs the satellite platform flight speed, and VgThe beam ground travel speed when the load is operating in strip mode. Theta can be obtained by the formula (10)maxAnd the starting position of the satellite
Finally, from the starting position of the satelliteAnd intersection of beam center and groundThe direction of the antenna beam center at the starting moment can be determinedBecause the antenna and the satellite are fixedly connected, the beam center pointing of the antenna is the Z-axis pointing of the satelliteY-axis pointing of satellitePerpendicular toWith satellite velocityForm a slant plane, the X-axis of the satellite pointingAndbecome the right hand rule.
At this time, the initial position of the satellite, the direction of the antenna beam at the initial time and the attitude information of the whole satellite in the multi-azimuth observation mode are obtained. In fig. 4, a flow chart for determining the satellite platform position and attitude at the time of the initial imaging is given.
Next, the coordinates of the beam pointing to the virtual center during the first imaging process need to be determined. During first sliding spotlight imaging, the flying distance L of the platform iss1As shown in equation (6), in combination with the satellite flight velocity VsThe working time T of the first sliding spotlight imaging can be obtainedacq。
WhereinCombining the satellite initial position obtained in the early stageThe position of the stage at the time of termination of the first imaging can be obtainedAs shown in formula (12)
During sliding beamforming imaging, the antenna beam center will always point at a certain rotational center. Therefore, the rotation point is necessarily at the vectorAnd the center of rotationThe following relationship must be satisfied:
based on formula (13) and the result obtained in step 1Vector, the center of rotation can be determinedThe position of (a). Fig. 5 shows the determination of the position of the beam pointing virtual point during the first imaging and the position of the beam pointing at the end of the imaging.
(2) Step 2: determining a fast maneuvering adjustment mode of a platform between two adjacent images
As can be seen from fig. 2, the maneuvering angular velocity of the satellite platform during the maneuvering process between adjacent images is significantly greater than the angular velocity during the imaging. The course of the angular speed during the maneuver is thus an angular speed ω from the moment of termination of the imaging phase 11eAccelerating and then decelerating to the angular velocity omega at the starting moment of the imaging stage 22sTo ensure the continuity and the smoothness of the working process. The relevant parameters that the process needs to determine are listed below:
beam rotation angular velocity acceleration time: t is1；
Acceleration during beam rotational acceleration: a is1；
Beam rotation acceleration deceleration time: t is2；
Acceleration during beam rotation deceleration: a is2
The allowable time of attitude maneuver of the satellite platform is T based on the given satellite platform in the attached table 1mAfter the attitude quick maneuver adjustment, the beam center squint angle theta at the initial moment of the second sliding beam bunching imaging2sIs composed of
The beam rotation angular velocity ω at that time2sIs composed of
Based on the beam center squint angle at the second starting time obtained by the equations (14) and (15) and the required beam rotation angular velocity, a system of equations for parameters such as acceleration time, acceleration, deceleration time, acceleration during deceleration, etc. of the maneuvering process can be established, as listed in equation (16)
The first equation represents that the sum of the acceleration time and the deceleration of the rotation of the satellite platform is the total maneuvering time T allowed by the platformm(ii) a The second equation shows that the beam rotation angular velocity at the initial time of the imaging phase 2 is exactly ω after acceleration and deceleration processes from the beam rotation velocity at the end time of the imaging phase 12s(ii) a The third equation shows that the angle through which the antenna beam turns throughout the maneuver is exactly the difference between the azimuthal squint angle from the end of phase 1 and the start of phase 2. The above three equations contain 4 unknowns, so it is necessary to design an acceleration time T first1And the other three parameters are determined. The corresponding analytical solution is:
T2＝Tm-T1 (17)
wherein T is used for ensuring the stable maneuvering process of the whole satellite1Is to ensure a1And a2The difference between them is sufficiently small.
Maneuvering time T based on platform between two adjacent imagingmIt needs to be determined at TmThe beam of the satellite platform is pointed during the maneuvering time. Fig. 6 shows the process of the satellite platform rapidly maneuvering to realize the antenna beam pointing from the scene terminal to the beginning.
In the context of figure 6 of the drawings,for the platform position at the end of the first sliding spotlight imaging,is the intersection of the antenna beam center and the ground at that moment.The platform position after the maneuvering process is completed (namely the platform position at the starting moment of the second sliding bunching imaging),the antenna beam center is pointed at the intersection point with the ground at that moment.
Andcan represent the relationship betweenIs composed of
By passingAndthe position of the satellite at the end of the maneuvering process can be obtainedIn obtainingAndafter three position vectors, the position vector needs to be determined on the groundTo ensure that the beam center vectors of the ending moment of the first sliding bunching and the starting moment of the second sliding bunching can have an intersection point (the antenna beam center always points to the point in the whole maneuvering process). If there is an intersection point between the two beam center vectors, thenAndfour points must be on one plane. Will be determined by these three points
First of all, the first step is to,must be associated with the sceneIs located at a distance from the front edge
Wherein theta is1eThe beam center squint angle at the end moment of the previous sliding beam bunching imaging is shown.
FIG. 7 shows the intersection point of the antenna beam center and the ground at two moments in timeAndis a positional relationship. From this figure, it can be seen that the determinationTwo steps are required: firstly, a vector parallel to the azimuth edge of the scene is made at a position which is at a distance D from the azimuth front edge of the sceneThe second step is to determine the vectorAnd plane surfaceThe intersection point of the antenna beam center pointing to the ground at the moment of termination of the maneuver
In thatAfter the determination, the intersection point of the beam center vector at the first sliding bunching termination time and the beam center vector at the second sliding bunching starting time can be determinedWhen in useAfter the determination, the beam center is always pointed in the whole platform maneuvering processTherefore, the direction of the Z axis of the platform can be obtainedNext, the orientations of the Y-axis and the X-axis of the stage can be obtained according to the same method as step 1. FIG. 8 shows the design flow of the fast maneuvering mode of the satellite platform between two adjacent sliding bunching imaging.
And (3) after the platform maneuvering process in the step (2) is finished, obtaining the beam center squint angle at the initial moment of the second sliding beam bunching imaging.
(3) And step 3: platform attitude maneuver adjustment mode for determining next sliding bunching imaging
After the squint angle and the rotation angular velocity at the initial stage of the 2 nd sliding bunching imaging are obtained in step 2, the following key parameters need to be determined in step 3:
the flying distance of the platform for the 2 nd sliding convergence imaging is L2s；
Beam center squint angle theta at 2 nd sliding convergence imaging termination moment2e；
Beam rotation angular velocity ω at the end time of the 2 nd sliding-focusing imaging2e。
Azimuth oblique viewing angle theta based on 2 nd sliding convergence imaging starting moment2sThe distance L traveled by the antenna beam on the ground at this stage can be derived2The sum of the azimuth length of the scene and the azimuth coverage length of the beam is:
then, by using the sliding bunching factor A, the flying distance L of the satellite platform in the sliding bunching stage 2 can be obtaineds2Is composed of
Oblique view angle theta at starting moment of 2 nd sliding bunching imaging phase2sA platform flight distance Ls2And the slant distance R of the sliding bunching rotation centertotThen the azimuthal squint angle theta of the stage at the end time can be obtained2eIs composed of
Finally, the beam squint angle theta at the termination time of the sliding bunching stage 2 is obtained according to the formula (25)2eThen, the beam rotation angular velocity ω corresponding to the time can be obtained2eIs composed of
And after the second sliding bunching imaging process is finished, the step 2 is performed again to determine the platform attitude maneuver process between the second sliding bunching imaging and the third sliding bunching imaging. The two steps are alternately carried out (until the azimuth oblique angle of the antenna beam exceeds the effective working range), and the satellite platform attitude requirements of each sliding beam-bunching imaging and maneuvering processes between two adjacent imaging can be obtained.
In the implementation process of the algorithm, the following process may be specifically performed:
1) according to the carrier frequency wavelength lambda of the SAR system and the flight speed of the satelliteDegree VsBeam footprint velocity V in stripe modegAnd SAR loaded antenna azimuth beam width thetaazAnd calculating the azimuth resolution rho in the strip modea_strip；
2) According to the azimuth resolution rho achieved under the multi-azimuth sliding bunching observation modeaAnd ρ obtained in the step 1)a_stripObtaining a resolution improvement factor A of the sliding bunching mode;
3) side-looking forward slope distance R according to scene centercAnd 2) obtaining a sliding bunching resolution improvement factor A obtained in the step 2), and obtaining a sliding bunching rotation center slope distance Rtot；
4) Starting squint angle theta according to current sliding bunching imaging1sAzimuth beam width theta of antennaaFront side view slope distance R of scene centercObtaining the length L of the area covered in the direction at the momentb；
5) From LbAnd the azimuth length L of the sceneaDetermining the travel distance L of the antenna beam footprint in the first sliding beam bunching imaging process1；
6) L obtained from step 5)1And determining the flight distance L of the platform in the first imaging process by using the sliding bunching improvement factor As1And a working time Tacq；
7) From LbAnd the azimuth length L of the sceneaDetermining the intersection point of the antenna beam center and the ground at the initial imaging moment of the satellite
8) Initial squint angle theta by sliding bunch imaging1sSatellite flight velocity VsAnd beam footprint velocity V in stripe modegAnd obtained in step 7)The position of the satellite imaging initial moment can be obtained
9) Starting time position by satellite imagingIntersection of antenna beam center and groundSatellite Z-axis pointing for determining starting time
10) Starting time position by satellite imagingSatellite Z-axis pointing at start timeAnd determining a virtual rotation center position by using the sliding bunching resolution improvement factor A
11) In the first sliding spotlight imaging process, the Z axis of the satellite platform always pointsThe Y-axis direction of the satellite is directed from the Z-axis direction and the satellite flying speed VsThe direction is determined by a right-hand rule, and the X-axis direction is determined by a Y-axis and a Z-axis through the right-hand rule;
12) according to the position of the starting moment of the first imaging of the satelliteLength of imaging time TacqAnd the flight orbit, the satellite position at the first imaging end time can be obtained(i.e. the)；
13) Obtained in step 12)And virtual center of rotation for a first imaging procedureThe attitude of the satellite at the termination time of the first imaging and the intersection point of the satellite and the ground are obtained
14) Obtained in step 12)And given the attitude maneuver time T between two imagingmDetermining the position of the satellite after maneuvering (the start time of the next imaging)
15) ByAnddetermining the plane determined by the three points and the distance D between the center of the beam and the scene after the maneuvering so as to determine the intersection point of the center of the beam and the ground after the maneuvering is finished
16) Satellite position after completion of maneuverAnd intersection of beam center and groundObtaining the Z-axis direction of the satellite at the maneuvering ending moment;
17) by vectorAndto determine the virtual center of rotation of the satellite platform during the maneuvering between two adjacent images
18) By virtual centre of rotation of a motorised processDetermining the three-axis attitude of the satellite in the whole maneuvering process;
19) byAnd determining the virtual rotation center of the second sliding bunching imaging by using the sliding bunching resolution improvement factor A
20) Entering the 11) step to the 14) step to carry out next sliding bunching imaging;
21) entering the steps 15) to 18) to carry out the next attitude maneuver adjustment process;
22) and (5) alternating the step (20) to the step (21) until the imaging work is finished.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the present invention, and those skilled in the art can make variations and modifications of the present invention without departing from the spirit and scope of the present invention by using the methods and technical contents disclosed above.

Claims (10)

1. A satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation is characterized by being realized in the following mode:
s1, determining the attitude maneuver mode of the satellite platform for the first imaging, namely determining the position of the satellite platform at the initial imaging moment, the satellite attitude and the virtual rotation center in the first imaging process;
s2, performing sliding bunching imaging, and determining the satellite position at the next imaging starting moment;
s3, determining satellite attitudes in the maneuvering process between two adjacent imaging and a virtual rotation center in the next imaging process according to the satellite position of the next imaging starting moment and the distance between the center of the beam after maneuvering and the scene, and adjusting the beam direction to the azimuth direction starting position of the scene in the attitude maneuvering adjustment process;
s4, repeating the steps s2 and s3 until the imaging work is finished;
in the sliding beamforming imaging process, the beam center of the satellite platform always points to the virtual rotation center determined in the imaging process.
2. The method according to claim 1, characterized in that the virtual center of rotation during the first imaging in step s1 is determined by:
initial squint angle theta by sliding bunch imaging1sSatellite platform flight speed VsAnd beam footprint velocity V in stripe modegAnd the intersection point of the antenna beam center and the ground at the initial imaging moment of the satelliteDetermining satellite imaging start time position
Starting time position by satellite imagingIntersection of antenna beam center and groundSatellite Z-axis pointing to determine starting imaging time
Starting time position by satellite imagingSatellite Z-axis pointing at start timeAnd determining a virtual rotation center position by using the sliding bunching resolution improvement factor A
3. The method of claim 2, wherein: the satellite imaging start time positionThe determination is made by:
determining the satellite imaging initial time, the intersection point of the antenna beam center and the groundAnd scene centerThe distance in the azimuth direction is further determined, and the intersection point of the antenna beam center and the ground is further determinedA location;
according to the intersection point of the antenna beam center and the groundPosition, determining the intersection pointSatellite position at zero doppler time
Determining connectionsAndvector and join of two pointsAndangle theta between vectors of two pointsmax；
According to the angle thetamaxDetermining the imaging start time position
Above, LaIs the azimuth length of the scene, RcPositive side-view skew distance, theta, at zero Doppler time from the center of the sceneazThe azimuth beamwidth of the antenna.
4. The method of claim 3, wherein: the included angle thetamaxAccording to the relationAnd (4) determining.
5. The method of claim 2, wherein: the virtual rotation center positionThe following relationship is satisfied:
to be started from a satellitePointing to first sliding spotlight imaging virtual center of rotationThe vector of (2).
6. The method of claim 1, wherein: step s2 is implemented by:
determining the satellite position of the imaging ending time according to the position of the imaging starting time of the satellite, the imaging time length and the flight orbit;
determining the satellite attitude at the imaging termination time and the intersection point of the satellite attitude and the ground according to the determined satellite position at the imaging termination time and the virtual rotation center of the imaging process;
and determining the satellite position at the moment of starting the next imaging after maneuvering according to the determined satellite position at the moment of ending the imaging and given attitude maneuvering time between two times of imaging.
7. The method of claim 1, wherein: step s3 is implemented by:
from the satellite position at the end of this imagingSatellite position at the start of next imagingAnd the intersection point of the beam center and the ground at the imaging termination momentDetermining the plane determined by the three points and the distance D between the center of the beam and the scene after maneuvering to determine the intersection point between the center of the beam and the ground after the attitude maneuvering is adjusted
Satellite position adjusted by attitude maneuverAnd intersection of beam center and groundObtaining the Z-axis direction of the satellite at the maneuvering ending moment;
by vectorAndto determine the virtual center of rotation of the satellite platform during the maneuvering between two adjacent images
By virtual centre of rotation of a motorised processDetermining the three-axis attitude of the satellite in the whole maneuvering process;
byAnd determining the virtual rotation center of the next sliding bunching imaging by using the sliding bunching resolution improvement factor A
8. The method of claim 7, wherein: the distance D is calculated according to the following formula:
wherein theta isendThe oblique angle of the beam center at the end of the sliding beam bunching imaging is RcPositive side-view skew distance, theta, at zero Doppler time from the center of the sceneazThe azimuth beamwidth of the antenna.
9. The method of claim 7, wherein: intersection of the beam center and the groundIs determined by the following steps:
firstly, a vector parallel to the azimuth edge of the scene is made at a position which is at a distance D from the azimuth front edge of the scene
Then, the vector is determinedAnd plane surfaceThe intersection point of the antenna beam center pointing to the ground at the moment of termination of the maneuver
As described aboveThe plane is the satellite position of the current imaging end timeSatellite position at the starting time of next imagingAnd the intersection point of the beam center and the ground at the current imaging end momentThe plane determined by the three points.
10. The method of claim 1, wherein: the beam rotation angular velocity in the attitude maneuver adjustment process is greater than that in the sliding beamforming imaging process.
CN201910765450.2A 2019-08-19 2019-08-19 Satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation Active CN110596704B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910765450.2A CN110596704B (en) 2019-08-19 2019-08-19 Satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910765450.2A CN110596704B (en) 2019-08-19 2019-08-19 Satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation

Publications (2)

Publication Number Publication Date
CN110596704A true CN110596704A (en) 2019-12-20
CN110596704B CN110596704B (en) 2021-10-01

Family

ID=68854916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910765450.2A Active CN110596704B (en) 2019-08-19 2019-08-19 Satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation

Country Status (1)

Country Link
CN (1) CN110596704B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111786087A (en) * 2020-08-07 2020-10-16 上海卫星工程研究所 Earth data transmission antenna layout method suitable for inter-satellite transmission
CN112014840A (en) * 2020-07-30 2020-12-01 西安空间无线电技术研究所 On-orbit implementation design method of satellite-borne SAR mosaic mode
CN112462365A (en) * 2020-09-21 2021-03-09 北京理工大学 Configuration optimization design method for acquiring satellite-borne scene matching SAR data
WO2023285432A1 (en) * 2021-07-14 2023-01-19 Iceye Oy Satellite with spot light mode for extended duration target imaging

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101513939B (en) * 2009-04-03 2011-01-05 北京航空航天大学 Two dimentional attitude control system of synthetic aperture radar satellite
CN103809178A (en) * 2014-01-17 2014-05-21 西安空间无线电技术研究所 Method for geosynchronous orbit synthetic aperture radar to realize continuous observation of coverage area
CN103076607B (en) * 2013-01-04 2014-07-30 北京航空航天大学 Method for realizing sliding spotlight mode based on SAR (Synthetic Aperture Radar) satellite attitude control
CN104090489A (en) * 2014-07-02 2014-10-08 中国科学院长春光学精密机械与物理研究所 Flexible agile satellite attitude maneuver rolling optimization control method
EP2816371A1 (en) * 2013-12-13 2014-12-24 Institute of Electronics, Chinese Academy of Sciences Method and device for steering attitude of satellite carrying synthetic aperture radar
CN106291557A (en) * 2016-08-30 2017-01-04 西安空间无线电技术研究所 A kind of satellite platform attitude maneuver method realizing satellite-borne SAR ultrahigh resolution slip beam bunching mode
CN107290961A (en) * 2017-06-29 2017-10-24 中国人民解放军国防科学技术大学 A kind of on-line scheduling method for quick satellite
CN107300699A (en) * 2016-04-15 2017-10-27 北京空间飞行器总体设计部 Mosaic mode implementation method based on quick Synthetic Aperture Radar satellite attitude maneuver
CN108267736A (en) * 2017-12-20 2018-07-10 西安空间无线电技术研究所 A kind of GEO SAR staring imagings mode orientation fuzziness determines method
CN109039418A (en) * 2018-06-15 2018-12-18 上海卫星工程研究所 Moonlet cluster network suitable for the real-time continuous monitoring in the space-based whole world
CN109507665A (en) * 2018-10-30 2019-03-22 北京空间飞行器总体设计部 It is a kind of based on spaceborne AIS real time information guidance star on autonomous imaging method

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101513939B (en) * 2009-04-03 2011-01-05 北京航空航天大学 Two dimentional attitude control system of synthetic aperture radar satellite
CN103076607B (en) * 2013-01-04 2014-07-30 北京航空航天大学 Method for realizing sliding spotlight mode based on SAR (Synthetic Aperture Radar) satellite attitude control
EP2816371A1 (en) * 2013-12-13 2014-12-24 Institute of Electronics, Chinese Academy of Sciences Method and device for steering attitude of satellite carrying synthetic aperture radar
CN103809178A (en) * 2014-01-17 2014-05-21 西安空间无线电技术研究所 Method for geosynchronous orbit synthetic aperture radar to realize continuous observation of coverage area
CN104090489A (en) * 2014-07-02 2014-10-08 中国科学院长春光学精密机械与物理研究所 Flexible agile satellite attitude maneuver rolling optimization control method
CN107300699A (en) * 2016-04-15 2017-10-27 北京空间飞行器总体设计部 Mosaic mode implementation method based on quick Synthetic Aperture Radar satellite attitude maneuver
CN106291557A (en) * 2016-08-30 2017-01-04 西安空间无线电技术研究所 A kind of satellite platform attitude maneuver method realizing satellite-borne SAR ultrahigh resolution slip beam bunching mode
CN107290961A (en) * 2017-06-29 2017-10-24 中国人民解放军国防科学技术大学 A kind of on-line scheduling method for quick satellite
CN108267736A (en) * 2017-12-20 2018-07-10 西安空间无线电技术研究所 A kind of GEO SAR staring imagings mode orientation fuzziness determines method
CN109039418A (en) * 2018-06-15 2018-12-18 上海卫星工程研究所 Moonlet cluster network suitable for the real-time continuous monitoring in the space-based whole world
CN109507665A (en) * 2018-10-30 2019-03-22 北京空间飞行器总体设计部 It is a kind of based on spaceborne AIS real time information guidance star on autonomous imaging method

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
LIQIANG XU，等: "The attitude tracking maneuvers of spaceborne spotlight SAR", 《2010 3RD INTERNATIONAL SYMPOSIUM ON SYSTEMS AND CONTROL IN AERONAUTICS AND ASTRONAUTICS》 *
LIU YINGYING，等: "Fuzzy attitude control for flexible satellite during orbit maneuver", 《2009 INTERNATIONAL CONFERENCE ON MECHATRONICS AND AUTOMATION》 *
T. LONG等: "A New Method of Zero-Doppler Centroid Control in GEO SAR", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 *
X. HAN,等: "Implementation method of Mosaic mode based on satellite attitude maneuver", 《2016 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM (IGARSS)》 *
唐文国，等: "一种遥感卫星宽幅无盲区拼接成像路径自主规划方法研究", 《上海航天》 *
张薇，等: "高轨SAR卫星在综合减灾中的应用潜力和工作模式需求", 《航天器工程》 *
李永昌: "敏捷卫星相机像移补偿关键技术研究", 《中国博士学位论文全文数据库 信息科技辑》 *
梁健，等: "敏捷SAR卫星序贯图像成像模式姿态机动策略研究", 《北京力学会.北京力学会第二十三届学术年会会议论文集》 *
章登义，等: "一种面向区域目标的敏捷成像卫星单轨调度方法", 《武汉大学学报(信息科学版)》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112014840A (en) * 2020-07-30 2020-12-01 西安空间无线电技术研究所 On-orbit implementation design method of satellite-borne SAR mosaic mode
CN111786087A (en) * 2020-08-07 2020-10-16 上海卫星工程研究所 Earth data transmission antenna layout method suitable for inter-satellite transmission
CN112462365A (en) * 2020-09-21 2021-03-09 北京理工大学 Configuration optimization design method for acquiring satellite-borne scene matching SAR data
CN112462365B (en) * 2020-09-21 2024-05-24 北京理工大学 Space-borne scene matching SAR data acquisition configuration optimization design method
WO2023285432A1 (en) * 2021-07-14 2023-01-19 Iceye Oy Satellite with spot light mode for extended duration target imaging
GB2608851B (en) * 2021-07-14 2024-04-10 Iceye Oy Satellite with spot light mode for extended duration target imaging

Also Published As

Publication number Publication date
CN110596704B (en) 2021-10-01

Similar Documents

Publication Publication Date Title
CN110596704B (en) Satellite platform attitude maneuver method for satellite-borne SAR multi-azimuth repeated observation
CN107132537B (en) A kind of SAR satellite on-orbit performance method for improving based on electromechanical combination scanning
CN110515078B (en) Wave position design method for airspace coverage
CN107300699B (en) Method for realizing mosaic mode based on agile synthetic aperture radar satellite attitude maneuver
CN107390181B (en) Radar high-resolution imaging method based on multi-beam scanning
CN113589285B (en) SAR real-time imaging method for aircraft
US5442364A (en) Alignment and beam spreading for ground radial airborne radar
CN110208797B (en) Quick-response SAR satellite high squint attitude maneuver method
CN115792907B (en) Method for designing azimuth imaging parameters of spaceborne SAR squint sliding bunching mode
CN110208801B (en) Universal SAR imaging PRF optimization design method
CN108613655B (en) Attitude adjustment method for imaging along inclined strip in agile satellite machine
CN112255606A (en) Method for calculating front side-view imaging attitude angle of Geo-SAR (synthetic aperture radar) satellite based on single reflector antenna
CN112014840B (en) On-orbit implementation design method of satellite-borne SAR mosaic mode
JP5298730B2 (en) Interference synthetic aperture radar system, processing method, directivity angle correction apparatus, directivity angle correction method, and program
WO2023117390A1 (en) High resolution wide swath sar imaging
CN116087953B (en) Satellite-borne SAR multi-target imaging parameter design method
CN118052081B (en) Parameter design method for high-orbit SAR system
RU2821956C1 (en) System and method for tracking antenna system of mobile satellite earth station
CN118401859A (en) High resolution wide swath SAR imaging
US20240150040A1 (en) Satellite with spot light mode for extended duration target imaging
CN111948650A (en) Satellite-borne bistatic SAR (synthetic Aperture Radar) combined Doppler guidance method based on electric scanning
AU2022420366A1 (en) Multi-spot imaging using synthetic aperture radar
CN117233757A (en) Satellite-borne SAR multi-angle scanning method with cooperative phased array electric scanning and attitude
CN118401860A (en) Multi-point imaging using synthetic aperture radar
KR20240117543A (en) Multi-point imaging using synthetic aperture radar

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant