爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 388|回复: 4

[求助] 两组经纬度矩阵如何求得最近的索引

[复制链接]

新浪微博达人勋

发表于 2024-2-29 14:36:10 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
求助一个问题,使用python语言。我有两组经纬度,分别为a文件的经纬度(维度为375, 335):lat_a, lon_a,b文件的经纬度(维度为2030, 1354):lat_b, lon_b.设计一个程序,使得对于a文件的每一个经纬度,都找出其在b文件经纬度的位置索引。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2024-2-29 15:02:17 | 显示全部楼层
本帖最后由 zhanghm2018 于 2024-2-29 15:04 编辑

提供经纬度矩阵在附件给大家

读取数据的方法:
import numpy as np
data = np.load("data.npz")
lat_a = data["lat_a"]
lon_a = data["lon_a"]
lat_b = data["lat_b"]
lon_b = data["lon_b"]

data.npz

21.93 MB, 下载次数: 2, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2024-2-29 15:04:46 | 显示全部楼层
zhanghm2018 发表于 2024-2-29 15:02
提供经纬度矩阵在附件给大家

读取数据的方法:

如果不想花费金币下载数据的加我的qq 704806663
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2024-2-29 15:24:38 | 显示全部楼层
本帖最后由 zhanghm2018 于 2024-2-29 15:28 编辑

最容易想到的方法是遍历a文件的经纬度
计算一次的消耗时间,从而得到计算一次全部的索引耗时1小时。(个人不能接受如此长时间的计算,想找办法加速计算)

method1.png
reult1.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2024-4-14 14:40:06 | 显示全部楼层
更新解法:
使用k-d树的算法可以很快的求出多个维度的点之间的最短距离。
distance 为最小距离,index为最小距离点的索引
final_method.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表