DIB（设备无关位图）旋转任意角度算法（单色位图）

      网上很多位图旋转的程序，但是一般都是8位、24位、32位位图的旋转，这些大于8位的位图每个像素都可以用整个字节表示，所以用char数组很容易实现对应像素复制。但是要对单色位图进行旋转的话，就涉及到按位复制，因为每个像素是用一个字节中的某一位表示的。

      我自己写了一个单色位图旋转的算法：

#include "math.h"
#define PI 3.14159
//角度到弧度转化的宏
#define RADIAN(angle) ((angle)*PI/180.0) 
HGLOBAL WINAPI RotateDIB(LPSTR lpDIB, int iRotateAngle)//单色位图旋转
{	
	LONG	lWidth,lHeight;// 源图像的宽度和高度	
	LONG	lNewWidth,lNewHeight;// 旋转后图像的宽度和高度

	LONG	lLineBytes;// 图像每行的字节数
	LONG	lNewLineBytes;//旋转后图像的宽度（lNewWidth'，必须是4的倍数）	
	
	LPSTR	lpDIBBits; //指向源图像的指针		
	LPSTR	lpNewDIBBits;
	
	LPSTR	lpSrc;    // 指向源象素的指针	
	LPSTR	lpDst;// 指向旋转图像对应象素的指针	
	
	HDIB	hDIB;// 旋转后新DIB句柄	
	LPSTR	lpNewDIB;// 指向旋转图像的指针(源图像指针是参数lpDIB)

	LPBITMAPINFOHEADER lpbmi;// 指向BITMAPINFO结构的指针（Win3.0）
	LPBITMAPCOREHEADER lpbmc;// 指向BITMAPCOREINFO结构的指针
		
	LONG	i,j;// 循环变量（象素在新DIB中的坐标）
	LONG	i0,j0;// 象素在源DIB中的坐标
	
	float	fRotateAngle;// 旋转角度（弧度）
	float	fSina, fCosa;// 旋转角度的正弦和余弦
	
	// 源图四个角的坐标（以图像中心为坐标系原点）
	float	fSrcX1,fSrcY1,fSrcX2,fSrcY2,fSrcX3,fSrcY3,fSrcX4,fSrcY4;	
	float	fDstX1,fDstY1,fDstX2,fDstY2,fDstX3,fDstY3,fDstX4,fDstY4;// 旋转后四个角的坐标（以图像中心为坐标系原点）
	
	// 两个中间常量
	float	f1,f2;
	
	lpDIBBits = ::FindDIBBits(lpDIB);// 找到源DIB图像象素起始位置		
	lWidth = ::DIBWidth(lpDIB);  // 获取图像的"宽度"（4的倍数）		
	lHeight = ::DIBHeight(lpDIB);// 获取图像的高度	
	lLineBytes = WIDTHBYTES(lWidth);//计算图像每行的字节数(单色位图)
	
	fRotateAngle = (float) RADIAN(iRotateAngle);// 将旋转角度从度转换到弧度		
	fSina = (float) sin((double)fRotateAngle);// 计算旋转角度的正弦
	fCosa = (float) cos((double)fRotateAngle);// 计算旋转角度的余弦
	
	// 计算原图的四个角的坐标（以图像中心为坐标系原点）
	fSrcX1 = (float) (- (lWidth  - 1) / 2);
	fSrcY1 = (float) (  (lHeight - 1) / 2);
	fSrcX2 = (float) (  (lWidth  - 1) / 2);
	fSrcY2 = (float) (  (lHeight - 1) / 2);
	fSrcX3 = (float) (- (lWidth  - 1) / 2);
	fSrcY3 = (float) (- (lHeight - 1) / 2);
	fSrcX4 = (float) (  (lWidth  - 1) / 2);
	fSrcY4 = (float) (- (lHeight - 1) / 2);
	
	// 计算新图四个角的坐标（以图像中心为坐标系原点）
	fDstX1 =  fCosa * fSrcX1 + fSina * fSrcY1;
	fDstY1 = -fSina * fSrcX1 + fCosa * fSrcY1;
	fDstX2 =  fCosa * fSrcX2 + fSina * fSrcY2;
	fDstY2 = -fSina * fSrcX2 + fCosa * fSrcY2;
	fDstX3 =  fCosa * fSrcX3 + fSina * fSrcY3;
	fDstY3 = -fSina * fSrcX3 + fCosa * fSrcY3;
	fDstX4 =  fCosa * fSrcX4 + fSina * fSrcY4;
	fDstY4 = -fSina * fSrcX4 + fCosa * fSrcY4;	
	
	lNewWidth  = (LONG) ( max( fabs(fDstX4 - fDstX1), fabs(fDstX3 - fDstX2) ) + 0.5);// 计算旋转后的图像实际宽度		
	lNewHeight = (LONG) ( max( fabs(fDstY4 - fDstY1), fabs(fDstY3 - fDstY2) )  + 0.5);// 计算旋转后的图像高度
	lNewLineBytes = WIDTHBYTES(lNewWidth);// 计算新图像每行的字节数(单色图)
	
	// 两个常数，这样不用以后每次都计算了
	f1 = (float) (-0.5 * (lNewWidth - 1) * fCosa - 0.5 * (lNewHeight - 1) * fSina + 0.5 * (lWidth  - 1));
	f2 = (float) ( 0.5 * (lNewWidth - 1) * fSina - 0.5 * (lNewHeight - 1) * fCosa + 0.5 * (lHeight - 1));
	
	
	// 分配内存，以保存新DIB
	hDIB = (HDIB) ::GlobalAlloc(GHND, lNewLineBytes * lNewHeight + *(LPDWORD)lpDIB + ::PaletteSize(lpDIB));	
	if (hDIB == NULL)// 判断是否内存分配失败
	{		
		return NULL;// 分配内存失败
	}
	// 锁定内存
	lpNewDIB =  (char * )::GlobalLock((HGLOBAL) hDIB);	
	memcpy(lpNewDIB, lpDIB, *(LPDWORD)lpDIB + ::PaletteSize(lpDIB));// 复制DIB信息头和调色板
	
	lpNewDIBBits = ::FindDIBBits(lpNewDIB);// 找到新DIB象素起始位置
	lpbmi = (LPBITMAPINFOHEADER)lpNewDIB;// 获取指针
	lpbmc = (LPBITMAPCOREHEADER)lpNewDIB;
	// 更新DIB中图像的高度和宽度
	if (IS_WIN30_DIB(lpNewDIB))// 对于Windows 3.0 DIB
	{	
		lpbmi->biWidth = lNewWidth;
		lpbmi->biHeight = lNewHeight;
	}
	else// 对于其它格式的DIB
	{		
		lpbmc->bcWidth = (unsigned short) lNewWidth;
		lpbmc->bcHeight = (unsigned short) lNewHeight;
	}

	for(i = 0; i < lNewHeight; i++)// 针对图像每行进行操作
	{	
		// 注意此处宽度和高度是新DIB的宽度和高度
		lpDst = (char *)lpNewDIBBits + lNewLineBytes * (lNewHeight - 1 - i);//第i行的数据（DIB是从下向上存储图像的）
		for(j = 0; j < lNewWidth; j++)// 针对图像每列进行操作
		{		 
		    int NewBytes = j/8;//一行中第几个字节
			int NewBits = j%8;//确定是在该字节中的第几位
			BYTE byNewByte = lpDst[NewBytes];//取一个字节，即8个象素	
			
			// 计算该象素在源DIB中的坐标
			i0 = (LONG) (-((float) j) * fSina + ((float) i) * fCosa + f2 + 0.5);//旋转前的目标行
			j0 = (LONG) ( ((float) j) * fCosa + ((float) i) * fSina + f1 + 0.5);//旋转前的目标列		
			// 判断是否在源图范围内
			if( (j0 >= 0) && (j0 < lWidth) && (i0 >= 0) && (i0 < lHeight))
			{				
				lpSrc = (char *)lpDIBBits + lLineBytes * (lHeight - 1 - i0);//指向源图像中第i0行
				int SrcBytes = j0/8;//该像素在这一行的第几个字节
				int SrcBits  = j0%8;//第几位
				BYTE bySrcByte = lpSrc[SrcBytes];//对应源图中的字节
				BYTE byTmp = (bySrcByte>>(7-SrcBits))&0x01;//取出源图像中的位,7-0		
			
				//注意：最高位代表最左边的象素!
				switch(NewBits)
				{
				case 0:
					byNewByte = (byTmp<<7)|byNewByte; //赋值给新DIB中对应的位
					break;
				case 1:
					byNewByte = (byTmp<<6)|byNewByte; 
					break;
				case 2:
					byNewByte = (byTmp<<5)|byNewByte; 
					break;
				case 3:
					byNewByte = (byTmp<<4)|byNewByte; 
					break;
				case 4:
					byNewByte = (byTmp<<3)|byNewByte; 
					break;
				case 5:
					byNewByte = (byTmp<<2)|byNewByte; 
					break;
				case 6:
					byNewByte = (byTmp<<1)|byNewByte; 
					break;
				case 7:
					byNewByte = byTmp|byNewByte; 
					break;
				}

				lpDst[NewBytes] = byNewByte;		
			}
			else
			{		
				lpDst[NewBytes] = 255;// 对于源图中没有的象素，直接赋值为255
			}
			
		}//列j结束	
	}//行i结束
	return hDIB;// 返回
}

具体调用可以参考我上传的资源 ：  VS2010实现单色位图旋转

上面这种算法边缘锯齿非常严重，下面是边缘处理稍微好一些的算法：

memset( lpNewDIBBits, 0xff, lNewLineBytes*lNewHeight );//先将每个像素设为白色
for(i = 0; i < lNewHeight; i++)// 针对图像每行进行操作
{	
	for(j = 0; j < lNewWidth; j++)// 针对图像每列进行操作
	{		 
		// 计算该象素在源DIB中的坐标
		i0 = (LONG) (-((float) j) * fSina + ((float) i) * fCosa + f2 + 0.5);//旋转前的目标行
		j0 = (LONG) ( ((float) j) * fCosa + ((float) i) * fSina + f1 + 0.5);//旋转前的目标列		
		// 判断是否在源图范围内
		if( (j0 >= 0) && (j0 < lWidth) && (i0 >= 0) && (i0 < lHeight))
		{				
			BYTE mask = *((char *)lpDIBBits+lLineBytes*i0+j0/8)&(0x80>>j0%8);
			mask = mask ? (0x80 >> j%8) : 0;
			//lpSrc = (char *)lpDIBBits + lLineBytes * (lHeight - 1 - i0);//指向源图像中第i0行
			*((char*)lpNewDIBBits + lNewLineBytes*(i) + (j/8)) &= ~(0x80 >> j%8);
			*((char*)lpNewDIBBits + lNewLineBytes*(i) + (j/8)) |= mask;
		}
	}//列j结束	
}//行i结束