## pixelData = TRUE
## The DICOM image data are read from dicomSeparate()
dcmList <- dicomSeparate(system.file("hk-40", package="oro.dicom"))
sliceLocation <- NULL
dcmImage <- create3D(dcmList, mode="integer")
graphics::image(dcmImage[,,1], col=grey(0:64/64), axes=FALSE,
xlab="", ylab="")
dSL <- abs(diff(sliceLocation))
plot(dSL, ylim=range(range(dSL) * 1.5, 0, 10), xlab="Index", ylab="mm",
main="Difference in Slice Location")
## pixelData = FALSE
## The DICOM image data are read from create3D()
## This may save on memory for large batches of DICOM data
dcmList <- dicomSeparate(system.file("hk-40", package="oro.dicom"),
pixelData=FALSE)
dcmImage <- create3D(dcmList, mode="integer", pixelData=FALSE)
graphics::image(dcmImage[,,1], col=grey(0:64/64), axes=FALSE,
xlab="", ylab="")Run the code above in your browser using DataLab