数字图像处理

旋转(rotation)有一个绕着什么转的问题,通常的做法是以图象的中心为圆心旋转,举个例子,图2.7旋转30度(顺时针方向)后如图2.8所示:
\
可以看出,旋转后图象变大了。另一种做法是不让图象变大,转出的部分被裁剪掉。如图2.9所示。
\
我们采用第一种做法,首先给出变换矩阵。在我们熟悉的坐标系中,将一个点顺时针旋转a角后的坐标变换公式,如图2.10所示,r为该点到原点的距离,在旋转过程中,r保持不变;b为r与x轴之间的夹角。
旋转前:x0=rcosb;y0=rsinb
旋转a角度后:
x1=rcos(b-a)=rcosbcosa+rsinbsina= x0cosa+y0sina;
y1=rsin(b-a)=rsinbcosa-rcosbsina= -x0sina+y0cosa;
以矩阵的形式表示: \
(2.5)
上面的公式中,坐标系xoy是以图象的中心为原点,向右为x轴正方向,向上为y轴正方向。它和以图象左上角点为原点o’,向右为x’轴正方向,向下为y’轴正方向的坐标系x’o’y’之间的转换关系如何呢?如图2.11所示。
图2.11 两种坐标系间的转换关系
设图象的宽为w,高为h,容易得到:
(2.6)
逆变换为:
(2.7)
有了上面的公式,我们可以把变换分成三步:
1.将坐标系o’变成o;
2.将该点顺时针旋转a角;
3.将坐标系o变回o’,这样,我们就得到了变换矩阵,是上面三个矩阵的级联。
\
(2.8)
要注意的是,因为新图变大,所以上面公式中出现了wold,hold,wnew,hnew,它们分别表示原图(old)和新图(new)的宽、高。我们从图2.8中容易看出:wnew=max(|x4-x1|,|x3-x2|);hnew=max(|y4-y1|,|y3-y2|)。
(2.8)的逆变换为
\(2.9)
这样,对于新图中的每一点,我们就可以根据公式(2.9)求出对应原图中的点,得到它的灰度。如果超出原图范围,则填成白色。要注意的是,由于有浮点运算,计算出来点的坐标可能不是整数,采用取整处理,即找最接近的点,这样会带来一些误差(图象可能会出现锯齿)。更精确的方法是采用插值,将在图象缩放时介绍。
源程序如下:
[cpp]
/**
* 程序名: Rotation.cpp
* 功 能: 实现灰度图像的旋转,如果超出原图范围,则用白色填充
* 测试位图为test.bmp放到工程目录下
*/
#include <iostream>
#include <fstream>
#include <cstring>
#include <cmath>
#include < windows.h>
using namespace std;
#define PI 3.1415926535
#define RADIAN(angle) (((angle)*PI)/180.0)
BITMAPFILEHEADER bmpFileHeader; //bmp文件头
BITMAPINFOHEADER bmpInfoHeader; //bmp信息头
RGBQUAD *pColorTable; //bmp颜色表
unsigned char *pBmpData; //bmp位图数据
unsigned char *pNewBmpData; //旋转后bmp位图数据
int newImgSize; //旋转后图像大小
/**
* 函数名: readBmp
* 参 数: fileName--要读取文件的文件名
* 功 能: 读取bmp位图数据,成功返回TRUE,否则返回FALSE
*/
BOOL readBmp(char *fileName)
{
FILE *fp = fopen(fileName,"rb"); //以二进制读方式打开
if (NULL == fp)
{
cout<<"The file is opened failure!"<<endl;
return FALSE;
}
//读取信息到相应的变量中
fread(&bmpFileHeader,sizeof(BITMAPFILEHEADER),1,fp);
fread(&bmpInfoHeader,sizeof(BITMAPINFOHEADER),1,fp);
pColorTable = new RGBQUAD[256];
fread(pColorTable,sizeof(RGBQUAD),256,fp);
int imgSize = bmpInfoHeader.biSizeImage; //文图数据的大小,4的倍数
pBmpData = new unsigned char[imgSize];
fread(pBmpData,sizeof(unsigned char),imgSize,fp);
fclose(fp);
return TRUE;
}
/**
* 函数名: rotation
* 参 数: rotAngle--要转换的角度
* 功 能: 实现bmp图像的旋转
*/
void rotation(int rotAngle)
{
double cosa,sina,srcX[4],srcY[4],dstX[4],dstY[4],rad;
int oldWidth,oldHeight,newWidth,newHeight;
rad = (double)RADIAN(rotAngle);
cosa = cos(rad);
sina = sin(rad);
//原图宽与高
oldWidth = bmpInfoHeader.biWidth;
oldHeight = bmpInfoHeader.biHeight;
//原图四个角的坐标
srcX[0] = -0.5 * oldWidth;
srcY[0] = 0.5 * oldHeight;
srcX[1] = 0.5 * oldWidth;
srcY[1] = 0.5 * oldHeight;
srcX[2] = 0.5 * oldWidth;
srcY[2] = -0.5 * oldHeight;
srcX[3] = -0.5 * oldWidth;
srcY[3] = -0.5 * oldHeight;
//新图四个角坐标
for(int i = 0; i < 4; i++ )
{
dstX[i] = cosa * srcX[i] + sina * srcY[i];
dstY[i] = -sina * srcX[i] + cosa * srcY[i];
// cout<<dstY[i]<<endl;
}
//新图的宽与高,向上取整
bmpInfoHeader.biWidth = newWidth = (int)(max(fabs(dstX[0] - dstX[2]),fabs(dstX[1] - dstX[3])) + 0.5);
bmpInfoHeader.biHeight = newHeight = (int)(max(fabs(dstY[0] - dstY[2]),fabs(dstY[1] - dstY[3])) + 0.5);
// cout<<newWidth<<newHeight<<endl;
//新图位图数据大小
bmpInfoHeader.biSizeImage = newImgSize = newHeight * ((newWidth * bmpInfoHeader.biBitCount + 31) / 32 * 4);
pNewBmpData = new unsigned char[newImgSize];
double temp1,temp2; //计算矩阵(2.9)中的两个常数,这样不用以后每次都计算了
temp1 = -0.5 * newWidth * cosa - 0.5 * newHeight * sina + 0.5 * oldWidth;
temp2 = 0.5 * newWidth * sina - 0.5 * newHeight * cosa + 0.5 * oldHeight;
memset(pNewBmpData,(BYTE)255,newImgSize); //先全部填充成白色
int x0,y0,x1,y1;
unsigned char *pOldTemp,*pNewTemp;
int oldLineByte,newLineByte;
oldLineByte = (oldWidth * bmpInfoHeader.biBitCount + 31) / 32 * 4;
newLineByte = (newWidth * bmpInfoHeader.biBitCount + 31) / 32 * 4;
//把旋转后的图像数据对应存储到pNewBmpData相应位置
for(y1 = 0; y1 < newHeight; y1++)
{
for(x1 = 0; x1 < newWidth; x1++ )
{
x0 = (int)(x1 * cosa + y1 * sina + temp1);
y0 = (int)(-x1 * sina + y1 * cosa + temp2);
if((x0 >= 0 && x0 < oldWidth) && (y0 >= 0 && y0 < oldHeight)) //这里不能为<=oldWidth或oldHeight
{
pOldTemp = pBmpData + (oldHeight - 1 - y0) * oldLineByte + x0;
pNewTemp = pNewBmpData + (newHeight - 1 - y1) * newLineByte + x1;
*pNewTemp = *pOldTemp;
}
}
}
}
/**
* 函数名: writeBmp
* 参 数: bmpName -- 旋转后的bmp文件名
* 功 能: 新建一个bmp文件,把旋转后的图像信息存入其中
*/
void writeBmp(char *bmpName)
{
FILE *fp = fopen(bmpName,"wb"); //以二进制写方式打开
if(NULL == fp)
cout<<"The file is opened failure"<<endl;
//写入选装后图像信息
fwrite(&bmpFileHeader,sizeof(BITMAPFILEHEADER),1,fp);
fwrite(&bmpInfoHeader,sizeof(BITMAPINFOHEADER),1,fp);
fwrite(pColorTable,sizeof(RGBQUAD),256,fp);
fwrite(pNewBmpData,sizeof(unsigned char),newImgSize,fp);
fclose(fp);
delete []pColorTable;
delete []pNewBmpData;
delete []pBmpData;
}
/**
* 函数名: work
* 参 数: 无
* 功 能: 实现处理工作
*/
void work()
{
char readBmpName[] = "test.bmp";
if ( !readBmp(readBmpName))
cout<<"The file "<<readBmpName<<"is read failure"<<endl;
cout<<"please input the angle to rotate(Clockwise):";
int rotAngle;
cin>>rotAngle;
rotation(rotAngle);
char writeBmpName[] = "test_new.bmp";
writeBmp(writeBmpName);
}
int main()
{
work();
return 0;
}
/**
* 程序名: Rotation.cpp
* 功 能: 实现灰度图像的旋转,如果超出原图范围,则用白色填充
* 测试位图为test.bmp放到工程目录下
*/
#include <iostream>
#include <fstream>
#include <cstring>
#include <cmath>
#include <windows.h>
using namespace std;
#define PI 3.1415926535
#define RADIAN(angle) (((angle)*PI)/180.0)
BITMAPFILEHEADER bmpFileHeader; //bmp文件头
BITMAPINFOHEADER bmpInfoHeader; //bmp信息头
RGBQUAD *pColorTable; //bmp颜色表
unsigned char *pBmpData; //bmp位图数据
unsigned char *pNewBmpData; //旋转后bmp位图数据
int newImgSize; //旋转后图像大小
/**
* 函数名: readBmp
* 参 数: fileName--要读取文件的文件名
* 功 能: 读取bmp位图数据,成功返回TRUE,否则返回FALSE
*/
BOOL readBmp(char *fileName)
{
FILE *fp = fopen(fileName,"rb"); //以二进制读方式打开
if (NULL == fp)
{
cout<<"The file is opened failure!"<<endl;
return FALSE;
}
//读取信息到相应的变量中
fread(&bmpFileHeader,sizeof(BITMAPFILEHEADER),1,fp);
fread(&bmpInfoHeader,sizeof(BITMAPINFOHEADER),1,fp);
pColorTable = new RGBQUAD[256];
fread(pColorTable,sizeof(RGBQUAD),256,fp);
int imgSize = bmpInfoHeader.biSizeImage; //文图数据的大小,4的倍数
pBmpData = new unsigned char[imgSize];
fread(pBmpData,sizeof(unsigned char),imgSize,fp);
fclose(fp);
return TRUE;
}
/**
* 函数名: rotation
* 参 数: rotAngle--要转换的角度
* 功 能: 实现bmp图像的旋转
*/
void rotation(int rotAngle)
{
double cosa,sina,srcX[4],srcY[4],dstX[4],dstY[4],rad;
int oldWidth,oldHeight,newWidth,newHeight;
rad = (double)RADIAN(rotAngle);
cosa = cos(rad);
sina = sin(rad);
//原图宽与高
oldWidth = bmpInfoHeader.biWidth;
oldHeight = bmpInfoHeader.biHeight;
//原图四个角的坐标
srcX[0] = -0.5 * oldWidth;
srcY[0] = 0.5 * oldHeight;
srcX[1] = 0.5 * oldWidth;
srcY[1] = 0.5 * oldHeight;
srcX[2] = 0.5 * oldWidth;
srcY[2] = -0.5 * oldHeight;
srcX[3] = -0.5 * oldWidth;
srcY[3] = -0.5 * oldHeight;
//新图四个角坐标
for(int i = 0; i < 4; i++ )
{
dstX[i] = cosa * srcX[i] + sina * srcY[i];
dstY[i] = -sina * srcX[i] + cosa * srcY[i];
// cout<<dstY[i]<<endl;
}
//新图的宽与高,向上取整
bmpInfoHeader.biWidth = newWidth = (int)(max(fabs(dstX[0] - dstX[2]),fabs(dstX[1] - dstX[3])) + 0.5);
bmpInfoHeader.biHeight = newHeight = (int)(max(fabs(dstY[0] - dstY[2]),fabs(dstY[1] - dstY[3])) + 0.5);
// cout<<newWidth<<newHeight<<endl;
//新图位图数据大小
bmpInfoHeader.biSizeImage = newImgSize = newHeight * ((newWidth * bmpInfoHeader.biBitCount + 31) / 32 * 4);
pNewBmpData = new unsigned char[newImgSize];
double temp1,temp2; //计算矩阵(2.9)中的两个常数,这样不用以后每次都计算了
temp1 = -0.5 * newWidth * cosa - 0.5 * newHeight * sina + 0.5 * oldWidth;
temp2 = 0.5 * newWidth * sina - 0.5 * newHeight * cosa + 0.5 * oldHeight;
memset(pNewBmpData,(BYTE)255,newImgSize); //先全部填充成白色
int x0,y0,x1,y1;
unsigned char *pOldTemp,*pNewTemp;
int oldLineByte,newLineByte;
oldLineByte = (oldWidth * bmpInfoHeader.biBitCount + 31) / 32 * 4;
newLineByte = (newWidth * bmpInfoHeader.biBitCount + 31) / 32 * 4;
//把旋转后的图像数据对应存储到pNewBmpData相应位置
for(y1 = 0; y1 < newHeight; y1++)
{
for(x1 = 0; x1 < newWidth; x1++ )
{
x0 = (int)(x1 * cosa + y1 * sina + temp1);
y0 = (int)(-x1 * sina + y1 * cosa + temp2);
if((x0 >= 0 && x0 < oldWidth) && (y0 >= 0 && y0 < oldHeight)) //这里不能为<=oldWidth或oldHeight
{
pOldTemp = pBmpData + (oldHeight - 1 - y0) * oldLineByte + x0;
pNewTemp = pNewBmpData + (newHeight - 1 - y1) * newLineByte + x1;
*pNewTemp = *pOldTemp;
}
}
}
}
/**
* 函数名: writeBmp
* 参 数: bmpName -- 旋转后的bmp文件名
* 功 能: 新建一个bmp文件,把旋转后的图像信息存入其中
*/
void writeBmp(char *bmpName)
{
FILE *fp = fopen(bmpName,"wb"); //以二进制写方式打开
if(NULL == fp)
cout<<"The file is opened failure"<<endl;
//写入选装后图像信息
fwrite(&bmpFileHeader,sizeof(BITMAPFILEHEADER),1,fp);
fwrite(&bmpInfoHeader,sizeof(BITMAPINFOHEADER),1,fp);
fwrite(pColorTable,sizeof(RGBQUAD),256,fp);
fwrite(pNewBmpData,sizeof(unsigned char),newImgSize,fp);
fclose(fp);
delete []pColorTable;
delete []pNewBmpData;
delete []pBmpData;
}
/**
* 函数名: work
* 参 数: 无
* 功 能: 实现处理工作
*/
void work()
{
char readBmpName[] = "test.bmp";
if ( !readBmp(readBmpName))
cout<<"The file "<<readBmpName<<"is read failure"<<endl;
cout<<"please input the angle to rotate(Clockwise):";
int rotAngle;
cin>>rotAngle;
rotation(rotAngle);
char writeBmpName[] = "test_new.bmp";
writeBmp(writeBmpName);
}
int main()
{
work();
return 0;
}
  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值