|
Size: 2902
Comment:
|
Size: 4445
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 12: | Line 12: |
| All the runs/series for your subject's visit will be in this directory. | All the runs/series for your subject's visit will be in this directory. You can do more elaborate searches with findsession. To get more info, run findsession with -help. |
| Line 16: | Line 16: |
| Unpacking is the processes through which you convert the data from DICOM format (which is not very useful) to a format that can be used by an analysis programe (eg, bhdr, nifti, analyze, mgh, mgz). The program that does this is unpacksdcmdir (this is only used to unpack Siemens DICOM files). mri_convert can also be used to convert individual runs, whereas unpacksdcmdir can be used to unpack any or all data from a directory. | Unpacking is the processes through which you convert the data from DICOM format (which is not very useful) to a format that can be used by an analysis programe (eg, bhdr, nifti, analyze, mgh, mgz). The program that does this is unpacksdcmdir (this is only used to unpack Siemens DICOM files). mri_convert can also be used to convert individual runs, whereas unpacksdcmdir can be used to unpack any or all data from a directory. Below we go through some simple unpacking. unpacksdcmdir can do a lot more. Run it with -help to get full documentation. |
| Line 18: | Line 18: |
| The first step in unpacking is to decide where you want your data to be stored. Usually, you choose a "parent" directory underwhich all of the visits for a study will be stored. Eg, /space/data/1/users/you/data-parent. | The first step in unpacking is to decide where you want your data to be stored. Usually, you choose a "parent" directory underwhich all of the visits for a study will be stored. Eg, /space/data/1/users/you/data-parent. The the actual target data for a given visit from a subject would be stored under something like /space/data/1/users/you/data-parent/yoursubject. |
| Line 22: | Line 22: |
| ["unpacksdcmdir"] -src dicomdir -targ /space/data/1/users/you/data-parent/yoursubject -scanonly /space/data/1/users/you/data-parent/yoursubject/scan.info | ["unpacksdcmdir"] -src dicomdir -targ targetdir -scanonly targetdir/scan.info |
| Line 24: | Line 24: |
| where dicomdir is the directory where the dicom data reside (as found above by findsession), targetdir is where you want the individual's data to go, and the scan.info file will contain a list of run/series numbers along with the scanning protocol that was used to acquire the data, and will look something like this: | |
| Line 25: | Line 26: |
| 1 circle_localizer ok 256 256 3 1 74407231[[BR]] 2 ge_functionals ok 64 64 10 146 74407258[[BR]] 3 ge_functionals ok 64 64 10 146 74406408[[BR]] 4 ge_functionals ok 64 64 10 146 74405551[[BR]] 5 ge_functionals ok 64 64 10 146 74401663[[BR]] 6 ge_functionals ok 64 64 25 3 74400790[[BR]] 7 ep2d_T1w ok 64 64 10 1 74400835[[BR]] 8 tfl3d1_ns T1_MPRAGE_sag ok 256 256 128 1 74401012[[BR]] 9 tfl3d1_ns T1_MPRAGE_sag ok 256 256 128 1 74402212[[BR]] |
|
| Line 26: | Line 36: |
| To actually do the unpacking, you will need to run unpacksdcmdir again with a slightly different set of arguments: | |
| Line 27: | Line 38: |
| NOTE: you must be in the targ dir for that command to work. This outputs the type of scans used. Now use a file at $FREESURFER_HOME/numaris4-protocols.unpackcfg |
["unpacksdcmdir"] -src dicomdir -targ targetdir -fsfast \ -run 2 bold bshort f \[[BR]] -run 3 bold bshort f \[[BR]] -run 4 bold bshort f \[[BR]] -run 5 bold bshort f \[[BR]] -run 6 bold bshort f \[[BR]] -run 8 3danat mgz 001.mgz \[[BR]] -run 9 3danat mgz 001.mgz |
| Line 31: | Line 47: |
| to correctly make a config file. If you cannot find the file, check here [http://www.nmr.mgh.harvard.edu/~daniel/links/how/numaris4-protocols.unpackcfg] |
The -fsfast argument tells unpacksdcmdir to use the FS-FAST directory structure (if you want to specify a structure explicitly, use -generic). Each -run flag is followed by 4 arguments: (1) the run number (as found in the scan.info file), (2) a directory, (3) format, and (4) file name. In fsfast directory structure, run 2 will be stored in targetdir/bold/002 in bshort format. bshort format will have many files that start with "f" (ie, the file name given as the 4th argument). Each of the fMRI runs (2-6) will be stored under the "bold" directory (this is also known as the "functional subdirectory" or fsd). Note that run 7 was not unpacked. Note also that the backslashes ("\") are unix command-line continuation characters and not part of the unpacksdcmdir arguments list. Runs 8 and 9 are an anatomicals which can be used in FreeSurfer processing; thay will be stored in targetdir/3danat/008/001.mgz and 009/001.mgz (this is compressed mgh format). You can launch FreeSurfer on these volumes with: |
| Line 34: | Line 49: |
| The format is: num folder format fileprefix Here is an example config file. ||1 scout bshort f|| ||3 bold bshort f|| ||4 bold bshort f|| ||5 bold bshort f|| ||6 bold bshort f|| ||7 bold bshort f|| ||8 bold bshort f|| The number is acquisition number. The folder is where that acquisition is stored (bold directory, scout directory etc) The format is the form of the data files. The fileprefix is how the filenames start. |
recon-all -s yoursubjectname -i targetdir/3danat/008/001.mgz -i targetdir/3danat/009/001.mgz -all |
| Line 49: | Line 51: |
| == Next run it again with your config file. == ["unpacksdcmdir"] -src <source> -targ <target> -cfg<configfile> -fsfast == If it is .ima format data... == Note: I have come across old data which was *.ima format. This can be unpacked using ["unpackimadir2"] which works just like unpacksdcmdir. == Check for TR bug == *** There is always a problem after unpacking that the TR is not set correctly. Please fix it. pico $SUBJECTS_DIR/$SUBJECT/bold/seq.info == Copy the 3danat == Goto the 3danat/ dir cd $SUBJECTS_DIR/$SUBJECT/ Copy contents to mri/orig. For example... cp 3danat/* mri/orig |
where "yoursubjectname" is the FreeSurfer name for the subject as found in $SUBJECTS_DIR. Alternatively, you can create yoursubjname/mri/orig can copy the mgz files into this directory as 001.mgz and 002.mgz and then launch recon-all. |
Unpack the data from bourget:
Find your subject...
To find your data, run the findsession command, passing it the "patient name", which is the name that you registered the subject under when scanning:
["findsession patientname"]
This will give you a list of all matching data sets, including paths to where each one can be found. Eg,
/space/chartres/6/siemens/TRIO-20501-20020416-124756-468000
All the runs/series for your subject's visit will be in this directory. You can do more elaborate searches with findsession. To get more info, run findsession with -help.
Unpack the data
Unpacking is the processes through which you convert the data from DICOM format (which is not very useful) to a format that can be used by an analysis programe (eg, bhdr, nifti, analyze, mgh, mgz). The program that does this is unpacksdcmdir (this is only used to unpack Siemens DICOM files). mri_convert can also be used to convert individual runs, whereas unpacksdcmdir can be used to unpack any or all data from a directory. Below we go through some simple unpacking. unpacksdcmdir can do a lot more. Run it with -help to get full documentation.
The first step in unpacking is to decide where you want your data to be stored. Usually, you choose a "parent" directory underwhich all of the visits for a study will be stored. Eg, /space/data/1/users/you/data-parent. The the actual target data for a given visit from a subject would be stored under something like /space/data/1/users/you/data-parent/yoursubject.
Next, you need to find out exactly what data are in the dicom directory. Specificially you need to know which run/series number corresponds to which acquisition. To do this run:
- ["unpacksdcmdir"] -src dicomdir -targ targetdir -scanonly targetdir/scan.info
where dicomdir is the directory where the dicom data reside (as found above by findsession), targetdir is where you want the individual's data to go, and the scan.info file will contain a list of run/series numbers along with the scanning protocol that was used to acquire the data, and will look something like this:
1 circle_localizer ok 256 256 3 1 74407231BR 2 ge_functionals ok 64 64 10 146 74407258BR 3 ge_functionals ok 64 64 10 146 74406408BR 4 ge_functionals ok 64 64 10 146 74405551BR 5 ge_functionals ok 64 64 10 146 74401663BR 6 ge_functionals ok 64 64 25 3 74400790BR 7 ep2d_T1w ok 64 64 10 1 74400835BR 8 tfl3d1_ns T1_MPRAGE_sag ok 256 256 128 1 74401012BR 9 tfl3d1_ns T1_MPRAGE_sag ok 256 256 128 1 74402212BR
To actually do the unpacking, you will need to run unpacksdcmdir again with a slightly different set of arguments:
["unpacksdcmdir"] -src dicomdir -targ targetdir -fsfast \
-run 2 bold bshort f \BR -run 3 bold bshort f \BR -run 4 bold bshort f \BR -run 5 bold bshort f \BR -run 6 bold bshort f \BR -run 8 3danat mgz 001.mgz \BR -run 9 3danat mgz 001.mgz
The -fsfast argument tells unpacksdcmdir to use the FS-FAST directory structure (if you want to specify a structure explicitly, use -generic). Each -run flag is followed by 4 arguments: (1) the run number (as found in the scan.info file), (2) a directory, (3) format, and (4) file name. In fsfast directory structure, run 2 will be stored in targetdir/bold/002 in bshort format. bshort format will have many files that start with "f" (ie, the file name given as the 4th argument). Each of the fMRI runs (2-6) will be stored under the "bold" directory (this is also known as the "functional subdirectory" or fsd). Note that run 7 was not unpacked. Note also that the backslashes ("\") are unix command-line continuation characters and not part of the unpacksdcmdir arguments list. Runs 8 and 9 are an anatomicals which can be used in FreeSurfer processing; thay will be stored in targetdir/3danat/008/001.mgz and 009/001.mgz (this is compressed mgh format). You can launch FreeSurfer on these volumes with:
recon-all -s yoursubjectname -i targetdir/3danat/008/001.mgz -i targetdir/3danat/009/001.mgz -all
where "yoursubjectname" is the FreeSurfer name for the subject as found in $SUBJECTS_DIR. Alternatively, you can create yoursubjname/mri/orig can copy the mgz files into this directory as 001.mgz and 002.mgz and then launch recon-all.
