A Rough Map of the Breast: Dose Calculation 22 Times Faster ✛

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 · 2 min
Monte Carlo dose estimates for a specific patient sped up from 27 hours to 72 minutes without losing accuracy.

A mapmaker drawing a route through a city doesn't sketch every brick in every wall. Streets, intersections, and distances are enough. The same goes for calculating radiation dose in breast tomography: you don't need every cell, just a faithful overall picture.

Previously, for each patient, a digital model of the breast was built from millions of tiny cubes (voxels), and the computer would churn through them for hours. The calculation took about 27 hours. That's like planning a trip across town by first studying the floor plan of every apartment in every building along the way.

The authors tested how much they could coarsen the map while keeping the route accurate.

It turned out that the cube size in the imaging plane affects the dose more than its thickness. Making the cubes 2.5–3 times larger in width and height changes the dose by 5–9%. But increasing only the slice thickness by the same 2.5–3 times shifts the dose by just 1–3%. So they reduced the resolution unevenly: more along the thickness, less in the plane. A final voxel size of 0.30 × 0.30 × 3.00 mm³ kept dose accuracy within 3.6–4.4% for all four samples. The number of cubes dropped by a factor of 20–31.

Next, they cut the number of simulated photons from 3 billion to 200 million. The statistical uncertainty stayed below 0.3%.

The total computation time shrank from 1600 minutes to 72 minutes. That's like a mapmaker swapping a detailed atlas for a tourist map and still not getting lost.

Previously, it was unclear whether such a coarse breast model could be used for synchrotron phase-contrast CT. Now it's clear: yes. The method will be useful in clinics where dose needs to be estimated quickly, before the main scan. The rough map turned out to be just as good as the detailed one, but 22 times faster.

The original phantom of one breast contained up to 95 million voxels, and after optimization only about 3 million remained. Yet the difference in the calculated dose was less than half a percent of the target radiation level.

MGD = DgN \times K_{air}
MGD is mean glandular dose, DgN is normalized dose coefficient, K_air is air kerma
MGD = \frac{\sum E_i}{\sum m_i}
E_i is the energy absorbed in glandular voxel i, m_i is the mass of that voxel
\sigma \propto \frac{1}{\sqrt{N}}
σ is statistical uncertainty, N is the number of primary photon histories
Concepts
tumor segmentation phantom computed tomography medical image segmentation