1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56
|
#!/bin/bash
if [ ! -e "$BART_TOOLBOX_PATH"/bart ] ; then
if [ -e "$TOOLBOX_PATH"/bart ] ; then
BART_TOOLBOX_PATH="$TOOLBOX_PATH"
else
echo "\$BART_TOOLBOX_PATH is not set correctly!" >&2
exit 1
fi
fi
export PATH="$BART_TOOLBOX_PATH:$PATH"
# oversampled radial trajectory
bart traj -r -y55 -x256 traj_tmp
bart scale 0.5 traj_tmp traj
# simulate k-space
bart phantom -t traj ksp
# compute Ram-Lak filter
bart rss 1 traj ramlak
# apply to data
bart fmac ksp ramlak ksp_filt
# adjoint nufft
bart nufft -a traj ksp img
bart nufft -a traj ksp_filt img_filt
# grid and degrid ones
bart ones 3 1 256 55 ones
bart nufft -a traj ones dens_tmp
bart nufft traj dens_tmp density
# sqrt
bart spow -- -1. density dcf
# inv sqrt
bart spow -- -0.5 density sqdcf
# adjoint nufft
bart fmac dcf ksp ksp_filt2
bart nufft -a traj ksp_filt2 img_filt2
# one channel all ones sensititty
bart ones 3 256 256 1 sens
# without dcf
bart pics -i30 -t traj ksp sens img_pics_i30
bart pics -i3 -t traj ksp sens img_pics_i3
# with dcf
bart pics -i30 -t traj -p sqdcf ksp sens img_pics_dcf_i30
bart pics -i3 -t traj -p sqdcf ksp sens img_pics_dcf_i3
|