根据经纬度自动寻找最近的点，并返回相应属性

最新推荐文章于 2024-06-29 21:27:33 发布

hydromonkey

最新推荐文章于 2024-06-29 21:27:33 发布

阅读量2.7k

点赞数

分类专栏：私人代码记录

本文链接：https://blog.csdn.net/hengcall/article/details/82807946

版权

本文介绍如何通过经纬度坐标找到最近的地理位置，并从中检索相关属性，适用于地理信息系统（GIS）和导航应用的技术实现。

摘要由CSDN通过智能技术生成

import pandas as pd
import numpy as np
from scipy.spatial.distance import cdist
import scipy.stats as stats
from pylab import *
STA_tmin = 'F:/AAA/STA-T-TMAX-TMIN-CSV/STA-TMEAN.csv'
PG_tmin = 'F:/AAA/GRID-TEM-CSV/IF-MEAN/PGTMEAN.csv'
sta_tmin = pd.read_csv(STA_tmin)
pg_tmin = pd.read_csv(PG_tmin)



# 两组数据
def closest_point(point, points): # 根据哪列，在哪里找
    """ Find closest point from a list of points. """
    return points[cdist([point], points).argmin()]  # 前者必须是个array

def match_value(df, col1, x, col2):  # 哪里找，根据什么找，和谁匹配,返回什么值。
    """ Match value x from col1 row to value in col2. """
    return df[df[col1] == x][col2].values[0]

pg_tmin['point'] = [(x, y) for x, y in zip(pg_tmin['lat'], pg_tmin['lon'])]
sta_tmin['point'] = [(x, y) for x, y in zip(sta_tmin['lat'], sta_tmin['lon&#