## pixelData = TRUE
dcmList <- dicomSeparate(system.file("hk-40", package="oro.dicom"))
dcmImage <- create3D(dcmList, mode="integer")
## Accumulate voxel dimensions
pS <- header2matrix(extractHeader(dcmList$hdr, "PixelSpacing", FALSE), 2)[1,]
sT <- extractHeader(dcmList$hdr, "SliceThickness")[1]
pixdim <- c(0, pS, sT, rep(1,4))
require("oro.nifti")
dcmNifti <- nifti(dcmImage, datatype=4, pixdim=pixdim)
par(bg="black")
image(dcmNifti)
orthographic(dcmNifti)Run the code above in your browser using DataLab