Нежелательный шумовой элемент во время дилатации

Виновником является ваш оператор if внутри вторых вложенных блоков цикла for:

for(int j = -limit; j <=limit ; j++) {
    for(int i = -limit; i <=limit ; i++) {
        if(*(SE + (j+1)*m + (i+1) ) > 0) {
            pixelvalue = (*(SE + (j+1)*m + ((2*limit)-(i+1)) ))  (Image.at<unsigned char>(y+j,x+i));
            if (pixelvalue>127)
                max=255.0;
        }
    }
}

То, как вы индексируете свой элемент структурирования, предполагает, что это маска 3 x 3, потому что вы ожидаете, что i и j колеблются от -1 до 1. Марк Рэнсом прав в своей оценке, у вас будут индексы с отрицательными значениями. Попробуйте установить m = 15, в этом и заключается ваша проблема. Это означает, что limit = 7, а затем начните с начала этого цикла, посмотрите, к какой строке и столбцу вы обращаетесь, когда пытаетесь получить доступ к маске с линейными индексами по строкам.

Если limit = 7, значит в начале j = -7, i = -7 и так (j+1)*m = (-7+1)*15 = -90. Также (2*limit) - (i+1) = (2*7) - (-7+1) = 20. В конечном итоге это означает, что вы пытаетесь получить доступ к ячейке памяти SE - 90 + 20 = SE - 70, поэтому вы фактически оцениваете ячейку памяти в 70 байтах (при условии, что unsigned char) от источника вашего элемента структурирования, а не в верхнем левом углу элемента структурирования. ... так вы уверены, что написали этот код правильно? Я удивлен, что вы еще не получили никаких ошибок сегментации. Возможно, вам повезло с 9 x 9, потому что смещение строки отрицательное, но смещение столбца будет положительным, и они могут быть равны по величине, так что вам, вероятно, повезло с доступом ко всем элементам структурирования тогда, но это, конечно, не будет быть верным для более высоких размеров маски.

Чтобы обеспечить охват всех коэффициентов маски, необходимо изменить смещение на +limit, а не +1. Я также не понимаю 2*limit внутри оператора if. Вы можете просто получить доступ к элементу маски напрямую. Поэтому нужно сделать так:

for(int j = -limit; j <=limit ; j++) {
    for(int i = -limit; i <=limit ; i++) {
        if(*(SE + (j+limit)*m + (i+limit) ) > 0) {
            pixelvalue = (*(SE + (j+limit)*m + (i+limit) )) * (Image.at<unsigned char>(y+j,x+i));
            if (pixelvalue>127)
                max=255.0;
        }
    }
}

Однако, если я могу что-то предложить, поскольку вы выполняете морфологическое бинарное расширение, вам нужно только проверить, не соприкасаются ли какие-либо пиксели в вашем структурирующем элементе с белым пикселем. Поскольку это квадратный структурирующий элемент, это будет еще проще. По сути, в окрестности m x m пикселя, если какой-либо пиксель не равен нулю, вы должны установить вывод на 255 и выйти из цикла. Небольшая оптимизация вашего кода может быть:

for(int j = -limit; j <=limit ; j++) {
    for(int i = -limit; i <=limit ; i++) {
        pixelValue = (*(SE + (j+limit)*m + (i+limit) )); // Get structuring element value
        pixelValue *= (Image.at<unsigned char>(y+j,x+i)); // Obtain corresponding intensity
        if (pixelValue != 0) { // If non-zero, set output to max and break
            max = 255.0;
            i = limit + 1; j = limit + 1; // break out of loop
        }
    }
}

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