根据地球上任意两点的经纬度计算两点间的距离
地球是一个近乎标准的椭球体，它的赤道半径为6378.140千米，极半径为6356.755千米，平均半径6371.004千米。如果我们假设地球是一个完美的球体，那么它的半径就是地球的平均半径，记为R。如果以0度经线为基准，那么根据地球表面任意两点的经纬度就可以计算出这两点间的地表距离（这里忽略地球表面地形对计算带来的误差，仅仅是理论上的估算值）。设第一点A的经纬度为(LonA, LatA)，第二点B的经纬度为(LonB, LatB)，按照0度经线的基准，东经取经度的正值(Longitude)，西经取经度负值(-Longitude)，北纬取90-纬度值(90-Latitude)，南纬取90+纬度值(90+Latitude)，则经过上述处理过后的两点被计为(MLonA, MLatA)和(MLonB, MLatB)。那么根据三角推导，可以得到计算两点距离的如下公式：
C = sin(MLatA)*sin(MLatB)*cos(MLonA-MLonB) + cos(MLatA)*cos(MLatB)
Distance = R*Arccos(C)*Pi/180
这里，R和Distance单位是相同，如果是采用6371.004千米作为半径，那么Distance就是千米为单位，如果要使用其他单位，比如mile，还需要做单位换算，1千米=0.621371192mile
如果仅对经度作正负的处理，而不对纬度作90-Latitude(假设都是北半球，南半球只有澳洲具有应用意义)的处理，那么公式将是：
C = sin(LatA)*sin(LatB) + cos(LatA)*cos(LatB)*cos(MLonA-MLonB)
Distance = R*Arccos(C)*Pi/180
以上通过简单的三角变换就可以推出。
如果三角函数的输入和输出都采用弧度值，那么公式还可以写作：
C = sin(LatA*Pi/180)*sin(LatB*Pi/180) + cos(LatA*Pi/180)*cos(LatB*Pi/180)*cos((MLonA-MLonB)*Pi/180)
Distance = R*Arccos(C)*Pi/180
也就是：
C = sin(LatA/57.2958)*sin(LatB/57.2958) + cos(LatA/57.2958)*cos(LatB/57.2958)*cos((MLonA-MLonB)/57.2958)
Distance = R*Arccos(C) = 6371.004*Arccos(C) kilometer = 0.621371192*6371.004*Arccos(C) mile = 3958.758349716768*Arccos(C) mile
在实际应用当中，一般是通过一个个体的邮政编码来查找该邮政编码对应的地区中心的经纬度，然后再根据这些经纬度来计算彼此的距离，从而估算出某些群体之间的大致距离范围(比如酒店旅客的分布范围-各个旅客的邮政编码对应的经纬度和酒店的经纬度所计算的距离范围-等等)，所以，通过邮政编码查询经纬度这样一个数据库是一个很有用的资源。
?
基本转换：
NTU = 度*100000
度 = NTU/100000
例如：
经度 = 116.21345°
纬度 = 39.445875°
NTU经度：116.21345 * 100000 = 11621345（NTU）
NTU纬度：39.445874*100000=3944587(NTU)
度分转换：
将度分单位数据转换为度单位数据
度=度+分/60
例如：
经度 = 116°20.12’
纬度 = 39°12.34’
经度 = 116 + 20.12 / 60 = 116.33533°
纬度 = 39 + 12.34 / 60 = 39.20567°
NTU经度 = 116.33533 * 100000 = 11633533（NTU）
NTU纬度 = 39.20567 * 100000 = 3920567（NTU）
度分秒转换：
将度分秒单位数据转换为度单位数据
度 = 度 + 分 / 60 + 秒 / 60 / 60
例如：
经度 = 116°20’43”
纬度 = 39°12’37”
经度 = 116 + 20 / 60 + 43 / 60 / 60 = 116.34528°
纬度 = 39 + 12 / 60 + 37 / 60 / 60 = 39.21028°
NTU经度 = 116.34528 * 100000 = 11634528（NTU）
NTU纬度 = 39.21028 * 100000 = 3921028（NTU）
其格式的经纬值先转换度，再采用上面的运算。
如何计算两点间距离
比如：点A的经度为11695400，纬度为3995400。点B的经度为11695300，纬度为3995300。
公式:两点间距离 = [ (A点经度 - B点经度)^2 + (A点纬度 - B点纬度)^2 ] ^ (1/2) = [ (11695400 - 11695300)^2
+ (3995400 - 3995300)^2 ] ^(1/2) =(10000+10000) ^ (1/2) =141米
公式说明:加法 + ，减法 - ，乘法 *，除法 /，幂运算 ^(1/2)表示平方根 ^2表示平方。
?
返回顶部
?
使用的比例级别总共有0-12级共13个级别，他们和实际的比例尺的对应关系分别是：
0级 -------------------- 1:5 000
1级 -------------------- 1:10 000
2级 -------------------- 1:20 000
3级 -------------------- 1:40 000
4级 -------------------- 1:80 000
5级 -------------------- 1:160 000
6级 -------------------- 1:320 000
7级 -------------------- 1:640 000
8级 -------------------- 1:1 280 000
9级 -------------------- 1:2 560 000
10级 ------------------- 1:5 120 000
11级 ------------------- 1:10 240 000
12级 ------------------- 1:20 480 000
13级 ------------------- 1:40 960 000
14级 ------------------- 1:81 920 000
简单的说，0级地图最详细，12级地图最简单。
返回顶部
?
比例级别定义表?
城市代码是某个城市的英文简称，标注在城市名称后面，便于用户记忆，点击查看城市代码表：
?
城市代码表
2 计算偏移角。
?
?
?
?
?
完整的列子如下：
?
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" >
<head runat="server">
??? <title>无标题页</title>
???
??? <script type="text/javascript">
???
??? function distanceByLnglat(lng1,lat1,lng2,lat2)
{
??? var radLat1 = Rad(lat1);
??? var radLat2 = Rad(lat2);
??? var a = radLat1 - radLat2;
??? var b = Rad(lng1) - Rad(lng2);
??? var 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 * 6378137.0;// 取WGS84标准参考椭球中的地球长半径(单位:m)
??? s = Math.round(s * 10000) / 10000;
alert(s);
???
// //下面为两点间空间距离（非球面体）
// var value= Math.pow(Math.pow(lng1-lng2,2)+Math.pow(lat1-lat2,2),1/2);
// alert(value);
}
function Rad(d)
{
??? return d * Math.PI / 180.0;
}
???
??? </script>
</head>
<body>
??? <form id="form1" runat="server">
??? <div onclick ="distanceByLnglat(116.95400,39.95400,116.95300,39.95300);">
??? test me
??? </div>
??? </form>
</body>
</html>
