Subject: | |
From: | |
Reply To: | |
Date: | Thu, 16 May 1996 12:42:07 EDT |
Content-Type: | text/plain |
Parts/Attachments: |
|
|
----------------------------Original message----------------------------
Distance computations between pairs of coordinates are complicated by the fact
that the earth is not flat. A good approximation is found in
47 CFR Sec. 73.208 which is used by radio engineers.
Sec. 73.208 Reference points and distance computations.
(c) The method given in this paragraph shall be used to compute
the distance between two reference points, except that, for
computation of distance involving stations in Canada and Mexico,
the method for distance computation specified in the applicable
international agreement shall be used instead. The method set
forth in this paragraph is valid only for distances not exceeding
475 km (295 miles).
(1) Convert the latitudes and longitudes of each reference point
from degree-minute-second format to degree-decimal format by
dividing minutes by 60 and seconds by 3600, then adding the results
to degrees.
(2) Calculate the middle latitude between the two reference
points by averaging the two latitudes as follows:
ML=(LAT1 (INFERIOR dd)+LAT2 (INFERIOR dd)) () 2
(3) Calculate the number of kilometers per degree latitude
difference for the middle latitude calculated in paragraph (c)(2)
as follows:
KPD (INFERIOR lat)=111.13209-0.56605 cos(2ML)+0.00120 cos(4ML)
(4) Calculate the number of kilometers per degree longitude
difference for the middle latitude calculated in paragraph (c)(2)
as follows:
KPD (INFERIOR lon)=111.41513 cos(ML)-0.09455 cos(3ML)+0.00012
cos(5ML)
(5) Calculate the North-South distance in kilometers as follows:
NS=KPD (INFERIOR lat) (LAT1 (INFERIOR dd)-LAT2 (INFERIOR dd))
(6) Calculate the East-West distance in kilometers as follows:
EW=KPD (INFERIOR lon) (LON1 (INFERIOR dd)-LON2 (INFERIOR dd))
(7) Calculate the distance between the two reference points by
taking the square root of the sum of the squares of the East-West
and North-South distances as follows:
DIST=(NS (SUPERSCRIPT 2+)EW (SUPERSCRIPT 2)) (SUPERSCRIPT 0.5)
(8) Round the distance to the nearest kilometer.
(9) Terms used in this section are defined as follows:
(i) LAT1 (INFERIOR dd) and LON1 (INFERIOR dd)=the coordinates of
the first reference point in degree-decimal format.
(ii) LAT2 (INFERIOR dd) and LON2 (INFERIOR dd)=the coordinates of
the second reference point in degree-decimal format.
(iii) ML=the middle latitude in degree-decimal format.
(iv) KPD (INFERIOR lat)=the number of kilometers per degree of
latitude at a given middle latitude.
(v) KPD (INFERIOR lon)=the number of kilometers per degree of
longitude at a given middle latitude.
(vi) NS=the North-South distance in kilometers.
(vii) EW=the East-West distance in kilometers.
(viii) DIST=the distance between the two reference points, in
kilometers.
Simple, eh?
I have placed a Q-BASIC Text file, which is simple txt that will run on
Q-BASIC with this very formula on our server, which you may access
at:
http://www.jmu.edu/wmra/dandb.txt
enjoy.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Bill Fawcett [log in to unmask]
James Madison University adm_wfawcett@jmuvax (bitnet)
821 S. Main Street, Harrisonburg, VA 22807 (703) 568-3809 voice
N. Lat. 38-26-27 W. Long. 78-52-33 (703) 568-3814 fax
- - Virginia- The home of Matthew Fontaine Maury - -
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
formula is found on
|
|
|