数字图像处理各种算法 C语言实现

包括图像空域滤波、频域滤波、几何变换(旋转、缩放、裁剪、镜像、平移)、傅里叶变换、余弦变换等等算法

#define _CRT_SECURE_NO_DEPRECATE
#include “gray_enhance.h“
#include “SpatialFiltering.h“
#include “FrequencyFiltering.h“

int main（）
{
	BITMAPFILEHEADER fileHeader;
	BITMAPINFOHEADER infoHeader;
	BYTE *extra = NULL;
	RGBQUAD *Quad8 = NULL *Quad4 = NULL;   //调色板

	RGB *img24 = NULL;
	RGB *img24_Out = NULL;
	BYTE *img8 = NULL;
	BYTE *img8_Out = NULL;

	FILE* pfin = NULL;
	FILE* pfout = NULL;
	/*FILE* txtfp = NULL;*/

	char *temp = NULL;

	int height width;
	int height_out width_out;
	int size;
	int i;
	int dif;
	int MirrorMode;

	int mode;				  //计算模式
	int intermode;			  //插值方式
	int template_size;		  //卷积模板大小
	int Spatial_filtermode;	  //空域滤波类型
	int trans_mode;			  //彩色转灰度转换方式

	int Freq_Filter_Mode;     //频域滤波类型
	int Freq_n;				  //滤波阶数
	int Freq_filter_class;    //频域滤波器选择

	int Up_Down Left_Right;  //图像平移位数
	int Crop_left Crop_right;//图像裁剪左右端
	int Crop_down Crop_up;	  //图像裁剪上下端

	double LogC;
	double ExpC ExpGamma;
	double LinearK LinearB;
	double RotAngle;
	double WidthScale HeightScale;
	double D0;			//频域滤波截止频率
	double D1;

	char file_in[800];
	char file_out[800];

	strcpy（file_in “C:\\zkyfile\\image_processing_test\\lenna.bmp“）;
	strcpy（file_out “C:\\zkyfile\\image_processing_test\\lenna_out.bmp“）;

	//double c;
	//int a b;
	//int d;
	////a = 103;

	//for （a = 1; a < 257; a++）
	//{
	//	b = 256;
	//	c = （a + 0.0） / （b + 0.0）;
	//	d = （int）（c*b）;
	//	printf（“%lf %d %d %d\n“ cada-d）;
	//}
	
	if （strcmp（file_in file_out） == 0）
	{
		printf（“输入文件和输出文件不能相等！\n“）;
		getchar（）;
		system（“pause“）;
		return 0;
	}

	pfin = fopen（file_in “rb“）;
	if （pfin == NULL）
	{
		printf（“打开文件失败！\n“）;
		getchar（）;
		system（“pause“）;
		return 0;
	}

	while （1）
	{
		printf（“\n\n\t\t欢迎来到图像处理系统！\n\n\n“）;
		printf（“线性：1    指数：2    对数：3    旋转：4    缩放：5    镜像：6\n“）;
		printf（“直方图均衡：7    均值滤波：8   彩色转灰度图（彩色图）：9   中值滤波：10\n  “）;
		printf（“空域滤波集成函数：11    傅里叶变换图：12    频域滤波集成函数：13\n“）;
		printf（“图像平移：14    图像裁剪：15   离散余弦变换:16    DCT与IDCT测试：17\n“）;
		printf（“退出：其他\n\n“）;
		printf（“请输入你要进行的操作：\n“）;
		scanf（“%d“ &mode）;

		ExpC = 1;
		ExpGamma = 3;
		LogC = 1;
		LinearK = 2;
		LinearB = 20;
		MirrorMode = 2;

		WidthScale = 1.0;
		HeightScale = 1.0;
		//输出图像长宽与输入长宽之比

		RotAngle = 45;

		template_size = 3;
		Spatial_filtermode = 1;
		trans_mode = 1;

		Freq_Filter_Mode = 1;
		Freq_filter_class = 1;
		Freq_n = 3;
		D0 = 30.0;
		D1 = D0 + 10.0;

		Up_Down = 0;
		Left_Right = 0;

		Crop_left = 1;
		Crop_right = 10;
		Crop_down = 1;
		Crop_up = 10;

		intermode = 1;

		switch （mode）
		{
		case Linear_Trans:
			printf（“请输入斜率: “）;
			scanf（“%lf“ &LinearK）;
			//必须用“%lf“，用“%f“会出问题
			printf（“请输入截距: “）;
			scanf（“%lf“ &LinearB）;
			break;
		case Exp_Trans:
			printf（“请输入常数项：“）;
			scanf（“%lf“ &ExpC）;
			printf（“请输入指数项：“）;
			scanf（“%lf“ &ExpGamma）;
			break;
		case Log_Trans:
			printf（“请输入常数c：“）;
			scanf（“%lf“ &LogC）;
		

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    ..A..H.     78336  2018-12-08 16:08  image_test1\.vs\image_test1\v15\.suo

     文件   43921408  2018-12-08 16:08  image_test1\.vs\image_test1\v15\Browse.VC.db

     文件   46661632  2018-12-08 16:08  image_test1\.vs\image_test1\v15\ipch\AutoPCH\3bf1d8971df8d01b\IMAGE_TEST1.ipch

     文件   42926080  2018-10-14 13:55  image_test1\.vs\image_test1\v15\ipch\AutoPCH\8c7c64b799a6e370\IMAGE_TEST1.ipch

     文件   46596096  2018-11-12 20:54  image_test1\.vs\image_test1\v15\ipch\AutoPCH\a0f1fa04430cab65\SPATIALFILTERING.ipch

     文件   46530560  2018-11-19 20:17  image_test1\.vs\image_test1\v15\ipch\AutoPCH\b5e694f83c549983\MAIN_HEAD.ipch

     文件   46465024  2018-11-01 21:09  image_test1\.vs\image_test1\v15\ipch\AutoPCH\e8ce77ca57a86b3b\IMAGE_TEST1.ipch

     文件   46530560  2018-11-12 20:48  image_test1\.vs\image_test1\v15\ipch\AutoPCH\e9af32aaa836a557\GRAY_ENHANCE.ipch

     文件   69140480  2018-10-13 13:00  image_test1\.vs\image_test1\v15\ipch\AutoPCH\eebd00148f5a872e\READ_IMAGE.ipch

     文件   46530560  2018-11-12 20:47  image_test1\.vs\image_test1\v15\ipch\AutoPCH\fd68c4f233e24ed2\GRAY.ipch

     文件     380928  2018-10-14 23:05  image_test1\.vs\image_test1\v15\Solution.VC.db

     文件      42496  2018-10-14 13:55  image_test1\Debug\image_test1.exe

     文件     495212  2018-10-14 13:55  image_test1\Debug\image_test1.ilk

     文件     643072  2018-10-14 13:55  image_test1\Debug\image_test1.pdb

     文件       1135  2018-10-14 13:55  image_test1\image_test1\Debug\image_test1.log

     文件      45546  2018-10-14 13:55  image_test1\image_test1\Debug\image_test1.obj

     文件        808  2018-10-14 13:55  image_test1\image_test1\Debug\image_test1.tlog\CL.command.1.tlog

     文件      26158  2018-10-14 13:55  image_test1\image_test1\Debug\image_test1.tlog\CL.read.1.tlog

     文件        616  2018-10-14 13:55  image_test1\image_test1\Debug\image_test1.tlog\CL.write.1.tlog

     文件        221  2018-10-14 13:55  image_test1\image_test1\Debug\image_test1.tlog\image_test1.lastbuildstate

     文件       1266  2018-10-14 13:55  image_test1\image_test1\Debug\image_test1.tlog\link.command.1.tlog

     文件       3122  2018-10-14 13:55  image_test1\image_test1\Debug\image_test1.tlog\link.read.1.tlog

     文件        580  2018-10-14 13:55  image_test1\image_test1\Debug\image_test1.tlog\link.write.1.tlog

     文件     429056  2018-10-14 13:55  image_test1\image_test1\Debug\vc141.idb

     文件     438272  2018-10-14 13:55  image_test1\image_test1\Debug\vc141.pdb

     文件      19982  2018-11-27 15:50  image_test1\image_test1\FrequencyFiltering.h

     文件      24361  2018-11-27 15:50  image_test1\image_test1\gray_enhance.h

     文件      14426  2018-11-29 21:29  image_test1\image_test1\image_test1.cpp

     文件       6473  2018-11-19 20:26  image_test1\image_test1\image_test1.vcxproj

     文件       1370  2018-11-19 20:26  image_test1\image_test1\image_test1.vcxproj.filters

............此处省略64个文件信息