Hi Dr. Yeh,
When I tried to correct and reconstruct the data with the command line, I got an unknown error. Do you know what's going on? To clarify my problem, I attached my command line and the output below. Thank you so much!
Sincerely,
Xitong
Command line I used:
subjects=(20220329_F019Y_10172 20220406_F019Y_10174 20220405_F024Y_10175 20220414_F022Y_10176 20220513_F020Y_10179)
for sub in ${subjects[@]}
do
input=/Shared/lss_kahwang_hpc/ThalHi_Diffusion/SRC_folder/${sub}_s06_HCPDTI.src.gz
reversed=/Shared/lss_kahwang_hpc/ThalHi_Diffusion/SRC_folder/${sub}_s05_DWIB0.nii.gz
singularity exec /Shared/lss_kahwang_hpc/opt/dsi_studio_latest/dsistudio_latest.sif dsi_studio \
--action=rec --source=${input} --rev_pe=${reversed} --method=7 --param0=1.25 --check_btable=1 \
--output=/Shared/lss_kahwang_hpc/ThalHi_Diffusion/Corrected_folder/
done
The output of one subject. I got the same error for all of them.
Running on compute node: argon-itf-bx47-40.hpc.
In directory: /Shared/lss_kahwang_hpc/ThalHi_Diffusion
Starting on: Sun May 22 11:47:37 CDT 2022
DSI Studio "Chen" May 15 2022
source=/Shared/lss_kahwang_hpc/ThalHi_Diffusion/SRC_folder/20220329_F019Y_10172_s06_HCPDTI.src.gz
action=rec
loading source...
prepare index file for future accelerated loading
SRC file loaded
src loaded
rev_pe=/Shared/lss_kahwang_hpc/ThalHi_Diffusion/SRC_folder/20220329_F019Y_10172_s05_DWIB0.nii.gz
topup/eddy
projected correction:0.938253 0.983177 0
source com:57.4211 55.5789 23.2632
rev pe com:65.8333 30.6111 16.6111
source and reverse phase encoding: AP_PA
create acq params at 20220329_F019Y_10172_s06_HCPDTI.topup.acqparams.txt
0 -1 0 0.05
0 1 0 0.05
get the bounding box for speeding up topup/eddy
create topup needed b0 nii.gz file from AP_PA b0
FSL installation found at /opt/fsl
run /opt/fsl/bin/topup at /Shared/lss_kahwang_hpc/ThalHi_Diffusion/SRC_folder with
--warpres=20,16,14,12,10,6,4,4,4
--subsamp=2,2,2,2,2,1,1,1,1
--fwhm=8,6,4,3,3,2,1,0,0
--miter=5,5,5,5,5,10,10,20,20
--lambda=0.005,0.001,0.0001,0.000015,0.000005,0.0000005,0.00000005,0.0000000005,0.00000000001
--estmov=1,1,1,1,1,0,0,0,0
--minmet=0,0,0,0,0,1,1,1,1
--scale=1
--imain=20220329_F019Y_10172_s06_HCPDTI.topup.AP_PA.nii.gz
--datain=20220329_F019Y_10172_s06_HCPDTI.topup.acqparams.txt
--out=20220329_F019Y_10172_s06_HCPDTI
--iout=20220329_F019Y_10172_s06_HCPDTI.topup.check_result
--verbose=1
ERROR:Unknown error
Warning: --method is not used. Please check command line syntax.
Warning: --param0 is not used. Please check command line syntax.
Warning: --check_btable is not used. Please check command line syntax.
Warning: --output is not used. Please check command line syntax.