Based on 4 MR sequences (multi-parametric MRI scans), classify each voxel within the tumor volume into high and low categories via kmeans clustering. With 16 (24) signal combinations across the 4 sequences (i.e. 0000-1111), every voxel in the tumor volume can be identified uniquely. The resultant habitat map shows the spatial heterogeneity within tumor.