#include
#include???

using?namespace?cv;
using?namespace?std;
float?correct（float?*sinogram?int?width?int?angle?float?SD?int?rotate）;//修正函數(shù)，包括正弦圖，寬度，角度（高度），源到探測器距離，旋轉中心
float?convfiltering（float?*sinogram?int?width?int?angle?int?start?int?end）;//卷積濾波函數(shù)，包括正弦圖，寬度，角度，開始與結束
float?inprojection（float?*sinogram?int?width?int?angle?float?anglestep?float?SD?int?rotate?float?*ROI?int?xmin?int?xmax?int?ymin?int?ymax?int?rotateangle）;//加權反投影函數(shù)，包括正弦圖，寬度，角度，步進角，源到探測器距離，旋轉中心，重建結果，ROI區(qū)域x與y最小最大坐標，旋轉角度

void?main（）
{
	Mat?angleimage?=?Mat（360?1920CV_32FC1Scalar（0.0））;
	Mat?result（401401?CV_32FC1）;
	Mat?picture;
	char?filename[256];
	double?start_t?=（double）cvGetTickCount（）;
	double?run_time;
	for（int?j?=500;?j<600;?j?=?j+100）????//重建的層號
	{
		for?（int?i?=?0;?i?		{
			sprintf_s（filename?“D:\\根系實驗數(shù)據(jù)\\CTyushiyan\\2018.3.8-1youji\\%d.png“i）;
????????????picture?=?imread（filename?IMREAD_ANYDEPTH）;//圖片必須添加到工程目錄
			Mat（picture.col（j）.t（））.copyTo（angleimage.row（i））;
		}

		correct（（float?*）angleimage.data?1920?360?4750.0730）;
		convfiltering（（float?*）angleimage.data?1920?36001900）;
		inprojection（（float?*）angleimage.data?1920?360?1.0?4750.0730?（float?*）result.data?-200200?-200?200?0）;
		normalize（result?result?1.0?0.0?NORM_MINMAX）;
		result?=?result?*?255;
		sprintf_s（filename?“D:\\根系實驗數(shù)據(jù)\\CTyushiyan\\有機斷面\\%d.png“j）;
		imwrite（filename?result）;
	}

	waitKey（0）;
}

float?correct（float?*sinogram?int?width?int?angle?float?SD?int?rotate）//修正函數(shù)，包括正弦圖，寬度，角度（高度），源到探測器距離，旋轉中心
{
	clock_t?t_start?t_end;//計時函數(shù)
	t_start?=?clock（）;

	int?i?j;
	float?k;
	for?（i?=?0;?i?	{
		k?=?SD?/?sqrt（SD*SD?+?（i?-?rotate）*（i?-?rotate））;//修正表達式
		for?（j?=?0;?j?		{
			sinogram[j*width?+?i]?*=?k;//將原正弦圖乘以修正表達式，得到新的正弦圖
		}
	}

	t_end?=?clock（）;
	float?time?=?（float）（t_end?-?t_start）?*?1000?/?CLOCKS_PER_SEC;//計時函數(shù)，計算耗時
	return?time;
}


float?convfiltering（float?*sinogram?int?width?int?angle?int?start?int?end）//卷積濾波函數(shù)，包括正弦圖，寬度，角度，開始與結束
{
	clock_t?t_start?t_end;//計時函數(shù)
	t_start?=?clock（）;

	if?（start?>=?end?||?start?=?width）//糾錯，防止出現(xiàn)錯誤
	{
		start?=?0;
		end?=?width?-?1;
	}
	float?*b?=?new?float[2?*?width?-?1];
	float?*h?=?new?float[width];
	float?*dst?=?new?float[width];
	int?i?j;
	float?a;
	for?（i?=?0;?i?	const?float?w?=?-0.20264236728f;//S-L濾波常數(shù)
	for?（i?=?0;?i?	{
		a?=?w?/?（4?*?i?*?i?-?1）;//S-L濾波函數(shù)
		b[width?-?1?-?i]?=?a;
		b[width?-?1?+?i]?=?a;
	}
	for?（j?=?0;?j?	{
		for?（i?=?0;?i?		{
			h[i]?=?sinogram[j*width?+?i];//將修正后的數(shù)據(jù)傳遞給需要進行卷積運算
		}
		int?minradius?=?（2?*?width?-?1）?/?2;
		int?maxradius?=?2?*?width?-?1?-?minradius;
		int?d?e?f;
		for?（d?=?start;?d?<=?end;?d++）
		{
			dst[d]?=?0;
			for?（e?=?-maxradius?+?1;?e?<=?minradius;?e++）
			{
				f?=?d?+?e;
				if?（f?>=?0?&&?f?					dst

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-03-18?21:13??FBPtest1\
?????目錄???????????0??2018-03-18?19:13??FBPtest1\Debug\
?????目錄???????????0??2018-03-18?21:13??FBPtest1\FBPtest1\
?????文件????26869760??2018-03-18?21:13??FBPtest1\FBPtest1.sdf
?????文件????????1397??2018-03-18?19:14??FBPtest1\FBPtest1.sln
?????文件???????33792??2018-03-18?21:13??FBPtest1\FBPtest1.v12.suo
?????目錄???????????0??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\
?????文件??????461620??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\fbp.obj
?????文件????????1534??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.log
?????目錄???????????0??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.tlog\
?????文件?????????612??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.tlog\cl.command.1.tlog
?????文件???????31936??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.tlog\CL.read.1.tlog
?????文件?????????408??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.tlog\CL.write.1.tlog
?????文件?????????160??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.tlog\FBPtest1.lastbuildstate
?????文件???????????2??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.tlog\link.command.1.tlog
?????文件???????????2??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.tlog\link.read.1.tlog
?????文件???????????2??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.tlog\link.write.1.tlog
?????文件???????????0??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\FBPtest1.tlog\unsuccessfulbuild
?????文件??????879616??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\vc120.idb
?????文件?????1970176??2018-03-18?19:13??FBPtest1\FBPtest1\Debug\vc120.pdb
?????文件????????5082??2018-03-18?21:13??FBPtest1\FBPtest1\fbp.cpp
?????文件????????7416??2018-03-18?19:13??FBPtest1\FBPtest1\FBPtest1.vcxproj
?????文件?????????944??2017-08-28?10:14??FBPtest1\FBPtest1\FBPtest1.vcxproj.filters
?????文件?????????165??2017-08-28?10:18??FBPtest1\FBPtest1\FBPtest1.vcxproj.user
?????目錄???????????0??2017-08-28?10:14??FBPtest1\FBPtest1\x64\
?????目錄???????????0??2018-03-18?19:15??FBPtest1\FBPtest1\x64\Debug\
?????文件??????579722??2018-03-18?19:15??FBPtest1\FBPtest1\x64\Debug\fbp.obj
?????文件?????????792??2018-03-18?19:15??FBPtest1\FBPtest1\x64\Debug\FBPtest1.Build.CppClean.log
?????文件????????1845??2018-03-18?19:15??FBPtest1\FBPtest1\x64\Debug\FBPtest1.log
?????目錄???????????0??2018-03-18?19:15??FBPtest1\FBPtest1\x64\Debug\FBPtest1.tlog\
?????文件?????????598??2018-03-18?19:15??FBPtest1\FBPtest1\x64\Debug\FBPtest1.tlog\cl.command.1.tlog
............此處省略13個文件信息