Interpolating 3D Gridded Data in Matlab

5 次查看（过去 30 天）

Mark 2022-11-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1857663-interpolating-3d-gridded-data-in-matlab

评论： Mark 2022-11-22

I have a 3D gridded data and I want to interpolate to increase the size of the data. Current resolution is 5 deg long x 5 deg latitude and I want to change the resolution to 2 deg long x 2 deg latitude as I am interested in Lat = 2 deg. The dimension of the data is (Longitude x Latitude x Height) and current size is (72x36x40).

Long= 72x36x40
Lat = 72x36x40
Height =72x36x40
Temperature =72x36x40

I tried to use the interpn function but at the end I am only getting NaN values. Below are my codes.

file='data.nc';
Lat1 = ncread(file,'Latitude'); % Size is 72x36x40
Long1 = ncread(file1, 'Longitude'); % Size is 72x36x40
Height = ncread(file1,'Height'); % Size is 72x36x40
Temp1 = ncread(file,'Temperature'); % Size is 72x36x40
[LongI,LatI,HeightI] = ndgrid(0:5:360,-90:5:90,0:10:600);
TempI = interpn(Long1,Lat1,Height,Temp1,LongI,LatI,HeightI,'linear');`

My goal is to convert 3D data at 5 long x 5 lat resolution to 2 long x 2 lat resolution.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Matt J 2022-11-21

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1857663-interpolating-3d-gridded-data-in-matlab#answer_1106448

编辑：Matt J 2022-11-21

在 MATLAB Online 中打开

This might be easier:

LUT= griddedInterpolant(Long1,Lat1,Height,Temp1,'linear','linear'); %create interp object
newgrid=LUT.GridVectors;
newgrid(1:2)=cellfun(@(z) z(1):2:z(end) , newgrid(1:2),'uni',0); %change the sampling lattice in the first two coordinates to be increments of 2
TempI = LUT(newgrid); %interpolate on newgrid

11 个评论
显示 9更早的评论隐藏 9更早的评论

Matt J 2022-11-22

编辑：Matt J 2022-11-22

That was given in my original answer. After the upsampling to 2 long x 2 lat, the upsampled {long,lat,height} are in newgrid, while the upsampled temperatures are in TempQ.

Mark 2022-11-22