Is there an established algorithm to determine the subset of the orbital parameters that will place an earth satellite in a given cone during some time window?

Say I have a half-cone with its apex on the surface of the earth and opening outwards towards space. Given some range of time $[t_0,t_1]$, I'd like to determine the subset of orbital parameters that would place a vehicle in this cone during this window.

I should add that I am interested in initial orbit parameters for two-body system and that the finer points such as precession and nutation can be ignored.

  • $\begingroup$ In theory, one could do this with equations, but it would really end up quite complex. In practice, it's usually easier to figure out where the satellite is, and if it's visible from the point on Earth. $\endgroup$
    – PearsonArtPhoto
    Jan 23, 2014 at 16:23
  • $\begingroup$ @PearsonArtPhoto But that's a fundamentally different problem. What I am interested in is not whether a vehicle with a given set parameters is visible; rather, given that a vehicle is visible, what is the subset of orbital parameters that it must belong to. $\endgroup$ Jan 23, 2014 at 17:22
  • $\begingroup$ What do you mean "half cone"? $\endgroup$
    – Mark Adler
    Jan 26, 2014 at 5:00
  • $\begingroup$ @MarkAdler It's a normal cone; I just used the term half cone so that it would not be interpreted as a double cone. $\endgroup$ Jan 26, 2014 at 18:15

2 Answers 2


Your best bet is to start with Gauss's method of orbit determination, which allows you to determine a two body orbit from three angular observations. Unfortunately, I don't have a good online reference handy, but it is discussed in Fundamentals of Astrodynamics by Bate, Mueller, and White.


no i did not find any algorithm. but their is this article on Wikipedia which shows that we don't have some kind of algorithm, instead the half cone that you mentioned has a different value at different position.

" the orbital state vectors (en.wikipedia.org/wiki/Orbital_state_vectors) of an orbit are cartesian vectors of position and velocity that together with their time (epoch ) uniquely determine the trajectory of the orbiting body in space. The body does not actually have to be in orbit for its state vector to determine its trajectory; it only has to move ballistically , i.e., solely under the effects of its own inertia and gravity.

Because even satellites in low earth orbit experience significant perturbations (primarily from the non-spherical shape of the earth), the Keplerian elements computed from the state vector at any moment are only valid at that time. Such element sets are known as osculating elements (en.wikipedia.org/wiki/Osculating_elements) because they coincide with the actual orbit only at that moment.

An osculating orbit and the object's position upon it can be fully described by the six standard Keplerian orbital elements (osculating elements), which are easy to calculate as long as one knows the object's position and velocity relative to the central body. The osculating elements would remain constant in the absence of perturbations . However, real astronomical orbits experience perturbations that cause the osculating elements to evolve, sometimes very quickly. In cases where general celestial mechanical analyses of the motion have been carried out, the orbit can be described by a set of mean elements with secular and periodic terms.

Put in more general terms, a perturbed trajectory can be analysed as if assembled of points, each of which is contributed by a curve out of a sequence of curves. Variables parameterising the curves within this family can be called orbital elements . Typically (though not necessarily), these curves are chosen as Keplerian conics, all of which share one focus. In most situations, it is convenient to set each of these curves tangent to the trajectory at the point of intersection. Curves that obey this condition (and also the further condition that they have the same curvature at the point of tangency as would be produced by the object's gravity towards the central body in the absence of perturbing forces) are called osculating, while the variables parameterising these curves are called osculating elements. In some situations, description of orbital motion can be simplified and approximated by choosing orbital elements that are not osculating. Also, in some situations, the standard (Lagrange-type or Delaunay-type) equations furnish orbital elements that turn out to be non-osculating."

My answer : So, as per the above statements you have to make many no. of cones for many instants. Hence you cannot have a single algoritm for complete orbit but at each instant you can apply some formula. I did not calculate it, you can. And may be those computer programmings for orbital calculations also use this method.

plz look these articles on wikipedia, then you will understand better, i have mentioned the links.

  • $\begingroup$ Thanks, but this is not I am looking for. A two-body model would be just fine. $\endgroup$ Jan 25, 2014 at 21:22
  • $\begingroup$ yes use the two body model. $\endgroup$
    – Niket
    Jan 26, 2014 at 9:13

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.