15
$\begingroup$

What I'm asking for

I stress that I am not asking for the formula---I know the formula, and how to derive it. Several different versions of it are reproduced near the end of the post. In fact, someone else has not only derived it as well, but also nicely presented one of the derivations here.

What I need is a reputable source for the formula so that, for example, one could put it on Wikipedia without violating its ban on reporting original research. [People have actually tried... But the relevant article has some very conscientious editor who deleted the section on the grounds of it being original research... and, well, unfortunately, the editor is correct, so there's not much point in trying to fight it.]

The reason I'm posting in Computer Graphics stackexchange

Since someone here might have modeled the way the Earth looks form orbit, perhaps he or she might know if this formula (or, more likely, some generalization of it) is published in some book, or journal, or conference proceedings, or class notes, etc.

I've done the "due googling"

Please understand that I'm not asking anyone to go searching for the answer on my behalf. I've done lots of googling already, and am only posting here as a last resort. My (far-fetched) hope is that someone here will simply know a reference right off the bat; if not... well, I hope at least you enjoyed the pretty picture below (if I do say so myself, with full awareness I'm talking to people interested in computer graphics of all things) before you moved on to bigger and better things.

Two sources that come close

  1. D. K. Lynch, "Visually discerning the curvature of the Earth," Applied Optics vol. 47, H39 (2008). It is freely available here. Unfortunately, instead of doing it the right way (which is not that hard), the author opted for a hack, which (a) I don't completely understand, and (b) which doesn't agree with what I know to be the correct formula.

  2. R. Hartley and A. Zisserman, Multiple View Geometry in Computer Vision, 2nd ed. (Cambridge University Press, Cambridge UK, 2004). In Sec. 8.3, "Action of a projective camera on quadrics," we read:

Suppose the quadric is a sphere, then the cone of rays between the camera centre and quadric is right-circular, i.e. the contour generator is a circle, with the plane of the circle orthogonal to the line joining the camera and sphere centres. This can be seen from the rotational symmetry of the geometry about this line. The image of the sphere is obtained by intersecting the cone with the image plane. It is clear that this is a classical conic section, so that the apparent contour of a sphere is a conic.

In principle, this would be exactly what is needed, if only just a bit more information were included---at least an expression for the eccentricity of the conic as a function of the distance to the sphere and the sphere radius (in the case when the image plane is perpendicular to a generatrix of the cone, as is the case when the pinhole camera is directed at a point on the horizon).

Details about the formula for which I need a scholarly reference

We assume a perfectly spherical, perfectly smooth Earth with no atmosphere. We point an idealized pinhole camera at the horizon, and, using straightforward central projection, compute the shape of the image of the horizon on the back of the camera (i.e. the shape it will have on film---the "film plane"). Here's a graphic (made in Asymptote, for those interested) that should make this clearer:

enter image description here

As we saw above, the image of the horizon is a portion of a conic section. Let $\varepsilon$ be the eccentricity of the conic; the derivation I mentioned above instead uses a parameter $k$, which is just the inverse eccentricity: $k=1/\varepsilon$. The eccentricity itself is given as $\varepsilon=1/\sqrt{\epsilon(2+\epsilon)}$, where $\epsilon=h/R$ is the ratio of the altitude $h$ of the pinhole above the surface of the Earth and the Earth radius $R$. [Instead of using $\epsilon$, which is the ratio of the altitude to $R$, it may be useful to use $\eta$, the ratio of the pinhole's distance to the Earth's center, $h+R$, to the Earth's radius: $\eta=(R+h)/R=1+\epsilon$. In terms of $\eta$, we have $\varepsilon=1/\sqrt{\eta^{2}-1}$.]

The distance from the pinhole (point $P$ in the graphic) to the film plane is taken to be one unit length.

The $y$-axis in the film plane is chosen to be parallel to the line joining the center of the Earth $C$ (not shown in the image) and the point on the horizon (labeled $V$ in the image) at which the camera is trained. This choice is well-defined because the line $CV$ must be parallel to the film plane. The reason for this is that both $CV$ and the film plane are perpendicular to the line of sight $PV$ (the line joining $P$ and $V$). And that is because 1. the line $PV$ is tangent to the Earth at $V$, thus perpendicular to $CV$, and 2. $PV$ is perpendicular to the film plane because the camera is trained at $V$. The $x$ axis is of course perpendicular to the $y$ axis and lies in the film plane, and the origin is chosen as the projection of the point $V$.

With these definitions out of the way, we are ready to write down a representation of the conic section that is the image of the Earth's horizon. This can be written in many ways, some of which are given below. What I need is a reputable reference for any one of these formulas, or for a formula equivalent to them.

1. The explicit formula given in the derivation mentioned above

The derivation I mentioned above gives this as the final version:

$[\,y\,(1/\varepsilon-\varepsilon)-1\,]^{2}+x^{2}(1/\varepsilon^{2}-1)=1.$

Let's represent this in a couple of additional ways.

2. Expression in terms of the canonical equation of a conic section

In this case, the equation takes the following form:

$x^2=2\mu y-(1-\varepsilon^{2}) y^2$,

where, in our case, $\mu=\varepsilon$.

The advantage of the canonical form is that it can deal with all conics on an equal footing, including in particular the case of the parabola, $\varepsilon=1$. In the ``standard'' formulation (see below), the case of the parabola can only be dealt with by taking the limit $\varepsilon\to 1$.

Details: the above formula holds in the case of a right circular cone, whose sides subtend an angle of $2\theta$, being intersected---at a distance $d$ from the vertex of the cone---by a plane at an angle $\omega$ relative to the cone axis. (To clarify: $d$ is the distance from the cone vertex to the point on the ellipse that is the closest to the cone vertex; that point is always one of the ends of the major axis of the ellipse). In this general case, the eccentricity is given as $\varepsilon=\cos \omega/\cos \theta$, while $\mu=d(\varepsilon-\cos|\omega+\theta|)$.

In terms of the above graphic: $d$ is the distance from $P$ to the film plane (i.e., the distance along the dotted red line); $\theta$ is the angle between the dotted red line and the axis of the cone (which is the line joining $P$ and the center of the Earth---the extension of the black line labeled $h$ in the graphic); the angle $\omega$ is the angle between the axis of the cone and the film plane.

Given that the film plane is perpendicular to the dotted red line, we have $\omega+\theta=\pi/2$; in addition, we take $d=1$, which then together give that $\mu=\varepsilon$.

3. Expression in terms of the ``standard form' of a conic section

This form is perhaps the most familiar:

$\frac{(x-x_{0})^{2}}{p^2}+\frac{(y-y_{0})^{2}}{q^2}=1$.

It is related to the parameters entering the canonical equation (see 2., above) as follows:

$x_{0}=0$;

$y_{0}=q=\frac{\mu}{1-\varepsilon^{2}}$ (which is $\frac{\varepsilon}{1-\varepsilon^{2}}$ in our case—note that $y_{0}=q$ follows from the fact that the ellipse passes through the origin); and

$p=q\smash[b]{\sqrt{|1-\varepsilon^{2}|}}=\frac{\mu}{\smash[b]{\sqrt{|1-\varepsilon^{2}|}}}$ (which is $\frac{\varepsilon}{\smash[b]{\sqrt{|1-\varepsilon^{2}|}}}$ in our case).

It is obvious that the parabolic case, $\varepsilon=1$, will create problems; as mentioned above, that case must be dealt with through taking the limit $\varepsilon\to 1$.

4. Expression in terms of a parametric curve

$x=\frac{(\epsilon+1) \cos (\alpha )}{\sin (\alpha )+\epsilon(\epsilon+2)}$

$y=\frac{\sqrt{\epsilon (\epsilon+2)} (\sin (\alpha )-1)}{\sin (\alpha )+\epsilon (\epsilon+2)},$

where $\alpha$ is the longitude of a point on the horizon, defined so that $\alpha=\pi/2$ corresponds to the point $V$ in the image above (i.e. to the point at which the pinhole camera is trained).

For how one might use these formulas, see this.

In conclusion...

Has anyone seen the formulas above in some reputable source, possibly in the context of modeling how the Earth looks from space? If so, could you let me know what this source was?

Thanks!

$\endgroup$
  • $\begingroup$ @trichoplax 1. Why now, after all this time? 2. Moreover, the response to the meta-question, while not firm, was leaning towards allowing the question. 3. Finally, as I have explained there, the argument about "opinionated answers and spam" is completely inapplicable in the case of requests for reputable sources for a specific fact. Either there is such a source or not. $\endgroup$ – linguisticturn Mar 13 '16 at 20:29
  • $\begingroup$ If anyone disagrees that off site resource requests should be off topic, they can have their say in Are questions asking for off site resources on topic? If anyone thinks requests for reputable sources should be an exception to the off site resources rule, they can have their say in Are requests for reputable sources on topic? $\endgroup$ – trichoplax Mar 14 '16 at 0:10
  • 1
    $\begingroup$ @trichoplax You yourself were uncertain whether this question is really a request for off-site resources, and you say that uncertainty remains. Is it a standard practice on this stackexchange that, when there is uncertainty about whether a question is off-topic, to err on the side of it being off-topic? $\endgroup$ – linguisticturn Mar 20 '16 at 14:23
  • $\begingroup$ You make a good point, which I've taken some time to think about. In the absence of any community support for excluding such questions, I have reopened this one. $\endgroup$ – trichoplax Mar 23 '16 at 15:34
  • $\begingroup$ This is allmost certainly never going to be answered. Besides what does it have to do with 3d graphics. $\endgroup$ – joojaa Mar 24 '16 at 5:29
3
$\begingroup$

The curve you are seeking is just the intersection of a plane (the back of the camera) and a right circular cone. This is not really a question about the earth, or views of planets from space; it's just plain simple 3D coordinate geometry. To find a reference, I'd recommend searching for "intersection of a plane and a cone" or "plane section of a cone", or "plane section of quadric", something like that.

I'd expect that you can find the relevant formulae (and derivations) in any standard text on 3D coordinate geometry. Some likely places are:

  • Salmon -- A Treatise on the Analytic Geometry of Three Dimensions
  • Somerville -- Analytical Geometry of Three Dimensions
  • Snyder and Sisam -- Analytic Geometry of Space

These are all pretty old books, and you may have trouble finding them.

You might also try asking on Math.StackExchange.

Calling the derivation "original research" seems absurd, to me. It's a high-school homework problem in analytic geometry.

$\endgroup$
  • $\begingroup$ Thank you for the answer! I will try to locate these sources. As far as your statements that this is not really original research or even about the Earth: to see why as far as Wikipedia goes both of these are disputable, see here and also e.g. here. Many Wikipedia editors will ultimately agree with you, but some will create problems. The easiest way to deal with the latter ones is to show an appropriate source. $\endgroup$ – linguisticturn Mar 25 '18 at 18:47

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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