zoukankan      html  css  js  c++  java

  • Geographical distance

    Introduction

    Calculating the distance between geographical coordinates is based on some level of abstraction; it does not provide an exact distance, which is unattainable if one attempted to account for every irregularity in the surface of the earth.[1] Common abstractions for the surface between two geographic points are:

    • Flat surface;
    • Spherical surface;
    • Ellipsoidal surface.

    All abstractions above ignore changes in elevation. Calculation of distances which account for changes in elevation relative to the idealized surface are not discussed in this article.

    Nomenclature

    Distance, D,\,! is calculated between two points, P_1\,! and P_2\,!. The geographical coordinates of the two points, as (latitude, longitude) pairs, are (phi_1,lambda_1)\,! and (phi_2,lambda_2),\,! respectively. Which of the two points is designated as P_1\,! is not important for the calculation of distance.

    Latitude and longitude coordinates on maps are usually expressed in degrees. In the given forms of the formulae below, one or more values must be expressed in the specified units to obtain the correct result. Where geographic coordinates are used as the argument of a trigonometric function, the values may be expressed in any angular units compatible with the method used to determine the value of the trigonometric function. Many electronic calculators allow calculations of trigonometric functions in either degrees or radians. The calculator mode must be compatible with the units used for geometric coordinates.

    Differences in latitude and longitude are labeled and calculated as follows:

    egin{align} Deltaphi&=phi_2-phi_1;\ Deltalambda&=lambda_2-lambda_1. end{align} \,!

    It is not important whether the result is positive or negative when used in the formulae below.

    "Mean latitude" is labeled and calculated as follows:

    phi_m=frac{phi_1+phi_2}{2}.\,!

    Colatitude is labeled and calculated as follows:

    For latitudes expressed in radians:
    heta=frac{pi}{2}-phi;\,!
    For latitudes expressed in degrees:
    heta=90^circ-phi.\,!

    Unless specified otherwise, the radius of the earth for the calculations below is:

    R\,! = 6,371.009 kilometers = 3,958.761 statute miles = 3,440.069 nautical miles.

    D_\,! = Distance between the two points, as measured along the surface of the earth and in the same units as the value used for radius unless specified otherwise.

    Singularities and discontinuity of latitude/longitude

    Longitude has singularities at the Poles (longitude is undefined) and a discontinuity at the ±180° meridian. Also, planar projections of the circles of constant latitude are highly curved near the Poles. Hence, the above equations for delta latitude/longitude (Deltaphi!, Deltalambda!) and mean latitude (phi_m!) may not give the expected answer for positions near the Poles or the ±180° meridian. Consider e.g. the value of Deltalambda! (“east displacement”) when lambda_1! and lambda_2! are on either side of the ±180° meridian, or the value of phi_m! (“mean latitude”) for the two positions (phi_1!=89°, lambda_1!=45°) and (phi_2!=89°, lambda_2!=−135°).

    If a calculation based on latitude/longitude should be valid for all Earth positions, it should be verified that the discontinuity and the Poles are handled correctly. Another solution is to use n-vector instead of latitude/longitude, since this representation does not have discontinuities or singularities.

    Flat-surface formulae

    A planar approximation for the surface of the earth may be useful over small distances. The accuracy of distance calculations using this approximation become increasingly inaccurate as:

    • The separation between the points becomes greater;
    • A point becomes closer to a geographic pole.

    The shortest distance between two points in plane is a straight line. The Pythagorean theorem is used to calculate the distance between points in a plane.

    Even over short distances, the accuracy of geographic distance calculations which assume a flat Earth depend on the method by which the latitude and longitude coordinates have been projected onto the plane. The projection of latitude and longitude coordinates onto a plane is the realm of cartography.

    The formulae presented in this section provide varying degrees of accuracy.

    Spherical Earth projected to a plane

    This formula takes into account the variation in distance between meridians with latitude:

    D=Rsqrt{(Deltaphi)^2+(cos(phi_m)Deltalambda)^2},
    where:
    Deltaphi\,! and Deltalambda\,! are in radians;
    phi_m\,! must be in units compatible with the method used for determining cos(phi_m).\,!
    To convert latitude or longitude to radians use
    1^circ = (pi/180)\,mathrm{radians}.
    Note: This approximation is very fast and produces fairly accurate result for small distances[citation needed]. Also, when ordering locations by distance, such as in a database query, it is much faster to order by squared distance, eliminating the need for computing the square root.

    Ellipsoidal Earth projected to a plane

    The FCC prescribes essentially the following formulae in 47 CFR 73.208 for distances not exceeding 475 km /295 miles:[2]

    D=sqrt{(K_1Deltaphi)^2+(K_2Deltalambda)^2},
    where
    D\,! = Distance in kilometers;
    Deltaphi\,! and Deltalambda\,! are in degrees;
    phi_m\,! must be in units compatible with the method used for determining cos(phi_m);\,!
    egin{align} K_1&=111.13209-0.56605cos(2phi_m)+0.00120cos(4phi_m);\ K_2&=111.41513cos(phi_m)-0.09455cos(3phi_m)+0.00012cos(5phi_m).end{align}\,!
    It may be interesting to note that:
    K_1=Mfrac{pi}{180}\,! = kilometers per degree of latitude difference;
    K_2=cos(phi_m)Nfrac{pi}{180}\,! = kilometers per degree of longitude difference;
    where M\,! and N\,! are the meridional and its perpendicular, or "normal", radii of curvature (the expressions in the FCC formula are derived from the binomial series expansion form of M\,! and N\,!, set to the Clarke 1866 reference ellipsoid).

    Polar coordinate flat-Earth formula

    D=Rsqrt{ heta^2_1;oldsymbol{+}; heta^2_2;mathbf{-};2 heta_1 heta_2cos(Deltalambda)},
    where the colatitude values are in radians. For a latitude measured in degrees, the colatitude in radians may be calculated as follows: heta=frac{pi}{180}(90^circ-phi).\,!

    Spherical-surface formulae

    Main article: Great-circle distance

    If we are willing to accept a possible error of 0.5%, we can use formulas of spherical trigonometry on the sphere that best approximates the surface of the earth.

    The shortest distance along the surface of a sphere between two points on the surface is along the great-circle which contains the two points.

    The great-circle distance article gives the formula for calculating the distance along a great-circle on a sphere about the size of the Earth. That article includes an example of the calculation.

    Tunnel distance

    A tunnel between points on Earth is defined by a line through three-dimensional space between the points of interest. The great circle chord length may be calculated as follows for the corresponding unit sphere:

    egin{align} &Delta{X}=cos(phi_2)cos(lambda_2) - cos(phi_1)cos(lambda_1);\ &Delta{Y}=cos(phi_2)sin(lambda_2) - cos(phi_1)sin(lambda_1);\ &Delta{Z}=sin(phi_2) - sin(phi_1);\ &C_h=sqrt{(Delta{X})^2 + (Delta{Y})^2 + (Delta{Z})^2}.end{align}

    The tunnel distance between points on the surface of a spherical Earth is D = R C_h. For short distances (Dll R), this underestimates the great circle distance by D(D/R)^2/24.

    Ellipsoidal-surface formulae

    Main article: Geodesics on an ellipsoid
     
    Geodesic on an oblate ellipsoid

    An ellipsoid approximates the surface of the earth much better than a sphere or a flat surface does. The shortest distance along the surface of an ellipsoid between two points on the surface is along the geodesic. Geodesics follow more complicated paths than great circles and in particular, they usually don't return to their starting positions after one circuit of the earth. This is illustrated in the figure on the right where f is taken to be 1/50 to accentuate the effect. Finding the geodesic between two points on the earth, the so-called inverse geodetic problem, was the focus of many mathematicians and geodesists over the course of the 18th and 19th centuries with major contributions by Clairaut,[3] Legendre,[4] Bessel,[5] and Helmert.[6] Rapp[7] provides a good summary of this work.

    Methods for computing the geodesic distance are widely available in geographical information systems, software libraries, standalone utilities, and online tools. The most widely used algorithm is by Vincenty,[8] who uses a series which is accurate to third order in the flattening of the ellipsoid, i.e., about 0.5 mm; however, the algorithm fails to converge for points that are nearly anti-podal. (For details, see Vincenty's formulae.) This defect is cured in the algorithm given by Karney,[9] who employs series which are accurate to sixth order in the flattening. This results in an algorithm which is accurate to full double precision and which converges for arbitrary pairs of points on the earth. This algorithm is implemented in GeographicLib.[10]

    The exact methods above are feasible when carrying out calculations on a computer. They are intended to give millimeter accuracy on lines of any length; we can use simpler formulas if we don't need millimeter accuracy, or if we do need millimeter accuracy but the line is short. Rapp,[11] Chap. 6, describes the Puissant method, the Gauss mid-latitude method, and the Bowring method.[12]

    Lambert's formula for long lines

    Lambert's formulae[13] give accuracy on the order of 10 meters over thousands of kilometers. First convert the latitudes scriptstyle phi_1, scriptstyle phi_2 of the two points to reduced latitudes scriptstyle eta_1, scriptstyle eta_2

    an eta = (1 - f) an phi,

    where f is the flattening. Then calculate the central angle sigma in radians between two points (eta_1 , ; lambda_1) and (eta_2 , ; lambda_2) on a sphere in the usual way (law of cosines or haversine formula), with longitudes lambda_1 ; and lambda_2 ; being the same on the sphere as on the spheroid.

    P = frac { eta_1 + eta_2 }{2} qquad Q = frac {eta_2 - eta_1}{2}

    X = ( sigma - sin sigma) frac {sin^2 P cos^2 Q}{ cos^2 frac { sigma}{2}} qquad qquad Y = ( sigma + sin sigma) frac {cos^2 P sin^2 Q}{ sin^2 frac { sigma}{2}}

    mathrm{distance} = a igl( sigma - frac f2 (X + Y) igr)

    where a is the equatorial radius of the chosen spheroid.

    On the GRS 80 spheroid Lambert's formula is off by

    0 North 0 West to 40 North 120 West, 12.6 meters
    0N 0W to 40N 60W, 6.6 meters
    40N 0W to 40N 60W, 0.85 meter

    Bowring's method for short lines

    Bowring maps the points to a sphere of radius R′, with latitude and longitude represented as φ′ and λ′. Define

    A = sqrt{1 + e'^2cos^4 phi_1}, quad B = sqrt{1 + e'^2cos^2 phi_1},

    where the second eccentricity squared is

    e'^2 = frac{a^2 - b^2}{b^2} = frac{f(2-f)}{(1-f)^2}.

    The spherical radius is

    R' = frac{sqrt{1 + e'^2 }}{B^2} a.

    (The Gaussian curvature of the ellipsoid at φ1 is 1/R′2.) The spherical coordinates are given by

    egin{align} anphi_1' &= frac{ anphi}B,\ Deltaphi' &= frac{Delta phi}{B}iggl[1 + frac{3 e'^2 }{4 B^2}(Delta phi) sin (2 phi_1 + frac23 Delta phi )iggr],\ Deltalambda' &= ADeltalambda, end{align}

    where Deltaphi=phi_2-phi_1, Deltaphi'=phi_2'-phi_1', Deltalambda=lambda_2-lambda_1, Deltalambda'=lambda_2'-lambda_1'. The resulting problem on the sphere may be solved using the techniques for great-circle navigation to give approximations for the spheroidal distance and bearing. Detailed formulas are given by Rapp,[11] §6.5 and Bowring.[12]

    See also

    References

     
    1. Lambert, W. D (1942). "The distance between two widely separated points on the surface of the earth". J. Washington Academy of Sciences 32 (5): 125–130.

    External links
  • 相关阅读:
    aspscheduler+uwsgi定时任务执行多次
    django定时任务
    django记录用户操作模块
    python缩小图片
    pymysql同时执行多条语句时报错
    MySQL8.0 修改密码
    linux中安装python3
    Mysql高性能
    Mysql高性能
    Mysql高性能
  • 原文地址：https://www.cnblogs.com/zhoug2020/p/5414023.html
Copyright © 2011-2022 走看看