84 Cerr(
"***** vDeformation::CheckParameters() -> No image dimensions provided !" << endl);
91 Cerr(
"***** vDeformation::CheckParameters() -> Wrong verbosity level provided !" << endl);
98 Cerr(
"***** vDeformation::CheckParameters() -> Number of transformations in the deformation has not been initialized !" << endl);
105 Cerr(
"***** vDeformation::CheckParameters() -> An error occurred while checking parameters of the child dynamic class !" << endl);
141 Cerr(
"***** vDeformation::ApplyDeformationsToBackwardImage() -> Called while not initialized !" << endl);
156 Cerr(
"***** vDeformation::ApplyDeformationsToBackwardImage()-> An error occurred while performing backward deformation of the backward image !" << endl);
157 Cerr(
"***** frame index " << a_fr <<
" respiratory image index " << rimg<<
" cardiac image index " << cimg<<
" !" << endl);
193 Cerr(
"***** vDeformation::ApplyDeformationsToHistoSensitivityImage() -> Called while not initialized !" << endl);
207 Cerr(
"***** vDeformation::ApplyDeformationsToHistoSensitivityImage()-> An error occurred while performing backward deformation of the backward image !" << endl);
208 Cerr(
"***** frame index " << a_fr <<
" respiratory image index " << rimg<<
" cardiac image index " << cimg<<
" !" << endl);
248 Cerr(
"***** vDeformation::PerformDeformation() -> Called while not initialized !" << endl);
270 Cerr(
"***** vDeformation::ApplyDeformations()-> An error occurred while performing backward deformation of the backward image !" << endl);
271 Cerr(
"***** frame index " << a_fr <<
" respiratory image index " << a_rimg<<
" cardiac image index " << a_cimg<<
" !" << endl);
293 Cerr(
"***** vDeformation::ApplyDeformations()-> An error occurred while performing forward deformation of the forward image !" << endl);
294 Cerr(
"***** frame index " << a_fr <<
" respiratory image index " << a_rimg<<
" cardiac image index " << a_cimg<<
" !" << endl);
328 Cerr(
"***** vDeformation::PerformHistoSensitivityDeformation() -> Called while not initialized !" << endl);
348 Cerr(
"***** vDeformation::PerformHistoSensitivityDeformation()-> An error occurred while performing backward deformation of the sensitivity image !" << endl);
349 Cerr(
"***** frame index " << fr <<
" respiratory image index " << rimg<<
" cardiac image index " << cimg<<
" !" << endl);
393 Cerr(
"***** vDeformation::PerformSensitivityDeformation() -> Called while not initialized !" << endl);
405 Cerr(
"***** vDeformation::PerformSensitivityDeformation -> An error occurred while performing forward deformation !" << endl);
406 Cerr(
"***** frame index " << fr <<
" respiratory gate index " << rg<<
" cardiac gate index " << cg<<
" !" << endl);
417 Cerr(
"***** vDeformation::PerformSensitivityDeformation -> An error occurred while performing backward deformation !" << endl);
418 Cerr(
"***** frame index " << fr <<
" respiratory gate index " << rg<<
" cardiac gate index " << cg<<
" !" << endl);
424 Cerr(
"***** vDeformation::PerformDeformation -> Unknown type of deformation !" << endl);
451 int incX = dX>0 ? 1 : -1;
455 dX = (dX>=0) ? dX : -dX;
456 dY = (dY>=0) ? dY : -dY;
457 dZ = (dZ>=0) ? dZ : -dZ;
461 if ((iiv < nb_tot_vox)
463 ap_outputImage[iov] += ap_inputImage[iiv] * (1.-dX)*(1.-dY)*(1.-dZ);
465 if ((iiv+incX < nb_tot_vox)
466 && (((int)iiv)+incX >= 0))
467 ap_outputImage[iov] += ap_inputImage[iiv+incX] * dX*(1.-dY)*(1.-dZ);
469 if ((iiv+incY < nb_tot_vox )
470 && (((
int)iiv)+incY >= 0))
471 ap_outputImage[iov] += ap_inputImage[iiv +incY] * (1.-dX)*dY*(1.-dZ);
473 if ((iiv+incX+incY < nb_tot_vox)
474 && ( ((
int)iiv)+incX+incY >= 0))
475 ap_outputImage[iov] += ap_inputImage[iiv+incX+incY] * dX*dY*(1.-dZ);
477 if ((iiv+incZ < nb_tot_vox)
478 && (((int)iiv)+incZ) >= 0)
479 ap_outputImage[iov] += ap_inputImage[iiv+incZ] * (1.-dX)*(1.-dY)*dZ;
481 if ((iiv+incX+incZ) < nb_tot_vox
482 && (((int)iiv)+incX+incZ) >=0 )
483 ap_outputImage[iov] += ap_inputImage[iiv+incX+incZ] * dX*(1.-dY)*dZ;
485 if ((iiv+incY+incZ < nb_tot_vox)
486 && (((int)iiv)+incY+incZ >= 0))
487 ap_outputImage[iov] += ap_inputImage[iiv+incY+incZ] * (1.-dX)*dY*dZ;
489 if ((iiv+incX+incY+incZ < nb_tot_vox)
490 && (((int)iiv)+incX+incY+incZ >= 0))
491 ap_outputImage[iov] += ap_inputImage[iiv+incX+incY+incZ] *dX*dY*dZ;
FLTNB **** m4p_forwardImage
Declaration of class oImageDimensionsAndQuantification.
FLTNB **** m4p_refDynForwardImage
FLTNB ***** m5p_refDynBackwardImage
FLTNB **** m4p_refDynSensitivityImage
FLTNB ****** m6p_backwardImage
int GetNbBackwardImages()
INTNB GetNbVoxXY()
Get the number of voxels in a slice.
Declaration of class oImageSpace.
int GetNbCardGates()
Get the number of cardiac gates.
This class holds all the matrices in the image domain that can be used in the algorithm: image...
INTNB GetNbVoxXYZ()
Get the total number of voxels.
int GetNbThreadsForProjection()
Get the number of threads used for projections.
INTNB GetNbVoxX()
Get the number of voxels along the X axis.
int GetNbRespGates()
Get the number of respiratory gates.
FLTNB ***** m5p_sensitivity
#define VERBOSE_DEBUG_NORMAL