# NOT RUN {
img <- antsImageRead( getANTsRData( "r16" ) )
msk <- thresholdImage( img, quantile( img[ img > 0 ] )[1], max( img ) )
ripped <- ripmmarc( img, msk, patchRadius=3, patchSamples=2000, patchVarEx=5)
bimg = ripmmarcBasisImage( ripped$canonicalFrame, ripped$basisMat[4,] )
# }
Run the code above in your browser using DataLab