GLT.calculateDVHs: Returns the DVH calculated from data stored into a geoLet...

In kbolab/moddicom: Package for handling DICOM RT files

Description

This function calculates DVH starting from DICOM RT-DOSE objects, previously loaded into a geoLet object.

N.B: if you are reading this help, your life is going to end in 10, 9, 8, 7, ...

Usage

GLT.calculateDVHs(obj.geoLet, ROIName, newPixelSpacing = NA,
  justTheDVH = TRUE, verbose = FALSE, forceReCalculus = FALSE,
  fastEngine = TRUE, decimation = FALSE, decimation.percentage = 0.8,
  smoothing = FALSE, smoothing.iterations = 10)

Arguments

obj.geoLet	the object geoLet to load the DICOM serie in
ROIName	the name of the ROI you are interested in having the DVH
newPixelSpacing	Optional. By default it uses same pixelValues of the main CT scan but if you want to interpolate the space in order to have a bigger detail you can increase the precision. Pay attention: it can take a lot o memory!!!
justTheDVH	Optional. TRUE by default. Because of this function is also able to calculate other ancillary data structures, if you want you can set this parameter to FALSE and get also the IMAGE voxel cube and te DOSE voxel cube (sampled at the same points in the 3D space)
verbose	Optional. FALSE by default. Set to TRUE it allow to see some logs during the computation
forceReCalculus	Optional. FALSE by default. Due to a caching system, implemented to improve performances, If a previous calculus has been interrupted the afterwards computations (on the same ROI) could have problems. If a computation has been interrupted it is a good practice to set this parameter to TRUE in the next computation in order to let the algorithm to re-build anything.
fastEngine	TRUE by default. By default it uses the vcgClostKD() function which should be quicker but probably a bit less accurate (?). If you want to try something different you can set this parameter on FALSE and it will use the most classical vcgClost
decimation	Optional. FALSE by default. Set it to TRUE if you want to enable decimation to the mesh structure
smoothing	Optional, FALSE by default. This parameter, if se to TRUE allow to smooth the mesh by a numerical factor indicated in the parameter smoothing.iterations.
decimationpercentage.	Optional. 0.8 by default. If the parameter decimation is set to TRUE this parameter indicates the percentage of triangles that should be sacrified in the first (and biggest) mesh-model. Killing triangles allows to have an easier to handle mesh and less memory-consumer but reduce the quality of the approximation. In most cases the optimal value has to be defined empirically. If decimation is set to FALSE dont's waste your time in tuning this parameter: decimation will not be performed!
smoothing.terations	Optional, 10by default. Setting this parameter makes sense only if smoothing was previously set to TRUE. Increasing the number of iterations we can improve the "smoothing" applied to the original mesh.

Value

Depending on the status of the parameter justTheDVH. If it is set to TRUE the function returns a dvhmatrix object. If set to FALSE it returns a list of three elements:

• voxelCube.CT this is a 3D array which contains the voxel cube of the original CT scan

• voxelCube.Dose this is a 3D array which contains the voxel cube of the DOSES. The dimension is the same of the voxelCube because the dose value in the space has been interpolated (via triliear interpolation) in order to overlap the centroids of the voxels between voxelCube.CT and voxelCube.Dose. This should allow to quickly plot the expected dose value on a specific pixel/voxel in the CT scan.

• DVHobj a dvhmatrix object

Examples

## Not run: 

obj<-geoLet()
obj$openDICOMFolder(pathToOpen='./DICOMSeries/pat001' );
a<-obj$GLT.calculateDVH( Structure = "GTV" )


## End(Not run)

kbolab/moddicom documentation built on Nov. 29, 2020, 9:11 p.m.