[Freesurfer] mri_watershed
We use an anatomical analysis program called Brain Image in our lab. We would like to speed up our analysis in Brain Image by replacing manual skull stripping with FreeSurfer's automated skull strip program. For this we would need to do a skull strip alone, without adjusting intensity values or contrast. Is this possible? Is this a reasonable thing to do? I have run mri_watershed with the input from orig on several subjects, and the output volume looks much brighter than the orig volume. Does mri_watershed normalize intensities or change the contrast in order to strip the skull? Thanks, Sasha Sasha Wolosin Research Assistant Developmental Cognitive Neurology Kennedy Krieger Institute 707 N. Broadway Baltimore, MD 21205 ph: (443) 923-9270 Disclaimer: The materials in this e-mail are private and may contain Protected Health Information. Please note that e-mail is not necessarily confidential or secure. Your use of e-mail constitutes your acknowledgment of these confidentiality and security limitations. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying, distribution, or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this e-mail in error, please immediately notify the sender via telephone or return e-mail. ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
[Freesurfer] process volume
hey guys, Im trying to process an orig scan from the terminal line: [EMAIL PROTECTED] orig]$ recon-all -subjid 1768-4 -segmentation INFO: SUBJECTS_DIR is /home/amc/freesurfer_new/freesurfer/subjects program versions used $Id: recon-all-nmr,v 1.72 2005/07/12 15:56:32 greve Exp $ $Id: mri_motion_correct2,v 1.9 2004/12/14 03:37:42 greve Exp $ mri_convert --version $ $Header: /software/source/minc/cvsroot/minc/progs/rawtominc/rawtominc.c,v 6.3 1999/10/19 14:45:31 neelin Exp $ minctracc.c $Header: /software/source/registration/mni_autoreg/minctracc/Main/minctracc.c,v 96.3 2000/03/15 08:42:41 stever Exp $ $Header: /software/source/minc/cvsroot/minc/progs/mincresample/mincresample.c,v 6.2 1999/10/19 14:45:27 neelin Exp $ $Header: /software/source/minc/cvsroot/minc/progs/mincaverage/mincaverage.c,v 6.1 1999/10/19 14:45:18 neelin Exp $ $Id: talairach2,v 1.7 2004/08/31 18:07:37 tosa Exp $ mri_convert --version $ # $Id: mritotal.in,v 1.10 2000/06/14 12:55:01 louis Exp $ # $Id: nu_estimate.in,v 1.15 1999/08/10 18:05:42 jgsled Exp $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ ### ~/freesurfer_new/freesurfer/subjects/1768-4 Using COR format - Segmentation Thu Aug 18 14:30:56 EDT 2005 ERROR: cannot find brain/COR-.info recon-all exited with errors at Thu Aug 18 14:30:56 EDT 2005 Notice how it says using COR format. The file I am trying to convert is in mgz format, using the new freesurfer version. I'm confused as to why it is doing this. Thanks, Dave David A. Soscia Neurosciences Advanced Imaging Research Center 47 New Scotland Ave. Albany, NY 12208-3479 Building S Rm. 132 P- (518)-262-0826 C- (518)-469-8317 [EMAIL PROTECTED] ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
[Freesurfer] segmentation fault during mri_convert
Hi, Is mri_convert setup yet to convert to mgz format? I have been trying it in both the std and dev environment from dicom format and spm format and i keep getting a segmentation fault. it still converts to cor format fine but i thought we were trying to get away from that. i dont know exactly whats causing it, but here is my command line and the running information. troy:mduff[84] mri_convert -ot mgz -it dicom 33749275 $SUBJECTS_DIR/SM2C1_autoreconall/mri/orig/002 mri_convert -ot mgz -it dicom 33749275 /space/troy/1/users/FREESURFER/SM2C1_autoreconall/mri/orig/002 reading from 33749275... Getting Series No Scanning Directory INFO: Found 131 files in /autofs/space/troy_002/users/fmridata/sm2c1_041118/3danat/dicom/028 INFO: Scanning for Series Number 28 INFO: found 128 files in series INFO: loading series header info. INFO: sorting. RunNo = 27 INFO: (256 256 128), nframes = 1, ismosaic=0 FileName /autofs/space/troy_002/users/fmridata/sm2c1_041118/3danat/dicom/028/33747925 Identification NumarisVersyngo MR 2004V 4VB11D ScannerModel Avanto PatientName sm2c1 Date and time StudyDate 20041118 StudyTime 180537.015000 SeriesTime190942.968000 AcqTime 190056.292500 Acquisition parameters PulseSeq tfl3d1_ns Protocol T1_MPRAGE_sag PhEncDir ROW EchoNo0 FlipAngle 7 EchoTime 3.31 InversionTime 1000 RepetitionTime2730 PhEncFOV 256 ReadoutFOV256 Image information RunNo 27 SeriesNo 28 ImageNo 1 NImageRows256 NImageCols256 NFrames 1 SliceArraylSize 1 IsMosaic 0 ImgPos 71.7733 130.0894 131.5385 VolRes 1. 1. 1.3300 VolDim256 256 128 Vc 0.1238 -0.9687 -0.2150 Vr -0.0454 0.2110 -0.9764 Vs -0.9913 -0.1306 0.0178 VolCenter -2.5670 21.9742 -19.4516 TransferSyntaxUID 1.2.840.10008.1.2.1 INFO: no Siemens slice order reversal detected (good!). 128 DICOM 3.0 files in list Found 128 DICOM Files WARNING: NumberOfFrames 0 != Found Count of slices 128. reading DICOM image... - DICOM meta-header file name /autofs/space/troy_002/users/fmridata/sm2c1_041118/3danat/dicom/028/33747925 Date and time study date 20041118 study time 180537.015000 series time 190942.968000 acquisition time190942.968000 Identification patient namesm2c1 manufacturerSIEMENS Dimensions number of rows 256 number of columns 256 number of frames128 pixel width 1 pixel height1 slice thickness 1.33 field of view 256 image number1 (might be not reliable) transfer syntax UID 1.2.840.10008.1.2.1 Acquisition parameters echo time 3.31 repetition time 2730 inversion time 1000 echo number 0 flip angle 0.122173 bits allocated 16 Spatial information first image position-71.7733 -130.089 131.539 last image position 95.662 -108.024 134.554 image orientation -0.123794 0.968728 -0.215038 0.0453834 -0.210951 -0.976442 - TR=2730.00, TE=3.31, TI=1000.00, flip angle=7.00 i_ras = (0.123794, -0.968728, -0.215038) j_ras = (-0.0453834, 0.210951, -0.976442) k_ras = (-0.99127, -0.130637, 0.0178497) writing to /space/troy/1/users/FREESURFER/SM2C1_autoreconall/mri/orig/002... Segmentation fault i saw a couple other problems on the list but nothing exactly like this. thanks for the help. margaret ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] mri_watershed
Dear Sasha, mri_watershed was initially designed to skull strip orig images, without any preprocessing step, such as intensity normalization or contrast modification. So, This certainly seems like a reasonable thing to do. However, mri_watershed might normalize intensities during the process but should NOT generate skull-stripped images with different intensities. We'll quickly look into this and let you know if this is a bug that needs to be fixed. Cheers, Florent On Fri, 19 Aug 2005, Sasha Wolosin wrote: We use an anatomical analysis program called Brain Image in our lab. We would like to speed up our analysis in Brain Image by replacing manual skull stripping with FreeSurfer's automated skull strip program. For this we would need to do a skull strip alone, without adjusting intensity values or contrast. Is this possible? Is this a reasonable thing to do? I have run mri_watershed with the input from orig on several subjects, and the output volume looks much brighter than the orig volume. Does mri_watershed normalize intensities or change the contrast in order to strip the skull? Thanks, Sasha Sasha Wolosin Research Assistant Developmental Cognitive Neurology Kennedy Krieger Institute 707 N. Broadway Baltimore, MD 21205 ph: (443) 923-9270 Disclaimer: The materials in this e-mail are private and may contain Protected Health Information. Please note that e-mail is not necessarily confidential or secure. Your use of e-mail constitutes your acknowledgment of these confidentiality and security limitations. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying, distribution, or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this e-mail in error, please immediately notify the sender via telephone or return e-mail. ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] process volume
It automatically tries to figure out what format you use for your analysis. It looks for mri/orig/XXX.mgz (which XXX is anything between 0-9). If it does not find that, it assumes COR. Specify -mgz to force mgz. David Soscia wrote: hey guys, Im trying to process an orig scan from the terminal line: [EMAIL PROTECTED] orig]$ recon-all -subjid 1768-4 -segmentation INFO: SUBJECTS_DIR is /home/amc/freesurfer_new/freesurfer/subjects program versions used $Id: recon-all-nmr,v 1.72 2005/07/12 15:56:32 greve Exp $ $Id: mri_motion_correct2,v 1.9 2004/12/14 03:37:42 greve Exp $ mri_convert --version $ $Header: /software/source/minc/cvsroot/minc/progs/rawtominc/rawtominc.c,v 6.3 1999/10/19 14:45:31 neelin Exp $ minctracc.c $Header: /software/source/registration/mni_autoreg/minctracc/Main/minctracc.c,v 96.3 2000/03/15 08:42:41 stever Exp $ $Header: /software/source/minc/cvsroot/minc/progs/mincresample/mincresample.c,v 6.2 1999/10/19 14:45:27 neelin Exp $ $Header: /software/source/minc/cvsroot/minc/progs/mincaverage/mincaverage.c,v 6.1 1999/10/19 14:45:18 neelin Exp $ $Id: talairach2,v 1.7 2004/08/31 18:07:37 tosa Exp $ mri_convert --version $ # $Id: mritotal.in,v 1.10 2000/06/14 12:55:01 louis Exp $ # $Id: nu_estimate.in,v 1.15 1999/08/10 18:05:42 jgsled Exp $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ ### ~/freesurfer_new/freesurfer/subjects/1768-4 Using COR format - Segmentation Thu Aug 18 14:30:56 EDT 2005 ERROR: cannot find brain/COR-.info recon-all exited with errors at Thu Aug 18 14:30:56 EDT 2005 Notice how it says using COR format. The file I am trying to convert is in mgz format, using the new freesurfer version. I'm confused as to why it is doing this. Thanks, Dave David A. Soscia Neurosciences Advanced Imaging Research Center 47 New Scotland Ave. Albany, NY 12208-3479 Building S Rm. 132 P- (518)-262-0826 C- (518)-469-8317 [EMAIL PROTECTED] ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer -- Douglas N. Greve, Ph.D. MGH-NMR Center [EMAIL PROTECTED] Phone Number: 617-724-2358 Fax: 617-726-7422 ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] segmentation fault during mri_convert
yes, this worked. thanks so much, margaret On Fri, 19 Aug 2005, Nick Schmansky wrote: Margaret, The specified output filename needs to have .mgz at the end. In your case, it would be: $SUBJECTS_DIR/SM2C1_autoreconall/mri/orig/002.mgz I noticed that $SUBJECTS_DIR/SM2C1_autoreconall/mri/orig/002 is a directory, which might be the source of the segfault (which, of course, shouldn't happen, and I will investigate this). So adding .mgz to your output file should work (I just tried it myself with your specified dicom file and it worked). Nick On Fri, 2005-08-19 at 12:42 -0400, Margaret Duff wrote: Hi, Is mri_convert setup yet to convert to mgz format? I have been trying it in both the std and dev environment from dicom format and spm format and i keep getting a segmentation fault. it still converts to cor format fine but i thought we were trying to get away from that. i dont know exactly whats causing it, but here is my command line and the running information. troy:mduff[84] mri_convert -ot mgz -it dicom 33749275 $SUBJECTS_DIR/SM2C1_autoreconall/mri/orig/002 mri_convert -ot mgz -it dicom 33749275 /space/troy/1/users/FREESURFER/SM2C1_autoreconall/mri/orig/002 reading from 33749275... Getting Series No Scanning Directory INFO: Found 131 files in /autofs/space/troy_002/users/fmridata/sm2c1_041118/3danat/dicom/028 INFO: Scanning for Series Number 28 INFO: found 128 files in series INFO: loading series header info. INFO: sorting. RunNo = 27 INFO: (256 256 128), nframes = 1, ismosaic=0 FileName /autofs/space/troy_002/users/fmridata/sm2c1_041118/3danat/dicom/028/33747925 Identification NumarisVersyngo MR 2004V 4VB11D ScannerModel Avanto PatientName sm2c1 Date and time StudyDate 20041118 StudyTime 180537.015000 SeriesTime190942.968000 AcqTime 190056.292500 Acquisition parameters PulseSeq tfl3d1_ns Protocol T1_MPRAGE_sag PhEncDir ROW EchoNo0 FlipAngle 7 EchoTime 3.31 InversionTime 1000 RepetitionTime2730 PhEncFOV 256 ReadoutFOV256 Image information RunNo 27 SeriesNo 28 ImageNo 1 NImageRows256 NImageCols256 NFrames 1 SliceArraylSize 1 IsMosaic 0 ImgPos 71.7733 130.0894 131.5385 VolRes 1. 1. 1.3300 VolDim256 256 128 Vc 0.1238 -0.9687 -0.2150 Vr -0.0454 0.2110 -0.9764 Vs -0.9913 -0.1306 0.0178 VolCenter -2.5670 21.9742 -19.4516 TransferSyntaxUID 1.2.840.10008.1.2.1 INFO: no Siemens slice order reversal detected (good!). 128 DICOM 3.0 files in list Found 128 DICOM Files WARNING: NumberOfFrames 0 != Found Count of slices 128. reading DICOM image... - DICOM meta-header file name /autofs/space/troy_002/users/fmridata/sm2c1_041118/3danat/dicom/028/33747925 Date and time study date 20041118 study time 180537.015000 series time 190942.968000 acquisition time190942.968000 Identification patient namesm2c1 manufacturerSIEMENS Dimensions number of rows 256 number of columns 256 number of frames128 pixel width 1 pixel height1 slice thickness 1.33 field of view 256 image number1 (might be not reliable) transfer syntax UID 1.2.840.10008.1.2.1 Acquisition parameters echo time 3.31 repetition time 2730 inversion time 1000 echo number 0 flip angle 0.122173 bits allocated 16 Spatial information first image position-71.7733 -130.089 131.539 last image position 95.662 -108.024 134.554 image orientation -0.123794 0.968728 -0.215038 0.0453834 -0.210951 -0.976442 - TR=2730.00, TE=3.31, TI=1000.00, flip angle=7.00 i_ras = (0.123794, -0.968728, -0.215038) j_ras = (-0.0453834, 0.210951, -0.976442) k_ras = (-0.99127, -0.130637, 0.0178497) writing to /space/troy/1/users/FREESURFER/SM2C1_autoreconall/mri/orig/002... Segmentation fault i saw a couple other problems on the list but nothing exactly like this. thanks for the help. margaret ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu
Re: [Freesurfer] mri_watershed
Hi Sasha, If you use tkmedit to visualize the results, then the contrast change you saw may be just due to the display problem (not exactly a problem) of tkmedit. You can move the mouse around the image to actually see whether the image intensity values were changed or not. I used mri_watershed a lot, and never saw it changed the original image intensity (in the brain region). But tkmedit will adjust its display contrast according to the dynamic range of the original image. So once skull is removed, the max and min intensity values of the image are very likely to be different, and the image will then look differently in tkmedit. -Xiao On Fri, 19 Aug 2005, florent segonne wrote: Dear Sasha, mri_watershed was initially designed to skull strip orig images, without any preprocessing step, such as intensity normalization or contrast modification. So, This certainly seems like a reasonable thing to do. However, mri_watershed might normalize intensities during the process but should NOT generate skull-stripped images with different intensities. We'll quickly look into this and let you know if this is a bug that needs to be fixed. Cheers, Florent On Fri, 19 Aug 2005, Sasha Wolosin wrote: We use an anatomical analysis program called Brain Image in our lab. We would like to speed up our analysis in Brain Image by replacing manual skull stripping with FreeSurfer's automated skull strip program. For this we would need to do a skull strip alone, without adjusting intensity values or contrast. Is this possible? Is this a reasonable thing to do? I have run mri_watershed with the input from orig on several subjects, and the output volume looks much brighter than the orig volume. Does mri_watershed normalize intensities or change the contrast in order to strip the skull? Thanks, Sasha Sasha Wolosin Research Assistant Developmental Cognitive Neurology Kennedy Krieger Institute 707 N. Broadway Baltimore, MD 21205 ph: (443) 923-9270 Disclaimer: The materials in this e-mail are private and may contain Protected Health Information. Please note that e-mail is not necessarily confidential or secure. Your use of e-mail constitutes your acknowledgment of these confidentiality and security limitations. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying, distribution, or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this e-mail in error, please immediately notify the sender via telephone or return e-mail. ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] coregistration
You may try FLIRT from FSL. http://www.fmrib.ox.ac.uk/fsl/flirt/ -Xiao On Fri, 19 Aug 2005 [EMAIL PROTECTED] wrote: I always have a problem to coregistrate T1 images to Talairach, Is there any way to improve the accuracy of coregistration besides tkregister2? Thanks, Antao ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
[Freesurfer] Permission denied
Hello, Recently moved freesurfer from one system to another, but I am having trouble executing it. ./csurf: Permission denied. I have license in the right folder. Has anyone experienced the same problem. Thanks, Arvind ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] Permission denied
Recently moved freesurfer from one system to another, but I am having trouble executing it. ./csurf: Permission denied. I have license in the right folder. Has anyone experienced the same problem. Make sure it is executable. Try: chmod a+x csurf -- Kevin Teich ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
[Freesurfer] process volume again
hey guys, Now I am getting this error when trying to process the volume: [EMAIL PROTECTED] 1768-4]$ recon-all -subjid 1768-4 -mgz -segmentation INFO: SUBJECTS_DIR is /home/amc/freesurfer_new/freesurfer/subjects program versions used $Id: recon-all-nmr,v 1.72 2005/07/12 15:56:32 greve Exp $ $Id: mri_motion_correct2,v 1.9 2004/12/14 03:37:42 greve Exp $ mri_convert --version $ $Header: /software/source/minc/cvsroot/minc/progs/rawtominc/rawtominc.c,v 6.3 1999/10/19 14:45:31 neelin Exp $ minctracc.c $Header: /software/source/registration/mni_autoreg/minctracc/Main/minctracc.c,v 96.3 2000/03/15 08:42:41 stever Exp $ $Header: /software/source/minc/cvsroot/minc/progs/mincresample/mincresample.c,v 6.2 1999/10/19 14:45:27 neelin Exp $ $Header: /software/source/minc/cvsroot/minc/progs/mincaverage/mincaverage.c,v 6.1 1999/10/19 14:45:18 neelin Exp $ $Id: talairach2,v 1.7 2004/08/31 18:07:37 tosa Exp $ mri_convert --version $ # $Id: mritotal.in,v 1.10 2000/06/14 12:55:01 louis Exp $ # $Id: nu_estimate.in,v 1.15 1999/08/10 18:05:42 jgsled Exp $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ ### Using mgz format - Segmentation Fri Aug 19 13:30:42 EDT 2005 ERROR: cannot find brain.mgz recon-all exited with errors at Fri Aug 19 13:30:42 EDT 2005 i dont have a brain.mgz in the brain directory because I thought the purpose of this process was to create it. Any suggestions? Dave David A. Soscia Neurosciences Advanced Imaging Research Center 47 New Scotland Ave. Albany, NY 12208-3479 Building S Rm. 132 P- (518)-262-0826 C- (518)-469-8317 [EMAIL PROTECTED] ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
[Freesurfer] new freesurfer atlases for dev
For those in the NMR center using the freesurfer dev environment, there will re a change in recon-all to use "new" atlaes by default. This change will be live on Monday. The new cortical atlases are based on automatically filling ventricles and subcortical structures using the automatic subcortical segmentations (asegs). This is part of the "autorecon" stream accessed with the -autorecon flags (ie, -autorecon-all, -autorecon1, etc). The new subcortical atlas is based on 40 subjects. To continue to use the old atlases, run recon-all with -old-atlas. If you use any of the "stage" flags (ie, -all, -stage1, -stage2, etc), recon-all will automatically use the old atlases to maintain consistency. Specifications of individual steps (eg, -surfreg) will use the new atlases. To help keep track of which atlases were used, recon-all now stores he names of the atlases used for a subject in the following files located in the subject's scripts directory: mri_em_register.atlas.dat - subcortical atlas (new=RB40_talairach_2005-07-20.gca, old=talairach_mixed.gca) mris_register.atlas.dat - cortical registration atlas (new=average.curvature.filled.buckner40.tif, old=average.tif) mris_ca_label.atlas.dat - cortical label atlas (new=curvature.buckner40.filled.desikan_killiany.gcs, old=atlas2002_simple.gcs) If these files do not exist, then the old atlases are implied. doug -- Douglas N. Greve, Ph.D. MGH-NMR Center [EMAIL PROTECTED] Phone Number: 617-724-2358 Fax: 617-726-7422 ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] mri_watershed
no, it should just strip. Is that not happening? If not, you can always use mri_mask to use the skull stripped volume as a mask to apply to whatever volume (in register with the orig) you want. cheers, Bruce On Fri, 19 Aug 2005, Sasha Wolosin wrote: We use an anatomical analysis program called Brain Image in our lab. We would like to speed up our analysis in Brain Image by replacing manual skull stripping with FreeSurfer's automated skull strip program. For this we would need to do a skull strip alone, without adjusting intensity values or contrast. Is this possible? Is this a reasonable thing to do? I have run mri_watershed with the input from orig on several subjects, and the output volume looks much brighter than the orig volume. Does mri_watershed normalize intensities or change the contrast in order to strip the skull? Thanks, Sasha Sasha Wolosin Research Assistant Developmental Cognitive Neurology Kennedy Krieger Institute 707 N. Broadway Baltimore, MD 21205 ph: (443) 923-9270 Disclaimer: The materials in this e-mail are private and may contain Protected Health Information. Please note that e-mail is not necessarily confidential or secure. Your use of e-mail constitutes your acknowledgment of these confidentiality and security limitations. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying, distribution, or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this e-mail in error, please immediately notify the sender via telephone or return e-mail. ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] segmentation fault during mri_convert
you need to specify a file with the extension .mgz as the output name On Fri, 19 Aug 2005, Margaret Duff wrote: Hi, Is mri_convert setup yet to convert to mgz format? I have been trying it in both the std and dev environment from dicom format and spm format and i keep getting a segmentation fault. it still converts to cor format fine but i thought we were trying to get away from that. i dont know exactly whats causing it, but here is my command line and the running information. troy:mduff[84] mri_convert -ot mgz -it dicom 33749275 $SUBJECTS_DIR/SM2C1_autoreconall/mri/orig/002 mri_convert -ot mgz -it dicom 33749275 /space/troy/1/users/FREESURFER/SM2C1_autoreconall/mri/orig/002 reading from 33749275... Getting Series No Scanning Directory INFO: Found 131 files in /autofs/space/troy_002/users/fmridata/sm2c1_041118/3danat/dicom/028 INFO: Scanning for Series Number 28 INFO: found 128 files in series INFO: loading series header info. INFO: sorting. RunNo = 27 INFO: (256 256 128), nframes = 1, ismosaic=0 FileName /autofs/space/troy_002/users/fmridata/sm2c1_041118/3danat/dicom/028/33747925 Identification NumarisVersyngo MR 2004V 4VB11D ScannerModel Avanto PatientName sm2c1 Date and time StudyDate 20041118 StudyTime 180537.015000 SeriesTime190942.968000 AcqTime 190056.292500 Acquisition parameters PulseSeq tfl3d1_ns Protocol T1_MPRAGE_sag PhEncDir ROW EchoNo0 FlipAngle 7 EchoTime 3.31 InversionTime 1000 RepetitionTime2730 PhEncFOV 256 ReadoutFOV256 Image information RunNo 27 SeriesNo 28 ImageNo 1 NImageRows256 NImageCols256 NFrames 1 SliceArraylSize 1 IsMosaic 0 ImgPos 71.7733 130.0894 131.5385 VolRes 1. 1. 1.3300 VolDim256 256 128 Vc 0.1238 -0.9687 -0.2150 Vr -0.0454 0.2110 -0.9764 Vs -0.9913 -0.1306 0.0178 VolCenter -2.5670 21.9742 -19.4516 TransferSyntaxUID 1.2.840.10008.1.2.1 INFO: no Siemens slice order reversal detected (good!). 128 DICOM 3.0 files in list Found 128 DICOM Files WARNING: NumberOfFrames 0 != Found Count of slices 128. reading DICOM image... - DICOM meta-header file name /autofs/space/troy_002/users/fmridata/sm2c1_041118/3danat/dicom/028/33747925 Date and time study date 20041118 study time 180537.015000 series time 190942.968000 acquisition time190942.968000 Identification patient namesm2c1 manufacturerSIEMENS Dimensions number of rows 256 number of columns 256 number of frames128 pixel width 1 pixel height1 slice thickness 1.33 field of view 256 image number1 (might be not reliable) transfer syntax UID 1.2.840.10008.1.2.1 Acquisition parameters echo time 3.31 repetition time 2730 inversion time 1000 echo number 0 flip angle 0.122173 bits allocated 16 Spatial information first image position-71.7733 -130.089 131.539 last image position 95.662 -108.024 134.554 image orientation -0.123794 0.968728 -0.215038 0.0453834 -0.210951 -0.976442 - TR=2730.00, TE=3.31, TI=1000.00, flip angle=7.00 i_ras = (0.123794, -0.968728, -0.215038) j_ras = (-0.0453834, 0.210951, -0.976442) k_ras = (-0.99127, -0.130637, 0.0178497) writing to /space/troy/1/users/FREESURFER/SM2C1_autoreconall/mri/orig/002... Segmentation fault i saw a couple other problems on the list but nothing exactly like this. thanks for the help. margaret ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] mri_watershed
Alternatively, depending on what version of tkmedit you have, you can use the -mm-main min max option to specify an out-of-range min and max visible value on the command line. As Xiao said, normally the min and max visible values that determine the color scale are preset to the actual min and max values in the volume. In the View-Configure-Brightness / Constrast dialog box, you cannot specify values outside of this range. However, if on the command, you do: tkmedit bert T1 -mm-main 0 500 This will set the min and max visible values to 0 and 500, even though 500 is above the max value in the volume. Note that your version of tkmedit may not have this option, so if you get an error when trying this, ignore this email. :) On Fri, 19 Aug 2005, Xiao Han wrote: Hi Sasha, If you use tkmedit to visualize the results, then the contrast change you saw may be just due to the display problem (not exactly a problem) of tkmedit. You can move the mouse around the image to actually see whether the image intensity values were changed or not. I used mri_watershed a lot, and never saw it changed the original image intensity (in the brain region). But tkmedit will adjust its display contrast according to the dynamic range of the original image. So once skull is removed, the max and min intensity values of the image are very likely to be different, and the image will then look differently in tkmedit. -Xiao On Fri, 19 Aug 2005, florent segonne wrote: Dear Sasha, mri_watershed was initially designed to skull strip orig images, without any preprocessing step, such as intensity normalization or contrast modification. So, This certainly seems like a reasonable thing to do. However, mri_watershed might normalize intensities during the process but should NOT generate skull-stripped images with different intensities. We'll quickly look into this and let you know if this is a bug that needs to be fixed. Cheers, Florent On Fri, 19 Aug 2005, Sasha Wolosin wrote: We use an anatomical analysis program called Brain Image in our lab. We would like to speed up our analysis in Brain Image by replacing manual skull stripping with FreeSurfer's automated skull strip program. For this we would need to do a skull strip alone, without adjusting intensity values or contrast. Is this possible? Is this a reasonable thing to do? I have run mri_watershed with the input from orig on several subjects, and the output volume looks much brighter than the orig volume. Does mri_watershed normalize intensities or change the contrast in order to strip the skull? Thanks, Sasha Sasha Wolosin Research Assistant Developmental Cognitive Neurology Kennedy Krieger Institute 707 N. Broadway Baltimore, MD 21205 ph: (443) 923-9270 Disclaimer: The materials in this e-mail are private and may contain Protected Health Information. Please note that e-mail is not necessarily confidential or secure. Your use of e-mail constitutes your acknowledgment of these confidentiality and security limitations. If you are not the intended recipient, be advised that any unauthorized use, disclosure, copying, distribution, or the taking of any action in reliance on the contents of this information is strictly prohibited. If you have received this e-mail in error, please immediately notify the sender via telephone or return e-mail. ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer -- Kevin Teich ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Re: [Freesurfer] process volume again
Not sure what's going on there. You may have gotten into a situation where some of the processing was done in COR and some im MGZ. It will assume to do things in mgz unless you specify -COR. One thing you can do to get things synched back up is to convert all your CORs to mgzs, eg, cd mri mri_convert brain brain.mgz doug David Soscia wrote: hey guys, Now I am getting this error when trying to process the volume: [EMAIL PROTECTED] 1768-4]$ recon-all -subjid 1768-4 -mgz -segmentation INFO: SUBJECTS_DIR is /home/amc/freesurfer_new/freesurfer/subjects program versions used $Id: recon-all-nmr,v 1.72 2005/07/12 15:56:32 greve Exp $ $Id: mri_motion_correct2,v 1.9 2004/12/14 03:37:42 greve Exp $ mri_convert --version $ $Header: /software/source/minc/cvsroot/minc/progs/rawtominc/rawtominc.c,v 6.3 1999/10/19 14:45:31 neelin Exp $ minctracc.c $Header: /software/source/registration/mni_autoreg/minctracc/Main/minctracc.c,v 96.3 2000/03/15 08:42:41 stever Exp $ $Header: /software/source/minc/cvsroot/minc/progs/mincresample/mincresample.c,v 6.2 1999/10/19 14:45:27 neelin Exp $ $Header: /software/source/minc/cvsroot/minc/progs/mincaverage/mincaverage.c,v 6.1 1999/10/19 14:45:18 neelin Exp $ $Id: talairach2,v 1.7 2004/08/31 18:07:37 tosa Exp $ mri_convert --version $ # $Id: mritotal.in,v 1.10 2000/06/14 12:55:01 louis Exp $ # $Id: nu_estimate.in,v 1.15 1999/08/10 18:05:42 jgsled Exp $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ ### Using mgz format - Segmentation Fri Aug 19 13:30:42 EDT 2005 ERROR: cannot find brain.mgz recon-all exited with errors at Fri Aug 19 13:30:42 EDT 2005 i dont have a brain.mgz in the brain directory because I thought the purpose of this process was to create it. Any suggestions? Dave David A. Soscia Neurosciences Advanced Imaging Research Center 47 New Scotland Ave. Albany, NY 12208-3479 Building S Rm. 132 P- (518)-262-0826 C- (518)-469-8317 [EMAIL PROTECTED] ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer -- Douglas N. Greve, Ph.D. MGH-NMR Center [EMAIL PROTECTED] Phone Number: 617-724-2358 Fax: 617-726-7422 ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer