A method for detecting heterotopic ossification, HO, in a 3D image of a subject, the method comprising: (a) receiving, by a processor of a computing device, a 3D image of a subject; (b) applying, by the processor, a global thresholding operation to the received 3D image to produce an initial bone mask that identifies an initial region of interest within the image comprising a graphical representation of bone; (c) determining, by the processor, a boundary value map using a 3D edge detection operation applied to the initial region of interest of the image identified by the initial bone mask, wherein the boundary value map identifies and includes intensity values of voxels of the 3D image that correspond to boundaries where bone meets soft tissue; (d) determining, by the processor, a bone threshold map using the initial bone mask and the boundary value map, wherein the bone threshold map comprises, for each voxel of the initial bone mask, a threshold value determined by extrapolating values of the boundary value map to voxels within the initial bone mask; (e) determining, by the processor, a final bone mask using the bone threshold map and the received 3D image; and (f) determining, by the processor, a distance map by applying a distance transform to the final bone mask; (g) applying, by the processor, a watershed segmentation operation to the distance map to identify a set of catchment basins and/or watershed lines within the distance map; (h) generating, by the processor, a first split bone mask using the final bone mask and the identified catchment basins and/or watershed lines from step (g); (i) applying, by the processor, one or more second derivative splitting filters to the 3D image to identify a set of split line voxels within the 3D image; (j) removing, by the processor, voxels corresponding to the set of split line voxels from the first split bone mask, thereby generating the second split bone mask; (k) determining, by the processor, a plurality of labeled split binary components of the second split bone mask via one or more morphological processing operations; (l) performing, by the processor, a region growing operation within the final bone mask using the plurality of labeled spit binary components of the second split bone mask as seeds, thereby producing a labeled final bone map; and (m) rendering, by the processor, a graphical representation of the labeled final bone map. (a) receiving, by a processor of a computing device, a 3D image of a subject; (b) applying, by the processor, a global thresholding operation to the received 3D image to produce an initial bone mask that identifies an initial region of interest within the image comprising a graphical representation of bone; (c) determining, by the processor, a boundary value map using a 3D edge detection operation applied to the initial region of interest of the image identified by the initial bone mask, wherein the boundary value map identifies and includes intensity values of voxels of the 3D image that correspond to boundaries where bone meets soft tissue; (d) determining, by the processor, a bone threshold map using the initial bone mask and the boundary value map, wherein the bone threshold map comprises, for each voxel of the initial bone mask, a threshold value determined by extrapolating values of the boundary value map to voxels within the initial bone mask; (e) determining, by the processor, a final bone mask using the bone threshold map and the received 3D image; and (f) determining, by the processor, a distance map by applying a distance transform to the final bone mask; (g) applying, by the processor, a watershed segmentation operation to the distance map to identify a set of catchment basins and/or watershed lines within the distance map; (h) generating, by the processor, a first split bone mask using the final bone mask and the identified catchment basins and/or watershed lines from step (g); (i) applying, by the processor, one or more second derivative splitting filters to the 3D image to identify a set of split line voxels within the 3D image; (j) removing, by the processor, voxels corresponding to the set of split line voxels from the first split bone mask, thereby generating the second split bone mask; (k) determining, by the processor, a plurality of labeled split binary components of the second split bone mask via one or more morphological processing operations; (l) performing, by the processor, a region growing operation within the final bone mask using the plurality of labeled spit binary components of the second split bone mask as seeds, thereby producing a labeled final bone map; and (m) rendering, by the processor, a graphical representation of the labeled final bone map.