Siemens provided the MindCORE Neuroimaging Facility witha demo of their built-in motion correction software for structural scans (BioMatrix Motion Sensor – SAMER). Below are examples of T1 MPRAGE scans collected from a participant performing various commonly observed movements while being scanned. Whole brain MPRAGEs were collected with the following parameters: matrix 256 x 256; 192 sagittal slices; 0.9mm slice thickness; FOV 240 x 240mm; resolution .9 x .938 x .938; TR 2300ms; TE 2.32ms; flip angle 8°; acceleration factor 2. The scans were collected on a 3T Siemens Cima.X scanner (software version XA61) with a 32-channel head coil and with normal padding placed around the participant’s head to reduce extraneous head movements. All images are presented in native space and with no additional processing. The default sequence (sequence name: tfl3d1_228ns) produces an uncorrected and a motion-corrected version of each scan, allowing the user to evaluate the effectiveness of the motion correction algorithm.
Each of the scans below is presented with key raw scan slices as well as with motion correction enabled.
Baseline Scan
No intentional participant motion.
Heavy Breathing
Participant breathing slowly and heavily with repetitive chest and some chin motion.
Chewing
Participant moved jaw and swallowed many times during the scan.
Movement
Participant moved substantially during the scan by shifting their weight, and purposely pressing their head against the sides of the head coil and compressing the padding.
Freesurfer segmentation
Each of the above scans were processed using Freesurfer’s recon-all command with default parameters. Presented below are representative coronal slices showing the gray-matter/white-matter boundary (blue), the pial surface (red), and major subcortical structures shaded. Segmentation errors (such as can be seen in the right hippocampus) were not manually corrected.
As can be clearly seen, segmentation largely fails in cases where there are major motion events (i.e., in the “move” condition) and these failures are not overcome with the motion correction procedure. Interestingly, segmentation appears worse with motion correction in the “breathe” condition. The motion correction procedure did not produce any differences in the “baseline” and “chew” conditions between the corrected and uncorrected output (see values in table below), meaning motion artifacts were not detected in these conditions.
The following table presents volumes (in mm3) of segmented structures for each scan with and without motion correction.
StructName | baseline | baseline moco |
breathe | breathe moco |
chew | chew moco |
move | move moco |
Right-Cerebellum-Cortex | 64134.3 | 64134.3 | 63658.9 | 63027.3 | 64270.9 | 64270.9 | 65980.4 | 68368.1 |
Left-Cerebellum-Cortex | 62126.3 | 62126.3 | 61662.7 | 61369.9 | 63205.6 | 63205.6 | 72474.5 | 67019.5 |
Brain-Stem | 24158.1 | 24158.1 | 22968.3 | 23893.8 | 23464.9 | 23464.9 | 20636.4 | 24524.3 |
Left-Cerebellum-White-Matter | 16321.4 | 16321.4 | 16196.8 | 17867.4 | 16105.6 | 16105.6 | 18215.4 | 21899.8 |
Right-Cerebellum-White-Matter | 15453.9 | 15453.9 | 16046 | 15236.2 | 15443.5 | 15443.5 | 18450.4 | 20366.7 |
Left-Lateral-Ventricle | 14535.2 | 14535.2 | 14351.6 | 14625.5 | 14838.2 | 14838.2 | 13581.6 | 14297 |
Right-Lateral-Ventricle | 9986.4 | 9986.4 | 10100.2 | 10062.5 | 10027.8 | 10027.8 | 11497.4 | 9784 |
Right-Thalamus | 8781.2 | 8781.2 | 9520.7 | 9147.9 | 8925.4 | 8925.4 | 9113.8 | 8251.2 |
Left-Thalamus | 8293.9 | 8293.9 | 8796.5 | 9022.8 | 9057 | 9057 | 9255.3 | 8376.7 |
Left-Putamen | 5898.4 | 5898.4 | 4702.5 | 4709.9 | 5970.9 | 5970.9 | 4393.3 | 3790.7 |
Right-Putamen | 5793.8 | 5793.8 | 5809.3 | 5909.2 | 5804.7 | 5804.7 | 4286.4 | 5420.1 |
Left-Hippocampus | 4926.9 | 4926.9 | 4777.6 | 4709.5 | 4923.3 | 4923.3 | 2843.8 | 5428.7 |
Right-VentralDC | 4672.8 | 4672.8 | 4758.5 | 4474.5 | 4895.7 | 4895.7 | 5196.2 | 4195.8 |
Right-Caudate | 4632.5 | 4632.5 | 4509.4 | 4765.8 | 4819.9 | 4819.9 | 4397.2 | 4480.4 |
Right-Hippocampus | 4649.6 | 4649.6 | 4689.5 | 4360.2 | 4664.6 | 4664.6 | 4984.3 | 3953 |
Left-VentralDC | 4676.8 | 4676.8 | 5108.4 | 4675.7 | 4818.4 | 4818.4 | 4648.7 | 4373.9 |
Left-Caudate | 4585.1 | 4585.1 | 4514 | 4293.8 | 4570.8 | 4570.8 | 4126.9 | 4236.1 |
Left-Pallidum | 2444.6 | 2444.6 | 2273.5 | 2464.2 | 2155.1 | 2155.1 | 1990.1 | 2326.6 |
Right-Amygdala | 2257 | 2257 | 2332.4 | 2386.2 | 2281.9 | 2281.9 | 2014.5 | 2141.4 |
Right-Pallidum | 2225.3 | 2225.3 | 2347.5 | 2440.4 | 2303.7 | 2303.7 | 1891.8 | 2559.8 |
4th-Ventricle | 2297 | 2297 | 2584 | 2410.6 | 2354.6 | 2354.6 | 2746.5 | 2217 |
WM-hypointensities | 1794.9 | 1794.9 | 3031.9 | 2686.5 | 2094 | 2094 | 10610.2 | 7157.5 |
Left-Amygdala | 1993.7 | 1993.7 | 1731.9 | 1930.4 | 1924.7 | 1924.7 | 889 | 1493.5 |
Left-choroid-plexus | 1127.6 | 1127.6 | 1065.1 | 1070.1 | 1052.5 | 1052.5 | 990.5 | 1097.3 |
Right-choroid-plexus | 1064.3 | 1064.3 | 913.2 | 960.2 | 941 | 941 | 603 | 756.2 |
CSF | 1199.8 | 1199.8 | 1171.4 | 1161.9 | 1254.4 | 1254.4 | 1352.5 | 1330.8 |
CC_Anterior | 995.6 | 995.6 | 1123.1 | 1108.3 | 1154.8 | 1154.8 | 1156.8 | 1006.6 |
CC_Posterior | 995.3 | 995.3 | 978.4 | 976.3 | 986.2 | 986.2 | 1124.3 | 1013.6 |
3rd-Ventricle | 1116 | 1116 | 1162.6 | 1147 | 1115.9 | 1115.9 | 1461.8 | 1394.8 |
Left-Accumbens-area | 744.7 | 744.7 | 792.4 | 666.7 | 713.9 | 713.9 | 656.9 | 557.6 |
Right-Accumbens-area | 694.6 | 694.6 | 745.6 | 728 | 749 | 749 | 797.1 | 744.3 |
Left-Inf-Lat-Vent | 739.2 | 739.2 | 790.2 | 1136.6 | 661.1 | 661.1 | 275.5 | 650 |
CC_Mid_Posterior | 579.3 | 579.3 | 579 | 552.9 | 580.5 | 580.5 | 776.6 | 617.3 |
CC_Mid_Anterior | 566.8 | 566.8 | 589.3 | 568.5 | 546.9 | 546.9 | 597.4 | 527.9 |
CC_Central | 495 | 495 | 481.2 | 489.4 | 495.2 | 495.2 | 543.8 | 525.4 |
Right-Inf-Lat-Vent | 423.8 | 423.8 | 575.2 | 413.6 | 464.9 | 464.9 | 744.8 | 408.3 |
Optic-Chiasm | 282 | 282 | 304.6 | 301.5 | 298.3 | 298.3 | 253.4 | 325.1 |
Left-vessel | 40.5 | 40.5 | 8.8 | 22.6 | 54.7 | 54.7 | 3.4 | 34.7 |
Right-vessel | 18.4 | 18.4 | 51.8 | 78.6 | 36 | 36 | 0 | 27.9 |
Recommendation
The Siemens BioMatrix Motion Sensor procedure appears to produce qualitatively better output when there are moderate amounts of movement, as in the “breathe” condition above. These apparent improvements may not be beneficial in automated segmentation pipelines, however. In cases where there is little motion (“baseline” and “chew” conditions), the motion correction algorithm does not appear to be engaged. The motion correction algorithm is not able to overcome artifacts in cases where there are large motion events (“move” condition). If such large motion artifacts are apparent on reconstructed scans, it is recommended to repeat the scan, perhaps with additional coaching or other support for the participant.
Acknowledgements: Thank you to Alex Gordienko for developing and implementing the behavioral protocol here, as well as for the reconstruction and display of the scan quality.