LAUSR

Abstract In the context of faulted or occluded regions in seismic images, seismic horizon reconstruction requires the use of fast interactive approaches with respect to any bounding domain and number of passing points. The only algorithms which respect these constraints and provide a result in reasonable time in comparison with the interaction are based on the solution of a partial derivative equation either on juxtaposed quadrangular regions or using a binary mask on the whole data. While the first method requires the use of many passing points and its local nature leads to different results depending on the number and location of constraints, the second is not compatible with fast reconstruction for more than one passing point. In this paper, we propose a global and fast reconstruction method on the polygon defined by constraint points. Direct and inverse Schwarz-Christoffel space transformations lead to global reconstructions in respect to all the constraints. Experiments both on synthetic and real seismic images exhibit better results than the conventional quadrangle and mask methods. Index terms Seismic horizon reconstruction, faulted seismic image, partial derivative equation, space transformation, Schwarz-Christoffel transformation.