MATLAB流带边缘颜色

Question

在流带图中，如何绘制丝带本身的边缘，而不是每个面的边缘？

例如：

figure
load wind
[sx sy sz] = meshgrid(80,20:10:50,0:5:15);
h = streamribbon(x,y,z,u,v,w,sx,sy,sz);
axis tight
shading interp;
view(3);
camlight
lighting gouraud

产生：

如果我添加set(h,'edgecolor','k')

我得到：

而我正在寻找的是这样的东西（图像的顶部和左侧已经用 gimp 进行了编辑，以表明我正在追逐的东西）：

Answer 1

正如评论中所说，没有简单的方法可以照亮色带表面的某些边缘而不是其他边缘。

但是，带状表面仅由边缘上的点定义，因此只需检索这些点并重新组织（内联）它们即可获得 3D 线的完美坐标集。
要获得一个封闭配置文件，我们只需要复制配置文件末尾的第一个点。

我将它打包到一个函数中，该函数将返回一组与功能区句柄数组大小相同的handles，并将您可以发送的任何其他参数传输到内部plot3 函数。所以用法很简单：

所以使用与问题完全相同的代码（检索变量h中的色带表面句柄数组），添加：

hs = StreamEdges( h , 'k','Linewidth',1 ) ;

将产生：

当然，你总是可以疯狂地在组中设置你的线属性：

set(hs,'LineWidth',1,'Color','c')
%// or
set(hs,'LineWidth',2,'Color','r','Marker','none')

---

函数StreamEdges.m的代码：

function hs = StreamEdges(hribbon,varargin)

wasOnHold = ishold ;            %// save hold state
hold on

hs = zeros( size(hribbon) ) ;   %// initialize output handle array
for ih=1:numel( hribbon ) ;
    %// retrieve X, Y, and Z data of each surface
    hsurf = handle( hribbon(ih)) ;  %// for pre-HG2 (2014a) versions
    xx = hsurf.XData ;
    yy = hsurf.YData ;
    zz = hsurf.ZData ;

    %// reoder coordinates to linearise surface profile (edge)
    xx = [xx(:,1) ; flipud(xx(:,2)) ; xx(1)] ;
    yy = [yy(:,1) ; flipud(yy(:,2)) ; yy(1)] ;
    zz = [zz(:,1) ; flipud(zz(:,2)) ; zz(1)] ;

    %// return double to stay compatible with pre-HG2
    hs(ih) = double( plot3(xx(:),yy(:),zz(:), varargin{:} ) ) ;
end

if ~wasOnHold ; hold off ; end      %// restore hold state
hs = handle(hs) ;                   %// convert "double" handle back to HG2 handles

---

该函数应适用于 HG2 之前的版本（=2014b 之前）或 HG2 之后的版本。它将返回的句柄数组将采用 HG2 格式（= 真正的句柄，而不仅仅是它们的 double 表示形式）。这使您可以访问图形对象上的点符号，即使对于 HG2 之前的版本也是如此。如果您对此不满意，只需注释函数的最后一行（或者要获得旧的句柄样式，您可以随时使用该函数 hs = double(hs)）

---

目标axes

除了保留hold 状态外，该函数的行为类似于大多数高级绘图函数，这意味着它们将直接转到当前的axes (gca)，或者如果不存在则创建一个新的.
您可以轻松地修改函数以接受axes 句柄作为目标，或者直接在不修改的情况下在参数中显式调用它：

hs = StreamEdges( h , 'Parent','target_axes_handle' , 'Color','k','Linewidth',1) ;