#include<stdio.h>

#define  Y_SIZE 480  /*縦が480ドット（この値は画像のサイズによってかえる）*/
#define  X_SIZE 640  /*横が640ドット（この値は画像のサイズによってかえる）*/

unsigned char image_in[Y_SIZE][X_SIZE][3];   /*入力画像配列*/
unsigned char image_out[Y_SIZE][X_SIZE][3];  /*出力画像配列*/
unsigned char bmph[54];  /*ビットマップヘッダー*/
unsigned char m[9];

int i,j,k;
unsigned char r[10];
unsigned char w[10];

int median();



main()
{
	FILE *fp;

	/*ファイルを入力*/
	printf("read_filename?");
	scanf("%s",r);
	fp=fopen(r,"rb");
	fread(bmph,sizeof(bmph),1,fp);
	fread(image_in,sizeof(image_in),1,fp);
	fclose(fp);

	printf("point1\n");

	/*image on median*/
	for(i=0;i<Y_SIZE;i++){
		for(j=0;j<X_SIZE;j++){
			for(k=0;k<3;k++){

				m[0]=image_in[i+1][j+1][k];
				m[1]=image_in[i+1][j-1][k];
				m[2]=image_in[i-1][j+1][k];
				m[3]=image_in[i+1][j-1][k];
				m[4]=image_in[i  ][j+1][k];
				m[6]=image_in[i+1][j-1][k];
				m[5]=image_in[i+1][j  ][k];
				m[7]=image_in[i-1][j  ][k];
				m[8]=image_in[i][j][k];
				image_out[i][j][k]=median(m);  /*メディアンフィルタ*/
			
			}
		}
	}
	printf("over median...........\n");
	printf("point2\n");

	/*ファイルを出力*/
	printf("write name?");
	scanf("%s",w);
	fp=fopen(w,"wb");
	fwrite(bmph,sizeof(bmph),1,fp);
	fwrite(image_out,sizeof(image_out),1,fp);
	fclose(fp);
	printf("point3\n");
	return(0);

}
/*9つの画素の中央値（メディアン）を求める*/
int median(unsigned char m[9]){

	int i,j,buf;

	for(j=0;j<8;j++){
		for(i=0;i<8;i++){
			if(m[i+1]<m[i]){
				buf=m[i+1];
				m[i+1]=m[i];
				m[i]=buf;
			}
		}
	}
	return m[4];
}