作业程序源码及说明及指南.docVIP

作业程序源码及说明 冯斌 计算机科学与技术一班 matlab中的加噪处理 说明：在matlab command window中直接处理。lena.bmp为原始图象，经过matlab加高斯白噪声得到lena20.bmp，该高斯白噪声标准差为20。程序如下： I=imread(d:\lena.bmp); I=double(I); max_I=max(max(I)); min_I=min(min(I)); II=(I-min_I)/(max_I-min_I); %图像归一化 III=imnoise(II,gaussian,0,20*20/(255*255)); %加入方差为20的高斯白噪声 im_noise=min_I+III*(max_I-min_I); %图像反归一化 imshow(III); 二、核心程序myMain.c /*对一个经过matlab加高斯白噪声的图象并手工选取规则大小lena20.bmp进行3*3区域算术均值滤波 */ /*该图象的原始图象是一个128*128的bmp：lenabmp*/ #include stdlib.h #include stdio.h #include math.h #define SIZEPALETTE 256 /* 位图彩色表个数256，占用空间大小256*4 */ //#define SIZEBMP 26062 /* 验证数据位图阵列数据大小166*157 该大小生成不可绘制的位图*/ #define SIZEBMP 16384 /*位图阵列数据大小128*128*/ typedef unsigned short int WORD; typedef unsigned long int DWORD; typedef unsigned char BYTE ; typedef struct DDWORD { WORD highhigh; WORD highlow; WORD lowhigh; WORD lowlow; }DDWORD; /*宏定义bmp相关数据结构*/ typedef struct tagBITMAPFILEHEADER { DWORD bfType; DDWORD bfSize; DWORD bfReserved1; DWORD bfReserved2; DDWORD bfOffBits; } BITMAPFILEHEADER; /*位图文件头*/ typedef struct tagBITMAPINFOHEADER{ DDWORD biSize; DDWORD biWidth; DDWORD biHeight; DWORD biPlanes; DWORD biBitCount; DDWORD biCompression; DDWORD biSizeImage; DDWORD biXPelsPerMeter; DDWORD biYPelsPerMeter; DDWORD biClrUsed; DDWORD biClrImportant; } BITMAPINFOHEADER; /*位图信息头格式定义*/ typedef struct tagRGBQUAD{ BYTE rgbBlue; BYTE rgbGreen; BYTE rgbRed; BYTE rgbReserved; } RGBQUAD;/*位图彩色表定义，构成包括一个3字节的RGBTRIPLE和保留字节 */ /*bmp数据结构宏定义结束*/ BITMAPFILEHEADER bmfHeader; /*位图文件头*/ BITMAPINFOHEADER bmiHeader; /*位图信息头 */ RGBQUAD m_Palette[SIZEPALETTE]; /*位图彩色表数组*/ BYTE m_ImageOld[SIZEBMP]; /*含噪图象阵列数据*/ BYTE m_ImageNew[SIZEBMP]; /*滤波图象阵列数据*/ main() { int i,j; int m,n; int width=128; int height=128; FILE *fpin; FILE *fpout; float coffi=1.0/9.0; double sum; /*读取文件fpin-lena20.bmp，将彩色表信息读进m_Palette中，阵列数据读入m_ImageOld中*/ fpin=fopen(E:\\myccspro\\