** Next:** Tests on synthetic data:
** Up:** Guitton: IRLS vs. Huber
** Previous:** The Iteratively Reweighted Least

We now compare IRLS and the Huber solver doing velocity-stack inversion.
This method was first introduced by Thorson and Claerbout (1985)
and completed with an IRLS algorithm by Nichols (1994).
Guitton and Symes (1999) introduced the Huber misfit function with a specific non-linear
solver as an alternative to IRLS.
The velocity estimation is made using inverse methods to minimize the function
where *E* is an error measure function, **H** is the modeling operator that
transforms the model space (velocity domain) into the data space
(CMP gathers). The adjoint operator (**H^**) is the well-known hyperbolic Radon transform
operator. The forward operation is
| |
(7) |

and the adjoint transformation becomes
| |
(8) |

where *w*_{o} is a weighting function Claerbout and Black (1997).
We now solve the inverse problem using both the Huber solver and IRLS on synthetic and real data.
Note that in the following sections, I assume that when
, we nearly solve the *l*^{1} problem for both the Huber solver and IRLS.

** Next:** Tests on synthetic data:
** Up:** Guitton: IRLS vs. Huber
** Previous:** The Iteratively Reweighted Least
Stanford Exploration Project

4/27/2000