Multiple sclerosis white matter (WM) lesions can affect brain tissue volume measurements of voxel-wise segmentation methods if these lesions are included in the segmentation process. Several authors have presented different techniques to improve brain tissue volume estimations by filling WM lesions before segmentation with intensities similar to those of WM. Here, we propose a new method to refill WM lesions, where contrary to similar approaches, lesion voxel intensities are replaced by random values of a normal distribution generated from the mean WM signal intensity of each two-dimensional slice. We test the performance of our method by estimating the deviation in tissue volume between a set of 30 T1-w 1.5 T and 30 T1-w 3 T images of healthy subjects and the same images where: WM lesions have been previously registered and afterwards replaced their voxel intensities to those between gray matter (GM) and WM tissue. Tissue volume is computed independently using FAST and SPM8. When compared with the state-of-the-art methods, on 1.5 T data our method yields the lowest deviation in WM between original and filled images, independently of the segmentation method used. It also performs the lowest differences in GM when FAST is used and equals to the best method when SPM8 is employed. On 3 T data, our method also outperforms the state-of-the-art methods when FAST is used while performs similar to the best method when SPM8 is used. The proposed technique is currently available to researchers as a stand-alone program and as an SPM extension

Sergi Valverde holds a FI-GDR2013 grant from the Generalitat de Catalunya