problematic case ABBE-R
seems the BiLaplacian split found a wrong threshold:

ABBE-R.  BiLaplacian split.  Residual = 2549

However, when we re-segment the BiGaussian model, we obtain much lower residual.

BiGaussian fit.  Residual = 645

ABBE-R:   Split using BiGaussian

However, for Phantoms the BiLaplacian was much better than Gaussian at least for most cases and all phantoms were processed using exclusively BiLaplacian.


Proposal:  We evaluate histogram fit using both {BiLaplace,BiGauss}.  We choose the method that provides smaller fitting residual.