数字图像处理冈萨雷斯英文版Chapter05图像复原.pptVIP

Algorithm: Level A: A1= zmedian – zmin A2= zmedian – zmax If A1 0 and A2 0, goto level B Else increase window size If window size = Smax repeat level A Else return zxy Level B: B1= zxy – zmin B2= zxy – zmax If B1 0 and B2 0, return zxy Else return zmedian Adaptive Median Filter zmin = minimum gray level value in Sxy zmax = maximum gray level value in Sxy zmedian = median of gray levels in Sxy zxy = gray level value at pixel (x,y) Smax = maximum allowed size of Sxy where Purpose: want to remove impulse noise while preserving edges Level A: A1= zmedian – zmin A2= zmedian – zmax Else ? Window is not big enough increase window size If window size = Smax repeat level A Else return zxy ? zmedian is not an impulse B1= zxy – zmin B2= zxy – zmax If B1 0 and B2 0, ? zxy is not an impulse return zxy ? to preserve original details Else return zmedian ? to remove impulse Adaptive Median Filter: How it works If A1 0 and A2 0, goto level B Level B: Determine whether zmedian is an impulse or not Determine whether zxy is an impulse or not (Images from Rafael C. Gonzalez and Richard E. Wood, Digital Image Processing, 2nd Edition. Adaptive Median Filter: Example Image corrupted by salt-and-pepper noise with pa=pb= 0.25 Image obtained using a 7x7 median filter Image obtained using an adaptive median filter with Smax = 7 More small details are preserved (Images from Rafael C. Gonzalez and Richard E. Wood, Digital Image Processing, 2nd Edition. Estimation of Degradation Model Degradation model: Purpose: to estimate h(x,y) or H(u,v) Methods: 1. Estimation by Image Observation 2. Estimation by Experiment 3. Estimation by Modeling or Why? If we know exactly h(x,y), regardless of noise, we can do deconvolution to get f(x,y) back from g(x,y). Estimation by Image Observation f(x,y) f(x,y)*h(x,y) g(x,y) Subimage Reconstructed Subimage DFT DFT Restoration process by estimation Original