`
kalogen
  • 浏览: 878229 次
  • 性别: Icon_minigender_1
  • 来自: 杭州
社区版块
存档分类
最新评论

根据地球上任意两点的经纬度计算两点间的距离

阅读更多

地球是一个近乎标准的椭球体,它的赤道半径为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

在实际应用当中,一般是通过一个个体的邮政编码来查找该邮政编码对应的地区中心的经纬度,然后再根据这些经纬度来计算彼此的距离,从而估算出某些群体之间的大致距离范围(比如酒店旅客的分布范围-各个旅客的邮政编码对应的经纬度和酒店的经纬度所计算的距离范围-等等),所以,通过邮政编码查询经纬度这样一个数据库是一个很有用的资源。

 

根据经纬度计算该位置点与附近某一经纬度位置点的位置偏移角
2008-05-18 13:44

如何计算?

设某一城市有一个城市地标,可以视为城市的中心点 ,其经纬度已知。

分散在城市中心分别
   (1)由多层环形交通公路(环线),环城市地标为基准,在不同层的环线上分布着重要的设施。
    (2)分布着多个卫星城、主要居民区和学校、基础设施。

现控制无人机分别自城市地标出发,依次探访和航拍环城公路上的重要设施。

无人机配置有GPS导航和电子罗盘。

自城市中心地标处控制无人机起飞。

1 计算任意经纬度点间的距离

用GPS测出两个点的经纬度后,如何计算这两个点之间的距离呢。用的是便携式GPS,只有经纬度显示,无坐标显示和输出。


设两点A、B的经、纬度分别为(jA,wA)(jB,wB),则半径为R的球面上两点间的最短距离(大圆弧)为:
弧AB=R*arccos[sin(wA)sin(wB)+cos(wA)cos(wB)*cos(jA-jB)]


google maps的脚本里代码。

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;
}


;该函数用IDL语言编写,利用GIS中根据两点经纬度计算距离公式

;计算经纬度两点间的距离
Function Length2LatLon,pts0,pts1
    radius=6378.14d;地球平均半径,单位公里
;    print,pts0,pts1

;    length=radius * acos(sin(pts0[1]*!PI/180.0)*sin(pts1[1]*!PI/180.0)+cos(pts0[1]*!PI/180.0)*cos(pts1[1]*!PI/180.0)*cos(pts0[0]*!PI/180.0-pts1[0]*!PI/180.0))
    length=radius * acos(sin(pts0[1]*!DTOR)*sin(pts1[1]*!DTOR)+cos(pts0[1]*!DTOR)*cos(pts1[1]*!DTOR)*cos(pts0[0]*!DTOR-pts1[0]*!DTOR))
    return,length
end
pro test
pts0=[117.548802d,38.625238d]
pts1=[117.548625d,38.625139d]
print,length2latlon(pts0,pts1)


end

 

另外在IDL中提供了计算两点距离的方法map_2points可以方便的计算不同单位(Miles,Meters)的距离,还可以通过设置不同的参数来得到不同的距离、角度值,

Syntax:
Result = MAP_2POINTS( lon0, lat0, lon1, lat1 [, DPATH=value | , /METERS | , /MILES | , NPATH=integer{2 or greater} | , /PARAMETERS | , RADIANS=value] [, /RADIUS] [, /RHUMB] )

具体可以参考IDL 的Help里面讲解

如果要算的距离是椭球面的距离,就难些。

因为大地线的方程是一个微分方程,所以大地线长度是对一个微分方程的开方的积分,无法写成简单的公式.通常是直接用差分代替微分,用求和代替积分算出,这是很严格的.

下面是近似地由大地经纬度求椭球面距离的逼近方法:

如两点在同一纬线上,距离为N*cos(phi)*(lambda2-lambda1)。

如两点在同一经线上,距离为M*dphi从phi1到phi2的积分。

其中N=a/sqrt(d), d=1-e*e*sin(phi)*sin(phi),

          M=a*(1-e*e)/[d*sqrt(d)],

其中长短轴a,b与扁率偏心率的关系是:a/b=1-f=sqrt(1-e*e).

把积分区域分成等分n小段,积分可以写成求和,既可算出。这也是严格的。

如两点不在同一经纬线上,且两点距离很近,则小区域的椭球面可以视为平面,于是两点的距离近似地是上面算出的距离的平方和的开方。如两点距离较远,把两点间的经纬差分成等分n小段,每小段的距离的和就是所求的距离。

地图基础知识 http://www.51ditu.com

  • 地图应用接口中使用的经纬度单位是NTU,您需要把其它单位的经纬度值转化过来,以下为常用的转换关系:

基本转换:
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)
其格式的经纬值先转换度,再采用上面的运算。

与经纬度有关的基础知识:

NTU 相当于十万分之一度。

如何计算经度方向距离

比如: 点A的纬度为3995400,点B的纬度为3995300,则这两个点在沿着纬线的方向相差100米。
点A的经度为11695400,点B的经度为11695300,则这两个点在沿着经线的方向相差大约为77米,
该值的计算方法可以用近似公式:
经度方向距离 = 经度差 * cos(纬度值) = 100 * cos(39) = 77米。

如何计算两点间距离

比如:点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>

来自: http://hi.baidu.com/liuliangzhou/blog/item/dcad4913326f2c045baf53e8.html
分享到:
评论
1 楼 hzw2312 2016-10-21  
C = sin(MLatA)*sin(MLatB)*cos(MLonA-MLonB) + cos(MLatA)*cos(MLatB)

这里面的C的值是什么?

相关推荐

    经纬度计算距离,经纬度计算距离公式,matlab

    计算地球上两点间距离通常采用Haversine公式,这是一个基于球面三角学的公式,适用于近似地球的球体。Haversine公式如下: sin²(Δlat/2) + cos(lat1) * cos(lat2) * sin²(Δlon/2) = c 其中,lat1和lat2分别是...

    两点经纬度计算距离工具

    这个小工具可以帮助用户快速、准确地计算地球上任意两个地理位置之间的直线距离。下面将详细解释相关知识点。 1. **经纬度系统**:地球被划分为经度和纬度网格,用于定位地球表面的任何位置。经度(Longitude)是从...

    通过两点经纬度,计算他们之间的距离 - JAVA

    这个任务通常使用哈弗辛公式(Haversine formula)来实现,该公式基于球面三角学,可以计算地球上任意两点间的大圆距离。 哈弗辛公式的基本思路是将地球视为一个完美的球体,然后利用经纬度坐标来计算两点之间的...

    WGS84坐标下根据经纬度计算两点之间的距离

    通过这些经纬度值,可以计算出地球上两点间的距离。通常,这种距离计算基于球面三角学原理,利用地球的平均半径来近似计算两点间的大圆弧长度。 #### 三、计算公式详解 在给定的部分内容中,提供了一个计算两点间...

    经纬度计算距离、方位角等java源码

    使用Haversine公式可以计算地球上任意两点间的直线距离(也称为大圆距离),这是基于球面三角学的。Haversine公式考虑了地球半径,并转换了经纬度的角度值。在Java代码中,这通常会涉及到以下步骤: - 将经纬度从...

    根据定位经纬度计算距离.txt

    该函数使用球面三角公式来计算地球上两点间的最短路径距离,即所谓的“大圆距离”。 ##### 参数解释 - `lon`: 经度(Longitude) - `lat`: 纬度(Latitude) - `core`: 第一个参考点的经度 - `corn`: 第一个参考点...

    mysql函数-根据经纬度坐标计算距离

    在二维平面上,我们可以使用余弦定理来计算两点间的欧氏距离,但在地球上,需要使用哈弗辛公式(Haversine Formula): 哈弗辛公式如下: \[ a = \sin^2\left(\frac{\Delta \phi}{2}\right) + \cos(\phi_1) \cdot \...

    WGS84计算两点经纬度距离

    总结来说,WGS84坐标系统在地理计算中起着关键作用,通过Haversine公式,我们可以准确地计算出地球上任意两点间的距离。利用编程语言实现这个公式,可以方便地集成到各种应用中。而"jwd.exe"这样的工具则为用户提供...

    Android-依据两点间经纬度坐标,计算两点间距离

    根据提供的文件信息,本文将详细解释如何在Android应用开发中基于两点间的经纬度坐标来计算这两点之间的地理距离。此方法对于实现地图定位、路线规划等功能具有重要意义。 ### 核心概念 #### 地理坐标系统 地球...

    两经纬度之间的距离

    在IT行业中,尤其是在地理信息系统(GIS)或者基于位置服务(LBS)的应用中,计算两个经纬度之间的距离是一项常见的...通过Haversine公式,我们可以准确地计算出地球上任意两点间的距离,进而用于各种实际应用场景。

    易语言经纬度计算

    在地理坐标系统中,地球被抽象为一个球体,经纬度则是描述地球上任意位置的坐标。经度(longitude)是从本初子午线(0度经线)出发,向东西两侧各分180度,东经用"+"表示,西经用"-"表示;纬度(latitude)是从赤道...

    用经纬度计算两点距离软件

    总的来说,这个软件通过处理经纬度数据,利用数学公式如Haversine来计算地球上任意两点之间的距离。这在物流、导航、气象学等领域都有广泛应用。用户只需提供起点和终点的经纬度,软件就能快速给出准确的距离。

    Excel计算两经纬度之间的距离

    首先,我们要知道计算两点间距离的方法通常是基于球面三角学,这里使用的是哈弗辛公式(Haversine formula),它能够计算地球上任意两点之间的大圆距离。公式如下: \[ a = \sin^2\left(\frac{\Delta \phi}{2}\...

    根据经纬度计算距离_计算_根据计算两点的长度_

    Vincenty's Formula是一种迭代方法,适用于任意两点间的最短距离。虽然比哈弗辛公式复杂,但能提供更高的精度。其公式涉及了更多的三角函数和坐标转换,需要多次迭代直到达到预定的精度阈值。 在编程实现中,你可以...

    经纬度距离计算+坐标转换

    经纬度是用来确定地球上任意位置的坐标系统,而距离计算则涉及到地理空间数据的处理,这对于地图应用、定位服务以及物流规划等领域至关重要。以下是关于这两个主题的详细解释: 一、经纬度系统 经纬度是一种全球性...

    Java地图工具类计算地球上两点之间的距离

    在Java编程中,地图工具类通常用于处理与地理坐标系统相关的...通过这样的工具类,开发者可以方便地在Java应用程序中获取地球上任意两点的距离信息,从而应用于各种场景,如导航、地理信息系统(GIS)或物流配送等。

    经纬度距离计算工具绿色版.rar

    经纬度距离计算工具是一款专为地理定位和导航需求设计的小型绿色软件,它基于地球的椭球体模型,准确地计算地球上任意两点之间的直线距离。这个工具尤其适用于那些需要进行地图分析、户外活动规划或者物流路线优化的...

    经纬度计算问题c++写的

    这些坐标定义了地球上任意一点的位置。 在C++中,进行经纬度计算通常涉及以下方面: 1. 数据类型选择:由于经纬度可能涉及到小数部分,所以一般使用浮点数类型如`float`或`double`来存储。考虑到精度问题,推荐...

    经纬度计算距离 源码

    这个源码可能是一个程序或者库,用于在地球上任意两点的经纬度坐标之间计算出它们的球面距离。接下来,我们将深入探讨这个主题。 首先,我们来理解经纬度坐标系统。经纬度是地球表面位置的一种表达方式,经度表示...

    经纬度测距函数,传入参数为地球上两点的经纬度,输出为两点的距离。

    在IT领域,尤其是在地理信息系统(GIS)和导航应用中,计算地球上两点之间的距离是一个常见的需求。这个需求可以通过使用经纬度测距函数来实现。这里,我们深入探讨一下这个主题。 首先,地球并不是一个完美的球体...

Global site tag (gtag.js) - Google Analytics