ants_compcorr_group <- function( moco_img = "" , bm_img = "" , cortmask_img = "" , mocoparams_csv = "" , dosvd = 0 )
{
if( nchar( moco_img ) == 0 )
{
print( " please set a valid moco image " ) ;
print( " usage : ants_compcorr_group( <moco-image> , <brain-mask-image> , <cortmask-image> , <mocoparams-csv> , <dosvd> ) " ) ;
return ;
}
if( nchar( bm_img ) == 0 )
{
print( " please set a valid brain-mask image " ) ;
print( " usage : ants_compcorr_group( <moco-image> , <brain-mask-image> , <cortmask-image> , <mocoparams-csv> , <dosvd> ) " ) ;
return ;
}
if( nchar( cortmask_img ) == 0 )
{
print( " please set a valid cortmask image " ) ;
print( " usage : ants_compcorr_group( <moco-image> , <brain-mask-image> , <cortmask-image> , <mocoparams-csv> , <dosvd> ) " ) ;
return ;
}
if( nchar( mocoparams_csv ) == 0 )
{
print( " please set a valid Moco-params csv file " ) ;
print( " usage : ants_compcorr_group( <moco-image> , <brain-mask-image> , <cortmask-image> , <mocoparams-csv> , <dosvd> ) " ) ;
return ;
}
# split the string into filename and extension
split_img = strsplit( moco_img , "." , fixed = TRUE )[[1]] ;
filename = split_img[1] ;
if( length( split_img ) == 2 )
{
extension = paste( "" , split_img[2] , sep = "." ) ;
}
else if( length( split_img ) == 3 )
{
extension = paste( "" , split_img[2] , split_img[3] , sep = "." ) ;
}
compcorr_img = paste( filename , "_compcorr" , extension , sep = "" ) ;
ImageMath( 4 , compcorr_img , "CompCorrAuto" , moco_img , bm_img , 3 ) ;
csv = paste( filename , ".csv" , sep = "" ) ;
sccan( "--timeseriesimage-to-matrix" , paste( "[" , paste( moco_img , cortmask_img , 3.0 , 0.0 , sep = "," ) , "]" , sep = "" ) , "-o" , csv ) ;
compcorr_csv = paste( filename , "_compcorr_compcorr" , ".csv" , sep = "" ) ;
filt_csv = paste( filename , "_filt" , ".csv" , sep = "" ) ;
antsr_frequency_filter( csv , filt_csv , 3.0 , 0.1 , 0.25 , compcorr_csv ) ;
if( dosvd == 1 )
{
RSF_Networks_img = paste( filename , "_RSF_Networks" , extension , sep = "" ) ;
sccan( "--sparse-svd" , paste( "[" , paste( filt_csv , cortmask_img , -0.05 , sep = "," ) , "]" , sep = "" ) , "-n" , 40 , "-i" , 40 , "--PClusterThresh" , 100 , "-o" , RSF_Networks_img ) ;
ea_rsf = paste( filename , "_ea_rsf" , sep = "" ) ;
RSF_NetworksprojectionsView1vec_csv = paste( filename , "_RSF_NetworksprojectionsView1vec" , ".csv" , sep = "" ) ;
antsr_resting_state_corr_eigenanat( ea_rsf , RSF_NetworksprojectionsView1vec_csv , RSF_NetworksprojectionsView1vec_csv , compcorr_csv, mocoparams_csv ) ;
}
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.