Calcul rapide de la dose : un modèle grossier vaut bien un modèle précis ✛

Original: "Fast Patient-Specific Breast CT Dosimetry: 22-Fold Acceleration of Monte Carlo MGD Estimation"
arXiv:2609.03263v1 · 2026-09-03 · CC BY 4.0 · 1 min
Le calcul de la dose d'irradiation pour le sein est passé de 27 heures à 72 minutes, avec une précision quasi inchangée.

Quand vous vous rendez dans un quartier inconnu, vous n'avez pas besoin d'une carte où chaque maison est dessinée. Les rues et les intersections suffisent. Il en va de même pour le calcul de la dose d'irradiation lors d'une tomographie mammaire : chaque détail n'est pas nécessaire, seule compte la vue d'ensemble.

Auparavant, pour chaque patiente, on construisait un modèle du sein composé de millions de minuscules cubes (voxels). L'ordinateur les parcourait tous, et le calcul prenait environ 27 heures. C'est comme si, avant un voyage, vous étudiiez le plan de chaque appartement de chaque immeuble.

Les chercheurs ont testé ce qui se passe lorsque l'on agrandit les cubes.

Il s'est avéré que la taille du cube dans le plan de l'image influence davantage la dose que son épaisseur. On a donc augmenté davantage l'épaisseur, et moins la largeur et la hauteur. Résultat : le nombre de cubes a été réduit d'un facteur de 20 à 31, tandis que la dose n'a varié que de quelques pour cent. Ensuite, on a également réduit le nombre de particules simulées. Le calcul est passé de 27 heures à 72 minutes, avec une précision quasi inchangée.

Cette méthode sera utile dans les hôpitaux : la dose pourra être estimée rapidement, avant même le scan principal. Une carte grossière s'est avérée aussi bonne qu'une carte détaillée, mais 22 fois plus rapide.

Le fantôme initial d'un sein contenait jusqu'à 95 millions de voxels, et après optimisation, il n'en restait qu'environ 3 millions. La différence dans la dose calculée était inférieure à un demi pour cent du niveau d'exposition cible.

MGD = DgN \times K_{air}
MGD est la dose glandulaire moyenne, DgN est le coefficient de dose normalisé, K_air est la kerma dans l'air
MGD = \frac{\sum E_i}{\sum m_i}
E_i est l'énergie absorbée dans le voxel glandulaire i, m_i est la masse de ce voxel
\sigma \propto \frac{1}{\sqrt{N}}
σ est l'incertitude statistique, N est le nombre d'histoires photoniques primaires
Notions
segmentation tumorale fantôme tomodensitométrie