如何在表面的脊上绘制点？

问题描述

我找不到在表面的上下脊上绘制一些点的方法。 我为表面编写了这段代码：

%initialize parameters
x   = linspace(-2,2,100);
y   = linspace(-2,80);
fxy = abs(log(x'.*y));


% it looks like this 
figure(1),clf
surf(x,y,fxy')
shading interp
axis square,rotate3d on
xlabel('X'),ylabel('y'),zlabel('f(x,y)')

为了找到较低的山脊，我写了类似下面的内容：

% find min

minval = min(min(fxy));
[xi,yi] = find( minval == fxy);


idx = sub2ind(size(fxy),xi,yi);

% plot the minimum as a red ball
hold on
plot3(x(xi),y(yi),fxy(idx)','ro','MarkerFaceColor','r','MarkerSize',12 )

然后针对下限（较低的山​​脊）设置阈值0.1，然后：

% finding lower ridge : points below a threshold of .01

[xi,yi] = find( fxy < minval+.01);

% to get the values,convert from matrix indices to linear indices
idx = sub2ind(size(fxy),yi);

% plot the close-to minimum points
plot3(x(xi),'markerfacecolor','markersize',12)

现在最大程度地做到了：

% finding maximum
maxval = max(max(fxy));
[xi,yi] = find( maxval == fxy);
% plot the maximum as a red ball
plot3(x(xi),fxy(xi,yi)',12 )

现在如何找到上脊？我将上限设置为5，结果是这样的：

% finding upper ridge : points below a threshold of 5

[xi,yi] = find( maxval- fxy <5);

% to get the values,12)

仅获得上脊的最佳阈值是什么？还有其他方法可以只在上边缘绘制点吗？

解决方法

对于顶部山脊，通常情况下，我会遵循Adriaan的评论并看一下派生词（您可以使用功能gradient）。

但是，在您的情况下，还有另一种方法，它可以通过以下事实实现：

• 在表面X的任何fxy切片上，脊y坐标将与此X切片的最大值重合。
• 并且，在表面Y的任何fxy切片上，脊x坐标将与此Y切片的最大值重合。

基于此观察，max函数是您所需要的：

% Find maximum Z value and index for
[zxmax,idxzx] = max(fxy) ;          % All X slices (columns)
[zymax,idxzy] = max(fxy,[],2) ;     % All Y slices (rows)

xr = x(idxzx) ; % get the actual X coordinates from the column indices
yr = y(idxzy) ; % get the actual Y coordinates from the row indices

% Display
hold on
plot3(xr,y,zxmax,'r','LineWidth',4)
plot3(x,yr,'b',4)

% Or if you want them dotted
% plot3(xr,'or')
% plot3(x,zymax,'ob')

绘制在初始曲面显示上方的图形将呈现：