Is the 3D griddata example correct?

4 次查看（过去 30 天）

Matt 2013-5-1

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/74277-is-the-3d-griddata-example-correct

I'm trying to get my head around griddata so that I can plot an intensity as the colour on a 3D surface (a hemisphere). When I start with the example on the Mathworks website "Interpolate 3-D data set over over a grid in the x-y plane" I get the same error as I do in my own code. The code copied from the Mathworks website is as follows:

rng(0,'twister')
x = 2*rand(5000,1)-1; 
y = 2*rand(5000,1)-1; 
z = 2*rand(5000,1)-1;
v = x.^2 + y.^2 + z.^2;
d = -0.8:0.05:0.8;
[xq,yq,zq] = meshgrid(d,d,0);
vq = griddata(x,y,z,v,xq,yq,zq);
surf(xq,yq,vq);
set(gca,'XTick',[-1 -0.5 0 0.5 1]);
set(gca,'YTick',[-1 -0.5 0 0.5 1]);

The subsequent error is as follows:

Error using griddata (line 51)
XI and YI must be the same size or vectors of different orientations.
Error in griddatatest (line 11)
vq = griddata(x,y,z,v,xq,yq,zq);

I presume I have copied and pasted it correctly from Mathworks Central?! I thought it would be a good place to start for my own code...

5 个评论
显示 3更早的评论隐藏 3更早的评论

Eric Sampson 2013-5-20

They did some pretty significant updates to the internals of griddata over the years, I can't remember exactly what release though. I think it does explain what you're seeing; there is likely a way to make in work in R2010 if needed, but the old version isn't as robust as the new one so it would require some fiddling around with the various griddata options.

Brady Flinchum 2013-5-21

Ya I checked it on one of the newer university computers with the 2012 release and it worked fine. After analyzing this script in more detail some it wasn't exactly what I wanted to do anyways.

Thanks for the speedy reply!

请先登录，再进行评论。

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