PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
 
PLoS One. 2016; 11(6): e0157595.
Published online 2016 June 15. doi:  10.1371/journal.pone.0157595
PMCID: PMC4909204

A Computationally Efficient Filter for Reducing Shot Noise in Low S/N Data

Thomas Abraham, Editor

Abstract

Functional multineuron calcium imaging (fMCI) provides a useful experimental platform to simultaneously capture the spatiotemporal patterns of neuronal activity from a large cell population in situ. However, fMCI often suffers from low signal-to-noise ratios (S/N). The main factor that causes the low S/N is shot noise that arises from photon detectors. Here, we propose a new denoising procedure, termed the Okada filter, which is designed to reduce shot noise under low S/N conditions, such as fMCI. The core idea of the Okada filter is to replace the fluorescence intensity value of a given frame time with the average of two values at the preceding and following frames unless the focused value is the median among these three values. This process is iterated serially throughout a time-series vector. In fMCI data of hippocampal neurons, the Okada filter rapidly reduces background noise and significantly improves the S/N. The Okada filter is also applicable for reducing shot noise in electrophysiological data and photographs. Finally, the Okada filter can be described using a single continuous differentiable equation based on the logistic function and is thus mathematically tractable.

Introduction

Functional multineuron calcium imaging (fMCI) is a promising optical technique that simultaneously records suprathreshold activity from a large number of neurons [13]. An action potential evokes a transient elevation in intracellular calcium ions in the cell body of neurons. By taking advantage of this phenomenon, fMCI infers the times of action potentials from fluctuations in the fluorescence intensities in neurons that were loaded with calcium-sensitive fluorescent indicators [4]. Typically, an action potential elicits a rapid fluorescence increase that occurs within tens of millisecond, followed by a slow decrease that persists for up to a second [5]. The amplitude of the calcium transients is as small as 2–10% over the baseline intensity when Oregon green 488 BAPTA-1AM (OGB1-AM), the most widely used chemical calcium indicator, is used. Therefore, fMCI is vulnerable to signal-irrelevant noise and often suffers from low signal-to-noise ratios (S/N) [3].

Shot noise is dominant in under low-light intensity conditions during fMCI data acquisition. A major source of shot noise in fMCI is the photodetectors. fMCI mainly uses photomultiplier tubes, charge-coupled devices, or complementary metal-oxide-semiconductors to detect photons. Their signal is inevitably contaminated with internal noise due to the quantal nature of light and the probabilistic nature of photon detection [6]. Thus, the number of photons detected in a given time window could be larger and smaller than the true number, and these stochastic fluctuations cause shot noise.

In fMCI time-series data, shot noise is usually reduced through post hoc filtering using a moving window. The filters are classified into two types: linear and nonlinear filters. Most of the linear filters are moving-average low-pass filters with linearly weighted windows, such as rectangular, Gaussian, and Hann windows [79]. The basic idea of the linear filters is to replace the value at a focused time with the averaged value over the window. Unfortunately, the linear filters smooth sharply changing values, such as the rising phase of calcium transients, and hinder the precise detection of the onset times of the calcium transients. Moreover, a single outlier due to shot noise affects nearby values within the window. This undesirable effect is relevant in fMCI data that are acquired at sampling rates as low as tens of hertz.

Nonlinear filters may overcome these problems [10]. A widely used nonlinear filter is the median filter [11]. A median filter replaces the value at a focused time with the median value within a moving window. This procedure tends to preserve the structure of a signal; however, unlike linear filters, a median filter selects the replaced values from a pool of values that already existed in the original data, thereby resulting in less efficient denoising.

In this study, we propose a novel nonlinear filter named the Okada filter, which employs the compromise strategies of linear and median filters. The Okada filter replaces the value at a focused time with the average of two values at both ends of a moving window unless the focused value is the median among values in the window. Therefore, the filtered data contain "more centered" values that did not exist in the original values (as in linear filter) and have a greater likelihood of preserving the structure of the putative signal because any values during a consecutive change are unlikely replaced (as occurs when using a median filter).

Materials and Methods

Animal ethics

Experiments were performed with the approval of the animal experiment ethics committee at the University of Tokyo (approval number: P24-6) and according to the University of Tokyo guidelines for the care and use of laboratory animals.

Functional multineuron calcium imaging

Hippocampal organotypic slices were prepared from Wistar/ST rats (SLC), as described elsewhere [12]. Rat pups were anaesthetized via hypothermia and isoflurane and then decapitated. The brains were removed and placed in ice-cold oxygenated Gey’s balanced salt solution supplemented with 25 mM glucose. The brains were sliced horizontally (300 μm thick) using a vibratome (DTK-1500, Dosaka). Then, entorhinal-hippocampal regions were trimmed using a surgical microknife. The slices were placed on Omnipore membrane filters (JHWP02500, Millipore) and incubated in 5% CO2 at 35°C. The culture medium, which was composed of 50% minimal essential medium (Invitrogen), 25% Hanks’ balanced salt solution, 25% horse serum (Gibco), and antibiotics, was changed every 3.5 days. Experiments were performed at 11−19 days in vitro. Slices were transferred to a 35-mm dish that was filled with 2 ml of dye solution and incubated for 40 min in a humidified incubator at 37°C in 5% CO2 with 0.0005% OGB1-AM (Invitrogen), 0.01% Pluronic F-127 (Invitrogen), and 0.005% Cremophor EL (Sigma-Aldrich) [13]. The slices were then recovered for 30 min in oxygenated artificial cerebrospinal fluid (aCSF) that consisted of (in mM) 127 NaCl, 26 NaHCO3, 3.3 KCl, 1.24 KH2PO4, 1.2 MgSO4, 1.2 CaCl2, and 10 glucose, which were then bubbled with 95% O2 and 5% CO2. A slice was mounted in a recording chamber and perfused with aCSF at a rate of 1.5–2.0 ml/min for 15 min. The hippocampal CA3 pyramidal cell layer was imaged at 10 Hz using a Nipkow-disk confocal microscope (CSU-X1, Yokogawa Electric) equipped with a cooled CCD camera (iXonEM+DV897, Andor Technology) and an upright microscope with a water-immersion objective lens (16×, 0.8 numerical aperture, Nikon) [14]. Fluorophores were excited at 488 nm with a laser diode and visualized with a 507-nm long-pass emission filter. The fluorescence change was measured as (FtF0)/F0, where, Ft is the fluorescence intensity at time t, and F0 is the fluorescence intensity averaged from –10 to 10 s relative to t.

Electrical physiological recording

Recordings were performed in a submerged chamber perfused at 6–8 ml/min with oxygenated aCSF at 35°C. Whole-cell patch-clamp recordings were obtained from hippocampal CA3 pyramidal cells that were visually identified under an infrared-differential, interference-contrast microscope. For cell-attached recordings, glass pipettes (3–6 MΩ) were filled with aCSF [15]. For voltage clamp recordings, pipettes were filled with a cesium-based solution that consisted of (in mM) 130 CsMeSO4, 10 CsCl, 10 HEPES, 10 creatine phosphate, 4 Mg-ATP and 0.3 Na2-GTP [16]. Excitatory postsynaptic currents (EPSCs) were recorded at a clamped voltage of –70 mV. Signals were amplified and digitized at a sampling rate of 2 kHz using a MultiClamp 700B amplifier and a Digidata 1440A digitizer that were controlled by pCLAMP 10.4 software (Molecular Devices).

Okada filter

The Okada filter is given by the following equation:

xtxt+xt1+xt+12xt2(1+eα(xtxt1)(xtxt+1)),

where xt–1, xt, and xt+1 denote the ΔF/F values at frames t–1, t, and t+1, respectively. This equation is based on the logistic function. Therefore, when α is sufficiently large, it can perform two processes; (1) if (xtxt–1)(xtxt+1) ≤ 0, xt is substituted with xt itself, and (2) if (xtxt–1)(xtxt+1) > 0, xt is substituted with (xt–1 + xt+1)/2. This process is repeated subsequently from the second frame to the penultimate frame, whereas the values from the first frame and the last frame are not changed. All routines were written in MATLAB (The MathWorks).

Median filter

Median filter replaces xt with the median value within xt–1, xt and xt+1. This process is repeated subsequently from the second frame to the penultimate frame, whereas the values from the first frame and the last frame are not changed. Unlike the Okada filter, the median filter uses the original xt–1 before substitution as xt–1.

Binomial filter

Binomial filter replaces xt with a binomial weighted convolution of xt–1, xt and xt+1 (here, xt ← 0.25 xt–1 + 0.5 xt + 0.25 xt+1). This replacement is repeated subsequently from the second frame to the penultimate frame, whereas the values from the first frame and the last frame are not changed. Unlike the Okada filter, the binomial filter uses the original xt–1 before substitution as xt–1.

Savitzky-Golay filter

Savitzky-Golay filter replaces xt with the approximate value estimated by fitting a low order polynomial to a series of {xt–1, xt, xt+1} using the least squared method. This process is repeated subsequently from the second frame to the penultimate frame, whereas the values from the first frame and the last frame are not changed. Unlike the Okada filter, the Savitzky-Golay filter uses the original xt–1 before substitution as xt–1.

Photograph with shot noise

A portrait of Lenna was downloaded from the Standard Image Data-Base. Shot noise was added with a Poisson distribution as follows:

xi,j=1poissrnd(xi,j),

where xi,j represents the value of the pixel at row i and line j. This process was applied to all pixels of the photograph.

Peak signal to noise ratio

The peak signal-to-noise ratio (PSNR) is defined using the mean squared error (MSE) as follows:

PSNR=10logp2MSE,

where p is the maximum luminance value, i.e., 255 for an 8-bit image or 65,025 for a 16-bit image. When the m×n original image (I) was compared to its filtered image (K), MSE is given by the following equation:

MSE=1mni=1mj=1n(I(i,j)K(i,j))2.

Frequency response

The power spectrum is computed using the Fourier transformation function with a 10-s hamming window. The ratio of the power of the filtered trace to that of the original trace at a given frequency is defined as a frequency response.

Statistics

Summarized data are reported as the means ± standard deviations (SDs). For each cell, the logarithm of the mean amplitude of calcium transients divided by the SD of the fluorescence intensities during baseline periods, except ± 0.5 s relative to the peak of calcium transient, was calculated as the S/N. Wilcoxon signed rank test and paired t-test were performed to assess the significance of the differences. P<0.05 was considered statistically significant.

Results

Okada filter

The Okada filter is designed to remove outliers in a given time-series vector {xt} (t ranges from 1 to T). At each step at time t, the Okada filter compares xt to the immediately preceding and following values, xt–1 and xt+1, respectively. If xt is the median among xt–1, xt, and xt+1, that is, if (xt−xt–1)(xtxt+1) ≤0, then xt is not changed. If xt is not the median, that is, if (xtxt–1)(xt−xt+1) >0, then xt is substituted with the average of xt–1 and xt+1, i.e., (xt–1 + xt+1)/2 (Fig 1A). This procedure is conducted sequentially from t = 2 to T– 1 (Fig 1B). In each step, the value that was updated in the preceding step is used as xt–1. Notably, the substitution process can be expressed by a single equation using a logistic function (Fig 1C). When the coefficient α is sufficiently large, the equation takes either xt or (xt–1 + xt+1)/2, depending on the sign of (xt−xt–1)(xt−xt+1). Specifically, if (xt−xt–1)(xt−xt+1) is negative, the equation substitutes xt with the same value xt, otherwise the equation substitutes xt with (xt–1 + xt+1)/2. This conditional substitution is aimed to filter out a sudden outlier, such as shot noise, in the time-series vector {xt} without affecting a continuous increase or decrease in {xt} because the latter might reflect signal, such as the onset of a calcium transient. This conditional process is important to preserve the waveform of signal; note that wildly used filters, such binomial and Savitzky-Golay filters, inevitably smooth the waveforms of sharp signal (e.g., calcium transients in fMCI data) and may lead to erroneous detections of the exact timings of signal onsets (Fig 2). Nonlinear filters, such as the Okada filter and median filter preserve the true signal onsets (Fig 2).

Fig 1
Introduction of the Okada filter.
Fig 2
The Okada filter does not burr the signal onsets.

Filtering fMCI data

We applied the Okada filter to fMCI data. The fluorescence intensities were measured from the cell bodies from CA3 neurons in organotypic hippocampal slice cultures loaded with OGB1-AM. Because the amplitudes of calcium transients were maximally ~10%, the constant α was set to 100. The Okada filter reduced noise from the original calcium trace, embossing the waveforms of spike-evoked calcium transients (Fig 1D). Other examples are shown in Fig 3A and 3B, in which the calcium traces of 20 representative neurons were filtered using the Okada filter.

Fig 3
Application of the Okada filter to fMCI data.

We sought to quantify the denoising power of the Okada filter. We randomly sampled 100 cells from our laboratory data storage and applied the Okada filter. Then, the levels of signal and noise and their ratios (S/N) were compared before and after the filtering (Fig 3C–3E). The mean amplitudes of calcium transients (signal) were significantly reduced in the Okada-filtered traces compared to the original traces (Fig 3C, P = 3.9×10−18, Z = 8.68, Wilcoxon signed rank test, n = 100 cells). This undesirable signal reduction resulted from substitutions of the peak values xt of calcium transients with (xt–1 + xt+1)/2. Noise was quantified as the SD of the fluctuations in fluorescence during the baseline periods without calcium transients. Noise was significantly reduced in Okada-filtered trace (Fig 3D, P = 4.5×10−18, Z = 8.67). We calculated the S/N by dividing the mean amplitudes of calcium transients (signal) by the SD of the baseline fluorescence fluctuations (noise). The S/N of the Okada-filtered traces was significantly higher than that of the original traces (Fig 3E, P = 9.6×10−18, Z = 8.58), indicating that the Okada filter reduces noise to a greater degree than the signal.

Comparison of denoising power

To estimate the merit of the Okada filter, we repeated the same analysis for the median, binomial and Savitzky-Golay filters (Figs (Figs4A,4A, ,5A5A and and6A;6A; n = 100 cells). To facilitate comparisons, the identical datasets were filtered. All these filters also reduced the signal and noise (Fig 4B, P = 1.2×10−17, Z = 8.55; Fig 4C, P = 3.9×10−18, Z = 8.68; Fig 5B, P = 4.0×10−18, Z = 8.68; Fig 5C, P = 3.9×10−18, Z = 8.68; Fig 6B, P = 5.0×10−18, Z = 8.65; Fig 6C, P = 3.9×10−18, Z = 8.68) and increased the S/N (Fig 4D, P = 3.9×10−18, Z = 8.68; Fig 5D, P = 5.1×10−18, Z = 8.65; Fig 6D, P = 5.1×10−18, Z = 8.65).

Fig 4
Application of the median filter to fMCI data.
Fig 5
Application of the binomial filter to fMCI data.
Fig 6
Application of the Savitzky-Golay filter to fMCI data.

We compared the S/N of the Okada-filtered and median filtered traces (Fig 7A). The S/N varied among datasets; however, as a whole, the Okada filter produced significantly higher S/N than did the median filter (Fig 7A, P = 3.2×10−9, Z = 5.92). This difference appeared to be due to the greater ability of the Okada filter to remove baseline noise. Therefore, to examine whether the noise level affected the performance of these filters, we plotted the ratios of the S/N of Okada-filtered traces to those of median-filtered traces against the S/N of the original traces. These parameters exhibited a negative correlation (Fig 7B, Pearson’s r = −0.62, P = 3.6×10−71, t98 = 50), and the regression line crossed y = 1 at the S/N of 1.19 dB. Thus, the Okada filter removed noise more efficiently from the lower S/N data. We also compared the S/N of the Okada-filtered traces with the binomial-filtered and Savitzky-Golay-filtered traces in the same manner. The Okada filter improved the S/N to a degree that was comparable to the binomial filter (Fig 7C, t = 0.16, P = 0.87, paired t-test) but that was lower than the Savitzky-Golay filter (Fig 7E, t = 4.2, P = 5.7×10−5). The high performances of binomial and Savitzky-Golay filters were presumably because these filters preserve the signal amplitude more than the Okada filter. However, like the Fig 7B plot, the Fig 7D and 7F plots also exhibited a negative correlation (Fig 7D, Pearson’s r = −0.47, P = 1.4×10−69, t98 = 48; Fig 7F, Pearson’s r = −0.39, P = 9.1×10−69, t98 = 47). These results confirm that the Okada filter worked relatively efficiently in the lower S/N data.

Fig 7
Comparison of S/N improved by the Okada and the median, binomial, and Savitzky-Golay filters.

We computed the power spectra of the original traces and the traces that were filtered by the Okada, median, binomial and Savitzky-Golay filters (Fig 8A–8D, n = 100). Then we plotted the frequency responses, the ratios of the powers of the filtered traces to that of the original traces (Fig 8E). The Okada and median filters broadly reduced the power with a low-pass trend, whereas the binominal and Savitzky-Golay filters reduced the power in a more frequency-specific manner.

Fig 8
Comparisons of the frequency responses of the Okada, median, binomial, and Savitzky-Golay filters.

Comparison of computational speed

We next measured the computing speed of the Okada filter. We used a commercially available laptop computer equipped with Intel® Core™ i5-3337U Processor (64 bit, 1.80 GHz) and 8.0 GB RAM. Using the MATLAB code based the 'if/else' statement in Fig 9i), we filtered 100 calcium traces used above. These data were taken at 10 frames per second for 300 s and thus contained 3,000 data points. The computation time spent to filter a single trace was, on average, 247 ± 58 μs (mean ± SD of 100 cells). The Okada filter can be expressed using the logistic function (Fig 9ii). That the commands executed in Fig 9i) and 9ii) are theoretically identical, but the computation time of Fig 9ii) was 499 ± 51 μs and was twice longer than that of Fig 9i) (P = 2.1×10−17, Z = 8.49, Wilcoxon signed rank test, n = 100). This difference is probably because the exponential function demanded a higher computational cost. We also measured the computational speed of the median filter. Fig 9iii) uses the 'if' statement to examine whether xt is the median. Its computational time was 3.76 ± 0.85 ms and was approximately 14 times longer than that of Fig 9i) (P = 3.9×10−18, Z = 8.68). If we used a MATLAB function 'median', the computational time increased further to 23.7 ± 1.6 ms. Moreover, we examined the computational speed of the binomial and Savitzky-Golay filters using the MATLAB codes described in Fig 9v) and vi). Their computational times were 26.0 ± 3.4 ms and 1.17 ± 0.15 ms, respectively, and were both significantly longer than the Okada filter (both P = 3.9×10−18, Z = 8.68). Thus, the Okada filter is superior to the other filters in terms of computational speed. The rapid computation of the Okada filter is mainly due to the fact that it can skips the substituting process when the focused value xt is the median; note that the linear filter has to update all values in the vector.

Fig 9
Comparison of computational speed of the Okada and median, binomial, and Savitzky-Golay filters.

Effect of the window length of Okada filter on denoising performance

In most of the denoising filters, the window length is modifiable. Therefore, we compared the performance of the Okada filters with different window lengths. To preserve the amplitude of signal, we changed the substitutional procedure of the Okada filter. Specifically, the focused value xt is not changed if it is the median among the values within the window, whereas if xt is not the median, xt is substituted with the average of three points around the median, i.e., the median and two immediately higher and lower values (two nearest neighbor values) of the median within the window. We applied the Okada filters with different window lengths, three, five and seven frames, to fMCI data. In all cases, the filters reduced the noise of traces (Fig 10, left) and improved the S/N (Fig 10, middle, window = 3, P = 9.6×10−18, Z = 8.58, window = 5, P = 6.3×10−17, Z = 8.36, window = 7, P = 4.7×10−16, Z = 8.12, Wilcoxon signed rank test, n = 100). Compared to the window of three frames (the original Okada filter), the filter with longer window improved the S/N more effectively (Fig 10, right, window = 5, P = 1.8×10−10, Z = 6.38, window = 7, P = 1.3×10−8, Z = 5.68). Thus, the Okada filter with longer window improved the S/N of fMCI data more efficiently.

Fig 10
Application of the Okada filter with different window lengths to fMCI data.

Modified Okada filters

The mean filter can be modified to the binomial and Gaussian filters using appropriate weighting factors. Likewise, the Okada filter can also be modified by changing the coefficient β in the equation in Fig 11A; higher β conditions preserve xt more heavily in the replaced value. That is, the modified Okada filter is expected to work more similarly to a binomial filter when the value replacement occurs. We thus tested the performance of the modified Okada filter with different values of β > 2 (Fig 11B, left). Under all β conditions tested, the modified filter improved the S/N of fMCI data (Fig 11B, middle, β = 2, P = 1.3×10−17, Z = 8.54; β = 3, P = 8.8×10−18, Z = 8.59; β = 4, P = 8.3×10−18, Z = 8.60; β = 5, P = 9.3×10−18, Z = 8.58; β = 6, P = 6.1×10−18, Z = 8.63; β = 7, P = 5.3×10−18, Z = 8.63; β = 8, P = 6.1×10−18, Z = 8.63; β = 9, P = 6.1×10−18, Z = 8.63; β = 10, P = 6.1×10−18, Z = 8.63, Wilcoxon signed rank test, n = 100). Compared to the original Okada filter (β = 2), the effect of the modified Okada filters with higher β values were less sufficient to the low S/N data (Fig 11B, right, β = 3, Pearson’s r = 0.41, P = 2.4×10−76, t98 = 56; β = 4, r = 0.22, P = 1.8×10−65, t98 = 43; β = 5, r = 0.50, P = 5.5×10−57, t98 = 35; β = 6, r = 0.51, P = 2.0×10−52, t98 = 31; β = 7, r = 0.50, P = 1.6×10−47, t98 = 27; β = 8, r = 0.52, P = 4.0×10−46, t98 = 26; β = 9, r = 0.52, P = 3.6×10−44, t98 = 25; β = 10, r = 0.52, P = 1.2×10−42, t98 = 24). As a result, the S/N of the filtered traces were the highest at β = 2 (the Okada filter) and were decreased with larger β (Fig 11C).

Fig 11
Application of the modified Okada filters to fMCI data.

Further applications of the Okada filter

Although we originally developed the Okada filter to denoise vector sequences with low S/N and with low temporal resolution, such as fMCI data, we finally tried to apply it to other forms of data.

We first filtered electrophysiological data acquired at a high temporal resolution (Fig 12). Excitatory postsynaptic currents were recorded at 2 kHz from CA3 pyramidal cells in a hippocampal slice culture. The Okada filter and the median filter both reduced the baseline noise; however, there was still noise in the filtered traces. We thus repetitively applied the filters and found a critical difference in the denoising performance between the Okada filter and the median filter, particularly in later stages of the iterations. Noise was gradually reduced when the Okada filter was repeated, but the median filter did not show this effect. The Okada filter fundamentally substitutes xt with a new value that is different from xt–1, xt, or xt+1. Conversely, the median filter necessarily selects a substituted value from either xt–1, xt, or xt+1, and thus, the effect of iterative filtering is limited.

Fig 12
Application of the Okada and median filters to electrophysiological data.

We also investigated the effect of the Okada filter on photo images with shot noise. The Okada filter compares values of three neighboring points in a one-dimensional vector {xt}, and this idea cannot be applied to two-dimensional (2D) matrix {xi,j}. Therefore, we made a small modification. The 2D-modified Okada filter substitutes xi,j with the average of the vertically and horizontally scanned Okada-filtered values. This process is expressed by the following equation:

xi,jxi,j+[xi1,j+xi+1,j2xi,j2{1+eα(xi,jxi1,j)(xi,jxi+1,j)}+xi,j1+xi,j+12xi,j2{1+eα(xi,jxi,j1)(xi,jxi,j+1)}]/2.

One representative video frame of fMCI is shown in Fig 13A, which contains shot noise because it was taken under low-light intensity conditions. Because the data were presented in the form of a 16-bit digital image, α was set to 106. When the Okada filter and the median filter were applied to this image, shot noise was reduced and the contours of individual neurons could be more easily identified (Fig 13A). We also filtered a natural photograph (portrait) in which the shot noise was artificially added. The Okada filter and the median filter both successfully smoothed out shot noise (Fig 13C). In both images, we calculated the PSNR to evaluate the denoising effect of the Okada and median filters. The PSNR of the fMCI image in Fig 13A and the photograph in Fig 13C were both higher than the original image (Fig 13B, original versus Okada, original versus median, P = 1.8×10−5, Z = 4.29, Wilcoxon signed rank test, F23,23 = 2.13, 2.20, F-test of equality of variances, n = 24 subregions; Fig 13D, P = 1.03×10−20 and 2.44×10−21, t15 = 32.5 and 34.6, paired t-test, F15,15 = 1.43 and 1.08, respectively, n = 16 compartments.). In the fMCI image, the Okada filter was slightly inferior to the medial filter in terms of the PSNR evaluation (P = 4.5×10−12, t23 = 13.0, paired t-test, F23,23 = 1.03).

Fig 13
Application of the Okada and median filters to photographs.

Discussion

In this study, we devised an Okada filter, a novel nonlinear filter that is optimized to data with low S/N and with low sampling rates. Indeed, the Okada filter reduces noise in fMCI data and improves their S/N. We compared the performance of the Okada filter to that of the median filter, binomial filter and Savitzky-Golay filter. The Okada filter surpassed the median filter in denoising and computational speed. Moreover, the Okada filter showed a higher computational speed than the binomial and Savitzky-Golay filters. Therefore, the Okada filter is more suitable for practical use. The frequency response analysis revealed that the Okada filter broadly reduces high frequency components, suggesting that the Okada filter works more effectively for white noise-like background fluctuations, such as thermal noise. In addition, the modified Okada filter also improved the S/N of the fMCI data. One may be able to adjust the parameter β to a given purpose.

The S/N in fMCI has been improved by development of more sensitive calcium indicators and more optimized optical systems, such as objective lens and photodetectors [17]; however, their shot noise is stochastically emitted during the process of photon conversion and can come from cosmic rays. Thus, online removal of the shot noise during image acquisition is technically difficult. Therefore, fMCI data are usually denoised using offline filters that are employed in image processing and sound processing [10,18,19]. Each filter has advantages and disadvantages. Linear filters are usually effective in reducing high-frequency noise, including shot noise; however, because of their nature, they often blur the edges of the signal and thereby deform the shapes of the signal. Conversely, nonlinear filters have less influence on the edges, but their denoising powers may be insufficient. In this way, the Okada filter reduces noise by averaging values within a window, such as the linear filter, and leaves the putative signal, such as the median filter, thereby benefitting from the advantages of both types of filters. Furthermore, we found that the Okada filter with longer windows improved the S/N of fMCI traces more efficiently. Because the amplitudes of signal tended to be more severely reduced when using longer windows (relatively to the frame rates), the Okada filters with long windows may be suitable for fMCI data recorded at higher temporal resolution. In the Okada filter, the parameter β is modifiable to optimize its denoising power depending on the S/N of the original data. Some of the recently developed calcium indicators are sensitive enough to exhibit the signal amplitude larger than 50% over the baseline intensity. In such cases with excellent S/N, higher β values might be suitable.

Another disadvantage of nonlinear filters is, in general, the high cost for von Neumann computation; however, the Okada filter is a nonlinear filter that exhibits more rapid computation than the median filter. Its computational speed is almost a match for that of simple linear filters. This efficient computation may be applicable to an immediate filter that works online during data sampling. Moreover, we packed the data processing of the Okada filter into a single mathematical formula, while the median filter cannot be formulated in a single algebraic equation. The equation of Okada filter is a continuous, integrable, and differentiable function. Thus, it is not only arithmetically tractable but may be also applicable to development of electric circuit-based denoising devices.

Acknowledgments

We thank Mr. Nariaki Uemura for his initial trials of the Okada filter and Mr. Tetsuya Sakaguchi for optimizing MATLAB codes.

Funding Statement

This work was supported by Grants-in-Aid for Science Research on Innovative Areas, "Mesoscopic Neurocircuitry" (22115003).

Data Availability

Data Availability

All relevant data are within the paper.

References

1. Takahashi N, Sasaki T, Usami A, Matsuki N, Ikegaya Y (2007) Watching neuronal circuit dynamics through functional multineuron calcium imaging (fMCI). Neurosci Res 58: 219–225. [PubMed]
2. Namiki S, Ikegaya Y (2009) Current application and technology of functional multineuron calcium imaging. Biol Pharm Bull 32: 1–9. [PubMed]
3. Takahashi N, Oba S, Yukinawa N, Ujita S, Mizunuma M, Matsuki N, et al. (2011) High-speed multineuron calcium imaging using Nipkow-type confocal microscopy. Curr Protoc Neurosci 2: 14 doi: 10.1002/0471142301.ns0214s57 [PubMed]
4. Smetters D, Majewska A, Yuste R (1999) Detecting action potentials in neuronal populations with calcium imaging. Methods 18: 215–221. [PubMed]
5. Sasaki T, Takahashi N, Matsuki N, Ikegaya Y (2008) Fast and accurate detection of action potentials from somatic calcium fluctuations. J Neurophysiol 100: 1668–1676. doi: 10.1152/jn.00084.2008 [PubMed]
6. Sjulson L, Miesenbock G (2007) Optical recording of action potentials and other discrete physiological events: a perspective from signal detection theory. Physiology (Bethesda) 22: 47–55. [PubMed]
7. Savitzky A, Golay MJE (1964) Smoothing and differentiation of data by simplified least squares procedures. Anal Chem 36: 1627–1639.
8. Little MA, Jones NS (2013) Signal processing for molecular and cellular biological physics: an emerging field. Philos Trans A Math Phys Eng Sci 371: 20110546 doi: 10.1098/rsta.2011.0546 [PMC free article] [PubMed]
9. Ikegaya Y, Aaron G, Cossart R, Aronov D, Lampl I, Ferster D, et al. (2004) Synfire chains and cortical songs: temporal modules of cortical activity. Science 304: 559–564. [PubMed]
10. Yaksi E, Friedrich RW (2006) Reconstruction of firing rate changes across neuronal populations by temporally deconvolved Ca2+ imaging. Nat Methods 3: 377–383. [PubMed]
11. Arias-Castro E, Donoho DL (2009) Does median filtering truly preserve edges better than linear filtering? Annal Stat 37: 1172–1206.
12. Koyama R, Muramatsu R, Sasaki T, Kimura R, Ueyama C, Tamura M, et al. (2007) A low-cost method for brain slice cultures. J Pharmacol Sci 104: 191–194. [PubMed]
13. Ishikawa T, Okamoto K, Ikegaya Y (2015) Homeostatic hippocampal activity against reduced glutamatergic neurotransmission. Int J Pharmacol 11: 318–326.
14. Mizunuma M, Norimoto H, Tao K, Egawa T, Hanaoka K, Sakaguchi T, et al. (2014) Unbalanced excitability underlies offline reactivation of behaviorally activated neurons. Nat Neurosci 17: 503–505. doi: 10.1038/nn.3674 [PubMed]
15. Sasaki T, Matsuki N, Ikegaya Y (2011) Action-potential modulation during axonal conduction. Science 331: 599–601. doi: 10.1126/science.1197598 [PubMed]
16. Takahashi N, Kitamura K, Matsuo N, Mayford M, Kano M, Matsuki N, et al. (2012) Locally synchronized synaptic inputs. Science 335: 353–356. doi: 10.1126/science.1210362 [PubMed]
17. Grienberger C, Konnerth A (2012) Imaging calcium in neurons. Neuron 73: 862–885. doi: 10.1016/j.neuron.2012.02.011 [PubMed]
18. Fiore L, Corsini G, Geppetti L (1996) Application of non-linear filters based on the median filter to experimental and simulated multiunit neural recordings. J Neurosci Methods 70: 177–184. [PubMed]
19. Szymanska AF, Kobayashi C, Norimoto H, Ishikawa T, Ikegaya Y, Nenadic Z Accurate detection of low signal-to-noise ratio neuronal calcium transient waves using a matched filter. J Neurosci Meth: in press. [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science