Graph-based segmentation of signal voids in time series of diffusion-weighted images of musculature in the human lower leg
Martin Schwartz1,2, Günter Steidle1, Petros Martirosian1, Bin Yang2, and Fritz Schick1

1Section on Experimental Radiology, Department of Radiology, University of Tuebingen, Tuebingen, Germany, 2Institute of Signal Processing and System Theory, University of Stuttgart, Stuttgart, Germany

Synopsis

The segmentation of signal voids, which occur in time-series of single-shot diffusion-weighted images, is important for an accelerated evaluation providing larger studies on this phenomenon. The proposed segmentation is based on a two-stage detection and segmentation approach, which utilizes a graph-based representation with random walker optimization. It was demonstrated that the presented method enables a fast and accurate segmentation of signal voids in time-series of diffusion-weighted images.

Purpose

In time-series of single-shot diffusion weighted (DW) images signal voids in different muscle groups have been observed 1. Due to the random appearance of these signal voids in the time domain and a large number of images, an automatic segmentation is desired to expedite the time-consuming evaluation and thus enable larger studies for revealing the underlying processes of this phenomenon. To prevent an underestimation of the size of the signal voids caused by ignoring regions with partial volume effects, a two-class detection and segmentation approach is proposed. Furthermore, a-priori knowledge about the location of vessels, which are responsible for signal fluctuations due to blood flow and pulsation 1, is directly incorporated into the segmentation instead of using an extra processing step 2.

Methods

The DW images were acquired on a 3 T MR scanner (Magnetom Skyra, Siemens Healthcare, Erlangen, Germany) with a stimulated echo DW EPI sequence (matrix size: 64 x 64, FoV = 200 x 200 mm², TE = 31 ms, TR = 500 ms, BW = 2004 Hz/px, TM = 145 ms, 500 repetitions and 15 channel Tx/Rx knee coil). The segmentation procedure is depicted in Fig. 1 and was implemented in MATLAB® (The Mathworks, Inc., USA). In the first step, bias-field and long-term variations were estimated and corrected by the BCFCM 3 (bias-corrected fuzzy c-means) and RLOESS 1,4 (robust local regression with a span of 5 %) algorithm. Repetitions with high signal energy were averaged to calculate a reference image without signal voids as input image for the bias-field correction algorithm. The BCFCM algorithm was utilized to get the bias-field of the muscle tissue and to determine a tissue mask to reduce the computational load in subsequent calculations. After these processing steps, the image xi is converted in a time-difference representation by using two temporally adjacent images xi-1 and xi+1 according to (1).

$$I_t(x_i)=x_i-\frac{x_{i-1}+x_{i+1}}{2} \quad \quad \quad \quad (1)$$

This simple high pass filtering enhances the temporal changes in signal intensity and emphasizes the starting points and endpoints of the signal void events. The event detection step is based on a standard two class Fuzzy C-Means algorithm 5,6 with an additional grey-level based rejection term to suppress false positive detections in noisy images. For the segmentation the variation corrected DW images are transformed in an undirected graph-based representation, while the edge weighting is derived from the grey-level of the connected voxels 7, the distance to the seed points and from the time-difference representation. To prevent detections at pulsation sources like vessels, the vessel locations are estimated from a flow compensated GRE image and are integrated as an additional weighting term in the weighting function. For the intra-slice weighting (in one DW image), an 8-point, and for the inter-slice weighting (over time), a 2-point neighborhood connection is utilized as shown in Fig. 2. For an automatic seed point placement, the derivation of an Euclidean distance map is calculated 8,9 (Fig. 3 c) and d)) to iteratively shift the equidistantly placed seed points of the “no event” class away from the event seed points. This approach determines a region between already segmented points (from the event detection) and the shifted seed points to reduce again the computational load (Fig. 3 e)). With the weighting between the voxels, a graph-based optimization problem had to be solved for the voxel classification. This optimization problem is solved by the random walker algorithm 7 (Fig. 3 f)), which reformulates the classification task into the solution of a linear system based on a combinatorial Dirichlet problem.

Results

In Fig. 3 partial results along the entire whole segmentation workflow are shown. Segmentation steps from the original DW image in a) to the final segmentation result in f) are included. Three original DW images with overlaid segmentation results are depicted in Fig. 4 showing signal voids in the m. gastrocnemius medialis and in the m. soleus. It can be seen that the vessel region between the m. soleus and the m. tibialis posterior as well as the fibula is not considered as signal void (marked in the middle and right DW image).

Conclusion

The proposed segmentation procedure enables a fast and accurate segmentation of signal voids in musculature of the human calf with considerations of the partial volume affected as well as vessel and bone regions. Moreover, the segmentation can be extended to a 3D + t approach for the segmentation of time series data of 3D or multi-slice diffusion weighted images by an extension of the graph construction step.

Acknowledgements

No acknowledgement found.

References

[1]: Steidle et al., NMR Biomed.:28(7); [2]: Schwartz et al., Proc ESMRMB 2015; [3]: Ahmed et al., IEEE Trans. Med. Imaging 2002:21(3); [4]: Cleveland, W., Journal of the American Statistical Association 1979:74(368); [5]: Dunn, J., Journal of Cybernetics 1974; [6]: Bezdek, J., Plenum Press New York 1981; [7]: Grady, L., IEEE Trans. Pattern Analysis and Machine Intelligence 2006:28(11); [8]: Maurer et al., IEEE Trans. On Pattern Analysis and Machine Intelligence 2003:25(2); [9]: Kohlmann et al., International Journal on Computer Assisted Radiology and Surgery 2015:10(4)

Figures

Fig. 1: Two-stage detection and segmentation procedure with an incorporation of a-priori knowledge of the vessel locations.

Fig. 2: Graph construction for the segmentation task with intra- and inter-slice weighting to yield a time series segmentation (2D + t).

Fig. 3: Segmentation workflow - a) original DW image, b) highlighted event detection, c) Euclidean distance map with gradient vectors, d) original DW image with gradient vectors, e) resulting automatic seeding and f) segmentation result.

Fig. 4: Original DW images with overlaid segmentation result. Marked vessel and fibula are not segmented as voids.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1923