假彩色编码小学期论文

1、 课程设计报告 设计题目 基于MATLAB图像假彩色编码 专 业 电子信息工程 班 级 一 班 学生姓名 黄 涛 实验成绩 指导教师 平 子 良 5年6月基于MATLAB图像假彩色编码黄涛(通信系电子信息工程1班)摘要: 通过对原始图像进行假彩色编码,可以提高图像的可辨识灰度差异,再现图像细节.该项技术广泛应用在工程实践和数字图像处理。 本文介绍了一些经典的图像假彩色编码原理与思想,包含密度分层假彩色编码频率假彩色编码等。其中主要介绍了单色密度三级-RGB复色密度、频域编码并以此进行试验，并在几种常规方法后进行自定义的基于MATLAB图像的假彩色编码。关键字: 密度分层法假彩色编码； 频域假彩

2、色编码； MATLABFalse color coding based on MATLAB image黄涛(通信系电子信息工程1班)Huang Tao (Communication Electronic Information Engineering 1)摘要:Abstract:通过对原始图像进行假彩色编码,可以提高图像的可辨识灰度差异,再现图像细节.该项技术广泛应用在工程实践和数字图像处理。By using false color coding of the original image, the image can be improved to recognize the gray dif

3、ference and reproduce the image details. This technology is widely used in engineering practice and digital image processing.本文介绍了一些经典的图像假彩色编码原理与思想,包含密度分层假彩色编码频率假彩色编码等。其中主要介绍了单色密度三级-RGB复色密度、频域编码并以此进行试验，并在几种常规方法后进行自定义的基于MATLAB图像的假彩色编码。In this paper, some classical image false color coding is introduc

4、ed, including density, pseudo color coding, frequency false color coding and so on. The paper mainly introduces the monochrome density three RGB composite color density, frequency domain coding and thus test, and custom after several conventional methods of based on MATLAB image pseudo color coding.

5、关键字: 密度分层法假彩色处理 频域假彩色编码原理 MATLABKeywords: density stratification method false color processing frequency false color coding principle MATLAB1. 前言在实际生活中有着很多的黑白照片，比如B超照片、X射线底片、金相图片等，但人类视觉能能同时分辨的灰度最多不过1520个层次，这就影响了人们对这些图像的分析处理。然而人眼的视锥细胞对彩色的分辨力却相当敏锐，以至于可以辨认出数以千万记的不同颜色。假彩色编码是通过光学手段给一张黑白图片中的不同黑白层次赋予不同的颜色，

6、以提高人眼对黑白底片的辨识能力。但由于传统的光学手段操作繁琐，容易产生误差，因而借助于MATLAB7.1软件进行仿真实验。MATLAB是MathWorks公司于1982年推出的一套高性能的数值计算和可视化数学软件,它集数值分析、矩阵运算、信号处理和图形显示于一体,构成了一个方便的、界面友好的用户环境,被誉为“巨人肩上的工具”。MATLAB中除主包外，还包含许多功能各异的工具箱用于解决各个领域的特定问题。它的工具箱主要有通信、控制系统、滤波器设计、图像处理、非线性控制设计、系统识别、神经网络、最优化、模糊逻辑、信号处理、鲁棒控制、统计等。借助于这些工具箱，用户可以非常方便地迸行分析、计算及设计工

7、作。由于MATLAB有着上述种种的便捷之处，借助其进行假彩色编码试验将会很便利、快捷。因此借助其对单色密度、复色密度、频域编码将是高效便捷的。 随着信息技术的发展，图像信息被广泛应用于多媒体通信和计算机系统中。但是图像具有庞大的数据量，不便于储存和传输，必须采用合适的方法对其进压缩，因此有必要对图像压缩标准编码进行研究。在众多的图像编码标准中，假彩色编码广泛应用于静态图此案的压缩，在数字图像压缩领域得到广泛的应用。2.灰度图像假彩色编码原理2.1 密度分层假彩色编码密度分层法假彩色处理又称为灰度分层法，是假彩色处理技术中最基本、最简单的方法。在图像处理中，图像灰度变换属于点运算范畴。点运算（P

8、oint Operation）的概念是当算子T的作用域是以每一个单个像素为单位，图像的输出g(x,y)只与位于(x,y)处得输入f(x,y)有关，实现的是像素点到点的运算处理时，称这种运算为“点运算”。点运算的表达式为： s=T（r）或者 g(x,y)=Tf(x,y)式中：r和s分别为输入、输出像素的灰度级；T为灰度变换函数的映射关系。通过上述式子可以将原图像（x，y）处的灰度f(x,y)变为Tf(x,y),T算子描述了输入灰度级和输出灰度级之间的映射关系。 密度假彩色处理是对图像的灰度值动态范围进行分割，使分割后的每个灰度值区间甚至每个灰度值本身对应于某一颜色。具体而言，假定把一幅图像看作一

9、个二维的强度函数，用于一个平行于图像坐标平面的平面（称之为密度切割平面）去切割图像的强度函数，这样图像灰度在分割处被分为上、下两部分，即两个灰度值区间。如果再对每一个区间赋予某种颜色，就可以将原来的灰度图像变换成只有两种颜色的图像。更进一步，如果用多个密度切割平面对图像函数进行分割，那么就可以将图像的灰度动态范围切割成多个区间。应当指出的是，每一灰度值区间赋予何种颜色，是由具体应用决定的，并无规律可言。但相邻灰度值区间的颜色差别不宜太小也不宜太大。太小将无法反映细节上的差异，太大则会导致图像出现不连续性。实际应用中，密度切割平面之间可以使等间隔的，也可以是不等间隔的，而且切割平面的划分也应依据

10、具体的应用范围和研究对象而确定。结论：密度假彩色编码不改变图像的大小，不考虑质量的变化，只改变图像的颜色，突出人想要突出的部分。 2.2 频率假彩色编码在频域假彩色增强时，先把灰度图像 f(x，y)中的不同频率成分别经过 FFT傅里叶变换到频率域。在频率域内，经过三个不同传递特性的滤波器，f(x，y)被分离成三个独立分量，然后对它们进行 IFFT傅里叶逆变换，便得到三幅代表不同频率分量的单色图像。接着对这三幅图像作进一步的附加处理，最后将它们作为三基色分量分别加到彩色显示器的红、绿、蓝显示通道，从而实现频率域分段的假彩色增强。频率域滤波的假彩色增强处理原理如图 2-1图2_1频域滤波的假彩色增

11、强处理框图频域假彩色编码的意义图像经过变换后处理起来比变换前更加简单方便，由于这种变化是对图像函数而言的，所以称为图像变换。现在研究的图像变换基本上都是正交变换，正交变换可以减少数据的相关性，获取图像的整体特点，有利于用较少的数据量表示原始图像，这对图像的分析、存储以及图像的传输都是非常有意义的。图像的傅里叶变换傅里叶变换分为连续函数的傅里叶变换和离散函数的傅里叶变换两种，由于计算机只能处理离散数值，所以连续傅里叶变换在计算机上无法直接使用。为了在计算机上实现傅里叶变换计算，必须把连续的函数离散化，即将连续傅里叶变换转换为离散傅里叶变换。频域滤波频域滤波增强方法是将图像从空间域变换到频域，在图

12、像的频域空间对图像进行滤波处理。根据信号分析理论，傅里叶变换和卷积理论是频域滤波技术的基础，因此，在频域空间的滤波与在空间域滤波一样可以通过卷积运算实现。假定g（x,y）表示函数f（x,y）与线性不变算子h（x,y）进行卷积运算的结果，即 g（x,y）= f（x,y）* h（x,y） 式（2_2）因此可得 G（u,v）= F(u,v)H(u,v) 式（2_3）其中，H（u,v）为传递函数，为突出F（u,v）中的不同频率部分而设计出不同的滤波器。滤波器分为很多种，常见的有高通滤波器，低通滤波器，带通滤波器。低通滤波器，它能减弱或消除傅里叶空间的高频分量，但不影响低频分量。因为高频分量对应图像中的

13、区域边缘等灰度值具有较大较快变化的部分,滤波器将这些分量滤去可使图像平滑。高通滤波器，它能够减弱或消除傅里叶空间的低频分量，但不影响高频分量。因为低频分量对应图像中灰度值缓慢变化的区域，因而与图像的整体特性，如整体对比度和平均灰度值等有关，高通滤波器将这些分量滤去可使图像锐化。带通滤波器，它能够减弱或消除傅里叶空间的特定频率分量，既不是图像中的高频成分也不是低频成分，而是在一个有限的频带范围内。这三种滤波器的特性曲线是三维图形，其二维剖面示意图如图2_2 所示：图 2_2频域低通、高通、带通三种滤波器剖面示意图本实验用滤波器 常见的滤波器的类型有很多，如理想型、巴特沃斯型、指数型等。本次实验中

14、应用的是巴特沃斯型滤波器。巴特沃斯低通滤波器的传递函数为： 式（2_4）式中，D0为截止频率；n为滤波器的阶次。和理想滤波器一样，巴特沃斯低通滤波器的特性曲线同样为三维图型，其剖面示意图如图4-2所示。一般情况下，当H（u,v）下降到最大值1/2时，D（u,v）为截止频率D0。 图2_3 巴特沃斯低通滤波器剖面特性图巴特沃斯滤波器又称为最大平坦滤波器，其通带与带阻之间的过渡比较平坦。因此，巴特沃斯低通滤波器的特点是：在通过频率与截止频率之间没有明显的不连续性，不会出现“振铃”效应，其效果好于理想低通滤波器。这也是本实验采用其的原因。巴特沃斯高通滤波器的传递函数为：式（2_4）巴特沃斯高通滤波器

15、性能的剖面示意图如图4-3所示。截止频率的取值方法也是降到最大值的1/2时取得。该通过频率与截止频率之间也没有明显的不连续性，图像增强后，“振铃”现象不明显。图2_4 巴特沃斯高通滤波器剖面特性图3.试验模拟结果及分析3.1单色密度实验对图像进行了灰度分段，对不同的灰度段赋予了三基色中不同的单色。编码程序：I=imread('C:Documents and SettingsAdministrator桌面平野测试图像2.jpg');J=zeros(256,256,3);I=double(I);m,n=size(I);for i=1:256(m)for j=1:256(n) if

16、0<=I(i,j)&&I(i,j)<=80 R(i,j)=255;G(i,j)=0;B(i,j)=0; else if 80<=I(i,j)&&I(i,j)<=190 R(i,j)=0;G(i,j)=255;B(i,j)=0; else if 190<=I(i,j)&&I(i,j)<=255 R(i,j)=0;G(i,j)=0;B(i,j)=255; end end endendend J(:,:,1)=R;J(:,:,2)=G;J(:,:,3)=B; figure, imshow(uint8(J);程序所得图

17、像：图（3-1）程序修改部分:for i=1:256(m)for j=1:256(n) if 0<=I(i,j)&&I(i,j)<=80 R(i,j)=0; G(i,j)=255; B(i,j)=0; else if 80<=I(i,j)&&I(i,j)<=190 R(i,j)=0; G(i,j)=0;B(i,j)=255; else if 190<=I(i,j)&&I(i,j)<=255 R(i,j)= 255;G(i,j)=0;B(i,j)= 0; 修改后程序所得图像：图（3-2）分析:该程序的修改使得红绿

18、蓝三色的所在位置发生了互换但其色域及颜色并没有发生改变。程序自定义修改部分： m,n=size(I);for i=1:256(m)for j=1:256(n) if 0<=I(i,j)&&I(i,j)<=80 R(i,j)=0; G(i,j)=255; B(i,j)=0; else if 80<=I(i,j)&&I(i,j)<=190 R(i,j)=0; G(i,j)=0;B(i,j)=255; else if 190<=I(i,j)&&I(i,j)<=255 R(i,j)= 255;G(i,j)= 255;B

19、(i,j)= 0; 自定义修改后程序所得图像：图（3-3）分析：自定义修改后的程序所做出的图像与图（3-2）相比色域并未再次互换，但发生了混色现象，这是因为在该区域内程序设置了R(i,j)= 255;G(i,j)= 255，使得该区域中同时存在等量的红色和绿色，发生混色后得到了黄色。3.2 复色密度实验对图像进行了灰度分段，对不同的灰度段赋予了不同量的三基色。编码程序：I=imread('C:Documents and SettingsAdministrator桌面平野测试图像 cat256rgb.png');J=zeros(256,256,3);I=double(I);for

20、 i=1:256for j=1:256 if 0<=I(i,j)&&I(i,j)<=63 R(i,j)=0;G(i,j)=4*I(i,j);B(i,j)=255; else if 63<=I(i,j)&&I(i,j)<=127 R(i,j)=0;G(i,j)=255;B(i,j)=511-4*I(i,j); else if 128<=I(i,j)&&I(i,j)<=191 R(i,j)=4*I(i,j);G(i,j)=255;B(i,j)=0; else if 192<=I(i,j)&&I

21、(i,j)<=255R(i,j)=255;G(i,j)=1025-4*I(i,j);B(i,j)=0; end end end endendend J(:,:,1)=R;J(:,:,2)=G;J(:,:,3)=B; figure, imshow(uint8(J);运行后所得图像： 图（3-4）程序修改部分：for i=1:256for j=1:256 if 0<=I(i,j)&&I(i,j)<=63 R(i,j)=0;G(i,j)=4*I(i,j);B(i,j)=255; else if 63<=I(i,j)&&I(i,j)<=12

22、7 R(i,j)=0;G(i,j)=255;B(i,j)=511-8*I(i,j); else if 128<=I(i,j)&&I(i,j)<=191 R(i,j)=4*I(i,j);G(i,j)=255;B(i,j)=0; else if 192<=I(i,j)&&I(i,j)<=255R(i,j)=255;G(i,j)=1025-8*I(i,j);B(i,j)=0; 修改后程序所的图像： 图（3-5）分析:由于该句程序的改变else if 63<=I(i,j)&&I(i,j)<=127 R(i,j)=0;G

23、(i,j)=255;B(i,j)=511-8I(i,j);使得图像的色域发生改变。3.3 频域编码实验对图像先进行了转化为灰度图像的操作，然后将灰度图像分频处理，使高、中、低频通过各自的滤波器后产生高、中、低频各自的图片，最后将三图进行合成，得到频域假彩色编码图。编码程序：clcI1=imread(' C:Documents and SettingsAdministrator桌面平野测试图像 9.jpg');I=rgb2gray(I1);figure,subplot(231),imshow(I1);subplot(234),imshow(I);M,N=size(I);I=dou

24、ble(I);F=fft2(I);fftshift(F);REDcut=100; %红色高频截止频率GREENcut=200; %绿色低频截止频率BLUEcenter=150; %蓝色频带中心BLUEwidth=100; %蓝色频带宽度BLUEu0=10;BLUEv0=10;for u=1:Mfor v=1:N D(u,v)=sqrt(u2+v2); %频率至中心距离 REDH(u,v)=1/(1+(sqrt(2)-1)*(D(u,v)/REDcut)2); %低通 GREENH(u,v)=1/(1+(sqrt(2)-1)*(GREENcut/D(u,v)2); %高通 BLUED(u,v)=

25、sqrt(u-BLUEu0)2+(v-BLUEv0)2); %蓝色频率至中心 if(BLUED(u,v)=BLUEcenter)&&(BLUED(u,v)=-BLUEcenter) BLUED(u,v)=BLUED(u,v)+1, else BLUED(u,v)=BLUED(u,v); endBLUEH(u,v)=1-1/(1+BLUED(u,v)*BLUEwidth/(BLUED(u,v)2-(BLUEcenter)2)2); %蓝带通endendRED=REDH.* F; %图像红频率REDcolor=ifft2(RED); %图像红色成分,相应于低频GREEN=GREEN

26、H.* F; %图像绿频率GREENcolor=ifft2(GREEN); %图像绿色成分,相应于高频BLUE=BLUEH.* F; %图像蓝频率BLUEcolor=ifft2(BLUE); %图像蓝色成分,相应于中频REDcolor=real(REDcolor)/256; %转换成实数图像GREENcolor=real(GREENcolor)/256;BLUEcolor=real(BLUEcolor)/256;R=abs(REDcolor);B=abs(BLUEcolor);G=abs(GREENcolor);subplot(232),imshow(R);subplot(233),imsho

27、w(G);subplot(235),imshow(B); for i=1:Mfor j=1:N OUT(i,j,1)=REDcolor(i,j); OUT(i,j,2)=GREENcolor(i,j); OUT(i,j,3)=BLUEcolor(i,j);endendOUT=abs(OUT);subplot(236),imshow(OUT);编码程序所得图像：低频成分高频成分原图频域编码假彩色图像中频成分转换成灰度图像 图（3-6）修改原始编码程序：clcI1=imread(' C:Documents and SettingsAdministrator桌面平野测试图像 9.jpg

28、9;);I=rgb2gray(I1);figure,subplot(231),imshow(I1);subplot(234),imshow(I);M,N=size(I);I=double(I);F=fft2(I);fftshift(F);REDcut=1500; %红色高频截止频率GREENcut=2; %绿色低频截止频率BLUEcenter=150; %蓝色频带中心BLUEwidth=100; %蓝色频带宽度BLUEu0=10;BLUEv0=10;for u=1:Mfor v=1:N D(u,v)=sqrt(u2+v2); %频率至中心距离 REDH(u,v)=1/(1+(sqrt(2)-1

29、)*(D(u,v)/REDcut)2); %低通 GREENH(u,v)=1/(1+(sqrt(2)-1)*(GREENcut/D(u,v)2); %高通 BLUED(u,v)=sqrt(u-BLUEu0)2+(v-BLUEv0)2); %蓝色频率至中心 if(BLUED(u,v)=BLUEcenter)&&(BLUED(u,v)=-BLUEcenter) BLUED(u,v)=BLUED(u,v)+1, else BLUED(u,v)=BLUED(u,v); endBLUEH(u,v)=1-1/(1+BLUED(u,v)*BLUEwidth/(BLUED(u,v)2-(BLU

30、Ecenter)2)2); %蓝带通endendRED=REDH.* F; %图像红频率REDcolor=ifft2(RED); %图像红色成分,相应于低频GREEN=GREENH.* F; %图像绿频率GREENcolor=ifft2(GREEN); %图像绿色成分,相应于高频BLUE=BLUEH.* F; %图像蓝频率BLUEcolor=ifft2(BLUE); %图像蓝色成分,相应于中频REDcolor=real(REDcolor)/256; %转换成实数图像GREENcolor=real(GREENcolor)/256;BLUEcolor=real(BLUEcolor)/256;R=a

31、bs(REDcolor);B=abs(BLUEcolor);G=abs(GREENcolor);subplot(232),imshow(R);subplot(233),imshow(G);subplot(235),imshow(B); for i=1:Mfor j=1:N OUT(i,j,1)=REDcolor(i,j); OUT(i,j,2)=GREENcolor(i,j); OUT(i,j,3)=BLUEcolor(i,j);endendOUT=abs(OUT);subplot(236),imshow(OUT);修改原始编码程序后所得图像：低频成分高频成分原图频域码假彩色图像中频成分转换为

32、灰度图像 图（3-7）分析:修改的程序编码为fftshift(F);REDcut=1500; %红色高频截止频率GREENcut=2; %绿色低频截止频率BLUEcenter=150; %蓝色频带中心BLUEwidth=100; %蓝色频带宽度将红色高频截止频率和绿色低频截止频率进行了修改,致使红色绿色的通过率发生了改变,从而发生了图(3-7)的混色现象。4.实验结论随着技术的进步，数字图像处理相关技术的应用也越发的广泛。相对于传统的文字信息，数字图像拥有更大的信息量，更加便于人们接收信息，已经应用到了社会的方方面面。本文首先介绍了基于MATLAB图像的假彩色编码，接着通过MATLAB7.1编

33、程分别实现了单色密度三级-RGB复色密度、频域编码三项实验，对多幅图像进行了仿真实验，并对实验数据进行分析，给出了对应的实验图像。实验结果表明假彩色编码实现简单，在多数的情况下都能得到较好的实验效果。利用MATLAB进行仿真实验，方法简单，误差小，大大提高了假彩色编码的效率和精度。借助基于MATLAB进行假彩色编码仿真解决了人类视觉不能同时分辨出灰度超过1520个层次的图像的问题，使人们对这些图像的分析和处理变得容易。在黑白图片中的不同黑白层次赋予不同的颜色，大大增强了图像的辨识度，就如同在上述的三个实验中所得到的图片所示，清晰明了的将图像的层次分开。经过此次的MATLAB实验课程设计，使得我