MATLAB图像亮度变换和空间滤波
f=imread('你要处理的图片');
imshow(f);
g=imcomplement(f); %g是翻转后的图像
figure,imshow(g);
g2=imadjust(f,[0.5,0.7],[0.1,0.5],0.5);%突出显示中间灰度部分
figure,imshow(g2,[]);
g3=imadjust(f,[0.5,0.75],[0.1,0.5],2);%gamma=2的整体亮度拉伸，图像总体将变暗
figure,imshow(g3,[]);%对数变换（这个图2感觉好漂亮啊）
f=imread('你的文件名');
imshow(f);
flog=log(1+double(f));%flog是f做对数变换后的图像
flog1=mat2gray(flog);%值变换到[0，1]之间的图像
g=im2uint8(flog1); %最终对数变换后的频谱图像
figure,imshow(g); %对数变换后的图像拓展了动态范围
%亮度拉伸变幻
f=imread('你的文件'); %f
imshow(f);%fmean是图像f的均值，作为亮度拉伸的控制点
fmean=mean2(im2double(f));%利用公式（2-2）对图像进行亮度拉伸
g=intrans(f,'stretch',fmean,0.9);
figure,imshow(g); 
f=imread('你的文件名'); %f是花粉图像
imshow(f); %花粉图像灰度分布不均匀
figure,imhist(f); %显示直方图
ylim('auto'); %直方图纵轴坐标刻度自动调整
%g是直方图均衡化的图像，均匀分布到256级灰度
g=histeq(f,256);
figure,imshow(g);
figure,imhist(g);
ylim('auto');
%介绍几种边缘处理方式
f=[1 1 1 0 0;2 2 2 0 0;3 3 3 0 0;0 0 0 4 4;0 0 0 5 5]; %f是5x5矩阵
imshow(f,[]);
w=ones(3,3); %w=3x3平均滤波器
gd=imfilter(f,w); %默认选项滤波，采用0扩充边界
%显示滤波结果
figure,imshow(gd,[]);
gr=imfilter(f,w,'replicate'); %使用复制边界选项扩充边界的滤波结果
%显示滤波结果
figure,imshow(gr,[]);
gr=imfilter(f,w,'symmetric'); %采用镜像反射方式扩充边界的滤波结果
%显示滤波
figure,imshow(gr,[]);
gc=imfilter(f,w,'circular'); %采用周期方式扩展边界的滤波结果
%显示滤波结果
figure,imshow(gc,[]);
%空间滤波
f=imread('你的图');
imshow(f);
w=[0 1 0;1 -4 1;0 1 0]; %手工指定拉普拉斯滤波器
g1=imfilter(f,w,'replicate'); %g1是空间滤波后的图像，g1的负值被截掉
figure,imshow(g1,[]);
f2=im2double(f); %f被转换为double类型
g2=imfilter(f2,w,'replicate');%同样的滤波操作，但是g2中会保留负值
figure,imshow(g2,[]);
g=f2-g2; %原图像减去空间滤波后的图像
figure,imshow(g2,[]); %处理后的凸显比原图像要清晰
%手工指定空间滤波器
f=imread('moon.tif'); %f同上，是月球表面图像
w4=fspecial('laplacian',0);%w4是拉普拉斯滤波器，中心元素的值是-4
w8=[ 1 1 1;1 -8 1;1 1 1];%w8是自定义滤波器，中心元素的值是-8
f=im2double(f); %f转换到double类型
g4=f-imfilter(f,w4,'replicate'); %g4是原图像减去空间滤波的结果
g8=f-imfilter(f,w8,'replicate'); %g8是原图像减去空间滤波的结果
figure,imshow(f);
figure,imshow(g4);
%可以看到手工处理的结果图像更加清晰
figure,imshow(g8);
%中值滤波
f=imread('你的图片');
imshow(f);%fn时候被椒盐噪声污染了图像
fn= imnoise(f,'salt & pepper',0.2);
figure,imshow(fn);
gm=medfilt2(fn); %gm是经过中值滤波的图像
figure,imshow(gm);
