学习之路（二） 图像变换

题目：给定一幅数字图像，每个像元对应一个字节（灰度范围为0～255）。使用查找表技术编写程序实现对该图像的线性变换。要求结果影像的灰度范围为0～255，即将原始图像的灰度最小值变换为0，灰度最大值变换为255（通过这样对应关系可以求出变换系数A和B）。

数据：图像一幅，文件名为：test.raw。影像大小为1210行，1240列。

要求：输出查找表；输出变换后影像，变换后影像也按二进制方式存储。文件名为result.raw；

#include<stdio.h>
#include<stdlib.h>
#include<string.h>

#define height 1210
#define width 1240

typedef unsigned char  BYTE;

int getmax(BYTE Arr[height][width])
{
	int i,j;
	int m=Arr[0][0];
	for(i=0;i<height;i++)
		for(j=0;j<width;j++)
		{
			if(m<Arr[height][width])
			{
				m=Arr[height][width];
			}
		}
	return m;
}

int getmin(BYTE Arr[height][width])
{
	int i,j;
	int m=Arr[0][0];
	for(i=0;i<height;i++)
		for(j=0;j<width;j++)
		{
			if(m>Arr[height][width])
			{
				m=Arr[height][width];
			}
		}
		return m;
}

int main()
{
	FILE *fp=NULL;

	//**定义为全局变量，防止出现局部变量过大产生bug
	static BYTE x[height][width]; //原始图片像素二维数组
	static BYTE Arr[height][width]; //转换后图像像素二维数组
	static BYTE *ptr;

	int i,j;

	if( ( fp= fopen("test.raw","rb") ) == NULL )
	{
		printf("can not open the raw image " );
		return 0;
	}
	else
	{
		printf("read ok!\n");
	}

	//分配内存并将图像读到二维数组中
	ptr=(BYTE*)malloc(sizeof(BYTE)*width*height);
	for(i=0;i<height;i++)
		for(j=0;j<width;j++)
		{
			fread(ptr,1,1,fp);
			x[i][j]=*ptr; // 把图像输入到2维数组中,变成矩阵型式
			//printf("%d ",x[i][j]);
			ptr++;
		}

	fclose(fp);

	//寻找图像中的像素最大值和最小直
	int max=getmax(x);
	int min=getmin(x);

	printf("最大值:%d\n最小值:%d\n",max,min);

	double A=255.0/double(max-min);
	double B=-255.0*double(min)/double(max-min);

	printf("A:%lf\nB:%lf\n",A,B);

	//y存储变换后的对应像素值
	int y[256];
	for(i=0;i<256;i++)
	{
		y[i]=A*i+B;
	}

	//输出查找表
	if((fp=fopen("lut.txt","w"))==NULL)
	{
		printf("can not create the look up table\n");
		return 0;
	}

	for(i=0;i<256*2;i++)
	{
		if(i<256)
			fprintf(fp,"%-4d",i);
		else if(i==256)
			fprintf(fp,"\n%-4d",y[i-256]);
		else
			fprintf(fp,"%-4d",y[i-256]);
	}



	//根据查找表技术变换影像并输出变换后影像
	if( ( fp = fopen( "result.raw", "wb" ) ) == NULL )
	{
		printf("can not create the raw_image\n");
		return 0;
	}

	for( i = 0; i < height; i++ )
	{
		for( j = 0; j < width ; j ++ )
		{
			Arr[i][j]=y[x[i][j]];
			fwrite( &Arr[i][j], 1 , 1, fp );
		}
	}

	printf("end!\n");
	fclose(fp);

	return 0;
}