Rigid registration for volumetric meshes

Hello,
In my recent project, I’m working with volumetric meshes of scapulae. I have already rigidly aligned them, as you can see in the picture, but I now aim to align a template non-rigidly to each of my scapulae. As I work with volumetric meshes, I would like to find not only a way to deform the surface of my template on the other meshes, but also the points inside.

The idea is to deform the template on the other meshes in order to get a set of deformations of the template, which are in correspondence (the points are organised the same way). Thus, an anatomical point on one mesh would have the same ID on all other meshes.

The final goal would be to create a statistical shape & intensity model that could represent the variety of shapes and intensities of the bone.

From now on, I’ve only been able to perform non-rigid registration of the original CT scans using ITK, which doesn’t work to create a statistical shape and intensity model, as it needs to be performed on meshes.

Is there anyone who knows how I could deal with that?
Thank you in advance for your answers!
Nicolas.