根据经纬度计算该位置点与附近某一经纬度位置点的位置偏移角
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小段,每小段的距离的和就是所求的距离。
地图基础知识
* 地图应用接口中使用的经纬度单位是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表示平方。
适用范围:未加密经纬度/加密经纬度
- 大小: 118.5 KB
分享到:
相关推荐
1、已知两点经纬度求距离方位角 2、已知一点经纬度、距离、方位角求另一点经纬度。 java源码。本人手敲,亲测简单实用。可用于百度地图、Google地图开发,误差极小基本可忽略。
经纬度转换软件,驴友必备。能转换十进制经纬度,偏移经纬度转真实经纬度。
经纬度偏差计算表格-含公式经纬度偏差计算表格-含公式
GPS到地图会有偏差,应要加上偏移量才会正确
在百度、谷歌、soso、等地图上,对经纬度偏移进行正确的纠偏的最佳技术解决方案,内容包含技术解决与相关获取资料方式!
最新全国百度、高德、谷歌等各类地图经纬度偏移,经纬度纠偏,gps偏移,gps纠偏数据库。
所谓的地图数据加密偏移,其实就是用一个偏移算法对地图的经纬度做一个加减偏移量,从而达到与实际地图不一致。这个偏移算法本身是没有什么规律可言的,每家地图服务商都有一套自己的加密偏移算法,既然算法上没有...
GPS的经纬度在google地图上有一定的偏移,该算法,实现了这个纠正
Google baidu 地图 gps偏移(经纬度偏移)的技术纠偏解决方案,在gps实际开发中的应用,最权威的解决方案
gps开发和互联网项目开发中,gps偏移的权威技术解决(gps偏移,经纬度偏移,gps纠偏,经纬度纠偏),中国唯一的官方权威解决方案
GPS返回的经纬度是以火星坐标为准,以此经纬度在Google地球上和Google地图上查找,两者显示存在偏差,本示例解决了此问题
baidumap 经纬度位移偏差 非常精确,测试过的,给大家分享分享
输入AB两个地点的经纬度,结果为B相对A正北顺时针偏移角度
用于GPS经纬度,转换为百度经纬度的偏移算法,采用简单的JS
google地图经纬度偏移修正算法完美解决方案.doc
Google baidu 地图 gps偏移(经纬度偏移)的技术纠偏解决方案 讲解如何通过官方API实现0偏移 不要相信网上的数据库或者算法之类的,感觉都有些扯淡,并且超级影响系统性能,通过官方的解决方案更完美,更稳定
国内地图有所谓的偏移问题,这是一个可以自动修正所谓偏移的经纬度小工具,不用安装的。
对GSP经纬度数据进行纠偏计算,过滤掉偏移无效的数据。画出规整的轨迹。
纠正高德地图定位位置偏差案例,只需要自己在高德开发中心申请的key替换到我的key,一运行就会弹出alert(经纬度)
借鉴网络开源的经纬度转化代码,由VB6重写编写的百度、高德地图经纬度转化工具,可进行批量经纬度转化,有部分的偏移情况,但是基本能控制在500m以内,高纬度地区偏移在米级