Thanks to Eugenio, for the tip!
We did look at (orig,nu,T1,brainmask,norm,nu,noneck, aseg, wm,
filled).mgz. The only one showing a difference where our white matter
error occurs is filled.mgz.
In some of the other previous images are minor differences around that
area, but not overlapping the area. This starts with the brainmask.mgz.
The wm.mgz looks almost identical though and this is what is used in
making filled.mgz, right?
Cheers,
Mel and Vincent
P.S.: Here's a quick script to visualize all differences:
recon1='/data1/vbeliveau/atlas/proc/MR/recon/v0013_good'; %change this
recon2='/data1/vbeliveau/atlas/proc/MR/recon/v0013_bad'; %change this
dest='/data1/vbeliveau/atlas/test'; %change this
img1=MRIread([recon1 '/mri/orig/001.mgz']);
img2=MRIread([recon2 '/mri/orig/001.mgz']);
img1.vol=img1.vol-img2.vol;
MRIwrite(img1,[dest '/diff_raw.mgz']);
disp('raw');
sum(img1.vol(:)>10e-4)
vol={'orig','nu','T1','brainmask','norm','nu_noneck','aseg','brain','wm','filled'};
for n=1:numel(vol)
img1=MRIread([recon1 '/mri/' vol{n} '.mgz']);
img2=MRIread([recon2 '/mri/' vol{n} '.mgz']);
img1.vol=img1.vol-img2.vol;
MRIwrite(img1,[dest '/diff_' vol{n} '.mgz']);
disp(vol{n});
sum(img1.vol(:)>10e-4)
end
-------- Original Message --------
Subject: Re: [Freesurfer] success and failure of wm surfaces in
(almost) identical structural image
Date: Tue, 26 Aug 2014 14:05:19 +0200 (CEST)
From: Eugenio Iglesias <e.igles...@bcbl.eu>
To: Melanie Ganz <melanie.g...@nru.dk>
It'd be really helpful if you looked at where the two pipelines start to
diverge. It's apparent that the origs are the same, and that the surfaces are
different, but where does the divergence begin? Are the nu.mgz's similar? How
about the the brainmask.mgz's?
Here's the recon table:
https://surfer.nmr.mgh.harvard.edu/fswiki/ReconAllDevTable
Best,
/E
Juan Eugenio Iglesias
Postdoctoral researcher BCBL
www.jeiglesias.com
www.bcbl.eu
Legal disclaimer/Aviso legal/Lege-oharra: www.bcbl.eu/legal-disclaimer
----- Original Message -----
From: "Melanie Ganz" <melanie.g...@nru.dk>
To: freesurfer@nmr.mgh.harvard.edu
Sent: Tuesday, August 26, 2014 7:50:23 AM
Subject: [Freesurfer] success and failure of wm surfaces in (almost) identical
structural image
Dear Freesurfer community,
We've encountered a strange situation where the pial and wm surface
delineation is successful for one image but contains wm (and subsequent
pial) errors for another, almost identical structural image (see
attached image, red and yellow are for the successful surface and blue
and green are for the failed one). We've tried to identify the source
for this discrepancy but I'm at a loss. All voxels of the original input
images are identical to the 4th decimal and when converted to orig.mgz
only a few voxels appear to be different due to rounding, none of which
overlap with the wm surface errors. Both images were processed on the
same setup (hardware and software, Freesurfer stable 5.3).
Unfortunately, this error is quite systematic in my dataset; about 10%
of the images show this type of error. Any suggestions on how to
investigate this? I've uploaded example of good and bad at
recon ftp://surfer.nmr.mgh.harvard.edu/transfer/incoming/failed_wm_recon.tar.gz
in case you'd like to have a closer look.
Thanks for your help.
Vincent and Melanie
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
http://www.partners.org/complianceline . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
http://www.partners.org/complianceline . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.
