EDGE-PRESERVING ENHANCEMENT OF SEISMIC IMAGES BY NONLINEAR ANISOTROPIC DIFFUSION

摘要

rocessing a seismic image comprising the steps of(a) obtaining a two or threedimensional initial image data set, wherein each element u(m=0) of the dataset is the initial image intensity of a point of the image;(b) calculating foreach point the partial derivatives of u(m) in n directions to obtain nderivative data sets of the partial derivatives uxi; (c) calculating for eachpoint a symmetric nxn structural matrix S, wherein the elements spq equaluxp.uxq, and an nxn diffusion matrix D, wherein the elements dij are afunction 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(.epsilon..Dgrad(u(m))), whereinc is a predetermined constant, 0=c=1, and wherein .epsilon. is a scalar,0=.epsilon.=, wherein .epsilon. is close to zero when near an edge and closeto 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.