BrainVoyager v23.0

Intensity Inhomogeneity Correction

Background

Intensity inhomogeneities in MR images can substantially reduce the accuracy of segmentation and registration. Despite advances of correcting spatial intensity inhomogeneities of modern scanner software, the advent of multichannel phased array coils and 7T+ scanners have increased (again) the importance of this problem for post-scan processing since images may exhibit substantial intensity inhomogeneities across space (see image on left side in figure above). This poses a problem for segmentation tools that are based on region growing approaches since white matter voxels at one position of the image might have the same intensity as grey matter voxels at other locations in the image. One well-performing method of intensity inhomogeneity correction (IIHC) that compares well with other methods (e.g. Sled et al., 1997) is based on a "surface fitting" approach. In this approach, low-order polynomials are used to model low-frequency variations across 3D image space; these polynomials are fitted to a subset of voxels that have been labeled as belonging to white matter. After estimating the low-frequency intensity fluctuations (producing a "bias field", see image on lower right side above), they are removed from the data producing voxels with much more homogeneous intensities (see image on upper right side above) that improve visualization and also are better starting points for subsequent segmentation tools.

In BrainVoyager QX 2.1 and earlier it was necessary to select presumed white matter voxels using region growing in specific intensity ranges. This manual "pre-segmentation" step is often time-consuming since several attempts are usually required until good intensity range values are found; furthermore, the obtained quality depends on the experience of the user. Version 2.2 of BrainVoyager QX introduced a new IIHC tool that contains as its core the same polynomial "surface fitting" approach but includes new tools that allow to run IIHC without user intervention. Since bias field estimation works more precise when it's focused on brain tissue only, a background cleaning step followed by a brain extraction step are useful preparator steps. More specifically, BrainVoyager QX executes the following steps during automatic intensity inhomogeneity correction:

Running Automatic IIHC

The described steps are available in the 16 Bit 3D Tools dialog (see below) that can be invoked using the Volumes > Inhomogeneity Correction, V16 Tools menu item. Note that the .v16 data set belonging to the .vmr data set (i.e. a file with the same name as the current VMR but with extension .v16) is automatically loaded if it can be located in the same folder, otherwise the Load .V16 button can be used to load an appropriate V16 file.

In most cases, pressing the GO button (see snapshot above) in the Intensity inhomogeneity correction (IIHC) field is sufficient to create a homogeneous data set in a few seconds. Starting the automatic IIHC tool triggers the sequence of processing steps described above resulting usually in a segmented, homogeneous brain data set (see top right image in figure of Background section above and right image in snapshot below). The preparatory steps option, Include mask generation, can be turned on or off; it is automatically checked in case that this step has not yet been performed (i.e. if no "_MaskBrain.vmr" file is availalbe in the folder of the VMR/V16 file). If the Map GM and WM peaks to standard VMR values is turned on (default), the grey matter and white matter intensities will be centered around intensiy values of 100 and 160, respectively. Detection of white matter can be influenced by changing the Tissue range threshold and Intensity threshold (percentile of low intensity homogeneous blocks) values, for details, see description above. The order of the Legendre polynomials can be set in the Order of polynom for fit option (default: 3). When the Save resulting IIH corrected .V16 data to disk option is checked (default), the resulting IIHC corrected V16 file will be saved to disk.

The snapshot above shows the result of IIHC for a 3D data set with only rather small inhomogeneities (left), nonetheless the resulting data set indicates a gain in white and grey matter separation and within-tissue homogeneity. Note that BrainVoyager stores the resulting file under the name of the original VMR/V16 data set but with the substring "_IIHC" (see title bar of window on right side in snapshot above). In addition to an inhomogeneity corrected VMR/V16 data set, the automatic procedure also generates representations of the generated bias fields. A bias field represents the fitted model of the low-frequency intensity variations of the automatically labeled white matter voxels over space (see bias field image below).The estimated bias fields are also stored to disk using the name of the original file plus the substring "_IIH-BiasField[nr]" where the "[nr]" part refers to the iteration of the IIHC process. Since the first calculated bias field highlights the most severe inhomogeneities, it is shown as default in the secondary VMR of the corrected VMR view even when running multiple iterations. Bias field visualizations from subsequent iterations can be loaded as separate (primary) VMRs or loaded in the secondary VMR view. To visualize the produced bias field, the usual 1st / 2nd volume toggle key (F8) can be used or the Show primary VMR / Show secondary VMR options in the Two VMR display options field in the Spatial Transf tab of the 3D Volume Tools dialog (see snapshot below).

While a single run of the automatic IIHC process leads to very good results, the process is repeated once as default (No. of cycles: 2) in order to further increase the resulting homogeneity. More than two iterations do usually not lead to significant further improvements as can be assessed with the displayed intensity histograms that are generated for each iteration (see figure below). Each histogram shows the result before (blue line) and after (yellow line) application of one cycle through the processing pipeline. In the plot on the left side below (iteration 1), the blue line indicates that white and grey matter voxel intensities overlap substantially since there is only one broad peak in the histogram; the yellow line indicates the successful application of the automatic inhomogeneity correction process clearly revealing two peaks corresponding to grey and white matter. After the second iteration (snapshot below, right side) not much further improvement is visible.

While IIHC can be executed automatically clicking the GO button, all involved steps can also be used in isolation, including background cleaning (Clean button), brain extraction (Extract Brain button), white matter detection (Detect WM button) and bias field estimation and removal (Correct IIH button). Separate execution with non-standard parameters may be useful in rare cases where the automatic process does not lead to satisfactory results; while the default setting of these parameters are usually working well, it might be useful to change them in specific situations.

Tip: For visualization purposes, within-tissue homogeneity can be further enhanced by running the Sigma filter from the Filter, smoothing field in the Segmentation tab of the 3D Volume Tools dialog. If one plans to run (advanced) segmentation tools subsequently, this step should not be performed since it is executed as part of automtic segmentation routines.

References

Hou, Z, Huang, S, Hu, Q, Nowinski, WL (2006). A fast and automatic method to correct intensity inhomogeneity in MR brain images. In: R. Larsen, M. Nielsen, and J. Sporring (Eds.): MICCAI 2006, LNCS 4191, pp. 324–331, Springer: Berlin, Heidelberg.

Sled (1997). A comparison of retrospective intensity non-uniformity correction methods for MRI. Lecture Notes in Computer Science.


Copyright © 2023 Rainer Goebel. All rights reserved.