据我了解,您有一组表示椭圆的点,您可以直接在图像矩阵中绘制它们(而不仅仅是在屏幕上显示它们)。
为此,您可以使用POLY2MASK 函数将椭圆转换为二进制掩码。然后通过计算它的perimeter,这将给我们一个二进制掩码,只表示构成椭圆的像素,将其应用于图像以设置像素的颜色。
考虑下面的例子。我正在使用来自上一个问题的函数calculateEllipse.m SO:
%# some image
I = imread('pout.tif');
sz = size(I);
%# ellipse we would like to draw directly on image matrix
[x,y] = calculateEllipse(100,50, 150,50, 30, 100);
%# lets show the image, and plot the ellipse (overlayed).
%# note how ellipse have floating point coordinates,
%# and also have points outside the image boundary
figure, imshow(I)
hold on, plot(x,y, 'LineWidth',2)
axis([-50 250 -50 300]), axis on
%# create mask for image pixels inside the ellipse polygon
BW = poly2mask(x,y,sz(1),sz(2));
%# get the perimter of this mask
BW = bwperim(BW,8);
%# use the mask to index into image
II = I;
II(BW) = 255;
figure, imshow(II)
这应该会给您提供比简单地舍入x 和y 的坐标更好的结果(另外它可以为我们处理越界点)。请务必阅读 POLY2MASK 的算法部分,了解它在亚像素级别上的工作原理。
编辑:
如果您使用的是 RGB 图像(3D 矩阵),同样适用,您只需更改我们使用二进制掩码的最后一部分:
%# color of the ellipse (red)
clr = [255 0 0]; %# assuming UINT8 image data type
%# use the mask to index into image
II = I;
z = false(size(BW));
II( cat(3,BW,z,z) ) = clr(1); %# R channel
II( cat(3,z,BW,z) ) = clr(2); %# G channel
II( cat(3,z,z,BW) ) = clr(3); %# B channel
figure, imshow(II)
这是另一种方式:
%# use the mask to index into image
II = I;
BW_ind = bsxfun(@plus, find(BW), prod(sz(1:2)).*(0:2));
II(BW_ind) = repmat(clr, [size(BW_ind,1) 1]);
figure, imshow(II)