Удаление выделенных пикселей из небольшого бинарного изображения

Существует множество алгоритмов заливки, которые идентифицируют все связанные пиксели с заданной начальной точкой.

В качестве альтернативы распространенный способ удалить небольшие выбросы, подобные этим, которые возникают из-за шума, - это размыть изображение, а затем расширить его, чтобы вернуться к тому же размеру, хотя, если вы делаете чисто CoG, вам не обязательно нужен шаг расширения.

Как насчет псевдокода:

for( y = 0; y < rows; y++ )
{    
   for ( x = 0; x < cols; x++ )
   {
       if ( pixel( x, y ).isColor() )
       {
          int sum = 0;
          // forgetting about edge cases for clarity...
          if ( !pixel( x-1, y-1 ).isColor() ) sum++;
          if ( !pixel( x,   y-1 ).isColor() ) sum++;
          if ( !pixel( x+1, y-1 ).isColor() ) sum++;
          if ( !pixel( x-1, y   ).isColor() ) sum++;
          if ( !pixel( x+1, y   ).isColor() ) sum++;
          if ( !pixel( x-1, y+1 ).isColor() ) sum++;
          if ( !pixel( x,   y+1 ).isColor() ) sum++;
          if ( !pixel( x+1, y+1 ).isColor() ) sum++;
          if ( sum >= 7 )
          {
             pixel( x, y ).setBackground();
             x -= 1;
             y -= 1;
          }
       }
   }
}

То есть удалить любой пиксель, окруженный 7 фоновыми пикселями. Если вы измените цвет пикселя обратно на самый ранний пиксель, который теперь может быть затронут.

Ваша мера «выброса» может измениться, например. вы можете считать диагональные пиксели равными 1/2. Например. пиксели непосредственно выше, ниже, слева и справа считаются за 2. Затем используйте другое число для суммы.

Вы можете повысить точность, увеличив размер фильтра — например, до 5x5 вместо 3x3. В этом случае пиксели на расстоянии 2 должны иметь еще меньшее значение.

Я считаю, что ваш алгоритм неверен. m10 (извините, я не мог понять, как делать индексы) в двоичном изображении — это сумма x-координат пикселей переднего плана, поэтому ваш код должен быть примерно таким:

for (i = 0, i < rows, i++) {
    for (j = 0, j < cols, j++) {
        if (colorAt(i,j).isForeground()) {
            tempSumX += i;
            totalForeground++;
        }
    }
}
x_ = tempSumX/totalForeground; //where x_ represents the x coordinate of the weighted center of mass.

Предполагая, что это связано с вашим предыдущим сообщением Алгоритм поиска самого длинного участка значения под любым углом в 2D-матрице, вы должны вычислить другие моменты первого и второго порядка в том же цикле:

m01 += j;
m20 += i*i;
m02 += j*j;
m11 += i*j;

(tempSumX в вашем алгоритме — это просто m10, а totalForeground — m00)

Попробуйте, и если у вас все еще возникают проблемы с выбросами, вы можете использовать маркировку подключенных компонентов http://en.wikipedia.org/wiki/Connected_component_labeling, чтобы найти наибольшую массу. (Доступно в Matlab с использованием bwlabel или bwconncomp.)