For a simple 1-D test of B-spline regularization, I chose the function shown in Figure , but sampled at irregular locations. To create two different regimes for the inverse interpolation problem, I chose 50 and 500 random locations. I interpolated these two sets of points to 500 and 50 regular samples, respectively. The first test corresponds to an under-determined situation, while the second test is clearly over-determined. Figures and show the input data for the two test after normalized binning to the selected regular bins.

bin500
50 random points binned to 500 regular
grid points. The random data are used for testing inverse
interpolation in an under-determined situation.
Figure 23 |

bin50
500 random points binned to 50 regular
grid points. The random data are used for testing inverse
interpolation in an over-determined situation.
Figure 24 |

I solved system ()-() by the iterative conjugate-gradient method, utilizing a recursive filter preconditioning for faster convergence. To construct the regularization operator , I used the method of the previous subsection with the tension-spline differential equation that I will describe in Chapter .

The least-squares differences between the true and the estimated model are plotted in Figures and . Observing the behavior of the model misfit versus the number of iterations and comparing simple linear interpolation with the third-order B-spline interpolation, we discover that

- In the under-determined case, both methods converge to the same final estimate, but B-spline inverse interpolation does it faster (with fewer iterations). However, the total computational gain is not significant because each B-spline iteration is more expensive than the corresponding linear interpolation iteration.
- In the over-determined case, both methods converge similarly at early iterations, but B-spline inverse interpolation results in a more accurate final estimate.

norm500
Model convergence in the
under-determined case. Dashed line: using linear
interpolation. Solid line: using third-order B-spline.
Figure 25 |

norm50
Model convergence in the
over-determined case. Dashed line: using linear
interpolation. Solid line: using third-order B-spline.
Figure 26 |

12/28/2000