数字图像处理_邻域平均法和中值滤波处理_实验报告

fread(&hdr->compress_type, 4, 1, fp);
fread(&hdr->data_size, 4, 1, fp);
fread(&hdr->resol_hori, 4, 1, fp);
fread(&hdr->resol_vert, 4, 1, fp);
fread(&hdr->nr_colors, 4, 1, fp);
所遇问题
分析原因
解决方法
盐全部几种在图像的下面部分，如下图
temp=rand()%nr_pixels中temp是int型，而int型的数有效范围是-32767~32767，所以不能获得32768~65536之间的随机数
1、可采取“回避”法，换一种取随机中的算法
2、把temp换成长整型或浮点型。
运行了几次但每次生成的图像中噪声分布都一样
可以让系统获取当前时间，让其转换为int值，再从这个范围内通过对nr_pixels取整，获得一个从0~nr_pixels范围内的随机数
srand((unsigned int)time(NULL));
i = ( rand() %nr_pixels);
在做实验时，我通过反复测试，惊讶的发现在读取位图上每个像素的信息时，居然是从左下角开始读的（我之前一直以为是左上角）。但实验过程中我有个地方一直想不明白也费了很多时间去调试的地方是，在随机算法中，我先是用temp=rand()%nr_pixels来获取一个从0~65536的随机数，但只能对一半图像加上噪声，单步跟踪时也没发现什么问题，我只好换一种算法来获得随机数，即row=rand()%hdr->weight;column=rand()%hdr->height;bitmap[row*256+column] = 255;就一切正常了。可是这两种获取随机数的方法实质上应该是一样的啊……后天听室友做的时候也出现相同状况，经过分析后发现temp是int型，而int型的数有效范围是-32767~32767，所以不能获得32768~65536之间的随机数，才恍然大悟只能给图像的一半加噪声的原因。
其中,hdr.h用来定义位图头文件的结构体，hdr.c用来读取位图的头文件信息，ff.c用来对图像加上椒盐噪声。
5.3
位图头文件的数据结构如下：
/*以下为定义位图头文件的结构体*/
struct bmphdr
{
char signature[2];//文件类型
int size;//位图大小
short reserved[2];//两个保留字
fwrite(&hdr->data_size, 4, 1, fpnew);
fwrite(&hdr->resol_hori, 4, 1, fpnew);
fwrite(&hdr->resol_vert, 4, 1, fpnew);
fwrite(&hdr->nr_colors, 4, 1, fpnew);
}
hdr = (struct bmphdr *)malloc(sizeof(struct bmphdr));//分配空间
fread(hdr->signature, 2, 1, fp);
if (hdr->signature[0] != 'B' || hdr->signature[1] != 'M')//如果打开文件的扩展名不是.BMP，则说明不是
fread(&hdr->important_color, 4, 1, fp);
if (hdr->offset > 54) //获得位图有效信息，偏移量=54时，才读的是位图第一个像素
fread(&hdr->info, 1024, 1, fp);
fclose(fp);//关闭文件
return hdr;
row = rand()%hdr->height;
column = rand()%hdr->width;
bitmap[row*256+column] = (unsigned char)0;//椒为黑色，所以给该像素赋值为0
}
// 打开一个以输出文件名命名的文件,设为可写的二进制形式
fpnew = fopen(argv[2], "wb+");
column = rand()%hdr->width;
bitmap[row*256+column] =(unsigned char) 255;//盐为白色，所以给该像素赋值为255
}
for(i = 0; i < PEPPER; i++)//加入PEPPER“颗”椒
{
//在0~nr_pixels中随机生成一个数，即选择要加上椒的像素
}
6.
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <memory.h>
#include "hdr.h"
#define SALT 500 //加盐数量为500
#define PEPPER 500//加椒数量为500
struct bmphdr *hdr;
BMP位图文件
{
printf("Not a bmp file!\n");
return NULL;
}
//以下都是读位图头文件信息
fread(&hdr->size, 4, 1, fp);
fread(hdr->reserved, 4, 1, fp);
fread(&hdr->offset, 4, 1, fp);
fwrite(hdr->signature, 2, 1, fpnew);
fwrite(&hdr->size, 4, 1, fpnew);
fwrite(hdr->reserved, 4, 1, fpnew);
fwrite(&hdr->offset, 4, 1, fpnew);
fwrite(&hdr->hdr_size, 4, 1, fpnew);
FILE *fp, *fpnew;//定义两个文件指针分别用于提取原图像的数据和进行邻域平均后的图像
//定义主函数的参数包括：输入的位图文件名和输出的位图文件名
argc=3;
argv[1]="test.bmp";
argv[2]="testsaltpepper.bmp";
if (argc != 3)
{
printf("please input the name of input and out bitmap files\n");
exit(1);
}
// 获取位图文件相关信息
hdr = get_header(argv[1]);
if (!hdr) exit(1);
//以二进制可读方式打开输入位图文件
fp = fopen("test.bmp", "rb");
fwrite(&hdr->important_color, 4, 1, fpnew);
if (hdr->offset > 54)
fwrite(hdr->info, (hdr->offset - 54), 1, fpnew);
//直方图均衡化的数据(bitmap)赋值
fwrite(bitmap,nr_pixels,1,fpnew);
unsigned char *bitmap;
int main(int argc,char *argv[3])
{
int i, nr_pixels;//定义整数 i 用于函数循环时的，nr_pixels用于计算位图的像素个数
int sum = 0;//sum用来存放用模块后新的像素值
int row,column;//用来临时保存该产生噪声的像素
short nr_planes;//位平面数
short bits_per_pixel;//每像素的位数
int compress_type;//压缩类型
int data_size;//数据大小
int resol_hori;//水平
int resol_vert;//垂直
int nr_colors;
int important_color;
{
FILE *fp;
struct bmphdr *hdr;
fp = foቤተ መጻሕፍቲ ባይዱen(filename, "rb");//打开位图文件
if (!fp)//位图文件打开出错时给出提示
{
printf("File open error or such file does not exist!\n");
return NULL;
四、实验环境
操作系统：Windows XP
试验平台：VC++ 6.0
五、设计思路
5.1
5.1．1、
基本思想是:生成一些随机数，根据这些随机数把对应的像素灰度级设为255。
5.1．
基本思想是:生成一些随机数，根据这些随机数把对应的像素灰度级设为0。
5.2
该程序一共包含3个文件hdr.h,hdr.c,和ff.c
char info[1024];
};
struct bmphdr *get_header(char filename[]);