2D矢量场或3D空间中的2D矢量之间的插值

我的程序创建了在3D空间中模拟风数据的2D矢量工作表。我想知道如何在这些2D图纸之间进行插值。 X和Y值不会变化，因为它们对应于保持静态的lat/long值。矢量的W分量被设置为0，只留下矢量的U/V分量在与两个2D矢量片之间的Z距离对应的Z范围内插值。2D矢量场或3D空间中的2D矢量之间的插值

我已经读入interp3，但我不确定它会与我正在尝试做什么一起工作。如果你能用这个来帮忙，我会无限欣赏，谢谢。

注意：我删除了一个类似于我以前问过的问题，因为我相信我过分复杂了我的问题，从而为其他人解决问题提供了帮助。

感谢您的帮助，并告诉我是否可以提供更多信息！ enter image description here

来源

2013-11-24 Pierson Sargent

是否有可能为你提供一个.MAT文件中的数据被绘制以及用于生成/绘制数据的代码？ – Falimond

我为我的天真评论事先表示歉意。除了代码之外，我可以创建一个.mat文件，我认为如果你注释掉read_grib和mex文件的解压部分，你就可以做所有的事情，但我是“.mat”文件的新手，所以我不确定！让我知道！ **感谢您在无意间教我关于“.mat”文件:) 通过Dropbox下载：https://www.dropbox.com/s/5bgrwexesqm9a52/Code%26MATfile_WindSim.zip –

从您发布的图片看起来像是最终结果和插值图将看起来像箭头的波浪墙，其中上面显示的彩色片表示对于某个Z值的波浪壁的横截面。如果这是正确的，那么可能interp3不是你正在寻找的函数，而是interp2或interp1可能足以看到你的lat和long值总是固定的。然后，对于每个长/经度坐标，您只需插入一个维度（Z），然后选择多个查询点（相当于添加的横截面）。 – Falimond

这里是我的代码插入两个Z表之间的任意数量的不同站。我相信它可以被重写为在interp1（）函数中包含所有Z表单。希望这可以让你朝着正确的方向前进。

zLevels = 5; %number of interpolated points between z50 and z75 
nStation = 100;  %number of (lat,long) pairs to interpolate 
for i = 1:nStation %for nStation different (lat, long) pairs generate interp. values 
     % generate zQuery points between z50 and z75 for each station 
     zQuery = ((1:zLevels)/zLevels)*range([z50mb_vec(i) z75mb_vec(i)]) + z75mb_vec(i); 
     % use interp1 to interpolate about the Z axis for U component 
     U(i,1:zLevels) = interp1([z50mb_vec(i) z75mb_vec(i)],[UwindVec50mb(i) UwindVec75mb(i)],zQuery); 
     % and for V component 
     V(i,1:zLevels) = interp1([z50mb_vec(i) z75mb_vec(i)],[VwindVec50mb(i) VwindVec75mb(i)],zQuery); 
end 

% defining some color codes for each zLevel, otherwise the plot is a mess 
% of colors 
colorcode = ['r' 'g' 'b' 'm' 'c' 'r' 'g' 'b' 'm' 'c' 'r']; 
for j = 1:nStation 
    for i = 1:zLevels 
    quiver3(lat_value(j), long_value(j), zQuery(i), U(j,i), V(j,i), 0, colorcode(i)); 
    hold on; 
    end 
end

生成图（Z50数据是蓝绿色，Z75是红色，插片，否则）： Plot from the above code

来源

2013-11-24 06:07:11 Falimond

我想这正是我要找的人。我将不得不通过它，但我认为这太棒了。感谢您花时间浏览并帮助！非常感谢！！！ –

只有一个问题Falimond，N值代表什么？ –

应该是zLevels。原本写入时我使用N，忘记替换它。 – Falimond

2D矢量场或3D空间中的2D矢量之间的插值

回答

相关问题