Non Convex Operators for Electromagnetic Geosounding Noise

Abstract

The problem of model recovering in the presence of impulse noise on the data is considered for the magnetotelluric inverse problem. The application of total variation regularization along with L1-norm penalized data fitting (TVL1) is the usual approach for the impulse noise treatment in image recovery. This combination works poorly when a high level of impulse noise is present on the data. A nonconvex operator named smoothly clipped absolute deviation (TVSCA) was recently applied to the image recovery problem. This operator is solved using an sequence of TVL1 equivalent problems, providing a significant improvement over TVL1. In practice, TVSCAD requires the selection of several parameters, a task that can be very difficult to attain. A more simple approach to the presence of impulse noise in data is presented here. A non convex function is also considered in the data fitness operator, along with the total variation regularization operator. The nonconvex operator is solved by following a half-quadratic procedure of minimization. Results are presented for synthetic and also for field data, assessing the proposed algorithm's capacity in model recovering under the influence of impulse noise on data for the magnetotelluric problem.