conv_template1.rar_convtemplate_edge sharpening_matlab 锐化_平滑模板_边

%-------------------------------------------------------------------------- %函数名称 % conv_template1 %输入参数 % f 二维数组x 或者为一副图像 % tX 模板卷积时中心的 X 坐标 % tY 模板卷积时中心的 Y 坐标 %输出变量 % g 输出图像 f 和模板template卷积后的卷积结果 %说明 % 图像和模板自带高度和宽度信息 % 注意边界问题 % 图像坐标系 -----------> X % | % | % | % Y % 在matlab中下标均从1开始 %-------------------------------------------------------------------------- function g = conv_template1(f,tX,tY) %f = double(f); [height_f width_f] = size(f); height_t = 3; width_t = 3; %------------------------------给g赋值------------------------------------- %模板卷积达不到的地方用原像素值代替 f = double(f); g = double(f); %gg = g; %---------------------在有效的范围内使用模板卷积--------------------------- for i = tY:1:height_f - height_t + tY for j = tX:1:width_f - width_t + tX temp = 0; ll = 1; for k = 1:height_t for l = 1:width_t %与(k,l)对应的像素在原图的坐标为 % (i-tY+k,j-tX+l) %temp = temp + f(i-tY+k,j-tX+l) * template(k,l); ff(ll) = g(i-tY+k,j-tX+l); ll = ll+1; end end %fr = 4 %转化成细化形式 fa(1) = ff(6); fa(2) = ff(3); fa(3) = ff(2); fa(4) = ff(1); fa(5) = ff(4); fa(6) = ff(7); fa(7) = ff(6); fa(8) = ff(9); fa(9) = fa(1); %细化 a = 0; b = 0; for m = 1:8 a = a + abs(fa(m+1) - fa(m)); b = b + fa(m); end a; A = (fa(1)&fa(7)) & (fa(2)|fa(6)) & ((fa(3)|fa(4)|fa(5)|fa(8))==0); B = (fa(1)&fa(3)) & (fa(4)|fa(8)) & ((fa(2)|fa(5)|fa(6)|fa(7))==0); C = (fa(3)&fa(5)) & (fa(2)|fa(6)) & ((fa(1)|fa(4)|fa(7)|fa(8))==0); D = (fa(5)&fa(7)) & (fa(8)|fa(4)) & ((fa(1)|fa(6)|fa(3)|fa(2))==0); %gggg=9 if ((a==0)|(a==2)) if (b~=1) %aaa = 45 if ((fa(1)&fa(3)&fa(5))==0) if ((fa(1)&fa(3)&fa(7))==0) g(i,j) = 0; %aaa = 3 end end end end if (a==4)&(A|B) if (b~=1) if ((fa(1)&fa(3)&fa(5))==0)&((fa(1)&fa(3)&fa(7))==0) g(i,j) = 0; %aaa = 4 end end end %-------------------------------------------------- if ((a==0)|(a==2)) if (b~=1) %aaa = 45 if ((fa(7)&fa(3)&fa(5))==0) if ((fa(1)&fa(5)&fa(7))==0) g(i,j) = 0; %aaa = 3 end end end end if (a==4)&(C|D) if (b~=1) if ((fa(1)&fa(3)&fa(5))==0)&((fa(1)&fa(3)&fa(7))==0) g(i,j) = 0; %aaa = 4 end end end end end %-------------------------把结果转化为uint8------------------------------- %根据需要。。。 %g = abs(g); %g =uint8(g + 0.5); %imshow(g);