“…We can estimate ε i ( d,t ij ) by using
and concatenate them into a vector ε̂ ( d ) = ( ε̂ i ( d , t ij ) : i = 1, …, n ; j = 1, …, m i ) T for each voxel. Specifically, we use a penalized likelihood approach with an L 1 penalty function for a Gaussian mixture model to cluster all residual vectors { ε̂ ( d ) : d ∈ 𝒟} into K homogeneous regions [Pan and Shen, 2007, Huang et al, 2015]. …”