| Radii on the Ellipsoid
See Also, GIS Tips Archives

Exaggerated illustration showing the plane in which
the transverse radius of curvature is measured.
Last time we discussed the important role of the radius in converting from the angular
units of latitude and longitude (e.g. degrees) to the linear units of X and Y coordinates
of the cartesian coordinate system (e.g. meters or feet). However, the discussion was
limited to the sphere, and for good reason. In a sphere, there is a single radius one has
to deal with. We know from previous articles that the ellipsoid has at least two radii:
the polar and equatorial. It is reasonable to ask, then, which of the two radii should the
Casual Cartographer use.
In the following, we will present formulas for various radii. In the formulas:

Approximating Spheres
The idea of using the average between the two would be a quite reasonable approach.
That is, lets approximate the ellipsoid with a sphere whose radius is some sort of average
of the two radii that we have. There are many variations to this basic approach.
Arithmetic Mean
What most of us would consider the average is technically known as the arithmetic mean;
i.e. add em up and divide by 2. That is:

This produces a reasonable value which can be, but is seldom, used.
Geometric Mean
If one takes a moment to examine the following formula for the geometric mean, one can
see that it is yet another, rather interesting, way of taking an average:

This "average" is referred to as the geometric mean and the resulting radius
is used quite often in cartography; especially when the double projection technique is
used.
Sphere of Equal Volume
Serious cartographers will, at times, choose to use a radius which produces a sphere
whose total volume is the same as the ellipsoid in use. In this case, the radius would be:

Sphere of Equal Area
Then again, if one is seeking an equal area solution, one might choose a radius such
that the surface area of the resulting sphere is the same as the surface area of the
original ellipsoid. In this case:

This value is often used, in conjunction with the double projection approach, to
produce a projection of the equal area variety.
EquiDistant Sphere
Then again, if one is seeking an equidistant solution, one could choose a radius such
that the length of the meridians on the resulting sphere have the same length as the
meridians on the original ellipsoid. In this case:
where 
So, in addition to the original two radii, we now have another five. Had enough? These
are just some of the radii which can be used to approximate the ellipsoid. However, all of
the above yield only a sphere which approximates the size of the ellipsoid. In many cases,
the resulting approximation is sufficiently accurate for the objectives of a particular
map or application. However, the actual number of radii associated with a specific
ellipsoid is limitless, i.e. infinite.
Radius of Curvature
We arrived at this discussion of radii on the ellipsoid looking for a solution to the
basic problem of converting angles to linear measures. We noted that the radius is the
value which can be used to convert angles to linear measures. Since the radius on the
sphere was always the same, or constant, we could use this relationship for both small and
large angles.
What do we do when the radius is not constant, as is the case of the ellipsoid? The
calculus provides us with the answer; but well restrict ourselves to Casual
Cartographer level terminology to explore it. To cross this bridge, we need only recognize
that at any given point on the ellipsoid, there is a "magic radius" which
will do the job at that point. That is, at a specific point on the ellipsoid, there is a
specific value which can be used as the radius which will accurately convert small changes
of angle to small changes in linear distances. (Those who have made it through Calculus I
will recognize this as the familiar theme throughout differential calculus.) This magic
radius is referred to as the radius of curvature. In Casual Cartographer terms, we define
the radius of curvature as "the radius of the specific circle which fits the
curvature of the ellipsoid better than any other circle at the specific point of
interest."
Due to the ellipsoid being an ellipse of revolution, it is clear that the radius of
curvature of a specific point depends solely on the latitude of the point. That is, all
points with the same latitude will have the same radius of curvature; regardless of the
longitude. Even so, there are an infinite number of "specific points" which can
be defined between the equator and the north pole. As we shall see, the radius of
curvature at each such point will be different. There are, therefore, an infinite number
of "radii of curvature" on the ellipsoid.
As a result, if one is mapping a small region localized at, for example, 40 degrees of
latitude, one could indeed use the radius of curvature of the ellipsoid at 40 degrees of
latitude. The results would, one could reasonably expect, would be better than any of the
other radii which we discussed above.
Once again, what appears simple and straight forward is not quite as simple. Note that
we defined our radius of curvature as the radius of a circle which best fits the ellipsoid
at a given point. We did not say the sphere which best fits the ellipsoid at the
specific point. There are some good reasons for this.
Perhaps this because the nature of determining the radius of the "best
fitting" sphere is unnecessarily complex. I must admit Ive never seen the
formulae for such a calculation. Perhaps that is proof enough.
More likely, however, is the fact that as Casual Cartographers, we need to make two
distinct conversions of angles to linear measure. That is, we need to convert degrees of
longitude to linear units (e.g. meters of feet) of X, and degrees of latitude to linear
units of Y. Also important to note is that in the cartesian X and Y system, the distances
we come up with for X will be in a direction normal to (i.e. perpendicular to) the
direction in which the Y distances apply. Thus, two different radii of curvature can be of
value. One is used to converting longitude and one for converting latitude, each of which
would produce the best results possible. Thus, I suggest that by using two separate radii,
the results are better than if used just one. As we shall see, these two radii of
curvature are not the same.
Meridional Radius of Curvature
For converting degrees of latitude to linear units of Y, the radius used is known as
the meridional radius of curvature. That is, the radius of curvature at a specific point
in the north/south direction. Conceptually, it is the simpler of the two radii involved,
and it is, essentially, the radius of curvature, at a specific latitude, of the ellipse
used to generate the ellipsoid. The meridional radius of curvature can be calculated from
the following formula:

Notice that the effect of the ellipsoid is introduced by the eccentricity squared and
the latitude of the point of interest. That is, the actual radius of curvature is not only
varies from one ellipsoid to the next, but is also dependent upon the latitude of the
actual point of interest.
The meridional radius of curvature can be used to faithfully convert a small change of
latitude into a corresponding change of our Y coordinate. One cannot use the meridional
radius of curvature to compute the distance on the ellipsoid from the equator to an
arbitrary latitude. Since the radius of curvature changes as the latitude changes, which
value of latitude would you use to compute the radius of curvature? This type of
calculation requires the use of the Calculus, and well deal with that some other
time. However, for a small change in latitude, the meridional radius of curvature,
determined at the specific latitude of interest, can be used to calculate the resulting
change in Y. Hopefully you can see that the smaller the change in latitude, the more
faithful the calculated change in Y will be. (That calculus stuff again.)
Transverse Radius of Curvature
The second radius, the one which will convert a small change in longitude to an
equivalent change in X is known as the transverse radius of curvature of an ellipsoid. It
is also known as the radius of the first vertical section. In any case it is more
difficult to grasp, conceptually, as it is defined as the radius of curvature in the plane
which is normal to (i.e. perpendicular to) both the surface of the ellipsoid at, and the
meridian passing through, the specific point of interest. This should not be surprising as
we are converting longitude to X, and the X axis is normal to (i.e. perpendicular to) the
Y axis. The transverse radius of curvature can be computed from the following formula:

As with the meridional radius of curvature, the transverse radius of curvature not only
varies from one ellipsoid to the next, but also with the latitude of the point of
interest.
In our previous tip, we noted that the radius of a small circle of latitude on the
sphere could be computed as:

On the ellipsoid, therefore, the radius of a small circle of latitude is:

or, replacing with the above definition, the radius of a small
circle of latitude on the ellipsoid becomes:

This radius can be used to convert both small and large changes in longitude to
appropriate changes in X, assuming that the latitude does not change, as the radius
calculation does not depend on longitude. Thus, the value of the transverse radius of
curvature becomes evident. But notice that the resulting expression has an even more
complex dependency on the latitude.
Points of Interest
It should not be a big surprise to anyone that the radii of curvature on the ellipsoid
varies from one ellipsoid to the next. It may be new to some Casual Cartographers to
observe that the radii of curvature is also dependent on the latitude of the specific
point in question. That is, to perform the conversion of angles of latitude and longitude
to linear Xs and Ys, the best radius to use varies with the latitude.
Some additional points of interest regarding the meridional and transverse radii of
curvature are of note:
1) At either pole, that is when the sine of the latitude is either +1 or 1, both
radii of curvature are the same, specifically 
2) Only at the poles will the two radii of curvature be the same.
3) Both radii of curvature achieve their maximum values at the poles.
4) At the poles, both radii of curvature are not only greater than the polar radius,
but is also greater than the equatorial radius!!!
5) At the equator, both radii achieve their minimum values.
6) At the equator, the transverse radii of curvature is equal to the equatorial radius.
7) The meridional radius of curvature equals the equatorial radius at about 55 degrees
of latitude.
8) The transverse radius of curvature equals the polar radius at about 35 degrees of
latitude.
Making a Projection
As we saw last time, on the sphere, we have a single radius which can be used to
convert the angular measures (e.g. degrees) of the spherical coordinate system to the
linear units of the cartesian system (e.g. meters of feet). In the ellipsoidal case, we
can approximate the ellipsoid with a suitably sized sphere. However, to achieve the best
results possible in the ellipsoidal case, we need not only deal with two different radii,
but both radii will change significantly depending upon the latitude of the specific
features being mapped. Thus, as we will learn in future tips, the ellipsoidal form of a
projection is often very much more complex than the corresponding spherical form.
In the last several tips, we have discussed coordinate systems, both spherical and
cartesian. In this and the previous tip, we discussed the means of converting from the
angular units of the spherical coordinate systems to the linear units of the cartesian
system. Perhaps you can see where this is leading us. Next time, well use the
results of these discussions to develop some of the simpler projections.
Top | More GIS Tips Articles | Back
to The Casual Cartographer August 1999 |