Hi everyone,
I'm wondering whether there are multiple versions of recon-all v 1.313.2.6 (2010/08/04 15:50:56), and/or mri_ca_normalize stable5.
I have been testing recon-all and compared our output against the sample processed-data in FS/TUtorial. I found that the output obtained was different from the sample, despite the fact that I ran recon-all with the same options as the sample outputs, on the same input files.
recon-all -s bert -i /analysis/Tests/FSTest_Tossi/Samples2/buckner_data/tutorial_subjs/bert/mri/orig/001.mgz -i /analysis/Tests/FSTest_Tossi/Samples2/buckner_data/tutorial_subjs/bert/mri/orig/002.mgz -i /analysis/Tests/FSTest_Tossi/Samples2/buckner_data/tutorial_subjs/bert/mri/orig/003.mgz -all -debug -clean -norandomness -allowcoredump -time -subjid BertByTossi20110119A
I performed a line-by-line comparison of the recon-all log between the sample and one from my run and found that our log differs from the sample (buckner_data/tutorial_subjs/bert/scripts/recon-all.log) at line 1034, where mri_ca_normalize runs the first path of smoothing (with sigma=0.200). Whereas the sample has 38 recursions, my run had 18 recursions for this section. The outputs hereafter differed from each other.
I used the same version as the sample was processed, recon-all: v 1.313.2.6 (2010/08/04 15:50:56), and mri_ca_normalize: stable5. Not to mention, the input files are identical.
Also, I ran the exact same line of the command more than three times, and obtained the exact same result. Therefore, I believe this is not due to any randomization introduced (if any) in normalization/smoothing. I used an x86_64 processor, the same as the sample. Thus, I do not believe this is due to any floating point resolution difference between architectures, either.
The only possible reason I can think of that might account for the discrepancy is that there are multiple versions or compiles "stable5" mri_ca_normalize.
I obtained the sample input and result from; ftp://surfer.nmr.mgh.harvard.edu/pub/data/bert.recon.tgz Our download md5sum matches the one posted.
Here's a list of input files; $ls -l /analysis/Tests/FSTest_Tossi/Samples2/buckner_data/tutorial_subjs/bert/mri/orig/ total 17776 -rw-rw-r-- 1 tossi mrilab 5996899 2010-08-09 13:01 001.mgz -rw-rw-r-- 1 tossi mrilab 6077242 2010-08-09 13:01 002.mgz -rw-rw-r-- 1 tossi mrilab 6083164 2010-08-09 13:01 003.mgz
Thank you very much in advance.
Tossi
/********************************************* Tossi Ikuta, Ph.D. Research Scientist The Feinstein Institute for Medical Research The Zucker Hillside Hospital (718) 470-4486 tikuta@nshs.edu *********************************************/