python实现matlab中interp2方法(有误差（小数点后第5位）)

浮麈

已于 2022-06-29 17:51:20 修改

阅读量1k

点赞数

分类专栏： python 文章标签： python matlab

于 2022-06-29 11:35:20 首次发布

本文链接：https://blog.csdn.net/weixin_43865446/article/details/125517106

版权

python 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

在matlab中的interp2方法

lat = ds1.data('lat')'; % 1x451 90:-0.4:-90
lon = ds1.data('lon')'-360; % 1x900 -180:0.4:179.6

%% 插值,从0.4度分辨率插值到0.125度
% 转换后
%lat 90:-0.125:-90 1x1441
%lon 0:0.125:359.875 1x2880
%U: 2880x1441 lon x lat

lat1 = 90:-0.125:-90; %1x1441
lon1 = 0:0.125:359.875; %1x2880
[LAT,LON] = meshgrid(lat,lon);
[LAT1,LON1] = meshgrid(lat1,lon1);
U1 = interp2(LAT,LON,U,LAT1,LON1); %lon x lat 2880x1441

如图：

在pyhon中实现

import numpy as np
from scipy import interpolate

#将矩阵的行进行倒叙
def invert(list1):
    return np.squeeze([row[::-1] for row in list1])




#示例
def test
    U = []#矩阵
    lat = np.transpose(ds['lat'].data)# 1x451 90:-0.4:-90
    lon = np.transpose(ds['lon'].data) # 1x900 -180:0.4:179.6

    lat1 = np.arange(90, - 90 + - 0.125, - 0.125)
    lon1 = np.arange(0, 359.875 + 0.125, 0.125)
    #lat、lon、lat1、lon1在interp2d方法下实现插值，但是行的顺序是颠倒的
    f = interpolate.interp2d(lat, lon, U)
    U1 = f(lat1, lon1)
    #将矩阵行中数据颠倒过来，就得到和matlab中和interp2d保持一致的数据了
    re= invert(U1)