Fitting Intravoxel Incoherent Motion Model to Diffusion MR Signals of the Human Breast Tissue using Particle Swarm Optimization

Intravoxel incoherent motion (IVIM) modeling offers the parameters f, D and D* as biomarkers for different lesion types and cancer stages from diffusion MR signals. Challenges with the available optimization algorithms in fitting the model to the signals motive new studies for improved parameter estimations. In this study, one thousand value sets of f, D, D* for human breast are assembled and used to generate five thousand diffusion MR signals considering noise-free and noisy situations exhibiting signal-to-noise ratios (SNR) of 20, 40, 60 and 80. The estimates of f, D, D* are obtained using Levenberg-Marquardt (LM), trust-region (TR) and particle swarm (PS) algorithms. On average, the algorithms provide the highest fitting performance for the noise-free signals (R2adj=1.000) and great fitting performances on the noisy signals with SNR>20 (R2adj>0.988). TR algorithm performs slightly better for SNR=20 (R2adj=0.947). TR and PS algorithms achieve the highest parameter estimation performance for all the parameters while LM algorithm reveals the highest performance for f and D only on the noise-free signals (r=1.00). For the noisy signals, performances increase while SNR increases. All algorithms accomplish poor performances for D* (r=0.01-0.20) while TR and PS algorithms perform same for f (r=0.48-0.97) and D (r=0.85-0.99) but remarkably better than LM algorithm for f (r=0.08-0.97) and D (r=0.53-0.99). Overall, TR and PS algorithms demonstrate better but indistinguishable performances. Without requiring any user-given initial value, PS algorithm may facilitate improved estimation of IVIM parameters of the human breast tissue. Further studies are needed to determine its benefit in clinical practice.