【正文】 ma1Bvalue()。 selfUIviSigma2Cvalue()。 selfUIviSigma2Dvalue()。}PEHCorrection::~PEHCorrection(){}bool PEHCorrection::OnbtnOpenFile(Fl_Button *o,PEHCorrectionUI *v){ PEHCorrection * self = (PEHCorrection*)(vwnd)。 char * fileName。 BOOL isOK= FALSE。 while(! isOK) { fileName = fl_file_chooser(Open the 3D file, *.gz,)。 if(!fileName) return false。 if (_access(fileName,0)==1) { if (! fl_ask(The file does not exit, to select another one?)) return false。 } else { isOK = TRUE。 } } strcpy(selffileName, fileName)。 vtxtFileNamevalue(fileName)。 return true。}void PEHCorrection::OnbtnGo(Fl_Button *o,PEHCorrectionUI *v){ PEHCorrection * self = (PEHCorrection*)(vwnd)。 selfsigma1A = vviSigma1Avalue()。 selfsigma1B = vviSigma1Bvalue()。 selfsigma2C = vviSigma2Cvalue()。 selfsigma2D = vviSigma2Dvalue()。 double sigma1A(selfsigma1A),sigma1B(selfsigma1B),sigma2C(selfsigma2C),sigma2D(selfsigma2D)。 Datashort3D vqData,eqData。 ct_information ct_info。 int temp。 //clock_t StartTime, EndTime。 time_t StartTime, EndTime。 time(amp。StartTime)。 if ( !read_ct( selffileName, vqData, ct_info ) ) { printf( Error in read_ct()\n )。 return。 } int tq = 2100。 int max = 70。 int xSize = () , ySize = () , zSize = () 。 //zSize = max 。 (xSize,ySize,zSize)。pragma region to calculate eq for(int i=0。i xSize。i++) { for(int j=0。 jySize。j++) { for(int k=0。kzSize。k++) { if (vqData(i,j,k) tq) eqData(i,j,k) = vqData(i,j,k) tq。 else eqData(i,j,k) = 0。 } } }pragma endregion mist::array3short rData,rn_1Data,rnData,hatVqData。 (xSize,ySize,zSize)。 (0)。 (xSize,ySize,zSize)。 (0)。 (xSize,ySize,zSize)。 (0)。 (xSize,ySize,zSize)。pragma region To calculate r0 for(int k=0。kzSize。k++)//for every slice { for(int i=0。i xSize。i++) { for(int j=0。 jySize。j++) { if (eqData(i,j,k)!=0 ) { double sigma1 = sigma1A * (vqData(i,j,k)2000) + sigma1B。 double range = ceil(1+2*sigma1)。 //cout sigma1\n。 int iimin=(i range0) ? 0 : i range。 int iimax=(i + range xSize) ? xSize : i + range。 int jjmin=(j range0) ? 0 : j range。 int jjmax=(j + range ySize) ? ySize : j + range。 for(int ii=iimin。ii iimax 。ii++) { for(int jj=jjmin。 jjjjmax。jj++) { int kk = k。 if( ((iii)||(jjj)||(kkk))) { if (Math::distance(ii,jj,kk,i,j,k) range) rData(ii,jj,kk)+=eqData(i,j,k)*Math::gaussian(sigma1,i,j,ii,jj)。 } } } } } } }pragma endregion rn_1Data = rData。pragma region to calculate rn by iteration until energy bees negligible. bool isNegligible 。 int iteration = 0。 do { char fileName[128]。 sprintf(fileName,r%04dData,iteration)。 write_ct( fileName, rn_1Data , ct_info )。 isNegligible = true。 for(int k=0。kzSize。k++)//for every slice { for(int i=0。i xSize。i++) { for(int j=0。 jySize。j++) { if (rn_1Data(i,j,k)!=0 ) { double sigma2 = sigma2C * (rn_1Data(i,j,k)) + sigma2D。 double range = ceil(1+2*sigma2)。 //cout sigma1\n。 int iimin=(i range0) ? 0 : i range。 int iimax=(i + range xSize) ? xSize : i + range。 int jjmin=(j range0) ? 0 : j range。 int jjmax=(j + range ySize) ? ySize : j + range。 for(int ii=iimin。ii iimax 。ii++) { for(int jj=jjmin。 jjjjmax。jj++) { int kk = k。 if( ((iii)||(jjj)||(kkk))) { if (Math::distance(ii,jj,kk,i,j,k) range) { doubled istribute=rn_1Data(i,j,k)*Math::gaussian(sigma2,i,j,ii,jj)。 if(distribute ) isNegligible = false。 rData(ii,jj,kk)+=distribute。 rnData(ii,jj,kk)+=distribute。 } } } } } } } } rn_1Data = rnData。 (0)。 iteration++。 }while(! isNegligible)。pragma endregionpragma region to calculate hatVq for(int k=0。kzSize。k++)//for every slice { for(int i=0。i xSize。i++) { for(int j=0。 jySize。j++) { hatVqData(i,j,k) = vqData(i,j,k) rData(i,j,k)。 } } } write_ct( hatVqData, hatVqData , ct_info )。pragma endregion time(amp。EndTime)。 // printf(\nCalculation time is %f for ADC \n, difftime(EndTime, StartTime))。 //EndTime = clock()。 //printf(Time = % sec\n, (double)(EndTime StartTime )/CLOCKS_PER_SEC)。 write_ct( rData, rData , ct_info )。 write_ct( eqData, eqData , ct_info )。 write_ct( hatVqData, hatVqData , ct_info )。}附錄21．PLUTO瀏覽器的界面附錄211 PLUTO瀏覽界面jySIZE