Question

我在圆柱形网格上输入数据，并想在MATLAB中使用slice绘制它们。为此，我首先使用pol2cart将参考坐标转换为笛卡尔坐标。

r   = linspace(1,4,4);
phi = linspace(0,2*pi,10);
z   = linspace(1,3,3);

[rmesh,phimesh,zmesh]=meshgrid(r,phi,z)
[xmesh,ymesh,zmesh]=pol2cart(phimesh,rmesh,zmesh)

当我现在使用slice时（例如slice(xmesh,ymesh,zmesh,ones(10,4,3),2,2,2)）会抛出错误，因为坐标矩阵未正确排序（Input grid is not a valid MESHGRID.）

如何修改矩阵以获得可绘制的结果？

Answer 1

遗憾的是，您无法使用圆柱坐标中给出的数据来使用切片。

来自matlab文档：

slice(X,Y,Z,V,sx,sy,sz) draws slices of the volume V. X, Y, and Z are three-dimensional arrays specifying the coordinates for V.
X, Y, and Z must be monotonic and orthogonally spaced (as if produced by the function meshgrid).

您可以使用griddata。

这是一个例子：

    r = linspace(1,4,4);
    phi = linspace(0,2*pi,10);
    z   = linspace(1,3,3);
    data = repmat(linspace(1,0,4),[10,1,3]);

    [rmesh,phimesh,zmesh]=meshgrid(r,phi,z);
    [xmesh,ymesh,zmesh]=pol2cart(phimesh,rmesh,zmesh);

    [xg, yg, zg] = meshgrid(linspace(-4,4,50),linspace(-4,4,50),linspace(1,3,3));
    gdata = griddata(xmesh,ymesh,zmesh,data,xg,yg,zg);

    slice(xg,yg,zg,gdata,2,2,2)

根据您拥有的数据类型以及显示“超出范围”的数据（意味着，根据您的示例：半径小于1或大于4）的重要性，您可以添加以下内容：隐藏您感兴趣的区域之外的数据：

rg = sqrt(xg.^2+yg.^2);
gdataNaN = gdata;
gdataNaN(rg<min(r)) = NaN;
gdataNaN(rg>max(r)) = NaN;

figure
slice(xg,yg,zg,gdataNaN,2,2,2)

如果这还不够，你必须实现自己的slice方法（基本上使用griddata方法）或者查看matlab中心文件交换。我还没有测试过，但是用于分析MRI图像的工具可能会起作用（例如，检查：http://www.mathworks.com/matlabcentral/fileexchange/27983-3d-slicer）。

编辑： http://www.mathworks.com/matlabcentral/fileexchange/30147-smartslice-and-igslice 这似乎是由有同样问题的人开发的。

Answer 2

你能在笛卡尔坐标上调用meshgrid吗？

r   = linspace(1,4,4);
phi = linspace(0,2*pi,10);
z   = linspace(1,3,3);

[x, y] = pol2cart(r,phi);

[xmesh,ymesh,zmesh]=meshgrid(x, y, z);

有时您可能会发现[ymesh,xmesh,zmesh]=meshgrid(x, y, z);对您的应用程序更有意义。

如何使用切片在圆柱形网格上绘制数据？

2 个答案: