vol.abserror: Absolute error volume

In espadon: Easy Study of Patient DICOM Data in Oncology

Absolute error volume

Description

The function vol.abserror provides the absolute error volume between 2 volumes.

Usage

vol.abserror(vol, vol.ref, T.MAT = NULL, alias = "", description = NULL)

Arguments

vol, vol.ref	"volume" class objects.
T.MAT	"t.mat" class object to link the reference frames of vol and vol.ref. T.MAT can be created by load.patient.from.Rdcm or load.T.MAT. If T.MAT = NULL, vol$ref.pseudo must be equal to vol.ref$ref.pseudo.
alias	Character string, $alias of the created object.
description	Character string, describing the created object. If description = NULL (default value), it will be set to "absolute error".

Value

Returns a "volume" class object (see espadon.class for class definitions), with the same grid as vol.ref, and representing the absolute error between vol and vol.ref.

See Also

vol.error

Examples

# loading of toy-patient objects (decrease dxyz)
step <- 5
patient <- toy.load.patient (modality = c("ct", "sct","rtstruct"), 
                             roi.name = c("eye", "brain","gizzard"),
                             dxyz = rep (step, 3))

patient$ct[[1]]$description
patient$ct[[2]]$description
# Creation of the absolute error volume between ct and synthetic ct
vAE <- vol.abserror (patient$ct[[2]], patient$ct[[1]], T.MAT = patient$T.MAT)

# Display
palette_vAE <- colorRampPalette(c("#00005F", "#0000FF", "#00FFFF", "#00FF00",
                                  "#FFFF00", "#FF7F00", "#FF0000", "#7F0000",
                                  "#5F0000")) (100)
breaks_vAE <- seq(floor (vAE$min.pixel), ceiling (vAE$max.pixel), 
                  length.out = 101)
                  
layout (mat = matrix(c(rep(1,6),2,2), ncol=4))                                
plot (vAE, view.coord =61, view.type = "trans", 
      col = palette_vAE, breaks = breaks_vAE)
display.palette(palette_vAE, breaks = breaks_vAE, 
                cex.axis = 1.2, main = vAE$unit)
par(mfrow=c(1,1))

espadon documentation built on April 11, 2025, 5:57 p.m.