Hello, I am trying to apply a deformation field to a volume using the BRAINS resample image module but it ends up not doing anything and I am not sure why. Some context is needed:
We are looking at the AAPM TG-132 report and assessing the deformation algorithms at our center compared to gold standard data. As part of this, we are provided a volume, the deformed volume, and the deformation field that was used. I have attached the original dataset. The other dataset can be downloaded here: https://www.dropbox.com/s/jnyydbcie9bg9fu/TG132P2-1-1%20-%20Dataset%201%20deformed.zip?dl=0
The dropbox download contains the deformed image as well as 2 registration dicoms. It is unclear to me why there are two but I tried with both.
In any case my process was as follows:
- Import the Original dicom, the deformed dicom, and the 2 registrations
- Use the BRAINS resample image to apply the transform with the original image as input, a new volume as output, and either of the registrations as the transform
- Compare the new image and the original image.
Output: They are the same volume with no deformations.
I am wondering if I am running into a memory issue as I am running this on a fairly limited machine. I just didn’t expect it to return a volume in which nothing occurred and no error messages were specified. I also attempted the following:
in the python module in slicer I ran the following commands:
voxelarray = array(‘SpatialRe’)
which returned None (i.e., the spatial registration is empty?).
I also imported the registration in Matlab and looked at the vector grid header which contains the transform and it was not empty as far as I could tell.
Any help would be great, thanks