DOES MEDIAN FILTERING TRULY PRESERVE EDGES BETTER THAN LINEAR FILTERING?

被引:166
作者
Arias-Castro, Ery [1 ]
Donoho, David L. [2 ]
机构
[1] Univ Calif San Diego, Dept Math, La Jolla, CA 92093 USA
[2] Stanford Univ, Dept Stat, Stanford, CA 94305 USA
关键词
Linear filtering; kernel smoothing; median filtering; running median; image denoising; minimax estimation; nonparametric regression; MINIMAX ESTIMATION; WAVELET SHRINKAGE;
D O I
10.1214/08-AOS604
中图分类号
O21 [概率论与数理统计]; C8 [统计学];
学科分类号
020208 ; 070103 ; 0714 ;
摘要
Image processing researchers commonly assert that "median filtering is better than linear filtering for removing noise in the presence of edges." Using a straightforward large-n decision-theory framework, this folk-theorem is seen to be false in general. We show that median filtering and linear filtering have similar asymptotic worst-case mean-squared error (MSE) when the signal-to-noise ratio (SNR) is of order 1, which corresponds to the case of constant per-pixel noise level in a digital signal. To see dramatic benefits of median smoothing in an asymptotic setting, the per-pixel noise level should tend to zero (i.e., SNR should grow very large). We show that a two-stage median filtering using two very different window widths can dramatically outperform traditional linear and median filtering in settings where the underlying object has edges. In this two-stage procedure, the first pass, at a fine scale, aims at increasing the SNR. The second pass, at a coarser scale, correctly exploits the nonlinearity of the median. Image processing methods based on nonlinear partial differential equations (PDEs) are often said to improve on linear filtering in the presence of edges. Such methods seem difficult to analyze rigorously in a decision-theoretic framework. A popular example is mean curvature motion (MCM), which is formally a kind of iterated median filtering. Our results on iterated median filtering suggest that some PDE-based methods are candidates to rigorously outperform linear filtering in an asymptotic framework.
引用
收藏
页码:1172 / 1206
页数:35
相关论文
共 36 条
[1]  
[Anonymous], 1982, IMAGE ANAL MATH MORP
[2]  
BARNER K, 2003, NONLINEAR SIGNAL LIN
[3]   Removing noise and preserving details with relaxed median filters [J].
Ben Hamza, A ;
Luque-Escamilla, PL ;
Martínez-Aroza, J ;
Román-Roldán, R .
JOURNAL OF MATHEMATICAL IMAGING AND VISION, 1999, 11 (02) :161-177
[4]   DETERMINISTIC PROPERTIES OF ANALOG MEDIAN FILTERS [J].
BOTTEMA, MJ .
IEEE TRANSACTIONS ON INFORMATION THEORY, 1991, 37 (06) :1629-1640
[5]  
Brandt J, 1998, UTILITAS MATHEMATICA, V54, P111
[6]  
Candès EJ, 2002, ANN STAT, V30, P784
[7]   Partial differential equations and mathematical morphology [J].
Cao, F .
JOURNAL DE MATHEMATIQUES PURES ET APPLIQUEES, 1998, 77 (09) :909-941
[8]   Vector median filters, inf-sup operations, and coupled PDE's: Theoretical connections [J].
Caselles, V ;
Sapiro, G ;
Chung, DK .
JOURNAL OF MATHEMATICAL IMAGING AND VISION, 2000, 12 (02) :109-119
[9]   Wedgelets: Nearly minimax estimation of edges [J].
Donoho, DL .
ANNALS OF STATISTICS, 1999, 27 (03) :859-897
[10]   Nonlinear pyramid transforms based on median-interpolation [J].
Donoho, DL ;
Yu, TPY .
SIAM JOURNAL ON MATHEMATICAL ANALYSIS, 2000, 31 (05) :1030-1061