中南大學数学院数字图像处理报告1.docVIP

数据图像处理报告 学生姓名 学 号 学 院 数学科学与计算技术学院 专业班级 本科生院制 图像处理实验（一）图像的幅度谱图及频域滤波 快速傅氏变换（FFT），是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的能克服时间域与频率域之间相互转换的计算障碍，在光谱、大气波谱分析、数字信号处理等方面有广泛应用。 f1(x,y)（128×128大小，暗处=0，亮处=255） 同屏显示原图f1和FFT(f1)的幅度谱图 f1=ones(128)*0.; for i=25:100 for j=60:70 f1(i,j)=255; end end figure(1); subplot(2,3,1); imshow(f1); title(f1原始图); f3=fft2(f1); subplot(2,3,2); imshow(abs(f3)); title(FFT(f1)幅谱图); 令f2(x,y)=，重复以上过程,并简述二者幅度谱异同的理由 for i=1:128 for j=1:128 f2(i,j)=((-1)^(i+j))*f1(i,j); end end subplot(2,3,3); imshow(f2); title(f2的原始图); f=fft2(f2); subplot(2,3,4); imshow(abs(f)); title(FFT(f2)幅谱图); 理由：因为f2是f1通过平移得到，而且正好平移到f1中部，若再做快速傅里叶 变换，则f2相比于f1,变换后主要能量（低频分量）集中在频率平面的中心。 对该图进行动态范围压缩以增强图象 利用DCT变化进行动态压缩 figure; I = im2double(f1); T = dctmtx(8); //求8*8模板的dct变换矩阵 B = blkproc(I,[8 8],P1*x*P2,T,T); //施DCT变换B=TIT’ mask = [1 1 1 1 0 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]; B2 = blkproc(B,[8 8],P1.*x,mask); //只取10个有效点 I2 = blkproc(B2,[8 8],P1*x*P2,T,T); //IDCT变换 subplot(1,2,1); imshow(f1); title(原图像); subplot(1,2,2); imshow(I2); title(变化后原图像); //显示原图与idct变换后的图 运行结果 将f2(x,y)顺时针旋转45度得到f 3(x , y) ，试显示FFT(f 3 )的幅度谱，并与FFT(f2)的幅度谱进行比较 f4=imrotate(f1,315, bilinear, crop); //双线性插值法旋转图像45度,并剪切图像,使其和原图像大小一致 subplot(2,3,5) imshow(f4) title(顺时针旋转图像45^{0}) f5=fft2(f4); subplot(2,3,6); imshow(abs(f5)); title(旋转图像后的幅谱图); 以上程序运行的结果显示如下： 图像处理实验（二）数字图像增强 灰度变换是图像增强的一种重要手段，使图像对比度扩展，图像更加清晰，特征更加明显。灰度级的直方图给出了一幅图像概貌的描述，通过修改灰度直方图来得到图像增强。 由于本文采用的是张彩色图片，首先需通过算法将其变成黑白 a = imread(tree