open-discussion
open-discussion > RE: slomoco with pestica4
Sep 22, 2016 09:09 PM | Christian Kaufmann
RE: slomoco with pestica4
Finally, after upgrading to AFNI 16.2.17, slicemoco_newalgorithm.sh
-d epiFile -r worked with my setup (Ubuntu 14.04 LTS, Matlab
R2012b). Thank you for your help so far!
However, I encountered another issue while running second order moco (whole error output below).
You might have an idea what's going wrong especially where it reads:
===
Error in convert_1dmat_into_tranmatarray (line 3)
[zdim tdim veclen] = size(transmat1d_zt);
Output argument "tranmat_zt" (and maybe others) not assigned during call to
"/mrt/Applications/pestica4/slomoco/convert_1dmat_into_tranmatarray.m>convert_1dmat_into_tranmatarray".
Error in slicemoco_newalgorithm_input (line 25)
tranmat_zt = convert_1dmat_into_tranmatarray(transmat1d_zt);% tranmat_zt =
reshape(tranmat_zt,[zmbdim tdim 3 4]); % (v1) debugged
===
Thank you!
Christian
Whole error output during "slicemoco_newalgorithm.sh -d epiFile -r:
===
finished slicewise in-plane motion correction
Exiting
Running Secondorder Motion Correction using SLOMOCO output
cp ../pestica4/card_pestica4.dat card_pestica4.dat
cp ../pestica4/resp_pestica4.dat resp_pestica4.dat
matlab -nojvm -nosplash addpath /mrt/Applications/pestica4/slomoco; addpath /mrt/Applications/pestica4/afni_matlab; addpath /mrt/Applications/pestica4/pestica_matlab; CARD=load('card_pestica4.dat'); RESP=load('resp_pestica4.dat'); slicemoco_newalgorithm_input('OutBrick_run_022.slicemocoxy_afni+orig','OutBrick_run_022.brain+orig','tempslmoco_volslc_alg_vol_OutBrick_run_022.slicemocoxy_afni/motion.wholevol_zt',CARD,RESP); exit;
< M A T L A B (R) >
Copyright 1984-2012 The MathWorks, Inc.
R2012b (8.0.0.783) 64-bit (glnxa64)
August 22, 2012
To get started, type one of these: helpwin, helpdesk, or demo.
For product information, visit http://www.mathworks.com.
>> Note: SMS acquisition is applied (MB acc. fac = 6).
Error in convert_1dmat_into_tranmatarray (line 3)
[zdim tdim veclen] = size(transmat1d_zt);
Output argument "tranmat_zt" (and maybe others) not assigned during call to
"/mrt/Applications/pestica4/slomoco/convert_1dmat_into_tranmatarray.m>convert_1dmat_into_tranmatarray".
Error in slicemoco_newalgorithm_input (line 25)
tranmat_zt = convert_1dmat_into_tranmatarray(transmat1d_zt);% tranmat_zt =
reshape(tranmat_zt,[zmbdim tdim 3 4]); % (v1) debugged
>> matlab -nojvm -nosplash addpath /mrt/Applications/pestica4/slomoco; addpath /mrt/Applications/pestica4/pestica_matlab; addpath /mrt/Applications/pestica4/afni_matlab; qa_slomoco('OutBrick_run_022.slicemocoxy_afni+orig','tempslmoco_volslc_alg_vol_OutBrick_run_022.slicemocoxy_afni/motion.wholevol_zt','tempslmocoxy_afni_OutBrick_run_022'); exit;
< M A T L A B (R) >
Copyright 1984-2012 The MathWorks, Inc.
R2012b (8.0.0.783) 64-bit (glnxa64)
August 22, 2012
To get started, type one of these: helpwin, helpdesk, or demo.
For product information, visit http://www.mathworks.com.
>> Note: SMS acquisition is applied (MB acc. fac = 6).
Index exceeds matrix dimensions.
Error in qa_slomoco (line 61)
outofplane(acqtp,:) = rotmat_zt(z,t,:); % [tdim*zmbdim [dL dP dS
rotmatvector (1:9)]]
>> *+ WARNING: If you are performing spatial transformations on an oblique dset,
such as ./OutBrick_run_022+orig.BRIK,
or viewing/combining it with volumes of differing obliquity,
you should consider running:
3dWarp -deoblique
on this and other oblique datasets in the same session.
See 3dWarp -help for details.
++ Oblique dataset:./OutBrick_run_022+orig.BRIK is 4.199786 degrees from plumb.
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as ./OutBrick_run_022+orig.BRIK,
or viewing/combining it with volumes of differing obliquity,
you should consider running:
3dWarp -deoblique
on this and other oblique datasets in the same session.
See 3dWarp -help for details.
++ Oblique dataset:./OutBrick_run_022+orig.BRIK is 4.199786 degrees from plumb.
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as ./OutBrick_run_022+orig.BRIK,
or viewing/combining it with volumes of differing obliquity,
you should consider running:
3dWarp -deoblique
on this and other oblique datasets in the same session.
See 3dWarp -help for details.
++ Oblique dataset:./OutBrick_run_022+orig.BRIK is 4.199786 degrees from plumb.
3drefit -saveatr -atrint TAXIS_NUMS 350 72 77002 -999 -999 -999 -999 -999 OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit: AFNI version=AFNI_16.2.17 (Sep 16 2016) [64-bit]
++ Authored by: RW Cox
** ERROR: Can't open dataset OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit processed 0 datasets
3drefit -saveatr -atrfloat TAXIS_FLOATS 0 1.05 0 -31.762 2 -999999 -999999 -999999 OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit: AFNI version=AFNI_16.2.17 (Sep 16 2016) [64-bit]
++ Authored by: RW Cox
** ERROR: Can't open dataset OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit processed 0 datasets
3drefit -saveatr -atrfloat TAXIS_OFFSETS 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit: AFNI version=AFNI_16.2.17 (Sep 16 2016) [64-bit]
++ Authored by: RW Cox
** ERROR: Can't open dataset OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit processed 0 datasets
3dNotes -h slicewise_moco_inplane.sh OutBrick_run_022+orig OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
Error: Cannot open dataset
===
However, I encountered another issue while running second order moco (whole error output below).
You might have an idea what's going wrong especially where it reads:
===
Error in convert_1dmat_into_tranmatarray (line 3)
[zdim tdim veclen] = size(transmat1d_zt);
Output argument "tranmat_zt" (and maybe others) not assigned during call to
"/mrt/Applications/pestica4/slomoco/convert_1dmat_into_tranmatarray.m>convert_1dmat_into_tranmatarray".
Error in slicemoco_newalgorithm_input (line 25)
tranmat_zt = convert_1dmat_into_tranmatarray(transmat1d_zt);% tranmat_zt =
reshape(tranmat_zt,[zmbdim tdim 3 4]); % (v1) debugged
===
Thank you!
Christian
Whole error output during "slicemoco_newalgorithm.sh -d epiFile -r:
===
finished slicewise in-plane motion correction
Exiting
Running Secondorder Motion Correction using SLOMOCO output
cp ../pestica4/card_pestica4.dat card_pestica4.dat
cp ../pestica4/resp_pestica4.dat resp_pestica4.dat
matlab -nojvm -nosplash addpath /mrt/Applications/pestica4/slomoco; addpath /mrt/Applications/pestica4/afni_matlab; addpath /mrt/Applications/pestica4/pestica_matlab; CARD=load('card_pestica4.dat'); RESP=load('resp_pestica4.dat'); slicemoco_newalgorithm_input('OutBrick_run_022.slicemocoxy_afni+orig','OutBrick_run_022.brain+orig','tempslmoco_volslc_alg_vol_OutBrick_run_022.slicemocoxy_afni/motion.wholevol_zt',CARD,RESP); exit;
< M A T L A B (R) >
Copyright 1984-2012 The MathWorks, Inc.
R2012b (8.0.0.783) 64-bit (glnxa64)
August 22, 2012
To get started, type one of these: helpwin, helpdesk, or demo.
For product information, visit http://www.mathworks.com.
>> Note: SMS acquisition is applied (MB acc. fac = 6).
Error in convert_1dmat_into_tranmatarray (line 3)
[zdim tdim veclen] = size(transmat1d_zt);
Output argument "tranmat_zt" (and maybe others) not assigned during call to
"/mrt/Applications/pestica4/slomoco/convert_1dmat_into_tranmatarray.m>convert_1dmat_into_tranmatarray".
Error in slicemoco_newalgorithm_input (line 25)
tranmat_zt = convert_1dmat_into_tranmatarray(transmat1d_zt);% tranmat_zt =
reshape(tranmat_zt,[zmbdim tdim 3 4]); % (v1) debugged
>> matlab -nojvm -nosplash addpath /mrt/Applications/pestica4/slomoco; addpath /mrt/Applications/pestica4/pestica_matlab; addpath /mrt/Applications/pestica4/afni_matlab; qa_slomoco('OutBrick_run_022.slicemocoxy_afni+orig','tempslmoco_volslc_alg_vol_OutBrick_run_022.slicemocoxy_afni/motion.wholevol_zt','tempslmocoxy_afni_OutBrick_run_022'); exit;
< M A T L A B (R) >
Copyright 1984-2012 The MathWorks, Inc.
R2012b (8.0.0.783) 64-bit (glnxa64)
August 22, 2012
To get started, type one of these: helpwin, helpdesk, or demo.
For product information, visit http://www.mathworks.com.
>> Note: SMS acquisition is applied (MB acc. fac = 6).
Index exceeds matrix dimensions.
Error in qa_slomoco (line 61)
outofplane(acqtp,:) = rotmat_zt(z,t,:); % [tdim*zmbdim [dL dP dS
rotmatvector (1:9)]]
>> *+ WARNING: If you are performing spatial transformations on an oblique dset,
such as ./OutBrick_run_022+orig.BRIK,
or viewing/combining it with volumes of differing obliquity,
you should consider running:
3dWarp -deoblique
on this and other oblique datasets in the same session.
See 3dWarp -help for details.
++ Oblique dataset:./OutBrick_run_022+orig.BRIK is 4.199786 degrees from plumb.
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as ./OutBrick_run_022+orig.BRIK,
or viewing/combining it with volumes of differing obliquity,
you should consider running:
3dWarp -deoblique
on this and other oblique datasets in the same session.
See 3dWarp -help for details.
++ Oblique dataset:./OutBrick_run_022+orig.BRIK is 4.199786 degrees from plumb.
*+ WARNING: If you are performing spatial transformations on an oblique dset,
such as ./OutBrick_run_022+orig.BRIK,
or viewing/combining it with volumes of differing obliquity,
you should consider running:
3dWarp -deoblique
on this and other oblique datasets in the same session.
See 3dWarp -help for details.
++ Oblique dataset:./OutBrick_run_022+orig.BRIK is 4.199786 degrees from plumb.
3drefit -saveatr -atrint TAXIS_NUMS 350 72 77002 -999 -999 -999 -999 -999 OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit: AFNI version=AFNI_16.2.17 (Sep 16 2016) [64-bit]
++ Authored by: RW Cox
** ERROR: Can't open dataset OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit processed 0 datasets
3drefit -saveatr -atrfloat TAXIS_FLOATS 0 1.05 0 -31.762 2 -999999 -999999 -999999 OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit: AFNI version=AFNI_16.2.17 (Sep 16 2016) [64-bit]
++ Authored by: RW Cox
** ERROR: Can't open dataset OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit processed 0 datasets
3drefit -saveatr -atrfloat TAXIS_OFFSETS 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 0 0.44 0.8775 0.265 0.7025 0.0875 0.5275 0.9675 0.3525 0.79 0.175 0.615 OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit: AFNI version=AFNI_16.2.17 (Sep 16 2016) [64-bit]
++ Authored by: RW Cox
** ERROR: Can't open dataset OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
++ 3drefit processed 0 datasets
3dNotes -h slicewise_moco_inplane.sh OutBrick_run_022+orig OutBrick_run_022.slicemocoxy_afni.slomoco_pestica+orig
Error: Cannot open dataset
===
Threaded View
Title | Author | Date |
---|---|---|
Christian Kaufmann | Jul 18, 2016 | |
Christian Kaufmann | Nov 4, 2016 | |
Christian Kaufmann | Sep 22, 2016 | |
Wanyong Shin | Oct 31, 2016 | |
Wanyong Shin | Jul 18, 2016 | |
Christian Kaufmann | Jul 20, 2016 | |