Presented herein are systems and methods that facilitate automated segmentation of 3D images of subjects to distinguish between regions of heterotopic ossification (HO) normal skeleton, and soft tissue. In certain embodiments, the methods identify discrete, differentiable regions of a 3D image of subject (e.g., a CT or microCT image) that may then be either manually or automatically classified as either HO or normal skeleton.