经纬度转公里 matlab,两地经纬度如何转换成两地距离的程序

function r = geodistance( ci , cf , m )

% Calculates the distance in meters between two points on earth surface.

%

% SYNTAX: r = geodistance( coordinates1 , coordinates2 , method ) ;

%

%          Where coordinates1 = [longitude1,latitude1] defines the

%          initial position and coordinates2 = [longitude2,latitude2]

%          defines the final position.

%          Coordinates values should be specified in decimal degrees.

%          Method can be an integer between 1 and 23, default is m = 6.

%         Methods 1 and 2 are based on spherical trigonometry and a

%         spheroidal model for the earth, respectively.

%          Methods 3 to 23 use Vincenty's formulae, based on ellipsoid

%         parameters.

%         Here it follows the correspondence between m and the type of

%         ellipsoid:

%

%         m =  3 -> ANS ,        m =  4 -> GRS80,    m = 5 -> WGS72,

%         m =  6 -> WGS84,       m =  7 -> NSWC-9Z2,

%         m =  8 -> Clarke 1866, m =  9 -> Clarke 1880,

%         m = 10 -> Airy 1830,

%         m = 11 -> Bessel 1841 (Ethiopia,Indonesia,Japan,Korea),

%         m = 12 -> Bessel 1841 (Namibia),

%         m = 13 -> Sabah and Sarawak (Everest,Brunei,E.Malaysia),

%         m = 14 -> India 1830, m = 15 -> India 1956,

%         m = 16 -> W. Malaysia and Singapore 1948,

%         m = 17 -> W. Malaysia 1969,

%         m = 18 -> Helmert 1906,m = 19 -> Helmert 1960,

%         m = 20 ->