Samuel Learning

計算地球兩點間的距離

前提: 把地球當作標準球體計算
不妨令兩點的經緯度分別為A:(L1,B1),B:(L2,B2)，且B在A的東邊，它們和北極P構成一球面三角形PAB，設AB間的角距離為a，角PAB為A1，角PBA為A2，則有：

cosa=sin(B1)×sin(B2)＋cos(B1)×cos(B2)×cos(L2-L1);(邊的余弦定理)
sin(A1)=cosB2×sin(L2-L1)/sina;（正弦定理）
sin(A2)=cosB1×sin(L2-L1)/sina.（正弦定理）

public class TwoPoints{
  private static final double R_OF_EARTH=6.371E6;
  
  //public static void main(String[] args)
  //{
  //  System.out.println(TwoPoints.distance(45,30,45,-60));
  //  System.out.println(R_OF_EARTH*Math.PI/3);
  //}
  
  public static double distance(double lat1,double lon1,double lat2,double lon2)
  {
    return R_OF_EARTH*Math.acos(Math.sin(toRadians(lat1))*
      Math.sin(toRadians(lat2))+Math.cos(toRadians(lat1))*
      Math.cos(toRadians(lat2))*
      Math.cos(Math.abs(toRadians(lon1)-toRadians(lon2))));
  }
  
  private static double toRadians(double angle)
  {
    return angle*Math.PI/180;
  }
}

posted on 2005-08-30 09:47 MingIsMe 閱讀(408) 評論(0)  編輯  收藏 所屬分類: 06 J2EE