Haversine Formula:
http://www.movable-type.co.uk/scripts/gis-faq-5.1.html
My question is slightly different. The mean radius of the Earth is approx.
6,372 km, but what formula is used to calculate the radius of the Earth
for any given latitude?
I wish my trig was up to working this out but I'm just a humble computer
scientist with no degree in math.
Thank you all for any help in this.
R * cos(latitude)
As I understand your question.
Thanks, but I think I explained my question badly; I'll re-phrase it.
For any specific latitude, what is the formula for calculating the radius
of the Earth at that latitude.
So if I wish to calculate the radius of the Earth at the latitude of
N56:29:35 (56.492956), what formula can I use?
Thanks.
6372 * cos(56.492956/180*pi) = 3518 km.
As you move toward one of the poles, this will finally become 0.
Or do you have Earth flattening in mind? If not, the answer can also be that
the radius doens't change.
http://www.usenet-replayer.com/faq/comp.infosystems.gis.html
In article <tPooj.326$cJ...@text.news.blueyonder.co.uk>,
Actually, your approach calculates the radius of a circle of latitude.
The radius of curvature of the Earth is a bit more complicated.
This issue is discussed in rather excruciating detail in the FAQ 5.1
article I mentioned in my previous post. (Look for the section labeled
5.1a which is just short of half way through the article.) For
convenience, and because that post could easily get separated from this
one, here it is again:
http://www.usenet-replayer.com/faq/comp.infosystems.gis.html
- Bob
Thanks for the pointer Bob - in both posts.
Cheers.
eventhough your question might look clear, it is in fact very vague. If
you are a seaman, you are looking for a seaman answer.
If your are a geodesist, you are looking for a geodesist answer.
If you are a mathematician, you are looking for a mathematical answer.
Here is a possible answer
The geoid:
The hypothetical surface of the Earth that coincides everywhere with
mean sea level and is perpendicular, at every point, to the direction of
gravity. The geoid is used as a reference surface for astronomical
measurements and for the accurate measurement of elevations on the
Earth's surface.
From that point of view, there is no formula to calculate the radius of
Earth, just readings from surveys.
The reference ellipsoid:
In geodesy, a reference ellipsoid is a mathematically-defined surface
that approximates the geoid, the truer figure of the Earth, or other
planetary body. Because of their relative simplicity, reference
ellipsoids are used as a preferred surface on which geodetic network
computations are performed and point coordinates such as latitude,
longitude, and elevation are defined.
If you take WGS84 as a reference ellipsoid, the Earth is flattened by
1/298.257
Then you can calculate the radius for each latitude.
Seamen use the great circle for long planning routes. They assume,
without large errors, that the Earth is spherical.
As you see, a question can rise many more questions
Jacquelin Hardy
a seaman that likes maths.
Uffe Kousgaard a écrit :
Matt asked the question.
Regards
Uffe
And I think Uffe knows the answers!
Jacquelin is quite correct. I haven't bothered getting involved in
this discussion because as Jacquelin says, the answer is "it depends
on your application". The original poster was being confused by
attempting to compute distances in a projected coordinate system,
which is a necessarily inaccurate way of carrying out the measurement
(I know, it IS accurate if you are doing distance from the pole of
projection of an equidistant projection, but typically that isn't the
case, and wasn't the case for the OP). In fact I suspect that the
simple cosine rule formula was quite accurate enough for his or her
needs; the difference between spherical geometry, ellipsoidal geometry
and surveyed distance on a geoid is only likely to bother a surveyor
looking for closure in a survey. Over long distances, the difference
is less than a kilometre - and that's for distances between Antarctica
and Europe.
A practical solution: if you want a good geodetic distance calcluator,
then use geod which is bundled with proj. That is a) robust b) written
by an expert in the field of geodesy (Gerry Evenden), c) accurate and
d) free! As far as I'm concerned, this is the gold standard for great
circle calculations on the ellipsoid; it avoids all the
computationally nasty things that can happen if you roll your own! And
as Jacquelin has pointed out, you can't compute distance on the geoid
analytically - but I think the difference would be less than a
millimetre! It is worth noting that the maximum vertical distance
between sphere and ellipsoid is about 11km. The maximum vertical
distance between ellipsoid and geoid is +85 metres to -106 metres);
the ellipsoid is a fantastically accurate depiction of the shape of
the Earth.
Paul
Uffe Kousgaard wrote:
-------------------------------------
-----------------------------------------------------------------------------
function Radius(Latitude : Degrees;
Spheroid : Spheroid_Type) return Meters is
e2 : long_float renames Models(Spheroid).Eccentricity_Squared;
b : Meters renames Models(Spheroid).Semi_Minor_Axis;
f : long_float renames Models(Spheroid).Flattening;
boa : long_float := 1.0 - f;
Phi : Radians := Deg_to_Rad(Latitude);
Psi : Radians := Geocentric_Latitude(Phi, boa);
cl : long_float;
r : Meters;
begin
cl := cos(Psi);
r := b/sqrt(1.0 - e2*cl**2);
return r;
end Radius;
##-----------------------------------------------##
Delivered via http://www.psjournal.com/
GPS/GIS Community
Web and RSS access to your favorite newsgroup -
comp.infosystems.gis - 3593 messages and counting!
##-----------------------------------------------##