《图像超分辨率重建图像处理课程设计》由会员分享，可在线阅读，更多相关《图像超分辨率重建图像处理课程设计（14页珍藏版）》请在金锄头文库上搜索。
1、.目录1 课程设计目的12图像处理系统设计容及要求12.1设计容12.2设计要求13 设计方案14 功能模块的具体实现14.1 空域插值放大的方法14.1.1 最邻近插值算法14.1.2 双线性插值算法14.1.3 双三次插值算法14.2 频域重建的方法14.2.1 DCT变换的介绍14.2.2 DCT放大图像放大算法原理14.3 频域分块重建的方法14.4 同态滤波器滤波处理14.4.1 同态滤波器原理14.4.2 同态滤波函数确实定15 总结与体会1参考文献1附录1. . .1 课程设计目的MATLAB7.0软件。MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及
2、交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进展有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言的编辑模式，代表了当今国际科学计算软件的先进水平。通过用MATLAB对图像进展处理，以实现以下目的。1.培养严谨的科学态度，正确的设计思想，科学的设计方法和良好的工作作风。2.培养独立思考的能力，独立检索资料、阅读文献、综合分析、计算机应用、数据及文字处理等能力。3.培养综合运用根底理论、根本知识的能力。通过课程设计得到工程
3、设计的初步锻炼。2图像处理系统设计容及要求2.1设计容图像超分辨率重建是利用低质量或低分辨率图像来产生高质量或高分辨率图像的技术，重建包括空域方法和频域方法。本设计要求用插值技术提高图像的分辨率。1利用插值技术将原始图像在空域放大2倍。2在DCT域放大原始图像2倍，设计滤波器在DCT域增强图像的高频信息。3对图像分块进展DCT变换，在DCT域对子图像进展放大和滤波增强高频信息。4比拟上述三种图像重建结果，设计软件界面。5设计方案、编写代码实现上述功能。2.2设计要求1利用数字图像处理技术，以MATLAB为平台，建立一个实现设计主题的简易处理系统。2能显示输入图像、中间图像和重建的图像。3程序代
4、码要有注释说明，调用MATLAB函数要清楚并理解函数的功能、使用围，在设计说明书中要写清楚函数的功能和参数意义。4完成设计说明书一份。5刻苦钻研，勤于思考，勇于实践，独立完成课程设计任务。6遵守纪律，在指定地点进展课程设计。7掌握有关课程的根本理论和根本知识。概念清楚，方案合理，数据可靠，计算正确，运行良好，图纸图表符合标准，设计说明书论文撰写规，辩论中答复下列问题正确。3 设计方案根据课程设计题目的要求，设计界面如图3-1所示：图3-1 设计界面本设计分为空域放大和DCT域放大两局部，空域放大主要工作有：在空域对现有的传统插值算法分别进展了研究与仿真实验，包括最近邻域插值，双线性插值，双三次
5、插值等，这些插值方法均是通过低通滤波，滤除和过滤图像数据中的高频信息。所以这些插值基函数对边缘和纹理信息都比拟丰富的图像的插值效果不是特别理想。DCT域放大主要工作有：通过DCT变换实现了由空间域到频域的转换，通过对频域处理可以方便的实现空间域较难实现的处理。而空间域与频域又存在一定的联系，为数字图像的处理提供了另一种方法。该算法在对整块图像进展处理时，尽管采用了增强系数对图像亮度效果进展补充，但对整幅图像高频局部预测采用填零方式，在图像像素位数增大即图像信息量增大时这种预测精度不如对图像分块处理后高，且基于JPEG格式图像多采用分成88子块分块压缩编码，对上述算法进展改良。改良后的算法，将原
6、始图像数据切割成接近88大小子块，对每一子块分别实施DCT放大算法。改良后的算法如下：对原始图像进展分块，然后对每一子块运用DCT图像放大算法，最后合并处理所有的子块。系统整体设计如图3-2所示。空域重建频域分块重建原始图低分辨率图像频域重建观察、比拟三种图像重建的结果设计软件界面图3-2 系统整体设计图4 功能模块的具体实现4.1 空域插值放大的方法4.1.1 最邻近插值算法这是最简单的算法，每一个输出像素都赋给输入图象中与其最邻近的采样点的值。插值核函数是：h(*)=1 0*0.5h(*)=0 0.5*在所有的插值方法中，这种插值方法速度最快，早期的应用比拟普遍，然而当图像中包含灰度有变化
7、的细微构造时，最近邻插值法会在图像中产生人工的痕迹。图像的边缘阶梯失真现象比拟明显。其实现效果如图4-1所示。图4-1 最邻近插值算法实现效果图在程序中可以直接调用函数也可自行编写。如自行编写，局部源程序如下：width = K * nrows; height = K * ncols;J = uint8(zeros(width,height);widthScale = nrows/width;heightScale = ncols/height;for * = 5:width - 5 for y = 5:height - 5* = * * widthScale; yy = y * height
8、Scale;if (*/double(uint16(*) = 1.0) & (yy/double(uint16(yy) = 1.0) J(*,y) = I(int16(*),int16(yy);else % a or b is not integera = double(round(*); % (a,b) is the base-dotb = double(round(yy);J(*,y) = I(a,b); % calculate J(*,y)endendendimwrite(J, lena2.jpg, jpg);figure;imshow(J);4.1.2 双线性插值算法设f ( * ,
9、y ) 为2 个变量的函数, 其在单位正方形顶点的值。假设希望通过插值得到正方形任意点的f ( * , y ) 值。可以令由双线性方程来定义的一个双曲抛物面与4 个点拟合。利用公式实现插值：双线性插值是对待插值象素周围的4个邻近像素的灰度按照距离进展加权平均，实质上是属于一阶插值。双线性插值的平滑作用有可能会使图像的细节产生退化，当放大倍数比拟大的时候，这种现象更加明显。同时，双线性插值的斜率不连续也会产生不希望的结果。其实现效果如图4-2所示。图4-2 双线性插值算法实现效果图在程序中可以直接调用函数也可自行编写。如自行编写，局部源程序如下：width = K * nrows; height
10、 = K * ncols;J = uint8(zeros(width,height);widthScale = nrows/width;heightScale = ncols/height;for * = 5:width - 5 for y = 5:height - 5* = * * widthScale; yy = y * heightScale;if (*/double(uint16(*) = 1.0) & (yy/double(uint16(yy) = 1.0) % if a and b is integer,then J(*,y) - I(*,y)J(*,y) = I(int16(*)
11、,int16(yy);else a = double(uint16(*); % (a,b) is the base-dotb = double(uint16(yy);*11 = double(I(a,b); % *11 - I(a,b)*12 = double(I(a,b+1); % *12 - I(a,b+1)*21 = double(I(a+1,b); % *21 - I(a+1,b)*22 = double(I(a+1,b+1); % *22 - I(a+1,b+1) J(*,y) = uint8( (b+1-yy) * (*-a)*21 + (a+1-*)*11) + (yy-b) *
12、 (*-a)*22 +(a+1-*) * *12) ); % calculate J(*,y)endendendimwrite(J, lena2.jpg, jpg);figure;imshow(J);4.1.3 双三次插值算法双三次插值是用待插枝点周围16个点作为参考像素值的一种三阶插值方法，典型的双三次插值核函数是：这里参考值a不同文献中取值不同，a=一1，a=05，a=075等等。其实现效果如图4-3所示。图4-3 双三次插值算法实现效果图4.2 DCT域插值放大的方法4.2.1 DCT变换的介绍离散余弦变换是从二种特殊形式的傅里叶变换转化过来的，是一种性能很好的正交变换方式。离散余弦变换
13、本质上仍然是离散傅立叶变换，二者在频域本质上是一样的。离散余弦变换因其是一种实数变换，其变换矩阵的基向量很好地描述了人类视觉的相关性，接近于最正确变换。因而DCT在图像处理中有很广泛的应用，并成为一些静态图像和视频压缩国际标准的根本处理模块，因而采用DCT变换可以很方便地应用于压缩域图像和视频中。数字图像可以通过傅里叶变换、离散余弦变换等由空间域转换到频域中表示，通过对频域的处理可以方便实现空间域较难实现的处理。而空间域和频域又存在一定的联系，为数字图像的处理提供了另一种方法。尺度变换的物理含义就放大意义来讲，如果信号在时域进展扩展，即当Oa1时，则其频谱将在频域进展压缩，信息更集中于低频，高频局部所含信息很少，同时幅值也会增大。根据这一原理，可以对图像频域进展处理以到达图像空间域放大的操作，即将源图像的频域数据作为目标图像频域数据的低频局部，而对于高频进展填零预测，幅值乘以相应的增强系数，再反变换回空间域可实现图像的放大。傅里叶变换由于频域复数运算，在处理及运算上有些复杂，而离散余弦变换是实数变换，可以方便地对应空间域图像。4.2.2 频域重建的方法根据图像的表示方法以及以上原理，我们可以仅在频域进展处理就可完成对图像的放大操作。对原图像进展DCT正变换得到图像频域数据F
