where is the well-sampled seismic data (or model), denotes the sampling operator, and denotes the observed data Chen et al. (2014a). Because of the missing traces, the operator is highly singular and thus makes equation 1 highly underdetermined.

There have been many different algorithms for solving equation 1 by imposing different constraints. In the seismic data processing field, the POCS algorithm Abma and Kabir (2006) is one of the most widely used methods to interpolate seismic data with irregularly missing traces, and has the following iterative expression:

where denotes the reconstructed data after th iteration, and are the forward and inverse sparsity-promoting transforms, is an identity matrix and is a thresholding operator. However, for regularly missing traces, POCS can not obtain satisfied results because of the strong aliasing noise in the sparse domain. Gan et al. (2015) has shown that the seislet-based POCS algorithm can still be effective for interpolating regularly missing traces provided that an accurate local slope map can be obtained. All the examples in this paper are based on the iterative framework 2. The thresholding operator is chosen as a soft-thresholding operator in the seislet domain:

where corresponds to the soft thresholding operator with a threshold value . denotes the coordinates vector. denotes the value the spatial point .

2020-02-28