中国地质大学数字图像处理报告.docxVIP

中国地质大学（武汉）数字图像处理实习报告书 姓 名: 专 业: 班 级: 学 号:指导老师: 第1题： 实验内容及目的: 根据灰度图象得到一副彩色图像（变换函数自定） ,分别显示1幅彩色图象的 R G B H, S, H, S, I分量（每个分量也各用 8 bit 分量（每个分量用 8 bit表示），和这幅彩色图象的 表木）。 通过本题学会灰度图和彩色图之间的转换，了解 RG陛间和HIS空间之间的关系。 实验原理: 图像由若干个像素点组成。灰度图像矩阵的每个元素对应一个灰度值， 伪彩色图像矩阵 的每个元素对应 R、G、B分量的值。将一幅灰度图变为一幅彩色图相当于在灰度值与 RGB 分量之间建立函数关系。HIS彩色空间由H （色调）、S （饱和度）、I （辉度）分量构成。RGB 分量和HIS分量的对应关系如下： I KR G B)/3; S -1 -ymin(R,G,B); 1 - -[(R-G) (R-B)] 【- cosq-2 ], .(R — G)2 (R — B)(G -B) e,G 至 B;2n -6,G B. 根据对应函数将 RGB分量转化为HIS分量并显示，即可得到所需的 HIS分量图。 思路分析： 用MATLAB处理图像首先要将灰度图像读入。读入后将图像矩阵中的元素转化为 uint8 类型处理。将灰度图变为彩色图就是在灰度值与 R、G、B分量之间建立关系，所以要先生 成R、G、B分量的矩阵，这就要知道原灰度图的大小（用 size（）函数实现）。建立R、G、 B分量的矩阵后，采用循环根据原图像的灰度值对其进行赋值。最后用 cat 0函数合并三 个分量得到生成的彩色图。 读入彩色图像并转为 uint8类型，根据彩色图像的构成将 R、G、B分量从图像矩阵中 分量出来并显示。 显示图像的HIS分量时，首先将图像转化为 double类型处理。然后分离出 RGB分量， 根据RGB分量和HIS分量之间的转化关系得到 HIS分量，此时的 HIS分量时double类型， 所以将这三个分量合并后转化为 8bit类型。将转化后的矩B^再分离出 HIS分量即可得到8bit 类型的HIS分量图。 代码: clear all close all %读入图像% 显示原始图像 %读入图像 % 显示原始图像 %等图像车专化为int8类型 %求得矩阵的维度 subplot(121); imshow(X);title( 灰度图像’)； X=im2uint8(X); [M,N]=size(X); %生成RGB矩阵 R=zeros(M,N); G=zeros(M,N); B=zeros(M,N); %建立灰度值与RGB分量之间的变换关系 for p=1:M for k=1:N if (X(p,k)128) R(p,k)=128; elseif (X(p,k)192) R(p,k)=4*X(p,k)-511; else R(p,k)=255; end if (X(p,k)64) B(p,k)=255; elseif (X(p,k)128) B(p,k)=511-4*X(p,k); else B(p,k)=0; end if (X(p,k)64) G(p,k)=4*X(p,k); elseif (X(p,k)192) G(p,k)=255; else G(p,k)=1023-4*X(p,k); end end end Z=cat(3,R,G,B); %合并三个分量 Subplot(122); imshow(Z); %显示伪彩色图 title(‘伪彩色图像’)； %%%%%%%%%%%%% %显示一巾g图的RGB分量和HIS分量 clear all ; close all ; X=imread( tower.jpg ); X=im2uint8(X); %等图像车专化为int8 类型 figure(1)； title(‘原始图像’)；％显示原始图像 imshow(X); XR=X(:,:,1); %^离出 RGB分量 XG=X(:,:,2); XB=X(:,:,3); figure(2); subplot(2,3,1); imshow(XR); %显示RGB分量图 title( R 分量’)； subplot(2,3,2); imshow(XG); title( G 分量’)； subplot(2,3,3); imshow(XB); title( B 分量’)； %HSI分量 X=im2double(X); XR=X(:,:,1); XG=X(:,:,2); XB=X(:,:,3); %根据RGB分量和HIS分量间的变换关系求得 HIS分量 num=0.5*((XR-XG)+(XR-XB)); den=sqrt((XR-XG).A2+(XR-XB).*