Figure 3

Figure 4

Since these data are collected in two series of one-dimensional tracks, it would be easiest to estimate a pair of one-dimensional PEFs on these two sets of tracks, as shown in the top halfs of Figures and .

We now have two PEFs which have been estimated in a data space, but the model which we wish to regularize with these PEFs is in a different space. This requires the introduction of two additional linear interpolation operators and , which pull bins from the model space into the ascending and descending track data spaces, respectively. The mappings used for these operators are shown in Figure .

Figure 5

Now that we have both two prediction-error filters for regularizations operators as well as linear interpolation operators that pull model points into the data space, we can put everthing together in the following fitting goals,

(3) |

The 1D PEFs can also be replaced by 2D PEFs that are estimated by scaling the filter so that it covers multiple sparse tracks. If this approach is taken, the interpolation can occur in the data space where the PEFs are estimated (using a single PEF for each of the two track spaces), or in the model space (using both PEFs simultaneously). The more straightforward data-space interpolation is shown in Figure .

Figure 6

Figures a and b are simply the interpolation of the data space with 2D PEFs estimated on the sparse tracks. Figures c and d are those interpolated results mapped back to the model space by using fitting goals (1), where the input data are now the interpolated sets of tracks in the first two panels. Since fitting goals (1) were applied when generating the new tracks, the track derivative is no longer necessary.

The results are mixed, as Figure c shows that the trend of the ridge was correctly identified by the PEF estimated on the ascending tracks. The PEF estimated on the descending tracks did not fare so well, as the direction of the ridge in the interpolated tracks of Figure d does not match the densely sampled tracks in Figure d. This is because the descending tracks are oblique to the structure, so the structure is aliased beyond the point where a spaced PEF can interpolate accurately. In either case, the result is better than that obtained with a Laplacian, and in the case of the ascending tracks is not that far from the PEF estiamted on a fully-sampled model space shown in Figure c.

5/3/2005