百度地圖開發：如何計算兩點之間距離  

2012-11-26 22:04:27|  分類： 百度地圖開發 |  标簽：android開發  百度地圖開發   |舉報 |字号大中小 訂閱

路線規劃提供了擷取路線距離的方法，見MKRoutePlan 類的 getDistance 方法。如果是計算任意兩點的距離，有兩種方法：一種利用勾股定理計算，适用于兩點距離很近的情況；一種按标準的球面大圓劣弧長度計算，适用于距離較遠的情況。

 
 1.           static           double DEF_PI           =           3.14159265359          ;           // PI       
      
 2.           static           double DEF_2PI          =           6.28318530712          ;           // 2*PI       
      
 3.           static           double DEF_PI180          =           0.01745329252          ;           // PI/180.0       
      
 4.           static           double DEF_R           =          6370693.5          ;           // radius of earth       
      
 5.           public           double GetShortDistance          (          double lon1,           double lat1,           double lon2,           double lat2          )       
      
 6.           {       
      
 7.            double ew1, ns1, ew2, ns2          ;       
      
 8.            double dx, dy, dew          ;       
      
 9.            double distance          ;       
      
 10.            // 角度轉換為弧度       
      
 11.            ew1           = lon1           * DEF_PI180          ;       
      
 12.            ns1           = lat1           * DEF_PI180          ;       
      
 13.            ew2           = lon2           * DEF_PI180          ;       
      
 14.            ns2           = lat2           * DEF_PI180          ;       
      
 15.            // 經度差       
      
 16.            dew           = ew1           - ew2          ;       
      
 17.            // 若跨東經和西經180 度，進行調整       
      
 18.            if           (dew           > DEF_PI          )       
      
 19.            dew           = DEF_2PI           - dew          ;       
      
 20.            else           if           (dew           <           -DEF_PI          )       
      
 21.            dew           = DEF_2PI           + dew          ;       
      
 22.            dx           = DEF_R           *           Math.          cos          (ns1          )           * dew          ;           // 東西方向長度(在緯度圈上的投影長度)       
      
 23.            dy           = DEF_R           *           (ns1           - ns2          )          ;           // 南北方向長度(在經度圈上的投影長度)       
      
 24.            // 勾股定理求斜邊長       
      
 25.            distance           =           Math.          sqrt          (dx           * dx           + dy           * dy          )          ;       
      
 26.            return distance          ;       
      
 27.           }       
      
 28.           public           double GetLongDistance          (          double lon1,           double lat1,           double lon2,           double lat2          )       
      
 29.           {       
      
 30.            double ew1, ns1, ew2, ns2          ;       
      
 31.            double distance          ;       
      
 32.            // 角度轉換為弧度       
      
 33.            ew1           = lon1           * DEF_PI180          ;       
      
 34.            ns1           = lat1           * DEF_PI180          ;       
      
 35.            ew2           = lon2           * DEF_PI180          ;       
      
 36.            ns2           = lat2           * DEF_PI180          ;       
      
 37.            // 求大圓劣弧與球心所夾的角(弧度)       
      
 38.            distance           =           Math.          sin          (ns1          )           *           Math.          sin          (ns2          )           +           Math.          cos          (ns1          )           *           Math.          cos          (ns2          )           *           Math.          cos          (ew1           - ew2          )          ;       
      
 39.            // 調整到[-1..1]範圍内，避免溢出       
      
 40.            if           (distance           >           1.0          )       
      
 41.                 distance           =           1.0          ;       
      
 42.            else           if           (distance           <           -          1.0          )       
      
 43.                  distance           =           -          1.0          ;       
      
 44.            // 求大圓劣弧長度       
      
 45.            distance           = DEF_R           *           Math.          acos          (distance          )          ;       
      
 46.            return distance          ;       
      
 47.           }       
      
 48.           double mLat1           =           39.90923          ;           // point1緯度       
      
 49.           double mLon1           =           116.357428          ;           // point1經度       
      
 50.           double mLat2           =           39.90923          ;          // point2緯度       
      
 51.           double mLon2           =           116.397428          ;          // point2經度       
      
 52.           double distance           = GetShortDistance          (mLon1, mLat1, mLon2, mLat2          )          ;