经纬度算距离

下面以54坐标系为例，介绍具体的计算方法。

完成经纬度BL 到平面直角坐标XY 的换算。在EXCEL 中，选择输入公式的起始单元格，例如：第2行第1列（A2格）为起始单元格，各单元格的格式如下：单元格；单元格内容；说明

A2；输入中央子午线，以度. 分秒形式输入，如115度30分则输入115.30；起算数据L0

B2；=INT(A2)+(INT(A2*100)-INT(A2)*100)/60+(A2*10000-INT(A2*100)*100)/3600；把L0化成度

C2；以度小数形式输入纬度值，如38°14′20″则输入38.1420；起算数据B D2；以度小数形式输入经度值；起算数据L

E2；=INT(C2)+(INT(C2*100)-INT(C2)*100)/60+(C2*10000-INT(C2*100)*100)/3600；把B 化成度

F2；=INT(D2)+(INT(D2*100)-INT(D2)*100)/60+(D2*10000-INT(D2*100)*100)/3600；把L 化成度

G2；=F2-B2；L-L0

H2；=G2/57.[1**********]23；化作弧度

I2；=TAN(RADIANS(E2))；Tan(B)

J2；=COS(RADIANS(E2))；COS(B)

K2；=0.[1**********]5*J2*J2

L2；=I2*I2

M2；=1+K2

N2；=6399698.9018/SQRT(M2)

O2；=H2*H2*J2*J2

P2；=I2*J2

Q2；=P2*P2

R2；=(32005.78006+Q2*(133.92133+Q2*0.7031))

S2；=6367558.49686*E2/57.[1**********]-P2*J2*R2+((((L2-58)*L2+61)*

O2/30+(4*K2+5)*M2-L2)*O2/12+1)*N2*I2*O2/2

计算结果X

T2；=((((L2-18)*L2-(58*L2-14)*K2+5)*O2/20+M2-L2)*O2/6+1)*N2*(H2*J2) 计算结果Y

按上面表格中的公式输入到相应单元格后，就可方便地由经纬度求得平面直角坐标。当输入完所有的经纬度后，用鼠标下拉即可得到所有的计算结果。经纬度算距离的算法两个：

设地球上某点的经度为lambda, 纬度为phi ，

则这点的空间坐标是

x=cos(phi)*cos(lambda)

y=cos(phi)*sin(lambda)

z=sin(phi)

设地球上两点的空间坐标分别为(x1,y1,z1),(x2,y2,z2)

则它们的夹角为

A=acos(x1*x2+y1*y2+z1*z2)，A 是角度

则两地距离为

A/180*pi*R，其中R 为地球平均半径6371

误差不超过1%

CREATE OR REPLACE FUNCTION Get_Distance(

bw1 IN NUMBER, dj1 IN NUMBER,

bw2 IN NUMBER, dj2 IN NUMBER) RETURN NUMBER

n_size NUMBER;

BEGIN

SELECT 3.1415926 * (12756/2) / 180 * ACOS( SIN(bw1)*SIN(bw2) +

COS(bw1)*COS(bw2)*COS(dj1-dj2)) INTO n_size

FROM dual;

RETURN n_size;

EXCEPTION

WHEN NO_DATA_FOUND THEN

RETURN 0;

END;

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;

}据说是google maps的代码