The purpose of this work is to propose a method to characterize tumour heterogeneity on MRI, using probabilistic classification based on a reference tissue. The method uses maps of the apparent diffusion coefficient (ADC), T relaxation, and a calculated map representing high-b-value diffusion-weighted MRI (denoted simDWI) to identify up to five habitats (i.e. sub-regions) of tumours. In this classification method, the parameter values (ADC, T , and simDWI) from each tumour voxel are compared against the corresponding parameter probability distributions in a reference tissue. The probability that a tumour voxel belongs to a specific habitat is the joint probability for all parameters. The classification can be visualized using a custom colour scheme. The proposed method was applied to data from seven patients with biopsy-confirmed soft tissue sarcoma, at three time-points over the course of pre-operative radiotherapy. Fast-spin-echo images with two different echo times and diffusion MRI with three b-values were obtained and used as inputs to the method. Imaging findings were compared with pathology reports from pre-radiotherapy biopsy and post-surgical resection. Regions of hypercellularity, high-T proteinaceous fluid, necrosis, collagenous stroma, and fibrosis were identified within soft tissue sarcoma. The classifications were qualitatively consistent with pathological observations. The percentage of necrosis on imaging correlated strongly with necrosis estimated from FDG-PET before radiotherapy (R = 0.97) and after radiotherapy (R = 0.96). The probabilistic classification method identifies realistic habitats and reflects the complex microenvironment of tumours, as demonstrated in soft tissue sarcoma.