# NOT RUN {
# }
# NOT RUN {
asl<-antsImageRead("PEDS012_20131101_pcasl_1.nii.gz",4)
# image available at http://files.figshare.com/1701182/PEDS012_20131101.zip
pcasl.processing <- aslPerfusion( asl, moreaccurate=T, interpolation="linear",
dorobust=0.95, useDenoiser=1:6 )
# user might compare to useDenoiser=FALSE
pcasl.parameters <- list( sequence="pcasl", m0=pcasl.processing$m0 )
cbf <- quantifyCBF( pcasl.processing$perfusion, pcasl.processing$mask, pcasl.parameters )
meancbf <- cbf$kmeancbf
print(mean(meancbf[ pcasl.processing$mask==1 ]))
antsImageWrite( meancbf , 'temp.nii.gz')
plotANTsImage( meancbf, slices='1x50x1')
# }
Run the code above in your browser using DataLab