Рекурсивная разметка связанных областей

Image = [0 0 0 0 0 0 0 1 0 0 0;
     0 0 0 1 1 0 0 1 0 0 0;
     0 0 0 1 1 0 0 1 1 0 0;
     0 1 1 1 1 0 0 1 1 1 1;
     0 1 1 1 1 0 0 0 0 0 0;
     0 1 1 1 0 0 0 0 0 0 0;
     0 0 0 0 0 0 0 0 0 0 0;
     0 0 0 0 0 0 0 1 1 0 0;
     0 0 0 0 0 0 1 1 1 0 0;
     0 1 1 0 0 1 1 1 1 0 0;
     0 1 1 0 0 1 1 1 1 1 1;
     0 1 1 1 0 1 1 0 1 1 0;
     0 0 0 0 0 0 0 0 0 0 0]
Matplot(Image*255) // Посмотрим на нашу матрицу как на картинке
[m,n]=size(Image); // Узнаем горизонтальный и вертикальные размеры матрицы
km = 0; kn = 0; // Они нам еще пригодятся
cur = 1; // Переменная для подсчета объектов
// Цикл по пикселям изображения
for i = 1:1:m
  for j = 1:1:n
    kn = j - 1;
    if kn <= 0 then
      kn = 1;
      B = 0;
    else
      B = Image(i,kn); // Смотри рисунок 3 в статье
    end
    km = i - 1;
    if km <= 0 then
      km = 1;
      C = 0;
    else
      C = Image(km,j); // Смотри рисунок 3 в статье
    end
    A = Image(i,j); // Смотри рисунок 3 в статье
    if A == 0 then // Если в текущем пикселе пусто - то ничего не делаем
    elseif B == 0 & C == 0 then // Если вокруг нашего пикселя пусто а он не пуст - то это повод подумать о новом объекте
        cur = cur + 1;
        Image(i,j) = cur;
    elseif B ~=0 & C == 0 then
        Image(i,j) = B;
    elseif B == 0 & C ~= 0 then
        Image(i,j) = C;
    elseif B ~= 0 & C ~= 0 then        
        if B == C then
          Image(i,j) = B;
        else
          Image(i,j) = B;
          Image(Image == C) = B;
        end
    end 
  end
end