比较dfs与最近的经纬度(Python,Pandas)
创始人
2024-12-13 12:31:32
0

下面是一个使用深度优先搜索(DFS)算法比较最近经纬度的示例代码,使用Python和Pandas库:

import pandas as pd
import math

def calc_distance(lat1, lon1, lat2, lon2):
    # 将经纬度转换为弧度
    lat1 = math.radians(lat1)
    lon1 = math.radians(lon1)
    lat2 = math.radians(lat2)
    lon2 = math.radians(lon2)

    # 使用球面三角法计算两个经纬度之间的距离
    dlon = lon2 - lon1
    dlat = lat2 - lat1
    a = math.sin(dlat/2)**2 + math.cos(lat1) * math.cos(lat2) * math.sin(dlon/2)**2
    c = 2 * math.atan2(math.sqrt(a), math.sqrt(1-a))
    distance = 6371 * c  # 地球半径为6371千米
    return distance

def dfs_compare_distance(df, lat, lon, visited=set()):
    visited.add((lat, lon))  # 将当前经纬度添加到已访问集合

    nearest_distance = float('inf')
    nearest_lat = None
    nearest_lon = None

    # 遍历数据框中的每个经纬度
    for index, row in df.iterrows():
        curr_lat = row['latitude']
        curr_lon = row['longitude']

        # 计算当前经纬度与目标经纬度的距离
        distance = calc_distance(lat, lon, curr_lat, curr_lon)
        
        # 如果距离更近,则更新最近经纬度和距离
        if distance < nearest_distance:
            nearest_distance = distance
            nearest_lat = curr_lat
            nearest_lon = curr_lon

        # 如果当前经纬度未访问过,则递归进行DFS比较
        if (curr_lat, curr_lon) not in visited:
            dfs_compare_distance(df, curr_lat, curr_lon, visited)

    # 打印最近经纬度和距离
    print(f"最近的经纬度:{nearest_lat}, {nearest_lon}")
    print(f"距离:{nearest_distance}千米")

# 创建一个示例数据框
data = {'latitude': [39.9042, 40.7128, 37.7749],
        'longitude': [116.4074, -74.0060, -122.4194]}
df = pd.DataFrame(data)

# 比较最近经纬度
dfs_compare_distance(df, 39.9042, 116.4074)

这段代码首先定义了一个calc_distance函数,用于计算两个经纬度之间的距离。然后定义了一个dfs_compare_distance函数,用于使用DFS算法比较最近的经纬度。该函数以给定的经纬度作为起点,遍历数据框中的每个经纬度,计算距离并更新最近经纬度和距离。同时,使用一个visited集合来记录已访问过的经纬度,避免重复访问。最后,打印出最近经纬度和距离。

在示例代码中,我们创建了一个包含三个经纬度的数据框,然后调用dfs_compare_distance函数来比较最近的经纬度。起点经纬度为(39.9042, 116.4074)(对应北京的经纬度)。程序将打印出最近的经纬度和距离。

相关内容

热门资讯

Android Recycle... 要在Android RecyclerView中实现滑动卡片效果,可以按照以下步骤进行操作:首先,在项...
安装apache-beam==... 出现此错误可能是因为用户的Python版本太低,而apache-beam==2.34.0需要更高的P...
Android - 无法确定任... 这个错误通常发生在Android项目中,表示编译Debug版本的Java代码时出现了依赖关系问题。下...
Android - NDK 预... 在Android NDK的构建过程中,LOCAL_SRC_FILES只能包含一个项目。如果需要在ND...
Akka生成Actor问题 在Akka框架中,可以使用ActorSystem对象生成Actor。但是,当我们在Actor类中尝试...
Agora-RTC-React... 出现这个错误原因是因为在 React 组件中使用,import AgoraRTC from “ago...
Alertmanager在pr... 首先,在Prometheus配置文件中,确保Alertmanager URL已正确配置。例如:ale...
Aksnginxdomainb... 在AKS集群中,可以使用Nginx代理服务器实现根据域名进行路由。以下是具体步骤:部署Nginx i...
AddSingleton在.N... 在C#中创建Singleton对象通常是通过私有构造函数和静态属性来实现,例如:public cla...
Alertmanager中的基... Alertmanager中可以使用repeat_interval选项指定在一个告警重复发送前必须等待...