工作用到了这方面的知识,记录一下,用大白话说明一下。
2023年3月30日新增:
BMP图片格式解析,快速求数组的中值,可以了解一下BMP图片的格式,快速求数组中值 在滤波中可以减少很多时间,也推荐看下!
一、BMP图片
1、一个BMP格式的图片可以分为三块,文件头信息,位图信息,位图数据(有的带调色板,可选)。
BMP文件头: 存储了该文件的类型,大小和存放信息等。
位图信息: 放BMP图像的高,宽,bpp等信息。
位图数据: 就是每个像素的值。
BMP图片的其余详细信息在代码的/* 头文件定义 */
下图是BMP格式图片图像的存储方式:从左到右,从下到上,此处存放的应该是位图数据
2、具体操作
其他知识:什么是中值滤波?
由上图可知,BMP图片的像素存放方式很像一个二维平面坐标系,把下图P7想象成坐标原点(0,0),在这个3*3的区域内,将这九个数排序,假设从小到大排序是:P2,P4,P8,P3,P9,P1,P5,P6,P7。中值就是P9,把图中的P5换成P9。一个像素的“滤波”就完成了,其余像素的操作跟这个一致。
有的小伙伴会发现,这样搞,那边界的,如P7,P8,P9怎么替换?这就用到边界处理的,边界拓展——复制图像边界来处理?啥意思,即用次外圈的像素替换最外圈的像素,看不懂没关系,举个例子。下图,P1,P4分别用P2,P5替换,同理P8,P9分别用P5,P6替换,这里P1,P4,P7,P8,P9都是最外圈像素,P2,P5,P6都是次内圈像素。那四个角呢,跟它直接挨着的都可以替换它,如P7,可以用P4,P5,P8替换,因为刚处理过一次后P4,P5,P8都一样了嘿。
思路:文件嘛,在计算机看来都是二进制的0,1。如何读取呢?用fread读取BMP图片的数据,读取到的数据存放到数组data中,从data中取出9个数据,排序后取中值替换即可。
代码大部分都是直接参考别人的,只有边界处理是自己加的。
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
/* 颜色表定义,注释掉7-13行代码亦可,因为代码未用到*/
typedef unsigned char BYTE;
typedef struct tagRGBQUAD{
BYTE rgbBlue;
BYTE rgbGreen;
BYTE rgbRed;
BYTE rgbReserved;
} RGBQUAD;
/* 定义头文件型 */
#pragma pack(1)//声明1字节对齐,因为BMP文件是1字节保存的文件,而系统默认4字节对齐。
typedef struct{
unsigned char id1;//位图文件的类型,必须为BM(占用0-1字节)
unsigned char id2;
unsigned int filesize;//位图文件的大小,以字节为单位(2-5字节)
unsigned int reserved;// 位图文件保留字,必须为0(6-9字节)
unsigned int bitmapdataoffset;//位图数据的起始位置,以相对于位图(10-13字节)
unsigned int bitmapheadersize;//BMP头的大小,固定为40(14-17字节)
unsigned int width;//图片宽度;以像素为单位(18-21字节)
unsigned int height;//图片高度;以像素为单位(22-25字节)
unsigned short planes;//图片位面数,必须为1(26-27字节)
unsigned short bitperpixel;//每个像素所需的位数,每个像素所需的位数(28-29字节)
//只能是以下几个数:1(双色),4(16色),8(256色)或24(真彩色) 灰度级
unsigned int compression;//是否压缩(30-33字节)
//只能是以下几个数:0(不压缩),1(BI_RLE8压缩类型),2(BI_RLE4压缩类型)
unsigned int bitmapdatasize;//位图的大小,以字节为单位(34-37字节)
unsigned int hresolution;//位图水平分辨率,每米像素数(38-41字节)
unsigned int vresolution;//位图垂直分辨率,每米像素数(42-45字节)
unsigned int colors;//位图实际使用的颜色表中的颜色数(46-49字节)
unsigned int importantcolors;//位图显示过程中重要的颜色数(50-53字节)
//unsigned int bmiColors[1];//调色板;(54 - 57字节)
unsigned char palette[256][4];//调色板 占256*4=1024字节
}BMPheader;//总大小40+14+1024=1078字节
typedef struct{
BMPheader* bmpheader ;
unsigned char* bitmapdata;//图片数据;
}BMPheaderfile;
/*
求文件长度的函数
*/
long getfilesize(FILE *f)
{
long pos,len;
pos=ftell(f);//ftell函数用于得到文件指针当前位置相对于文件首的偏移字节数
fseek(f,0,SEEK_END);//fseek函数用于移动文件指针相对于SEEK_END的偏移量为0,相当于移动到文件最后
len=ftell(f);//len就是文件的长度
fseek(f,pos,SEEK_SET);//将文件指针移动到原来的地方
return len;
}
/*
主函数
*/
int main()
{
BMPheaderfile *output=(BMPheaderfile*)malloc(sizeof(BMPheaderfile));//定义一个输出指针
unsigned char *data=NULL;
FILE *fpr,*fpw;
/*
打开文件
*/
if((fpr=fopen("start1.bmp","rb"))==NULL)
{
printf("cannot open this file");
exit(0);
}
if((fpw=fopen("end1.bmp","wb"))==NULL)
{
printf("cannot wirte this file");
exit(0);
}
long length=getfilesize(fpr);
printf("文件的长度为%ld\n",length);
printf("文件的头长度为%ld\n",sizeof(BMPheader));
data=(unsigned char*)malloc(length*sizeof(char));//分配空间
/* 读文件,从fpr指向的文件中读出1个length长度(即图片的大小)的数据到data所指的内存空间去 */
if(0==fread(data,1,length,fpr))
{
printf("read failed\n");
exit(0);
}
fclose(fpr);//释放指针
output->bmpheader=(BMPheader*)malloc(sizeof(BMPheader));
/* 从data中拷贝sizeof(BMPheader)大小到output->bmpheader */
memcpy(output->bmpheader,data,sizeof(BMPheader));
/*
打印出图像中头文件的信息
*/
int height=output->bmpheader->height;
int width=output->bmpheader->width;
printf("filesize is %d\n",output->bmpheader->filesize);
printf("该图像每个像素所需要的位数:%d\n",output->bmpheader->bitperpixel);
printf("height is %d\n",output->bmpheader->height);
printf("width is %d\n",output->bmpheader->width);
data=data+sizeof(BMPheader);
output->bitmapdata=data;
/*
中值滤波算法(选择3×3的滑动窗口)
*/
unsigned char pixel[9]={0};//滑动窗口的像素值,初始为0
unsigned char mid;//中值
unsigned char temp;//中间变量
unsigned char border;
int flag;
int m,i,j,x,h,w,y;
for(j=1;j<height-1;j++)
{
for(i=1;i<width-1;i++)
{
//将3×3滑动窗口中的所有像素值放入pixel[m]
m=0;
for(y=j-1;y<=j+1;y++)
{
for(x=i-1;x<=i+1;x++)
{
pixel[m]=data[y*width+x];
m=m+1;
}
}
//让一位数组pixel[9]进行降序排列
do
{
flag=0;//循环结束的标志
for(m=0;m<9;m++)
{
if(pixel[m]<pixel[m+1])
{
temp = pixel[m];
pixel[m] = pixel[m+1];
pixel[m+1] = temp;
flag=1;
}//if
}//for
} while (flag==1);
mid=pixel[4];
output->bitmapdata[width*j+i]=mid;
}
}
/* 边界问题:拓展边界——复制图像边界 */
/* 下边界 */
for(int q = 2; q < width; q++)
{
output->bitmapdata[q]=output->bitmapdata[width+q];
}
/* 上边界 */
for(int w = width*(width-2)+2; w < width*(width-1)-1; w++)
{
output->bitmapdata[w+width]=output->bitmapdata[w];
}
/* 左边界 */
for(int e = 0; e < height; e++)
{
output->bitmapdata[e*height]=output->bitmapdata[e*height+1];
}
/* 右边界 */
for(int r = 0; r < height; r++)
{
output->bitmapdata[r*height-1]=output->bitmapdata[r*height-2];
}
/*
N=5的中值滤波线段状
*/
/*
for(i=0;i<height;i++)
{
for(j=2;j<width-1;j++)
{
m=0;
for(x=j-2;x<=j+2;x++)
pixel[m++]=data[i*width+x];
for(h=0;h<5;h++)
for(w=h+1;w<5;w++)
{
if(pixel[h]>pixel[w])
{
temp=pixel[w];
pixel[w]=pixel[h];
pixel[h]=temp;
}
}
output->bitmapdata[i*width+j]=pixel[2];
}
}
*/
/*
保存图像文件
*/
fseek(fpw,0,0); //fseek(fpw,0,SEEK_SET)
fwrite(output->bmpheader,1,sizeof(BMPheader),fpw);//写入图像的头文件
fwrite(output->bitmapdata,1,length-sizeof(BMPheader),fpw);//写入图像的数据信息
fclose(fpw);//释放指针
return 0;
}
原图:无边界处理:
有边界处理: