# NOT RUN {
fi <- antsImageRead(getANTsRData("r16") ) %>%
resampleImage(c(60,60),1,0)
mi <- antsImageRead(getANTsRData("r64") ) %>%
resampleImage(c(60,60),1,0)
mytx <- antsRegistration(fixed=fi, moving=mi, typeofTransform = c('SyN') )
mcimg = antsImageRead( mytx$fwd[1] )
msk = getMask( fi )
vv=multichanneltovector(mcimg,msk)
mcimg2=vectortomultichannel( vv, msk )
vv2=multichanneltovector(mcimg2,msk)
cor.test(vv2,vv)
stopifnot( all( mcimg2[30,30] == mcimg[30,30] ) )
# }
Run the code above in your browser using DataLab