Abstract: A seismic data set is processed by selecting a set of gathers in the seismic data set. A frequency independent constraint matrix is calculated in a first gather. A model solution in at least one second gather is calculated by using the constraint matrix in a high-resolution Radon transform.