How to use ITK to segment bones in osteoporotic patients

Osteoporosis reduces bone density and thus the strength of bone in the CT data set (shown in the figure), making it impossible to use traditional thresholding techniques to handle segmentation.

I noticed that the Morphon registration algorithm mentioned in this paper
seems to work for this case. There is a corresponding ITK implementation, but the source code it provides has a version mismatch with the current ITK.

I wonder if there is any other possible way in ITK to deal with this situation? Is DisplacementFieldTransform a viable approach? Displacement fields are also mentioned in Morphon, but I didn’t find an example of applying DisplacementFieldTransform, so I turn to you for help.

You could use elastix, possibly via its ITKElastix interface.

Thank you! I will try.