Home Geology and Mineralogy Seismic data filtering using non-local means algorithm based on structure tensor
Article Open Access

Seismic data filtering using non-local means algorithm based on structure tensor

  • Shuai Yang , Anqing Chen EMAIL logo and Hongde Chen
Published/Copyright: May 25, 2017
Download Article (PDF)
Become an author with De Gruyter Brill
Submit Manuscript Author Information
Open Geosciences
From the journal Open Geosciences Volume 9 Issue 1

Abstract

Non-Local means algorithm is a new and effective filtering method. It calculates weights of all similar neighborhoods’ center points relative to filtering point within searching range by Gaussian weighted Euclidean distance between neighborhoods, then gets filtering point’s value by weighted average to complete the filtering operation. In this paper, geometric distance of neighborhood’s center point is taken into account in the distance measure calculation, making the non-local means algorithm more reasonable. Furthermore, in order to better protect the geometry structure information of seismic data, we introduce structure tensor that can depict the local geometrical features of seismic data. The coherence measure, which reflects image local contrast, is extracted from the structure tensor, is integrated into the non-local means algorithm to participate in the weight calculation, the control factor of geometry structure similarity is added to form a non-local means filtering algorithm based on structure tensor. The experimental results prove that the algorithm can effectively restrain noise, with strong anti-noise and amplitude preservation effect, improving PSNR and protecting structure information of seismic image. The method has been successfully applied in seismic data processing, indicating that it is a new and effective technique to conduct the structure-preserved filtering of seismic data.

Keywords: Non-Local means; Structure tensor; Filter; Similarity; Interpretative processing

1 Introduction

Filtering is a key step in seismic data processing [1]. Due to increasing demands for high-quality data, interpretation-oriented processing becomes more and more important. Improved signal-to-noise ratio (S/N) by filtering processing may facilitate subsequent studies of seismic geomorphology, seismic facies and micro-structures as well as reservoir prediction.

Seismic data filtering may be realized through Gaussian filtering, median filtering, mean filtering and various transform-domain approaches, e.g. Fourier transform-based f-k and f-x filtering, wavelet transform-based filtering [2,3], curvelet transform-based filtering [4] and Radon transform-based filtering. Additionally, there is the filtering technique based on a geostatistical method which makes use of the variogram function after decomposition to do geostatistical estimation respectively to get the estimated value of the various components of the original data, including the effective information and noise values. It has been applied to filtering of seismic data [5]. These techniques are effective for seismic data filtering, but most of them have the disadvantages of blurring the details and edges of an image in the process of noise reduction. Marginal reflections or seismic discontinuities may be caused by faults, fractures, channels, lenticular geobodies or reefs. Such features are closely related to the interpretation of hydrocarbon reservoirs. Consequently some new techniques, e.g. filtering based on mathematical morphology [6], structure constrained edge-preserved filtering [7], and anisotropic diffusion filtering based on partial differential equation [8], have been developed to preserve edge structures or the contacts between reflections and formations.

Buades et al. made a comparative study of some typical filtering techniques and developed a non-local means (NLM) algorithm for noise removal [9-11]. This algorithm was demonstrated to have better performance than other typical denoising methods. David Bonar et al. successfully applied the NLM algorithm in seismic data denoising with the idea of neighborhood matching and correlation as well as weighted smoothing [12, 13], but they did not take account of seismic geometric properties such as directivity and local contrast. Apparently, if a seismic profile or attribute slice is considered as a digital image, each pixel would be characterized by two properties, i.e. numerical property of the pixel itself and geometric structure of neighboring pixels. In this paper, we develop a new filtering algorithm for seismic data (image) processing. It integrates similarity evaluation by local contrast with NLM filtering based on structure tensor, which captures the geometric structure of an image [8]. The calculation of neighborhood similarity distance involves numerical similarity and neighborhood geometric distance. A coherence attribute, designed to describe local contrast, is added into weighting factor estimation to evaluate similarity. A control factor is also included to adjust the weight of geometric property. As per model data processing and residual error images, this algorithm achieves a better result of edge preservation. In addition, a large control factor reduces the weight of geometric property for better noise attenuation. The results of field data processing also show that this algorithm is capable of edge preservation while noise suppression.

2 Theoretical studies

2.1 NLM filtering

Previous denoising algorithms were developed under the assumption that the raw image is regular. Moreover many details of the image may also be smoothed in denoising. These two disadvantages have been remedied by NLM filtering developed by Buades et al. [14]. The basic idea of this algorithm is that an image usually contains a mass of noises, thus an image ν is defined as

v=u+n(1)

As per Equation (1), the image v is composed of the original noise-free image u and random noise n. At the pixel i, the NLM filtering result ν̂(i) is simply the weighted average of all of pixels within the noisy image v.

v^(i)=∑jw(i,j)v(j)(2)

In Equation (2), the weighting factor w(i, j) is dependent on the similarity between pixels i and j, in addition must satisfy the conditions 0 ≤ w(i, j) ≤ 1 and ∑jw(i,j)=1. The similarity between the pixels i and j is represented by Gaussian weighted Euclidean distance D(i, j) from the pixel i and its neighborhood Ni to the pixel j and its neighborhood Nj. Note that each pixel i has its own independent weighting factor of the other pixels j within the image, which is calculated by Equation (3).

w(i,j)=1Z(i)exp−D2(i,j)h2(3)

where Z(i) is the normalizing factor to ensure ∑jw(i,j)=1,and defined by

Z(i)=∑jexp−D2(i,j)h2(4)

The filtering parameter h is a constant which controls the decay rate of the exponential function and thus determines the degree of filtering. For example, a large value for h will provide very similar weight for all pixels j, thus the image would be blurred. On the contrary, a small value for h will provide a significant weight for only a few of the pixels j, thus noises cannot be attenuated sufficiently.

The Gaussian weighted Euclidean distance D2(i, j) in Equation (3) is defined by the following expression.

D2(i,j)=|v(Ni)−v(Nj)|2,a2=∑lnlGa(xl,yl)(v(Ni(l))−v(Nj(l)))2(5)

where the operator |∙|2,a2 denotes the squared factor of Gaussian weighted Euclidean distance D2(i, j), Ni represents a neighborhood with the center at the pixel i, which is usually a square domain, Ga represents the Gaussian kernel with the standard deviation a, and l represents one of the total nl elements within a neighborhood.

For a 2D image, the Gaussian kernel can be defined by,

Ga(x,y)=exp−(x−x0)2+(y−y0)22a2(6)

where X0 and y0 represent the center of Gaussian kernel with x and y corresponding to the coordinates of the element l in equation (5).

For three neighborhoods (in red, green and blue, respectively) around a square image (in Figure 1a), the weights of the centers in respect to other all of elements within the image are calculated by Equation (3). We obviously can see, for red neighborhood (in Figure 1a), the centers of all red neighborhoods in Figure 1b have larger weighting factors due to good similarity. Similarly for green neighborhood (in Figure 1a), the centers of all green neighborhoods in Figure 1c have larger weighting factors. For blue neighborhood (in Figure 1a), the centers of all blue neighborhoods in Figure 1d have larger weighting factors.

Figure 1: The weights schematic diagrams of Non-Local means algorithm
Figure 1:

The weights schematic diagrams of Non-Local means algorithm

As shown in the Figure 1 and Equation (5), for two neighborhoods with identical numerical distribution and structures, the square factor of Gaussian weighted Euclidean distance between them is equal to zero. Therefore, for the center of red neighborhood (in Figure 1a), the centers of all red neighborhoods in Figures 1b have identical weighting factors, which is obviously unreasonable. Hence, D2(i, j) in Equation (3) is corrected to be D¯2 (i, j) by the idea of inverse distance weighted interpolation.

D¯2(i,j)=D2(i,j)+|i−j|2(7)

where |i−j| is the geometric distance between the pixels i and j, for the pixel i with the coordinates of xi and yi and the pixel j with the coordinates of xj and yj, |i—j|2=(xi−xj )2 + ( yi −yj )2.

In addition, there is a defect in neighborhood distance measurement in the Equation (5). The center of the neighborhood has a much larger weighting factor than other all of elements within the region of search, which means if i=j, the contribution from adjacent elements to similarity is reduced by over-weighting. For example, if the center X is a noisy pixel and surrounding pixels are less noisy, the distance calculated by Equation (5) with the large weighting factor for X is unfavorable for denoising. In order to solve this problem, we adopt that 0≤w(i, i)< 1 when i=j. Then, the impact of the weighting factor at the center may be somewhat mitigated. In numerical calculation, w(i, i)=0.5 or w(i, i) = max(w(i, j)∀i≠j). In following calculation, w(i, i) = 0.5.

Each pixel in the image would be compared with all pixels in the process of NLM filtering, which results in mass computation and low efficiency. For an image with M×N pixels, M×N weighting factors would be calculated repeatedly for each pixel. Therefore altogether (M×N)2 factors would be calculated in the filtering. It is difficult to popularize such an inefficient algorithm. More importantly, many irrelevant neighborhoods with large distance are assigned weighting factors, this would impair the result of filtering.

Buades [14] proposed that the region of search is limited within a square area, i.e. a neighborhood, around the pixel to be handled so as to improve efficiency of filtering. Thus for a region of S×S and an image with M×N pixels, the computational complexity is M×N×s(S×S-1) (because it is unnecessary to calculate the weighting factor for the center of the neighborhood itself). As a result, the computational complexity is greatly reduced compared with the original NLM algorithm.

2.2 Structure tensor and local contrast

The definition of structure tensor proposed by J. Weickert et al. [8] is very useful in image analysis to estimate the direction field and extract local structures of the image.

Sρ=Gρ∗(∇vσ⋅∇vσT)=∂vσ∂x2∗Gρ∂vσ∂x∂vσ∂y∗Gρ∂vσ∂x∂vσ∂y∗Gρ∂vσ∂y2∗Gρ(8)

where Sρ is the structure tensor, ∇v is the gradient of the image, vσ = v * Gσ is the image of v after Gaussian filtering, Gσ is the Gaussian kernel defined by σ. A procedure of Gaussian filtering before gradient estimation could eliminate the impacts of noises, ρ is the variance of the Gaussian filter Gρ. More information involved by the convolution of Gρ may facilitate the delineation of such geometric structures as linear laminae, fracture boundaries and angular regions.

The structure tensor is a symmetric positive semidefinite matrix as per the definition. In accordance with the principles of matrix eigenvector decomposition, the structure tensor is decomposed to be

Sρ=s11s12s12s22=υ1υ2.λ100λ2.υ1Tυ2T(9)

For two eigenvectors solved, assume λ1 > λ2, i.e. a positive λ1 and negative λ2.

λ1=12s11+s22+(s11−s22)2+4s122λ2=12s11+s22−(s11−s22)2+4s122(10)

Two corresponding unit orthogonal eigenvectors ν1 and ν2 are

{ω=(2S12S22−S11+(S22−S11)2+4S122)ν1=ν2=ω|ω|ν|2(11)

In the local coordinate system composed of two eigenvectors at (x, y), v1 represents the direction in parallel to ∇ν with the maximum rate of change or the direction with the maximum contrast in the geometric structure, which corresponds to the vertical direction of the structure. ν2 stands for the direction perpendicular to ∇v with the minimum rate of change or the direction with the minimum contrast, which corresponds to the direction of the structure. There are three scenarios on seismic images.

  1. λ1≈λ2 ≈0 indicates blank reflections or similar reflection strength in a fat region.

  2. λ1 > > λ2 ≈0 indicates marginal reflections caused by faults, pinch-outs, etc.

  3. λ1≥λ2 > > 0 indicates reflection termination in angular regions.

Local features of each pixel in the image are described by the variations in two directions. Here we define a coherence parameter for quantitative description of local contrast.

H=(λ1−λ1)2=(s11−s12)2+4s122(12)
H indicates local contrast of pixel.

2.3 NLM filtering based on structure tensor

A seismic profile or attribute slice has two properties, numerical property of the pixel itself and geometric structure of neighboring pixels. Large local contrast usually occurs at the marginal structures and faults, while small local contrast occurs in the flat regions. Hence, it is inadequate to involve only numerical similarity in weighting factor calculation for neighborhood similarity. According to the definition, structure tensor describes local geometric properties of the image. Therefore, the coherence parameter H which describes local contrast should be considered in weighting factor calculation.

In accordance with the idea of bilateral filter [15], the weighting factor w(i, j) for NLM filtering algorithm based on structure tensor is defined as follows.

w(i,j)=1Z(i)exp−D2(i,j)h2.exp−H2(i,j)h2(13)

The normalized factor is

Z(i)=∑jexp−D2(i,j)h2.exp−H2(i,j)h2(14)

Similar to D2(i, j) in Equation (5), H2(i, j) is defined by the following expression.

H2(i,j)=H(Ni)−H(Nj)2,a2=∑lnlGa(l)(H(Ni(l))−H(Nj(l)))2(15)

In practical application, a geometric similarity control factor δ [16] is designed in weighting factor calculation to adjust H2(i, j), then w(i, j) is defined to be

w(i,j)=1Z(i)exp−D2(i,j)h2.exp−δH2(i,j)h2(16)

Note that δ∈[0,1]. When the δ is equal to 0, the algorithm takes the original form of NLM filtering. An exponential factor ∈ (0, 1] is added compared with a NLM algorithm. As per Equation (15), the exponential factor is equal to 1 when two neighborhoods have identical H, i.e. identical geometric structures (H2(i, j)=0), when weighting factors are mainly dependent on the numerical similarity of the neighborhoods.

A seismic data volume differs greatly from a conventional image in data range. The order of amplitude variations of seismic data processed by different operators may reach n-th power of 10. Consequently, the coherence parameter H also has large data range. The filtering parameter h and control factor δ should be adjusted as per the scale and coherence parameter of seismic data. The control factor δ is usually very small to offset the impacts of large H2(i, j). The factor δ adjusts the contribution of geometric properties to weighting factors to generate different results of filtering. An increase in δ makes the exponential factor, as well as weighting factors decrease. Such process may preserve geometric properties but cannot eliminate noises sufficiently. This parameter should be adjusted as per the requirements of data processing, which would be discussed in the following model tests.

According to Equation (16), neighborhood similarity involves numerical similarity and geometric similarity. If such a neighborhood is searched to match the original neighborhood, it would be assigned with a large weighting factor. Thus geometric properties are better preserved after processing.

The steps are detailed as follows.

  1. Define the size of S×S, e.g. 21×21, for a region of search and the neighborhood size (for example, 3×3 or 5×5). Input a proper Gaussian function standard deviation a, filtering parameter h and control factor δ.

  2. For boundary data filtering, use mirror reflections to extend the original boundary of an image as per the size of neighborhood. For example, if the neighborhood is of 5×5, the boundary would be extended with 2 data.

  3. Conduct Gaussian smoothing for all elements in the image with the operator G(0, σ2) in accordance with neighborhood size.

  4. Use the Sobel operator or central difference to calculate partial derivative.

  5. Convolve the Gaussian filter G(0, ρ2) and ∂vσ∂x2,∂vσ∂y2 and ∂vσ∂x.∂vσ∂y, respectively to obtain three elements S11, S12 and S22 of the structure tensor matrix.

  6. Calculate coherence parameter H as per Equation (12) to describe geometric similarity.

  7. For any element to be evaluated in the image, calculate the distance D¯2 (i, j) between the neighborhood of the element to be evaluated and that of match element and local contrast H2(i, j) as per neighborhood size.

  8. In terms of the filtering parameter h and control factor δ of geometric similarity, use Equation (16) to calculate weighting factor w(i, j) of all neighborhood centers within the region of search in respect to the element to be evaluated.

  9. Use Equation (2) to calculate the weighted average of all neighborhood centers within the region of search to obtain the value of the element to be evaluated.

  10. Fulfill the filtering of all elements. If the result is not acceptable, return to step (1) to adjust parameters for recalculation until the result is satisfactory.

In this workflow, the coherence parameter H of each element is calculated first because the structure tensor is only related to data distribution in its neighborhood.

3 Case studies

In order to evaluate the results of NLM filtering based on structure tensor, a layered model was designed with faults and wedge-like pinch-outs (Figure 3). Interval velocities were also defined. The underground structure was gridded to be a discrete model with 128 traces and 200 sampling points per trace. The synthetic seismic profile (Figure 4) was made with a Ricker wavelet designed with the dominant frequency of 30 Hz, 64 sampling points and sampling interval of 2 ms. Figure 5 shows the record with 10% random noises. Figure 6 shows the coherence calculated from the noisy model data (Figure 5). The coherence attribute is reconciled with the geometric properties of model data and thus could be used to describe local contrast.

Figure 2: The schematic diagram of elliptic model about structure tensor
Figure 2:

The schematic diagram of elliptic model about structure tensor

Figure 3: Multi-layer geological model contained fault and wedge pinch out
Figure 3:

Multi-layer geological model contained fault and wedge pinch out

Figure 4: Synthetic seismic record
Figure 4:

Synthetic seismic record

Figure 5: Seismic record of added 10% random noise
Figure 5:

Seismic record of added 10% random noise

Figure 6: Result of coherence measure
Figure 6:

Result of coherence measure

Figure 7 shows the result after Gaussian filtering of noisy model data (Figure 5). The neighborhood size is 11×11, and variance is 1. Figure 9 shows the result of NLM filtering with the region of search 21×21, neighborhood 11×11, Gaussian weighted variance 1 and filtering parameter h=50. Figure 11 shows the result of NLM filtering based on structure tensor with the region of search 21×21, neighborhood 11×11, Gaussian weighted variance 1, filtering parameter h=50 and δ = 2.5 × 10−12. Figures 8, 10 and 12 show noise images derived from three algorithms respectively. Gaussian filtering has eliminated noises but also over smoothed events and faults and pinch-out reflections, so this algorithm is not an ideal method for seismic data filtering. NLM filtering and NLM filtering based on structure tensor have generated very similar results (Figure 9 and Figure 11). The former has eliminated noises effectively (Figure 10), however, the latter has also effectively protected events and fault and pinch-out reflections with apparent geometric properties in addition to denoising (Figure 12). Furthermore, the result of the latter (Figure 11) has a higher peak S/N of 59.2443 compared with the value of 57.3538 for the former (Figure 9). In conclusion, in comparison with NLM filtering, the algorithm of NLM filtering based on structure tensor has better performance of noise suppression and edge preservation.

Figure 7: Gauss filtering's result
Figure 7:

Gauss filtering's result

Figure 8: Noise of Gauss filtering
Figure 8:

Noise of Gauss filtering

Figure 9: Filtering result of NLM algorithm
Figure 9:

Filtering result of NLM algorithm

Figure 10: Noise of NLM algorithm
Figure 10:

Noise of NLM algorithm

Figure 11: Filtering result of NLM algorithm based on structure tensor(δ = 2.5 × 10−12)
Figure 11:

Filtering result of NLM algorithm based on structure tensor(δ = 2.5 × 10−12)

Figure 12: Noise of NLM algorithm based on structure tensor(δ =2.5 × 10−12)
Figure 12:

Noise of NLM algorithm based on structure tensor(δ =2.5 × 10−12)

In field data processing by NLM filtering based on structure tensor (Figure 13), the seismic profile has 551 traces and 1000 sampling points per trace. The sampling interval is 2 ms. Processing parameters are: (1) the region of search 11×11, neighborhood 5×5; (2) Gaussian weighted variance 1; (3) filtering parameter h = 1000 and δ = 10−18. The results of filtering and noise image are shown in Figures 14 and 15, respectively. It can be seen that after filtering the noise is depressed effectively, while the texture of seismic profile is well preserved, especially fault reflections have been improved significantly, with clearer fault boundaries and fault combination features (white circle area in Figure 13 and Figure 14), and event continuity has been enhanced. The part (white rectangular area in Figure 13) that can be interpreted falsely as a weakcontinuous reflection in the noise profile, shows the feature of lens shape after filtering, which is closer to the geometry structure of the geological body (white rectangular area in Figure 14). It is favorable for the identification of river and beach bars.

Figure 13: Practical seismic profile contained noise
Figure 13:

Practical seismic profile contained noise

Figure 14: Filtering result of seismic profile
Figure 14:

Filtering result of seismic profile

Figure 15: Noise of NLM algorithm based on structure tensor
Figure 15:

Noise of NLM algorithm based on structure tensor

Figure 16 shows a seismic amplitude slice with two suspect channels, in which a large meandering river channel is shown on the left of the center and a small meandering river channel is shown in the lower left, and this figure is contaminated by some noises, which is not good for a study on characteristics of seismic topography. After NLM filtering based on structure tensor with the region of search 7×7, a neighborhood of 5×5, the Gaussian weighted variance 1, filtering parameter of h=1000, and δ = 10−9. The result (Figure 17) and noise image (Figure 18) show improved S/N and enhanced channel boundaries and structures. Such result provides more evidence for further studying the geomorphic features and seismic facies.

Figure 16: Seismic amplitude slice contained noise
Figure 16:

Seismic amplitude slice contained noise

Figure 17: Filtering result of seismic amplitude slice
Figure 17:

Filtering result of seismic amplitude slice

Figure 18: Noise of NLM algorithm based on structure tensor
Figure 18:

Noise of NLM algorithm based on structure tensor

Buades [10] demonstrated that NLM filtering with neighborhood size of 5×5, 7×7, 9x9 and 11×11 had a good anti-noise capacity and also preserved the details or geometric properties of the original image. For seismic data filtering, the filtering parameters should be selected according to the range of amplitude variations due to the very large amplitude variations of seismic data.

4 Conclusions

The calculation of neighborhood similarity distance involves numerical similarity and neighborhood geometric distance. If some neighborhoods are searched to match the original neighborhood, the one with the smallest geometric distance would be assigned with the largest weighting factor. The coherence attribute, designed to describe local contrast, is adopted to enhance the sensitivity of distance to noises in the region with small contrast. Hence, the weighting factor could be adjusted by NLM denoising algorithm for the center of each similar neighborhood in a noisy region with similar reflection strength. Our test of model data prove that the algorithm could preserve the structures of the image and generate high peak S/N, especially for the seismic data with significant geometric properties. In addition, this algorithm has been demonstrated practically and credibly in processing field data.

Tel.: +86-028-84073252; Fax: +86-028-84078992

Acknowledgement:

The authors are grateful to the anonymous reviewers whose suggestions have been helpful in preparation of the manuscript, and thank editor for carefully revising the manuscript.

References

[1] Yilmaz, Ö., Seismic data analysis: Processing, inversion, and interpretation of seismic data, 2nd ed.: SEG., 200110.1190/1.9781560801580Search in Google Scholar

[2] Li Z.C., Zhang Y.C., Wang Q.Z., Joint denosing method based on wavelet transform and multi-level median filtering. Geophysical Prospecting for petroleum, 2009, 48(5), 470-474 (in Chinese with English abstract)Search in Google Scholar

[3] Chen X.P., Cao S.Y., Second-generation wavelet transform and its application in the de-noising of seismic data. Geophysical Prospecting for petroleum, 2004, 43(6), 547-550 (in Chinese with English abstract)Search in Google Scholar

[4] Peng C., Chang Z., Zhu S.J., Noise elimination method based on curvelet transform. Geophysical Prospecting for petroleum, 2008, 47(5), 461-464(in Chinese with English abstract)Search in Google Scholar

[5] Henning H., Thierry C., David L.M., Erika A., Pierre L. and Didier L., On the use of Geostatistical Filetring Techniques in seismic processing. 73rd SEG International Exposition & Annual Meeting, Dallas, 2003, 7-11Search in Google Scholar

[6] Serra J., Image Analysis and Mathematical Morphology. Academic Press, New York, 1982Search in Google Scholar

[7] Wang J., Chen Y.H., Xu D.H., Qiao Y.L., Structure-oriented edge-preserving smoothing based on accurate estimation of orientation and edges. Applied Geophysics, 2009, 6(eq4), 367-37610.1007/s11770-009-0038-6Search in Google Scholar

[8] Weichert J., Anisotropic diffusion in image processing. PhD thesis, University of Kaiserslautern, Germany, 1996Search in Google Scholar

[9] Buades A., Coll B., Morel J.M., A review of image denoising algorithms, with a new one. Multiscale Modeling and Simulation, 2005, 4, 490-53010.1137/040616024Search in Google Scholar

[10] Buades A., Coll B., Morel J.M., Image denoising methods. a new nonlocal principle. SIAM Review, 2010, 52, 113-14710.1137/090773908Search in Google Scholar

[11] Buades A., Chien A., Morel J.M., Osher S., Topology preserving linear filtering applied to medical imaging. SIAM Journal on Imaging Science, 2008, 1(1), 26-5010.1137/07070022XSearch in Google Scholar

[12] David B., Mauricio S., Denoising seismic data using the nonlocal means algorithm. Geophysics, 2012, 77(1), A5-A810.1190/geo2011-0235.1Search in Google Scholar

[13] Wang W.H., Guo X.B., Shi Y., Non-local means filtering denoising approach and its application in VSP reverse time migration seismic data. Geophysical Prospecting for petroleum, 2015, 54(2), 165-171(in Chinese with English abstract)Search in Google Scholar

[14] Buades A., Coll B., Morel J.M, Nonlocal image and movie denoising. International Journal of Computer Vision, 2008, 76(2), 123-13910.1007/s11263-007-0052-1Search in Google Scholar

[15] Tomasi C., Manduchi R., Bilateral filtering for gray and color images. Sixth International Conference on Computer Vision, Bombay, 1998, 839-846Search in Google Scholar

[16] Xu J., Research of image denoising using Non-Local Means algorithm. Master thesis, Nanjing University of Science and Technology, China, 2009(in Chinese)Search in Google Scholar

Received: 2016-9-22
Accepted: 2017-2-3
Published Online: 2017-5-25

© 2017 S. Yang et al.

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.

Articles in the same Issue

  1. Regular Articles
  2. Two types of gabbroic xenoliths from rhyolite dominated Niijima volcano, northern part of Izu-Bonin arc: petrological and geochemical constraints
  3. Regular Articles
  4. CRSP, numerical results for an electrical resistivity array to detect underground cavities
  5. Regular Articles
  6. Magma evolution inside the 1631 Vesuvius magma chamber and eruption triggering
  7. Regular Articles
  8. Probabilities of Earthquake Occurrences along the Sumatra-Andaman Subduction Zone
  9. Regular Articles
  10. Modelling of carrying capacity in National Park - Fruška Gora (Serbia) case study
  11. Regular Articles
  12. Knickzone Extraction Tool (KET) – A new ArcGIS toolset for automatic extraction of knickzones from a DEM based on multi-scale stream gradients
  13. Regular Articles
  14. When the display matters: A multifaceted perspective on 3D geovisualizations
  15. Regular Articles
  16. Dependence of Gully Networks on Faults and Lineaments Networks, Case Study from Hronska Pahorkatina Hill Land
  17. Regular Articles
  18. Environmental Geochemistry of Geophagic Materials from Free State Province in South Africa
  19. Regular Articles
  20. Neotectonic interpretations and PS-InSAR monitoring of crustal deformations in the Fujian area of China
  21. Regular Articles
  22. Dual-shale-content method for total organic carbon content evaluation from wireline logs in organic shale
  23. Regular Articles
  24. The dolerite dyke swarm of Mongo, Guéra Massif (Chad, Central Africa): Geological setting, petrography and geochemistry
  25. Regular Articles
  26. Seismic data filtering using non-local means algorithm based on structure tensor
  27. Regular Articles
  28. Pore Distribution Characteristics of the Igneous Reservoirs in the Eastern Sag of the Liaohe Depression
  29. Regular Articles
  30. Three-dimensional structural model of the Qaidam basin: Implications for crustal shortening and growth of the northeast Tibet
  31. Regular Articles
  32. Failure Mode of the Water-filled Fractures under Hydraulic Pressure in Karst Tunnels
  33. Regular Articles
  34. Creating a low carbon tourism community by public cognition, intention and behaviour change analysisa case study of a heritage site (Tianshan Tianchi, China)
  35. Regular Articles
  36. Mapping Mangrove Density from Rapideye Data in Central America
  37. Regular Articles
  38. Marine sediment cores database for the Mediterranean Basin: a tool for past climatic and environmental studies
  39. Regular Articles
  40. Retrofitting the Low Impact Development Practices into Developed Urban areas Including Barriers and Potential Solution
  41. Regular Articles
  42. Spatial uncertainty of a geoid undulation model in Guayaquil, Ecuador
  43. Regular Articles
  44. Structure and Filling Characteristics of Paleokarst Reservoirs in the Northern Tarim Basin, Revealed by Outcrop, Core and Borehole Images
  45. Regular Articles
  46. Ground volume assessment using ’Structure from Motion’ photogrammetry with a smartphone and a compact camera
  47. Regular Articles
  48. Classification of coal seam outburst hazards and evaluation of the importance of influencing factors
  49. Regular Articles
  50. Geochemical characterization of Neogene sediments from onshore West Baram Delta Province, Sarawak: paleoenvironment, source input and thermal maturity
  51. Regular Articles
  52. Influence of Social-economic Activities on Air Pollutants in Beijing, China
  53. Regular Articles
  54. Spectral properties of weathered and fresh rock surfaces in the Xiemisitai metallogenic belt, NW Xinjiang, China
  55. Regular Articles
  56. Geochemistry of sandstones and shales from the Ecca Group, Karoo Supergroup, in the Eastern Cape Province of South Africa: Implications for provenance, weathering and tectonic setting
  57. Regular Articles
  58. Petrology and geochemistry of meta-ultramafic rocks in the Paleozoic Granjeno Schist, northeastern Mexico: Remnants of Pangaea ocean floor
  59. Regular Articles
  60. Distal turbidite fan/lobe succession of the Late Oligocene Zuberec Fm. – architecture and hierarchy (Central Western Carpathians, Orava–Podhale basin)
  61. Regular Articles
  62. Fourier Transform Infrared Spectroscopy of Clay Size Fraction of Cretaceous-Tertiary Kaolins in the Douala Sub-Basin, Cameroon
  63. Regular Articles
  64. Optimized AVHRR land surface temperature downscaling method for local scale observations: case study for the coastal area of the Gulf of Gdańsk
  65. Regular Articles
  66. New non-linear model of groundwater recharge: Inclusion of memory, heterogeneity and visco-elasticity
  67. Regular Articles
  68. “Urban geosites” as an alternative geotourism destination - evidence from Belgrade
  69. Regular Articles
  70. A customized resistivity system for monitoring saturation and seepage in earthen levees: installation and validation
  71. Regular Articles
  72. Consideration of Landsat-8 Spectral Band Combination in Typical Mediterranean Forest Classification in Halkidiki, Greece
  73. Regular Articles
  74. Coda Wave Attenuation Characteristics for North Anatolian Fault Zone, Turkey
  75. Regular Articles
  76. Modal composition and tectonic provenance of the sandstones of Ecca Group, Karoo Supergroup in the Eastern Cape Province, South Africa
  77. Regular Articles
  78. Quantitative studies of the morphology of the south Poland using Relief Index (RI)
  79. Regular Articles
  80. Interpretation of sedimentological processes of coarse-grained deposits applying a novel combined cluster and discriminant analysis
  81. Regular Articles
  82. Utilizing borehole electrical images to interpret lithofacies of fan-delta: A case study of Lower Triassic Baikouquan Formation in Mahu Depression, Junggar Basin, China
  83. Regular Articles
  84. Grain size statistics and depositional pattern of the Ecca Group sandstones, Karoo Supergroup in the Eastern Cape Province, South Africa
  85. Regular Articles
  86. Carbonate stable isotope constraints on sources of arsenic contamination in Neogene tufas and travertines of Attica, Greece
  87. Regular Articles
  88. Appreciation of landscape aesthetic values in Slovakia assessed by social media photographs
  89. Regular Articles
  90. Geochemistry of Selected Kaolins from Cameroon and Nigeria
  91. Regular Articles
  92. Spatial pattern of ASG-EUPOS sites
  93. Regular Articles
  94. A Stream Tilling Approach to Surface Area Estimation for Large Scale Spatial Data in a Shared Memory System
  95. Regular Articles
  96. A location-based multiple point statistics method: modelling the reservoir with non-stationary characteristics
  97. Regular Articles
  98. Water Inrush Analysis of the Longmen Mountain Tunnel Based on a 3D Simulation of the Discrete Fracture Network
  99. Regular Articles
  100. A Computer Program for Practical Semivariogram Modeling and Ordinary Kriging: A Case Study of Porosity Distribution in an Oil Field
  101. Regular Articles
  102. Imaging and locating paleo-channels using geophysical data from meandering system of the Mun River, Khorat Plateau, Northeastern Thailand
  103. Regular Articles
  104. Rare earth element contents of the Lusi mud: An attempt to identify the environmental origin of the hot mudflow in East Java – Indonesia
  105. Regular Articles
  106. Is Nigeria losing its natural vegetation and landscape? Assessing the landuse-landcover change trajectories and effects in Onitsha using remote sensing and GIS
  107. Regular Articles
  108. Methodological approach for the estimation of a new velocity model for continental Ecuador
https://doi.org/10.1515/geo-2017-0013
Keywords for this article
Non-Local means; Structure tensor; Filter; Similarity; Interpretative processing
Creative Commons
BY-NC-ND 3.0

Articles in the same Issue

  1. Regular Articles
  2. Two types of gabbroic xenoliths from rhyolite dominated Niijima volcano, northern part of Izu-Bonin arc: petrological and geochemical constraints
  3. Regular Articles
  4. CRSP, numerical results for an electrical resistivity array to detect underground cavities
  5. Regular Articles
  6. Magma evolution inside the 1631 Vesuvius magma chamber and eruption triggering
  7. Regular Articles
  8. Probabilities of Earthquake Occurrences along the Sumatra-Andaman Subduction Zone
  9. Regular Articles
  10. Modelling of carrying capacity in National Park - Fruška Gora (Serbia) case study
  11. Regular Articles
  12. Knickzone Extraction Tool (KET) – A new ArcGIS toolset for automatic extraction of knickzones from a DEM based on multi-scale stream gradients
  13. Regular Articles
  14. When the display matters: A multifaceted perspective on 3D geovisualizations
  15. Regular Articles
  16. Dependence of Gully Networks on Faults and Lineaments Networks, Case Study from Hronska Pahorkatina Hill Land
  17. Regular Articles
  18. Environmental Geochemistry of Geophagic Materials from Free State Province in South Africa
  19. Regular Articles
  20. Neotectonic interpretations and PS-InSAR monitoring of crustal deformations in the Fujian area of China
  21. Regular Articles
  22. Dual-shale-content method for total organic carbon content evaluation from wireline logs in organic shale
  23. Regular Articles
  24. The dolerite dyke swarm of Mongo, Guéra Massif (Chad, Central Africa): Geological setting, petrography and geochemistry
  25. Regular Articles
  26. Seismic data filtering using non-local means algorithm based on structure tensor
  27. Regular Articles
  28. Pore Distribution Characteristics of the Igneous Reservoirs in the Eastern Sag of the Liaohe Depression
  29. Regular Articles
  30. Three-dimensional structural model of the Qaidam basin: Implications for crustal shortening and growth of the northeast Tibet
  31. Regular Articles
  32. Failure Mode of the Water-filled Fractures under Hydraulic Pressure in Karst Tunnels
  33. Regular Articles
  34. Creating a low carbon tourism community by public cognition, intention and behaviour change analysisa case study of a heritage site (Tianshan Tianchi, China)
  35. Regular Articles
  36. Mapping Mangrove Density from Rapideye Data in Central America
  37. Regular Articles
  38. Marine sediment cores database for the Mediterranean Basin: a tool for past climatic and environmental studies
  39. Regular Articles
  40. Retrofitting the Low Impact Development Practices into Developed Urban areas Including Barriers and Potential Solution
  41. Regular Articles
  42. Spatial uncertainty of a geoid undulation model in Guayaquil, Ecuador
  43. Regular Articles
  44. Structure and Filling Characteristics of Paleokarst Reservoirs in the Northern Tarim Basin, Revealed by Outcrop, Core and Borehole Images
  45. Regular Articles
  46. Ground volume assessment using ’Structure from Motion’ photogrammetry with a smartphone and a compact camera
  47. Regular Articles
  48. Classification of coal seam outburst hazards and evaluation of the importance of influencing factors
  49. Regular Articles
  50. Geochemical characterization of Neogene sediments from onshore West Baram Delta Province, Sarawak: paleoenvironment, source input and thermal maturity
  51. Regular Articles
  52. Influence of Social-economic Activities on Air Pollutants in Beijing, China
  53. Regular Articles
  54. Spectral properties of weathered and fresh rock surfaces in the Xiemisitai metallogenic belt, NW Xinjiang, China
  55. Regular Articles
  56. Geochemistry of sandstones and shales from the Ecca Group, Karoo Supergroup, in the Eastern Cape Province of South Africa: Implications for provenance, weathering and tectonic setting
  57. Regular Articles
  58. Petrology and geochemistry of meta-ultramafic rocks in the Paleozoic Granjeno Schist, northeastern Mexico: Remnants of Pangaea ocean floor
  59. Regular Articles
  60. Distal turbidite fan/lobe succession of the Late Oligocene Zuberec Fm. – architecture and hierarchy (Central Western Carpathians, Orava–Podhale basin)
  61. Regular Articles
  62. Fourier Transform Infrared Spectroscopy of Clay Size Fraction of Cretaceous-Tertiary Kaolins in the Douala Sub-Basin, Cameroon
  63. Regular Articles
  64. Optimized AVHRR land surface temperature downscaling method for local scale observations: case study for the coastal area of the Gulf of Gdańsk
  65. Regular Articles
  66. New non-linear model of groundwater recharge: Inclusion of memory, heterogeneity and visco-elasticity
  67. Regular Articles
  68. “Urban geosites” as an alternative geotourism destination - evidence from Belgrade
  69. Regular Articles
  70. A customized resistivity system for monitoring saturation and seepage in earthen levees: installation and validation
  71. Regular Articles
  72. Consideration of Landsat-8 Spectral Band Combination in Typical Mediterranean Forest Classification in Halkidiki, Greece
  73. Regular Articles
  74. Coda Wave Attenuation Characteristics for North Anatolian Fault Zone, Turkey
  75. Regular Articles
  76. Modal composition and tectonic provenance of the sandstones of Ecca Group, Karoo Supergroup in the Eastern Cape Province, South Africa
  77. Regular Articles
  78. Quantitative studies of the morphology of the south Poland using Relief Index (RI)
  79. Regular Articles
  80. Interpretation of sedimentological processes of coarse-grained deposits applying a novel combined cluster and discriminant analysis
  81. Regular Articles
  82. Utilizing borehole electrical images to interpret lithofacies of fan-delta: A case study of Lower Triassic Baikouquan Formation in Mahu Depression, Junggar Basin, China
  83. Regular Articles
  84. Grain size statistics and depositional pattern of the Ecca Group sandstones, Karoo Supergroup in the Eastern Cape Province, South Africa
  85. Regular Articles
  86. Carbonate stable isotope constraints on sources of arsenic contamination in Neogene tufas and travertines of Attica, Greece
  87. Regular Articles
  88. Appreciation of landscape aesthetic values in Slovakia assessed by social media photographs
  89. Regular Articles
  90. Geochemistry of Selected Kaolins from Cameroon and Nigeria
  91. Regular Articles
  92. Spatial pattern of ASG-EUPOS sites
  93. Regular Articles
  94. A Stream Tilling Approach to Surface Area Estimation for Large Scale Spatial Data in a Shared Memory System
  95. Regular Articles
  96. A location-based multiple point statistics method: modelling the reservoir with non-stationary characteristics
  97. Regular Articles
  98. Water Inrush Analysis of the Longmen Mountain Tunnel Based on a 3D Simulation of the Discrete Fracture Network
  99. Regular Articles
  100. A Computer Program for Practical Semivariogram Modeling and Ordinary Kriging: A Case Study of Porosity Distribution in an Oil Field
  101. Regular Articles
  102. Imaging and locating paleo-channels using geophysical data from meandering system of the Mun River, Khorat Plateau, Northeastern Thailand
  103. Regular Articles
  104. Rare earth element contents of the Lusi mud: An attempt to identify the environmental origin of the hot mudflow in East Java – Indonesia
  105. Regular Articles
  106. Is Nigeria losing its natural vegetation and landscape? Assessing the landuse-landcover change trajectories and effects in Onitsha using remote sensing and GIS
  107. Regular Articles
  108. Methodological approach for the estimation of a new velocity model for continental Ecuador
Sign up now to receive a 20% welcome discount
Institutional Access
How does access work?
Have an idea on how to improve our website?
Please write us.
© 2026 De Gruyter Brill
Downloaded on 14.1.2026 from https://www.degruyterbrill.com/document/doi/10.1515/geo-2017-0013/html
Scroll to top button