# loading of toy-patient objects (decrease dxyz and increase beam.nb for
# better result)
step <- 5
patient <- toy.load.patient (modality = c("ct","mr", "rtdose"),
dxyz = rep (step, 3), beam.nb = 3)
MR <- patient$mr[[1]]
CT <- patient$ct[[1]]
D <- patient$rtdose[[1]]
# display
display.kplane (CT)
display.kplane (MR, k = floor (length(MR$k.idx)*5/8),
col = grey.colors (256, start = 0, end = 1),
breaks = seq (0, 500, length.out = 257), bg = "darkblue")
display.kplane (D, k = floor (length(D$k.idx)*3/8),
col = rainbow (256, s = seq (1, 0, length.out = 256),
start = 0, end = 4/6,
alpha = seq (0.8, 0, length.out=256),
rev = TRUE),
bg = "darkblue", ord.flip = TRUE, sat.transp = FALSE,
interpolate = FALSE)
display.kplane (CT, k = floor (length(CT$k.idx)/3), col = pal.RVV (1000),
bg = "darkblue", ord.flip = TRUE, interpolate = FALSE)
Run the code above in your browser using DataLab