Reconstruction


The reconstruction step processes an SRC file from Step T1 to generate the FIB file, which can be further used in fiber tracking.

Step T2: Open SRC File(s)

Click on [Step T2 Reconstruction] button in the main project window to select one or multiple SRC files.

DSI Studio will present a reconstruction window as shown in the figure to the right.

Tip: you can select multiple SRC files here, and DSI Studio will reconstruct each of them respectively. Even if you have additional preprocessing such as flipping images or rotating them to the MNI space, DSI Studio will apply the same procedure to each of the SRC files.

image

Visual Quality Inspection (Optional)

Switch the first tab [Source Images].

  1. Check eddy current distortion and motion: click on the b-table list on the left and use the keyboard arrow key to scroll down the list. If the brain is distorted across DWI, it is likely caused by the eddy current. Subject head movement can also be seen when scrolling through the DWIs.

  2. Check bad slices: click on the [Show bad slices] button (available after 10/2/2019 version) and DSI Studio will mark bad slices and their belonging DWI volume in red (see below). I recommend using sagittal view to check the bad slices (click on the sagittal view button on the right)

  1. Checking other problems:

Step T2a: Setup a Mask

image

The purpose of the mask is to filter out the background region, increase the reconstruction efficacy, and facilitate further visualization. In the mask selection window, there is a track bar at the bottom for slice selection. The mask can be generated using several built-in functions provided by DSI Studio:

Recommendation steps include thresholding, smoothing, and defragment. You can open/save the mask to a txt file or nifti file.

Preprocessing (Human Data)

If you have reverse-phase encoding acquisition:

  1. On the top menu, click on [Corrections][TOPUP/EDDY]
  2. Select the reverse phase encoding acquisition (nii.gz or src.gz). If you are not sure how to identify it, please refer to the video tutorial
  3. DSI Studio will call FSL’s topup and eddy to handle susceptibility artifact and eddy current distortion.

If you do not have reverse-phase encoding acquisition:

  1. On the top menu, click on [Corrections][EDDY]
  2. DSI Studio will call FSL’s eddy to handle susceptibility artifact and eddy current distortion.

This preprocessing may take a long time (in hours)

Preprocessing (Animal Data)

  1. Correct image orientation

    Most animal scans are likely to be acquired at an orientation different from the animal template in DSI Studio. You will need to rotate the image volume to match the template so that the atlas function can be used in Step T3. The following is the template orientation for mouse/marmoset/rhesus data.

    Moving the sliding bar to the “right” should go to the “top” of the brain. For rhesus and marmoset scans, the slices should be oriented in axial view, which top slice position on the right of the sliding bar.

    image

    Use [Edit][Image Flip] or [Image Swap] to make sure your image volume matches that of the template. Note that either “Image Flip” or “Image Swap” will create a mirrored image, and thus you may need either 2 or 4 flip/swap operations to avoid creating mirrow images. you can always add an additional left-right flip to ensure an even number of operations.

    For correcting oblique angles, use [Edit][Rotate…]

  2. remove background signals or crop image volume (optional, recommended)

    Animal scans often include signals from tissues outside the brain, and they can cause a problem in the atlas or template-related analysis.

    This can be tackled by using [Edit][Erase Background Signals] to eliminate the signals outside the mask.

    Also use [Edit][Crop Background] to reduce image volume,

  3. Reduce FOV (optional)

    After removing background signals, there could be a large empty background. The FOV can be reduced to fit the mask by [Edit][Crop background].

  4. Make isotropic (optional)

    Many animal scans may have a non-isotropic resolution, and sometimes it can cause registration problems.

    The data can be interpolated using [Edit][Make isotropic]

  5. Check [Check b-table] at Step T2b(2)

    If your analysis does not consider fiber orientation (e.g. just analyze FA, ADC…etc.) You can ignore the b-table flipping problem. Otherwise, you will need to pay attention to the reconstructed fiber orientations.

    The b-table from animal scans are often flipped in the image dimension, and DSI Studio can help correct it if users check the [Step T2b(2)][Check b-table]. This b-table checking function will examine a total of 24 different flip and swap conditions and figure out the one that gives the best fiber coherence. However, the b-table checking function may fail if the SNR is not good, and the reconstructed fiber orientations may be flipped. If this happens, you may need to uncheck this function, and manually flip or swap b-table using the [B-table][flip] function.

Step T2b(1): Select a Reconstruction Method

Diffusion Tensor Imaging (DTI)

DTI analysis will generate one fiber orientation per voxel and associated anisotropy and diffusivity measure. GQI and QSDR are recommended because DSI Studio will also calculate DTI metrics for GQI and QSDR. The only condition to choose DTI here is when GQI fails to produce correct results.

Generalized Q-sampling Imaging (GQI)(Recommended)

GQI is a model-free method for resolving fiber orientations and quantifying the anisotropy of diffusing water. DSI Studio will also calculate DTI metrics even if GQI is used.

The diffusion sampling length ratio controls GQI’s sensitivity, and it is highly recommended to optimize the ratio using a preliminary scan following these steps:

  1. Reconstruct FIB files using ratios of 0.3, 0.4, 0.5, …2.0 etc.
  2. Open one FIB file in [Step T3 Fiber Tracking]
  3. Switch to the coronal view [View Coronal View] and adjust the [Tracking Threshold] option at the right upper corner under [Step T3c: Options][Tracking Parameters] so that the value covers only the white matter:

image

  1. Locate non-crossing regions (e.g. mid corpus callosum marked by the light red circle) and crossing regions (e.g. lateral corpus callosum marked by light blue). The optimal value should resolve crossing patterns in crossing regions (light blue) and still have clean fiber direction at mid corpus callosum (light red).

The general principle is to “select the highest L that has an acceptable amount of spurious fibers.

image

Q-Space Diffeomorphic Reconstruction (QSDR)

QSDR [11] is the MNI version of GQI. Choose QSDR if you would like to have all the following analyses in the MNI space (e.g. correlation tractography and related connectometry analysis).

To carry out other images modalities to the MNI space with QSDR, click on the [Attach Images…] button at the right side of the output resolution box. To carry out other images modalities to the MNI space with QSDR, click on the [Attach Images…] button at the right side of the output resolution box. If you have T1W-based ROI to be transformed along with the T1W, add the T1W first before adding the ROI. If you have T1W-based ROI to be transformed along with the T1W, add the T1W first before adding the ROI.

The QSDR reconstruction requires the assignment of a template (e.g. human, monkey, rat, mouse). Please make sure that your data match the template.

Step T2b(2): Specify Parameters

Output ODF

Checking this check box will export all ODF information and allow DSI Studio to run connectometry analysis or visualize the ODFs for inspection. Note that enabling this option will greatly increase the size of the .fib file.

Check b-table

DSI Studio allows for checking and correcting whether the b-table is flipped and/or swapped in x, y, or z-direction. The function is realized by using a fiber coherence index to check which settings a better result [14]. The checkbox for “check b-table” is located on the bottom of the reconstruction dialog after selecting the reconstruction method (e.g. DTI, GQI, …etc.).

*If the diffusion signals are corrupted, this b-table checking may result in a random configuration that can be easily identified from the file name of the FIB file.

Advanced Options: Scheme balance

Check this option if you are using a shell acquisition to make sure that the orientational sampling is homogeneous.

Advanced Options: No high b for DTI

Check this if your tensor estimation uses only b-values lower than 4,000.

Advanced Options: Other output metrics

Specify the name of the metrics (separated by comma) to be included in the output FIB file: