圖像平滑化處理鄰域平均法C語言實現

1、圖像平滑化處理 鄰域平均法C語言實現鄰域平均法是簡單的空域圖像平滑化處理方法，對圖像噪聲的去除有一定的效果。#ifndef BMP_H_4_INCLUDED#define BMP_H_4_INCLUDED typedef unsigned short WORD；typedef unsigned long DWORD；typedef long LONG；typedef unsigned char BYTE；typedef struct tagBITMAPFILEHEADER/bmfh WORD bfType；DWORD bfSize；WORD bfReserved1；WORD bfReserve

2、d2；DWORD bfOffBits；BITMAPFILEHEADER；typedef struct tagBITMAPINFOHEADER/bmih DWORD biSize；LONG biWidth；LONG biHeight；WORD biPlanes；WORD biBitCount；DWORD biCompression；DWORD biSizeImage；LONG biXPelsPerMeter；LONG biYPelsPerMeter；DWORD biClrUsed；DWORD biClrImportant；BITMAPINFOHEADER；typedef struct tagRG

3、BQUAD/rgbq BYTE rgbBlue；BYTE rgbGreen；BYTE rgbRed；BYTE rgbReserved；RGBQUAD；typedef struct tagBITMAPINFOBITMAPINFOHEADER bmiHeader；RGBQUAD bmiColors1；BITMAPINFO；#include stdio.h#include stdlib.h#include string.h#include malloc.h#include ctype.h#include process.h#include"BMP_4.h"BITMAPFILEHE

4、ADER bmfh；BITMAPINFOHEADER bmih；BYTE*imgData；void readData()；int ReadFileHeader(BITMAPFILEHEADER*)；int ReadInfoHeader(BITMAPINFOHEADER*)；int CreatePalette(RGBQUAD)；int ReadPixelData(BYTE*)；LONG GetLineBytes(int,int)；void NeighborAverage(double)；int SaveAsImage(char*)；int main()int i；char saveasfilep

5、ath256；DWORD dwLineBytes；double Error；readData()；dwLineBytes=GetLineBytes(bmih.biWidth,bmih.biBitCount)；printf("Input error！n")；scanf("%f",&Error)；NeighborAverage(Error)；printf("Save as another path(ex.d：/poon.bmp)n")；scanf("%s",saveasfilepath)；i=SaveAsIma

6、ge(saveasfilepath)；if(i=-1)printf("Error：failed to save the image.n")；return 0；void readData()int i,k,h；DWORD dwLineBytes；i=ReadFileHeader(&bmfh)；if(i=0)printf("Read file header successful！n")；k=ReadInfoHeader(&bmih)；if(k=0)printf("Read info header successful！n"

7、)；dwLineBytes=GetLineBytes(bmih.biWidth,bmih.biBitCount)；imgData=(BYTE*)malloc(dwLineBytes*bmih.biHeight*sizeof(BYTE)；h=ReadPixelData(imgData)；if(h=0)printf("Read pixel data successful！n")；if(i=0&&k=0&&h=0)printf("Read datas successful！n")；int ReadFileHeader(BITMA

8、PFILEHEADER*bmfh)FILE*dataFile；dataFile=fopen("shuaige.bmp","r")；if(fread(&bmfh-bfType,sizeof(WORD),1,dataFile)！=1)printf("Can not read bfType in the file header.n")；fclose(dataFile)；return-1；if(fread(&bmfh-bfSize,sizeof(DWORD),1,dataFile)！=1)printf("Can no

9、t read bfSize in the file header.n")；fclose(dataFile)；return-1；if(fread(&bmfh-bfReserved1,sizeof(WORD),1,dataFile)！=1)printf("Can not read bfReserved1 in the file header.n")；fclose(dataFile)；return-1；if(fread(&bmfh-bfReserved2,sizeof(WORD),1,dataFile)！=1)printf("Can not r

10、ead bfReserved2 in the file header.n")；fclose(dataFile)；return-1；if(fread(&bmfh-bfOffBits,sizeof(DWORD),1,dataFile)！=1)printf("Can not read bfOffBits in the file header.n")；fclose(dataFile)；return-1；fclose(dataFile)；return 0；int ReadInfoHeader(BITMAPINFOHEADER*bmih)FILE*dataFile；d

11、ataFile=fopen("shuaige.bmp","r")；fseek(dataFile,14,SEEK_SET)；if(fread(&bmih-biSize,sizeof(DWORD),1,dataFile)！=1)printf("Can not read biSize in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biWidth,sizeof(LONG),1,dataFile)！=1)printf("Can not r

12、ead biWidth in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biHeight,sizeof(LONG),1,dataFile)！=1)printf("Can not read biHeight in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biPlanes,sizeof(WORD),1,dataFile)！=1)printf("Can not read biPlane

13、s in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biBitCount,sizeof(WORD),1,dataFile)！=1)printf("Can not read biBitCount in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biCompression,sizeof(DWORD),1,dataFile)！=1)printf("Can not read biCompr

14、ession in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biSizeImage,sizeof(DWORD),1,dataFile)！=1)printf("Can not read biSizeImage in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biXPelsPerMeter,sizeof(LONG),1,dataFile)！=1)printf("Can not rea

15、d biXPelsPerMeter in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biYPelsPerMeter,sizeof(LONG),1,dataFile)！=1)printf("Can not read biYPelsPerMeter in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biClrUsed,sizeof(DWORD),1,dataFile)！=1)printf(&quo

16、t;Can not read biClrUsed in the info header.n")；fclose(dataFile)；return-1；if(fread(&bmih-biClrImportant,sizeof(DWORD),1,dataFile)！=1)printf("Can not read biClrImportant in the info header.n")；fclose(dataFile)；return-1；fclose(dataFile)；return 0；int CreatePalette(RGBQUAD pal)int i；i

17、f(sizeof(pal)/sizeof(RGBQUAD)！=256)printf("The size of the palette must be 256.n")；return-1；for(i=0；i 256；i+)pali.rgbBlue=i；pali.rgbGreen=i；pali.rgbRed=i；pali.rgbReserved=0；return 0；int ReadPixelData(BYTE*imgData)BYTE*data；FILE*dataFile；DWORD dwLineBytes；dwLineBytes=GetLineBytes(bmih.biWid

18、th,bmih.biBitCount)；data=(BYTE*)malloc(dwLineBytes*bmih.biHeight*sizeof(BYTE)；dataFile=fopen("shuaige.bmp","rb")；fseek(dataFile,bmfh.bfOffBits,SEEK_SET)；if(fread(data,dwLineBytes*bmih.biHeight*sizeof(BYTE),1,dataFile)！=1)printf("Can not read the pixel data.n")；free(data

19、)；fclose(dataFile)；return-1；memcpy(imgData,data,dwLineBytes*bmih.biHeight*sizeof(BYTE)；free(data)；fclose(dataFile)；return 0；LONG GetLineBytes(int imgWidth,int bitCount)return(imgWidth*bitCount+31)/32*4；void NeighborAverage(double Error)int dwLineBytes；int NewdwLineBytes；int NewHeight；BYTE*NewData；do

20、uble gray；int i,j；dwLineBytes=GetLineBytes(bmih.biWidth,bmih.biBitCount)；NewdwLineBytes=dwLineBytes+2；NewHeight=bmih.biHeight+2；NewData=(BYTE*)malloc(NewdwLineBytes*NewHeight*sizeof(BYTE)；for(i=0；i NewHeight；i+)for(j=0；j NewdwLineBytes；j+)NewDataNewdwLineBytes*i+j=0；for(i=0；i bmih.biHeight；i+)for(j=

21、0；j dwLineBytes；j+)NewDataNewdwLineBytes*(i+1)+j+1=imgDatadwLineBytes*i+j；for(i=1；i NewHeight-1；i+)for(j=1；j NewdwLineBytes-1；j+)gray=(double)(NewDataNewdwLineBytes*(i-1)+j-1+NewDataNewdwLineBytes*i+j-1+NewDataNewdwLineBytes*(i+1)+j-1+NewDataNewdwLineBytes*(i-1)+j+NewDataNewdwLineBytes*(i+1)+j+NewDa

22、taNewdwLineBytes*(i-1)+j+1+NewDataNewdwLineBytes*i+j+1+NewDataNewdwLineBytes*(i+1)+j+1)/8；if(abs(double)(NewDataNewdwLineBytes*i+j)-gray)Error)imgDatadwLineBytes*(i-1)+j-1=(int)gray；int SaveAsImage(char*filepath)FILE*dataFile；RGBQUAD pal256；int i；DWORD dwLineBytes；dwLineBytes=GetLineBytes(bmih.biWid

23、th,bmih.biBitCount)；dataFile=fopen(filepath,"wb")；if(！dataFile)printf("Error：Can not open the file：%sn",filepath)；return-1；for(i=0；i 256；i+)pali.rgbReserved=0；pali.rgbBlue=i；pali.rgbGreen=i；pali.rgbRed=i；if(fwrite(&bmfh.bfType,sizeof(WORD),1,dataFile)！=1)printf("Can not

24、write bfType in the file header.n")；fclose(dataFile)；return-1；if(fwrite(&bmfh.bfSize,sizeof(DWORD),1,dataFile)！=1)printf("Can not write bfSize in the file header.n")；fclose(dataFile)；return-1；if(fwrite(&bmfh.bfReserved1,sizeof(WORD),1,dataFile)！=1)printf("Can not write bf

25、Reserved1 in the file header.n")；fclose(dataFile)；return-1；if(fwrite(&bmfh.bfReserved2,sizeof(WORD),1,dataFile)！=1)printf("Can not write bfReserved2 in the file header.n")；fclose(dataFile)；return-1；if(fwrite(&bmfh.bfOffBits,sizeof(DWORD),1,dataFile)！=1)printf("Can not wri

26、te bfOffBits in the file header.n")；fclose(dataFile)；return-1；if(fwrite(&bmih.biSize,sizeof(DWORD),1,dataFile)！=1)printf("Can not write biSize in the info header.n")；fclose(dataFile)；return-1；if(fwrite(&bmih.biWidth,sizeof(LONG),1,dataFile)！=1)printf("Can not write biWidt

27、h in the info header.n")；fclose(dataFile)；return-1；if(fwrite(&bmih.biHeight,sizeof(LONG),1,dataFile)！=1)printf("Can not write biHeight in the info header.n")；fclose(dataFile)；return-1；if(fwrite(&bmih.biPlanes,sizeof(WORD),1,dataFile)！=1)printf("Can not write biPlanes in t

28、he info header.n")；fclose(dataFile)；return-1；if(fwrite(&bmih.biBitCount,sizeof(WORD),1,dataFile)！=1)printf("Can not write biBitCount in the info header.n")；fclose(dataFile)；return-1；if(fwrite(&bmih.biCompression,sizeof(DWORD),1,dataFile)！=1)printf("Can not write biCompression in the info header.n")；fclose(dataFile)；return-1；if(fwrite(&bmih.biSizeImage,sizeof(DWORD),1,dataFile)！=1)printf("Can not write biSizeImage in the info header.n")；fclose(dataFile)；return-1；if(fwrite(&bmih.biXPelsPerMeter,sizeof(LONG),1,dataFile)！=1)pri