# NOT RUN {
img <- antsImageRead( getANTsRData("r16") , 2 )
labels = getMask(img,1,150) + getMask(img,151,225) * 2
labelsDilated = multiLabelMorphology(labels, 'MD', 2)
# should see original label regions preserved in dilated version
# label N should have mean N and 0 variance
labelStats(labelsDilated, labels)
# }
Run the code above in your browser using DataLab