地图两个坐标点之间的距离

一次程序测试过程中发现之前的Oracle存储过程计算地图两个坐标点之间的距离有误差，之后在百度中完善更新了一下，具体过程如下：

这是之前的函数，有误差

1. CREATE OR REPLACE FUNCTION GetDistance(lat1 number,
2.                                        lng1 number,
3.                                        lat2 number,
4.                                        lng2 number) RETURN NUMBER is
5. earth_padius number := 6378.137;
6. radLat1      number := Radian(lat1);
7. radLat2      number := Radian(lat2); 
8. a            number := radLat1 - radLat2;
9. b            number := Radian(lng1) - Radian(lng2);
10. s            number := 0;
11. begin
12.   s := 2 *Asin(Sqrt(power(sin(a / 2), 2) + cos(radLat1) * cos(radLat2) * power(sin(b / 2), 2)));
13.   s := s * earth_padius;
    
14.   s := Round(s * 10000) / 10;
15. return s;
16. end;  

后来在网上找的一个存储过程：

    oracle存储过程：

1. create or replace procedure SP_GET_DISTANCE  
2. (cx in number,cy in number,sx in number, sy in number,distance out varchar2)   
3. is  
4. d number;  
5. x number;  
6. y number;  
7. r number;  
8. pi number;  
9. begin  
10. --开始计算  
11.   r:=6371229;--地球半径  
12.   pi:=3.14159265358979323;--圆周率  
13.   x:=(sx-cx)*pi*r*cos((sy+cy)/2*pi/180)/180;  
14.   y:=(sy-cy)*pi*r/180;  
15.   d:=SQRT(power(x,2)+power(y,2));  
16.   distance:=to_char(d,9999999999999.99);  
17. end SP_GET_DISTANCE;  

    java代码：

1. package com.wpn.web.util;  
2.   
3. public class Distance {  
4.     private final static double PI = 3.14159265358979323;// 圆周率  
5.     private final static double R = 6371229;  // 地球的半径  
6.   
7.     private Distance() {  
8.     }  
9.   
10.     /** 
11.      * 纬度lat 经度lon 
12.      * @param longt1 
13.      * @param lat1 
14.      * @param longt2 
15.      * @param lat2 
16.      * @return 
17.      */  
18.     public static double getDistance(double longt1, double lat1, double longt2, double lat2) {  
19.         double x, y, distance;  
20.         x = (longt2 - longt1) * PI * R * Math.cos(((lat1 + lat2) / 2) * PI / 180) / 180;  
21.         y = (lat2 - lat1) * PI * R / 180;  
22.         distance = Math.hypot(x, y);  
23.         return distance;  
24.     }  
25.   
26.     /*public enum GaussSphere { 
27.         Beijing54, Xian80, WGS84, 
28.     } 
29.  
30.     private static double Rad(double d) { 
31.         return d * Math.PI / 180.0; 
32.     } 
33.     public static double DistanceOfTwoPoints(double lng1, double lat1, double lng2, double lat2, GaussSphere gs) { 
34.         double radLat1 = Rad(lat1); 
35.         double radLat2 = Rad(lat2); 
36.         double a = radLat1 - radLat2; 
37.         double b = Rad(lng1) - Rad(lng2); 
38.         double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a / 2), 2) + Math.cos(radLat1) * Math.cos(radLat2) * Math.pow(Math.sin(b / 2), 2))); 
39.         s = s * (gs == GaussSphere.WGS84 ? 6378137.0 : (gs == GaussSphere.Xian80 ? 6378140.0 : 6378245.0)); 
40.         s = Math.round(s * 10000) / 10000; 
41.         return s; 
42.     }*/  
43.       
44.     public static void main(String[] arg){  
45.         double longt1 = 116.515502;  
46.         double lat1 = 39.863898;  
47.         double longt2 = 116.304187;  
48.         double lat2 = 40.052584;  
49.         System.out.println(getDistance(longt1,lat1,longt2,lat2));  
50.     }  
51. }  

摘自https://www.cnblogs.com/feizhuliu/archive/2012/12/07/6312112.html

把Oracle存储过程改为我所需要的函数：

  CREATE OR REPLACE FUNCTION GetDistance_new(cx number,
                                           cy number,
                                           sx number,
                                           sy number) RETURN NUMBER is
  r  number;
  pi number;
  x  number;
  y  number;
  d  number;
  s  number;
begin
  r  := 6371229; --地球半径  
  pi := 3.14159265358979323; --圆周率  
  x  := (sx - cx) * pi * r * cos((sy + cy) / 2 * pi / 180) / 180;
  y  := (sy - cy) * pi * r / 180;
  d  := SQRT(power(x, 2) + power(y, 2));
  s  := Round(d, 2);
  return s;

end;

完成。

分享一个测试两个地理坐标直线距离的百度API：http://developer.baidu.com/map/jsdemo.htm#a6_1

高德地图测试两个地理坐标直线距离：

http://lbs.amap.com/api/javascript-api/example/calcutation/calculate-distance-between-two-markers