|
Size: 2273
Comment:
|
Size: 4136
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 3: | Line 3: |
| The description of the current longitudinal stream can be found here: [[LongitudinalProcessing]] | A detailed description of the current longitudinal stream can be found here: [[LongitudinalProcessing]] == FS Version 5.1 Preview == * Will include the fixed binaries * Will resample all time points to the base space (in the motion correction / conformalize step to avoid additional resampling) * Will include revised/improved editing We recommend to wait for this release. Initial studies show improved results due to the common voxel space. Furthermore, manual inspection is significantly simplified. == FS Version 5.0 == '''!!! Do not use FS 5.0 for longitudinal processing, before updating some binaries !!!''' * A non-related, last minute, incorrect 'bug-fix' in the transformation library introduced a new bug that breaks mri_robust_template, mri_robust_register (half way space). * You can find fixed versions of these binaries for your system [[ http://surfer.nmr.mgh.harvard.edu/pub/dist/freesurfer/misc | here ]]. * Also grab a modified recon-all from [[http://surfer.nmr.mgh.harvard.edu/pub/dist/freesurfer/misc/allplatforms | here]] that fixes the partial matching in the sanity check. Improvements/Fixes: * More robust base space estimation: keeps the base image centered * Added sanity checks (if longitudinal TP is in base and if baseid is different from the time points id's) No modifications to the data flow in the base or long stream. MartinReuter |
| Line 7: | Line 30: |
| Improvements to the new longitudinal stream (only on volume related issues, surface improvements to come). | Improvements to the new longitudinal stream (modified base stream, fixes/improvements in the long stream). |
| Line 10: | Line 33: |
| * only '''one''' template estimation on the norm.mgz files of all timepoints * a new block '''-base-init''' added to recon-all for the initial template creation * this creates the maps from each TP to base and the '''norm_template''' * then create an orig_template as orig/001.mgz to initialize the base run * the norm_template is used as brainmask * in gcareg and canorm the norm_template is used instead of the nu * this setup reduced the difficulties of dealing with two spaces (two template estimations) * should lead to significant run time improvements especially with several time points |
* Only '''one''' template estimation on the norm.mgz files of all timepoints, to prevent difficulties with two possibly different template locations (T1 and norm as in 4.4). * A new block '''-base-init''' added to recon-all for the initial template creation, to be able to call later parts of the script individually without re-running the initialization everytime. * The -base-init block creates the maps from each TP to base and the '''norm_template''' (using the norm.mgz of all TPs). * Then the orig.mgz of all TPs are resampled and averaged as orig/001.mgz to initialize the base run. * The norm_template is used as brainmask (as it contains only brain). * In -gcareg and -canorm the norm_template is used instead of the nu.mgz. * This setup (with only one estimation) should lead to significant run time improvements especially with several time points. |
| Line 20: | Line 42: |
| * mri_ca_label now uses the correct intensity scaling factors of the base * the talairach.lta is now created by concatenation (tpN -> base -> talairach) |
* The talairach.lta is now created by concatenation (tpN -> base -> talairach) * The brainmask is copied/mapped from the base (which is basically an OR of all TPs) * mri_ca_normalize has a new -long algorithm using the base aseg as init, fixing the bias towards smaller volumes in subcortical structures, as seen in 4.4. * mri_ca_label now uses the correct intensity scaling factors of the base. MartinReuter |
| Line 25: | Line 51: |
| First version of a working longitudinal stream, without optimizing each step. | First version of a working longitudinal stream, without optimizing each step. |
| Line 28: | Line 54: |
| * We have one '''-base run''' and '''-long runs''' for each TP * A base template (median) is created and used to initialize the longitudinal runs * The base is '''unbiased''' and can be viewed as an initial guess where things are |
* We have one '''-base run''' and '''-long runs''' for each TP. * A base template (median) is created and used to initialize the longitudinal runs. * The base is '''unbiased''' and can be viewed as an initial guess where things are. |
| Line 32: | Line 58: |
| * and [[mri_robust_template]] (unbiased robust template estimation) * Probabilistic fusion was added [[mri_fuse_segmentations]] to incorporate label information from the other TPs at a specific location * All TPs have to be processed cross sectionally (independently) first |
* and [[mri_robust_template]] (unbiased robust template estimation). * Probabilistic fusion was added [[mri_fuse_segmentations]] to incorporate label information from other TPs at a specific location. * All TPs have to be processed cross sectionally (independently) first. |
| Line 36: | Line 62: |
| == FS Version pre 4.4 == | MartinReuter == FS pre Version 4.4 == |
| Line 40: | Line 68: |
| The stream is described here: [[LongitudinalProcessingPreV4.4]]. | The old stream is described here: [[LongitudinalProcessingPreV4.4]]. |
| Line 43: | Line 71: |
| * One of the TPs was used to initialize the others (TP1 by default) * This lead to a bias wrt to TP1 * selecting a different TP for the initialization completely changed the results * A fix to init TP1 with itself improved things, but did not remove the bias |
* One of the TPs was used to initialize the others (TP1 by default). * This lead to a bias wrt to TP1. * Selecting a different TP for the initialization completely changed the results. * A quick fix to init TP1 with itself improved things, but did not remove the bias. |
Longitudinal Stream Change Log
A detailed description of the current longitudinal stream can be found here: LongitudinalProcessing
FS Version 5.1 Preview
- Will include the fixed binaries
- Will resample all time points to the base space (in the motion correction / conformalize step to avoid additional resampling)
- Will include revised/improved editing
We recommend to wait for this release. Initial studies show improved results due to the common voxel space. Furthermore, manual inspection is significantly simplified.
FS Version 5.0
!!! Do not use FS 5.0 for longitudinal processing, before updating some binaries !!!
- A non-related, last minute, incorrect 'bug-fix' in the transformation library introduced a new bug that breaks mri_robust_template, mri_robust_register (half way space).
You can find fixed versions of these binaries for your system here.
Also grab a modified recon-all from here that fixes the partial matching in the sanity check.
Improvements/Fixes:
- More robust base space estimation: keeps the base image centered
- Added sanity checks (if longitudinal TP is in base and if baseid is different from the time points id's)
No modifications to the data flow in the base or long stream.
FS Version 4.5
Improvements to the new longitudinal stream (modified base stream, fixes/improvements in the long stream).
Base Stream:
Only one template estimation on the norm.mgz files of all timepoints, to prevent difficulties with two possibly different template locations (T1 and norm as in 4.4).
A new block -base-init added to recon-all for the initial template creation, to be able to call later parts of the script individually without re-running the initialization everytime.
The -base-init block creates the maps from each TP to base and the norm_template (using the norm.mgz of all TPs).
- Then the orig.mgz of all TPs are resampled and averaged as orig/001.mgz to initialize the base run.
- The norm_template is used as brainmask (as it contains only brain).
- In -gcareg and -canorm the norm_template is used instead of the nu.mgz.
- This setup (with only one estimation) should lead to significant run time improvements especially with several time points.
Long Stream:
The talairach.lta is now created by concatenation (tpN -> base -> talairach)
- The brainmask is copied/mapped from the base (which is basically an OR of all TPs)
- mri_ca_normalize has a new -long algorithm using the base aseg as init, fixing the bias towards smaller volumes in subcortical structures, as seen in 4.4.
- mri_ca_label now uses the correct intensity scaling factors of the base.
FS Version 4.4
First version of a working longitudinal stream, without optimizing each step.
Main Differences:
We have one -base run and -long runs for each TP.
- A base template (median) is created and used to initialize the longitudinal runs.
The base is unbiased and can be viewed as an initial guess where things are.
New tools: mri_robust_register (symmetric registration)
and mri_robust_template (unbiased robust template estimation).
Probabilistic fusion was added mri_fuse_segmentations to incorporate label information from other TPs at a specific location.
- All TPs have to be processed cross sectionally (independently) first.
FS pre Version 4.4
These old versions of the longitudinal stream should not be used!
The old stream is described here: LongitudinalProcessingPreV4.4.
Short Info:
- One of the TPs was used to initialize the others (TP1 by default).
- This lead to a bias wrt to TP1.
- Selecting a different TP for the initialization completely changed the results.
- A quick fix to init TP1 with itself improved things, but did not remove the bias.
Original Author: MartinReuter