img <- makeImage(c(4,4), rnorm(4*4))
mask <- makeImage(c(4,4),
as.matrix(c(0,0,0,0,
0,1,1,0,
0,1,1,0,
0,0,0,0), nrow=4))
plot(img, list(mask))
## Not run:
# mnit<-getANTsRData('mni')
# mnit<-antsImageRead(mnit)
# mniafn<-getANTsRData('mnia')
# mnia<-antsImageRead(mniafn)
# mnia<-thresholdImage(mnia,22,25)
# mnia<-smoothImage(mnia,1.5)
# mnia2<-antsImageRead(mniafn)
# mnia2<-thresholdImage(mnia2,1,4)
# mnia2<-smoothImage(mnia2,1.5)
# ofn<-paste(tempfile(),'.png',sep='')
# # write directly to a file
# plot( mnit, list(mnia,mnia2), slices=seq(50, 140, by=5),
# window.overlay = c(0.25,1), axis=2,
# color.overlay=c('red','blue'), outname = ofn )
# ## End(Not run)
Run the code above in your browser using DataLab