資源簡介
基于C++編寫的CT圖像重建例子,包含CT原始數(shù)據(jù)
代碼片段和文件信息
//?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
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件??????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.xm
?????文件???????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個(gè)文件信息
評論
共有 條評論