matlab利用经纬度计算距离_【MATLAB】经纬度换算距离

文章来源于我的B站专栏:用经纬度算距离​www.bilibili.com

突然对于经纬度与距离感兴趣了(公选课讲到东风系列弹道导弹射程),就想了解一下如何通过经纬度来计算距离。百度了一下,觉得不满意就自己尝试做一下,都是些基本的数学知识,感兴趣玩玩。(如果有错漏之处,欢迎指正!)

首先说思路:经纬度转换为空间直角坐标系,由此得到两个向量,求出向量夹角,由向量夹角和地球半径求出弧长,即距离。

首先说一下经纬度,经度,是地球上一个地点离一根被称为本初子午线的南北方向走线以东或以西的度数。本初子午线的经度是0°,地球上其它地点的经度是向东到180°或向西到180°。

纬度,是指某点与地球球心的连线和地球赤道面所成的线面角,其数值在0至90度之间。位于赤道以北的点的纬度叫北纬,记为N,位于赤道以南的点的纬度称南纬,记为S。

这里定义:X轴为地心到(经度:0°、纬度:0°)的向量。

Y轴为地心到(经度:90°、纬度:0°)的向量。

Z轴为地心到(纬度:90°)的向量。经纬度与空间直角坐标系的对应关系(C为球面上一点)

由此可以得到由经纬度到空间直角坐标系的对应关系:(注意正弦函数的输入为弧度)单位向量

由角度到弧度对应关系:度-->弧度

  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
### 回答1: 在MATLAB中,可以使用haversine公式来计算两个经纬度坐标之间的球面距离。具体步骤如下: 1. 首先,将经纬度坐标转换为弧度制。这可以通过将经纬度值除以180再乘以π来实现。 2. 然后,使用haversine公式来计算球面距离。haversine公式如下: d = 2 * R * asin(sqrt(sin((lat2-lat1)/2)^2 + cos(lat1) * cos(lat2) * sin((lon2-lon1)/2)^2)) 其中,d是距离,R是地球的半径(一般情况下为6371千米),lat1和lon1是第一个点的纬度和经度,lat2和lon2是第二个点的纬度和经度。 3. 最后,根据需要可以将距离换算为其他单位。 下面是一个MATLAB代码示例,可以根据上述步骤计算两个经纬度坐标之间的球面距离: ```matlab function distance = calculateDistance(lat1, lon1, lat2, lon2) R = 6371; % 地球半径,单位:千米 % 将经纬度转换为弧度制 lat1 = deg2rad(lat1); lon1 = deg2rad(lon1); lat2 = deg2rad(lat2); lon2 = deg2rad(lon2); % 使用haversine公式计算球面距离 dlat = lat2 - lat1; dlon = lon2 - lon1; a = sin(dlat/2)^2 + cos(lat1) * cos(lat2) * sin(dlon/2)^2; c = 2 * atan2(sqrt(a), sqrt(1-a)); distance = R * c; end ``` 使用这个函数,可以通过传入经纬度坐标调用该函数来计算两点之间的球面距离。例如: ```matlab lat1 = 31.21563; lon1 = 121.50891; lat2 = 39.90420; lon2 = 116.40740; distance = calculateDistance(lat1, lon1, lat2, lon2); disp(distance); % 输出球面距离,单位:千米 ``` 运行以上代码,将会输出上海和北京之间的球面距离约为1030千米。 ### 回答2: 在Matlab中,计算两点之间的球面距离可以利用Haversine公式。该公式基于经纬度坐标系,以一个球体来近似地描述地球的形状,计算两点之间的弧长。以下是一个实现该功能的简单示例代码: ```matlab function distance = computeDistance(lat1, lon1, lat2, lon2) % 地球的平均半径(单位:千米) radius = 6371; % 将角度转换为弧度 lat1 = deg2rad(lat1); lon1 = deg2rad(lon1); lat2 = deg2rad(lat2); lon2 = deg2rad(lon2); % 计算两点间的差值 dlat = lat2 - lat1; dlon = lon2 - lon1; % 使用Haversine公式计算球面距离 a = sin(dlat/2)^2 + cos(lat1) * cos(lat2) * sin(dlon/2)^2; c = 2 * atan2(sqrt(a), sqrt(1-a)); distance = radius * c; end ``` 在上述代码中,`lat1`和`lon1`表示第一个坐标点的纬度和经度,`lat2`和`lon2`表示第二个坐标点的纬度和经度。最终,函数返回两点之间的球面距离。注意,这里使用的是地球的平均半径,该值为6371公里。 你可以调用这个函数并传入相应的经纬度数据,以计算得到两点之间的球面距离。 ### 回答3: 在Matlab中,我们可以使用Haversine公式来计算两点之间的球面距离,该公式适用于两点之间的直线距离小于200km的情况。首先,我们需要获取两个点的经纬度坐标。 假设点A的经纬度坐标为(A_lat, A_lon),点B的经纬度坐标为(B_lat, B_lon),则可以使用以下代码计算两点之间的球面距离: ```matlab function distance = calculateDistance(A_lat, A_lon, B_lat, B_lon) R = 6371; % 地球半径(单位:km) % 将角度转换为弧度 lat1 = deg2rad(A_lat); lon1 = deg2rad(A_lon); lat2 = deg2rad(B_lat); lon2 = deg2rad(B_lon); % 使用Haversine公式计算球面距离 delta_lat = lat2 - lat1; delta_lon = lon2 - lon1; a = sin(delta_lat/2)^2 + cos(lat1) * cos(lat2) * sin(delta_lon/2)^2; c = 2 * atan2(sqrt(a), sqrt(1-a)); distance = R * c; end ``` 使用上述函数即可计算两个经纬度坐标点之间的球面距离。函数参数A_lat、A_lon分别表示点A的纬度和经度,B_lat、B_lon表示点B的纬度和经度。最后,返回的距离单位为千米(km)。 例如,我们可以调用该函数进行实际计算: ```matlab A_lat = 39.9042; A_lon = 116.4074; B_lat = 31.2304; B_lon = 121.4737; distance = calculateDistance(A_lat, A_lon, B_lat, B_lon); disp(distance); % 输出两个点之间的球面距离(单位:km) ``` 通过以上示例,我们可以用Matlab计算任意两点间的球面距离。注意,如果两点之间的直线距离大于200km,我们应该使用其他更精确的球面距离计算方法。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值