[计算机软件及应用]matlab高级编程与应用-图像处理实验报告.docVIP

图像处理实验报告 自03 张驰昱 2010012028 基础知识 以测试图像的中心为圆心，图像的长和宽中较小值的一半为半径画一个红颜色的圆 问题分析： 构造一个与原图像大小相同的矩阵并按照要求在矩阵中“画”出形状，然后点乘原图像矩阵就可以了。新的构造矩阵在需要保留原图像的地方须置为1。 具体实现： load(hall.mat); hall=im2double(hall_color); %把8位无符号整形线性转化为0到1之间的数 [xx,yy,~]=size(hall); %得到图像长宽 R=min(xx,yy)/2; %设半径R为长宽中的较小值的一半 circle=ones(xx,yy); %申请一个与原图像相同大小的1矩阵 r=[0:2*pi/400:2*pi]; %一个完整圆周的弧度向量 x=floor(xx/2+1+R*cos(r)); %完整圆周的x坐标向量 y=floor(yy/2+1+R*sin(r)); %完整圆周的y坐标向量，配合之前x坐标值就是圆周坐标了 circle(xx*(y-1)+x)=0; %设圆周上每个坐标点的值为0，至此一个画着黑圆的矩阵构造完成 %下面是把这个构造好的矩阵点乘原图像矩阵，???在R层矩阵加上一个“白圆矩阵”就是红圆了 hall_color_a(:,:,1)=hall(:,:,1).*circle+~circle; hall_color_a(:,:,2)=hall(:,:,2).*circle; hall_color_a(:,:,3)=hall(:,:,3).*circle; imwrite(hall_color_a,hall_color_a.jpg,JPEG); 代码优化： 本题花了我一点脑筋把整个圆周的矩阵坐标构造的两个向量中去，免去了循环，大大提高程序效率，基本是瞬间运行完成。 运行结果： 将测试图像涂成国际象棋状的“黑白格”的样子，其中“黑”即黑色，“白”则意味着保留原图。 问题分析： 继续前面的构造新矩阵点乘原图像矩阵的思想。我用零矩阵和一矩阵拼出一个大的棋盘矩阵。0元素的地方就是涂黑。 具体实现： load(hall.mat); hall=im2double(hall_color); %把8位无符号整形线性转化为0到1之间的数 [xx,yy,~]=size(hall); %得到图像长宽 x0=xx/8;y0=yy/8; %每个小格的长宽，除以8是因为棋盘是8x8的 a=[ones(x0,y0),zeros(x0,y0)]; %黑与白是棋盘的基本单元，黑与白即0和1 chess=repmat([a;fliplr(a)],4,4); %拼接加重复，很容易构造出一个棋盘型的01矩阵 hall_color_b(:,:,1)=hall(:,:,1).*chess; %RGB三层都乘上这个棋盘矩阵 hall_color_b(:,:,2)=hall(:,:,2).*chess; hall_color_b(:,:,3)=hall(:,:,3).*chess; imwrite(hall_color_b,hall_color_b.jpg,JPEG); 代码优化： 这道问题的逻辑并不复杂，但是我始终想用矩阵做进一步简化，脱离一个个数值判断再运算的苦海。这个程序我简化了很多遍，一遍遍的优化逻辑和循环，最后发现这个二维加判断的问题在matlab中可以一层循环都不用，程序效率极高，也是瞬间运行完成。 运行结果：  图像压缩编码 图像的预处理是将每个像素灰度值减去128 ，这个步骤是否可以在变换域进行？请在测试图像中截取一块验证你的结论。 问题分析： 可以。DCT变换就是矩阵乘法，矩阵乘法有分配率，先乘再加和先加再乘是一样的。 具体实现： clear;clc; load hall.mat; load JpegCoeff.mat;%文件读入过程在以下代码展示中均省略，因为都一样 s=hall_gray(1:8,1:8); %取一小块 c1=dct2(s-128); %先减再变换 c2=dct2(s)-dct2(128*ones(8,8)); %先变换再减 err1=c1-c2 %差值 运行结果： 差值err1是10-12量级的，可以忽略 请编程实现二维DCT ，并和MATLAB 自带的库函数dct2比较是否一致 问题分析： 一维DCT变换如下： 二维情况只要再右乘转置矩阵即可：C=DPTDT 所以问题的关键即构造D矩阵，而D矩阵可以用kron函数快速构造出。 具体实现： function c=my_dct2(p) [N,~]=size(p); D(1,:)=ones(1,N)/(N^0.5); %第一行要分开构造 D(2:N,:)=(2/N)^0.5*cos(pi/2/N*kron([1:N