Error using interp3 "Error using griddedInterpolant. Grid arrays must have NDGRID structure"

21 次查看（过去 30 天）

KCE 2023-8-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2013727-error-using-interp3-error-using-griddedinterpolant-grid-arrays-must-have-ndgrid-structure

编辑： Bruno Luong 2023-8-29

I am trying to interpolate observation heights from pressures in 3D space. I am using interp3 because it is generally faster than griddata, however I return this error, and I am not really sure what NDGRID structure means. I am running it in a loop as well. The following sizes of variables in interp3(X,Y,Z,V,Xq,Yq,Zq):

for ns=1:length(obs_is)
    obs_hts(ns)=interp3(repmat(y_pts,[1 1 size(h,3)]),repmat(x_pts,[1 1 size(h,3)]),double(p),double(h),obs_lats(ns),obs_lons(ns),obs_pres(ns));
end

X: latitudes in 3D model space (199x199x50)

Y: longitudes in 3D model space (199x199x50)

Z: pressure in 3D model space (199x199x50)

V: heights in 3D model space (199x199x50)

Xq: observation latitude (scalar)

Vq: observation longitude (scalar)

Zq: observation pressure (scalar)

And yes, these are all double (not single) variables.

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

请先登录，再进行评论。

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

回答（2 个）

Star Strider 2023-8-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2013727-error-using-interp3-error-using-griddedinterpolant-grid-arrays-must-have-ndgrid-structure#answer_1294542

The ndgrid function creates gridded data (in 2, 3, ..., n) dimehsions from its argument vectors.

Considering the nature of your problem, the interpn function might be more appropriate.

10 个评论
显示 8更早的评论隐藏 8更早的评论

Walter Roberson 2023-8-28

interp3() has a try/catch that tests the error generated by griddedInterpolant() and conditionally substitutes a different error message referencing meshgrid(). But perhaps the error codes got updated while the interp3() error code test did not get updated?

Paul 2023-8-28

在 MATLAB Online 中打开

That sounds like an odd way to handle input argument checking to interp3. Seems like that should be done at the interp3 level before the call to griddedInterpolant.

At least calling interp3 with ndgrid inputs generates the expected error

[X,Y,Z] = ndgrid(1:4);
V = X + Y + Z;
interp3(X,Y,Z,V,1.5,2.5,3.5)
Error using interp3
Input grid is not a valid MESHGRID.

请先登录，再进行评论。

Bruno Luong 2023-8-28

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2013727-error-using-interp3-error-using-griddedinterpolant-grid-arrays-must-have-ndgrid-structure#answer_1295266

编辑：Bruno Luong 2023-8-29

在 MATLAB Online 中打开

Check with this

for ns=1:length(obs_is)
    X = repmat(y_pts,[1 1 size(h,3)]);
    Y = repmat(x_pts,[1 1 size(h,3)]);
    Z = double(p);
    if ~all(X == X(1,:,1), 'all')
        if all(X == X(:,1,1), 'all') && all(Y == Y(1,:,1), 'all')
            warning('You might swap X and Y')
        else
            error('X for ns=%d X is not meshgrid generated', ns);
        end
    end
    if ~all(Y == Y(:,1,1), 'all')
        error('Y for ns=%d Y is not meshgrid generated', ns);
    end
    if ~all(Z == Z(1,1,:), 'all')
        error('Z for ns=%d Z is not meshgrid generated', ns);
    end
end