《数字图像处理》实验课上老师让尽量给出自己的源码,查阅资料初步了解了卷积的基本原理,所以针对实验课上的例子给出了自己的卷积变换代码。总之,卷积变换代码并不复杂,前提是了解了卷积的基本原理,这里有一个参考链接感觉讲的很棒,请查阅过后再阅读我的代码。
博客链接:https://blog.youkuaiyun.com/jinv5/article/details/52874880
%%拉普拉斯算子增强(卷积在下面)
clear all
load imdemos circuit%在图像处理工具箱中加载'circuit'图像
I = circuit;
I = double(I);
imshow(I, []);
h = [0 1 0;1 -4 1; 0 1 0];
J = conv2(I, h, 'same');%以h为卷积核对I进行卷积,返回与I相同大小的中心部分
K = I - J;
figure, imshow(J, []);
figure, imshow(K, []);
%%自己编码实现卷积运算
%先补0
myF = zeros(2, 4 + size(I, 2));
myF(3:(size(I, 1) + 2), 1:2) = 0;
myF(3:(size(I, 1) + 2),3:(size(I, 1) + 2)) = I;
myF(3:(size(I, 1) + 2),size(I, 1) + 3:size(I, 1) + 4) = 0;
myF(size(I, 1) + 3:size(I, 1) + 4, 1:size(I, 1) + 4) = 0;
%进行卷积 卷积顺序改变不会影响卷积的结果
P = ones(4 + size(I, 1), 4 + size(I, 2));
for i = 3:(2 + size(I, 1))
for j = 3:(2 + size(I, 2))
for m = 1:3
for n = 1:3
P(i, j) = P(i, j) + myF(i + m-2, j + n-2)*h(m, n);
end
end
end
end
P = P(3:size(I, 1) + 2, 3:size(I, 2));
figure, imshow(P, []);
效果展示:
原图I: 利用matlab卷积函数conv2(I, h, 'same')求得的结果: 自己编码求得的卷积变换结果: