在Python中将纬度经度数据reshape 为二维矩阵

发布于08月08日

我有一个包含regularly spaced个数据的数组，格式为lats,lons,values.我的lats的步长为1，从-90到90不等，我的经度从-180到180，步长为1.

我想生成一个矩阵，其中除了在纬度和经度上定义的值之外，所有值都是0.做这件事最好的方法是什么？我的try 是:

lats=np.array([-85,-43,2])
lons=np.array([-44,-66,24])
values=np.array([1,2,3]) #some random data and points

matrix=np.zeros((181,361))

for i in lats:
  for j in lons:
    matrix[i][j]=values[]

但我不确定要对变量values使用什么索引.

推荐答案

我们知道指数(0，0)是纬度-90，经度-180.要将经度和纬度值转换为索引，请向纬度值添加90，向经度值添加180.这违反了您的规范，即坐标都将是整数.

示例代码:

import numpy as np

lats = np.array([-85,-43,2])
lons = np.array([-44,-66,24])
values = np.array([1,2,3])

coords = np.stack((lats, lons))
offset = np.array([90, 180])
indices = coords + offset[:, None]

matrix = np.zeros((181,361))
matrix[*indices] = values