# NOT RUN {
fi = antsImageRead( getANTsRData( "r16" ) )
seg = kmeansSegmentation( fi, 3 )
msk = thresholdImage( seg$segmentation, 1, Inf )
pseg = priorBasedSegmentation( fi, seg$probabilityimages, msk, 0.25, 0.1, 3 )
# }
Run the code above in your browser using DataLab