我正在尝试此代码生成基于https://www.crisluengo.net/archives/324中的代码的自由链代码,但它使用DIPimage。因此,有人知道如何通过dip_array函数吗?
代码:
clc;
clear all;
Image = rgb2gray(imread('https://upload-icon.s3.us-east-2.amazonaws.com/uploads/icons/png/1606078271536061993-512.png'));
BW = imbinarize(Image);
BW = imfill(BW,'holes');
BW = bwareaopen(BW, 100);
BW = padarray(BW,60,60,'both')
BW = imcomplement(BW);
imshow(BW)
[B,L] = bwboundaries(BW,'noholes');
%%%%%%%https://www.crisluengo.net/archives/324%%%%%
directions = [ 1, 0
1,-1
0,-1
-1,-1
-1, 0
-1, 1
0, 1
1, 1];
indx = find(dip_array(img),1)-1;
sz = imsize(img);
start = [floor(indx/sz(2)),0];
start(2) = indx-(start(1)*sz(2));
cc = []; % The chain code
coord = start; % Coordinates of the current pixel
dir = 1; % The starting direction
while 1
newcoord = coord + directions(dir+1,:);
if all(newcoord>=0) && all(newcoord<sz) ...
&& img(newcoord(1),newcoord(2))
cc = [cc,dir];
coord = newcoord;
dir = mod(dir+2,8);
else
dir = mod(dir-1,8);
end
if all(coord==start) && dir==1 % back to starting situation
break;
end
end
我不想翻译整个代码,我现在没有时间,但我可以给出一些指针:
-
dip_array(img)
提取具有dip_image
对象img
内部像素值的MATLAB数组。如果您在这里使用BW
作为输入图像,您可以简单地删除对dip_array
:indx = find(BW,1)-1
的调用。 -
imsize(img)
返回图像img
的大小。MATLAB函数size
是等价的(在这种特殊情况下)。 -
dip_image
对象的维度与MATLAB数组的维度不同:它们被索引为img(x,y)
,而MATLAB数组被索引为BW(y,x)
。 -
dip_image
对象的索引从0开始,而不是像MATLAB数组那样从1开始。
最后两点改变了计算start
的方式。我想应该是这样的:
indx = find(BW,1);
sz = size(BW);
start = [1,floor((indx-1)/sz(2))+1];
start(1) = indx-((start(2)-1)*sz(1));
但是使用ind2sub
更容易(不确定为什么我在博客文章中做了显式计算):
indx = find(BW,1);
sz = size(BW);
start = ind2sub(sz,indx);
出于同样的原因,您可能还想交换directions
的两列,并将all(newcoord>=0) && all(newcoord<sz)
更改为all(newcoord>0) && all(newcoord<=sz)
。