经度/纬度计算两点间的距离

环境：struts+hibernate，有一个article(文章类)类，文章中有三个字段、进度、纬度、距离(其中距离在是求出来的).
目标是:1.查询出最新发布的1000个文章。
2.根据经度纬度求出article中的距离。
3.在根据距离排序.
问题是：1000条记录存在临时表中,还是存在像Hashtable的列表中，之后再对1000条记录当中的每个article修改，求出距离根据经度和纬度。
考虑到效率问题，或者有其他什么好方法没有？可以存储1000条记录，进行我想的操作。
或者有其他好方法，希望顶一下。

 

解决方法:

1.查询出1000条记录，新的记录。

2.之后计算出经度/纬度之间的距离。存放在临时表中。

3.排序，加分页

计算两点间经度/纬度的距离公式如

java 代码

package com.hipengyou.apply.util;
/**
 * 关于地图类
 * @author
 *
 */
public class GoogleMap {
 private static final double EARTH_RADIUS = 6378.137;
 private static double rad(double d)
 {
    return d * Math.PI / 180.0;
 }
/**
 * 计算两点间距离（关于经度/纬度）
 * @param lat1
 * @param lng1
 * @param lat2
 * @param lng2
 * @return
 */
 public static double GetDistance(double lat1, double lng1, double lat2, double lng2)
 {
    double radLat1 = rad(lat1);
    double radLat2 = rad(lat2);
    double a = radLat1 - radLat2;
    double b = rad(lng1) - rad(lng2);
    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)));
    s = s * EARTH_RADIUS;
    s = Math.round(s * 10000) / 10000;
    return s;
 }
 public static double D_jw(double wd1,double jd1,double wd2,double jd2)
 {
     double x,y,out;
     double PI=3.14159265;
     double R=6.371229*1e6;

     x=(jd2-jd1)*PI*R*Math.cos( ((wd1+wd2)/2) *PI/180)/180;
     y=(wd2-wd1)*PI*R/180;
     out=Math.hypot(x,y);
     return out/1000;
 }

public static void main(String args[]){
  System.out.println(GoogleMap.D_jw(39.89945,116.40969,39.14307,117.20010));
 }
}

 

js代码

 

private const double EARTH_RADIUS = 6378.137;
private static double rad(double d)
{
   return d * Math.PI / 180.0;
}

public static double GetDistance(double lat1, double lng1, double lat2, double lng2)
{
   double radLat1 = rad(lat1);
   double radLat2 = rad(lat2);
   double a = radLat1 - radLat2;
   double b = rad(lng1) - rad(lng2);
   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)));
   s = s * EARTH_RADIUS;
   s = Math.Round(s * 10000) / 10000;
   return s;
}

非常感谢，帮了我大忙了：）
虽然我也没看明白到底原理是什么，但验算了A（60，30），B（60，90）两点之间，此段代码和我用余弦定理算出来的结果很一致。
余弦定理的步骤是：1、算A、B弦长：地球半径R*cos（经度差60）=R/2；
2、算角AOB,O为地球圆心，利用余弦定理，
cosAOB=(2R*R-(R/2)^2) /2*R*R=7/8;
3、弧AB的长为：R*arc cos(7/8);求毕 

 

以下地址是www.51ditu.com上计算l两个间经度纬度的距离地址

http://api.51ditu.com/docs/mapsapi/examples/LTFunction/getPointsDistance.htm