matlab图像处理效果展示22课件

1、图像处理效果展示2常看看探索RGB图像的像素特点原图程序 A=imread(山水3.jpg); B=A; B(:,:,1)=0; figure image(A) figure image(B) C=A; C(:,:,2)=0; figure image(C) D=A; D(:,:,3)=0; figure image(D)探索RGB图像的像素特点取出红色探索RGB图像的像素特点取出绿色探索RGB图像的像素特点取出蓝色对比两种取色方法A=imread(山水3.jpg);C=A;C(:,:,2)=0;G=imsubtract(A,C);Figuresubplot(121);imshow(G)sub

2、plot(122);imshow(A(:,:,2)对比两种取色方法A=imread(山水3.jpg);D=A;D(:,:,3)=0;B=imsubtract(A,D);Figuresubplot(121);imshow(B)subplot(122);imshow(A(:,:,3)图像乘法运算 A=imread(海豚1.jpg)； C=immultiply(A,0.5);图像乘（除）法运算 A=imread(海豚1.jpg); D=imdivide(A,2);图像增强之中值滤波 clear A=imread(山水.jpg); B=rgb2gray(A); a,b=size(B); C=B; fo

3、r i=2:a-1for j=2:b-1G=B(i-1,j) B(i+1,j) B(i,j) B(i,j-1) B(i,j+1);F=sort(G);C(i,j)=F(3);endend figure imshow(C)前后图像增强之邻域平均 clear A=imread(山水.jpg); B=rgb2gray(A); a,b=size(B); K=B; for i=2:a-1for j=2:b-1R=B(i-1,j) B(i+1,j) B(i,j-1) B(i,j+1);K(i,j)=ceil(sum(R)/4);endend figure imshow(K)前前前后 图像增强之大值滤波 A

4、=imread(刘翔.jpg); B=rgb2gray(A); a,b=size(B); C=B; for i=2:a-1for j=2:b-1if B(i,j)200C(i,j)=B(i,j);elseC(i,j)=0;endendend figure subplot(121);imshow(B) title() title(处理前); subplot(122);imshow(C) title(大值滤波器滤波后);图像增强之锐化滤波 A=imread(傅.jpg); a,b=size(A); B=A; for i=2:a-1for j=2:b-1x1=abs(A(i,j)-A(i,j-1);

5、x2=abs(A(i,j)-A(i,j+1);y1=abs(A(i,j)-A(i-1,j);y2=abs(A(i,j)-A(i+1,j);d1=abs(A(i,j)-A(i-1,j-1);d2=abs(A(i,j)-A(i+1,j-1);d3=abs(A(i,j)-A(i-1,j+1);d4=abs(A(i,j)-A(i+1,j+1);H=x1 x2 y1 y2 d1 d2 d3 d4;t=max(H);if t0 B(i,j)=t;else B(i,j)=0;endendend figure imshow(A) figure imshow(B)图像增强之锐化滤波x1=abs(B(i,j)-B

6、(i,j-1);x2=abs(B(i,j)-B(i,j+1);y1=abs(B(i,j)-B(i+1,j);y2=abs(B(i,j)-B(i-1,j);图像增强之锐化滤波d1=abs(B(i,j)-B(i-1,j-1);d2=abs(B(i,j)-B(i-1,j+1);d3=abs(B(i,j)-B(i+1,j-1);d4=abs(B(i,j)-B(i+1,j+1);考虑去噪过程中减少边缘模糊化原灰度图去噪不顾虑边缘模糊加权平均去噪减少边缘模糊程序实现（续前图） A=imread(山水2.jpg); B=rgb2gray(A); C=B; a,b=size(B); for i=2:a-1fo

7、r j=2:b-1T=B(i,j) B(i-1,j) B(i+1,j) B(i,j+1) B(i,j-1);t=1/5*sum(T(:);C(i,j)=t;endend figure imshow(C) figure imshow(B) D=B; for i=2:a-1for j=2:b-1T=B(i,j) B(i-1,j) B(i+1,j) B(i,j+1) B(i,j-1);t=1/5*sum(T(:);if B(i,j)uint8(t+20)D(i,j)=t;elseD(i,j)=B(i,j);endendend figure imshow(D)傅里叶变换展示图傅里叶变换 clear c

8、lose allA=imread(fft1.jpg); B=double(A); C=fft2(B); D=fftshift(C); figure subplot(121);imshow(A) subplot(122);imshow(log(1+abs(D),)滤去高频部分显示效果 A=imread(山水2.jpg); B=rgb2gray(A); C=double(B); D=fft2(C); E=fftshift(D); figure subplot(121);imshow(A)subplot(122);imshow(log(1+abs(E),) D(abs(D) figure subpl

9、ot(121);imshow(log(1+abs(D),) F=fftshift(D); subplot(122);imshow(log(1+abs(F),)去掉高频部分后回复灰度图像效果展示ifftshift指令；ifft2指令；uint8指令的运用； 冈萨雷斯版里面的解释非常形象：一个恰当的比喻是将傅里叶变换比作一个玻璃棱镜。棱镜是可以将光分解为不同颜色的物理仪器，每个成分的颜色由波长（或频率）来决定。 傅里叶变换可以看作是数学上的棱镜，将函数基于频率分解为不同的成分。当我们考虑光时,讨论它的光谱或频率谱。同样, 傅立叶变换使我们能通过频率成分来分析一个函数。 图像傅立叶变换的物理意义 图

10、像的频率是表征图像中灰度变化剧烈程度的指标，是灰度在平面空间上的梯度。如：大面积的沙漠在图像中是一片灰度变化缓慢的区域，对应的频率值很低；而对于地表属性变换剧烈的边缘区域在图像中是一片灰度变化剧烈的区域，对应的频率值较高。傅立叶变换在实际中有非常明显的物理意义，设f是一个能量有限的模拟信号，则其傅立叶变换就表示f的谱。从纯粹的数学意义上看，傅立叶变换是将一个函数转换为一系列周期函数来处理的。从物理效果看，傅立叶变换是将图像从空间域转换到频率域，其逆变换是将图像从频率域转换到空间域。换句话说，傅立叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数，傅立叶逆变换是将图像的频率分布函数变

11、换为灰度分布函数 傅立叶变换以前，图像（未压缩的位图）是由对在连续空间（现实空间）上的采样得到一系列点的集合，我们习惯用一个二维矩阵表示空间上各点，则图像可由z=f(x,y)来表示。由于空间是三维的，图像是二维的，因此空间中物体在另一个维度上的关系就由梯度来表示，这样我们可以通过观察图像得知物体在三维空间中的对应关系。为什么要提梯度？因为实际上对图像进行二维傅立叶变换得到频谱图，就是图像梯度的分布图，当然频谱图上的各点与图像上各点并不存在一一对应的关系，即使在不移频的情况下也是没有。我现在是将一幅图做傅里叶变换，取其幅值和频率分别作图，我想问下，这两幅图代表了图像的什么东西？是每个像素之间的差

12、还是什么？ 傅立叶频谱图上我们看到的明暗不一的亮点，实际上图像上某一点与邻域点差异的强弱，即梯度的大小，也即该点的频率的大小（可以这么理解，图像中的低频部分指低梯度的点，高频部分相反）。一般来讲，梯度大则该点的亮度强，否则该点亮度弱。这样通过观察傅立叶变换后的频谱图，也叫功率图，我们首先就可以看出，图像的能量分布，如果频谱图中暗的点数更多，那么实际图像是比较柔和的（因为各点与邻域差异都不大，梯度相对较小），反之，如果频谱图中亮的点数多，那么实际图像一定是尖锐的，边界分明且边界两边像素差异较大的。对频谱移频到原点以后，可以看出图像的频率分布是以原点为圆心，对称分布的。将频谱移频到圆心除了可以清晰

13、地看出图像频率分布以外，还有一个好处，它可以分离出有周期性规律的干扰信号，比如正弦干扰，一副带有正弦干扰，移频到原点的频谱图上可以看出除了中心以外还存在以某一点为中心，对称分布的亮点集合，这个集合就是干扰噪音产生的，这时可以很直观的通过在该位置放置带阻滤波器消除干扰 另外我还想说明以下几点： 1、图像经过二维傅立叶变换后，其变换系数矩阵表明： 若变换矩阵Fn原点设在中心，其频谱能量集中分布在变换系数短阵的中心附近(图中阴影区)。若所用的二维傅立叶变换矩阵Fn的原点设在左上角，那么图像信号能量将集中在系数矩阵的四个角上。这是由二维傅立叶变换本身性质决定的。同时也表明一股图像能量集中低频区域。 2

14、 、变换之后的图像在原点平移之前四角是低频，最亮，平移之后中间部分是低频，最亮，亮度大说明低频的能量大（幅角比较大）第二个问题是我想绘制图像的功率谱，和白噪声做对比，是画FFT后的幅值的平方是吧？还有图像的功率代表了什么意义呢? 彩色模型RGB：M*N*3R=M*N*1=（：，：，1）G=M*N*1=（：，：，2）B=M*N*1=（：，：，3）彩色模型CMY模型：大多数在纸上沉积的设备都需这种格式，印刷模式。C青色，M洋红色，Y黄色。我们看物体时使用了减色模式，这种格式也使用了减色模式。CMY三色混合生成纯黑色，但是染料中杂质的存在，是混合颜色为灰色，我们需加入黑色颜料调色，所以有成为CMYK

15、模式对称问题 A=imread(山水3.jpg); B=A; a,b,c=size(A); a1=floor(a/2); b1=floor(b/2); c1=floor(c/2); B(1:a1,1:b,1:c)=A(a:-1:a-a1+1,1:b,1:c); figure subplot(121);imshow(A) subplot(122);imshow(B)上下对称图相对称 C=A; C(1:a,1:b1,1:c)=A(1:a,b:-1:b+1-b1,1:c); figure subplot(121);imshow(A) subplot(122);imshow(C)左右对称图像对称 D=

16、A; if a figure subplot(121);imshow(A) subplot(122);imshow(D)斜线对称图像镜像(jingxaing.m)123456789321654987789456123直方图均衡化直方图均衡化的效果图像直方图规定化直方图规定化的效果 A=imread(海豚1.jpg); B=rgb2gray(A); A1=imread(刘翔.jpg); B1=rgb2gray(A1); figure subplot(221);imshow(B) subplot(222);imshow(B1) subplot(223);imhist(B) subplot(224)

17、;imhist(B1) figure subplot(221); C=histeq(B,imhist(B1); imhist(C) title(按刘翔.jpg规定直方图转化的海豚直方图); subplot(222);imhist(B) title(海豚的直方图); subplot(223); C1=histeq(B1,imhist(B); imhist(C1) title(按海豚.jpg规定直方图转化的刘翔直方图); subplot(224);imhist(B1) 通过删去行和烈将高级的图像变换为低级的图像。牡丹图像不同的大小展示牡丹图像不同的大小展示牡丹图像不同的大小展示牡丹图像不同的大小展

18、示牡丹图像不同的大小展示将图像恢复原来大小将原图像隔行隔列删减后，恢复为原来大小时图像的效果（复制删减后存在的像素点的灰度值）将图像恢复原来大小512将图像恢复原来大小将图像恢复原来大小灰度级变化256-8clear A=imread(牡丹3.jpg);B=rgb2gray(A);a,b=size(B);for i=1:afor j=1:bif rem(B(i,j),32)=16D(i,j)=B(i,j)+(32-rem(B(i,j),32);elseif rem(B(i,j),32)0&rem(B(i,j),32) H=uint8(D); figure subplot(121);imshow(B) subplot(122);imshow(H)有问题的转换 A=imread(伞.jpg);B=rgb2gray(A); a,b=size(B);for i=1:afor j=1:bif B(i,j)125C(i,j)=255;elseC(i,j)=0;endendend D=uint8(C); figure imshow(D) E=D; E=D; for i=2:a-1for j=2:b-1x1=abs(B(i,j)-B(i,j-1);x2=abs(B(i,j)-B(i,j+1);y1=abs(B(i,j)-B(i-1,j);y2=abs(B(i,j)-B(i+1,j);H=x1 x