3

我正在制作一个泛化圆柱体功能的函数，以便圆柱体有盖，可以是任何尺寸和方向。然而，从气缸的外观来看，我遇到了堵塞。为了使盖子看起来正确，弯曲部分需要一组阴影，而盖子需要另一组。（在您要求制作 3 个表面之前不是一个选项）

以下是相关代码：

    surface(xSurf,ySurf,zSurf,c,'EdgeColor','none','FaceLighting','phong');

如果您想查看整个代码。

感谢您的帮助，

约翰

    function varargout = DrawCylinder(x,y,z,r,h,aVec,bVec,cVec,ccolor, npts)
    % DrawCylinder Generate a three-dimensional cylinder
    %
    % DrawCylinder(x,y,z,a,b,c,aVec,bVec,CVec,ccolor, npts) 
    % creates a surface plot of a cylinder whose center is at (x,y,z), has
    % semiaxes of length a, b, and c.  The unit vectors associated with each
    % semixis are aVec, bVec, and cVec and must be size 3 x 1 (column vector)
    % with size of npts + 1.
    %
    % H = DrawCylinder(...) creates the surface plot and returns the handle H to each
    % graphical object created.
    %
    % [X Y Z] = DrawCylinder(...) does not generate the surface plot and returns
    % the data necessary to create the surface using:
    % SURF(X,Y,Z);
    %
    % [X Y Z C] = DrawCylinder(...) does not generate the surface plot and returns
    % the data necessary to create the surface using:
    % SURF(X,Y,Z,C,'cdataMapping','direct');

    %CREATE SURFACE FOR CYLINDER
    [xCyl,yCyl,zCyl]=cylinder(1,npts);

    xSurf=[zeros(1,max(size(xCyl)));xCyl;zeros(1,max(size(xCyl)))];
    ySurf=[zeros(1,max(size(yCyl)));yCyl;zeros(1,max(size(yCyl)))];
    zSurf=[zeros(1,max(size(zCyl)));zCyl;ones(1,max(size(zCyl)))] - 0.5;

    xSurf = xSurf*r;
    ySurf = ySurf*r;
    zSurf = zSurf*h;

    %ROTATE CYLINDER
    %Make sure aVec,bVec, and cVec are column unit vectors:
    if all(size(aVec)==[1,3])
        aVec=aVec';
    end
    if all(size(bVec)==[1,3])
        bVec=bVec';
    end
    if all(size(cVec)==[1,3])
        cVec=cVec';
    end
    aVec=aVec/norm(aVec); %Make unit vectors
    bVec=bVec/norm(bVec);
    cVec=cVec/norm(cVec);

    rot = [aVec,bVec,cVec]; %The rotation matrix

    [iMax, jMax] = size(xSurf);

    for i=1:iMax
        for j=1:jMax
            rotatedPt = rot*[xSurf(i,j);ySurf(i,j);zSurf(i,j)];
            xSurf(i,j) = rotatedPt(1);
            ySurf(i,j) = rotatedPt(2);
            zSurf(i,j) = rotatedPt(3);
        end
    end

    %TRANSLATE CYLINDER
    xSurf = xSurf + x;
    ySurf = ySurf + y;
    zSurf = zSurf + z;

    c = ccolor*ones(size(xSurf));
    if nargout == 0
        surface(xSurf,ySurf,zSurf,c,'EdgeColor','none','FaceLighting','phong');
    elseif nargout == 1
        varargout = {surface(xSurf,ySurf,zSurf,c,'EdgeColor','none','FaceLighting','phong');};
    elseif nargout == 3
        varargout = {xSurf ySurf zSurf};
    elseif nargout == 4
        varargout = {xSurf ySurf zSurf c};
    end

    end

2012 年 8 月 18 日编辑：

只是为了让你能看到。

这就是我得到的...

这就是我得到的

这就是我想要的...

这就是我想要的

4

1 回答 1

1

我很确定“阴影”只是指端盖的颜色，而不是更复杂的效果。如果是这样，那很简单，我将举一个灰度示例

改变

c = ccolor*ones(size(xSurf));

到

c = (ccolor/255)*ones(size(xSurf));
c([1 3],:)=max(0,(ccolor-10))/255;

第一行c用归一化的 ccolor 初始化矩阵（期望 8 位灰度 ccolor 输入，它归一化为 0..1）。第二行将帽（第 1 行和第 3 行）更改为稍深的颜色，在 0 处触底，并单独保留圆柱体表面（第 2 行和第 4 行）。

为了确保正确查看结果，您需要更改 nargout==0 条件，使其看起来像这样

surface(xSurf,ySurf,zSurf,c,'EdgeColor','none','FaceLighting','phong');
colormap(gray(256));
caxis([0 1]);

颜色图只是设置颜色图，类似于 8 位灰度。caxis 命令相当关键。根据 Matlab 的表面文档

MATLAB 对此数据执行线性变换以从当前颜色图中获取颜色

为了我们的目的，这是不好的。由于我们只有两个值，最低值将更改为 0，最高值将更改为 1。这实际上忽略了我们的 ccolor 输入，并给出了一个带有两个黑色大写字母的白色圆柱体。使用caxis([0 1])保留完整比例和 ccolor 在其中的位置。

更新：

听起来我误解了你想要什么，实现非常接近你想要的效果的最简单方法是将'MeshStyle'设置为'row'，如下所示：

surface(xSurf,ySurf,zSurf,c,'EdgeColor','k','FaceLighting','phong','MeshStyle','row');

这将为您提供以下结果：加盖气缸

仍然有一个中心点，但它是迄今为止产生这种效果的最简单的方法。

于 2012-08-18T17:08:53.283 回答