To preserve the edges of the model Clapp et al. (1998) propose adding a weight to zero the model residual at the edges. Lizarralde and Swift (1999) implement a similar approach for the inversion of VSP data for interval velocity. This approach requires human intervention for reflector picking. We want to design a weight which de-weights edges in the model residual, but which is estimated automatically.

The 2-D gradient magnitude is a good isotropic edge-detection operator that can be used to calculate the diagonal weights. As we show in the deblurring problem Valenciano et al. (2003), using the gradient magnitude we can iteratively obtain sharp edges.

We perform the following non linear iterations: starting with , at the *k*^{th} iteration the algorithm solves

(8) |

(9) |

10/14/2003