数字图像处理图像缩放的VC实现

图像缩放的VC+实现实验要求1. 试用最邻近插值和双线性插值对所给图像进行缩小和放大操作。先缩小12倍，再放大12倍。2. 描述图像的变化，分析变化的原因，并给出对两种插值方法的评价。实验原理1. 最邻近插值（近邻取样法）最邻近插值是一种简单的插值算法，也称为零阶插值。它输出的像素灰度值就等于距离它映射到的位置最近的输入像素的灰度值。对于通过反向变换得到的的一个浮点坐标，对其进行简单的取整，得到一个整数型坐标，这个整数型坐标对应的像素值就是目的像素的像素值，也就是说，取浮点坐标最邻近的左上角点（对于DIB是右上角，因为它的扫描行是逆序存储的）对应的像素值。可见，最邻近插值简单且直观，但得到的图像质量不高；当图像中包含像素之间灰度级有变化的细微结构时，最邻近查值法会在图像中产生人为的加工痕迹。2. 双线性插值 对于一个目的像素，设置坐标通过反向变换得到的浮点坐标为(i+u,j+v)，其中i、j均为非负整数，u、v为0,1)区间的浮点数，则这个像素得值 f(i+u,j+v) 可由原图像中坐标为 (i,j)、(i+1,j)、(i,j+1)、(i+1,j+1)所对应的周围四个像素的值决定，即：f(i+u,j+v) = (1-u)(1-v)f(i,j) + (1-u)vf(i,j+1) + u(1-v)f(i+1,j) + uvf(i+1,j+1)其中f(i,j)表示源图像(i,j)处的的像素值，以此类推。这就是双线性内插值法。双线性内插值法计算量大，但缩放后图像质量高，不会出现像素值不连续的的情况。由于双线性插值具有低通滤波器的性质，使高频分量受损，所以可能会使图像轮廓在一定程度上变得模糊。概要设计1. 本实验的系统软件整体流程图如下： 本实验使用的平台是VC6.0，使用VC+实现的图形用户界面。是否退出程序是否继续处理其他图片保存处理后的图片最近邻插值放大/缩小双线性插值放大/缩小设置放大缩小倍数打开图片2功能模块的划分和描述由上面的程序整体流程图可以看出，本程序主要分为下面几个模块：1. 读取灰度图像模块：该模块主要实现读入一幅灰度图像。2. 设置放大缩小倍数模块：可以自由设置图片要放大和缩小的比例，且设置了最大只能放大15倍，最小缩小0.001倍，不过可以多次放大和缩小。3. 最近邻插值模块：实现最近邻插值放大和缩小。4. 双线性插值模块：实现双线性插值放大和缩小。5. 保存处理后的图片模块：实现图片的保存。详细设计由于读取灰度图像模块、保存处理后的图片模块可以由系统自带的类CFileDialog类实现，因此本实验的重点是放在设置放大缩小倍数模块、最近邻插值模块、双线性插值模块这三个模块，以下是关于这三个模块的实现：1 设置放大和缩小比例系数本程序中新建了一个继承自对话框类CDialog的类SetClass用来接收用户输入的放大缩小比例的值，用户通过菜单的设置项来设置该系数，该对话框的界面如下：其中的X和Y分别是用来设置图片的横向和纵向放大系数。用户设置的变量的值存放在文档类的四个浮点型变量m_BXZoomRatio（横向放大倍数，初值为1.25），m_BYZoomRatio（纵向放大倍数，初值为1.25），m_SXZoomRatio（横向缩小系数，初值为0.8），m_SYZoomRatio（纵向缩小系数，初值是0.8）中。其程序流程图如下：CMyDIPDoc:OnZoombig()放大倍数菜单-设置缩小倍数CMyDIPDoc:OnZoomsmall()其代码实现如下：/放大系数设置菜单项的响应函数void CMyDIPDoc:OnZoombig() / TODO: Add your command handler code hereSetClass setclass; /设置对话框if (setclass.DoModal() = IDOK) /将用户设置值存到文档类中的相应变量中m_BXZoomRatio = setclass.m_getx;m_BYZoomRatio = setclass.m_gety;/缩小系数设置菜单项的响应函数void CMyDIPDoc:OnZoomsmall() / TODO: Add your command handler code hereSetClass setclass;if (setclass.DoModal() = IDOK) /将用户设置值存到文档类中的相应变量中m_SXZoomRatio = setclass.m_getx;m_SYZoomRatio = setclass.m_gety;2. 最近邻插值模块和双线性插值模块本程序实现这两个模块功能代码是在一个全局函数中实现的，其原型为：HGLOBAL WINAPI ZoomDIB(LPSTR lpDIB, float fXZoomRatio, float fYZoomRatio, BOOL tag)；其中lpDIB是指向需处理的设备无关位图的指针，fXZoomRatio和fYZoomRatio分别为横向和纵向缩放系数，tag为一标记变量（tag=1时执行最近邻插值操作，tag=0时执行双线性插值操作）。程序流程图如下：OnMenuitem32779()放大最近邻插值OnMenuitem32778()缩小ZoomDIB()OnMenuitem32780()放大双线性插值OnMenuitem32781()缩小其实现代码如下：HGLOBAL WINAPI ZoomDIB(LPSTR lpDIB, float fXZoomRatio, float fYZoomRatio, BOOL tag)/ 源图像的宽度和高度LONGlWidth;LONGlHeight;/ 缩放后图像的宽度和高度LONGlNewWidth;LONGlNewHeight;/ 缩放后图像的宽度（lNewWidth，必须是4的倍数）LONGlNewLineBytes;/ 指向源图像的指针LPSTRlpDIBBits;/ 指向源象素的指针LPSTRlpSrc;/ 缩放后新DIB句柄HDIBhDIB;/ 指向缩放图像对应象素的指针LPSTRlpDst;/ 指向缩放图像的指针LPSTRlpNewDIB;LPSTRlpNewDIBBits;/ 指向BITMAPINFO结构的指针（Win3.0）LPBITMAPINFOHEADER lpbmi;/ 指向BITMAPCOREINFO结构的指针LPBITMAPCOREHEADER lpbmc;/ 循环变量（象素在新DIB中的坐标）LONGi;LONGj;/ 象素在源DIB中的坐标LONGi0;LONGj0;doublex,y;/ 图像每行的字节数LONG lLineBytes;/ 找到源DIB图像象素起始位置lpDIBBits = :FindDIBBits(lpDIB);/ 获取图像的宽度lWidth = :DIBWidth(lpDIB);/ 计算图像每行的字节数lLineBytes = WIDTHBYTES(lWidth * 8);/ 获取图像的高度lHeight = :DIBHeight(lpDIB);/ 计算缩放后的图像实际宽度/ 此处直接加0.5是由于强制类型转换时不四舍五入，而是直接截去小数部分lNewWidth = (LONG) (:DIBWidth(lpDIB) * fXZoomRatio + 0.5);/ 计算新图像每行的字节数lNewLineBytes = WIDTHBYTES(lNewWidth * 8);/ 计算缩放后的图像高度lNewHeight = (LONG) (lHeight * fYZoomRatio + 0.5);/ 分配内存，以保存新DIBhDIB = (HDIB) :GlobalAlloc(GHND, lNewLineBytes * lNewHeight + *(LPDWORD)lpDIB + :PaletteSize(lpDIB);/ 判断是否内存分配失败if (hDIB = NULL)/ 分配内存失败return NULL;/ 锁定内存lpNewDIB = (char * ):GlobalLock(HGLOBAL) hDIB);/ 复制DIB信息头和调色板memcpy(lpNewDIB, lpDIB, *(LPDWORD)lpDIB + :PaletteSize(lpDIB);/ 找到新DIB象素起始位置lpNewDIBBits = :FindDIBBits(lpNewDIB);/ 获取指针lpbmi = (LPBITMAPINFOHEADER)lpNewDIB;lpbmc = (LPBITMAPCOREHEADER)lpNewDIB;/ 更新DIB中图像的高度和宽度if (IS_WIN30_DIB(lpNewDIB)/ 对于Windows 3.0 DIBlpbmi-biWidth = lNewWidth;lpbmi-biHeight = lNewHeight;else/ 对于其它格式的DIBlpbmc-bcWidth = (unsigned short) lNewWidth;lpbmc-bcHeight = (unsigned short) lNewHeight;/ 针对图像每行进行操作for(i = 0; i lNewHeight; i+)/ 针对图像每列进行操作for(j = 0; j = 0) & (j0 = 0) & (i0 lHeight)/ 指向源DIB第i0行，第j0个象素的指针lpSrc = (char *)lpDIBBits + lLineBytes * (lHeight - 1 - i0) + j0;/ 复制象素*lpDst = *lpSrc;else/ 对于源图中没有的象素，直接赋值为255* (unsigned char*)lpDst) = 255;else/否则，执行下面的双线性插值代码double m = x - LONG(x);/X方向的小数部分double n = y - LONG(y);/Y方向的小数部分unsigned char rd, ld, lu, ru;/r：右；l：左；d：下；u：上unsigned char pix; ld = *(char *)lpDIBBits + lLineBytes * (lHeight - 1 - LONG(y) + LONG(x);/左下角点的像素值rd = *(char *)lpDIBBits + lLineBytes * (lHeight - 1 - LONG(y) + LONG(x) + 1);/ 右下角点的像素值lu = *(char *)lpDIBBits + lLineBytes * (lHeight - 1 - LONG(y) - 1) + LONG(x);/ 左上角点的像素值ru = *(char *)lpDIBBits + lLineBytes * (lHeight - 1 - LONG(y) - 1) + LONG(x) + 1);/ 右上角点的像素值pix = (1 - m) * (1 - n) * ld + (1 - m) * n * lu + m * n * ru + m * (1 - n) * rd;/双线性插值，得到变换后的像素* (unsigned char*)lpDst) = pix;/将像素值赋给目的图片/ 返回return hDIB;该函数分别在视图类CMyDIPView中下面的函数（也就是相应菜单项的响应函数）中调用：void OnMenuitem32778();/最近邻插值缩小图片void OnMenuitem32779();/最近邻插值放大图片void OnMenuitem32780();/双线性插值放大图片void OnMenuitem32781();/双线性插值缩小图片实验结果运行程序，各模块实验结果如下：读入灰度图片和设置放大缩小系数的结果读入图片：设置放大缩小系数：最近邻插值放大缩小图片结果图片使用最近邻插值缩小12倍：再使用最近邻插值放大12倍：双线性插值放大缩小图片结果缩小12倍：再放大12倍：结果分析上面实验结果分别用最邻近插值和双线性插值对钟表图片进行缩小和放大操作。先缩小12倍，再放大12倍。我们可以看到用最近邻插值方法处理后的图片产生了严重的疵点，边界出现了毛疵，图像质量很差。而用双线性插值处理后的图片边界较为平滑，图像质量较好，但是图中的一些细节也被模糊了。产生上面结果的原因：最近邻插值是取距离源图像像素点最近的点的像素值赋给新图像相应的像素点，因此在新图像中各像素点之间的方差比较大，因此处理后的图像就产生了很多疵点；而双线性插值取的是源图像中对应的四个点像素的平均值作为新图像像素点的值，因此在新图像中各像素点之间的方差比较小，图像就比较平滑，但是这种方法也会产生一些源图像中所不具有的像素值，会导致图像的一些细节的丢失，并且是不可逆的。参考文献1 冈萨雷斯. 数字图像处理（第二版）. 阮秋琦等译. 电子工业出版社. 2009.6(5)2 谢凤英 赵丹培. Visual C+ 数字图像处理. 电子工业出版社. 2008.9(1)