EDGE-PRESERVING ENHANCEMENT OF SEISMIC IMAGES BY NONLINEAR ANISOTROPIC DIFFUSION

摘要

rocessing a seismic image comprising the steps of(a) obtaining a two or three dimensional initial image data set, wherein each element u(m=0) of the data set is the initial image intensity of a point of the image;(b) calculating for each point the partial derivatives of u(m) in n directions to obtain n derivative data sets of the partial derivatives uxi; (c) calculating for each point a symmetric nxn structural matrix S, wherein the elements spq equal uxp.uxq, and an nxn diffusion matrix D, wherein the elements dij are a function of the elements spq; (d) calculating for each point u(m+1) from u(m) using the following equation u(m+1)=u(m)+c.div(ε.Dgrad(u(m))), wherein c is a predetermined constant, 0≤c≤1, and wherein ε is a scalar, 0≤ε≤, wherein ε is close to zero when near an edge and close to 1 when far away from an edge (10, 11, 12, 14, 15); and(e) repeating steps (b) through (e) M times to obtain the processed image.