N-D фильтрация многомерных изображений
Прочтите цветное изображение в рабочую область и отобразите его.
originalRGB = imread('peppers.png');
imshow(originalRGB)
Создайте фильтр размытия движения с помощью fspecial функция.
h = fspecial('motion', 50, 45);Примените фильтр к оригинальному изображению, чтобы создать изображение с размытием движения. Обратите внимание, что imfilter является более эффективной памятью, чем некоторые другие функции фильтрации, поскольку она выводит массив совпадающего типа данных, что и входной массив изображений. В этом примере выход является массивом uint8.
filteredRGB = imfilter(originalRGB, h); figure, imshow(filteredRGB)
Снова отфильтровайте изображение, на этот раз задав опцию репликации границы.
boundaryReplicateRGB = imfilter(originalRGB, h, 'replicate');
figure, imshow(boundaryReplicateRGB)
По умолчанию imfilter использует корреляцию, потому что функции создания фильтра тулбокса производят ядра корреляции. Используйте необязательный параметр для свертки.
Создайте выборку.
A = magic(5)
A = 5×5
17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9
Создайте фильтр.
h = [-1 0 1];
Фильтр с помощью корреляции, по умолчанию.
imfilter(A,h)
ans = 5×5
24 -16 -16 14 -8
5 -16 9 9 -14
6 9 14 9 -20
12 9 9 -16 -21
18 14 -16 -16 -2
Фильтруйте с помощью свертки, задавая imfilter с необязательным параметром.
imfilter(A,h,'conv')ans = 5×5
-24 16 16 -14 8
-5 16 -9 -9 14
-6 -9 -14 -9 20
-12 -9 -9 16 21
-18 -14 16 16 2
В этом примере выход imfilter имеет отрицательные значения, когда вход имеет класс double. Чтобы избежать отрицательных значений, преобразуйте изображение в различные типы данных перед вызовом imfilter. Для примера, когда тип входа uint8, imfilter обрезает выходные значения до 0. Также может быть целесообразно преобразовать изображение в целый тип со знаком.
A = magic(5)
A = 5×5
17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9
Фильтрация изображения с помощью imfilter.
h = [-1 0 1]; imfilter(A,h)
ans = 5×5
24 -16 -16 14 -8
5 -16 9 9 -14
6 9 14 9 -20
12 9 9 -16 -21
18 14 -16 -16 -2
Заметьте, что результат имеет отрицательные значения. Чтобы избежать отрицательных значений в выходном изображении, преобразуйте вход изображение в uint8 перед выполнением фильтрации. С момента входа в imfilter является классом uint8, выходы также имеют класс uint8, и imfilter обрезает отрицательные значения до 0.
A = uint8(magic(5)); imfilter(A,h)
ans = 5x5 uint8 matrix
24 0 0 14 0
5 0 9 9 0
6 9 14 9 0
12 9 9 0 0
18 14 0 0 0
Эта функция может использовать преимущества аппаратной оптимизации для типов данных uint8, uint16, int16, single, и double чтобы запустить быстрее.
The imfilter функция вычисляет значение каждого выходного пикселя, используя арифметику с двойной точностью и плавающей точностью. Если результат превышает область значений типов данных, то imfilter обрезает результат до допустимой области значений типа данных. Если это целочисленный тип данных, то imfilter округляет дробные значения.
Если вы задаете равномерное ядро h, тогда центр ядра floor((size(h) + 1)/2).
Для примера - центр фильтра с 4 элементами [0.25 0.75 -0.75 -0.25] является вторым элементом, 0.75. Этот фильтр дает те же результаты, что и фильтрация с 5-элементным фильтром [0 0.25 0.75 -0.75 -0.25].