图像的缩放和数学计算.docx

1、图像的缩放和数学计算医学图像处理实验二：图像的缩放和数学运算摘要 本次实验的主要内容是通过了解学习图像缩放的原理和方法和图像的数学运算。 对一张图片用不同的方法进行缩放，其中包括最近邻内插法，双线性内插法和双三次内插法三种。改变缩放的倍数，比较不同缩放方法的异同； 了解学习图像进行加减乘除四种数学运算的原理以及两个图片进行数学运算的条件。并对符合条件的两幅图进行数学运算处理，用运算结果帮助理解图片数学运算的原理。一、技术讨论1.1实验原理图像的缩放实验的原理： 图片是由像素矩阵表示的； 最近邻内插法（Pixel Replication）是把最邻近的灰度赋予每个新位置，利用opencv调用函数c

2、vResize(img1, img3, CV_INTER_NN)实现最近邻插值以及cvResize；重新调整图像； 双线性内插（Bilinear Interpolation）是用4个最近邻点去估计给定位置的灰度，令v(x,y)为原图像(x,y)位置的像素点的灰度值v(X,Y)=ax+by+cxy+d，其中a,b,c,d可由四个系数通过（x，y）最近四个邻点的方程得出，之后利用opencv调用函数cvResize(img1, img3, CV_INTER_LINEAR)实现双线性插值； 双三次内插（Bicubic interpolation）是用16个最近邻点去估计给定位置的灰度，然后用open

3、cv调用函数cvResize(img1, img3, CV_INTER_CUBIC)实现双三次内插。图像的数学运算实验的原理： 图片是由像素矩阵表示的； 图像的算数操作是像素阵列操作，具体如下： 做完运算的灰度值会超出0-255范围，这时再做一个映射映射到0-255。1.2实验方法图像的缩放实验主要运用了三种方法：最近邻内插法（Pixel Replication），双线性内插（Bilinear Interpolation）和双三次内插（Bicubic interpolation）。图像的数学运算实验只是简单的对符合数学运算大小条件的两张图片像素矩阵进行加减乘除，若做完运算的灰度值超出0-255

4、范围，则再做一个映射映射到0-255。没有特别的方法。二、结果与讨论图片的缩放实验一（此时scale=2） 1-1. Original image1-2.Pixel Replication result1-3.Bilinear Interpolation result1-4.Bicubic_Interpolation_result图片的缩放实验二（此时scale=5）2-1.Original image2-2.Pixel Replication result2-3.Bilinear Interpolation result2-4.Bicubic_Interpolation_result讨论：

5、最近邻内插：放大后的图像有很严重的马赛克，缩小后的图像有很严重的失真，现象由1-2,2-2可知； 对于双线性内插值法，由计算公式可知它计算量大，但缩放后图像质量高，不会出现像素值不连续的的情况。由于双线性插值具有低通滤波器的性质，使高频分量受损，所以图像轮廓在一定程度上变得模糊，由1-3,2-3可知； 双三次内插比计算量更大，因此计算精度也更高，计算新增加的像素的灰度值时，需要考虑周围的16个邻点，因此此方法能保留图片更好的细节质量。由实验结果1-4,2-4可对比另外两种实验结果可得； 在放大倍数增加的时候，三种缩放方法的差异更明显，这一点可通过对比scale分别为2和5时的对应图片信息可知。

6、注：在Qt软件上这些细节看的更清楚一些，在Word文档没那么清晰。图片的数学运算实验 3-1.original image1（详细信息：1200*2000） 3-2.original image2（详细信息：800*533） 3-3.加法结果 3-4.减法结果 3-5.乘法结果 3-6.除法结果讨论： 加法运算主要实现图像的叠加性，3-3明显比原图，即3-2色调暗，就是和3-1进行加法运算的结果，但由于3-1和3-2图像大小不一样，显示效果不太明显； 减法运算实现图像相减附录（分页写）项目一程序：.pro目录下的程序如下：QT += coreQT -= guiTARGET = Project-

7、ArithmeticOperationsCONFIG += consoleCONFIG -= app_bundleTEMPLATE = appSOURCES += main.cppINCLUDEPATH+=D:Qtopencv2.2includeopencvD:Qtopencv2.2includeopencv2D:Qtopencv2.2includeLIBS+=D:Qtopencv2.2liblibopencv_calib3d220.dll.aD:Qtopencv2.2liblibopencv_contrib220.dll.aD:Qtopencv2.2liblibopencv_core220.

8、dll.aD:Qtopencv2.2liblibopencv_features2d220.dll.aD:Qtopencv2.2liblibopencv_flann220.dll.aD:Qtopencv2.2liblibopencv_gpu220.dll.aD:Qtopencv2.2liblibopencv_highgui220.dll.aD:Qtopencv2.2liblibopencv_imgproc220.dll.aD:Qtopencv2.2liblibopencv_legacy220.dll.aD:Qtopencv2.2liblibopencv_ml220.dll.aD:Qtopencv

9、2.2liblibopencv_objdetect220.dll.aD:Qtopencv2.2liblibopencv_video220.dll.a.cpp目录下的程序如下：#include #include #include opencv2/imgproc/imgproc.hpp#include opencv2/imgproc/imgproc_c.hint main() /载入图片 IplImage* img0 = cvLoadImage(d:/1.jpg); IplImage* img1 = cvCreateImage( cvSize(img0-width, img0-height), 8

10、, 1 );cvZero(img1); IplImage* img2 = cvCreateImage( cvSize(3*img0-width, 3*img0-height), 8, 1 );cvZero(img2); /转化为单通道 cvCvtColor(img0, img1, CV_BGR2GRAY ); cvNamedWindow(MyImage); cvShowImage(MyImage,img1); /缩放 float Scale =2; CvSize dst_cvsize; /目标图像尺寸 dst_cvsize.width = img1-width * Scale; /目标图像的宽

11、为源图象宽的scale倍 dst_cvsize.height = img1-height * Scale; /目标图像的高为源图象高的scale倍 IplImage* img3 = cvCreateImage( dst_cvsize, 8, 1 );cvZero(img3); IplImage* img4 = cvCreateImage( dst_cvsize, 8, 1 );cvZero(img4); IplImage* img5 = cvCreateImage( dst_cvsize, 8, 1 );cvZero(img5); cvResize(img1, img3, CV_INTER_N

12、N);/最近邻插值 cvSaveImage(Pixel Replication result.jpg,img3); cvNamedWindow(Pixel Replication result); cvShowImage(Pixel Replication result,img3); cvResize(img1, img4, CV_INTER_LINEAR);/双线性插值 cvSaveImage(Bilinear Interpolation result.jpg,img4); cvNamedWindow(Bilinear Interpolation result); cvShowImage(B

13、ilinear Interpolation result,img4); cvResize(img1, img5, CV_INTER_CUBIC); /双三次内插 cvSaveImage(Bicubic_Interpolation_result.jpg,img5); cvNamedWindow(Bicubic_Interpolation_result); cvShowImage(Bicubic_Interpolation_result,img5); /等待用户相应 cvWaitKey(0); /释放图像 cvReleaseImage(&img0); cvReleaseImage(&img1);

14、cvReleaseImage(&img2); cvReleaseImage(&img3); cvReleaseImage(&img4); cvReleaseImage(&img5); /毁掉图像显示窗口 cvDestroyWindow(MyImage); cvDestroyWindow(Pixel Replication result); cvDestroyWindow(Bilinear Interpolation result); cvDestroyWindow(Bicubic_Interpolation_result); return 0;程序流程图如下：项目二程序：.pro目录下的程序如下：QT += coreQT -= guiTARGET = Project-ArithmeticOperationsCONFIG += consoleCONFIG -= app_bundleTEMPLATE = appSOURCES += main.cppINCLUDEPATH+=D:Qtopencv2.2includeopencvD:Qtopencv2.2includeopencv2D:Qtopencv2.2includeLIBS+=D:Qtopencv2.2liblibopencv_calib3d220.dll.aD:Qtopencv2.2liblib