Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Eddy Current Sensor Array for Electromagnetic Sensing and Crack Reconstruction with High Lift-Off in Railway Tracks
Previous Article in Journal
When Trustworthiness Meets Face: Facial Design for Social Robots
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Efficacy of Nonlinear Filters in CMOS for 2-D Signal Processing for Image Quality Enhancement

by
Hector Bandala-Hernandez
1,
Alejandro Bautista-Castillo
1,*,
José Miguel Rocha-Pérez
1,2,
Victor Hugo Carbajal Gómez
1 and
Alejandro Díaz-Sánchez
1
1
National Institute for Astrophysics, Optics and Electrónics, Luis Enrique Erro #1, Sta María Tonanzintla, Puebla 72840, Mexico
2
Benemérita Universidad Autónoma de Puebla, Facultad de Ciencias de la Electronica, Ciudad Universitaria, Mexico City 72000, Mexico
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(13), 4213; https://doi.org/10.3390/s24134213
Submission received: 25 April 2024 / Revised: 1 June 2024 / Accepted: 18 June 2024 / Published: 28 June 2024
(This article belongs to the Section Sensing and Imaging)

Abstract

:
This study rigorously investigates the effectiveness of nonlinear filters in CMOS for 2-D signal processing to enhance image quality. We comprehensively compare traditional linear filters’ performance, which operate on the principle of linearity, with nonlinear filters, such as the median-median (Med-Med) approach, designed to handle nonlinear data. To ensure the validity of our findings, we use widely accepted metrics like normalized squared error (NSE), peak signal-to-noise ratio (PSNR), and structural similarity index (SSIM) to quantify the differences. Our simulations and experiments, conducted under controlled conditions, demonstrate that nonlinear filters in CMOS outperform linear filters in removing impulse noise and enhancing images. We also address the challenges of implementing these algorithms at the hardware level, focusing on power consumption and chip area optimization. Additionally, we propose a new architecture for the Med-Med filter and validate its functionality through experiments using a 9-pixel image sensor array. Our findings highlight the potential of nonlinear filters in CMOS for real-time image quality enhancement and their applicability in various real-world imaging applications. This research contributes to visual technology by combining theoretical insights with practical implementations, paving the way for more efficient and adaptable imaging systems.

1. Introduction

The convergence of various technologies, such as 2-D signal processing, image filters, CMOS processors, and image quality assessment, has sparked a vibrant discourse in visual technology research [1,2]. This intersection is a meeting point of technologies and a battleground of challenges and opportunities, from fundamental theoretical foundations to practical applications. Two-dimensional signal processing, for instance, equips us with valuable tools for eliminating unwanted components [3]. Within 2-D arrays, the strategic application of linear and nonlinear filters in the spatial domain becomes a dynamic solution to filter out high-frequency undesired elements [1,4,5]. However, the conflict presented by linear filters, which are applicable in specific situations but result in blurring and information loss, prompts the investigation of nonlinear alternatives. The theoretical foundations of two-dimensional CMOS processors, conceptualizing them as processing elements (PEs) intricately woven into kernels within 2-D sensor arrays, is another area of intense research. The primary focus is addressing nonlinear non-uniformities among PEs arising from systematic and random mismatches [6,7]. This scenario explores various operations, customizing each to address issues arising from the statistical distribution of non-uniformities in each PE’s response. In the expansive realm of visual technology, filters play a pivotal role in shaping and enhancing digital images. From subtle smoothing adjustments to advanced edge-detection techniques, these tools are the backbone of image processing. Each filter serves a unique function, including smoothing, eliminating imperfections, edge detection, and highlighting critical structures.
Noise introduces challenges within digital imagery that demand nuanced, thorough understanding [8]. Fixed pattern noise (FPN), or spatial non-uniformity, refers to a consistent and predictable deviation from expected pixel values [9]. It often emerges because of imperfections or non-uniformities in imaging sensors. As image quality becomes paramount, introducing figures of merit for image quality offers a systematic approach to evaluating fidelity, clarity, and perceptual accuracy [10,11]. Metrics like the structural similarity index (SSI), peak signal-to-noise ratio (PSNR), normalized squared error (NSE), mean squared error (MSE), perceptual image quality assessment (PIQA), and structural similarity map (SSIM map) serve as benchmarks for assessing various facets of image quality [12]. This exploration seeks to comprehend linear and nonlinear filters, the role of 2-D-CMOS pixels, and the importance of image quality assessment metrics [13]. By comprehensively understanding and critically evaluating these elements, we aim to empower diverse fields with tools for consistently delivering and appraising superior visual content.
This paper is structured as follows. Section 2 analyzes the balance between quality improvement and energy consumption for different filter implementations. Section 3 describes the design of a weighted median filter (WMF) circuit for these multidimensional filters. Section 4 compares the performance of the WMF circuit with previous median filters. Section 5 explains the implementation of the Med-Med filter using a 9-pixel image sensor array. Section 6, presents the characterization setup. This is followed by Section 7, which presents experimental results. Finally, Section 8 concludes.

2. High-Level Simulation Applying Multilevel 2-D Filters

Three crucial metrics for assessing 2-D data quality, namely NSE, PSNR, and SSIM [7,11], are employed. Nonlinear contamination is introduced to the sample to validate the hypothesis, which is then processed using a mean filter, median filter, Med-Med, and Mean-Med techniques. These filters are not randomly chosen but specifically selected to correspond to the multilevel approximations explored in this study, underscoring their significance in our research.
Simulations of the filters above were conducted on 509 × 383 pixels and 256 grayscale images, with the results depicted in Figure 1.
The original image is deliberately contaminated with 20 % impulse noise, a level carefully chosen to assess performance under various working conditions. This deliberate choice of noise level ensures that our research covers a wide range of scenarios. Figure 1 illustrates the outcomes of mean, median, median of median, and mean of median filters, respectively.
Table 1 compares specific performance metrics between processed images and their unaltered counterparts, focusing on NSE, PSNR, and SSIM. As depicted in Table 1, SSIM ideally reaches 1; however, the mean filter exhibits a notably low value, attributed to the blurring effect observed when applied to a 2-D array. Furthermore, human perception of image quality in terms of SSIM is notable; notably, the image produced by the mean filter displays a low resemblance to the original image, which is evident in its low SSIM. Conversely, employing the Med-Med filter yields the highest SSIM.
Moreover, normalized squared error (NSE) could significantly affect the perception of differences between the original and processed images. NSE necessitates the determination of the maximum pixel value. In our simulation, the variable M A X i , j is initialized to 256, representing the total number of grayscale levels within the image. The process depicted in Figure 2 illustrates the extraction of NSE, mirroring the procedure outlined in Figure 1.
Figure 2a demonstrates the results obtained by applying an ideal filter, resulting in an output image identical to the original, with a normalized squared error N S E = 0 . Following this, Figure 2b showcases a corrupted image with 20 % impulse noise, while Figure 2c–f illustrate the NSE values of the applied filters. Noticeably, the NSE of the Med-Med and Mean-Med filters is significantly lower than that of both the mean and median Filters. In certain instances, the scaling in Figure 2 is adjusted to ensure perceptibility to the human eye.
Figure 3 shows the MSE behavior with four filters applied to an image contaminated with impulsive noise, showing that the mean filter is ineffective at any noise level, maintaining a consistently high MSE of around 8965. In contrast, the median and Med-Med filters exhibit better performance at low-to-moderate noise levels. However, their effectiveness diminishes as the noise increases, with the MSE rising from 120.1 to 983.1 and 114.1 to 2013.6, respectively. The Mean-Med filter performs best at low noise levels, with an MSE starting at 67.7 and increasing to 1900.7 at very high noise levels. Median-based filters are more effective at handling impulsive noise, particularly at lower levels, but their efficacy decreases with higher noise percentages.

Power Constraints

Considerations must be taken into account when implementing a low-power integrated circuit (IC), especially regarding power consumption. Figure 4 illustrates a comparison of the estimated power consumption between Med-Med and Mean-Med filters. In previous research [14], a low-power median filter was assessed and serves as a reference benchmark. The comparison in Figure 4 is based on the number of cells within a 2-D array. It is assumed that each median filter within a multilevel filter setup consumes 90 nW of static power per cell. The maximum power consumption depicted in Figure 4 is computed for an array consisting of 32 × 32 processing cells.
The observed compromise between enhancing quality and conserving power suggests that while the Mean-Med scheme provides better quality enhancement than the Med-Med algorithm, implementing it on each cell increases power consumption as the number of cells grows. Consequently, when dealing with extensive 2-D arrays, opting for column-wise filtering over deploying a processing element on each cell is preferable to alleviate power usage. This method capitalizes on the inherent parallelism in column-wise operations, thereby reducing the overall power requirements.
On the other hand, the Med-Med algorithm emerges as a more fitting solution for achieving superior quality enhancement without sacrificing low-power characteristics. Its inherently localized nature enables efficient processing without substantially escalating power demands, rendering it an appealing choice for applications where power efficiency is crucial.

3. Hardware Implementations of the Algorithms

Initially, attention is directed towards the fundamental design of the Mean-Med and Med-Med filters, specifically the circuit responsible for extracting the median from a set of values [15,16,17]. This median extractor circuit is versatile and can operate as a weighted median extractor, enhancing its applicability across various scenarios.

Median Extractor and Weighted Median Filter

The weighted median filter (WMF) is defined by a function of the estimator β , which minimizes the weighted L 1 norm, as depicted in Equation (1) [18,19,20].
Σ i = 1 N w i x i β m i n i , w i N
In this context, N represents the number of inputs, x i represents the i-th entry of a continuous-valued input vector, and w i denotes the weight coefficient vector, where m i n represents the statistical minimum. Each element in the weight vector W = w 1 , w N must be a real, non-negative integer, and the sum of the components must be an odd number.
The depicted median circuit extractor in Figure 5 employs a transconductance comparator per datum. The proposed weighted median filter (WMF) topology can be realized by utilizing the tail current to establish different saturation levels and assign weights to each input. Each I b n represents a non-negative multiple of a base current I x = 10 nA. Thus, if I b 1 to I b 5 are 10 nA, the circuit operates as a standard median extractor. The sub-window schematic circuit of this non-weighted median extraction is:
For this circuit, operational transconductance amplifiers (OTAs) with input data exceeding v o u t will saturate positively, resulting in positive saturation currents + I o u t i flowing into node v o u t . Conversely, for OTAs with input data below V o u t , the output currents will yield negative output saturation currents I o u t j from node V o u t . Maintaining equilibrium, the number of weighted data below the median must be consistent, as illustrated in Equation (2).
Σ i I o u t i Σ j I o u t j 0
For all tail currents except that of the comparator whose input is the median of the set of inputs, the output voltage, V o u t , of the corresponding OTA closely approximates v i n . In any weighted median filter where the sum of the weights is an odd number, the output voltage tracks the corresponding input in the dataset. All comparators in the proposed filter operate in the subthreshold region and utilize I b n as the tail current when a non-weighted median is desired, resulting in a power consumption of 90 nW for the entire WMF cell. When a weighted median is required, the power consumption is calculated using Equation (3) [21,22,23].
P W M F = 18 n W Σ i n w i
Table 2 displays the assigned transistor sizes of the weighted median circuit in Figure 5, the bias current, and the supply voltage of the median circuit extractor.
The weight vector of the WMF can be represented as:
I b n = I x w n n = 1 , 2 5
where w n N and is the weight used for the n-th datum. Although the comparator transconductance could vary due to differences in the tail currents, the filter’s performance remains unaffected. A weighted median extraction example is depicted in Figure 6, where the mask W b = { 3 , 1 , 1 , 1 , 3 } is applied. Both median and weighted median extractors function as the core block for the implemented multilevel filters in this study. Subsequent sections detail the design of various multilevel filter implementations.

4. WMF Experimental and Simulation Results

The WMF circuit was fabricated using a 0.5 μ m technology provided by ON Semiconductors. Table 3 details the five inputs utilized for characterizing the chip. Figure 6b illustrates the experimental response of the median circuit, employing the weight mask W = { 1 , 1 , 1 , 1 , 1 } , with all comparators biased at a current of 10 nA. Discrepancies between the corners of the simulated output signal (pink line) and those of the experimental output (yellow line) are observed; the experimental results do not fully reach the corners corresponding to the ideal median due to a single comparator being assigned to each resulting median. However, in cases where the median aligns with such corners, the circuit needs to switch between the current and the subsequent comparator.
The five patterns of input window covers utilized during the weighting procedure are indicated as { w 1 , w 2 , w 3 , w 4 , w 5 } , where w n symbolizes the weight coefficient for the n-th input. The cover pattern W a = { 3 , 1 , 1 , 1 , 1 } was utilized with the WM filter. This assigns a superior weight to input V i n 1 when compared with the remaining inputs. Within this pattern, the digit 1 represents the standard value of 10 nA for the tail current of the respective comparator, while the digit 3 corresponds to a tail current of 30 nA for the weighted transconductance comparator. MATLAB simulation outcomes employing cover W a are illustrated in Figure 7a, where five waveforms (two sine waves, two triangular waves, and a square wave) are utilized for an ideal WMF; the yellow waveform portrays the weighted median outcome. Furthermore, Figure 7b demonstrates experimental outcomes of the identical cover application. In this evaluation, the yellow waveform represents the WMF outcome, where five 10 KHz/400 mVpp waveforms were employed for the simulation outcomes and 10 KHz/600 mVpp for the experimental outcomes, in alignment with the earlier scenario.
The current mode technique of the described circuit enable overcoming certain constraints of analog median filters reported previously [14,18,20]. This comparison is presented in Table 4, which assesses latency, silicon area, power consumption, supply voltage, operational frequency, technology, input count, comparator gain, accuracy, type (analog or digital), and input common mode range (ICMR).
Based on Table 4, this study necessitates minimal silicon area and significantly reduces power consumption compared with other research efforts, rendering it appropriate for ultra-low power applications. The proposed analog WMF underwent validation through simulations utilizing the BSIM3 Level 49 model for a one-dimensional signal and two weight masks. Discrepancies between simulation and experimental outcomes arise from constraints on the output signal’s corners, influenced by the switching between saturated comparators operating concurrently. To mitigate such discrepancies, increasing the transconductance of the comparators may be necessary.

5. 9-Pixels Image Sensors Array for Median of Medians

According to the high-level simulations performed in Figure 1 and Figure 2, the Med-Med filter demonstrates intermediate improvements in NSE, PSNR, and SSIM. Consequently, the hardware implementation of the Med-Med filter is shown in Figure 8 [24,25,26]. The filter comprises the 5-input median extractor circuit, as depicted in Figure 5. Therefore, the pixel array must be configured as a 5 × 5 grid and applied in the following manner. The first median filter processes pixels from the top-right corner to the bottom-left corner, the second median filter processes pixels from the top-left corner to the bottom-right corner, the third median filter processes the central column, and the fourth median filter processes the central row. The outputs of these four filters are connected to a fifth median filter, with the fifth input of the filter also processing the central pixel of the array. Consequently, this arrangement constitutes the Med-Med filter.
The VLSI implementation of this algorithm was fabricated using AMS 0.35 μ m Opto technology within a multi-project wafer. Due to design area limitations, the 5 × 5 array could not be implemented, and a 3 × 3 array was utilized instead. To maintain the functionality of the 5-input median extractor circuit, the I n 1 inputs were connected to V D D and the I n 5 inputs were connected to V S S . Each pixel consists of a 3 TAPS (3-transistor active pixel sensor). Since this array does not feature a correlated double sampling (CDS) system, a 4 TAPS (4-transistor active pixel sensor) was not implemented. The circuit employs a two-stage analog comparator, as illustrated in Figure 9 [27,28].
The low-frequency voltage gain of the comparator is show in Equation (5)
V o u t V i n = g m 1 g m 5 ( r o 4 r o 2 ) ( r o 6 r o 8 )
The initial phase comprises a differential pair employing NMOS transistors at its differential input and a PMOS transistor current mirror serving as an active load. The subsequent phase incorporates another differential pair utilizing PMOS transistors at its input and an NMOS transistor current mirror as an active load. The transistor type selection at each stage’s input aims to address the offset in the DC level voltage at the output. Equation (6) define the r o , n resistance.
r o , n = 1 λ I D
Subsequently, a tail current of 20 μ A is suggested for the initial stage, while 10 μ A is proposed for the second stage.
Table 5 shows the characteristics of the comparator.

6. Characterization Set-Up

Firstly, Figure 10 shows the layout of the 3 × 3 sensors array for the Med-Med filter accomplishment. Each section of the array is composed of a 5-input, subthreshold median filter used as sub-filter for the Med-Med filter, 3-transistors active pixel sensor, comparator, and a bias cell that supplies a current of 12.5 nA to each WMF. Figure 11 depicts a micro-photograph showcasing the 3 × 3 sensor array tailored for the Med-Med filter. The fabrication process leveraged AMI’s 0.35 μ m CMOS Opto Technology, ensuring precise and efficient functionality.
A meticulously designed and fabricated PCB (printed circuit board), as depicted in Figure 12, was utilized to complement this characterization of the Med-Med IC implementation. Concurrently, optic excitation parameters were configured to effectively apply mask patterns across the 3 × 3 array. For this purpose, a laser operating at 750 nm wavelength, coupled with an optic fiber, was employed to generate a spot size of 150 μ m (as illustrated in Figure 13).
The experimental arrangement, as depicted in Figure 13, comprises diverse elements. The left side features a micro-manipulator containing the optic fiber, alongside the PCB, as illustrated in Figure 13. This PCB connects to a Keysight 16851A logic analyzer, enabling the capture of associated output digital words. Conversely, the right side of the image provides a microscopic view, demonstrating the accurate positioning of the fiber onto the chip.

7. Experimental Results

The pixels of the 9-pixel array were illuminated using an optical fiber termination. The fiber was aligned manually, with the moving capability along the X and Y directions across the array. Due to alignment imperfections, the fiber was positioned as centrally as possible over the desired pixel, although some of the incident light reached adjacent pixels. Consequently, the pixel values were measured manually, one by one, to obtain the median sub-filter values, their outputs, and, ultimately, the output of the main median filter centered on the image below.

Results for M s k 2

Msk2 refers to the illumination pattern of the 3 × 3 pixel array centered at position (1,2), resulting in a HEX 35 value for this pixel, color red in Figure 14. Due to the lower intensity of illumination in the neighboring pixels, different values are obtained (Figure 14). In this instance, the output of the Med1 filter is HEX 90 (Figure 15).
M E D F 1 sub-filter has inputs: { 90 , C A , 90 } ; therefore, its output is ninety, as shown in Figure 15.
The M E D F 2 and M E D F 3 filters in M s k 2 have the same result as their respective sub-filter. M E D F 2 has inputs { 35 , C A , 3 E } ; therefore, its median is 3 E . Similarly, inputs for M E D F 3 are { 3 E , C A , 3 E } , resulting in the HEX 3 E value. This result is shown in Figure 16 highlighted in blue.
Inputs for M E D F 4 sub-filter applying M s k 2 are { C A , C A , 35 } ; hence, the corresponding output is C A , as shown in Figure 17.
Finally, inputs for the highest hierarchy median filter are { 90 , 3 E , 3 E , C A , 3 E } , resulting in a Med-Med of 3 E , shown in Figure 18.
Having accessed each of the values, the functionality of the Med-Med filter was verified. This confirms that the hardware-level implementation is performing real-time filtering with filters consuming 90 nW.
Table 6 provides a detailed comparison between our proposed methodology and three alternative approaches within the field. While each method presents unique features, our solution, developed using AMS 0.35 μ m technology, offers several significant advantages. Primarily, it showcases power efficiency, with power consumption as low as 450 nW during Med-Med operation when all five filters are active, and only 90 nW for a single median filter at 500 KHz with 1.8 Vpp, making it highly suitable for energy-constrained applications. Furthermore, our methodology has a dynamic range spanning from 0.1 V to 1.7 V, ensuring versatility in signal processing, coupled with an accuracy of 40 mV. Additionally, our design, featuring a WMF operating in the subthreshold regime, promises robust performance across varied scenarios. In contrast, the alternatives vary in their offerings; for instance, ref. [29] present a pixel array design with a single median filter, providing a different trade-off between power consumption and performance. Ref. [30] explore the use of comparators, while [17] implement a median filter, each with its strengths and limitations. However, our methodology sets itself apart by requiring no data sorting, simplifying the processing pipeline, and occupying a silicon area of 0.013 mm   2 , making it space-efficient and cost-effective. Notably, our implementation encompasses a holistic solution, incorporating a 9-pixel array with photosensors, readout components, comparators, median filters, and digital logic for reading, ensuring comprehensive signal processing in imaging applications.

8. Conclusions

This study comprehensively explores the intersection between 2-D signal processing and image quality assessment, focusing on the crucial role of linear and nonlinear filters alongside the design and implementation of two-dimensional CMOS processors. Through simulations, practical experiments, and theoretical analyses, we have showcased that these nonlinear filters, particularly the Med-Med approach, offer substantial improvements in NSE, PSNR, and SSIM compared with conventional linear filter-based methods.
Implementing these algorithms in integrated circuits poses additional challenges, including power management and chip area optimization, which are essential. While the Mean-Med approach may offer notable enhancements in image quality, its integration into low-power circuits could incur prohibitive energy consumption costs. Conversely, the Med-Med approach delivers significant image quality improvements without excessively compromising power efficiency.
Moreover, we introduced a novel Med-Med filter architecture realized in a 9-pixel image sensor array and confirmed its functionality through practical experiments. The obtained results affirm the efficacy and feasibility of this architecture in real-world applications and highlight its potential for integration into computer vision systems and imaging devices.
The conclusions drawn from our analysis show the positive impact of the Med-Med filter on image quality. With a notably low NSE of 115.2, the Med-Med filter significantly outperforms both mean and median filters regarding restoration accuracy, showing minimal discrepancy between the original and restored images, translating to superior visual fidelity.
Furthermore, the PSNR of 27.5 achieved with the Med-Med filter confirms its efficacy in preserving image quality, as higher values denote minimal information loss during restoration. Conversely, mean and median filters show inferior PSNR performance, suggesting a more substantial degradation of image quality.
The structural similarity index (SSIM) of 0.988 reached with the Med-Med filter further underscores its ability to preserve the structure and details of the original image, showing high similarity between the restored and reference pictures and confirming the effectiveness of the Med-Med filter in keeping crucial visual information.
In summary, the results highlight the superiority of the Med-Med filter over mean and median filters concerning restoration accuracy and preservation of image quality. Its ability to reduce NSE, increase PSNR, and keep high SSIM positions it as a highly effective choice for image restoration applications where visual quality is paramount. This study underscores the significance of interdisciplinary research in visual technology, emphasizing the necessity to address theoretical challenges and practical considerations for advancing more efficient, exact, and adaptable imaging systems with broad applications across various industries.

Author Contributions

Conceptualization, H.B.-H. and A.D.-S.; formal analysis, A.D.-S.; investigation, J.M.R.-P.; methodology, V.H.C.G.; writing—original draft preparation, A.B.-C. and V.H.C.G.; writing—review and editing, A.D.-S. and J.M.R.-P.; supervision, A.D.-S. All authors have read and agreed to the published version of the manuscript.

Funding

I would like to extend my deepest gratitude to the National Council of Humanities, Sciences, and Technologies (CONAHCYT) for their generous financial support under grant A1-S-43214. This funding has been instrumental in advancing my research, allowing me to achieve significant milestones and contribute valuable findings to the scientific community. Your support is profoundly appreciated and has made a considerable impact on the success of my project.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kim, H.J.; Shrestha, A.; Sapkota, E.; Pokharel, A.; Pandey, S.; Kim, C.S.; Shrestha, R. A study on the effectiveness of spatial filters on thermal image pre-processing and correlation technique for quantifying defect size. Sensors 2022, 22, 8965. [Google Scholar] [CrossRef] [PubMed]
  2. Mendez, T.; Nayak, S.G.; Kumar, P.V.; Kedlaya, K.V. Performance metric evaluation of error-tolerant adders for 2D image blending. Electronics 2022, 11, 2461. [Google Scholar] [CrossRef]
  3. Bigas, M.; Cabruja, E.; Forest, J.; Salvi, J. Review of CMOS image sensors. Microelectron. J. 2006, 37, 433–451. [Google Scholar] [CrossRef]
  4. Ohta, J. Smart CMOS Image Sensors and Applications; CRC Press: Boca Raton, FL, USA, 2020. [Google Scholar]
  5. Li, N.; Liu, T.; Li, H. An improved adaptive median filtering algorithm for radar image co-channel interference suppression. Sensors 2022, 22, 7573. [Google Scholar] [CrossRef] [PubMed]
  6. Cherniak, G.; Nemirovsky, J.; Nemirovsky, A.; Nemirovsky, Y. Modeling Signal-to-Noise Ratio of CMOS Image Sensors with a Stochastic Approach under Non-Stationary Conditions. Sensors 2023, 23, 7344. [Google Scholar] [CrossRef] [PubMed]
  7. El-Desouki, M.; Deen, M.J.; Fang, Q.; Liu, L.; Tse, F.; Armstrong, D. CMOS image sensors for high speed applications. Sensors 2009, 9, 430–444. [Google Scholar] [CrossRef] [PubMed]
  8. Yonemoto, K.; Sumi, H. A numerical analysis of a CMOS image sensor with a simple fixed-pattern-noise-reduction technology. IEEE Trans. Electron Devices 2002, 49, 746–753. [Google Scholar] [CrossRef]
  9. Klosowski, M.; Sun, Y. Fixed Pattern Noise Reduction and Linearity Improvement in Time-Mode CMOS Image Sensors. Sensors 2020, 20, 5921. [Google Scholar] [CrossRef] [PubMed]
  10. Ma, Z.; Yu, H.; Cui, K.; Yu, Y.; Tao, C. Design and Study of a Two-Dimensional (2D) All-Optical Spatial Mapping Module. Sensors 2024, 24, 2219. [Google Scholar] [CrossRef]
  11. Theuwissen, A.J. CMOS image sensors: State-of-the-art. Solid-State Electron. 2008, 52, 1401–1406. [Google Scholar] [CrossRef]
  12. Xiang, W.; Shen, J.; Zhang, L.; Zhang, Y. Infrared and Visual Image Fusion Based on a Local-Extrema-Driven Image Filter. Sensors 2024, 24, 2271. [Google Scholar] [CrossRef]
  13. Ungureanu, V.I.; Negirla, P.; Korodi, A. Image-Compression Techniques: Classical and “Region-of-Interest-Based” Approaches Presented in Recent Papers. Sensors 2024, 24, 791. [Google Scholar] [CrossRef]
  14. Bandala-Hernandez, H.; Rocha-Pérez, J.M.; Díaz-Sánchez, A.; Lemus-López, J.; Vázquez-Leal, H.; Díaz-Armendariz, A.; Ramírez-Angulo, J. Weighted median filters: An analog implementation. Integration 2016, 55, 227–231. [Google Scholar] [CrossRef]
  15. Díaz-Sánchez, A.; Ramírez-Angulo, J.; Lopez-Martin, A.; Sánchez-Sinencio, E. A fully parallel CMOS analog median filter. IEEE Trans. Circuits Syst. II Express Briefs 2004, 51, 116–123. [Google Scholar] [CrossRef]
  16. Hung, Y.C.; Shieh, S.H.; Tung, C.K. A real-time current-mode CMOS analog median filtering cell for system-on-chip applications. In Proceedings of the 2007 IEEE Conference on Electron Devices and Solid-State Circuits, Tainan, Taiwan, 20–22 December 2007; pp. 361–364. [Google Scholar]
  17. Yildirim, M. Analog circuit implementation based on median filter for salt and pepper noise reduction in image. Analog Integr. Circuits Signal Process. 2021, 107, 195–202. [Google Scholar] [CrossRef]
  18. Muñiz-Montero, C.; Ramírez-Salinas, M.A.; Villa-Vargas, L.A.; Molina-Lozano, H.; Ponce-Ponce, V.H.; Sánchez-Gaspariano, L.A.; Arellano-Gutiérrez, D. A compact CMOS Class-AB analog median filter. In Proceedings of the 2012 IEEE 3rd Latin American Symposium on Circuits and Systems (LASCAS), Playa del Carmen, Mexico, 29 February–2 March 2012; pp. 1–4. [Google Scholar]
  19. Jendernalik, W.; Blakiewicz, G.; Jakusz, J.; Szczepański, S. A nine-input 1.25 mW, 34 ns CMOS analog median filter for image processing in real time. Analog Integr. Circuits Signal Process. 2013, 76, 233–243. [Google Scholar] [CrossRef]
  20. Juneja, M.; Mohana, R. An improved adaptive median filtering method for impulse noise detection. Int. J. Recent Trends Eng. 2009, 1, 274. [Google Scholar]
  21. Li, H.; Mu, X.; Yang, Y.; Mason, A.J. Low power multimode electrochemical gas sensor array system for wearable health and safety monitoring. IEEE Sens. J. 2014, 14, 3391–3399. [Google Scholar] [CrossRef]
  22. Cevik, I.; Ay, S.U. An ultra-low power energy harvesting and imaging (EHI) type CMOS APS imager with self-power capability. IEEE Trans. Circuits Syst. I Regul. Pap. 2015, 62, 2177–2186. [Google Scholar] [CrossRef]
  23. Chen, Y.; Jaeger, R.C.; Suhling, J.C. CMOS sensor arrays for high resolution die stress mapping in packaged integrated circuits. IEEE Sens. J. 2013, 13, 2066–2076. [Google Scholar] [CrossRef]
  24. Wang, Z.; Bovik, A.C.; Sheikh, H.R.; Simoncelli, E.P. Image quality assessment: From error visibility to structural similarity. IEEE Trans. Image Process. 2004, 13, 600–612. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, Z.; Bovik, A.C. A universal image quality index. IEEE Signal Process. Lett. 2002, 9, 81–84. [Google Scholar] [CrossRef]
  26. Martin, A.J.; Alonso, A.H.; Ruiz, D.; Gude, I.; De Marziani, C.; Pérez, M.C.; Alvarez, F.J.; Gutierrez, C.; Urena, J. EMFi-based ultrasonic sensory array for 3D localization of reflectors using positioning algorithms. IEEE Sens. J. 2014, 15, 2951–2962. [Google Scholar] [CrossRef]
  27. Amaral, P.; Goes, J.; Paulino, N.F.; Steiger-Garção, A. An improved low-voltage low-power CMOS comparator to be used in high-speed pipeline ADCs. In Proceedings of the ISCAS (5), Phoenix-Scottsdale, AZ, USA, 26–29 May 2002; pp. 141–144. [Google Scholar]
  28. Razavi, B.; Wooley, B.A. Design techniques for high-speed, high-resolution comparators. IEEE J. Solid-State Circuits 1992, 27, 1916–1926. [Google Scholar] [CrossRef]
  29. Jendernalik, W.; Jakusz, J.; Blakiewicz, G.; Szczepański, S. CMOS implementation of an analogue median filter for image processing in real time. Bull. Pol. Acad. Sci.-Tech. Sci. 2013, 61, 725–730. [Google Scholar] [CrossRef]
  30. Jendernalik, W. On analog comparators for CMOS digital pixel applications. A comparative study. Bull. Pol. Acad. Sci.-Tech. Sci. 2016, 64, 271–278. [Google Scholar] [CrossRef]
Figure 1. Performance of different 2-D filters on a grayscale image.
Figure 1. Performance of different 2-D filters on a grayscale image.
Sensors 24 04213 g001
Figure 2. (a) Original image with focal plane, (b) Image with 20% impulse noise, (c) Filtered image using mean filter, (d) Filtered image using median filter, (e) Filtered image using median of medians filter, (f) Filtered image using mean of medians filter.
Figure 2. (a) Original image with focal plane, (b) Image with 20% impulse noise, (c) Filtered image using mean filter, (d) Filtered image using median filter, (e) Filtered image using median of medians filter, (f) Filtered image using mean of medians filter.
Sensors 24 04213 g002
Figure 3. MSE comparing the filtered image to the noise.
Figure 3. MSE comparing the filtered image to the noise.
Sensors 24 04213 g003
Figure 4. Power consumption of an IC implementation using Mean-Med and Med-Med algorithms: (a) logarithmic X-axis, (b) both logarithmic X- and Y-axes for improved visualization.
Figure 4. Power consumption of an IC implementation using Mean-Med and Med-Med algorithms: (a) logarithmic X-axis, (b) both logarithmic X- and Y-axes for improved visualization.
Sensors 24 04213 g004
Figure 5. The 5-input median extractor circuit.
Figure 5. The 5-input median extractor circuit.
Sensors 24 04213 g005
Figure 6. Transient response of the 5-input median filter applying the mask W = { 1 , 1 , 1 , 1 , 1 } : (a) simulation results, (b) experimental results, 400 mVpp signals at 10 KHz.
Figure 6. Transient response of the 5-input median filter applying the mask W = { 1 , 1 , 1 , 1 , 1 } : (a) simulation results, (b) experimental results, 400 mVpp signals at 10 KHz.
Sensors 24 04213 g006
Figure 7. WMF response applying mask W a = { 3 , 1 , 1 , 1 , 1 } : (a) MATLAB simulation and (b) experimental results, 400 mVpp signals at 10 KHz.
Figure 7. WMF response applying mask W a = { 3 , 1 , 1 , 1 , 1 } : (a) MATLAB simulation and (b) experimental results, 400 mVpp signals at 10 KHz.
Sensors 24 04213 g007
Figure 8. Scheme for the multilevel filter for median of medians extractor.
Figure 8. Scheme for the multilevel filter for median of medians extractor.
Sensors 24 04213 g008
Figure 9. Schematic circuit of the analog comparator in the Mean-Med implementation for image sensors.
Figure 9. Schematic circuit of the analog comparator in the Mean-Med implementation for image sensors.
Sensors 24 04213 g009
Figure 10. The 3 × 3 sensor array for the Med-Med filter accomplishment. CMOS Opto Technology.
Figure 10. The 3 × 3 sensor array for the Med-Med filter accomplishment. CMOS Opto Technology.
Sensors 24 04213 g010
Figure 11. Photography of 3 × 3 sensor array for the Med-Med filter accomplished using AMI 0.35 μ m CMOS Opto Technology.
Figure 11. Photography of 3 × 3 sensor array for the Med-Med filter accomplished using AMI 0.35 μ m CMOS Opto Technology.
Sensors 24 04213 g011
Figure 12. PCB for experimental measurements of Med-Med implementation.
Figure 12. PCB for experimental measurements of Med-Med implementation.
Sensors 24 04213 g012
Figure 13. Alignment of the optical fiber to Med-Med implementation.
Figure 13. Alignment of the optical fiber to Med-Med implementation.
Sensors 24 04213 g013
Figure 14. Mask 2 projected onto the 9-pixel array for Med-Med filter.
Figure 14. Mask 2 projected onto the 9-pixel array for Med-Med filter.
Sensors 24 04213 g014
Figure 15. Digital response for the Med-Med image sensor prototype for the M s k 2 exposure corresponding to M E D F 1 .
Figure 15. Digital response for the Med-Med image sensor prototype for the M s k 2 exposure corresponding to M E D F 1 .
Sensors 24 04213 g015
Figure 16. Digital response for the Med-Med image sensor prototype for the exposure corresponding to M s k 2 : M E D F 2 and M E D F 3 .
Figure 16. Digital response for the Med-Med image sensor prototype for the exposure corresponding to M s k 2 : M E D F 2 and M E D F 3 .
Sensors 24 04213 g016
Figure 17. Digital response for the Med-Med image sensor prototype for the M s k 2 exposure corresponding to M E D F 4 .
Figure 17. Digital response for the Med-Med image sensor prototype for the M s k 2 exposure corresponding to M E D F 4 .
Sensors 24 04213 g017
Figure 18. Digital response for the Med-Med image sensor prototype for the M s k 2 exposure corresponding to M E D F m e d .
Figure 18. Digital response for the Med-Med image sensor prototype for the M s k 2 exposure corresponding to M E D F m e d .
Sensors 24 04213 g018
Table 1. Mean squared error, peak signal-to-noise ratio, and structural similarity index are evaluated for a clean image compared with an image corrupted with 20 % impulse noise.
Table 1. Mean squared error, peak signal-to-noise ratio, and structural similarity index are evaluated for a clean image compared with an image corrupted with 20 % impulse noise.
Measure\FilterMeanMedianMed-MedMean-Med
MSE8965.3120.9115.271.2
PSNR8.627.327.529.6
SSIM0.0070.9770.9880.971
Table 2. Dimensions of 5-input median extractor circuit. 
Table 2. Dimensions of 5-input median extractor circuit. 
( W / L ) M 1 M 10 ( W / L ) M 11 M 12 I x V DD
40 μ m/1.2 μ m60 μ m/1.2 μ m10 nA1.8 V
Table 3. Input signals of the median extractor circuit.
Table 3. Input signals of the median extractor circuit.
WFM InputSignal (400 mVpp)Frequency
v i n 1 Sinusoidal10 kHz
v i n 2 Sinusoidal (90°)10 kHz
v i n 3 Triangular (45°)10 kHz
v i n 4 Triangular (45°)10 kHz
v i n 5 Square10 kHz
Table 4. Median filters comparison.
Table 4. Median filters comparison.
[18][16][19]This Work
Latency4.4 ns34 ns-5 μ s
Silicon Area0.03 mm   2 0.014 mm   2 0.004 mm   2 0.013 mm   2
Power13.5 μ W-1.25 mW90 nW
V-Supply3 V3.3 V1.8 V
F r e q @ C L 10 kHz @ 10 pF10 MHz @ ?-500 kHz @ 10 pF
Technology0.5 μ m0.35 μ m0.5 μ m
# of Inputs3935
Comparator Gain-3800 V/V-100 V/V
Accuracy-10 mV0.4 μ A discriminability40 mV
TypeAnalog VLSI
ICMR−0.9 V to 1.5 V0 V to 1.8 V-0.1 V to 1.7 V
Table 5. Comparator characteristics. 
Table 5. Comparator characteristics. 
( W / L ) 1 , 2 , 7 , 8 ( W / L ) 3 , 4 ( W / L ) 5 , 6 P o w e r d i s s S e t t t i m e V s u p p l y
9 μ /0.9 μ 3 μ /0.9 μ 16 μ /0.9 μ 385.63 μ W9 ns±1.65 V
Table 6. Array comparison.
Table 6. Array comparison.
This Work’s Median Filter[29] (Pixel Array, 1 Median Filter)[30] (Comparator)[17] (Median Filter)
TechnologyAMS 0.35 μ m0.35 μ mAMS 0.35 μ mTSMC 0.18 μ m
Comparator gain100 V/V3800 V/V>1000 V/VNot reported
Vsupply1.8 V3.3 V3.3 V2 V (vdd = vss = 1V)
Power450 nW @Med-Med operation
(5 filters active),
90 nW (single median filter)
@500 KHz, 1.8 Vpp
1.25 mW7–22 μ WNot reported
Dynamic range0.1 V–1.7 V0–1.8 VNot reported0–255 μ A
Type of filterWeighted analog median—SubthresholdAnalog median filter *N/AAnalog median filter *
number of inputsMed-Med mode = 21, Standalone = 59N/A9
Lab measurementsYesYesSimulationsSimulations
Data sorting requiredNoYes **N/AYes
Accuracy40 mV10 mVN/ANot reported
Silicon area0.013 mm   2 0.014 mm   2 Not reportedNot reported
Additional comments9-pixel array, including photosensors,
readout, comparator, median filters,
and digital logic for reading
32 × 32 pixel array with
external FPGA logic and
single median filter
for the whole array
Explores the use
of analog comparators
in CMOS imagers
A median circuit implementation
* MAX-MIN circuit required. ** Median operation is performed after readout via one median filter for whole array.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bandala-Hernandez, H.; Bautista-Castillo, A.; Rocha-Pérez, J.M.; Carbajal Gómez, V.H.; Díaz-Sánchez, A. Exploring the Efficacy of Nonlinear Filters in CMOS for 2-D Signal Processing for Image Quality Enhancement. Sensors 2024, 24, 4213. https://doi.org/10.3390/s24134213

AMA Style

Bandala-Hernandez H, Bautista-Castillo A, Rocha-Pérez JM, Carbajal Gómez VH, Díaz-Sánchez A. Exploring the Efficacy of Nonlinear Filters in CMOS for 2-D Signal Processing for Image Quality Enhancement. Sensors. 2024; 24(13):4213. https://doi.org/10.3390/s24134213

Chicago/Turabian Style

Bandala-Hernandez, Hector, Alejandro Bautista-Castillo, José Miguel Rocha-Pérez, Victor Hugo Carbajal Gómez, and Alejandro Díaz-Sánchez. 2024. "Exploring the Efficacy of Nonlinear Filters in CMOS for 2-D Signal Processing for Image Quality Enhancement" Sensors 24, no. 13: 4213. https://doi.org/10.3390/s24134213

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop