ADVANCES AND LIMITATIONS IN A PARAMETRIC GEOMETRIC CORRECTION OF CHRIS/PROBA DATA
|
|
- Tabitha Tyler
- 5 years ago
- Views:
Transcription
1 ADVANCES AND LIMITATIONS IN A PARAMETRIC GEOMETRIC CORRECTION OF CHRIS/PROBA DATA Luis Alonso, José Moreno LEO Group, Dpt. de Termodinamica, Univ. Valencia, Dr. Moliner 50, Burjassot, Spain, luis.alonso@uv.es ABSTRACT/RESUME The geometric correction of remote sensing images is a preprocessing step required for most applications and studies. The success of CHRIS/PROBA as scientific mission is producing a large archive of image sets that is continuously growing. Although the geometric correction problem is already solved by means of ground control points (GCP) or image corregistration techniques, these methods require intense human intervention. Therefore they are not convenient for processing large number of images. For this reason the development of a quasi-automatic parametric algorithm for geometric correction of CHRIS/PROBA data was undertaken by our group. This work has given rise to some questions about the image acquisition procedure. As a result new information about the operation of CHRIS/PROBA has recently become available for the user community, pointing out that the operation for acquisition is based on time, instead of fixed fly-by zenith angles (FZAs), as was documented previously. It also states that to determine each particular acquisition the system assumes a circular orbit with radius equal to the semi-major axis. This paper presents a study of the implication that the new information has on how the parametric algorithm must be implemented, as well as the impact on the acquisition geometry, the observation angles, and the pointing accuracy. 1. INTRODUCTION CHRIS/PROBA is the first hyperspectral sensor with pointing capabilities and high spatial resolution. Initially designed as a technology demonstrator, the uniqueness of the data obtained has probed to be very useful for scientific studies. In order to make full use of the CHRIS/PROBA multiangular capability, all five images from a given acquisition set must be accurately co-registered to make sure that the values measured from a given coordinate actually correspond to the same surface for the five angles. To achieve this, it is necessary an accurate geometric correction or geolocation algorithm. Every pixel of each image must have a geographic coordinate assigned. The current algorithms based on Ground Control Points GCP cannot assure this level of co-registration between acquisitions for different reasons: the need of a homogeneous distribution of GCPs throughout the image, the interpolation process between GCPs unable to compene the local non-linear distortions within the image, and in particular, the reduced across-track spatial resolution of the larger Fly-by Zenith Angle (FZA) images makes it more difficult to assign GCPs (subpixel precision would be required). For the Barrax test site in La Mancha (Spain), which is optimal for this type of approach due to its flatness, the number and size of fields, and landscape features that allow setting a large number of GCPs, there are portions of the images that have up to three pixels of error in the co-registration, corresponding to approximately 100m. Therefore, in those portions of the image, only fields larger than 200m will have enough number of pixels coregistered in the five observation angles. Also, the distortions introduced by the GCP correction are inhomogeneously distributed throughout the scene. This method is particularly problematic for water studies, due to the impossibility to set GCP within the water area; especially difficult is the case of coastal waters where GCP can only be assigned to a portion (sometimes small) of the image. Therefore, an algorithm for accurate geolocation is needed. Such algorithm was proposed by our group in the 2nd CHRIS/PROBA workshop [1], taking a parametric approach that would reproduce in detail the actual ellite orbit and the acquisition process, based on the knowledge available about the operation of the platform and the sensor. 2. PARAMETRIC ALGORITHM A parametric geometric correction would allow geolocation with higher accuracy, without introducing inhomogeneous distortions in the images. In flat study sites where topography has no impact in the geometry of the problem, a simple parametric method (without the use of a digital elevation model (DEM)) would provide almost direct corrections, with very little adjustments and supervision. The work done to develop this algorithm raised a number of questions about details of the operation of the Proc. of the 3rd ESA CHRIS/Proba Workshop, March, ESRIN, Frascati, Italy, (ESA SP-593, June 2005)
2 ellite. These questions were answered in a recent communication by ESA to the CHRIS/PROBA users [2]. In this communication the actual acquisition procedure was detailed, and it strongly differed from what was commonly known. The algorithm was based on the assumption that the images were acquired at fixed observation angles, depicted in Figure 1. According to the documentation available at the time [3], the observation angles of each image were defined by the Fly-by Zenith Angle (FZA) and the Minimum Zenith Angle (MZA) [1,3], and the images were acquired in such a way that the centre of the images corresponded to a FZA of 0º, ±36º and ±55º; and the acquisition time was 10s for each image. According to recently released documentation it seems that this is just a simplified description of the actual procedure, which in reality operates on a time basis, instead of a fixed angles basis. FZA = +55 FZA = +36 FZA = -55 orbit FZA = -36 ground track FZA = 0 MZA N target site Figure 1. Nominal FZA and its relation to azimuth and zenith observation angles. The first documentation on CHRIS/PROBA defined the image acquisition to be angle based. To determine and program the operation of the ellite for the acquisition procedure the system follows a series of steps: Using the ellite s position and velocity at 390s prior to acquisition start the semi-major axis of the orbit and the angular velocity are calculated. From this moment until the end of the acquisition the orbit is assumed to be circular, with radius equal to the calculated semi-major axis. Then the five acquisitions are defined to take place within a cone of 55º centred on the target. The total time of observation (the time elapsed from the start of the first image acquisition to the end of the last image acquisition) is the time that the platform needs to cross this cone. The total observation time (Tobs) is then divided into scan time (Tsc), slew time (Tsl) and margin time (Tmar). Scan time is fixed to 20s for the image acquisition, slew time is dedicated to reorient the ellite for the next acquisition in 12.5s, and margin time is added before and after scan time to damp transients and its duration is determined by the other elements of Tobs. Tobs = Tsc + Tmar + 1st acquisition Tsl + Tmar + Tsc + Tmar + 2nd acquisition Tsl + Tmar + Tsc + Tmar + 3rd acquisition Tsl + Tmar + Tsc + Tmar + 4th acquisition Tsl + Tmar + Tsc 5th acquisition Once the timing is determined, the??? and Control System calculates the necessary attitude change rates to keep the target on sight while scanning the area with a slowdown factor of 5 to increase the integration time. The actual onboard implementation of the pointing manoeuvre restricts the rotation to two axis: pitch and roll. The above procedure results in a set of different FZA angles for different altitudes of the ellite, according to Chris Data Format Issue4 [4] reflected in Table 1. Table 1. Change of FZA with ellite altitude Max altitude Min. altitude ±55 FZA angles ( ) ±49.52 ±56.18 ±36 FZA angles ( ) ±29.27 ±32.02 Inclination ( ) Period (min) Altitude (km) NEW APPROACH FOR PARAMETRIC MODEL Then new information about the acquisition procedure implies a radical change in the whole concept of how CHRIS/PROBA works. The parametric approach is still feasible but the complete structure of the model had to be adapted to simulate the actual acquisition procedure. One advantage of the new approach is that the inputs required by the model have been r reduced to: - Satellite position ( R ) at orbit fixation time. r - Velocity vector ( V ) at orbit fixation time. - Target coordinates ( R r tgt ) These data permits to calculate all the parameters that will determine the ellite evolution with respect to the target in the following minutes [5].
3 One of the most relevant parameters is the semi-major axis, obtained with Eq. 1, where GM is the gravitation constant of the Earth. 1 = 2 R r V r (1) 2 GM The angular momentum H defines the orbital plane and its rotation axis: r H = R r V r (2) Using the semi-major axis a it is possible to calculate the angular velocity ω of the ellite and the φ angle, necessary to calculate the total observation time T obs and the time evolution of the platform: ω = GM 3 φ = 55º asin cos(35º) r a max R tgt (3) (4) Finally, it is possible to obtain the observation time, that determines the whole acquisition process: T obs = 2 φ ω (5) 4. TEST CASE FOR THE MODEL To test the algorithm we have focused on the acquisition of the 12th of July, For this date we have available the telemetry data and the corresponding TLE. These data allow checking the consistency of the method, as both provide the necessary information to retrieve the ellite s position and velocity at a given time, as well as other orbital parameter while they are obtained with different methods. Using both types of data we have calculated the circular orbit that corresponds to the acquisition procedure, and have found a significant difference between the actual altitude of the platform and the corresponding semimajor axis at the moment of orbit fixation, shown in Table 2; obtaining similar results from telemetry and TLE propagation. Table 2. Difference between actual ellite height from Earth centre and calculated semi-major axis at the moment of orbit fixation, derived by the two methods. Units are in kilometers. Sat. height Difference: Telemetry TLE Estimation Difference: Ground height: The difference is of 49.9km in the case of telemetry data, and 47.6km in the case of TLE propagation. This difference that might seem small in comparison to the distance of the ellite to the Earth s centre turns very significant with respect to the altitude of the platform over the Earth s surface. In this test case it represents a decrement of 7.8% on the altitude above ground of the ellite, see Fig. 3. Because the whole determination of the ellite pointing and imaging programming depends on the assumption of a circular orbit of radius equal to the semi-major axis the difference between and the actual altitude can have a big impact in the accurate determination of the observation angles Altitude from Earth centre [km] actual altitude estimated Ground Telemetry sample Figure 2. Actual altitude of the ellite and the altitude of the corresponding semi-major axis. The vertical lines indicate the instants of orbit fixation (yellow), image acquisition start (green) and end (red). The situation worsens when we look at the evolution of the actual altitude of the ellite, which increases as it overpasses the target area, as illustrated in Fig. 3; the increment of altitude, with respect to the altitude at the moment of orbit fixation, is 15km at the start of image acquisition and 20km at the end of the acquisition process.
4 665 Sat height above ground Height [km] elliptical circular ~55º 55º km between 1st and 5th image 20km between orbit fix and acquisition Φ Telemtry Sample Figure 3. Variation in altitude of the ellite during the orbit determination and imaging process. These results are particular to this test case, but it is possible to determine how the situation would be for the extreme cases of the orbit in that date, that is, if the acquisition occurs at perigee or apogee. The orbit of PROBA had an excentricity of and a semimajor axis of km according to TLE, which has results in a perigee radius of 6914km, and an apogee radius of 7031km. In these cases the error in the orbit height would be of ±58km. This difference in height has a direct impact on the accuracy of the programmed observation angles, and could produce pointing problems during image generation, as can be seen in Fig. 4. Further analysis of archived telemetry will produce sounder results about the magnitude of the effects derived from the difference in the actual height and the estimated height used for operation. 5. ANALYSIS OF POINTING ACCURACY With the aim of checking if these differences in altitude were due to some computation error or they actually represent the present situation of CHRIS/PROBA acquisition the pointing accuracy was analyzed over whole collection of acquisitions over Barrax by determining the coordinates of the centre pixel of each image and their distance to the programmed target coordinates. Figure 4. Impact of using a circular orbit of radius for acquisition determination in an elliptical orbit of higher altitude. Diagram not at scale. We restricted the analysis to those acquisitions performed during 2004, after the enhancement in the pointing accuracy with the star-trackers, in order to include only those images with the best attitude control, i.e. highest pointing capability. We have also included in this analysis the images acquired over Rosarito test site (also in Spain) to increase the significance of the results. Fig. 6 (left) shows the position of the images centres acquired the 16/07/2004 over the test site Barrax; the blue circle indicates the programmed image centre, and the blue crosses indicate de limit of a ±0.01º area (the precision of the centre coordinates). In Fig. 6 (right) a graph shows the position of the image centers for all the Barrax acquisitions during The programmed centre is the red cross, the image centre of the FZA 0º acquisitions are distinguished by symbols with an outer frame, the ±36 are filled symbols, and the ±55 are empty symbols. The inclination of the orbit at the given latitude is displayed as a black dashed line as reference of the along-track direction. The centre image pixel coordinates of the Rosarito test site have been provided by the Center for Hydrographic Studies (CEDEX) to improve the analysis. These data are plotted in Fig. 7 in a similar way than the data from Barrax.
5 UTM Y [m] Programmed Centre 16/07/2004 MZA= -8º 15/07/2004 MZA= 19º 30/06/2004 MZA=-17º 27/05/2004 MZA= 14º 23/03/2004 MZA= 10º /07/ UTM X [m] Figure 6. Center of image position: Estimation of pointing accuracy. Barrax test site (Spain) May MZA=14 15Jun MZA=-20 27Aug MZA=4 23Oct MZA=19 26Nov MZA=-4 FZA=0 FZA=±55 Target 20/05/ UTM Easting [m] Figure 7. Center of image position: Estimation of pointing accuracy. Rosarito test site (Spain). As with all ellite imaging instruments, sources of uncertainty in CHRIS/PROBA pointing exist and lead to misalignments of the images with respect programmed pointing [4]. Typical sources of error include: Alignment (instrument and spacecraft) Datation, i.e. timing (error in time translate into errors on the computed earth rotation, hence on pointing) Limitation in the PROBA on-board feedback loop used to maintain maximum pointing accuracy during the acquisition (residual calibration errors on wheel axis and inertia ratio) The image centers alignment follow the orbit inclination rather well indicating a good across-track pointing, although they present a systematic deviation eastwards of the orbit track; this could be due to a problem of instrument alignment or a timing shift. It has also been detected a large along-track dispersion of the five images of each acquisition (FZA = 0º, ±36º
6 and ±55º), in some cases more than 7km, that is, more than half the size of the image. This result would be consistent with the effect of a large difference between the operational orbit height and the actual orbit height. It is also interesting to note that in many cases the centers of the ±55º images fall to the north of the 0º centre while the ±36º images centers are located to the south. This distribution cannot be explained with the arguments expressed before, which would produce an ordered dispersion. Data from other test sites would increase statistical significance of this problem, also allowing a thorough analysis of image centre coordinates, which can provide a deeper insight to improve PROBA s pointing accuracy. 6. ACTUAL TIME OF OBSERVATION Up to now we have only considered changes in altitude, in orbital mechanics a change in altitude implies a change in angular velocity, and this would affect the T obs required. It is possible to find a relationship between the observation time (T o obs) of a circular orbit that fits the actual orbit and the calculated T obs for image acquisition, using Eqs. 3 and 5 where T o obs m = = m T (6) R o implying that the time required in the actual orbit to travel an arc of 2 φ is larger than the T obs used for imaging. Hence the angular distance covered by PROBA in its orbit 2 φ in a time T obs would be obs 3/2 (7) 2 φ' = 2 φ / m (8) which is shorter than the distance needed to point the target correctly. This implies that the estimated attitude for the acquisition based on the reference orbit will be performed on an orbit that is actually at a different altitude, and which requires an observation time to cross the 55º cone that is different from the calculated T obs, illustrated in Fig.8. These two mis-estimations add up to produce an error in the pointing of the five angles that could be not necessarily consecutive (e.g. image centers located from north to south in a ±55º, 0º, ±36º, instead of -55º, -36º, 0º, +36º, +55º), contrary to the case of and appropriate T obs with wrong altitude, that would produce ordered pointing errors. circular estimated circular approximated R o 2 φ 55º Φ Figure 8. Diagram on the effect of an understimation of ellite height and time of observation. 7. CONCLUSIONS Image acquisition based on fixed time instead of fixed FZAs has serious implications in observation geometry assumptions, especially in the determination of observation angles. Parametric correction needs accurate information about location and velocity of the ellite at the orbit fixation before the image acquisition. Inputs can be retrieved from telemetry or from TLE (orbit propagation). It has been found that the assumption of a circular orbit with radius the semi-major axis could be inappropriate even for a near circular orbit like PROBA s, producing a mis-estimation of observation angles during the determination of the image acquisition procedure. Satellite pointing errors have been noticed with some systematic deviations consistent with the misestimation of the acquisition process. Further analysis could help to improve PROBA performance. Geometric correction and per pixel observation angles, our main goal, are still possible with the new approach. Despite the fact that the objective of our work is to obtain an algorithm for quasi-automatic geometric correction these results on the performance of CHRIS/PROBA acquisition procedure could be relevant for the improvement of the ellite. For this reason it now our intention to include in the model the possibility to allow changes in different parameters: - From the orbit (circular, elliptical) - From the Earth (rotation, oblanteness, etc) - From the procedure (integration time, slowdown factor, )
7 The nominal model could become a useful tool for sensitivity analysis on acquisition performance, or as simulator of alternative acquisition procedures. ACKNOWLEDGEMNT The authors want to thank R. Peña, J. A. Dominguez and A. Verdú for kindly providing the pointing data from the Rosarito CHRIS/PROBA test site. REFERENCES 1. Alonso L. and Moreno J., Quasi-Automatic Geometric Correction and Related Geometric Issues in the Exploitation of CHRIS/Proba Data, Proceedings of the second CHRIS/Proba Workshop, April 2004, ESA/ESRIN, Frascati, Italy (ESA SP-578, July 2004) 2. Davidson M. and Vuilleumier P., Note on CHRIS Acquisition Procedure and Image Geometry revision 1.3, ESA, Noordwijk, The Netherlands, Cutter M. A., CHRIS Data Format revision 2, SIRA, Kent, U.K., Cutter M. A., CHRIS Data Format revision 4, SIRA, Kent, U.K., Wertz J. R. and Larson W. J., Space mission analysis and design 3 rd ed., Dordrecht: Kluwer Academic Publishers, 1999
PROBA 1. F. Teston ESA/ESTEC D/TEC-EL
PROBA 1 F. Teston ESA/ESTEC D/TEC-EL Frederic.Teston@esa.int PROBA 1 launch PROBA 1 has been launched on 21 October 2001 Orbital parameters: Altitude: 681-561 km Near polar (inclination of 97.9 ) Sun-synchronous
More informationNEW EUMETSAT POLAR SYSTEM ATTITUDE MONITORING SOFTWARE
NEW EUMETSAT POLAR SYSTEM ATTITUDE MONITORING SOFTWARE Pablo García Sánchez (1), Antonio Pérez Cambriles (2), Jorge Eufrásio (3), Pier Luigi Righetti (4) (1) GMV Aerospace and Defence, S.A.U., Email: pgarcia@gmv.com,
More informationCircular vs. Elliptical Orbits for Persistent Communications
5th Responsive Space Conference RS5-2007-2005 Circular vs. Elliptical Orbits for Persistent Communications James R. Wertz Microcosm, Inc. 5th Responsive Space Conference April 23 26, 2007 Los Angeles,
More informationSatellite Communications
Satellite Communications Lecture (3) Chapter 2.1 1 Gravitational Force Newton s 2nd Law: r r F = m a Newton s Law Of Universal Gravitation (assuming point masses or spheres): Putting these together: r
More informationACHIEVING THE ERS-2 ENVISAT INTER-SATELLITE INTERFEROMETRY TANDEM CONSTELLATION.
ACHIEVING THE ERS-2 ENVISAT INTER-SATELLITE INTERFEROMETRY TANDEM CONSTELLATION M. A. Martín Serrano (1), M. A. García Matatoros (2), M. E. Engdahl (3) (1) VCS-SciSys at ESA/ESOC, Robert-Bosch-Strasse
More informationDETERMINATION OF SCIAMACHY LINE-OF-SIGHT MISALIGNMENTS
DETERMINATION OF SCIAMACHY LINE-OF-SIGHT MISALIGNMENTS Manfred Gottwald (1), Eckhart Krieg (1), Sander Slijkhuis (1), Christian von Savigny (2), Stefan Noël (2), Heinrich Bovensmann (2), Klaus Bramstedt
More informationMULTI PURPOSE MISSION ANALYSIS DEVELOPMENT FRAMEWORK MUPUMA
MULTI PURPOSE MISSION ANALYSIS DEVELOPMENT FRAMEWORK MUPUMA Felipe Jiménez (1), Francisco Javier Atapuerca (2), José María de Juana (3) (1) GMV AD., Isaac Newton 11, 28760 Tres Cantos, Spain, e-mail: fjimenez@gmv.com
More informationThe VENμS mission: Earth Observation with High Spatial and Temporal Resolution Capabilities
Ben Gurion University The VENμS mission: Earth Observation with High Spatial and Temporal Resolution Capabilities G. Dedieu 1, O. Hagolle 2, A. Karnieli 3, S. Cherchali 2 P. Ferrier 2 and Y. Yaniv 4 1
More informationOrbit Design Marcelo Suárez. 6th Science Meeting; Seattle, WA, USA July 2010
Orbit Design Marcelo Suárez Orbit Design Requirements The following Science Requirements provided drivers for Orbit Design: Global Coverage: the entire extent (100%) of the ice-free ocean surface to at
More informationAN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT
AAS 16-366 AN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT Jason A. Reiter * and David B. Spencer INTRODUCTION Collision avoidance maneuvers to prevent orbital
More informationCHAPTER 3 PERFORMANCE
PERFORMANCE The launch performance given in this chapter is based on the following assumptions: The LV system parameters being all nominal values; Mass of the LV adapter and the separation system are included
More informationFive Years Orbit Experience of a Small Satellite Hyperspectral Imaging Mission
SSC07-VIII-3 Five Years Orbit Experience of a Small Satellite Hyperspectral Imaging Mission Dr Mike A Cutter & Prof Sir Martin Sweeting Surrey Satellite Technology Ltd, Tycho House, 20 Stephenson Road,
More informationGuillem Sòria, José A. Sobrino, Juan C. Jiménez-Muñoz, Mónica Gómez, Juan Cuenca, Mireia Romaguera and Malena Zaragoza
Guillem Sòria, José A. Sobrino, Juan C. Jiménez-Muñoz, Mónica Gómez, Juan Cuenca, Mireia Romaguera and Malena Zaragoza University of Valencia, Spain MERIS (A)ATSR WORKSHOP ESRIN, Frascati, Italy 26-30
More informationAppendix B. A proposition for updating the environmental standards using real Earth Albedo and Earth IR Flux for Spacecraft Thermal Analysis
19 Appendix B A proposition for updating the environmental standards using real Earth Albedo and Earth IR Romain Peyrou-Lauga (ESA/ESTEC, The Netherlands) 31 st European Space Thermal Analysis Workshop
More informationAutonomous Vision Based Detection of Non-stellar Objects Flying in Formation with Camera Point of View
Autonomous Vision Based Detection of Non-stellar Objects Flying in Formation with Camera Point of View As.Prof. M. Benn (1), Prof. J. L. Jørgensen () (1) () DTU Space, Elektrovej 37, 4553438, mb@space.dtu.dk,
More informationANALYSIS OF THE MEGHA-TROPIQUES TRAJECTORY. DETERMINATION OF RENDEZ-VOUS CONDITIONS WITH THE TERRA SATELLITE.
ANALYSIS OF THE MEGHA-TROPIQUES TRAJECTORY. DETERMINATION OF RENDEZ-VOUS CONDITIONS WITH THE TERRA SATELLITE. Michel Capderou (1), Florent Deleflie (2), and Olivier Chomette (1) (1) LMD-IPSL Ecole Polytechnique,
More informationFundamentals of Satellite technology
Fundamentals of Satellite technology Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Orbital Plane All of the planets,
More informationSatellite meteorology
GPHS 422 Satellite meteorology GPHS 422 Satellite meteorology Lecture 1 6 July 2012 Course outline 2012 2 Course outline 2012 - continued 10:00 to 12:00 3 Course outline 2012 - continued 4 Some reading
More informationMODULE 2 LECTURE NOTES 1 SATELLITES AND ORBITS
MODULE 2 LECTURE NOTES 1 SATELLITES AND ORBITS 1. Introduction When a satellite is launched into the space, it moves in a well defined path around the Earth, which is called the orbit of the satellite.
More informationAttitude Control Simulator for the Small Satellite and Its Validation by On-orbit Data of QSAT-EOS
SSC17-P1-17 Attitude Control Simulator for the Small Satellite and Its Validation by On-orbit Data of QSAT-EOS Masayuki Katayama, Yuta Suzaki Mitsubishi Precision Company Limited 345 Kamikmachiya, Kamakura
More informationOrbits in Geographic Context. Instantaneous Time Solutions Orbit Fixing in Geographic Frame Classical Orbital Elements
Orbits in Geographic Context Instantaneous Time Solutions Orbit Fixing in Geographic Frame Classical Orbital Elements Instantaneous Time Solutions Solution of central force motion, described through two
More informationCHAPTER 3 PERFORMANCE
PERFORMANCE 3.1 Introduction The LM-3A performance figures given in this chapter are based on the following assumptions: Launching from XSLC (Xichang Satellite Launch Center, Sichuan Province, China),
More informationAUTONOMOUS AND ROBUST RENDEZVOUS GUIDANCE ON ELLIPTICAL ORBIT SUBJECT TO J 2 PERTURBATION.
AUTONOMOUS AND ROBUST RENDEZVOUS GUIDANCE ON ELLIPTICAL ORBIT SUBJECT TO J 2 PERTURBATION Emmanuel GOGIBUS (1), Hervé CHARBONNEL (2), Patrick DELPY (3) (1) Astrium Space Transportation, 66 route de Verneuil,
More informationCalibration of MERIS on ENVISAT Status at End of 2002
Calibration of MERIS on ENVISAT Status at End of 2002 Bourg L. a, Delwart S. b, Huot J-P. b a ACRI-ST, 260 route du Pin Montard, BP 234, 06904 Sophia-Antipolis Cedex, France b ESA/ESTEC, P.O. Box 299,
More informationDesign of Attitude Determination and Control Subsystem
Design of Attitude Determination and Control Subsystem 1) Control Modes and Requirements Control Modes: Control Modes Explanation 1 ) Spin-Up Mode - Acquisition of Stability through spin-up maneuver -
More informationAPPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING
APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING Orbit selection and sensor characteristics are closely related to the strategy required to achieve the desired results. Different types
More informationMODELLING OF PERTURBATIONS FOR PRECISE ORBIT DETERMINATION
MODELLING OF PERTURBATIONS FOR PRECISE ORBIT DETERMINATION 1 SHEN YU JUN, 2 TAN YAN QUAN, 3 TAN GUOXIAN 1,2,3 Raffles Science Institute, Raffles Institution, 1 Raffles Institution Lane, Singapore E-mail:
More informationProton Launch System Mission Planner s Guide SECTION 2. LV Performance
Proton Launch System Mission Planner s Guide SECTION 2 LV Performance 2. LV PERFORMANCE 2.1 OVERVIEW This section provides the information needed to make preliminary performance estimates for the Proton
More informationEnd of Life Re-orbiting The Meteosat-5 Experience
End of Life Re-orbiting The Meteosat-5 Experience Milan EUMETSAT, Darmstadt, Germany This article illustrates the orbit maneuver sequence performed during Meteosat- 5 End of Life (EOL) re-orbiting operations
More informationDEVIL CHAPTER 6 TEST REVIEW
IB PHYSICS Name: Period: Date: # Marks: 51 DEVIL PHYSICS BADDEST CLASS ON CAMPUS CHAPTER 6 TEST REVIEW 1. A cyclist rides around a circular track at a uniform speed. Which of the following correctly gives
More informationGeneration X. Attitude Control Systems (ACS) Aprille Ericsson Dave Olney Josephine San. July 27, 2000
Generation X Attitude Control Systems (ACS) Aprille Ericsson Dave Olney Josephine San July 27, 2000 ACS Overview Requirements Assumptions Disturbance Torque Assessment Component and Control Mode Recommendations
More informationPERFORMANCE OF THE FIRST ORBIT MANEUVER ON CBERS-1 SATELLITE BY INPE S CONTROL CENTER
Série Arquimedes, Volume 2, Anais do DINCON 2003, pp. 831-840 2º Congresso Temático de Aplicações de Dinâmica e Controle da Sociedade Brasileira de Matemática Aplicada e Computacional (SBMAC). São José
More information(ii) We have already found out the velocity of the satellite in orbit in part (i) (using equation (2.5)) to be km/s
Chapter 2 Questions and Solutions Question 1. Explain what the terms centrifugal and centripetal mean with regard to a satellite in orbit around the earth. A satellite is in a circular orbit around the
More informationGEOSTATIONARY/COLLOCATION MISSION ANALYSIS USING FOCUSGAT
GEOSTATIONARY/COLLOCATION MISSION ANALYSIS USING FOCUSGAT EUROEAN SACE ASTRONOMY CENTRE (ESA/ESAC), MADRID, SAIN 3 6 MAY 010 Isabel Vera Trallero (1), Javier Cuesta Cabanás (1), Ángel Luis Cerviño (1),
More informationAATSR derived Land Surface Temperature from heterogeneous areas.
AATSR derived Land Surface Temperature from heterogeneous areas. Guillem Sòria, José A. Sobrino Global Change Unit, Department of Thermodynamics, Faculty of Physics, University of Valencia, Av Dr. Moliner,
More informationAnalysis of Relative Motion of Collocated Geostationary Satellites with Geometric Constraints
www.dlr.de Chart 1 Analysis of Relative Motion of Collocated Geostationary Satellites with Geometric Constraints SFFMT2013, München F. de Bruijn & E. Gill 31 May 2013 www.dlr.de Chart 2 Presentation Outline
More informationCHAPTER 3 PERFORMANCE
PERFORMANCE 3.1 Introduction The LM-3B performance figures given in this chapter are based on the following assumptions: Launching from XSLC (Xichang Satellite Launch Center, Sichuan Province, China),
More informationAttitude Determination using Infrared Earth Horizon Sensors
SSC14-VIII-3 Attitude Determination using Infrared Earth Horizon Sensors Tam Nguyen Department of Aeronautics and Astronautics, Massachusetts Institute of Technology 77 Massachusetts Avenue, Cambridge,
More informationUsage of IGS TEC Maps to explain RF Link Degradations by Spread-F, observed on Cluster and other ESA Spacecraft
Usage of IGS TEC Maps to explain RF Link Degradations by Spread-F, observed on Cluster and other ESA Spacecraft J. Feltens, J. Dow, G. Billig, D. Fornarelli, S. Pallaschke, B. Smeds, H.-J. Volpp, P. Escoubet,
More informationIntroduction to Satellite Orbits
Introduction to Satellite Orbits Dinesh Manandhar Center for Spatial Information Science The University of Tokyo Contact Information: dinesh@iis.u-tokyo.ac.jp Slide : 1 Basics of Satellite Orbits The speed
More informationCOMPARISON OF ANGLES ONLY INITIAL ORBIT DETERMINATION ALGORITHMS FOR SPACE DEBRIS CATALOGUING
COMPARISON OF ANGLES ONLY INITIAL ORBIT DETERMINATION ALGORITHMS FOR SPACE DEBRIS CATALOGUING Fran Martinez Fadrique, Alberto Águeda Maté, Joan Jorquera Grau, Jaime Fernández Sánchez, Laura Aivar García
More informationRETRIEVAL OF VEGETATION BIOPHYSICAL VARIABLES FROM CHRIS/PROBA DATA IN THE SPARC CAMPAING
RETRIEVAL OF VEGETATION BIOPHYSICAL VARIABLES FROM /PROBA DATA IN THE SPARC CAMPAING S. Gandia, G. Fernández, J. C. García, J. Moreno Laboratory for Earth Observation Department of Thermodynamics. Faculty
More informationAATSR DERIVED LAND SURFACE TEMPERATURE OVER A HETEROGENEOUS REGION
AATSR DERIVED LAND SURFACE TEMPERATURE OVER A HETEROGENEOUS REGION José A. Sobrino (1), Guillem Sòria (1), Juan C. Jiménez- Muñoz (1), Belen Franch (1), Victoria Hidalgo (1) and Elizabeth Noyes (2). (1)
More informationPW-Sat two years on orbit.
13th of February 2014 is the second anniversary of launch of the first polish student-made satellite PW-Sat. Currently Students' Space Association on Warsaw University of Technology is working on another
More informationLecture D30 - Orbit Transfers
J. Peraire 16.07 Dynamics Fall 004 Version 1.1 Lecture D30 - Orbit Transfers In this lecture, we will consider how to transfer from one orbit, or trajectory, to another. One of the assumptions that we
More informationCreating Satellite Orbits
Exercises using Satellite ToolKit (STK) vivarad@ait.ac.th Creating Satellite Orbits 1. What You Will Do Create a low-earth orbit (LEO) satellite Create a medium-earth orbit (MEO) satellite Create a highly
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E S 8-9 Part B, Lectures 8-9 23 March, 2017 C O N T E N T S In this lecture, we will look at factors that cause an orbit to change over time orbital
More informationFIBER OPTIC GYRO-BASED ATTITUDE DETERMINATION FOR HIGH- PERFORMANCE TARGET TRACKING
FIBER OPTIC GYRO-BASED ATTITUDE DETERMINATION FOR HIGH- PERFORMANCE TARGET TRACKING Elias F. Solorzano University of Toronto (Space Flight Laboratory) Toronto, ON (Canada) August 10 th, 2016 30 th AIAA/USU
More informationCHRIS HYPERSPECTRAL MISSION FIVE YEARS SINCE LAUNCH
CHRIS HYPERSPECTRAL MISSION FIVE YEARS SINCE LAUNCH Paul Stephens Dr Mike Cutter Fordate Island, Indonesia Mission Overview Instrument Overview Observation Programme Applications Snapshots Future Plans
More informationSECTION 9 ORBIT DATA - LAUNCH TRAJECTORY
SECTION 9 ORBIT DATA - LAUNCH TRAJECTORY --~'- SECTION 9 LAUNCH TRAJECTORY 9.1 MISSION PROFILE IUE was launched by a three-stage Delta 2914 launch vehicle from Cape Kennedy on January 26, 1978 at l7 h
More informationCNES R&D and available software for Space Images based risk and disaster management
CNES R&D and available software for Space Images based risk and disaster management 1/21 Contributors: CNES (Centre National d Etudes Spatiales), Toulouse, France Hélène Vadon Jordi Inglada 2/21 Content
More information= o + t = ot + ½ t 2 = o + 2
Chapters 8-9 Rotational Kinematics and Dynamics Rotational motion Rotational motion refers to the motion of an object or system that spins about an axis. The axis of rotation is the line about which the
More informationAn experimental study of angular variations of brightness surface temperature for some natural surfaces
An experimental study of angular variations of brightness surface temperature for some natural surfaces Juan Cuenca, José A. Sobrino, and Guillem Soria University of Valencia, c./ Dr. Moliner 5, 46 Burjassot,
More informationReduction of the Response Time of Earth Observation Satellite Constellations using Inter-satellite Links
Reduction of the Response Time of Earth Observation Satellite Constellations using Inter-satellite Links S. De Florio Data quality of Earth observation satellites is often evaluated in terms of short system
More informationEUCLID SKY SCANNING TOOL
EUCLID SKY SCANNING TOOL 5TH INTERNATIONAL CONFERENCE ON ASTRODYNAMICS TOOLS AND TECHNIQUES ESA/ESTEC, NOORDWIJK, THE NETHERLANDS 29 MAY 1 JUNE 2012 Francesco Cacciatore (1), Mariano Sánchez (2), Alberto
More informationPrevious Lecture. Orbital maneuvers: general framework. Single-impulse maneuver: compatibility conditions
2 / 48 Previous Lecture Orbital maneuvers: general framework Single-impulse maneuver: compatibility conditions closed form expression for the impulsive velocity vector magnitude interpretation coplanar
More informationGeoreferencing. datum. projection. scale. The next few lectures will introduce you to these elements. on the Earth, you ll need to understand how
Georeferencing GOAL: To assign a location to all the features represented in our geographic information data In order to do so, we need to make use of the following elements: ellipsoid/geoid To determine
More informationQuaternion-Based Tracking Control Law Design For Tracking Mode
A. M. Elbeltagy Egyptian Armed forces Conference on small satellites. 2016 Logan, Utah, USA Paper objectives Introduction Presentation Agenda Spacecraft combined nonlinear model Proposed RW nonlinear attitude
More informationAttitude Determination and. Attitude Control
Attitude Determination and Placing the telescope in orbit is not the end of the story. It is necessary to point the telescope towards the selected targets, or to scan the selected sky area with the telescope.
More informationResults and Analyses of Debris Tracking from Mt Stromlo
Results and Analyses of Debris Tracking from Mt Stromlo Jizhang Sang 1, Ian Ritchie, Matt Pearson, and Craig Smith EOS Space Systems Pty Ltd, Australia Abstract In the last 2 years, EOS Space Systems has
More informationRAPID GEOSYNCHRONOUS TRANSFER ORBIT ASCENT PLAN GENERATION. Daniel X. Junker (1) Phone: ,
RAPID GEOSYNCHRONOUS TRANSFER ORBIT ASCENT PLAN GENERATION Daniel X. Junker (1) (1) LSE Space GmbH, Argelsrieder Feld 22, 82234 Wessling, Germany, Phone: +49 160 9111 6696, daniel.junker@lsespace.com Abstract:
More informationPhysics 201 Midterm Exam 3
Name: Date: _ Physics 201 Midterm Exam 3 Information and Instructions Student ID Number: Section Number: TA Name: Please fill in all the information above Please write and bubble your Name and Student
More informationPractice Test 3. Name: Date: ID: A. Multiple Choice Identify the choice that best completes the statement or answers the question.
Name: Date: _ Practice Test 3 Multiple Choice Identify the choice that best completes the statement or answers the question. 1. A wheel rotates about a fixed axis with an initial angular velocity of 20
More informationThe AEOLUS Mission - In Orbit Commissioning and Verification
The AEOLUS Mission - In Orbit Commissioning and Verification ADM-Aeolus CAL/VAL Workshop, ESRIN P McGoldrick, J Brewster, J Marshall, F Fabre Airbus Defence and Space, UK and France 12 th February 2015
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E 6 Part B, Lecture 6 17 March, 2017 C O N T E N T S In this lecture, we will look at various existing satellite tracking techniques. Recall that we
More informationDE-ORBITATION STUDIES AND OPERATIONS FOR SPIRALE GTO SATELLITES
DE-ORBITATION STUDIES AND OPERATIONS FOR SPIRALE GTO SATELLITES François BONAVENTURE (1), Slim LOCOCHE (2), Anne-Hélène GICQUEL (3) (1) Tel. (+33) (0)5 62 19 74 27, E-mail. francois.bonaventure@astrium.eads.net
More informationIAA-CU A Simulator for Robust Attitude Control of Cubesat Deploying Satellites
A Simulator for Robust Attitude Control of Cubesat Deploying Satellites Giovanni Mattei, George Georgiou, Angelo Pignatelli, Salvatore Monaco Abstract The paper deals with the development and testing of
More informationSatellite Position Determination of LEO Spacecraft
Satellite Position Determination of LEO Spacecraft S.C.Rathnakara ISRO Satellite Centre 5 th September, 2007 International Committee on Global Navigation Satellite Systems(ICG) Things to follow Brief background
More informationSCIAMACHY SOLAR OCCULTATION: OZONE AND NO 2 PROFILES
SCIAMACHY SOLAR OCCULTATION: OZONE AND NO 2 PROFILES Klaus Bramstedt, Astrid Bracher, Jerome Meyer, Alexej Rozanov, Heinrich Bovensmann, and John P. Burrows Inst. of Environmental Physics, University of
More informationResearch Article Study on Zero-Doppler Centroid Control for GEO SAR Ground Observation
Antennas and Propagation Article ID 549269 7 pages http://dx.doi.org/1.1155/214/549269 Research Article Study on Zero-Doppler Centroid Control for GEO SAR Ground Observation Yicheng Jiang Bin Hu Yun Zhang
More informationMAE 142 Homework #5 Due Friday, March 13, 2009
MAE 142 Homework #5 Due Friday, March 13, 2009 Please read through the entire homework set before beginning. Also, please label clearly your answers and summarize your findings as concisely as possible.
More informationAnalysis of Relative Motion of Collocated Geostationary Satellites with Geometric Constraints
Analysis of Relative Motion of Collocated Geostationary Satellites with Geometric Constraints F.J. de Bruijn, E. Gill Abstract This paper investigates the design of the relative motion orbits of geostationary
More informationSpace Surveillance with Star Trackers. Part II: Orbit Estimation
AAS -3 Space Surveillance with Star Trackers. Part II: Orbit Estimation Ossama Abdelkhalik, Daniele Mortari, and John L. Junkins Texas A&M University, College Station, Texas 7783-3 Abstract The problem
More informationF. Letterio, S. Tonetti, S. Cornara, G. Vicario presented by Mariano Sánchez Nogales
DESEO Design Engineering F. Letterio, S. Tonetti, S. Cornara, G. Vicario presented by Mariano Sánchez Nogales DEIMOS Space S.L.U., Spain - 1 - Table of Contents DESEO Overview Toolkit Heritage Software
More informationAccuracy of EOCFI Software functions in the Lib and Orbit libraries
estec European Space Research and Technology Centre Keplerlaan 1 2201 AZ Noordwijk The Netherlands T +31 (0)71 565 6565 F +31 (0)71 565 6040 www.esa.int Accuracy of EOCFI Software functions in the Lib
More informationANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS
ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.
More informationSpinning Satellites Examples. ACS: Gravity Gradient. ACS: Single Spin
Attitude Determination and Attitude Control Placing the telescope in orbit is not the end of the story. It is necessary to point the telescope towards the selected targets, or to scan the selected sky
More informationSub-millimeter size debris monitoring system with IDEA OSG 1
Sub-millimeter size debris monitoring system with IDEA OSG 1 Masahiko Uetsuhara, Mitsunobu Okada, Yasunori Yamazaki, Astroscale Pte. Ltd. Toshiya Hanada Kyushu University ABSTRACT The 20-kg class microsatellite
More informationOrbit and Transmit Characteristics of the CloudSat Cloud Profiling Radar (CPR) JPL Document No. D-29695
Orbit and Transmit Characteristics of the CloudSat Cloud Profiling Radar (CPR) JPL Document No. D-29695 Jet Propulsion Laboratory California Institute of Technology Pasadena, CA 91109 26 July 2004 Revised
More informationFigure 1. View of ALSAT-2A spacecraft
ALSAT-2A TRANSFER AND FIRST YEAR OPERATIONS M. Kameche (1), A.H. Gicquel (2), D. Joalland (3) (1) CTS/ASAL, 1 Avenue de la Palestine, BP 13, Arzew 31200 Oran, Algérie, email:mo_kameche@netcourrier.com
More informationERS-ENVISAT CROSS-INTERFEROMETRY SIGNATURES OVER DESERTS. Urs Wegmüller, Maurizio Santoro and Christian Mätzler
ERS-ENVISAT CROSS-INTERFEROMETRY SIGNATURES OVER DESERTS Urs Wegmüller, Maurizio Santoro and Christian Mätzler Gamma Remote Sensing AG, Worbstrasse 225, CH-3073 Gümligen, Switzerland, http://www.gamma-rs.ch,
More informationInnovative image geo-referencing tool for decision support in wildfire fighting
Modelling, Monitoring and Management of Forest Fires II 173 Innovative image geo-referencing tool for decision support in wildfire fighting A. Losso, L. Corgnati & G. Perona Turin Polytechnic, Italy Abstract
More informationSwimming in Newtonian Spacetime: Motion by Cyclic Changes in Body Shape
January 5, 004 Swimming in Newtonian Spacetime: Motion by yclic hanges in Body Shape Michael J. Longo University of Michigan, Ann Arbor, MI 48109 E-mail: mlongoumich.edu In a recent article in Science,
More informationClient RSK Environment Ltd Job No. J387 Date 19/11/2016. Project: Environmental Survey Location: Otranto Vessel: RV Atlante
FUGRO OCEANSISMICA S.p.A. GEODETIC PARAMETERS Client RSK Environment Ltd Job No. J387 Date 19/11/216 Project: Environmental Survey Location: Otranto Vessel: RV Atlante WGS-84 Geodetic Parameters Datum
More informationResearch Article The Impact on Geographic Location Accuracy due to Different Satellite Orbit Ephemerides
Mathematical Problems in Engineering Volume 2009, Article ID 856138, 9 pages doi:10.1155/2009/856138 Research Article The Impact on Geographic Location Accuracy due to Different Satellite Orbit Ephemerides
More informationADEN ALOS AVNIR-2 C YCLIC R EPORT C YCLIC R EPORT #29 28 J ULY 2009 TO 12 S EPTEMBER 2009
ADEN ALOS AVNIR-2 C YCLIC R EPORT C YCLIC R EPORT #29 28 J ULY 2009 TO 12 S EPTEMBER 2009 This RGB composite AVNIR-2 image has been taken on the 16th of September, 2008 from orbit 14095 over the desert
More informationExtending the Patched-Conic Approximation to the Restricted Four-Body Problem
Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean
More informationRELATIVE MISSION ANALYSIS FOR PROBA 3: SAFE ORBITS AND CAM
5TH INTERNATIONAL CONFERENCE ON SPACECRAFT FORMATION FLYING MISSIONS AND TECHNOLOGIES RELATIVE MISSION ANALYSIS FOR PROBA 3: SAFE ORBITS AND CAM Munich, 31 st May 2013 T. V. Peters D. Escorial Presented
More informationSatellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN
Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,
More informationStatus of the Gravitational Redshift Test with Eccentric Galileo Satellites
Status of the Gravitational Redshift Test with Eccentric Galileo Satellites J. Ventura-Traveset, R. Prieto-Cerdeira, D. Blonski ICG 11, Sochi, 7 November 2016 Contents 1. Galileo satellites 5 & 6 Status
More informationPhysics 201 Midterm Exam 3
Physics 201 Midterm Exam 3 Information and Instructions Student ID Number: Section Number: TA Name: Please fill in all the information above. Please write and bubble your Name and Student Id number on
More informationPhysics for Scientists and Engineers 4th Edition, 2017
A Correlation of Physics for Scientists and Engineers 4th Edition, 2017 To the AP Physics C: Mechanics Course Descriptions AP is a trademark registered and/or owned by the College Board, which was not
More informationSSTD = Standard deviation SMA = Semi Major Axis
- 1 C - EPC-95-212 NOVEL ORBT RASNG STRATEGY MAKES LOW THRUST COMMERCALLY VABLE. ARNON SPTZER* ABSTRACT A new technique for utilizing low thrust Electric Propulsion System (EPS) for orbit raising has been
More informationCURRENT STATUS OF SCIAMACHY POLARISATION MEASUREMENTS. J.M. Krijger 1 and L.G. Tilstra 2
% % CURRENT STATUS OF SCIAMACHY POLARISATION MEASUREMENTS JM Krijger 1 and LG Tilstra 2 1 SRON (National Institute for Space Research), Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands, krijger@sronnl
More informationGalileo Extended Slots Characterisation and Relation with the Nominal Constellation
Galileo Extended Slots Characterisation and Relation with the Nominal Constellation By Andrés AYALA, ) Rubén CASTRO, 2) Nityaporn SIRIKAN, 3) Daniel BLONSKI, 3) and Daniel NAVARRO 3) ) ESTEC, European
More informationSolar Orbiter Ballistic Transfer Mission Analysis Synthesis
European Space Agency Agence Spatiale Européenne directorate of operations and infrastructure ground systems engineering department mission analysis office MAO Working Paper No. 483 Issue 1, Rev. 0 Solar
More informationRosetta Optical Navigation
Rosetta Optical Navigation Mathias Lauer ESA / ESOC / Flight Dynamics Page 1 Overview Rosetta: Mission and Spacecraft Asteroid Flybys (Steins): - Scenario - Navigation Strategy - Image Processing - Autonomous
More informationPHYS 1303 Final Exam Example Questions
PHYS 1303 Final Exam Example Questions 1.Which quantity can be converted from the English system to the metric system by the conversion factor 5280 mi f 12 f in 2.54 cm 1 in 1 m 100 cm 1 3600 h? s a. feet
More informationEstimating Geometric Aspects of Relative Satellite Motion Using Angles-Only Measurements
Estimating Geometric Aspects of Relative Satellite Motion Using Angles-Only Measurements Jason Schmidt 1 Dept of Mechanical & Aerospace Engineering, Utah State University, Logan, UT, 84322 Thomas Alan
More informationHomework 2; due on Thursday, October 4 PY 502, Computational Physics, Fall 2018
Homework 2; due on Thursday, October 4 PY 502, Computational Physics, Fall 2018 Department of Physics, Boston University Instructor: Anders Sandvik PERTURBATION OF A GEO-STATIONARY ORBIT Here you will
More information