Image Sequence Restoration by Median Filtering

Total Page:16

File Type:pdf, Size:1020Kb

Image Sequence Restoration by Median Filtering Rochester Institute of Technology RIT Scholar Works Theses 2004 Image sequence restoration by median filtering Shawn R. Jackson Follow this and additional works at: https://scholarworks.rit.edu/theses Recommended Citation Jackson, Shawn R., "Image sequence restoration by median filtering" (2004). Thesis. Rochester Institute of Technology. Accessed from This Thesis is brought to you for free and open access by RIT Scholar Works. It has been accepted for inclusion in Theses by an authorized administrator of RIT Scholar Works. For more information, please contact [email protected]. Image Sequence Restoration by Median Filtering by Shawn R. Jackson A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of Master of Science in Computer Engineering Supervised by Dr. Andreas Savakis Department of Computer Engineering Kate Gleason College of Engineering Rochester Institute of Technology Rochester, NY November of 2004 Approved By: Andreas Savakis Dr. Andreas Savakis, Professor and Department Head Primary Advisor - R.l T. Dept. of Computer Engineering Kenneth Hsu Dr. Kenneth Hsu, Professor Secondary Advisor - R.l T. Dept. of Computer Engineering Greg P. Semeraro Dr. Greg Semeraro, Assistant Professor Secondary Advisor - R.l T. Dept. of Computer Engineering Thesis Release Permission Form Rochester Institute of Technology Kate Gleason College of Engineering Title: Image Sequence Restoration by Median Filtering 1, Shawn R. Jackson, hereby grant permission to the WaUace Memorial Library to reproduce my thesis in whole or part. Shawn R. Jackson Shawn R. Jackson 11 Dedication This work is dedicated to Mom, Dad, Nicole, and Kristina. Thanks for pushing me the whole way. 111 Acknowledgements I would like to show my appreciation to Dr. Savakis, Dr. Hsu, and Dr. Semeraro, Thank you for your guidance and insight. I would also like to thank Eastman Kodak Company for providing the necessary tools and hardware. IV Abstract Median filters are non-linear filters that fit in the generic category of order-statistic filters. Median filters are widely used for reducing random defects, commonly characterized by impulse pepper" or "salt and noise in a single image. Motion estimation is the process of estimating the displacement vector between like pixels in the current frame and the reference frame. When dealing with a motion sequence, the motion vectors are the key for operating on corresponding pixels in several frames. This work explores the use of various motion estimation algorithms in combination with various median filter algorithms to provide noise suppression. The results are compared using two sets of metrics: performance-based and objective image quality-based. These results are used to determine the best motion estimation / median filter combination for image sequence restoration. The primary goals of this work are to implement a motion estimation and median filter algorithm in hardware and develop and benchmark a flexible software alternative restoration process. There are two unique median filter algorithms to this work. The first filter is a modification to a single frame adaptive median filter. The modification applied motion compensation and temporal concepts. The other is an adaptive extension to the multi-level (ML3D) filter, called adaptive multi-level (AML3D) filter. The extension provides adaptable filter window sizes to the multiple filter sets that comprise the ML3D filter. The adaptive median filter is capable of filtering an image in 26.88 seconds per frame and results in a PSNR improvement of 5.452dB. The AML3D is capable of filtering an image in 14.73 seconds per frame and results in a PSNR improvement of 6.273dB. The AML3D is a suitable alternative to the other median filters. Table of Contents Chapter 1: Introduction 1 Chapter 2: Techniquesfor Motion Estimation ofImage Sequences 6 2.1 Motion Estimation 9 2.2 Block Matching Methods 9 2.3 Gradient Based Methods 14 2.4 Pel-Recursive Based Methods 20 2.5 Selected Motion Estimation Algorithms 22 2.5.1 Simplified Black's Motion Estimation Algorithm 23 2.5.2 Horn-Schunck Motion Estimation Algorithm 24 Chapter 3: Techniquesfor Spado-Temporal Filtering ofImage Sequences. 30 3.1 Temporal Filtering 31 3.2 Spatio-Temporal Filtering 33 3.3 Median Filtering 34 3.4 Selected Median Filtering Algorithms 37 3.4.1 Motion Compensated Median Filter Algorithm 40 3.4.2 ML3D Filter Algorithm 41 3.4.3 ML3Dex Filter Algorithm 43 3.4.4 Arce Bi-directional Multi-level Median Filter Algorithm 46 3.4.5 Adaptive Median Filter Algorithm 47 3.4.6 Adaptive ML3D Filter Algorithm 50 3.4.7 Summary of Selected Median Filter Algorithms 53 Chapter 4: Software Implementation ofRestoration Process 55 4.1 Software Interaction Diagram 59 Chapter 5: Hardware Implementation ofRestoration Process 65 5.1 Description of Hardware Platform 65 vi 5.1.1 Overview ofWildStar PCI Hardware 65 5.1.2 Diagram ofWildStar PCI Hardware 68 5.2 Detail Investigation of Hardware Algorithms 70 5.2.1 Data Flow Analysis ofHardware Design 70 5.2.2 Detailed Investigation ofPE1 Design 73 5.2.3 Detailed Investigation of PEO Design 76 5.3 Host Software Roles and Responsibilities 81 6.1 Image Sequence Selection and Preparation 84 6.2 Performance and Timing Based Results 88 6.3 Image Quality Based Results 94 6.4 General Conclusions 98 6.5 Path Forward 103 6.6 Summary 106 Bibliography 6 Appendix: Example Restoration Results 112 vn List of Figures Figure 2.1: Sample Motion Between Frames 8 Figure 2.2: Full Search Block Matching Window Size 13 Figure 2.3: OSA with Displacement of 7 and Minimum Distortion Function at (-1, -1) 15 Figure 2.4: Hierarchical Framework for Efficient Motion Estimation 18 Figure 2.5: Aperture Problem 19 Figure 2.6: Windows Area Selection in Frames 25 Figure 3.1: Example 2-Stage MMF Topology with 3 Sub-Filter Windows 37 Figure 3.2: Original Image ofPeppers 38 Figure 3.3: Corrupted Image with 10% Salt and Pepper Noise 38 Figure 3.4: Corrupted Image Filtered with 3x3 Median Filter 39 Figure 3.5: Corrupted Image Filtered with MMF 39 Figure 3.6: Diagram Depicting Motion Compensated Median Filtering 41 Figure 3.7: Diagram Depicting ML3D Filtering 43 Figure 3.8: Diagram Depicting ML3Dex Filtering 45 Figure 3.9: Diagram Depicting the First Stage of Arce Bi-directional MMF 48 Figure 3.10: Diagram Depicting the AML3D Filter 53 Figure 4.1: Software Interaction Diagram 61 Figure 4.2: Software Interaction Diagram Focused on Image Reading and Writing 62 Figure 4.3: Software Interaction Diagram Featuring Motion Estimation Components 63 Figure 4.4: Software Interaction Diagram Featuring Filtering Components 64 Figure 5.1: WildStar PCI Test Board 68 Figure 5.2: WildStar PCI Mezzanine Memory Module 69 Figure 5.3: Block Diagram ofWildStar PCI Connections 69 Figure 5.4: Overall Image Flow and Design Partitioning 72 Figure 5.5: Memory Flow and Processing Element Responsibilities 72 Figure 5.6: Block Diagram ofWildStar PCI PE1 75 Figure 5.7: Register Map for WildStar PCI PE1 75 Figure 5.8: Block Diagram of WildStar PCI PEO 77 Figure 5.9: Register Map for WildStar PCI PEO 78 Figure 6.1: Original Image from 13 Days Sequence 85 Figure 6.2: Original Image from BW Sequence 86 Figure 6.3: Original Image from Martial Sequence 86 Figure 6.4: Original Image from NASA Sequence 87 Figure 6.5: Original Image from Walk Sequence 87 Figure 6.6: More Sample Resultant Images from BW Sequence 99 Vlll Figure 6.7: Time / Quality Plot for Simplified Black's Motion Estimator 102 Figure 6.8: Time /Quality Plot for Horn-Schunck Motion Estimator 102 IX List of Tables Table 6.1: Table ofProperties ofthe 13 Days Image Sequences 85 Table 6.2: Table ofTotal Times for BW Sequence on Hardware Implementation 90 Table 6.3: Table ofTotal Times for BW Sequence on Software Implementation Using Black's Method for Motion Estimation and Adaptive Median Filter 90 Table 6.4: Table ofTotal Times for Averages ofEach Sequence 91 Table 6.5: Table ofAverage Times for Averages of Each Sequence 91 Table 6.6: Table ofRanked Overall Performance for Filter / Motion Estimation Combination 93 Table 6.7: Table ofImage Quality Statistics for BW Image Sequence using Adaptive and AML3D Median Filters using Black's Method for Motion Estimator 95 Table 6.8: Table of Image Quality Improvements for BW Image Sequence using Adaptive and AML3D Median Filters using Black's Method for Motion Estimator 95 Table 6.9: Table of Image Quality Statistics for Grayscale Sequences Averaged 96 Table 6.10: Table of Image Quality Statistics for Color Sequences Averaged 96 Table 6.11: Table of Image Quality Statistics for All Sequences Averaged 97 Table 6.12: Table ofRanked Overall Improvement for Filter / Motion Estimation Combination 98 Glossary AML3D Adaptive ML3D Filter - an adaptive MMF proposed in this work BM Block Matching motion estimation algorithm Code Value Intensity level of a pixel DFD Displaced Frame Difference DMA Direct Memory Access FPGA Field Programmable Gate Array LUT Look-Up Table MAE Mean Absolute Error ML3DEx Extended ML3D MMF MMF Multi-level Median Filter MSE Mean Squared Error Motion Field Set ofmotion vectors for an image Motion Vector Displacement vector from a pixel in one frame to where it is in the next frame PDC Pixel Difference Classification Pixel Picture element RMSE Root Mean Squared Error VHDL VHSIC Hardware Description Language - High level VLSI design language VHSIC Very High Speed Integrated Circuit VLSI Very Large Scale Integration XI Chapter 1 : Introduction As the demand for DVD and digital content increases, the need for digital film restoration also increases. The desire to preserve old films from degradation is one example of the need for efficient digital film restoration. The cellulose nitrate based film used during the silent picture era is highly volatile and susceptible to age degradation.
Recommended publications
  • A Convolutional Neural Networks Denoising Approach for Salt and Pepper Noise
    A Convolutional Neural Networks Denoising Approach for Salt and Pepper Noise Bo Fu, Xiao-Yang Zhao, Yi Li, Xiang-Hai Wang and Yong-Gong Ren School of Computer and Information Technology, Liaoning Normal University, China [email protected] Abstract. The salt and pepper noise, especially the one with extremely high percentage of impulses, brings a significant challenge to image denoising. In this paper, we propose a non-local switching filter convolutional neural network denoising algorithm, named NLSF-CNN, for salt and pepper noise. As its name suggested, our NLSF-CNN consists of two steps, i.e., a NLSF processing step and a CNN training step. First, we develop a NLSF pre-processing step for noisy images using non-local information. Then, the pre-processed images are divided into patches and used for CNN training, leading to a CNN denoising model for future noisy images. We conduct a number of experiments to evaluate the effectiveness of NLSF-CNN. Experimental results show that NLSF-CNN outperforms the state-of-the-art denoising algorithms with a few training images. Keywords: Image denoising, Salt and pepper noise, Convolutional neural networks, Non-local switching filter. 1 Introduction Images are always polluted by noise during image acquisition and transmission, resulting in low-quality images. The removal of such noise is a must before subsequent image analysis tasks. To achieve this, researchers have proposed many denoising algorithms [1-7], where the representative ones include non-local means (NLM) [1] and block-matching 3D filtering (BM3D) [2].In this paper, we focus on the salt and pepper noise, where pixels exhibit maximal or minimal values with some pre-defined probability.
    [Show full text]
  • 22Nd International Congress on Acoustics ICA 2016
    Page intentionaly left blank 22nd International Congress on Acoustics ICA 2016 PROCEEDINGS Editors: Federico Miyara Ernesto Accolti Vivian Pasch Nilda Vechiatti X Congreso Iberoamericano de Acústica XIV Congreso Argentino de Acústica XXVI Encontro da Sociedade Brasileira de Acústica 22nd International Congress on Acoustics ICA 2016 : Proceedings / Federico Miyara ... [et al.] ; compilado por Federico Miyara ; Ernesto Accolti. - 1a ed . - Gonnet : Asociación de Acústicos Argentinos, 2016. Libro digital, PDF Archivo Digital: descarga y online ISBN 978-987-24713-6-1 1. Acústica. 2. Acústica Arquitectónica. 3. Electroacústica. I. Miyara, Federico II. Miyara, Federico, comp. III. Accolti, Ernesto, comp. CDD 690.22 ISBN 978-987-24713-6-1 © Asociación de Acústicos Argentinos Hecho el depósito que marca la ley 11.723 Disclaimer: The material, information, results, opinions, and/or views in this publication, as well as the claim for authorship and originality, are the sole responsibility of the respective author(s) of each paper, not the International Commission for Acoustics, the Federación Iberoamaricana de Acústica, the Asociación de Acústicos Argentinos or any of their employees, members, authorities, or editors. Except for the cases in which it is expressly stated, the papers have not been subject to peer review. The editors have attempted to accomplish a uniform presentation for all papers and the authors have been given the opportunity to correct detected formatting non-compliances Hecho en Argentina Made in Argentina Asociación de Acústicos Argentinos, AdAA Camino Centenario y 5006, Gonnet, Buenos Aires, Argentina http://www.adaa.org.ar Proceedings of the 22th International Congress on Acoustics ICA 2016 5-9 September 2016 Catholic University of Argentina, Buenos Aires, Argentina ICA 2016 has been organised by the Ibero-american Federation of Acoustics (FIA) and the Argentinian Acousticians Association (AdAA) on behalf of the International Commission for Acoustics.
    [Show full text]
  • Analysis of the Entropy-Guided Switching Trimmed Mean Deviation-Based Anisotropic Diffusion Filter
    Analysis of the Entropy-guided Switching Trimmed Mean Deviation-based Anisotropic Diffusion filter U. A. Nnolim Department of Electronic Engineering, University of Nigeria Nsukka, Enugu, Nigeria Abstract This report describes the experimental analysis of a proposed switching filter-anisotropic diffusion hybrid for the filtering of the fixed value (salt and pepper) impulse noise (FVIN). The filter works well at both low and high noise densities though it was specifically designed for high noise density levels. The filter combines the switching mechanism of decision-based filters and the partial differential equation-based formulation to yield a powerful system capable of recovering the image signals at very high noise levels. Experimental results indicate that the filter surpasses other filters, especially at very high noise levels. Additionally, its adaptive nature ensures that the performance is guided by the metrics obtained from the noisy input image. The filter algorithm is of both global and local nature, where the former is chosen to reduce computation time and complexity, while the latter is used for best results. 1. Introduction Filtering is an imperative and unavoidable process, due to the need for sending or retrieving signals, which are usually corrupted by noise intrinsic to the transmission or reception device or medium. Filtering enables the separation of required and unwanted signals, (which are classified as noise). As a result, filtering is the mainstay of both analog and digital signal and image processing with convolution as the engine of linear filtering processes involving deterministic signals [1] [2]. For nonlinear processes involving non-deterministic signals, order statistics are employed to extract meaningful information [2] [3].
    [Show full text]
  • Total Variation Denoising for Optical Coherence Tomography
    Total Variation Denoising for Optical Coherence Tomography Michael Shamouilian, Ivan Selesnick Department of Electrical and Computer Engineering, Tandon School of Engineering, New York University, New York, New York, USA fmike.sha, [email protected] Abstract—This paper introduces a new method of combining total variation denoising (TVD) and median filtering to reduce noise in optical coherence tomography (OCT) image volumes. Both noise from image acquisition and digital processing severely degrade the quality of the OCT volumes. The OCT volume consists of the anatomical structures of interest and speckle noise. For denoising purposes we model speckle noise as a combination of additive white Gaussian noise (AWGN) and sparse salt and pepper noise. The proposed method recovers the anatomical structures of interest by using a Median filter to remove the sparse salt and pepper noise and by using TVD to remove the AWGN while preserving the edges in the image. The proposed method reduces noise without much loss in structural detail. When compared to other leading methods, our method produces similar results significantly faster. Index Terms—Total Variation Denoising, Median Filtering, Optical Coherence Tomography I. INTRODUCTION Fig. 1. 2D slices of a sample retinal OCT volume. Optical Coherence Tomography (OCT), a non-invasive in vivo imaging technique, provides images of internal tissue mi- crostructures of the eye [1]. OCT imaging works by directing light beams at a target tissue and then capturing and processing the backscattered light [2]. To create a 3D volume image of the eye, the OCT device scans the eye laterally [2]. However, current OCT imaging techniques inherently introduce speckle noise [3].
    [Show full text]
  • Transportation Noise and Annoyance Related to Road
    Méline et al. International Journal of Health Geographics 2013, 12:44 INTERNATIONAL JOURNAL http://www.ij-healthgeographics.com/content/12/1/44 OF HEALTH GEOGRAPHICS RESEARCH Open Access Transportation noise and annoyance related to road traffic in the French RECORD study Julie Méline1,2*, Andraea Van Hulst3,4, Frédérique Thomas5, Noëlla Karusisi1,2 and Basile Chaix1,2 Abstract Road traffic and related noise is a major source of annoyance and impairment to health in urban areas. Many areas exposed to road traffic noise are also exposed to rail and air traffic noise. The resulting annoyance may depend on individual/neighborhood socio-demographic factors. Nevertheless, few studies have taken into account the confounding or modifying factors in the relationship between transportation noise and annoyance due to road traffic. In this study, we address these issues by combining Geographic Information Systems and epidemiologic methods. Street network buffers with a radius of 500 m were defined around the place of residence of the 7290 participants of the RECORD Cohort in Ile-de-France. Estimated outdoor traffic noise levels (road, rail, and air separately) were assessed at each place of residence and in each of these buffers. Higher levels of exposure to noise were documented in low educated neighborhoods. Multilevel logistic regression models documented positive associations between road traffic noise and annoyance due to road traffic, after adjusting for individual/ neighborhood socioeconomic conditions. There was no evidence that the association was of different magnitude when noise was measured at the place of residence or in the residential neighborhood. However, the strength of the association between neighborhood noise exposure and annoyance increased when considering a higher percentile in the distribution of noise in each neighborhood.
    [Show full text]
  • A Spatial Median Filter for Noise Removal in Digital Images
    A Spatial Median Filter for Noise Removal in Digital Images James Church, Dr. Yixin Chen, and Dr. Stephen Rice Computer Science and Information System, University of Mississippi [email protected],{ychen,rice}@cs.olemiss.edu Abstract Each of the algorithms coveredcan be applied to one di- mensional as well as two dimensional signals. Figure 1 In this paper, six different image filtering algorithms are demonstrates five common filtering algorithms applied compared based on their ability to reconstruct noise- to an original image. affected images. The purpose of these algorithms is to remove noise from a signal that might occur through the transmission of an image. A new algorithm, the Spatial Median Filter, is introduced and compared with the cur- rent image smoothing techniques. Experimental results demonstrate that the proposed algorithm is comparable (a) (b) (c) to popular image smoothing algorithms. In addition, a modification to this algorithm is introduced to achieve more accurate reconstructions over other popular tech- niques. (d) (e) (f) Figure 1. Examples of common filtering ap- 1. Smoothing Algorithms proaches. (a) Original Image (b) Mean Filtering (c) Median Filtering (d) Root signal of Median The inexpensiveness and simplicity of point-and- Filtering (e) Component-wise Median Filtering shoot cameras, combined with the speed at which bud- (f) Vector Median Filtering ding photographerscan send their photos over the Inter- net to be viewed by the world, makes digital photogra- The simplest of smoothing algorithms is the Mean phy a popular hobby. With each snap of a digital pho- Filter as defined in (1). The Mean Filter is a linear filter tograph, a signal is transmitted from a photon sensor which uses a mask over each pixel in the signal.
    [Show full text]
  • 2007 Registration Document
    2007 REGISTRATION DOCUMENT (www.renault.com) REGISTRATION DOCUMENT REGISTRATION 2007 Photos cre dits: cover: Thomas Von Salomon - p. 3 : R. Kalvar - p. 4, 8, 22, 30 : BLM Studio, S. de Bourgies S. BLM Studio, 30 : 22, 8, 4, Kalvar - p. R. 3 : Salomon - p. Von Thomas cover: dits: Photos cre 2007 REGISTRATION DOCUMENT INCLUDING THE MANAGEMENT REPORT APPROVED BY THE BOARD OF DIRECTORS ON FEBRUARY 12, 2008 This Registration Document is on line on the website www .renault.com (French and English versions) and on the AMF website www .amf- france.org (French version only). TABLE OF CONTENTS 0 1 05 RENAULT AND THE GROUP 5 RENAULT AND ITS SHAREHOLDERS 157 1.1 Presentation of Renault and the Group 6 5.1 General information 158 1.2 Risk factors 24 5.2 General information about Renault’s share capital 160 1.3 The Renault-Nissan Alliance 25 5.3 Market for Renault shares 163 5.4 Investor relations policy 167 02 MANAGEMENT REPORT 43 06 2.1 Earnings report 44 MIXED GENERAL MEETING 2.2 Research and development 62 OF APRIL 29, 2008: PRESENTATION 2.3 Risk management 66 OF THE RESOLUTIONS 171 The Board first of all proposes the adoption of eleven resolutions by the Ordinary General Meeting 172 Next, six resolutions are within the powers of 03 the Extraordinary General Meeting 174 SUSTAINABLE DEVELOPMENT 79 Finally, the Board proposes the adoption of two resolutions by the Ordinary General Meeting 176 3.1 Employee-relations performance 80 3.2 Environmental performance 94 3.3 Social performance 109 3.4 Table of objectives (employee relations, environmental
    [Show full text]
  • Computational Photography Filtering: Smoothing Objectives Example
    Filtering: Smoothing • Point-process and Neighboring Pixels Computations on an Image for Image Smoothing Computational Photography CS 4475/6475 Maria Hybinette 1 Maria Hybinette 2 Maria Hybinette Objectives From Pixel/Point Operations to Groups of Pixels • Smooth an image while considering a neighborhood n of pixels – Average filtering (use the average of n) – Median filtering (use the median of n) • special non-linear filtering and smoothing • Previously we used pixel to pixel arithmetic (e.g., add(), sub())) approach • How to Smooth a Signal? • Approach: Use a group of pixels : Neighboring Values of a neighborhood. 1. Moving Average [1 1 1 1 1] X 1/5 2. Weighted Moving Average [1 4 6 4 1] X 1/16 3 Maria Hybinette 4 Maria Hybinette Example: 3x3 Kernel Example: 3x3 Kernel • Smoothing Process over an Image using Averages • 1*90/9=10, 2*90/9=20, …. 5 Maria Hybinette 6 Maria Hybinette Results in Shades of Gray • What about the edges? 7 Maria Hybinette 8 Maria Hybinette Edges Mathematical Notation, Observations, and Computing the Result Matrix • A Small image ‘h’ is “rubbed” over (larger) image ‘F’ that we are smoothing • Kernel: h – Example: 3x3 area around each original pixel used – Neighborhood size – k=1, window size = 3. – k >1, • Window size is 2k+1 – k=1, 2(1)+1 = 3 – Our kernel is 3x3 – Recall : A pixel is referenced by (i, j) • Expand the image we are smoothing: • Want Result Matrix ‘G’ – Wrap around pixels. (3x3 filter needs one more pixel – F à G using h. on each border, 5x5 needs 2) – G = h (operation) F – 9 Copy edges Maria Hybinette 10 Maria Hybinette Computing the Values ofA Result Mathematical Matrix RepresentationSmoothing when for a- iSmoothing are all 1s.
    [Show full text]
  • Noise Reduction by Vector Median Filtering
    GEOPHYSICS, VOL. 78, NO. 3 (MAY-JUNE 2013); P. V79–V86, 18 FIGS., 1 TABLE. 10.1190/GEO2012-0232.1 Noise reduction by vector median filtering Yike Liu1 1984; Marfurt, 2006), Radon transform (Sacchi and Porsani, 1999; ABSTRACT Sacchi et al., 2004), edge-preserving smoothing (Luo et al., 2002), and scalar median filter (SMF) (Mi and Margrave, 2000). Compared The scalar median filter (SMF) is often used to reduce to other methods, the SMF filtering technique often produces less noise in scalar geophysical data. We present an extension smearing among adjacent samples after noise attenuation. For ex- of the SMF to a vector median filter (VMF) for suppressing ample, the SMF can remove an abnormal impulse from seismic re- noise contained in geophysical data represented by multidi- cords without smearing the impulse into its nearby samples as the mensional, multicomponent vector fields. Although the mean or f-k filters do. The SMF can also be used to separate the up- SMF can be applied to each component of a vector field in- and downgoing wavefields in vertical seismic profile (VSP) data dividually, the VMF is applied to all components simulta- because it often reduces interference between these two wavefields. neously. Like the SMF, the VMF intends to suppress Despite the fact that numerous types of geophysical data are random noise while preserving discontinuities in the vector naturally represented by vector fields, the vector median filter fields. Preserving such discontinuities is essential for ex- (VMF) is seldom employed in data processing. On the other hand, ploration geophysics because discontinuities often manifest SMF is commonly used in exploration geophysics.
    [Show full text]
  • Denoising of X-Ray Images Using the Adaptive Algorithm Based on the LPA-RICI Algorithm
    Journal of Imaging Article Denoising of X-ray Images Using the Adaptive Algorithm Based on the LPA-RICI Algorithm Ivica Mandi´c,Hajdi Pei´c,Jonatan Lerga * ID and Ivan Štajduhar Faculty of Engineering, University of Rijeka, Vukovarska 58, 51000 Rijeka, Croatia; [email protected] (I.M.); [email protected] (H.P.); [email protected] (I.Š.) * Correspondence: [email protected]; Tel.: +385-51-651-583 Received: 5 December 2017; Accepted: 2 February 2018; Published: 5 February 2018 Abstract: Diagnostics and treatments of numerous diseases are highly dependent on the quality of captured medical images. However, noise (during both acquisition and transmission) is one of the main factors that reduce their quality. This paper proposes an adaptive image denoising algorithm applied to enhance X-ray images. The algorithm is based on the modification of the intersection of confidence intervals (ICI) rule, called relative intersection of confidence intervals (RICI) rule. For each image pixel apart, a 2D mask of adaptive size and shape is calculated and used in designing the 2D local polynomial approximation (LPA) filters for noise removal. One of the advantages of the proposed method is the fact that the estimation of the noise free pixel is performed independently for each image pixel and thus, the method is applicable for easy parallelization in order to improve its computational efficiency. The proposed method was compared to the Gaussian smoothing filters, total variation denoising and fixed size median filtering and was shown to outperform them both visually and in terms of the peak signal-to-noise ratio (PSNR) by up to 7.99 dB.
    [Show full text]
  • Transportation Noise and Annoyance Related to Road Traffic
    Transportation noise and annoyance related to road traffic in the French RECORD study. Julie Méline, Andraea van Hulst, Frédérique Thomas, Noëlla Karusisi, Basile Chaix To cite this version: Julie Méline, Andraea van Hulst, Frédérique Thomas, Noëlla Karusisi, Basile Chaix. Transportation noise and annoyance related to road traffic in the French RECORD study.. International Journal of Health Geographics, BioMed Central, 2013, 12 (1), pp.44. 10.1186/1476-072X-12-44. inserm- 00870153 HAL Id: inserm-00870153 https://www.hal.inserm.fr/inserm-00870153 Submitted on 5 Oct 2013 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Méline et al. International Journal of Health Geographics 2013, 12:44 INTERNATIONAL JOURNAL http://www.ij-healthgeographics.com/content/12/1/44 OF HEALTH GEOGRAPHICS RESEARCH Open Access Transportation noise and annoyance related to road traffic in the French RECORD study Julie Méline1,2*, Andraea Van Hulst3,4, Frédérique Thomas5, Noëlla Karusisi1,2 and Basile Chaix1,2 Abstract Road traffic and related noise is a major source of annoyance and impairment to health in urban areas. Many areas exposed to road traffic noise are also exposed to rail and air traffic noise.
    [Show full text]
  • Median Filtering by Threshold Decomposition: Induction Proof
    Median Filtering by Threshold Decomposition: Induction Proof by Connor Bramham, Brandon Mayle, Maura Gallagher, Douglas Magruder, and Cooper Reep Math Club Socrates Preparatory School www.socprep.org [email protected] Sponsors: Saul Rubin and Neal C. Gallagher, Ph.D, Socrates Preparatory School, 3955 Red Bug Lake Rd., Casselberry, FL 32707 [email protected] [email protected] Abstract In building a robot for the FTC competition, our team needed to remove motor noise from our sensor signals. So we settled on using a median filter because of the medians superior removal of impulsive noise. For us, however, the foundational publications that describe these filters were challenging to understand. Having learned the concept of proof by induction from the MIT OpenCourseWare course, “Mathematics for Computer Science” (MIT Course Number 6.042J / 18.062J), we developed an original proof for the principle of median filter threshold decomposition in order to better understand their operation. The induction is over the number of quantized threshold levels for the sequence of input values as applied to both the standard and recursive median filter. I Introduction Many of this paper’s authors participated in the 2017 First Tech Challenge robotics competition. Our completed robot contained many motors and sensors, including ultrasound sensors for distance measurements, gyroscopic sensors for direction measurements, color sensors to sense the brightness and color of beacons, and more. The robot’s electric motors generated a considerable amount of noise that found its way onto the signals from our sensors to our robot controller. One method used to address the noise problem is to use an aluminum sheet to shield robot electronics from the drive motors and servos, but more was needed.
    [Show full text]