数字图像处理-图像缩放的VC++实现.doc

1、图像缩放的 VC+实现实验要求1. 试用最邻近插值和双线性插值对所给图像进行缩小和放大操作。先缩小 12倍，再放大 12 倍。2. 描述图像的变化，分析变化的原因，并给出对两种插值方法的评价。实验原理1. 最邻近插值（近邻取样法）最邻近插值是一种简单的插值算法，也称为零阶插值。它输出的像素灰度值就等于距离它映射到的位置最近的输入像素的灰度值。对于通过反向变换得到的的一个浮点坐标，对其进行简单的取整，得到一个整数型坐标，这个整数型坐标对应的像素值就是目的像素的像素值，也就是说，取浮点坐标最邻近的左上角点（对于 DIB 是右上角，因为它的扫描行是逆序存储的）对应的像素值。可见，最邻近插值简单且直观

2、，但得到的图像质量不高；当图像中包含像素之间灰度级有变化的细微结构时，最邻近查值法会在图像中产生人为的加工痕迹。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)表示源

3、图像(i,j)处的的像素值，以此类推。这就是双线性内插值法。双线性内插值法计算量大，但缩放后图像质量高，不会出现像素值不连续的的情况。由于双线性插值具有低通滤波器的性质，使高频分量受损，所以可能会使图像轮廓在一定程度上变得模糊。概要设计1. 本实验的系统软件整体流程图如下：本实验使用的平台是 VC6.0，使用 VC+实现的图形用户界面。打开图片设置放大缩小倍数最近邻插值放大/缩小双线性插值放大/缩小保存处理后的图片是否继续处理其他图片退出程序否是2功能模块的划分和描述由上面的程序整体流程图可以看出，本程序主要分为下面几个模块：1. 读取灰度图像模块：该模块主要实现读入一幅灰度图像。2. 设置放

4、大缩小倍数模块：可以自由设置图片要放大和缩小的比例，且设置了最大只能放大 15 倍，最小缩小 0.001 倍，不过可以多次放大和缩小。3. 最近邻插值模块：实现最近邻插值放大和缩小。4. 双线性插值模块：实现双线性插值放大和缩小。5. 保存处理后的图片模块：实现图片的保存。详细设计由于读取灰度图像模块、保存处理后的图片模块可以由系统自带的类CFileDialog 类实现，因此本实验的重点是放在设置放大缩小倍数模块、最近邻插值模块、双线性插值模块这三个模块，以下是关于这三个模块的实现：1 设置放大和缩小比例系数本程序中新建了一个继承自对话框类 CDialog 的类 SetClass 用来接收用户

5、输入的放大缩小比例的值，用户通过菜单的设置项来设置该系数，该对话框的界面如下：其中的 X 和 Y 分别是用来设置图片的横向和纵向放大系数。用户设置的变量的值存放在文档类的四个浮点型变量 m_BXZoomRatio（横向放大倍数，初值为 1.25） ，m_BYZoomRatio（纵向放大倍数，初值为 1.25） ，m_SXZoomRatio（横向缩小系数，初值为 0.8） ，m_SYZoomRatio（纵向缩小系数，初值是 0.8）中。其程序流程图如下：其代码实现如下：/放大系数设置菜单项的响应函数void CMyDIPDoc:OnZoombig() / TODO: Add your comma

6、nd handler code hereSetClass setclass; /设置对话框if (setclass.DoModal() = IDOK)/将用户设置值存到文档类中的相应变量中m_BXZoomRatio = setclass.m_getx;m_BYZoomRatio = setclass.m_gety;/缩小系数设置菜单项的响应函数void CMyDIPDoc:OnZoomsmall() 菜单-设置放大倍数缩小倍数CMyDIPDoc:OnZoombig()CMyDIPDoc:OnZoomsmall()/ TODO: Add your command handler code her

7、eSetClass 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 是指向需处理的设备无关位图的指针，fXZoomRa

8、tio 和fYZoomRatio 分别为横向和纵向缩放系数，tag 为一标记变量（ tag=1 时执行最近邻插值操作，tag=0 时执行双线性插值操作） 。程序流程图如下：其实现代码如下：HGLOBAL WINAPI ZoomDIB(LPSTR lpDIB, float fXZoomRatio, float fYZoomRatio, BOOL tag)/ 源图像的宽度和高度LONG lWidth;LONG lHeight;/ 缩放后图像的宽度和高度LONG lNewWidth;LONG lNewHeight;最近邻插值双线性插值放大缩小放大缩小OnMenuitem32779()OnMenuit

9、em32778()OnMenuitem32780()OnMenuitem32781()ZoomDIB()/ 缩放后图像的宽度（lNewWidth，必须是 4 的倍数）LONG lNewLineBytes;/ 指向源图像的指针LPSTR lpDIBBits;/ 指向源象素的指针LPSTR lpSrc;/ 缩放后新 DIB 句柄HDIB hDIB;/ 指向缩放图像对应象素的指针LPSTR lpDst;/ 指向缩放图像的指针LPSTR lpNewDIB;LPSTR lpNewDIBBits;/ 指向 BITMAPINFO 结构的指针（ Win3.0）LPBITMAPINFOHEADER lpbmi;

10、/ 指向 BITMAPCOREINFO 结构的指针LPBITMAPCOREHEADER lpbmc;/ 循环变量（象素在新 DIB 中的坐标）LONG i;LONG j;/ 象素在源 DIB 中的坐标LONG i0;LONG j0;double x,y;/ 图像每行的字节数LONG lLineBytes;/ 找到源 DIB 图像象素起始位置lpDIBBits = :FindDIBBits(lpDIB);/ 获取图像的宽度lWidth = :DIBWidth(lpDIB);/ 计算图像每行的字节数lLineBytes = WIDTHBYTES(lWidth * 8);/ 获取图像的高度lHeig

11、ht = :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, lN

12、ewLineBytes * 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);/

13、获取指针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

14、;/ 针对图像每行进行操作for(i = 0; i = 0) / 复制象素*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 *

15、(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);/

16、右上角点的像素值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();/双线性插值缩小图片实验结果运行程序，各模块实验结果如下：读入灰度图片和设置放大缩小系数的结果读入图片：设置放大缩小系数：