CT图像重建

基于C++编写的CT图像重建例子，包含CT原始数据

// Plotting.cpp : Defines the entry point for the console application.
//

#include “stdafx.h“
#include 
#include 
#include 
#include 
#include 
#include 
#include 

using namespace std;

//the grey value 0~255 change to attenuation -4~2;
double grey2attenu（CvScalar newcolor）
{
	return （newcolor.val[0]/255*6-4）;
}

//the attenuation -4~2 change to grey value 0~255;
CvScalar attenu2grey（double attenu）
{
	double a;
	a=（attenu+4）/6*255;
	return cvScalar（a）;
}


int main（）
{
	const char* saveimg1=“S-L Head Model.bmp“;
	const char* saveimg1rotate1=“S-L Head Model Rotate.bmp“;
	const char* saveimg1proj1=“S-L Head Model Back Projection Reconstruction.bmp“;

	const int imagedepth=IPL_DEPTH_8U;

	int SNR;		//the signal noise rate

	srand（（unsigned）time（ NULL ））;		//set the random number seed system time
	
	int tijk;//loop indexes

	//the R-L convolution sequence （sharp frequence domain） 
	const unsigned int RLlen=800; 
	double rlconv[2*RLlen+1];
	for （i=0;i<2*RLlen+1;i++）
	{
		if （i==RLlen） rlconv[i]=0.25;
		if （（（i%2）==0）&&（i!=RLlen）） rlconv[i]=0.;
		if （（（i%2）!=0）&&（i!=RLlen）） rlconv[i]=-1.0/CV_PI/CV_PI/（（i-RLlen）*（i-RLlen））;
	}

	//the S-L convolution sequence （smooth frequence domain）
	const unsigned int SLlen=800;
	double slconv[2*SLlen+1];
	for （i=0;i<2*SLlen+1;i++） slconv[i]=-2./CV_PI/CV_PI/（4*（i-SLlen）*（i-SLlen）-1.）;

	//my convolutions with different frequency respond
	const unsigned int MYCONVLEN=60;
	const double convsharphighpass[2*MYCONVLEN+1]={-0.000240858140043-0.000133223717409-0.0000051541704450.0001369392936370.000284263415617
0.0004249274624080.0005437727031510.0006230332806720.0006439782802440.000589494970690
0.0004473598599180.000213744039480-0.000103649544388-0.000483483296816-0.000890339445913
-0.001276395732730-0.001585422616489-0.001758962282615-0.001744177246757-0.001502505408111
-0.001017986023455-0.0003039789510530.0005929741423560.0015931853246580.002588399771881
0.0034510289709450.0040475229611550.0042547304172400.0039775329690010.003165636034711
0.0018272240470750.000037291132652-0.002061145017234-0.004264068980025-0.006321971804775
-0.007963063485987-0.008922404125941-0.008974258625376-0.007964308634454-0.005838020030214
-0.0026615322953880.0013680679949710.0059253291905860.0105728860305390.014791359623433
0.0180207909559990.0197103054129780.0193715874560480.0166310359127320.011275272117329
0.003285033288548-0.007146610768674-0.019614402582766-0.033520369782713-0.048111703314857
-0.062534256227342-0.075897240654455-0.087343467603161-0.096118785174237-0.101634309392224
0.896484375000000-0.101634309392238-0.096118785174230-0.087343467603166-0.075897240654451
-0.062534256227345-0.048111703314854-0.033520369782715-0.019614402582764-0.007146610768675
0.0032850332885500.0112752721173270.0166310359127330.0193715874560480.01

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

     文件      30771  2007-11-20 21:08  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\CTconv_reconstuct.cpp

     文件   17104896  2014-03-06 15:12  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\CTconv_reconstuct.sdf

     文件        900  2011-07-19 17:33  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\CTconv_reconstuct.sln

    ..A..H.     12288  2011-07-19 17:33  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\CTconv_reconstuct.suo

     文件       4615  2007-11-17 14:52  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\CTconv_reconstuct.vcproj

     文件       1411  2007-11-20 21:12  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\CTconv_reconstuct.vcproj.SIMONXU.Simon Xu.user

     文件       6380  2014-03-06 15:09  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\CTconv_reconstuct.vcxproj

     文件       1234  2011-07-19 17:25  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\CTconv_reconstuct.vcxproj.filters

     文件        143  2011-07-19 17:25  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\CTconv_reconstuct.vcxproj.user

     文件    2359296  2011-07-19 17:27  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\ipch\ctconv_reconstuct-d9859939\ctconv_reconstuct-e53b8102.ipch

     文件       1366  2007-11-17 14:50  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\ReadMe.txt

     文件        304  2007-11-17 14:50  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\stdafx.cpp

     文件        376  2007-11-17 14:50  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\stdafx.h

     文件       2659  2011-07-19 17:25  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\UpgradeLog.xml

     文件       3348  2011-07-19 17:25  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\_UpgradeReport_Files\UpgradeReport.css

     文件      12579  2009-08-31 01:49  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\_UpgradeReport_Files\UpgradeReport.xslt

     文件         69  2011-07-19 17:25  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\_UpgradeReport_Files\UpgradeReport_Minus.gif

     文件         71  2011-07-19 17:25  C++CT重建\CT_Reconstruction\CT重建仿真\CTconv_reconstuct\_UpgradeReport_Files\UpgradeReport_Plus.gif

     文件     641078  2007-11-20 00:46  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=10） S-L Head Model R-L Convolution Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-20 00:44  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=10） S-L Head Model S-L Convolution Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-20 00:30  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=25） S-L Head Model Direct Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-20 20:56  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=25） S-L Head Model My Convolution（121P） Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-20 15:16  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=25） S-L Head Model My Convolution3（121P） Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-20 00:32  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=25） S-L Head Model R-L Convolution Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-20 00:34  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=25） S-L Head Model S-L Convolution Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-19 23:17  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=5） S-L Head Model Direct Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-20 20:59  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=5） S-L Head Model My Convolution（121P） Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-20 21:01  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=5） S-L Head Model R-L Convolution Back Projection_1degree_Scan.bmp

     文件     641078  2007-11-20 21:04  C++CT重建\CT_Reconstruction\CT重建仿真\Noisy（SNR=5） S-L Head Model S-L Convolution Back Projection_1degree_Scan.bmp

     文件        492  2008-02-25 18:13  C++CT重建\CT_Reconstruction\CT重建仿真\Readme.txt

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