Блок схема по коду

clc, clear
close all
addpath D:\
 I = imread('12345.bmp');
img=im2double(rgb2gray(I));
 %I = gpuArray(imread('kurs.jpg'));
%imshow(I);
[Gx, Gy] = imgradientxy(img);
[Gmag, Gdir] = imgradient(Gx, Gy);
% figure, imshow(Gmag, []), title('Gradient magnitude')
% figure, imshow(Gdir, []), title('Gradient direction')
% figure, imshow(Gx, []), title('Directional gradient: X axis')
% figure, imshow(Gy, []), title('Directional gradient: Y axis')
 
%операторы собела
 
h=fspecial('sobel');
h2 = h';
gx=imfilter(double(I), [1 -1]);
gy=imfilter(double(I), [ 1;-1]);
g = abs(gx)- abs(gy);
%figure, imshow(uint8(g)), title('sobel');
 
% Добавление шума и очистка изображения
g = imnoise(I,'salt & pepper',0.01);
 
%Перевод в двоичный код
level = graythresh(g);
BW = im2bw(g,0.7); 
 
%Удаляем шум
L = medfilt2(BW,[5 5]);
%figure, imshow(L);
 
%сглаживаем высокочастотный шум
H = fspecial('average');
average = imfilter(L,H,'replicate');
%figure, imshow(L);
 
se=strel('disk',9);
Ldil = imclose(average, se);
 
 
se=strel('square',5);%параметры эрозии 
Ldil=imerode(Ldil,se);%делаем эрозию изображения
[Ldil, num]=bwlabel(Ldil,4);
feats=regionprops (Ldil,'Area','BoundingBox');
big=feats(1);%самый крупный объект
big_idx = 1;
for i=1:1:num;
    if((feats(i).Area)>big.Area)%поиск более крупного из всех
        big=feats(i);
        big_idx = i;
    end
end 
[x,y]=size(Ldil)
for i=1:1:x
   for j=1:1:y
      if (Ldil(i,j)~= 1)
          Ldil(i,j)=0;
      end
   end
end
 
bbox = feats(big_idx).BoundingBox
figure, imshow(Ldil);
imshow(img);
rectangle('Position', bbox, 'EdgeColor', 'green');
line=imcrop(img, [bbox(1),bbox(2)+bbox(4)/2,bbox(3),1]);
rectangle('Position', [bbox(1),bbox(2)+bbox(4)/2,bbox(3),1], 'EdgeColor', 'red');
for i=1:1:bbox(3)
   if(line(1,i)<0.6)
       bin(i)=0;
   else
       bin(i)=1;
   end
end
 
 
bin2 = []
for i=1:100
    bin2 = vertcat(bin2, bin);
end
figure
colormap gray
imagesc(bin2)