# DeStripe: frequency-based algorithm for removing stripe noises from AFM images

Shu-Wen W Chen, Jean-Luc Pellequer

► To cite this version:

Shu-Wen W Chen, Jean-Luc Pellequer. DeStripe: frequency-based algorithm for removing stripe noises from AFM images. BMC Structural Biology, 2011, 11 (1), pp.7. 10.1186/1472-6807-11-7. hal-04268279

HAL Id: hal-04268279

https://hal.science/hal-04268279v1

Submitted on 2 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

# DeStripe: frequency-based algorithm for removing stripe noises from AFM images

Chen and Pellequer
SOFTWARE

Open Access

# DeStripe: frequency-based algorithm for removing stripe noises from AFM images

Shu-wen W Chen, Jean-Luc Pellequer*

## Abstract

**Background:** Atomic force microscopy (AFM) is a relatively recently developed technique that shows a promising impact in the field of structural biology and biophysics. It has been used to image the molecular surface of membrane proteins at a lateral resolution of one nanometer or less. An immediate obstacle of characterizing surface features in AFM images is stripe noise. To better interpret structures at a sub-domain level, pre-processing of AFM images for removing stripe noises is necessary. Noise removal can be performed in either spatial or frequency domain. However, denoising processing in the frequency domain is a better solution for preserving edge sharpness.

**Results:** We have developed a denoising protocol, called *DeStripe*, for AFM bio-molecular images that are contaminated with heavy and fine stripes. This program adopts a divide-and-conquer approach by dividing the Fourier spectrum of the image into central and off-center regions for noisy pixels detection and intensity restoration; it is also applicable to other images interfered with high-density stripes such as those acquired by the scanning electron microscope. The denoising effect brought by *DeStripe* provides better visualization for image objects without introducing additional artifacts into the restored image.

**Conclusions:** The DeStripe denoising effect on AFM images is illustrated in the present work. It allows extracting extended information from the topographic measurements and implicitly enhances the molecular features in the image. All the presented images were processed by *DeStripe* with the raw image as the only input without any requirement for other prior information. A web service, http://biodev.cea.fr/destripe, is available for running *DeStripe*.

## Background

Unlike other optical-based microscopes, atomic force microscope (AFM) is a sensing instrument [1]. In brief, a nano-sized tip located beneath a micro-cantilever scans across a field of deposited molecules. The cantilever deflection can be detected by a laser beam that reflects off the back of the cantilever. With a set of piezoelectric ceramics connected to the cantilever, the so-called height image of molecules can be made at a constant applied force. Because of its exceptional high signal/noise ratio, AFM is able to measure the topography of a single isolated molecule with a lateral resolution of a few nanometers and a vertical resolution of a few Angstroms [2].

Acquisition of high-resolution images of macromolecules in aqueous solution using AFM does not require sample staining [3]. Development of AFM imaging techniques in life sciences is progressing [4] including imaging single isolated molecules at high speed [5]. To date very high-resolution imaging by AFM has been obtained on membrane proteins 2 D crystals [6,7] as well as on densely packed proteins in native membranes [8-10]. Although AFM is primarily an imaging tool, it also allows measuring inter- and intra-molecular interactions on the pico-Newton scale [11-14] and force-probing the surfaces of living cells at the single molecule level [15,16] in order to map protein receptors for example [17].

Image quality is highly involved in feature interpretation and extraction for sampled objects in all kinds of imaging systems. Noise is a critical artifact that influences image quality and is mainly produced during image acquisition. It needs to be removed or reduced

<footnote>* Correspondence: jlpellequer@cea.fr
CEA, iBEB, Service de Biochimie et Toxicologie Nucléaire, F-30207 Bagnols sur Cèze, France</footnote>
for further data processing to acquire desirable information or feature interpretation. There is no exception for AFM images. Among all types of noises, stripe noise is the notorious one that profoundly degrades image quality acquired by AFM, and is a consequence of the scanning pattern. During the scanning along one line of image, the sample surface height is acquired by an oscillating motion of the cantilever tip in the perpendicular direction to the substrate plate. Stripe noise may occur, for example using the tapping mode of AFM [18], from a loss or inadequate acquisition of height information. Any abrupt increase of the force from samples exerted on the tip would make a dramatic change in the tip vibration such that the noise cannot stay constant during the scanning; moreover these noise errors cannot be averaged off. A critical factor causing this change is the interaction between the sample and the tip. Thereby, these stripes vary in intensity, length, local density and frequency range, subjective to a variety of factors, such as sample preparation conditions and the constituents therein on the substrate plate. In particular, sharp and irregular boundaries or protuberances in the object distribution usually produce serious stripes across the image.

Noise is usually modeled as an additive term in the intensity distribution with a Gaussian or Poisson form of zero mean and constant variance [19]. This assumption becomes inapplicable if the image intensity is non-randomly distributed. From the noise origin described above, the high-density and fine stripes observed in AFM images do not occur randomly (neither in time nor in space). An image frequency spectrum can depict the noise characteristics [20]. In traditional strategies of Fourier transform, a filter with high or low pass is often set up to gate the frequencies composed of stripe noises. However, this method sometimes removes some image details when stripe noises are mixed with components from object textures. Recently, a combination of wavelet and Fourier transform has been developed to remove both stripe and ring artifacts from images [21]. Alternatively, a method based on the heterogeneity of image frequency spectrum has been proposed to remove periodic and quasi-periodic stripes [22]. However, as shown in this paper, the stripe noises observed in AFM images are not in a unique or specific form. Although we adopted a similar approach for developing the denoising protocol, some new designs are highlighted for the removal of non-uniform and high-density stripe noises from AFM images.

## Implementation

To be visualized by the human eye, the denoising procedure was not performed straightforwardly on the spectrum amplitude but on its logarithmic scale, the corresponding image of which is called LogF. We

proposed a divide-and-conquer strategy to proceed the denoising in two separate regions of the frequency domain, i.e., the central and the off-center regions, as most of the high intensities are concentrated in the neighborhood of the origin. In order to best preserve the original data and reduce the computational task, we employed a decision-based algorithm [23] to select pixels for the variance test [24] and intensity reconstruction. Heterogeneity measurements were performed for detecting a spectral pixel potentially responsible for stripe noise in the raw AFM image.

Regarding the image restoration, the intensity replacement was also performed in the image spectrum, and the processed image was obtained by the inverse Fourier transformation. In general, a smoothing procedure applied in the image spatial domain can be exactly used in its frequency domain, e.g., a median-like filter [22]. Yet, processing in the frequency domain is a better approach than in the spatial domain for that the edge sharpness can be better preserved. Moreover, we attempted choosing a size for the local window to restore the image intensity as small as possible. Consequently, the accuracy of image restoration is only jeopardized at a minimal degree in the case where pixels are falsely picked up as noisy. The design of our denoising protocol is outlined in Figure 1 and details are described below.

## Heterogeneity function

The heterogeneity of LogF determines whether a pixel is noisy or not. The heterogeneity measurement ranges from 0 to 1; the larger the value, the more heterogeneous is the intensity. There are two components in the heterogeneous function, i.e., abrupt change in intensity and intensity itself. The former is represented by the Laplacian of LogF. Table 1 presents the discrete Laplacian operator used in DeStripe. For each component, the values were offseted and normalized over the entire image such that they are confined within [0, 1]. Accordingly, the heterogeneity function $H$ at the pixel $(i, j)$ is expressed as

$$H(i, j) = \frac{L(i, j) - L_{\min}}{L_{\max} - L_{\min}} \times \frac{I(i, j) - I_{\min}}{I_{\max} - I_{\min}},$$

where $L(i, j)$ and $I(i, j)$ represent the Laplacian and intensity values at pixel $(i, j)$, respectively, $L_{\min}$, $L_{\max}$, $I_{\min}$ and $I_{\max}$ correspond to their maxima and minima. As a result, the pixels with higher $H$ values possess higher intensity and experience more dramatic change in intensity.

## Global sampling of pixels

A preliminary sampling of noisy pixels was done by thresholding the $H$ value based on an internally determined value ($H_{ref}$) extracted from the heterogeneity
<figure><img src="image_3.png" /><figcaption>**Figure 1** The flow chart of DeStripe, see details in Implementation. The dash-lined box describes the steps of reducing the number of potentially noisy pixels. The denotation of pixel sets at various steps is stated in Implementation.</figcaption></figure>

histogram which is analog to the homogeneity level in removal of impulse noises [25]. For convenience, we denoted this set of sampled pixels as P<sub>n</sub>1 and the algorithm is described below:

1. Calculate *H* values of LogF.

2. Form the histogram of *H* with 20 bins such that the variation of *H* in each bin is within 0.05.

<table><caption>**Table 1** The discrete Laplacian operator used in DeStripe</caption><tr><td>-1</td><td>-1</td><td>-1</td></tr><tr><td>-1</td><td>8</td><td>-1</td></tr><tr><td>-1</td><td>-1</td><td>-1</td></tr></table>

3. Find the most spreading group of consecutive bins with non-zero populations in the histogram.

4. From the most spreading group, find the threshold bin in the direction of increasing heterogeneity. The threshold bin is defined as the first bin encountered such that the ratio of its population relative to that of the most populated bin is ≤ 0.5.

5. Calculate the *H* threshold as $H_{ref} = 0.5 (H_{up} + H_{low})$, where $H_{up}$ and $H_{low}$ correspond to the upper and lower limits of the threshold bin.

6. Calculate the intensity threshold as $I_{ref} = 0.5 (max_{ref} + ave_{ref})$, where $ave_{ref}$ and $max_{ref}$ are the average and the maximum of intensities over the
pixels whose $H \le H_{ref}$ respectively. A relevant parameter, variance (var$_{ref}$), is also calculated.

$$7. P_{n1} = \{(i, j) | H(i, j) > H_{ref} \ \& \ I(i, j) > I_{ref}\}.$$

## Divide-and-Conquer Strategy

Due to dramatic variations in intensity in the central region of the frequency domain, we divided Pn1 into two groups; one was referred to as the central region and the other as the off-center region.

### Formation of the central region

The central region was considered as a circular disk. The initial radius of the disk was derived from the moment of inertia tensor of Pn1, where the mass magnitude was replaced by the intensity value in the tensor array [26]. We calculated the eigenvalues ($\sigma_x, \sigma_y$) and eigenvectors ($\hat{e}_x, \hat{e}_y$) of the tensor array. The initial radius equals the square root of $\sigma_x + \sigma_y$. Starting from the center of intensity distribution, ($i_0, j_0$), we expanded stepwise the region outwardly with an increment of 1/10 of the initial radius value. At each expansion step, we counted the Pn1 pixels and the total within the newly expanded region; if the ratio of the two numbers was $\le 0.85$, then the expansion was stopped and all the visited pixels $\in P_{n1}$ were included as members in the central region; we denoted it as $C_0$ and Pn1 - $C_0$ as Pn2.

### Sampling of noisy pixels in the central region

In order to avoid vain data treatments and reduce the number of false noisy pixels recruitment, we modeled the intensity distribution of $C_0$ by an anisotropic Gaussian function. We used the Levenberg-Marquadt algorithm [27] to fit the intensity nonlinearly into the model function. In other words, we minimized the objective

function, $\sum_{i,j} f(i, j)^2$, where $f(i, j)$ was expressed as

$$f(i, j) = \frac{I(i, j)}{I(i_0, j_0)} - e^{-\left\{ -c_1 \frac{(i-i_0)^2}{\sigma_x^2} - c_2 \frac{(j-j_0)^2}{\sigma_y^2} \right\}},$$

where $I(i_0, j_0)$, $c_1$ and $c_2$ are fitting parameters; $I(i_0, j_0)$ is the restored intensity at ($i_0, j_0$), and $c_1$ and $c_2$ shape the anisotropic breadth of the Gaussian fitting; $I(i, j)$ is the intensity at ($i, j$). The relationship between ($\sigma_x', \sigma_y'$) and ($\sigma_x, \sigma_y$) can be obtained through ($\hat{e}_x, \hat{e}_y$) as [28]

$$\sigma_x' = \sigma_x \times \cos^2 \theta + \sigma_y \times \sin^2 \theta, \text{ and}$$

$$\sigma_y' = \frac{\sigma_x \times \sigma_y}{\sigma_x'}$$

where $\theta$ is the angle that rotates counterclockwise the Cartesian x, y-coordinate axes aligned with ($\hat{e}_x, \hat{e}_y$). Note that ($\sigma_x, \sigma_y$) and ($\hat{e}_x, \hat{e}_y$) were computed based on the right-handed rule while ($\sigma_x', \sigma_y'$) were defined according to the convention of image presentation, i.e., the origin is located at the left-top corner of the image, the image j row and i column represent the x- and y-axes, respectively. For a pixel ($i, j) \in C_0$, if $f(i, j) \le 0$, then the pixel was not considered as corrupted with noise, otherwise it was included in a set denoted by $C_{n1}$.

We selected data points in $C_{n1}$ with a 10-bin histogram using the same thresholding method as described previously for global sampling. The collected pixels were further screened and clustered based on two criteria: (1) if the pixel distribution in the form of horizontal or vertical line was > 2/3 of the region of interest; (2) if the length of consecutive pixels was $\ge 4$. The qualified pixels with either condition were considered as noisy and formed the set, $C_{n2}$. The same procedure was also applied to the off-center region, Pn2; the resulting pixel set is denoted by Pn3. Maps of these sets are presented in additional file 1 Figure S1, for each study image.

### Constrained variance (CVAR) test and intensity restoration

For a processed pixel ($i_c, j_c$), the local mean and variance for the CVAR test were calculated as:

$$ave = \frac{1}{N_{\text{var}}} \sum_{i=-N_s}^{N_s} \sum_{j=-N_s}^{N_s} I(i + i_c, j + j_c),$$

$$\sigma_{\text{var}} = \frac{1}{N_{\text{var}}} \sum_{i=-N_s}^{N_s} \sum_{j=-N_s}^{N_s} (I(i + i_c, j + j_c) - ave)^2.$$

The imposed constraint is that only non-noisy pixels were counted for $N_{var}$ within the $(2N_s+1) \times (2N_s+1)$ local window, i.e., ($i+i_c, j+j_c$) neither $\in C_{n2}$ nor $\in P_{n2}$. In the present work, $N_s = 1$. Define $\sigma_{std} = \sqrt{\sigma_{var}}$, if $I(i_c, j_c) - ave > \sigma_{std}$, then $I(i_c, j_c)$ was replaced by ave, otherwise the original $I(i_c, j_c)$ was reserved. For Pn2, the CVAR test was not performed point by point; instead, we clustered the connected pixels and the test was performed starting at the boundary pixels of each cluster.

### Filter function

The major component of DeStripe is the filter function that turns a noisy image into a clean one. Consider $S(i, j)$ and $S(i, j)'$ as the measured and restored intensity values at pixel ($i, j$) in the LogF image, respectively, then the

filter function was calculated as $\Phi(i, j) = \frac{e^{S(i, j)}}{e^{S(i, j)}}$; the
value range is (0, 1], i.e. $0 < \Phi(i, j) \le 1$. Accordingly, the restored and the associated noise images were obtained by the inverse FFT of the products of $\exp[S(i, j)]$ with $\Phi(i, j)$ and $1-\Phi(i, j)$, respectively. The image formed by the $\Phi(i, j)$ values is henceforth called $\Phi$-image.

## Results and Discussion

We present two biological systems of which the topographic images were measured by AFM. One contains biomembranes and the other is constituted of proteins belonging to a large family of hydrolase enzymes, GTPase. We applied *DeStripe* to these raw AFM images for stripe noise reduction. The main purpose of the denoising is to better reveal molecular features in images that are distinguishable to human vision. Lastly, the applicability of *DeStripe* was also tested for stripe removal from images acquired by a scanning electron microscope, SEM [29].

## Denoising vs. high-resolution AFM imaging

In high-resolution AFM imaging, one goal is to measure molecular topographies down to a sub-nanometer scale. In most occasions, one may seek an AFM tip as small as possible for probing the sample surface; the usual radius for AFM tips is 5-10 nm [30]. However, we found that the AFM image quality also profoundly depended on the study system. In other words, one set of instrumentation parameters may not guarantee to gain similar quality for images with different biological constituents, as the sample-related stripe noise is always a major hindrance.

The denoising results of AFM imaging on biomembrane surfaces are shown in Figure 2. Two sets of images are arranged in two columns to indicate that the topographic measurements were independently performed with different experimental preparations and AFM setups. The representation of restored and noise images sharing the same intensity range with raw images is referred to additional file 1 Figure S2. We present here no details in AFM instrumentation or sample preparation, implying no need to know how this image was obtained for denoising. The top row presents the AFM raw images while the second and third rows illustrate the corresponding restored and noise images. The last row shows the $\Phi$-image. Similar presentation for GTPase enzymes is presented in Figure 3.

We found that the Gaussian model was appropriate for fitting the intensity distribution of the central region, and for assigning an appropriate value to the center pixel, ($i_0, j_0$). A strong noise occurred at the center pixel usually leads to the appearance of heavy or notorious stripes in the image. The Gaussian model greatly reduced the number of pixels processed; one may notice much less pixels

restored in the central region than in the off-center region. Naturally, one may speculate that there is more chance to recruit false noisy pixels in the off-center region. Recall that the FFT amplitude in the off-center region is much smaller and more homogeneous than that in the central region. The replacement of intensity values at false noisy pixels may not dramatically affect the image quality; that is the underlined rationale for dividing the candidate pixels into two regions.

In structural biology, the distribution of molecules is one feature of concern for AFM imaging. In general, observable individuality is a prerequisite prior to interpretation of molecular features for imaged objects. One may overlook important details simply because the entire image looks so dim due to the presence of stripe noises with very high intensity. The consequence

<figure><img src="image_4.png" /></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 3** AFM imaging on GTPase enzymes. Numbers at the bottom indicate the number of pixel composing the image. The scan size is 0.25 × 0.25 µm² for the left column and 0.125 × 0.125 µm² for the right column, while the intensity unit for both is 0.1 nm.</figcaption></figure>

of stripe removal can be evidenced by comparison of raw and processed images. First, the visibility of fine structures in the image was enhanced by the protocole. The left system in Figure 2 clearly demonstrates this effect. The restored image revealed that the dark areas in between bright segmental regions in the raw image were in fact distributed with granular particles. Second, the particle shape or cluster form became better observed. After noise reduction, stripes were eliminated or separated into short segments, and some brightly fused regions or islands were resolved into assemblies of individual granular particles; even the individual shape of background particles was noticeable.

Recently, enhancing image contrast by subtraction of a smoothed image from the raw image has been used for better visualizing the image objects [31]. The subtraction method is equivalent to extraction of edge features. In the present study, the purpose of denoising is to remove the noise prior to any other image processing and

enhanced visibility on existent features in the image is a natural result.

## Noise vs. the image quality of AFM

For the same biological molecule but prepared in different conditions, the occurring pattern of noise is dissimilar as illustrated in Figures 2 and 3. Note that stripes can be due to horizontal or vertical noises [32]. In AFM surface measurements, vertical stripe is the major pattern of noise that intensely affects image quality. To a worse degree, their intensity distribution is not uniform, as seen in the raw image on the left of Figure 2. This complicates noise characterization and estimation. Likewise, these inhomogeneous stripes are also observable from the right noise image of Figure 3. In the right image of Figure 3 two observed dark bands are however due to scan lines misalignments along the height ordinate. *DeStripe* does not aim to remove such type of noise, yet it can be eliminated or reduced by flattening the image [31].

The noise images in Figures 2 and 3 present the noise components peeled off from the raw image. Comparing AFM images on different systems, the degree of noise removal by *DeStripe* was found to vary. This reflects in the distribution of potentially noisy pixels and the values in the Φ-image. From the intensity distribution of Φ-image, one may perceive the noisy pixels identified and the degree of intensity modification. The smaller intensity value in the Φ-image, the greater proportion of the spectrum amplitude is removed. The performance of *DeStripe* in denoising can be evaluated by visual inspection of the noise image. Our results show that the stripe noise pattern is very different from one system to the other. It implies that *DeStripe* is able to automatically tune the denoising performance. It is noteworthy that here all the images processed by *DeStripe* use the same set of parameter values; these values were chosen by trial and error such that the denoising can be effective for various images. Consequently, the only user-provided input is the raw image.

Taken together with restored images, we found that *DeStripe* under-denoises somewhere but over-denoises elsewhere in the raw image. On the one hand, it is trivial to diagnose a case of under-denoising if there is any stripe noise visible in the restored image. By lowering the heterogeneity threshold to include more candidate pixels in C<sub>n</sub>2, we found that the visible stripe noises can be further removed from the restored image (data not shown). We have also attempted to use lower values for restored intensities and some stripes diminished. On the other hand, if there is any non-stripe structure pattern in the noise image correlated with object features in the restored image, it may imply that *DeStripe* cut off too much the image intensity. In order to obtain the true
surface measured by AFM, the structure pattern of noise
image can be used as a guide for judging whether the
denoising is appropriate. All the noise images presented
here exhibit almost purely stripe noises, mainly in linear
form. Indeed, a compromise between noise reduction
and preservation of structural features remains
challenging.

By comparing intensity range between restored and
raw images, our results yield a smaller size for restored
images than for the raw ones. Nevertheless, except for
the left system in Figure 2 all other systems show com-
parable sizes for both images. From the results of the
left system in Figure 2 the range of restored intensities
is almost half that of the raw ones; it implies that near
50% of the intensity magnitude measured at some pixels
by AFM is attributed to the noise. Certainly, those topo-
graphic measurements cannot represent a true value for
the molecular surface; one can no longer consider these
topographic measurements as true surface heights if
heavy and bright stripe noises contaminate the image at
such level. Consequently, the denoising is inevitable for
better characterizing the surface feature of AFM images
for a realistic biological system.

Comparisons of denoising effects with Gwyddion

Gwyddion is open-source software for AFM image pro-
cessing [33]. We applied the “correct lines” tool for
comparing the denoising effects on high-density stripe
removal. The results from Gwyddion are presented in
Figure 4. By visual inspection, Gwyddion scarcely
removed any noise from the image presented in row 4;
this is observable either by comparing the restored and
the raw images or from the intensity range of the noise
image. As a matter of fact, the intensity values are essen-
tially zero in the noise image on the right of Figure 3. The
comparison reveals that DeStripe is more effective than
Gwyddion on stripe reduction for AFM images. As a
result, DeStripe manifests better the individuality of objects
observed from the AFM image.

For topographic measurements, we assume that the
noise is a positive quantity overlaid over the molecular
surface; therefore the denoising effect should not yield
restored images with greater intensity values than the
raw measurements. In the present work, the processed
images by DeStripe satisfy this important feature, reveal-
ing no other additive component imposed on the original
image during the denoising procedure. This also reflects
in the positive range of intensity in the noise images. In
contrast, some noise intensities fall within the negative
range in the denoising results from Gwyddion, see Figure 4.
This indicates that Gwyddion may yield some restored
intensities greater than the original data, thus creating new
artifacts in the raw AFM image. Consequently, the two
advantageous attributes, i.e., effective denoising and no

extra artifact, make *DeStripe* superior to *Gwyddion* in
improving the AFM image quality.

Application to SEM image

DeStripe performs denoising using an anisotropic Gaus-
sian function to fit the spectral amplitude distribution
close to the origin of the frequency domain. This mod-
eling has been shown to provide an appropriate estimate
of intensity for the origin of the frequency domain that
is corrupted with noise. To further explore this aspect,
we ran the program for a SEM micrograph. The results
are presented in Figure 5. The originally measured
image is seriously interfered with stripe noises, and we
found that these severe stripes were mainly attributed to
the contribution of the amplitude value at the origin.
This ascertains that DeStripe strategy design is useful
for images corrupted with heavy stripe noises. This
image was processed using the same set of parameter
values as for the AFM images.

Conclusions

AFM emerges as a new nanotechnology for investigating
the structure of biological systems with mesoscopic
sizes, ranging from cellular morphology to single protein

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5** SEM imaging on a sintered specimen of CeO₂ (cerium oxide) at a sintering temperature of 1400°C for 2 hours and 50 minutes. The image is composed of 512 × 512 pixels with each side measuring 70.55 µm; the intensity unit is 0.1 µm.</figcaption></figure>

topography. Although this technique can probe bio-
molecular surface down to a sub-nanometer scale, noise
artifacts produced along with the measurements are also
visible. We aimed to unveil the layer of noise from the
image in order to witness the true topography beneath,
and thus developed the automatic program *DeStripe*. By
ripping off the noise contribution from the measured
intensity, one obtains the surface heights closer to the
true value and is able to inspect the surface features
more easily; implicitly, the denoising processing
enhances the molecular features of concern. *DeStripe*
involves several methods for reducing the number of
pixels for intensity restoration, including global sam-
pling, non-local variance test, Gaussian fitting, and local
CVAR test. Gaussian modeling revealed its superiority
over conventional variance test in identifying the noisy

pixels in the central region where intensity drastically
changes. Moreover, the criterion of non-negative noise
was insistently used for the reason that AFM imaging
was directed for probing surfaces unexplored before;
denoising simply for pleasing human eye is a false
pursuit.

Availability and requirements

• Project name: DeStripe

• Project server: http://biodev.cea.fr/destripe

• Operating system: Linux

• Programming language: Fortran 77 and C

* Any restriction to use by non-academics: license needed

**Additional material**

**Additional file 1: DeStripe: supplementary figures.** The file contains two figures.

List of abbreviations

AFM: atomic force microscope; SEM: scanning electron microscope; SPM: scanning probe microscope; FFT: fast Fourier transform; CVAR: constrained variance; GTP: guanosine tri-phosphate.

Acknowledgements

The authors would like to thank M. Odorico and Prof. P. Parot (CEA Marcoule) for providing AFM images, Dr. J. Armengaud for GTPases molecules, and R. Podor (ICS Marcoule) for providing the SEM image. They also thank A. Martel (CEA GIPSI Saclay) for web site developments. The authors thank Dr. E. Darrouzet for editing and improving the manuscript. This work was partially supported by the French ANR agency (Grant ANR-07- PCVI-0002-01) and by the French Atomic and alternative Energies Commission.

Authors' contributions

SWWC participated in the design of the study, carried out the software
development, and drafted the manuscript. JLP participated in the design
and coordinated the study. Both authors read and approved the final
manuscript.

Received: 4 October 2010 Accepted: 1 February 2011

Published: 1 February 2011

References

1. Binnig G, Quate CF, Gerber C: *Atomic force microscope*. *Phys Rev Lett* 1986, 56:930-933.

2. Fechner P, Boudier T, Mangenot S, Jaroslawski S, Sturgis JN, Scheuring S: Structural information, resolution, and noise in high-resolution atomic force microscopy topographs. *Biophys J* 2009, 96:3822-3831.

3. Hörber JKH, Miles MJ: Scanning probe evolution in biology. *Science* 2003, 302:1002-1005.

4. Parot P, Dufrene YF, Hinterdorfer P, Le Grimellec C, Navajas D, Pellequer JL, Scheuring S: *Past, present and future of atomic force microscopy in life sciences and medicine*. *J Mol Recogn* 2007, 20:418-431.

5. Ando T, Uchihashi T, Kodera N, Yamamoto D, Taniguchi M, Miyagi A, Yamashita H: *High-speed atomic force microscopy for observing dynamic biomolecular processes*. *J Mol Recognit* 2007, 20:448-458.

6. Engel A, Gaub HE: Structure and mechanics of membrane proteins. *Annu Rev Biochem* 2008, 77:127-148.

7. Oesterhelt F, Scheuring S: *High-resolution imaging and force measurement of individual membrane proteins by AFM*. *Curr Nanosci* 2006, 2:329-335.
8. Hoh JH, Lal R, John SA, Revel JP, Arnsdorf MF: Atomic force microscopy and dissection of gap junctions. *Science* 1991, **253**:1405-1408.

9. Muller DJ, Schabert FA, Buldt G, Engel A: Imaging purple membranes in aqueous solutions at sub-nanometer resolution by atomic force microscopy. *Biophys J* 1995, **68**:1681-1686.

10. Scheuring S, Seguin J, Marco S, Levy D, Robert B, Rigaud JL: Nanodissection and high-resolution imaging of the Rhodopseudomonas viridis photosynthetic core complex in native membranes by AFM. *Atomic force microscopy. Proc Natl Acad Sci USA* 2003, **100**:1690-1693.

11. Odorico M, Teulon JM, Bellanger L, Vidaud C, Bessou T, Chen S-wW, Quemeneur E, Parot P, Pellequer JL: Energy landscape of chelated uranyl - antibody interactions by Dynamic Force Spectroscopy. *Biophys J* 2007, **93**:645-654.

12. Robert P, Benoliel AM, Pierres A, Bongrand P: What is the biological relevance of the specific bond properties revealed by single-molecule studies? *J Mol Recognit* 2007, **20**:432-447.

13. Neuman KC, Nagy A: Single-molecule force spectroscopy: optical tweezers, magnetic tweezers and atomic force microscopy. *Nat Meth* 2008, **5**:491-505.

14. Teulon JM, Parot P, Odorico M, Pellequer JL: Deciphering the energy landscape of the interaction uranyl-DCP with antibodies using dynamic force spectroscopy. *Biophys J* 2008, **95**:L63-L65.

15. Gaboriaud F, Parcha BS, Gee ML, Holden JA, Strugnell RA: Spatially resolved force spectroscopy of bacterial surfaces using force-volume imaging. *Colloids Surf B Biointer faces* 2008, **62**:206-213.

16. Walter NG, Huang CY, Manzo AJ, Sobhy MA: Do-it-yourself guide: how to use the modern single-molecule toolkit. *Nat Meth* 2008, **5**:475-489.

17. Müller DJ, Helenius J, Alsteens D, Dufrêne YF: Force probing surfaces of living cells to molecular resolution. *Nat Chem Biol* 2009, **5**:383-390.

18. Zhong Q, Inniss D, Kjoller K, Elings V: Fractured polymer/silica fiber surface studied by tapping mode atomic force microscopy. *Surf Sci Lett* 1993, **290**:L688-L692.

19. Gravel P, Beaudoin G, De Guise JA: A method for modeling noise in medical images. *IEEE Trans Medical Imaging* 2004, **22**:1221-1232.

20. Gonzalez RC, Woods RE: Digital image processing. New Jersey: Pearson Education; 3 2008.

21. Munch B, Trtik P, Marone F, Stampanoni M: Stripe and ring artifact removal with combined wavelet-Fourier filtering. *Opt Express* 2009, **17**:8568-8591.

22. Aizenberg I, Butakoff C: Frequency domain median-like filter for periodic and quasi-periodic noise removal. *SPIE Proc Image Process* 2002, **4667**:181-191.

23. Chan RH, Ho C-W, Nikolova M: Salt-and-pepper noise removal by median-type noise detectors and detailpreserving regularization. *IEEE Trans Image Process* 2005, **14**(10):1479-1485.

24. Koziol JA: *Intuitive biostatistics* New York: Oxford University Press; 1995.

25. Pok G, Liu J-C, Nair AS: Selective removal of impulse noise based on homogeneity level information. *IEEE Trans Image Process* 2003, **12**:85-91.

26. Baierlein R: *Newtonian dynamics* McGraw-Hill; 1983.

27. Madsen K, Nielsen HB, Tingleff O: *Methods for non-linear least squares problems*. 2 edition. Technical University of Denmark; 2004.

28. Lampert CH, Wirjadi O: An optimal nonorthogonal separation of the anisotropic Gaussian convolution filter. *IEEE Trans Image Process* 2006, **15**:3501-3513.

29. Reimer L, Hawkes PW: Scanning electron microscopy: physics of image formation and microanalysis. Springer, 2 1998.

30. Trinh M-H, Odorico M, Bellanger L, Jacquemond M, Parot P, Pellequer J-L: Tobacco mosaic virus as an AFM tip calibrator. *J Mol Recognit* 2011, **24**.

31. Kienberger F, Pastushenko VP, Kada G, Puntheeranurak T, Chtcheglova L, Riethmueller C, Rankl C, Ebner A, Hinterdorfer P: Improving the contrast of topographical AFM images by a simple averaging filter. *Ultramicroscopy* 2006, **106**:822-828.

32. Chi C-y, Zhang J-x, Liu Z-j: Study on methods on noise reduction in a stripped image. *XXI ISPRS Congress, Youth Forum: 2008; Beijing The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences*; 2008, 213-216.

33. Gwyddion - Free SPM (AFM, SNOM/NSOM, STM, MFM, ...) data analysis software. [http://gwyddion.net](http://gwyddion.net).

doi:10.1186/1472-6807-11-7

Cite this article as: Chen and Pellequer: DeStripe: frequency-based algorithm for removing stripe noises from AFM images. *BMC Structural Biology* 2011 11:7.

Submit your next manuscript to BioMed Central
and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central