I use 3D slicer5.0.3 and in ANTs there are three different methods(stages) which are QuickSyN, Rigid, Rigid plus Affine.
Now I have many CT images of different subjects head. And I want to registrate the different CT images to the same reference before I extract the landmark location information.
I compared these three methods that for Rigid and Rigid+Affine, the new volume after registration looks similar. But for QuickSyN, it is different. And I am not sure whether the QuickSyN will change the shape and the CT value of the original CT. Just looks very different and the max CT value is also changed.
I would like to ask and also double check whether the method of QuickSyN will change the shape of CT images. Because I would like to keep the shape and if it is true, I probably will choose from Rigid and Rigid+Affine.
I guess you are looking at the SlicerANTs extension, right? QuickSyN appears to be a SlicerANTs setting that carries out multiple steps of ANTs registration: Rigid, then Affine, then SyN. Learn about the meanings of these here.
And I am not sure whether the QuickSyN will change the shape and the CT value of the original CT.
Registration will not change CT value; it doesn’t know about your imaging modality and will not touch the voxel values.
The question about changing “shape” depends on what you mean by shape, i.e. what you’re interested in preserving. SyN is a type of deformable (“non-parametric”) registration using velocity fields, so probably doesn’t preserve shape by any geometric definition of the word, but it should preserve topology. Affine will not necessarily preserve distances, but the transform will be uniform across the image (looks the same in each local region; affine just means linear with a possible translation) and it should preserve angles (typically an affine transform need not preserve angles, but the way ANTs seems to use the word it should). Rigid preserves both angles and distances.
To be clear, that new volume is your original volume after the resultant deformation field calculated by ANTs applied to it, regardless whether the registration is Rigid+affine only or SyN. The intensity values are likely to be different from the original CT, due to necessary interpolation to warp the image in all cases. Just because max value remained the same doesn’t mean that the individual voxel intensities are not modified. The effect of the affine registration might be low, if the reference and target is already in the same orientation, and of approximate size.
If your targets and reference are similar in shape, affine might be sufficient. For more precise local deformation, you need to use SyN.
Hi, thank you a lot! Now I am clear that the CT value will change.
The targets and reference are similar in shape because all of them are head but they are different in size.
Yesterday I employed the Rigid method on all targets because I dont want to change the size. But for some cases, the result is not good. So I used ‘IGT Fiducial registration Wizard’ to rigidly registrate my targets. For example, I used three anotamy points (which represents the same point for different individual) to build my rigid transfor matrix.
But I am confused that for those cases with different sizes, how did the target position itself? And also whether the 3D slicer rendering model can represent the actual size in the real world?
I am still not clear on what you are trying to do with the image registration. There are typically two ways you can use the registration.
There are a set of canonical landmarks in some sort of reference sample, and you deform your new target to this reference, and then use the inverse of the transformation (using a metric like SyN that does provide invertible transforms) to transfer the canonical landmarks back to the new subject.
You can register a number of new samples into an existing reference sample, and then do the statistics on the resultant warp fields (e.g., compare whether warp fields are statistically significant across groups, or visualize the expansion/contraction of individual voxels etc).
Whether you will use affine or deformable registration (or the total warp field, which is combining these two) depends on the questions you are asking.
Newborn to 9 months old is a very difficult registration problem due to closure of the fontanelles, matching the size of the skull will not be sufficient, as the shape of the skull is changing allometrically (in an non-uniform) way. The choice of your registration parameters (similarity metric and the geometric transformation model) will affect the output, and you have no real basis to prefer one over the other, if you do not have constrains or an objective measure to assess the outcome (e.g., which method gives you the landmark estimates closest to your expert landmarks).
Below, the second row are outputs of antsRegistration using the same pair of images in the top row (new born to 6 months of age) using different similarity and transformation models. A uses a large deformation model based on mean squared metric and velocity fields, B is the SyNCC (SyN with cross-correlation). Both of them does a good job of reducing the global size difference due to growth, but the local deformations are very different. Without having an objective method of choosing one over the other, you are sort of stuck.
Landmarks are usually the way to constrain this (or use them to choose one outcome over the other), but sounds like you are trying to use the image registration to transfer the landmarks perhaps?