#include?“subfunction.h“

#define?MU?50000
#define?BETA_MAX?128
#define?ITERMAX?5

Mat?Image?=?imread（“E:/pic/lena.png“?0）;
Mat?D1?=?（Mat_（1?2）?<Mat?D2?=?（Mat_（2?1）?<Size?sizeg?=?Image.size（）;
Mat?K?=?Mat::ones（15?15?CV_32FC1）/225;
Mat?noise（sizeg?CV_32FC1）;
Mat?Dx（sizeg?CV_32FC1）;
Mat?Dy（sizeg?CV_32FC1）;
int?beta?=?1;
int?beta_rate?=?2;
int?iter?=?1;
int?outIter?=?1;
Mat?g?Fg?u?w1?w2?Fk?FabsDx?FabsDy?temp[2]?conjDx?conjDy?conjFk?Dux?Duytemp1temp2temp3FabsFkFu;
int?main（）{
	system（“color?3f“）;
/*************初始化*************/
	if?（Image.empty（））{
		cout?<		return?-1;
	}
	Image.convertTo（Image?CV_32FC1）;
	randn（noise?0?1）;
	g?=?blurcircul（Image15Image.rowsImage.cols）;
	g?=?g?+?0.001*noise;
	Dx?=?psf2otf（D1?sizeg）;
	Dy?=?psf2otf（D2?sizeg）;
	Fk?=?psf2otf（K?sizeg）;
	conjDx?=?conj（Dx）;
	conjDy?=?conj（Dy）;
	conjFk?=?conj（Fk）;
	Fg?=?fft2（g）;
	split（Dx?temp）;
	magnitude（temp[0]?temp[1]?FabsDx）;
	split（Dy?temp）;
	magnitude（temp[0]?temp[1]?FabsDy）;
	split（Fk?temp）;
	magnitude（temp[0]?temp[1]?FabsFk）;
	pow（FabsDx?2?FabsDx）;//取模值的平方
	pow（FabsDy?2?FabsDy）;//取模值的平方
	pow（FabsFk?2?FabsFk）;//取模值的平方
	g.copyTo（u）;
/*******************************/
	while?（beta		iter?=?1;
		while?（iter	/*******updata?w1?and?w2*******/
			Dux?=?diffX（u）;
			Duy?=?diffY（u）;
			two_D_shrinkage（Dux?Duy?beta?w1?w2）;

	/***********updata?u***********/
			float?gamma?=?MU?/?beta;
			temp1?=?fft2（w1）;
			temp2?=?fft2（w2）;
			temp1?=?complexMul（conjDx?temp1）;
			temp2?=?complexMul（conjDy?temp2）;
			temp3?=?complexMul（conjFk?Fg）;
			temp1?=?temp1?+?temp2?+?gamma*temp3;
			temp2?=?FabsDx?+?FabsDy?+?gamma*FabsFk;
			split（temp1?temp）;
			temp[0]?=?temp[0]?/?temp2;
			temp[1]?=?temp[1]?/?temp2;
			merge（temp?2?Fu）;
			dft（Fu?Fu?DFT_INVERSE?+?DFT_SCALE）;
			split（Fu?temp）;
			u?=?temp[0];
	/******************************/
			iter?=?iter?+?1;
		}
		cout?<		outIter?=?outIter?+?1;
		beta?=?beta*beta_rate;
	}

	displayImage（g?“Blurr?Image“）;
	displayImage（u?“Deblurr?Image“）;
	waitKey（0）;
	system（“pause“）;
	return?0;
}

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2016-05-18?22:00??FTVD\
?????目錄???????????0??2016-05-17?16:15??FTVD\Debug\
?????文件??????127488??2016-05-18?21:54??FTVD\Debug\FTVD.exe
?????文件?????2264600??2016-05-18?21:54??FTVD\Debug\FTVD.ilk
?????文件?????2223104??2016-05-18?21:54??FTVD\Debug\FTVD.pdb
?????目錄???????????0??2016-05-18?21:54??FTVD\FTVD\
?????目錄???????????0??2016-05-18?21:54??FTVD\FTVD\Debug\
?????文件????????2509??2016-05-18?21:54??FTVD\FTVD\Debug\FTVD.log
?????目錄???????????0??2016-05-18?21:54??FTVD\FTVD\Debug\FTVD.tlog\
?????文件???????26586??2016-05-18?21:54??FTVD\FTVD\Debug\FTVD.tlog\CL.read.1.tlog
?????文件????????1314??2016-05-18?21:54??FTVD\FTVD\Debug\FTVD.tlog\CL.write.1.tlog
?????文件?????????158??2016-05-18?21:54??FTVD\FTVD\Debug\FTVD.tlog\FTVD.lastbuildstate
?????文件????????1210??2016-05-18?21:54??FTVD\FTVD\Debug\FTVD.tlog\cl.command.1.tlog
?????文件????????4176??2016-05-18?21:54??FTVD\FTVD\Debug\FTVD.tlog\link.command.1.tlog
?????文件????????5296??2016-05-18?21:54??FTVD\FTVD\Debug\FTVD.tlog\link.read.1.tlog
?????文件?????????446??2016-05-18?21:54??FTVD\FTVD\Debug\FTVD.tlog\link.write.1.tlog
?????文件??????307260??2016-05-18?21:54??FTVD\FTVD\Debug\main.obj
?????文件??????301221??2016-05-17?16:34??FTVD\FTVD\Debug\subfunction.obj
?????文件??????658432??2016-05-18?21:54??FTVD\FTVD\Debug\vc120.idb
?????文件??????946176??2016-05-18?21:54??FTVD\FTVD\Debug\vc120.pdb
?????文件????????4376??2016-05-17?15:37??FTVD\FTVD\FTVD.vcxproj
?????文件????????1164??2016-05-17?15:37??FTVD\FTVD\FTVD.vcxproj.filters
?????文件????????2294??2016-05-18?21:54??FTVD\FTVD\main.cpp
?????文件????????5296??2016-05-17?16:34??FTVD\FTVD\subfunction.cpp
?????文件?????????542??2016-05-17?16:33??FTVD\FTVD\subfunction.h
?????文件????15269888??2016-05-18?22:00??FTVD\FTVD.sdf
?????文件?????????958??2016-05-17?14:35??FTVD\FTVD.sln
?????文件???????21504??2016-05-18?22:00??FTVD\FTVD.v12.suo