OsteoFlow: Lyapunov-Guided Flow Distillation for Predicting Bone Remodeling after Mandibular Reconstruction
Organizations: Department of Electrical and Computer Engineering, University of British Columbia, Vancouver, British Columbia, Canada · Sunnybrook Research Institute, University of Toronto, Toronto, Ontario, Canada · Department of Mechanical Engineering, University of British Columbia, Vancouver, British Columbia, Canada · Department of Surgery, University of British Columbia, Vancouver, British Columbia, Canada
Abstract
Predicting long-term bone remodeling after mandibular reconstruction would be of great clinical benefit, yet standard generative models struggle to maintain trajectory-level consistency and anatomical fidelity over long horizons, particularly in low-data regimes. We introduce OsteoFlow, a flow-based framework predicting Year-1 post-operative CT scans from Day-5 scans. Our core contribution is Lyapunov-guided trajectory distillation: Unlike one-step distillation, our method distills a continuous trajectory over transport time from a registration-derived stationary velocity field teacher. Combined with a resection-aware image loss, this enforces geometric correspondence without sacrificing generative capacity. Evaluated on 344 paired regions of interest, OsteoFlow significantly outperforms state-of-the-art baselines, reducing mean absolute error in the surgical resection zone by ~20%. This highlights the promise of trajectory distillation for long-term prediction in low-data clinical settings. Code is available on GitHub: https://github.com/hamidreza-aftabi/OsteoFlow.
Figures & tables
| Method | Mid slab | Full volume | ||||||
| Dice (%) | MS-SSIM (%) | MAE All HU | MAE Bone HU | Dice (%) | MS-SSIM (%) | MAE All HU | MAE Bone HU | |
| Day5 Year1 | 85.60 | 64.18 | 95.89 | 177.94 | 94.36 | 87.13 | 41.78 | 70.32 |
| ResUNet++ [ 25 ] | 92.65 | 78.10 | 56.85 | 102.67 | 96.57 | 91.82 | 27.76 | 45.80 |
| Pix2Pix [ 26 ] | 92.51 | 78.26 | 56.96 | 103.55 | 96.45 | 91.86 | 27.75 | 46.38 |
| GRIT [ 27 ] | 92.45 | 77.05 | 60.07 | 100.63 | 96.48 | 91.32 | 30.31 | 45.78 |
| cDDPM( ) [ 28 , 29 ] | 93.64 | 80.02 | 50.45 | 87.15 | 97.01 | 92.21 | 26.14 | 39.71 |
| Mid slab | Full volume | |||||||||||||
| Dice (%) | MS-SSIM (%) | MAE All HU | MAE Bone HU | Dice (%) | MS-SSIM (%) | MAE All HU | MAE Bone HU | |||||||
| Teacher Effect | ✓ | 88.97 | 66.76 | 82.47 | 139.84 | 95.25 | 87.45 | 40.07 | 63.89 | |||||
| ✓ | 90.79 | 71.23 | 77.00 | 123.85 | 95.72 | 89.05 | 37.08 | 55.91 | ||||||
| ✓ | ✓ | 91.54 | 73.20 | 69.79 | 115.66 | 96.00 | 89.68 | 34.51 | 52.91 | |||||
| ✓ | ✓ | 89.67 | 69.22 | 79.29 | 137.79 | 95.49 | 88.48 | 37.60 | 60.03 | |||||
| Full Pipeline | ✓ | ✓ | 92.86 | 79.25 | 55.24 | 97.28 | 96.70 | 92.32 | 26.98 | 45.06 | ||||
| Backbone | Mid slab | Full volume | ||||||||
| Dice (%) | MS-SSIM (%) | MAE All HU | MAE Bone HU | Dice (%) | MS-SSIM (%) | MAE All HU | MAE Bone HU | |||
| SwinUNETR | ✓ | 92.30 | 77.79 | 58.83 | 103.92 | 96.58 | 91.44 | 28.28 | 46.82 | |
| Res-UNet | ✓ | 92.86 | 79.25 | 55.24 | 97.28 | 96.70 | 92.32 | 26.98 | 45.06 | |
| SwinUNETR | ✓ | ✓ | 93.27 | 80.24 | 51.95 | 89.21 | 96.86 | 92.74 | 25.49 | 41.24 |
| Res-UNet | ✓ | ✓ | 93.65 | 81.39 | 49.47 | 83.65 | 97.00 | 93.09 | 24.95 | 40.02 |
Longitudinal Medical Image Analysis
Longitudinal medical image analysis studies how scans change over time to track disease, predict outcomes, and guide clinical decisions.
See all 100 papers