1、Matlab7.0 与数字图像处理绪论作业 1作业要求：用 Matlab 编程实现：1）将附件图片 1 进行颜色转换，得到新的颜色。2）程获得附件图片 2、3 所示的颜色表。3）画出图 4 的像素灰度统计图。参考知识林宇星的参考程序Matlab 矩阵操作1、构造数组的方法：增量发和 linspace(first,last,num)first 和 last 为起始和终止数，num 为需要的数组元素个数。2、构造矩阵的方法：可以直接用 来输入数组，也可以用以下提供的函数来生成矩阵。ones() 创建一个所有元素都为 1 的矩阵，其中可以制定维数，1，2.个变量zeros() 创建一个所有元素都为
2、0 的矩阵eye() 创建对角元素为 1，其他元素为 0 的矩阵diag() 根据向量创建对角矩阵，即以向量的元素为对角元素magic() 创建魔方矩阵rand() 创建随机矩阵，服从均匀分布randn() 创建随机矩阵，服从正态分布randperm() 创建随机行向量horcat C=A,B，水平聚合矩阵，还可以用 cat(1,A,B)vercat C=A;B，垂直聚合矩阵, 还可以用 cat(2,A,B)repmat(M,v,h) 将矩阵 M 在垂直方向上聚合 v 次，在水平方向上聚合 h 次blkdiag(A,B) 以 A，和 B 为块创建块对角矩阵length 返回矩阵最长维的的长度n
3、dims 返回维数numel 返回矩阵元素个数size 返回每一维的长度，rows,cols=size(A)reshape 重塑矩阵， reshape(A,2,6),将 A 变为 26 的矩阵，按列排列rot90 旋转矩阵 90 度，逆时针方向fliplr 沿垂轴翻转矩阵flipud 沿水平轴翻转矩阵transpose 沿主对角线翻转矩阵ctranspose 转置矩阵，也可用 A或 A.，这仅当矩阵为复数矩阵时才有区别inv 矩阵的逆det 矩阵的行列式值trace 矩阵对角元素的和norm 矩阵或矢量的范数，norm(a，1)，norm(a，Inf).normest 估计矩阵的最大范数矢量c
4、hol 矩阵的 cholesky 分解cholinc 不完全 cholesky 分解lu LU 分解luinc 不完全 LU 分解qr 正交分解kron(A,B) A 为 mn，B 为 pq，则生成 mpnq 的矩阵，A 的每一个元素都会乘上 B，并占据 pq 大小的空间rank 求出矩阵的刺pinv 求伪逆矩阵Ap 对 A 进行操作A.P 对 A 中的每一个元素进行操作matlab 冒号的用法总结：http:/ 如何读入多幅 bmp 格式的图片for i=1:100for j=1:25img=imread(d:,num2str(i),filesep,num2str(j),.bmp); %文件
5、位于 d 盘下面组号文件夹下面的文件，注意命名的时候，如果有补零的话，num2str 那里要加参数。%处理endendBMP 数据格式http:/ 格式的数据存储在一个 3 维矩阵中，此 3 维矩阵为 3 层，每层由二维像素矩阵组成。第一层为红色 Red 的灰度第二层为绿色 Green 的灰度第三层为蓝色 Blue 的灰度绘图函数将矩阵以表格输出xlswrite Write Microsoft Excel spreadsheet file (.xls)xlswrite(filename, M)绘制矩阵的三维数据统计图http:/ = meshgrid(-3:.125:3);Z = peaks(
6、X,Y);meshc(X,Y,Z);axis(-3 3 -3 3 -10 5)作业源代码1) 将附件图片 1 进行颜色转换，得到新的颜色。function Homework_of_introductionImage1=imread(E:研究生课题 研究生课程学习数字图像处理绪论作业11.bmp);Red_layer=Image1(:,:,1);%提取 Red 层灰度值Green_layer=Image1(:,:,2);%提取 Green 层灰度值Blue_layer=Image1(:,:,3);%提取 Blue 层灰度值imwrite(Red_layer,E:研究生课题研究生课程学习 数字图像
7、处理绪论作业1Red 层灰度值.bmp);imwrite(Green_layer,E:研究生课题 研究生课程学习数字图像处理 绪论作业 1Green 层灰度值 .bmp);imwrite(Blue_layer,E:研究生课题研究生课程学习数字图像处理 绪论作业1Blue 层灰度值.bmp);RG_mix=Image1;RG_mix(:,:,3)=0;%RG 混合图RB_mix=Image1;RB_mix(:,:,2)=0;%RB 混合图BG_mix=Image1;BG_mix(:,:,1)=0;%BG 混合图imwrite(RG_mix,E:研究生课题研究生课程学习数字图像处理绪论作业1RG
8、混合图.bmp);imwrite(RB_mix,E:研究生课题研究生课程学习数字图像处理 绪论作业1RB 混合图 .bmp);imwrite(BG_mix,E:研究生课题研究生课程学习数字图像处理绪论作业1BG 混合图.bmp);2) 程获得附件图片 2、3 所示的颜色表function Homework_of_introduction2Image2=imread(E:研究生课题 研究生课程学习数字图像处理绪论作业 12.bmp);Image3=imread(E:研究生课题 研究生课程学习数字图像处理绪论作业 13.bmp);layer1=Image2(:,:,1);%提取第一层layer2=
9、Image2(:,:,2);%提取第二层layer3=Image2(:,:,3);%提取第三层imwrite(layer1,E:研究生课题研究生课程学习数字图像处理 绪论作业 12 提取第一层.bmp);%以图像形式输出imwrite(layer2,E:研究生课题研究生课程学习数字图像处理 绪论作业 12 提取第二层.bmp);%以图像形式输出imwrite(layer3,E:研究生课题研究生课程学习数字图像处理 绪论作业 12 提取第三层.bmp);%以图像形式输出xlswrite(E:研究生课题研究生课程学习数字图像处理 绪论作业 12 提取第一层.xls,layer1);%以表格数据形式
10、输出图像数据xlswrite(E:研究生课题研究生课程学习数字图像处理 绪论作业 12 提取第二层.xls,layer2);%以表格数据形式输出图像数据xlswrite(E:研究生课题研究生课程学习数字图像处理 绪论作业 12 提取第三层.xls,layer3);%以表格数据形式输出图像数据layer1=Image3(:,:,1);%提取第一层layer2=Image3(:,:,2);%提取第二层layer3=Image3(:,:,3);%提取第三层imwrite(layer1,E:研究生课题研究生课程学习数字图像处理 绪论作业 13 提取第一层.bmp);%以图像形式输出imwrite(la
11、yer2,E:研究生课题研究生课程学习数字图像处理 绪论作业 13 提取第二层.bmp);%以图像形式输出imwrite(layer3,E:研究生课题研究生课程学习数字图像处理 绪论作业 13 提取第三层.bmp);%以图像形式输出xlswrite(E:研究生课题研究生课程学习数字图像处理 绪论作业 13 提取第一层.xls,layer1);%以表格数据形式输出图像数据xlswrite(E:研究生课题研究生课程学习数字图像处理 绪论作业 13 提取第二层.xls,layer2);%以表格数据形式输出图像数据xlswrite(E:研究生课题研究生课程学习数字图像处理 绪论作业 13 提取第三层.
12、xls,layer3);%以表格数据形式输出图像数据3) 画出图 4 的像素灰度统计图function Homework_of_introduction3Image4=imread(E:研究生课题 研究生课程学习数字图像处理绪论作业 14.bmp);layer1=Image4(:,:,1);%提取第一层layer2=Image4(:,:,2);%提取第二层layer3=Image4(:,:,3);%提取第三层imwrite(layer1,E:研究生课题研究生课程学习数字图像处理 绪论作业 14 提取第一层.bmp);%以图像形式输出imwrite(layer2,E:研究生课题研究生课程学习数字图像处理 绪论作业 14 提取第二层.bmp);%以图像形式输出imwrite(layer3,E:研究生课题研究生课程学习数字图像处理 绪论作业 14 提取第三层.bmp);%以图像形式输出statis_table=zeros(65536,3);k=1;for i=1:256;for j=1:256;statis_table(k,1)=i;statis_table(k,2)=j;statis_table(k,3)=layer1(i,j);k=k+1;endend%灰度值图像输出meshc(statis_table)
