本文档最早发布于
1.256级灰度BMP文件读写的源代码!
首先要明白256级灰度BMP文件的格式
1.首先是一个14个字节的文件头,定义如下
typedef struct tagBITMAPFILEHEADER{
WORD bfType;
DWORD bfSize;
WORD bfReserved1;
WORD bfReserved2;
DWORD bfOffBits;
} BITMAPFILEHEADER, *PBITMAPFILEHEADER;
bfType是表明BMP文件类型的数据,在这里我们填入的是0x4d42,其实就是BM两个字,bfSize是文件大小,bfOffBits是文件头到数据块的偏移量,对于256级灰度图,就是1078个字节,后面会做描述
2.接下来是40个字节的是描述位图属性的40个字节
typedef struct tagBITMAPINFOHEADER{
DWORD biSize;
LONG biWidth;
LONG biHeight;
WORD biPlanes;
WORD biBitCount;
DWORD biCompression;
DWORD biSizeImage;
LONG biXPelsPerMeter;
LONG biYPelsPerMeter;
DWORD biClrUsed;
DWORD biClrImportant;
} BITMAPINFOHEADER, *PBITMAPINFOHEADER;
这里面只有biWidth表示宽度,biPlanes表示高度,biBitCount对于256级灰度正好是8
3.由于是256级灰度图,那么有256个调色板数据,每个调色板是如下定义的
typedef struct tagRGBQUAD {
BYTE rgbBlue;
BYTE rgbGreen;
BYTE rgbRed;
BYTE rgbReserved;
}RGBQUAD, *PRGBQUAD;
调色板数据其实告诉了显示器实际显示的时候的具体颜色,所以调色板长度是1024字节
4.最后是按行组织的图像数据,但这些数据并不是简单的按照图像的高度宽度w*h的数组数据这些数据最重要的特点是
a.按行组织,每行宽度是w,但是要进行4个字节的对齐。比如如果是图像宽度是253,那么数据对齐后一行还是有256个字节。对齐可以用下面的宏来计算
#define GET_ALIGN(x) (((x+3)/4)*4)
b.图像数据是倒行的,也就是数据第一行对应图像最后一行,最后一行数据对应第一行
图像的实际数据之前的偏移量是:sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+256*sizeof(RGBQUAD)=14+40+1024=1078个字节。
下面是实际的BMP文件输入输出函数代码:
/*********************************************/
从文件读入内存代码:
pw返回宽度;
ph返回高度;
函数返回内存指针;
/*********************************************/
unsigned char* read_bmp(const char* pszFileName, int* pw, int* ph)
{
BITMAPFILEHEADER bfh;
BITMAPINFOHEADER bmh;
FILE *fp;
unsigned char* pImg = NULL;
int i;
fp=fopen(pszFileName, "rb");//二进制打开
if(fp==NULL) return NULL;
fread(&bfh, sizeof(BITMAPFILEHEADER),1, fp);
fread(&bmh, sizeof(BITMAPINFOHEADER),1, fp);
//判断是否8bit的图像
if(bfh.bfType!=0x4d42 && bmh.biBitCount!=8) return NULL;
pImg=(unsigned char*)malloc(bmh.biWidth*bmh.biHeight);//开辟空间
*pw=bmh.biWidth;
*ph=bmh.biHeight;
for(i=0; iw-1 || y0<0 || y0>h-1) return -1;
for(i=y0; iw-3 || y>h-3) return;
pImg[(y-1)*w+x-2]=1;
pImg[ y *w+x-2]=1;
pImg[(y+1)*w+x-2]=1;
pImg[(y+2)*w+x-1]=1;
pImg[(y+2)*w+ x ]=1;
pImg[(y+2)*w+x+1]=1;
pImg[(y+1)*w+x+2]=1;
pImg[ y *w+x+2]=1;
pImg[(y-1)*w+x+2]=1;
pImg[(y-2)*w+x+1]=1;
pImg[(y-2)*w+ x ]=1;
pImg[(y-2)*w+x-1]=1;
}
/**************************************************/
画矩形函数,这个函数描绘了一个像素值为color的函数
矩形宽度为x0到x1,高度为y0到y1;*pImg指向位图数据区。
/**************************************************/
void draw_rect(unsigned char *pImg, int x0, int y0, int x1, int y1, int w, int h, int color)
{
int i, j;
if(x0<0 || x1>w-1 || y0<0 || y1>h-1) return;
for(i=y0; i>10;
}
同样的办法,插值两次就变成了二位图像插值算法。注意,我们临时的图像数组大小是512*512,插值后的图像
大小不能超过这个,否则会产生内存越界的错误。
void scale_img(unsigned char* pImgIn, int w1, int h1, unsigned char* pImgOut, int w2, int h2)
{
int i, j;
int i0, j0;
int f1, f2;
unsigned char pImg[512*512];
for(j=w2-1; j>=0; j--)
{
j0 = j*w1/w2;
f1 = 1024*j*w1/w2-1024*j0;
f2 = 1024-f1;
for(i=0; i>10);
}
}//宽度插值
for(i=h2-1; i>=0; i--)
{
i0 = i*h1/h2;
f1 = 1024*i*h1/h2-1024*i0;
f2 = 1024-f1;
for(j=0; j>10);
}
}//高度插值
}
下面是把我们的图像插值成原来的88%,125%
int step5(void)
{
int w, h, w2, h2, w3, h3;
unsigned char *pImg, *pImg2, *pImg3;
printf("step5: scale image test!\n");
pImg= read_bmp("test.bmp", &w, &h);
if(pImg==NULL) return -1;
w2=(int)(0.88*w); h2=(int)(0.88*h);
w3=(int)(1.25*w); h3=(int)(1.25*h);
pImg2=(unsigned char*)malloc(256*256);
pImg3=(unsigned char*)malloc(w3*h3);
scale_img(pImg, w, h, pImg2, w2, h2);
write_bmp("test5a.bmp", pImg2, w2, h2);
scale_img(pImg, w, h, pImg3, w3, h3);
write_bmp("test5b.bmp", pImg3, w3, h3);
free(pImg3);
free(pImg2);
return 1;
}
这是插值后的效果图,源代码已经同步更新