# NOT RUN {
data(powers_areal_mni_itk)
# }
# NOT RUN {
fixed <- antsImageRead( getANTsRData("ch2") )
moving <- antsImageRead( getANTsRData("mni") )
mytx <- antsRegistration( fixed=fixed, moving=moving, typeofTransform = c("SyN") )
data( "powers_areal_mni_itk", package = "ANTsR", envir = environment() )
coords = powers_areal_mni_itk[,1:3]
ch2reg = antsRegistration( fixed, moving, typeofTransform = "SyN" )
coordsw <- antsApplyTransformsToPoints( dim=3, points=coords,
transformlist=ch2reg$fwdtransforms, whichtoinvert=c(FALSE,FALSE) )
ptrd = 3
powersLabels = makePointsImage( coordsw, moving, radius = pard )
plot( moving, powersLabels, axis=3, nslices=30 )
# }
Run the code above in your browser using DataLab