冈萨雷斯 数字图像处理(第三版)matlab代码 图3.43 混合空间增强

clear all; clc ; %%%%混合空间增强 I= imread('Fig0343(a)(skeleton_orig).bmp'); I=rgb2gray(I); %将彩色图像转化为灰度图像 %%%拉普拉斯算子 I1=fspecial('laplacian',0); %掩膜采用标准laplacian掩膜fspecial('laplacian',0) I2=imfilter(I,I1,'replicate'); %滤波 %%%%%%效果并不理想，因为原图是一幅uint8类图像，输出结果仍为uint8类，所有像素均为正值，而拉普拉斯滤波模板中存在负值，变换结果中的所有负值被截掉了。 %%%%%%解决此问题，须先进行类型变换 J2=im2double(I); %将f转换为归一化的double类图像，然后再进行滤波： J3=imfilter(J2,I1,'replicate'); %滤波 % g1=I-I2; g2=J2+J3; %拉普拉斯算子 %%%sobel算子 K1=fspecial('sobel'); K2=filter2(K1,I); %%%均值滤波 M1=filter2(fspecial('average',5),K2)/255; %%%相乘 N1=immultiply(M1,g2); %%%求和 O1=imadd(double(I),N1); %%%%幂律变换(伽马变换) P=uint8(O1); P1=imadjust(P,[],[],0.5); subplot(241);imshow(I);title('图(a),原始灰度图像'); subplot(242);imshow(J3,[]);title('图(b),拉普拉斯操作后的结果'); subplot(243);imshow(g2,[]);title('图(c),图a与图b相加得到的锐化图像'); subplot(244);imshow(K2,[0,90]);title('图(d),sobel算子'); subplot(245);imshow(M1);title('图(e),用5×5模板均值滤波器平滑后的sobel图像'); subplot(246);imshow(N1);title('图(f),图(c)和图(e)相乘得到的掩蔽图像'); subplot(247);imshow(O1,[0,255]);title('图(g),图(a)和图(f)求和图像'); subplot(248);imshow(P1);title('图(h),对图(g)幂律变换(伽马变换)图像');