如何根据经纬度精确计算两点之间的距离

以下代码用python2运行

#!/usr/bin/python
#coding=utf-8
from __future__ import print_function
import sys
reload(sys)
sys.setdefaultencoding("utf-8")
import mathdef haversine_distance(lon1, lat1, lon2, lat2):# 将纬度和经度从度转换为弧度lat1_rad = math.radians(lat1)lon1_rad = math.radians(lon1)lat2_rad = math.radians(lat2)lon2_rad = math.radians(lon2)# 计算差值dlon = lon2_rad - lon1_raddlat = lat2_rad - lat1_rad# 应用 Haversine 公式a = math.sin(dlat / 2)**2 + math.cos(lat1_rad) * math.cos(lat2_rad) * math.sin(dlon / 2)**2c = 2 * math.atan2(math.sqrt(a), math.sqrt(1 - a))distance = 6371 * c * 1000  # 地球平均半径约为 6371 千米return distancedef vincenty(lon1, lat1, lon2, lat2):# 定义椭球体参数a = 6378137.0  # 长半轴(单位:米)f = 1 / 298.257223563  # 扁率b = (1 - f) * a  # 短半轴(单位:米)# 将经纬度从度数转换为弧度lon1, lat1, lon2, lat2 = map(math.radians, [lon1, lat1, lon2, lat2])# 差值L = lon2 - lon1U1 = math.atan((1 - f) * math.tan(lat1))U2 = math.atan((1 - f) * math.tan(lat2))sinU1, cosU1 = math.sin(U1), math.cos(U1)sinU2, cosU2 = math.sin(U2), math.cos(U2)# 初始设定lamb = Liter_limit = 100  # 最大迭代次数for _ in range(iter_limit):sin_lamb = math.sin(lamb)cos_lamb = math.cos(lamb)sin_sigma = math.sqrt((cosU2 * sin_lamb) ** 2 + (cosU1 * sinU2 - sinU1 * cosU2 * cos_lamb) ** 2)if sin_sigma == 0:return 0  # 重合点cos_sigma = sinU1 * sinU2 + cosU1 * cosU2 * cos_lambsigma = math.atan2(sin_sigma, cos_sigma)sin_alpha = cosU1 * cosU2 * sin_lamb / sin_sigmacos2_alpha = 1 - sin_alpha ** 2cos2_sigma_m = cos_sigma - 2 * sinU1 * sinU2 / cos2_alphaif math.isnan(cos2_sigma_m):cos2_sigma_m = 0  # 赤道情况C = f / 16 * cos2_alpha * (4 + f * (4 - 3 * cos2_alpha))lamb_prev = lamblamb = L + (1 - C) * f * sin_alpha * (sigma + C * sin_sigma * (cos2_sigma_m + C * cos_sigma * (-1 + 2 * cos2_sigma_m ** 2)))if abs(lamb - lamb_prev) < 1e-12:breakelse:raise ValueError("Vincenty formula failed to converge")u2 = cos2_alpha * (a ** 2 - b ** 2) / (b ** 2)A = 1 + u2 / 16384 * (4096 + u2 * (-768 + u2 * (320 - 175 * u2)))B = u2 / 1024 * (256 + u2 * (-128 + u2 * (74 - 47 * u2)))delta_sigma = B * sin_sigma * (cos2_sigma_m + B / 4 * (cos_sigma * (-1 + 2 * cos2_sigma_m ** 2) -B / 6 * cos2_sigma_m * (-3 + 4 * sin_sigma ** 2) *(-3 + 4 * cos2_sigma_m ** 2)))s = b * A * (sigma - delta_sigma)return s  # 距离,单位:米# 示例:计算北京和上海之间的距离
beijing_lon, beijing_lat = 116.407396, 39.904200
shanghai_lon, shanghai_lat = 121.473701, 31.230416distance = vincenty(beijing_lon, beijing_lat, shanghai_lon, shanghai_lat)
print("北京和上海之间的距离是 {:.2f} 米".format(distance))

haversine算法与vincenty的差距很小,对于500多米的近距离,误差在1米多。

将haversine用sql语句表达结果为:

select 6371*1000 * acos(sin(radians(a.latitude)) * sin(radians(b.latitude)) + cos(radians(a.latitude)) * cos(radians( b.latitude)) * cos(radians(b.longitude) - radians(a.longitude))) AS distance 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/pingmian/40864.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

【C#】函数方法、属性分文件编写

1.思想 分文件编写是面向对象编程的重要思想&#xff0c;没有实际项目作为支撑很难理解该思想的精髓&#xff0c;换言之&#xff0c;一两个函数代码量因为太少无法体现分文件编写减少大量重复代码的优势。 2.项目结构介绍 整项目的名称叫AutoMetadata&#xff0c;是一个基于W…

JAVA:文件防重设计指南

1、简述 在现代应用程序中&#xff0c;处理文件上传是一个常见的需求。为了保证文件存储的高效性和一致性&#xff0c;避免重复存储相同的文件是一个重要的优化点。本文将介绍一种基于哈希值的文件防重设计&#xff0c;并详细列出实现步骤。 2、设计原理 文件防重的基本思路…

Dns被莫名篡改的逆向分析定位(笔记)

引言&#xff1a;最近发现用户的多台机器上出现了Dns被莫名修改的问题&#xff0c;从系统事件上看并未能正常确定到是那个具体软件所为&#xff0c;现在的需求就是确定和定位哪个软件具体所为。 解决思路&#xff1a; 首先到IPv4设置页面对Dns进行设置&#xff1a;通过ProcExp…

setjmp和longjmp函数使用

这里用最简单直接的描述&#xff1a;这两组函数是用于实现类似vscode全局的标签跳转功能&#xff0c;setjmp负责埋下标签&#xff0c;longjmp负责标签跳转。 #include <stdio.h> #include <stdlib.h> #include <setjmp.h>jmp_buf envbuf1; jmp_buf envbuf2;…

HttpServer内存马

HttpServer内存马 基础知识 一些基础的方法和类 HttpServer&#xff1a;HttpServer主要是通过带参的create方法来创建&#xff0c;第一个参数InetSocketAddress表示绑定的ip地址和端口号。第二个参数为int类型&#xff0c;表示允许排队的最大TCP连接数&#xff0c;如果该值小…

【面试系列】软件工程师高频面试题及详细解答

欢迎来到我的博客&#xff0c;很高兴能够在这里和您见面&#xff01;欢迎订阅相关专栏&#xff1a; ⭐️ 全网最全IT互联网公司面试宝典&#xff1a;收集整理全网各大IT互联网公司技术、项目、HR面试真题. ⭐️ AIGC时代的创新与未来&#xff1a;详细讲解AIGC的概念、核心技术、…

JAVA每日作业day7.4

ok了家人们今天学习了Date类和simpleDateformat类&#xff0c;话不多说我们一起看看吧 一.Date类 类 java.util.Date 表示特定的瞬间 ( 日期和时间 ) &#xff0c;精确到毫秒。 1.2 Date类的构造方法 public Date(): 用来创建当前系统时间对应的日期对象。 public Date(long …

【java开发环境】多版本jdk 自由切换window和linux

win10 一、准备 各种版本的jdk&#xff0c;按自己的需要下载。 我这里是需要jdk17和jdk8。 1、jdk17 下载&#xff1a;Java Downloads | Oracle&#xff0c;选择exe后缀文件 2、jdk8下 载&#xff1a;Java Downloads | Oracle&#xff0c;选择exe后缀文件 二、详细步骤 1、…

Linux线程:编织并发的梦幻世界

目录 &#x1f6a9;引言 &#x1f6a9;听故事&#xff0c;引概念 &#x1f6a9;生产者消费者模型 &#x1f680;再次理解生产消费模型 &#x1f680;挖掘特点 &#x1f6a9;条件变量 &#x1f680;条件变量常用接口 &#x1f680;条件变量的原理 &#x1f6a9;引言 上一篇…

Redis分布式锁代码实现详解

引言 在分布式系统中&#xff0c;资源竞争和数据一致性问题常常需要通过锁机制来解决。Redis作为一个高性能的键值存储系统&#xff0c;因其提供的原子操作、丰富的数据结构以及网络延迟低等特点&#xff0c;成为了实现分布式锁的理想选择。本文将详细介绍如何使用Redis来实现…

打卡第2天----数组双指针,滑动窗口

今天是参与训练营第二天&#xff0c;这几道题我都看懂了&#xff0c;自己也能写出来了&#xff0c;实现思路很重要&#xff0c;万事开头难&#xff0c;希望我可以坚持下去。希望最后的结果是量变带来质变。 一、理解双指针思想 leetcode编号&#xff1a;977 不止是在卡尔这里…

一篇文章说清楚Filter(过滤器)、Interceptor(拦截器)和AOP(切面儿)

文章目录 前言一、Filter&#xff08;过滤器&#xff09;1.说明2.实现filterChain.doFilter() 3.order优先级4.解决跨域5.拦截返回错误信息JSON 二、Interceptor&#xff08;拦截器&#xff09;1.说明2.实现preHandlepostHandleafterCompletion 3.执行顺序图4.排除特定路径拦截…

【LinuxC语言】线程退出函数的实现

文章目录 前言线程退出函数的作用和实现流程函数实现线程池代码概况总结前言 在并发编程中,线程的创建和销毁是两个非常重要的环节。在Linux C语言环境下,我们可以使用POSIX线程库(pthread)来创建和控制线程。然而,线程的退出并不是一个简单的过程,它涉及到许多细节,如…

论文学习——基于类型检测的动态自适应多目标优化算法

论文题目&#xff1a;Dynamic adaptive multi-objective optimization algorithm based on type detection 基于类型检测的动态自适应多目标优化算法&#xff08;Xingjuan Cai a,b, Linjie Wu a,∗, Tianhao Zhao a, Di Wu c, Wensheng Zhang d, Jinjun Chen e&#xff09;Inf…

yum命令提示 错误:rpmdb: BDB0113 Thread/process 4153/139708200269632

一、报错信息 [rootDawn yum.repos.d]# yum clean all 错误&#xff1a;rpmdb: BDB0113 Thread/process 4153/139708200269632 failed: BDB1507 Thread died in Berkeley DB library 错误&#xff1a;db5 错误(-30973) 来自 dbenv->failchk&#xff1a;BDB0087 DB_RUNRECOVE…

HTTP与HTTPS协议区别及应用场景

HTTP&#xff08;超文本传输​​协议&#xff09;和 HTTPS&#xff08;安全超文本传输​​协议&#xff09;都是用于通过网络传输数据的协议。虽然它们有一些相似之处&#xff0c;但在安全性和数据保护方面也存在显著差异。 在这篇博文中&#xff0c;我们将探讨 HTTP 和 HTTPS…

websocket的工作原理

WebSocket 是一种在 Web 应用程序中实现全双工通信的协议。它允许客户端和服务器之间建立一个持久的连接&#xff0c;双方可以在这个连接上相互发送数据&#xff0c;而不需要通过传统的 HTTP 请求/响应模式。这使得 WebSocket 特别适用于实时应用&#xff0c;如聊天应用、在线游…

pytorch中的contiguous()

官方文档&#xff1a;https://pytorch.org/docs/stable/generated/torch.Tensor.contiguous.html 其描述contiguous为&#xff1a; Returns a contiguous in memory tensor containing the same data as self tensor. If self tensor is already in the specified memory forma…

自然语言处理与Transformer模型:革新语言理解的新时代

引言 自然语言处理&#xff08;NLP&#xff09;是人工智能和计算机科学的一个重要分支&#xff0c;旨在使计算机能够理解、生成和处理人类语言。随着互联网和数字化信息的爆炸性增长&#xff0c;NLP在许多领域中的应用变得越来越重要&#xff0c;包括&#xff1a; 搜索引擎&am…

domain 和 特征之间的关系

在计算机学术届&#xff0c;通常说的domain是什么意思&#xff1f; 例子1&#xff1a;在数学和函数的语境中&#xff0c;“domain” 通常指的是函数的定义域&#xff0c;即函数可以接受的输入值的集合。 假设这个时候定义域是{1,3,4,6,7,8,9,55}&#xff0c;那么定义域里面的…