# Create a clustered volume
mask <- read_vol(system.file("extdata", "global_mask_v4.nii", package="neuroim2"))
clusters <- rep(1:5, length.out=sum(mask))
cvol <- ClusteredNeuroVol(mask, clusters)
# Convert to DenseNeuroVol
dvol <- as(cvol, "DenseNeuroVol")
Run the code above in your browser using DataLab