谈谈从经纬度坐标到欧式距离的转化【matlab实现】

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

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

基于此，我用matlab实现编程：

最后我们不妨再编个程试试，计算全国各省会（含港澳台）之间的直线距离，

原始数据是这样：

编程如下：

输出我们的结果文件，省会城市欧式距离表，打开后如下，单位：/km

这样，任务就完成啦！

谈谈从经纬度坐标到欧式距离的转化【matlab实现】相关推荐

根据经纬度坐标计算实际距离
double hypot( double x, double y) ... { return sqrt(x * x + y * y); } double distance( d ...
经纬度计算距离html,js版谷歌地图计算两经纬度坐标点的距离
球面上两点间大圆劣弧的距离html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...
java版谷歌地图计算两经纬度坐标点的距离
球面上两点间大圆劣弧的距离 public class GoogleMapHelper {private static double EARTH_RADIUS = 6378.137; //地球半径//将 ...
经纬度坐标与实际距离的计算转换
在做数学建模题目的时候,经常遇到题目所给的数据是经纬度坐标,如城市的经纬度坐标等,在求解此类问题时,第一步往往是求解各点之间的距离矩阵. 求解距离矩阵的时候,我们通常会想到用欧氏距离公式直接求解两点之 ...
python计算坐标点欧式距离_Python计算一个点到所有点的欧式距离实现方法
Python计算一个点到所有点的欧式距离实现方法如下所示: distances = np.sqrt(np.sum(np.asarray(airportPosition - x_vals)**2, a ...
Hive计算两个经纬度坐标的直线距离
-- Hive计算两个经纬度坐标直线距离(单位:公里) -- 公式 -- ASIN( -- SQRT( -- POWER(SIN((lat1-lat2)*ACOS(-1)/360),2) + -- C ...
python计算坐标点欧式距离_计算机视觉课堂笔记-4
三维计算机视觉相机模型与多视几何相机标定与稀疏重建立体视觉三维重建 (1)人脑中的腹部通道进行物体视觉,目前主要采用基于学习的方法 (2)背部通道进行空间的视觉,目前主要采用基于几何的方法三维 ...
经纬度坐标点和距离之间的转换
1.纬度相同,经度不同在纬度相同的情况下: 经度每隔0.00001度,距离相差约1米: 每隔0.0001度,距离相差约10米: 每隔0.001度,距离相差约100米: 每隔0.01度,距离相差约10 ...
Java：计算地球上两个经纬度坐标之间的距离-geodesy和geotools实现
目录方式一:自定义公式计算方式二:geodesy计算距离方式三:geotools计算距离两个点的经纬度 latitude纬度 longitude经度地点 22.678611 113.8056 ...

谈谈从经纬度坐标到欧式距离的转化【matlab实现】

谈谈从经纬度坐标到欧式距离的转化【matlab实现】相关推荐

最新文章

热门文章