0

我正在尝试使用此代码根据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
4

1 回答 1

1

我不想翻译整个代码，我现在没有时间，但我可以指点一下：

  • dip_array(img)dip_image提取具有对象内部像素值的 MATLAB 数组img。如果您BW在此处用作输入图像，则可以简单地删除对dip_array:的调用indx = find(BW,1)-1。

  • imsize(img)返回图像的大小img。MATLAB 函数size是等效的（在这种特殊情况下）。

  • 对象的维度与dip_imageMATLAB 数组的维度不同：它们的索引为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).

于 2021-01-16T22:15:58.880 回答