marsEqCurrently available automated segmentation tools only provide results for brain tissues (gray matter, white matter, cerebrospinal fluid (CSF)), have a limited FOV, and do not guarantee continuity and smoothness of tissues, which is crucially important for accurate modeling of transcranial electrical or magnetic stimulation. MARS is a tool that addresses these needs.

MARS is an extended toolbox for SPM software. It is developed based on the SPM8 toolbox “New Segment”, which is an implementation of the Unified Segmentation algorithm (Ashburner & Friston 2005), where the image intensity model and anatomical prior (atlas) are combined into a Bayesian inference framework. MARS further added the morphological constraints on neighboring tissue voxels to the framework, which is encoded by Markov Random Field (MRF). This MRF information is added into the updating equation of the E-step of the algorithm, not just as a post-processing step. Compared to the New Segment, the resulting segmentation is shown to be more smooth and with less discontinuities, while at the same time the accuracy is not significantly changed. For more details on the algorithm, see Yu Huang, Lucas C. Parra, Fully Automated Whole-Head Segmentation with Improved Smoothness and Continuity, with Theory Reviewed. PLoS ONE 10(5): e0125477.

MARS is always initialized by the New Segment. In other words, this toolbox will run New Segment in the first step to get the optimal estimates on the bias field, the registration parameters between atlas and the image(s) to be segmented, and the initial estimates of segmentation posteriors and parameters of Gaussian mixture model. MARS will then continues to update the segmentation posteriors and Gaussian mixture model by using both atlas (tissue probability map, TPM) and tissue correlation map (TCM), while keep the bias field and registration parameters untouched. The toolbox is backwards compatible with the New Segment toolbox, as users can disable the MARS update so that only New Segment will be run.

The TCM is a way of representing MRF constraint. It indicates the probability of co-occurrence of specific two tissues in two neighboring voxels. For example, white matter cannot be adjacent to the air, and thus the probability of observing white matter as a neighbor of air voxel is 0, and this is encoded in the TCM. We developed three categories of TCM: the local, global and regional. The local TCM gives different co-occurrence probabilities depending on the location of the voxels, while the global version uses the same value for all the voxels in the image. The local TCM has nice theoretical properties, but not practical for implementation, as it needs too much memory to store and compute; while the global one is too rough to include any local information. The regional TCM is a trade-off between the two: it encodes locality information such as the interior of white matter and the boundary between grey and white matters, while at the same time uses the same co-occurrence value for one specific locality. Therefore, the regional TCM is highly recommended, and local TCM is not recommended unless you have a 64-bit computer with at least 50 GB memory.

Download the Matlab code and TPM/TCM here.

The software can be accessed using this link.

This is the future home of the graphical user interface for the analytical solution of the current flow induced during transcranial current stimulation (tCS) in a multi-shell spherical model of the human head.

Reference: J. P. Dmochowski, M. Bikson, L. C. Parra, “The point spread function of the human head and its implications for transcranial current stimulation.”  Physics in Medicine and Biology.

In collaboration with Dr. Chris Rorden at Georgia Tech, an improved atlas (tissue probability map) has been developed to be used with New Segment in the free software SPM8.

To reduce the intensive manual labor needed in correction of segmentation errors generated from SPM8, a simple and fast Matlab script has been developed for automatic post-processing of the MRI segmentation from SPM8 New Segment.

For individualized tDCS modeling, a Matlab function for fast automated high-density virtual electrode placement following 10/10 international standard (or any user-defined electrode arrangement) has also been developed, to put electrodes on each individual head model.

A script for automated solving of tDCS finite element model in Abaqus under all possible bipolar electrode configurations is also attached, which can save manual labor when running FEM study in Abaqus for different electrode configurations.


For details, please refer to the following paper:

Yu Huang, Jacek P. Dmochowski, Yuzhuo Su, Abhishek Datta, Christopher Rorden, Lucas C. Parra, Automated MRI Segmentation for Individualized Modeling of Current Flow in the Human Head, Journal of Neural Engineering, 10(6):066004, 2013.

Download Matlab codes and atlas here.

The BONSIA and SPHERES free modeling tools are live and FREE. The most advanced AND simple to use freely available modeling packages to help design TES montages including for tDCS, tACS, tRNS, and more. Check them out and give us your feedback on how to improve!

The paper explaining the principles is here. (PDF)

BONSAI is online here

SPHERES is online here

Screen Shot 2014-05-09 at 10.16.28 AM