Hyperostosis in Combination With Low Skull Density Ratio: A Potential Contraindication for Magnetic Resonance Imaging-Guided Focused Ultrasound Thalamotomy

A rapid element pressure field simulation method for transcranial phase correction in focused ultrasound therapy

Transcranial focused ultrasound ablation has emerged as a promising technique for treating neurological disorders. The clinical system exclusively employed the ray tracing method to compute phase aberrations induced by the human skull, taking into account computational time constraints. However, this method compromises slightly on accuracy compared to simulation-based methods. This study evaluates a fast simulation method that simulates the time-harmonic pressure field within the region of interest for effective phase correction. Experimental validation was carried out using a 512-element, 670 kHz hemispherical transducer for fourex vivoskulls. The ray tracing method achieved a restoration ratio of 64.81% ± 4.33% of acoustic intensity normalized to hydrophone measurements. In comparison, the rapid simulation method demonstrated improved results with a restoration ratio of 73.10% ± 7.46%, albeit slightly lower than the full-wave simulation which achieved a restoration ratio of 75.87% ± 5.40%. The rapid simulation methods exhibited computational times that were less than five minutes for parallel computation with 8 threads. The incident angle was calculated, and a maximum difference of 6.8 degrees was found when the fixed position of the skull was changed. Meanwhile, the restoration ratio of acoustic intensity was validated to be above 70% for different target positions away from the geometrical focus of the transducer. The favorable balance between time consumption and correction accuracy makes this method valuable for clinical treatment applications.

Comparison of the impact of skull density ratio with alternative skull metrics on magnetic resonance-guided focused ultrasound thalamotomy for tremor

OBJECTIVE

One of the key metrics that is used to predict the likelihood of success of MR-guided focused ultrasound (MRgFUS) thalamotomy is the overall calvarial skull density ratio (SDR). However, this measure does not fully predict the sonication parameters that would be required or the technical success rates. The authors aimed to assess other skull characteristics that may also contribute to technical success.

METHODS

The authors retrospectively studied consecutive patients with essential tremor who were treated by MRgFUS at their center between 2017 and 2021. They evaluated the correlation between the different treatment parameters, particularly maximum power and energy delivered, with a range of patients' skull metrics and demographics. Machine learning algorithms were applied to investigate whether sonication parameters could be predicted from skull density metrics alone and whether including combined local transducer SDRs with overall calvarial SDR would increase model accuracy.

RESULTS

A total of 62 patients were included in the study. The mean age was 77.1 (SD 9.2) years, and 78% of treatments (49/63) were performed in males. The mean SDR was 0.51 (SD 0.10). Among the evaluated metrics, SDR had the highest correlation with the maximum power used in treatment (ρ = -0.626, p < 0.001; proportion of local SDR values ≤ 0.8 group also had ρ = +0.626, p < 0.001) and maximum energy delivered (ρ = -0.680, p < 0.001). Machine learning algorithms achieved a moderate ability to predict maximum power and energy required from the local and overall SDRs (accuracy of approximately 80% for maximum power and approximately 55% for maximum energy), and high ability to predict average maximum temperature reached from the local and overall SDRs (approximately 95% accuracy).

CONCLUSIONS

The authors compared a number of skull metrics against SDR and showed that SDR was one of the best indicators of treatment parameters when used alone. In addition, a number of other machine learning algorithms are proposed that may be explored to improve its accuracy when additional data are obtained. Additional metrics related to eventual sonication parameters should also be identified and explored.