Matlab 创建单面距离贴图

我用这段代码创建了一个黑色的0平方，中间有一条白色的1垂直线，我想创建一个关于这条线的距离贴图，但只在垂直线的左侧。我该怎么做

下面的代码在直线的两侧生成一个距离贴图

c=zeros(500,500);
c(:,250)=1;
figure, imshow(c)

[D, idx]= bwdist(c,'euclidean')

---

您可以计算整个图像的距离贴图，然后将不感兴趣的一侧归零（或将值设置为

NaN

）

D = bwdist(c, 'euclidean'); 
D(:,251:end) = NaN;

一种更可靠的方法（无需对任何列进行硬编码）是在调用

bwdist

之前修改

c

，方法是将行右侧的任何内容设置为

1

，使每个像素的结果距离为

0

。您可以通过计算跨行的累积和来实现这一点

D = bwdist(cumsum(c, 2) > 0, 'euclidean');

---

事实上，这个问题本可以用另一种富有挑战性的方式提出。如果我有一个不同的输入形状（不是一条单行），如果我想在特定的角度间隔（不仅仅是左侧）中找到距离贴图，该怎么办。这是我的解决办法。找到输入形状的（x，y）坐标，然后为每个像素找到特定方向+方向间隔的遮罩。下面是函数

function mask = directional_mask(im, mainDirection, dirInterval)

% im: BW input image
% mainDirection : mask direction (in degrees)
% dirInterval:  mask direction interval (from dirMask-dirInterval to dirMask+dirInterval)

[Ys, Xs] = find(im);

[sizeX, sizeY] = size(im);

[X, Y] = meshgrid(1:sizeX,1:sizeY);

mask = im * 0;

for i = 1 : numel(Xs)

    refX = Xs(i);
    refY = Ys(i);

    [theta, ~] = cart2pol(X-refX, Y-refY);

    % adding pi/2 so that north angle is 0 degree
    dirmat = wrapTo2Pi(theta + pi/2);
    dirmat = rad2deg(dirmat);

    % upper and lower direction intervals
    thetaUp = mainDirection + dirInterval;
    thetaDown = mainDirection - dirInterval;

    dirmat2 = dirmat(:);

    % finding indices of the angle intervals
    thetas = [thetaUp, thetaDown];
    if thetaUp >= 360
        thetaUp = thetas(1) - 360;
        thetaDown = thetas(2);
        dirmat2((dirmat2>thetaUp) & (dirmat2<thetaDown)) = nan;
    elseif thetaDown < 0
        thetaUp = thetas(1);
        thetaDown = thetas(2) + 360;
        dirmat2((dirmat2>thetaUp) & (dirmat2<thetaDown)) = nan;
    else
        dirmat2((dirmat2>thetaUp) | (dirmat2<thetaDown)) = nan;
    end

    % final mask
    tmp = im*0;
    tmp(~isnan(dirmat2)) = 1;

    mask = mask | tmp;

end

mask = double(mask);

假设我想要180度的距离图。只要改变“主要方向”

让我们试一下不同的形状

im=zeros(100);
im(40:60,40:60)=1; im2=imrotate(im,30);im=imresize(im2,size(im)); im=im~=0;

mainDirection = 180;
dirInterval = 5;

如果你想要一个特定方向间隔（θ+-间隔）的遮罩，只需更改“dirInterval”。这里有一个例子

mainDirection = 180;
dirInterval = 30;

在最初的问题中，@Mac想屏蔽左侧的距离图。然后：

im=zeros(100); im(:,50)=1;

mainDirection = 270;
dirInterval = 90;

mainDirection = 180;
dirInterval = 30;

im=zeros(100); im(:,50)=1;

mainDirection = 270;
dirInterval = 90;