如何在Matlab中永久地在图像的两个坐标之间画一条线？ [复制]

Question

我有一组要按顺序连接的点。假设点是 (A1,A2,A3,...A9)；我想将A1连接到A2，A2连接到A3等等，最后将A9连接到A1。

我只需要知道一个可以帮助我将 A1 连接到 A2 的函数，其余的可以使用 for 循环来完成。

我知道连接两点之前已经在这里问过好几次了，但我找不到我需要的答案。一些解决方案建议使用“绘图”和“线条”，但这些函数会将结果叠加在图像上，实际上不会对原始图像进行任何更改。

我确实尝试过，并设法使用“saveas”和“print”函数保存了结果图，但是图像没有以正确的格式保存，并且使用这些函数的参数存在很多问题.此外，我真的不想保存图像，如果我能用线条获得所需的图像，这只是我愿意添加的不必要的开销。 我也尝试过“imline”来画线，但它似乎是交互式的。

This 特定问题完美地反映了我的问题，但是当我运行作为解决方案给出的代码 sn-ps 时，它们都在结果图像中给出了一组点。

我在this 找到here 的图片链接中尝试了上述代码。

虚线是上面链接中所有三个代码 sn-ps 的输出。 例如，我运行第一个代码是这样的：

I = imread('berries_copy.png'); 
grayImage=rgb2gray(I); 
img =false(size(grayImage,1), size(grayImage,2)); 

我写上面这段代码只是为了得到一个黑色的图像用于以下操作：

x = [500 470];        % x coordinates 
y = [300 310];        % y coordinates
nPoints = max(abs(diff(x)), abs(diff(y)))+1;    % Number of points in line
rIndex = round(linspace(y(1), y(2), nPoints));  % Row indices
cIndex = round(linspace(x(1), x(2), nPoints));  % Column indices
index = sub2ind(size(img), rIndex, cIndex);     % Linear indices
img(index) = 255;  % Set the line points to white
imshow(img);       % Display the image

This 是上述代码以及其他两个代码的结果图像，如您所见，它只是黑色背景上的几个点，这不是所需的输出。

我更改了代码并使用“绘图”函数来获得this 输出，这正是我想要的。无论如何我可以将虚线输出更改为实线吗？

或者，如果有人能提出一个简单的函数或方法来从 A1 到 A2 画一条线，并且实际上会改变输入图像，我将不胜感激。 （我真的希望这只是我是一个新手，而不是 Matlab 没有一个简单的功能来在图像中画一条线。）

附：我没有计算机视觉工具箱，如果可能的话，我想找到一个不涉及它的解决方案。

Answer 1

您的第一个问题是您正在创建一个 blank 图像，其大小与您使用此行加载的图像大小相同：

img =false(size(grayImage,1), size(grayImage,2));

当您向其中添加线条时，您会得到一个带有白色线条的黑色图像，正如预期的那样。

您的第二个问题是您尝试将grayscale intensity images 的解决方案应用于RGB (Truecolor) image，这需要您使用modify the data at the given indices for all three color planes（红色、绿色和蓝色）。以下是如何从my other answer修改灰度解决方案：

img = imread('berries_copy.png');  % Load image
[R, C, D] = size(img);             % Get dimension sizes, D should be 3
x = [500 470];                     % x coordinates
y = [300 310];                     % y coordinates
nPoints = max(abs(diff(x)), abs(diff(y)))+1;    % Number of points in line
rIndex = round(linspace(y(1), y(2), nPoints));  % Row indices
cIndex = round(linspace(x(1), x(2), nPoints));  % Column indices
index = sub2ind([R C], rIndex, cIndex);         % Linear indices
img(index) = 255;                  % Modify red plane
img(index+R*C) = 255;              % Modify green plane
img(index+2*R*C) = 255;            % Modify blue plane
imshow(img);                       % Display image
imwrite(img, 'berries_line.png');  % Save image, if desired

以及生成的图像（注意右下角浆果上方的白线）：

Answer 2

您可以使用 Bresenham 算法。当然它已经实现了，你可以在这里找到它：Bresenham optimized for Matlab。该算法选择近似于一条线的像素。

一个简单的例子，使用你的变量名可以是：

I = rgb2gray(imread('peppers.png'));
A1 = [1 1];
A2 = [40 40];
[x y] = bresenham(A1(1),A1(2),A2(1),A2(2));
ind = sub2ind(size(I),x,y);
I(ind) = 255;
imshow(I) 

Answer 3

您可以使用 imshow 显示图像，然后使用 plot 绘制线条并保存图形。检查以下代码：

I = imread('peppers.png') ;
imshow(I)
hold on
[m,n,p] = size(I) ;
%% Get random points A1, A2,..A10
N = 9 ;
x = (n-1)*rand(1,N)+1 ;
y = (m-1)*rand(1,N)+1 ;
P = [x; y]; % coordinates / points 
c = mean(P,2); % mean/ central point 
d = P-c ; % vectors connecting the central point and the given points 
th = atan2(d(2,:),d(1,:)); % angle above x axis
[th, idx] = sort(th);   % sorting the angles 
P = P(:,idx); % sorting the given points
P = [P P(:,1)]; % add the first at the end to close the polygon 
plot( P(1,:), P(2,:), 'k');
saveas(gcf,'image.png')