资源描述：
《数字图象 处理系列问题研究 matlab实现》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、数字图象处理系列问题研究matlab实现一、实验目的:掌握基本的图象增强和变换方法，观察图象增强的效果，加深对图像增强和一些基本变换的理解。熟悉如何用matlab实验平台对图像进行增强和变换处理。二、实验内容:1.编程实现打开一个灰度图像，能够显示图像。2.图像数字化以及相关操作(图像灰度的平均值，协方差矩阵，灰度标准差)。3.编程实现图像傅立叶变换和反变换，实现带噪小丑图像去除噪音。4.用直接灰度变换改变图像(求反，增强对比度，动态范围压缩，灰度切分)。三、实验过程:注:除小丑去噪外其他操作均用图1所示的Gem图像(
2、Gem是我的网名);这里面的所有图象我都先变换为灰度图判断是否为灰度图象:flag=isgray(A)判断是否为真彩图象:flag=isrgb(A)真彩图象→灰度图象:I=rgb2gray(A);1.打开并显示一幅灰度图像:在matlab命令窗口中输入如下代码:I=imread('d:gem.bmp');flag=isgray(I)%判断是否是灰度图片flag=1，ture是灰度图片;否则不是flag=1imshow(I)所显示结果如下:图1原始gem图像2.求图像灰度平均值，协方差矩阵和灰度标准差:(1)求灰度平均
3、值①实现方式1:因为实验过程1中打开灰度图像时已完成图像的数字化，即图像已经变成一个矩阵，因此直接调用matlab自带的求矩阵均值函数mean2()求图像的均值。代码如下:I=imread('d:gem.bmp');J=mean2(I)J=189.8383②实现方式2:先调用size函数求出矩阵的大小，然后求出矩阵所有值的平均值即是原始图像的灰度均值。代码如下:I=imread('d:gem.bmp');%读入原始图像[m,n]=size(I);%求对应矩阵的规格，即图像的像素总数J=double(I);a=0;%
4、变量a赋初值forj=1:m%循环求和fork=1:na=a+J(j,k);endendavg=a/(238*107)%求平均运行结果:avg=189.8383(2)灰度标准差:K=std2(I)K=88.4966(3)协方差矩阵cov(double(I))4.06074.24914.51134.82655.23136.14618.77918.64928.40313.92414.08984.33274.62675.00695.93018.64928.83908.69433.73213.89704.13784.43054
5、.80655.71938.40318.69438.83023.62353.74393.96344.23354.58075.47018.07348.40978.64833.46123.54493.76124.02904.36945.24207.78798.13238.39753.28473.33753.54813.81184.13974.98207.42377.76028.01923.15293.28373.51003.79204.13664.97237.39157.67527.90272.52512.65932.913
6、03.21283.60994.48037.01847.23457.4691..3.图像傅立叶变换和反变换直接调用matlab自带的函数实现，代码如下:I=imread('d:gem.bmp');C=fft2(double(I));%对图像进行傅立叶变换B=fftshift(fft2(double(I)));%对频谱图进行平移处理D=ifft2(B);%傅立叶反变换subplot(1,3,1),imshow(I)%画图像原始图subplot(1,3,2),imshow(log(abs(B)+1),)%画直接变换频谱图s
7、ubplot(1,3,3),imshow(abs(D),)%显示反傅立叶变换后的图像显示效果如下:图2图像的傅立叶变换及其反变换实现小丑图像去除噪音:小丑图象除噪，是我经过多次探索才辛苦的弄出来的A=imread('d:小丑.jpg');I=rgb2gray(A);subplot(2,2,1),imshow(I);B=fftshift(fft2(double(I)));subplot(2,2,2),imshow(log(abs(B)+1),);%如图对B里面的数据进行如下修改:(x,y)={(44,52),(107,
8、41),(22,89),(87,76)}"这些点是我用QQ抓图定位"对上面(x,y)这些点为中心N8领域D8=2的幅值全部置为0%图3显示的是对第一个点(44，52)处理数据subplot(2,2,3),imshow(log(abs(B)+1),);D=ifft2(B);subplot(2,2,4),imshow(abs(D),)
