这里是一个file exchange link到bresenham.m
改变了你的代码来获取所有的4个重心
%// read your input image
im = imread('http://i.stack.imgur.com/xeqe8.jpg');
BW = im>220;
CC = bwconncomp(BW);
stat = regionprops(CC,'Centroid');
figure; imshow(BW); hold on
for x = 1: numel(stat)
plot(stat(x).Centroid(1),stat(x).Centroid(2),'ro');
end
这里是输出:
进一步实施:
%// putting all the Centroid coordinates into corresponding x,y variable
x = [stat(1).Centroid(1),stat(2).Centroid(1),stat(3).Centroid(1),stat(4).Centroid(1)];
y = [stat(1).Centroid(2),stat(2).Centroid(2),stat(3).Centroid(2),stat(4).Centroid(2)];
%// obtain row and col dim
[r,c] = size(BW);
%// get all x,y values connecting the centroid points
[xAll{1},yAll{1}] = bresenham(x(1),y(1),x(4),y(4));
[xAll{2},yAll{2}] = bresenham(x(2),y(2),x(3),y(3));
[xAll{3},yAll{3}] = bresenham(x(3),y(3),x(4),y(4));
%// change row and col subs to linear index
for ii = 1:3
idx{ii} = sub2ind(size(BW),yAll{ii},xAll{ii});
end
%// change grayscale image to 3D (as you want red line)
out = repmat(im,[1,1,3]);
%// obtaining corresponding index of all 3 slices
for ii = 1:3
idxall{ii} = bsxfun(@plus, idx{ii},[0:2].*(r*c));
end
%// keep only the index of 1st slice to 255 and changing rest to 0 to obtain a red line.
%// Similar process for blue line except keep the index in the 3rd slice to 255
out(cat(1,idxall{:})) = 0;
out(idx{1}) = 255;
out(idx{2}) = 255;
out(idx{3}+2*(r*c)) = 255;
%// see what you have obtained
figure; imshow(out);hold on
for x = 1: numel(stat)
plot(stat(x).Centroid(1),stat(x).Centroid(2),'bo');
end
结果:
注:行可能看起来带点由于图片的尺寸较大,但其持续
最后一个数字放大看实线:
进一步说:
您可能必须采取@Spektre的意见,发现倾斜使用atan2
角度。有关更多说明,请参阅his answer。
你的任务不清楚。 “测量这个区域之间的角度”是什么意思? –
图片中您也在广场内圈了圈,但是当我运行您的代码时,我没有看到那个圈子。你是怎么得到这个的? –
这个问题和你[上一个问题](http://stackoverflow.com/q/30462831/2545927)有什么区别? – kkuilla