Python网格数据升尺度降尺度

处理地理数据经纬度精度不匹配问题

降尺度

from scipy.interpolate import interp2d
import numpy as np
import matplotlib.pyplot as plt

#数据准备，x y 为网格行列格点数，以1为一格范围
x=np.arange(0,10)
y=np.arange(0,10)
#生成10*10网格数据
data = np.arange(100).reshape((10,10))
plt.imshow(data)
plt.colorbar()

# 创建新的格点行列数，以0.5为一格范围
x1 = np.arange(0,10,0.5)
y1 = np.arange(0,10,0.5)

#线性插值 或更改kind采取非线性插值
f1 = interp2d(x,y,data,kind='linear')
data2 = f1(x1,y1)

#画图
plt.imshow(data2)
plt.colorbar()

升尺度

x1 = np.arange(0,10,2)
y1 = np.arange(0,10,2)
f1 = interp2d(x,y,data,kind='linear')
data2 = f1(x1,y1)
plt.imshow(data2)
plt.colorbar()