已知两点经纬度计算距离

2024-06-12 06:58

本文主要是介绍已知两点经纬度计算距离,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

已知两点经纬度计算距离

该博文来自于ieayoio的博客:http://www.ieayoio.com/

最近做了一些计算最近附近商店距离的工作,距离计算原以就是勾股定理那么简单,然而真去做的时候才知道太天真了,经线上纬度每差一度可以有一个确定的值,然而,因为纬线圈的长度随着纬度的改变而改变,不能单纯的计算每个单位经度上距离的一个确定值。
经过一番的查找,原来计算距离有一个公式如下:
S = 2 a r c s i n s i n 2 a 2 + c o s ( L a t 1 ) × c o s ( L a t 2 ) × s i n 2 b 2 × 6378.137 S = 2 arcsin\sqrt{sin^2\frac a2 + cos(Lat1) × cos(Lat2) × sin^2\frac b2} × 6378.137 S=2arcsinsin22a+cos(Lat1)×cos(Lat2)×sin22b ×6378.137

  1. Lat1 Lung1 表示A点经纬度,Lat2 Lung2 表示B点经纬度,这里的单位是弧度而不是角度;

  2. a=Lat1 – Lat2 为两点纬度之差 b=Lung1 -Lung2 为两点经度之差;

  3. 6378.137为地球半径,单位为千米;

从公式中可以看到两点的经度Lung1和Lung2没有包含在公式中,而只是包含在差值之中。

知道公式了,网上的各种教程无非就是套用公式,例如数据库要查询可以写一个自定义的函数,一下以mysql为例:

DROP FUNCTION IF EXISTS GetDistance;
DELIMITER $$
CREATE FUNCTION GetDistance(LatBegin DOUBLE, LngBegin DOUBLE, LatEnd DOUBLE, LngEnd DOUBLE)
RETURNS DOUBLE  
BEGIN  DECLARE Distance DOUBLE;DECLARE EARTH_RADIUS DOUBLE;DECLARE RadLatBegin,RadLatEnd,RadLatDiff,RadLngDiff DOUBLE; SET EARTH_RADIUS = 6378.137;SET RadLatBegin = LatBegin * PI()/180.0;SET RadLatEnd = LatEnd * PI()/180.0;SET RadLatDiff = RadLatBegin - RadLatEnd;SET RadLngDiff = LngBegin * PI()/180.0 - LngEnd * PI()/180.0;   SET Distance = 2 *ASIN(SQRT(POWER(SIN(RadLatDiff/2), 2)+COS(RadLatBegin)*COS(RadLatEnd)*POWER(SIN(RadLngDiff/2), 2)));SET Distance = Distance * EARTH_RADIUS;RETURN Distance;
END  $$
DELIMITER ;

这样定义一个自定义函数之后就可以直接利用这个函数进行数据库的查询和排序,例如:

SELECT * FROM `shop` WHERE GetDistance(`double_latitude`,`double_longitude`,34.54,234.4) < 5 ORDER BY GetDistance(`double_latitude`,`double_longitude`,34.54,234.4) DESC LIMIT 1,20

另外再提供一个Java版本:

import static java.lang.Math.*;/*** Created by ieayoio on 16-12-2.*/
public class DistanceUtils {public static double getDistance(double LatBegin, double LngBegin, double LatEnd, double LngEnd) {double EARTH_RADIUS = 6378.137;double RadLatBegin = LatBegin * PI / 180.0;double RadLatEnd = LatEnd * PI / 180.0;double RadLatDiff = RadLatBegin - RadLatEnd;double RadLngDiff = LngBegin * PI / 180.0 - LngEnd * PI / 180.0;double Distance = 2 * asin(sqrt(pow(sin(RadLatDiff / 2), 2) + cos(RadLatBegin) * cos(RadLatEnd) * pow(sin(RadLngDiff / 2), 2)));return Distance * EARTH_RADIUS;}
}

本文链接:http://www.ieayoio.com/2016/12/02/已知经纬度计算距离/

参考资料:
https://zhidao.baidu.com/question/401287294.html
http://blog.csdn.net/smartsmile2012/article/details/45339787
http://blog.csdn.net/xiyang_1990/article/details/16803735

这篇关于已知两点经纬度计算距离的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1053506

相关文章

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显

计算数组的斜率,偏移,R2

模拟Excel中的R2的计算。         public bool fnCheckRear_R2(List<double[]> lRear, int iMinRear, int iMaxRear, ref double dR2)         {             bool bResult = true;             int n = 0;             dou

GPU 计算 CMPS224 2021 学习笔记 02

并行类型 (1)任务并行 (2)数据并行 CPU & GPU CPU和GPU拥有相互独立的内存空间,需要在两者之间相互传输数据。 (1)分配GPU内存 (2)将CPU上的数据复制到GPU上 (3)在GPU上对数据进行计算操作 (4)将计算结果从GPU复制到CPU上 (5)释放GPU内存 CUDA内存管理API (1)分配内存 cudaErro

Java - BigDecimal 计算分位(百分位)

日常开发中,如果使用数据库来直接查询一组数据的分位数,就比较简单,直接使用对应的函数就可以了,例如:         PERCENT_RANK() OVER(PARTITION BY 分组列名 ORDER BY 目标列名) AS 目标列名_分位数         如果是需要在代码逻辑部分进行分位数的计算,就需要我们自己写一个工具类来支持计算了 import static ja

OpenStack离线Train版安装系列—2计算节点-环境准备

本系列文章包含从OpenStack离线源制作到完成OpenStack安装的全部过程。 在本系列教程中使用的OpenStack的安装版本为第20个版本Train(简称T版本),2020年5月13日,OpenStack社区发布了第21个版本Ussuri(简称U版本)。 OpenStack部署系列文章 OpenStack Victoria版 安装部署系列教程 OpenStack Ussuri版

新一代车载(E/E)架构下的中央计算载体---HPC软件架构简介

老规矩,分享一段喜欢的文字,避免自己成为高知识低文化的工程师: 屏蔽力是信息过载时代一个人的特殊竞争力,任何消耗你的人和事,多看一眼都是你的不对。非必要不费力证明自己,无利益不试图说服别人,是精神上的节能减排。 无人问津也好,技不如人也罢,你都要试着安静下来,去做自己该做的事.而不是让内心的烦躁、焦虑、毁掉你本就不多的热情和定力。 时间不知不觉中,快要来到夏末秋初。一年又过去了一大半,成