程序博客网 > 如何自学成为程序员

C/C++ BMP（24位真彩色）图像处理（4）------图像の旋转

来源：互联网 发布：如何自学成为程序员 编辑：程序博客网 时间：2024/06/05 18:57

    历经一个多月，CSDN貌似终于好像把文章列表阅读量信息归零BUG给修好了，于是乎放篇做期末大作业时写的文章上来测测效果，可别又像上次一样一发文章就又坑爹了啊！

    本篇谈的是图像的旋转，不算是什么新鲜的题目了。但是现在由于很多工具如MATLAB、OpenCV等都把算法写好给用户调用，导致大多用户只知其然不知其所以然，所以回顾一下也是好的。

    图像的旋转，说到底就是每个像素点绕着某个圆心旋转一定角度。如果是写代码的话，旋转的角度和圆心应该是已知的条件，我们第一个思路是根据已知条件求取出图像经过旋转后的新的像素点的坐标。这个思路显然是没有错的，然而有个问题就是旋转之后的图像有些像素点已经被转出到图像的大小之外，而另外有一些像素点无法跟原图像相对应。这就使得顺着这样一个思路编程变得复杂。因而我们反过来想，不如先将旋转后的图像规划好，然后通过一个方法去计算新图像每个像素点是否对应于原图像的像素点，如果对应那么该像素点处于原图像的什么位置，如果不对应那该像素点又该取值多少。通过这个方法我们遍历新图像的所有像素点则可得到旋转后的图像。

    OK，上面我们已经将图像旋转的思路理清了。现在我们还不知道的就是图像旋转前后的像素点对应关系。这个关系稍微复杂一点，本人参考了这个博客的计算方法，下面进行推导或者说重现。


    请看上图，经过上面的推导，我们是要从x’，y’去得到x，y的。这也就是说我们最终要得到的形式是：

    x=f(x’,y’)

    y=g(x’,y’)

    这样子的。OK,计算公式推导如下（用word编写的，风格略不同。。。）：

           

    OK，理论工作做完了，接下来就是编程了，需要注意的是由于旋转后图像的点对应于旋转前图像的点的位置基本上不大可能是正对着的，所以需要用到双线性插值（该理论的具体分析在上一篇文章详细提到）。

    主要代码如下：

[cpp] view plain copy
print?
  1. /*******************图像旋转******************/  
  2.     int centerofrotation_x = MYDRAW_WIDTH / 2;//旋转中心x坐标  
  3.     int centerofrotation_y = MYDRAW_HEIGHT / 2;//旋转中心y坐标  
  4.     double degree = 60;//逆时针旋转角度  
  5.     double radian = degree / 180 * pi;//角度转弧度  
  6.     for (int hnum = 0; hnum < MYDRAW_HEIGHT; hnum++)  
  7.     for (int wnum = 0; wnum < MYDRAW_WIDTH; wnum++)  
  8.     {  
  9.         int pixel_point = hnum*write_width + wnum * 3;//映射尺度变换图像数组位置偏移量  
  10.   
  11.         double d_original_img_wnum = (wnum - centerofrotation_x)*cos(radian) - (hnum - centerofrotation_y)*sin(radian) + centerofrotation_x;  
  12.         double d_original_img_hnum = (wnum - centerofrotation_x)*sin(radian) + (hnum - centerofrotation_y)*cos(radian) + centerofrotation_y;  
  13.   
  14.         if (d_original_img_wnum<0 || d_original_img_wnum>width || d_original_img_hnum<0 || d_original_img_hnum>height)//找不到与原图像的对于关系  
  15.         {  
  16.             pColorDataMid[pixel_point] = pColorDataMid[pixel_point + 1] = pColorDataMid[pixel_point + 2] = 0;  
  17.             continue;  
  18.         }  
  19.   
  20.         int i_original_img_hnum = d_original_img_hnum;  
  21.         int i_original_img_wnum = d_original_img_wnum;  
  22.         double distance_to_a_x = d_original_img_wnum - i_original_img_wnum;//在原图像中与a点的水平距离  
  23.         double distance_to_a_y = d_original_img_hnum - i_original_img_hnum;//在原图像中与a点的垂直距离  
  24.   
  25.         int original_point_a = i_original_img_hnum*l_width + i_original_img_wnum * 3;//数组位置偏移量，对应于图像的各像素点RGB的起点,相当于点A    
  26.         int original_point_b = i_original_img_hnum*l_width + (i_original_img_wnum + 1) * 3;//数组位置偏移量，对应于图像的各像素点RGB的起点,相当于点B  
  27.         int original_point_c = (i_original_img_hnum + 1)*l_width + i_original_img_wnum * 3;//数组位置偏移量，对应于图像的各像素点RGB的起点,相当于点C   
  28.         int original_point_d = (i_original_img_hnum + 1)*l_width + (i_original_img_wnum + 1) * 3;//数组位置偏移量，对应于图像的各像素点RGB的起点,相当于点D   
  29.         if (i_original_img_hnum == MYDRAW_HEIGHT - 1)  
  30.         {  
  31.             original_point_c = original_point_a;  
  32.             original_point_d = original_point_b;  
  33.         }  
  34.         if (i_original_img_wnum == MYDRAW_WIDTH - 1)  
  35.         {  
  36.             original_point_a = original_point_b;  
  37.             original_point_c = original_point_d;  
  38.         }  
  39.   
  40.         pColorDataMid[pixel_point] =  
  41.             pColorData[original_point_a] * (1 - distance_to_a_x)*(1 - distance_to_a_y) +  
  42.             pColorData[original_point_b] * distance_to_a_x*(1 - distance_to_a_y) +  
  43.             pColorData[original_point_c] * distance_to_a_y*(1 - distance_to_a_x) +  
  44.             pColorData[original_point_c] * distance_to_a_y*distance_to_a_x;  
  45.         pColorDataMid[pixel_point + 1] =  
  46.             pColorData[original_point_a + 1] * (1 - distance_to_a_x)*(1 - distance_to_a_y) +  
  47.             pColorData[original_point_b + 1] * distance_to_a_x*(1 - distance_to_a_y) +  
  48.             pColorData[original_point_c + 1] * distance_to_a_y*(1 - distance_to_a_x) +  
  49.             pColorData[original_point_c + 1] * distance_to_a_y*distance_to_a_x;  
  50.         pColorDataMid[pixel_point + 2] =  
  51.             pColorData[original_point_a + 2] * (1 - distance_to_a_x)*(1 - distance_to_a_y) +  
  52.             pColorData[original_point_b + 2] * distance_to_a_x*(1 - distance_to_a_y) +  
  53.             pColorData[original_point_c + 2] * distance_to_a_y*(1 - distance_to_a_x) +  
  54.             pColorData[original_point_c + 2] * distance_to_a_y*distance_to_a_x;  
  55.   
  56.     }  
  57.     /*******************图像旋转******************/  
    整个工程代码如下：
[cpp] view plain copy
print?
  1. #include <string.h>   
  2. #include <math.h>     
  3. #include <stdio.h>     
  4. #include <stdlib.h>     
  5. #include <malloc.h>  
  6.   
  7. #include<time.h>//时间相关头文件，可用其中函数计算图像处理速度  
  8.   
  9. #define   WIDTHBYTES(bits) (((bits)+31)/32*4)//用于使图像宽度所占字节数为4byte的倍数  
  10.   
  11. #define MYDRAW_HEIGHT 1054  //目标图像高度  
  12. #define MYDRAW_WIDTH 1500  //目标图像宽度  
  13.   
  14. #define pi 3.1415926535  
  15.   
  16. typedef unsigned char  BYTE;  
  17. typedef unsigned short WORD;  
  18. typedef unsigned long  DWORD;  
  19. typedef long LONG;  
  20.   
  21. //位图文件头信息结构定义  
  22. //其中不包含文件类型信息（由于结构体的内存结构决定，要是加了的话将不能正确读取文件信息）  
  23.   
  24. typedef struct tagBITMAPFILEHEADER {  
  25.     DWORD  bfSize;          //文件大小  
  26.     WORD   bfReserved1;     //保留字，不考虑  
  27.     WORD   bfReserved2;     //保留字，同上  
  28.     DWORD  bfOffBits;       //实际位图数据的偏移字节数，即前三个部分长度之和  
  29. } BITMAPFILEHEADER;  
  30.   
  31. //信息头BITMAPINFOHEADER，也是一个结构，其定义如下：  
  32.   
  33. typedef struct tagBITMAPINFOHEADER{  
  34.     //public:  
  35.     DWORD   biSize;             //指定此结构体的长度，为40  
  36.     LONG    biWidth;            //位图宽  
  37.     LONG    biHeight;           //位图高  
  38.     WORD    biPlanes;           //平面数，为1  
  39.     WORD    biBitCount;         //采用颜色位数，可以是1，2，4，8，16，24，新的可以是32  
  40.     DWORD   biCompression;      //压缩方式，可以是0，1，2，其中0表示不压缩  
  41.     DWORD   biSizeImage;        //实际位图数据占用的字节数  
  42.     LONG    biXPelsPerMeter;    //X方向分辨率  
  43.     LONG    biYPelsPerMeter;    //Y方向分辨率  
  44.     DWORD   biClrUsed;          //使用的颜色数，如果为0，则表示默认值(2^颜色位数)  
  45.     DWORD   biClrImportant;     //重要颜色数，如果为0，则表示所有颜色都是重要的  
  46. } BITMAPINFOHEADER;  
  47.   
  48. void main()  
  49. {  
  50.     long now = 0;  
  51.     now = clock();//存储图像处理开始时间  
  52.   
  53.     BITMAPFILEHEADER bitHead, writebitHead;  
  54.     BITMAPINFOHEADER bitInfoHead, writebitInfoHead;  
  55.     FILE* pfile;//输入文件  
  56.     FILE* wfile;//输出文件  
  57.   
  58.     char strFile[50] = "15.bmp";//打开图像路径，BMP图像必须为24位真彩色格式  
  59.     char strFilesave[50] = "16.bmp";//处理后图像存储路径  
  60.     fopen_s(&pfile, strFile, "rb");//文件打开图像  
  61.     fopen_s(&wfile, strFilesave, "wb");//打开文件为存储修改后图像做准备  
  62.     //读取位图文件头信息  
  63.     WORD fileType;  
  64.     fread(&fileType, 1, sizeof(WORD), pfile);  
  65.     fwrite(&fileType, 1, sizeof(WORD), wfile);  
  66.     if (fileType != 0x4d42)  
  67.     {  
  68.         printf("file is not .bmp file!");  
  69.         return;  
  70.     }  
  71.     //读取位图文件头信息  
  72.     fread(&bitHead, 1, sizeof(tagBITMAPFILEHEADER), pfile);  
  73.     writebitHead = bitHead;//由于截取图像头和源文件头相似，所以先将源文件头数据赋予截取文件头  
  74.     //读取位图信息头信息  
  75.     fread(&bitInfoHead, 1, sizeof(BITMAPINFOHEADER), pfile);  
  76.     writebitInfoHead = bitInfoHead;//同位图文件头相似  
  77.   
  78.     writebitInfoHead.biHeight = MYDRAW_HEIGHT;//为截取文件重写位图高度  
  79.     writebitInfoHead.biWidth = MYDRAW_WIDTH;//为截取文件重写位图宽度  
  80.     int mywritewidth = WIDTHBYTES(writebitInfoHead.biWidth*writebitInfoHead.biBitCount);//BMP图像实际位图数据区的宽度为4byte的倍数，在此计算实际数据区宽度  
  81.     writebitInfoHead.biSizeImage = mywritewidth*writebitInfoHead.biHeight;//计算位图实际数据区大小  
  82.   
  83.     writebitHead.bfSize = 54 + writebitInfoHead.biSizeImage;//位图文件头大小为位图数据区大小加上54byte  
  84.     fwrite(&writebitHead, 1, sizeof(tagBITMAPFILEHEADER), wfile);//写回位图文件头信息到输出文件  
  85.     fwrite(&writebitInfoHead, 1, sizeof(BITMAPINFOHEADER), wfile);//写回位图信息头信息到输出文件  
  86.   
  87.     int width = bitInfoHead.biWidth;  
  88.     int height = bitInfoHead.biHeight;  
  89.     //分配内存空间把源图存入内存     
  90.     int l_width = WIDTHBYTES(width*bitInfoHead.biBitCount);//计算位图的实际宽度并确保它为4byte的倍数  
  91.     int write_width = WIDTHBYTES(writebitInfoHead.biWidth*writebitInfoHead.biBitCount);//计算写位图的实际宽度并确保它为4byte的倍数  
  92.   
  93.     BYTE    *pColorData = (BYTE *)malloc(height*l_width);//开辟内存空间存储图像数据  
  94.     memset(pColorData, 0, height*l_width);  
  95.   
  96.     BYTE    *pColorDataMid = (BYTE *)malloc(mywritewidth*MYDRAW_HEIGHT);//开辟内存空间存储图像处理之后数据  
  97.     memset(pColorDataMid, 0, mywritewidth*MYDRAW_HEIGHT);  
  98.   
  99.     long nData = height*l_width;  
  100.     long write_nData = mywritewidth*MYDRAW_HEIGHT;//截取的位图数据区长度定义  
  101.   
  102.     //把位图数据信息读到数组里     
  103.     fread(pColorData, 1, nData, pfile);//图像处理可通过操作这部分数据加以实现  
  104.   
  105.     /*******************图像处理部分******************/  
  106.     /*******************图像旋转******************/  
  107.     int centerofrotation_x = MYDRAW_WIDTH / 2;//旋转中心x坐标  
  108.     int centerofrotation_y = MYDRAW_HEIGHT / 2;//旋转中心y坐标  
  109.     double degree = 60;//逆时针旋转角度  
  110.     double radian = degree / 180 * pi;//角度转弧度  
  111.     for (int hnum = 0; hnum < MYDRAW_HEIGHT; hnum++)  
  112.     for (int wnum = 0; wnum < MYDRAW_WIDTH; wnum++)  
  113.     {  
  114.         int pixel_point = hnum*write_width + wnum * 3;//映射尺度变换图像数组位置偏移量  
  115.   
  116.         double d_original_img_wnum = (wnum - centerofrotation_x)*cos(radian) - (hnum - centerofrotation_y)*sin(radian) + centerofrotation_x;  
  117.         double d_original_img_hnum = (wnum - centerofrotation_x)*sin(radian) + (hnum - centerofrotation_y)*cos(radian) + centerofrotation_y;  
  118.   
  119.         if (d_original_img_wnum<0 || d_original_img_wnum>width || d_original_img_hnum<0 || d_original_img_hnum>height)//找不到与原图像的对于关系  
  120.         {  
  121.             pColorDataMid[pixel_point] = pColorDataMid[pixel_point + 1] = pColorDataMid[pixel_point + 2] = 0;  
  122.             continue;  
  123.         }  
  124.   
  125.         int i_original_img_hnum = d_original_img_hnum;  
  126.         int i_original_img_wnum = d_original_img_wnum;  
  127.         double distance_to_a_x = d_original_img_wnum - i_original_img_wnum;//在原图像中与a点的水平距离  
  128.         double distance_to_a_y = d_original_img_hnum - i_original_img_hnum;//在原图像中与a点的垂直距离  
  129.   
  130.         int original_point_a = i_original_img_hnum*l_width + i_original_img_wnum * 3;//数组位置偏移量，对应于图像的各像素点RGB的起点,相当于点A    
  131.         int original_point_b = i_original_img_hnum*l_width + (i_original_img_wnum + 1) * 3;//数组位置偏移量，对应于图像的各像素点RGB的起点,相当于点B  
  132.         int original_point_c = (i_original_img_hnum + 1)*l_width + i_original_img_wnum * 3;//数组位置偏移量，对应于图像的各像素点RGB的起点,相当于点C   
  133.         int original_point_d = (i_original_img_hnum + 1)*l_width + (i_original_img_wnum + 1) * 3;//数组位置偏移量，对应于图像的各像素点RGB的起点,相当于点D   
  134.         if (i_original_img_hnum == MYDRAW_HEIGHT - 1)  
  135.         {  
  136.             original_point_c = original_point_a;  
  137.             original_point_d = original_point_b;  
  138.         }  
  139.         if (i_original_img_wnum == MYDRAW_WIDTH - 1)  
  140.         {  
  141.             original_point_a = original_point_b;  
  142.             original_point_c = original_point_d;  
  143.         }  
  144.   
  145.         pColorDataMid[pixel_point] =  
  146.             pColorData[original_point_a] * (1 - distance_to_a_x)*(1 - distance_to_a_y) +  
  147.             pColorData[original_point_b] * distance_to_a_x*(1 - distance_to_a_y) +  
  148.             pColorData[original_point_c] * distance_to_a_y*(1 - distance_to_a_x) +  
  149.             pColorData[original_point_c] * distance_to_a_y*distance_to_a_x;  
  150.         pColorDataMid[pixel_point + 1] =  
  151.             pColorData[original_point_a + 1] * (1 - distance_to_a_x)*(1 - distance_to_a_y) +  
  152.             pColorData[original_point_b + 1] * distance_to_a_x*(1 - distance_to_a_y) +  
  153.             pColorData[original_point_c + 1] * distance_to_a_y*(1 - distance_to_a_x) +  
  154.             pColorData[original_point_c + 1] * distance_to_a_y*distance_to_a_x;  
  155.         pColorDataMid[pixel_point + 2] =  
  156.             pColorData[original_point_a + 2] * (1 - distance_to_a_x)*(1 - distance_to_a_y) +  
  157.             pColorData[original_point_b + 2] * distance_to_a_x*(1 - distance_to_a_y) +  
  158.             pColorData[original_point_c + 2] * distance_to_a_y*(1 - distance_to_a_x) +  
  159.             pColorData[original_point_c + 2] * distance_to_a_y*distance_to_a_x;  
  160.   
  161.     }  
  162.     /*******************图像旋转******************/  
  163.     /*******************图像处理部分******************/  
  164.   
  165.     fwrite(pColorDataMid, 1, write_nData, wfile);   //将处理完图像数据区写回文件  
  166.     fclose(pfile);  
  167.     fclose(wfile);  
  168.   
  169.     printf("图像处理完成\n");  
  170.     printf("运行时间为：%dms\n", int(((double)(clock() - now)) / CLOCKS_PER_SEC * 1000));//输出图像处理花费时间信息  
  171. }  
    原图像以及得到的旋转图像如下：


原图像


旋转后图像





阅读全文
0 0
如何自学成为程序员 如何自学成为程序员
原创粉丝点击
热门IT博客
网络听课神器网络听课神器
从事数据行业从事数据行业
数据图
linux如何编辑文档
iphone6s 4g网络用不了
java短信发送例子
潞河国际教育学园 知乎
ajax遍历json数组对象
徐老师淘宝店铺
kool香烟淘宝上有卖吗
mac 常用软件
苹果mac怎么更新系统
时间序列数据的例子
移动4g是什么网络模式
淘宝网怎么开网店
面部识别软件动物
java编程思想pdf 目录
php开发用什么软件
js signature生成
怎么加淘宝刷流量群
热门问题 老师的惩罚 人脸识别 我在镇武司摸鱼那些年 重生之率土为王 我在大康的咸鱼生活 盘龙之生命进化 天生仙种 凡人之先天五行 春回大明朝 姑娘不必设防，我是瞎子 御宅书屋总榜 欲宅屋总点击榜 欲宅屋总收藏榜 都市之活了九世 不想上天榜 榜样三 学习榜样体会 2019榜样观后感或心得体会 榜样心得体会300字 榜样4的观后感 榜样4心得体会300字 榜样4什么时候播出 榜样2观后感心得体会 2018勤廉榜样心得体会2篇文库 榜组词 榜的组词是什么 榜字组词 榜怎么组词 榜的组词 榝 桑蚕丝与榨蚕丝的区别 榨干 榨干男配h 快穿榨干男配h 快穿榨干 小妖精榨干我 榨干我 破壁机榨果汁要加水吗 榨山茶油机 你所谓的爱是榨干我的理由 恩施榨广椒 榨广椒的做法 炸广椒的做法 榨水果汁的食谱大全 怎么榨水果汁 水果汁怎么榨 现榨橙汁 榨果汁食谱 九阳豆浆机榨果汁 自己榨果汁怎么搭配 现榨果汁能减肥吗 榨果汁的做法

程序博客网，程序员的互联网技术博客家园。csdn论坛精品 msdn技术资料都在这里