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