Difference between revisions of "Projects:DTI DWI QualityControl"

From NAMIC Wiki
Jump to: navigation, search
 
(17 intermediate revisions by 2 users not shown)
Line 4: Line 4:
 
=  Diffusion Tensor and Diffusion Weighted Imaging Quality Control =
 
=  Diffusion Tensor and Diffusion Weighted Imaging Quality Control =
  
DWIs data suffer from inherent low SNR, overall long scanning time of multiple directional encoding with correspondingly large risk to encounter several kinds of artifacts. These artifacts can be too severe for a correct and stable estimation of the diffusion tensor. Thus, a quality control (QC) procedure is absolutely necessary for DTI studies. We are developing a framework, called DTIPrep, for assessing and correcting DWIs and DTI.
+
[[Image:DTIPRep GUI.png|400px|thumb|right| DWI-based QC results using DTIPrep through three steps: 1) converting dicom to  nrrd format of DWI image, 2) loading the protocol and running the software and 3) potential needed of visual checking and final saving. In this example, gradient #11 suffers from intensity artifact and is excluded. The sphere shows 3D view of gradients distribution before and after running DTIPrep in blue and green colors respectively. ]]
 +
 
 +
[[Image:GradientsDistribution.png|400px|thumb|right|3D view of gradients distribution with different b-values. ]]
 +
 
 +
As theoretical work characterizing DTI grows, it is essential to increase its practical usability from a clinical environment perspective. Inherently, DWI images suffer from a vast variety of artifacts and the acquisition time for diffusion MRI is longer than conventional MRI due to the need for multiple acquisitions to obtain directionally encoded Diffusion Weighted Images (DWI). This leads to increased motion artifacts and reduced signal-to-noise ratio (SNR). Therefore, in a clinical environment, this imaging technique needs additional processes such as appropriate QC assessment methods to increase its practical usability. We are developing a framework, called DTIPrep, for assessing and correcting DWIs and DTI.
  
 
= Description =
 
= Description =
 +
 +
== Current framework for DWI QC ==
 +
DTIPrep is the first comprehensive and fully automatic pre-processing tool for DWI and DTI quality control can provide a crucial piece for robust DTI analysis studies. The protocoling, reporting, visual controlling and data correction capabilities are used to produce high consistence and inter-rater reliable QC results. This framework is organized by pipeline steps include: 1) Dicom to NRRD converting, 2) image info checking, 3) diffusion information checking, 4) rician LMMSE noise filter, 5) slice-wise intensity checking, 6) interlace-wise intensity checking, 7) Averaging baseline images, 8) Eddy current and motion correction, 9) gradient-wise checking of residual motion/deformations, 10) joint rician LMMSE noise filter, 11) brain masking, 12) DTI computing, 13) dominant direction artifact (vibration artifact) checking, 14) optional visual checking and 15) simulation-based bias analysis. 
  
  
[[Image:Screenshot.png|400px|thumb|right| Diffusion Weighted Imaging and Diffusion Tensor Imaging Quality Control_DTIPrep]]
+
== Entropy Based Diffusion Imaging Quality Control  ==
 +
We have proposed new QC step for detecting drop-out signal intensities which can be caused by mechanical vibration artifacts. This step detects and potentially removes these residual artifacts that are not commonly detected in the individual DWIs. The artifacts appear in color-FA images in either widespread or local dominant direction. In order to detect such artifacts, we proposed a new approach via the entropy of the Principal Direction (PD) histogram computed over the major region of the image (e.g. the full brain). Given a prior knowledge of expected entropy values for acceptable scans, the quality of the DTI image is categorized into acceptable, suspicious and highly suspicious/rejection categories using calculated the standard scores.
  
[[Image:753361684_3.png|400px|thumb|right|3D view of gradients before and after Quality Control procedures]]
 
  
== Current framework for DWI QC ==
+
[[Image:Vibration artifact.png|800px|thumb|center| Vibration artifacts. a) Artifact free scan. Top-left, a representative axial slice of the color-FA map. Top-right, spherical histogram of the PD distribution within the entire brain. Bottom-right, tractography of genu and splenium of the corpus callosum. Bottom-left, genu tract in more detail. (b) Vibration artifacts may manifest as localized (prefrontal region for this example) signal-loss in the DWI image or as dominant L-R (red) direction. (c) Vibration artifact in the absence of localize DWI signal loss. Spherical viewpoints chosen to show locations of highest histogram frequency.]]
DTIPrep is the first comprehensive and fully automatic pre-processing tool for DWI and DTI quality control can provide a crucial piece for robust DTI analysis studies. The protocoling, reporting, visual controlling and data correction capabilities are used to produce high consistence and inter-rater reliable QC results. This framework is organized by pipeline steps include: 1) image info checking, 2) diffusion info checking, 3) Slice-wise intensity checking, 4) Interlace-wise intensity checking, 5) Averaging baseline images, 6) Eddy current and motion correction, 7) gradient-wise checking of residual motion/deformations, 8) Computing DTI measurements and saving and 9) Optional visual QC.
+
 
 +
For our correction step, we employed a simple, iterative leave-one-out-strategy over all individual DWI images by recomputing DTI images and correspond- ing entropies. At each iteration, the DWI with maximal improvement is removed and all leave-one-out entropies are recomputed. This process is continued until either the z-score is in acceptable range or a maximum threshold for exclusion is reached.
 +
 
 +
[[Image:Correction.png|700px|thumb|center| The correction results show visible improvement in contrast within the cingulum and fornix tracts (left) and fiber tractography of splenium (right). ]]
  
== Ongoing extension using DTI QC ==
 
  
In our ongoing research, we are extending our QC procedures to include a QC step based on the derived DTI data. This step detects and potentially removes residual artifacts that are not commonly detected in the individual DWIs.  In our large scale population studies, we observed several such artifacts, most specifically an artifact of "dominating direction".Our experiments show that the residual artifacts presence after DTIPrep DWIs QC can be detected and corrected using knowledge of DTI. We introduce new approach in DTI QC including the ''detection'' and the ''corrections'' steps. Our approach use new entropy-based benchmark comes up with Principal Directions (PDs) histogram implemented by multi-level subdivision icosahedron within different regions of brain. Given training of the measurement, the quality of DTI is categorized into acceptable, suspicious and unacceptable groups using calculated the standard scores.
+
== Ongoing Research about Simulation-based Bias Analysis==
  
[[Image:Diagram2_copy.png|400px|thumb|right| Visualization of dominating direction artifact using the color FA DTI image (left column) and the corresponding spherical histogram of the principal directions within the brain. Top: Example of an acceptable DTI dataset. Middle: Example of a "green" direction (anterior-posterior) dominating artifact. Bottom: "Red" direction (left-right) artifact.]]
+
It is important to note that there are some pitfalls associated with these QC approaches in the result of QC after correcting these artifacts. The correction processes modify a configuration of gradient sampling from a scan either by changing the gradients directions or excluding individual DWI’s along a subset of the gradients due to artifacts. Conventionally, QC procedures exclude artifact-affected DWIs from subsequent computations leading to a cleaned, reduced set of DWIs, called DWI-QC. Often, a rejection threshold is heuristically/empirically chosen above which the entire DWI-QC data is rendered unacceptable and thus no DTI is computed. In this work, we propose a simulation-based automatic DTI QC to assess the resulting tensor properties from DWI-QC techniques. We define two error metrics based on the directional distribution of bias for FA and PD. These metrics can provide a promising benchmark for post-QC assessment of the remaining DWIs. Based on our simulation results, we introduce rejection metrics (automatic thresholds) with respect to magnitude and directional distribution of bias for FA and PD.
  
[[Image:AutismTable3_copy.png|400px|thumb|right| The z-score of orientational entropies for autism sub- jects. The detected images are indexed from 13 and all suffer from color artifac.]]
+
In experimental results, we applied our method on acquisition schemes and also individual scans post-QC. These results show that the proposed rejection metrics provide an effective assessment of post-QC individual scan and also acquisition protocols. Furthermore, our results confirm that higher degrees of uniformity in the sampling gradients results in lower overall bias. Thus, determination of diffusion properties with minimal error requires an evenly distributed gradient directions before and after QC. This method will be incorporated in DTIPrep, used for QC of DWI/DTI data.
  
We employ the correction step by excluding gradients which have the most contribution in the artifacts. We continue ex- cluding gradients till the z-score of whole brain of updated image will be close enough to the trained.
+
[[Image:Mc schemes.png|500px|thumb|center| Bias analysis of different DWI schemes via MC simulation. (a-c) Gradient direction distribution for 3 acquisition schemes: 42-direction quasi-uniform, Phillips 32-direction non-uniform, and 6-direction uniform. (d- f) Estimated error distribution in PD computation. (g-i) Estimated error distribution in FA computation as a percentage of true FA. Number of iteration 200,000, true FA value of 0.4 and SNR = 10. ]]
  
[[Image:Correction_Figure_copy.png|400px|thumb|right| The correction step result (the corrected image in right side) shows much more improvement in observing cc and fx tracts. The FA profiles of the genu and splenium are shown in bottom ( blue color: corrected image)]]
+
[[Image:Mc scans.png|600px|thumb|center| Scan-specific bias analysis via MC simulation. (a-b) The gradient direction distribution after excluding 20% of the total number of gradients due to artifacts, in two configurations: non-clustered and clustered exclusion (blue: excluded gradients, green: included gradients). (c) Estimated error distribution of PD computation given the gradient sampling schemes in (a) and (b). (d) Estimated error distribution in FA computation, shown as percentage of true FA. MC simulation was performed for 200,000 iterations for this experiment, with a true FA value of 0.4 and SNR = 10.]]
  
 
= Publications =
 
= Publications =

Latest revision as of 18:53, 2 January 2013

Home < Projects:DTI DWI QualityControl
Back to UNC Algorithms


Diffusion Tensor and Diffusion Weighted Imaging Quality Control

DWI-based QC results using DTIPrep through three steps: 1) converting dicom to nrrd format of DWI image, 2) loading the protocol and running the software and 3) potential needed of visual checking and final saving. In this example, gradient #11 suffers from intensity artifact and is excluded. The sphere shows 3D view of gradients distribution before and after running DTIPrep in blue and green colors respectively.
3D view of gradients distribution with different b-values.

As theoretical work characterizing DTI grows, it is essential to increase its practical usability from a clinical environment perspective. Inherently, DWI images suffer from a vast variety of artifacts and the acquisition time for diffusion MRI is longer than conventional MRI due to the need for multiple acquisitions to obtain directionally encoded Diffusion Weighted Images (DWI). This leads to increased motion artifacts and reduced signal-to-noise ratio (SNR). Therefore, in a clinical environment, this imaging technique needs additional processes such as appropriate QC assessment methods to increase its practical usability. We are developing a framework, called DTIPrep, for assessing and correcting DWIs and DTI.

Description

Current framework for DWI QC

DTIPrep is the first comprehensive and fully automatic pre-processing tool for DWI and DTI quality control can provide a crucial piece for robust DTI analysis studies. The protocoling, reporting, visual controlling and data correction capabilities are used to produce high consistence and inter-rater reliable QC results. This framework is organized by pipeline steps include: 1) Dicom to NRRD converting, 2) image info checking, 3) diffusion information checking, 4) rician LMMSE noise filter, 5) slice-wise intensity checking, 6) interlace-wise intensity checking, 7) Averaging baseline images, 8) Eddy current and motion correction, 9) gradient-wise checking of residual motion/deformations, 10) joint rician LMMSE noise filter, 11) brain masking, 12) DTI computing, 13) dominant direction artifact (vibration artifact) checking, 14) optional visual checking and 15) simulation-based bias analysis.


Entropy Based Diffusion Imaging Quality Control

We have proposed new QC step for detecting drop-out signal intensities which can be caused by mechanical vibration artifacts. This step detects and potentially removes these residual artifacts that are not commonly detected in the individual DWIs. The artifacts appear in color-FA images in either widespread or local dominant direction. In order to detect such artifacts, we proposed a new approach via the entropy of the Principal Direction (PD) histogram computed over the major region of the image (e.g. the full brain). Given a prior knowledge of expected entropy values for acceptable scans, the quality of the DTI image is categorized into acceptable, suspicious and highly suspicious/rejection categories using calculated the standard scores.


Vibration artifacts. a) Artifact free scan. Top-left, a representative axial slice of the color-FA map. Top-right, spherical histogram of the PD distribution within the entire brain. Bottom-right, tractography of genu and splenium of the corpus callosum. Bottom-left, genu tract in more detail. (b) Vibration artifacts may manifest as localized (prefrontal region for this example) signal-loss in the DWI image or as dominant L-R (red) direction. (c) Vibration artifact in the absence of localize DWI signal loss. Spherical viewpoints chosen to show locations of highest histogram frequency.

For our correction step, we employed a simple, iterative leave-one-out-strategy over all individual DWI images by recomputing DTI images and correspond- ing entropies. At each iteration, the DWI with maximal improvement is removed and all leave-one-out entropies are recomputed. This process is continued until either the z-score is in acceptable range or a maximum threshold for exclusion is reached.

The correction results show visible improvement in contrast within the cingulum and fornix tracts (left) and fiber tractography of splenium (right).


Ongoing Research about Simulation-based Bias Analysis

It is important to note that there are some pitfalls associated with these QC approaches in the result of QC after correcting these artifacts. The correction processes modify a configuration of gradient sampling from a scan either by changing the gradients directions or excluding individual DWI’s along a subset of the gradients due to artifacts. Conventionally, QC procedures exclude artifact-affected DWIs from subsequent computations leading to a cleaned, reduced set of DWIs, called DWI-QC. Often, a rejection threshold is heuristically/empirically chosen above which the entire DWI-QC data is rendered unacceptable and thus no DTI is computed. In this work, we propose a simulation-based automatic DTI QC to assess the resulting tensor properties from DWI-QC techniques. We define two error metrics based on the directional distribution of bias for FA and PD. These metrics can provide a promising benchmark for post-QC assessment of the remaining DWIs. Based on our simulation results, we introduce rejection metrics (automatic thresholds) with respect to magnitude and directional distribution of bias for FA and PD.

In experimental results, we applied our method on acquisition schemes and also individual scans post-QC. These results show that the proposed rejection metrics provide an effective assessment of post-QC individual scan and also acquisition protocols. Furthermore, our results confirm that higher degrees of uniformity in the sampling gradients results in lower overall bias. Thus, determination of diffusion properties with minimal error requires an evenly distributed gradient directions before and after QC. This method will be incorporated in DTIPrep, used for QC of DWI/DTI data.

Bias analysis of different DWI schemes via MC simulation. (a-c) Gradient direction distribution for 3 acquisition schemes: 42-direction quasi-uniform, Phillips 32-direction non-uniform, and 6-direction uniform. (d- f) Estimated error distribution in PD computation. (g-i) Estimated error distribution in FA computation as a percentage of true FA. Number of iteration 200,000, true FA value of 0.4 and SNR = 10.
Scan-specific bias analysis via MC simulation. (a-b) The gradient direction distribution after excluding 20% of the total number of gradients due to artifacts, in two configurations: non-clustered and clustered exclusion (blue: excluded gradients, green: included gradients). (c) Estimated error distribution of PD computation given the gradient sampling schemes in (a) and (b). (d) Estimated error distribution in FA computation, shown as percentage of true FA. MC simulation was performed for 200,000 iterations for this experiment, with a true FA value of 0.4 and SNR = 10.

Publications

Key Investigators

  • UNC Algorithms: Mahshid Farzinfar, Zhexing Liu, Martin Styner, Clement Vachet
  • Utah Algorithms: Tom Fletcher, Ross Whitaker, Guido Gerig, Sylvain Gouttard

Links