Image Enhancement by Edge-Preserving Filtering

Similar documents
Measurements of MeV Photon Flashes in Petawatt Laser Experiments

Argonne National Laboratory P.O. Box 2528 Idaho Falls, ID

UCRL-ID Broad-Band Characterization of the Complex Permittivity and Permeability of Materials. Carlos A. Avalle

The ACT External HEPA Push-Through Filter Assembly. A. A. Frigo, S. G. Wiedmeyer, D. E. Preuss, E. F. Bielick, and R. F. Malecha

Preparation of Random Phase Plates for Laser Beam Smoothing

Performance of Smoothing by Spectral Dispersion (SSD) with Frequency Conversion on the Beamlet Laser for the National Ignition Facility

Up-conversion Time Microscope Demonstrates 103x Magnification of an Ultrafast Waveforms with 300 fs Resolution. C. V. Bennett B. H.

Cascaded Wavelength Division Multiplexing for Byte-Wide Optical Interconnects

THE MEASURED PERFORMANCE OF A 170 GHz REMOTE STEERING LAUNCHER

Performance of Keck Adaptive Optics with Sodium Laser Guide Stars

STP-NU ROADMAP TO DEVELOP ASME CODE RULES FOR THE CONSTRUCTION OF HIGH TEMPERATURE GAS COOLED REACTORS (HTGRS)

Implementation of an Acoustic Emission Proximity Detector for Use in Generating Glass Optics. M. A. Piscotty, J. S. Taylor, K. L.

National Accelerator LaboratoryFERMILAB-TM-1966

High-]FrequencyElectric Field Measurement Using a Toroidal Antenna

Five-beam Fabry-Perot velocimeter

Microsecond-long Lasing Delays in Thin P-clad InGaAs QW Lasers

GA A22897 QUASI-OPTIC COMPONENTS IN OVERSIZED CORRUGATED WAVEGUIDE FOR MILLIMETER-WAVE TRANSMISSION SYSTEMS

Report on Ghosting in LL94 RAR Data

Accelerator and Fusion Research Division Lawrence Berkeley Laboratory University of California Berkeley, CA 94720

Recent advances in ALAMO

Lawrence Berkeley National Laboratory Recent Work

AIGaAs/InGaAIP Tunnel Junctions for Multifunction Solar Cells. Sharps, N. Y. Li, J. S. Hills, and H. Hou EMCORE Photovoltaics

GA A22574 ADVANTAGES OF TRAVELING WAVE RESONANT ANTENNAS FOR FAST WAVE HEATING SYSTEMS

MAPPING INDUCED POLARIZATION WITH NATURAL ELECTROMAGNETIC FIELDS FOR EXPLORATION AND RESOURCES CHARACTERIZATION BY THE MINING INDUSTRY

PEP-I11Magnet Power Conversion Systems:.

GA A23983 AN ADVANCED COLLABORATIVE ENVIRONMENT TO ENHANCE MAGNETIC FUSION RESEARCH

Nanosecond, pulsed, frequency-modulated optical parametric oscillator

STP-PT-032 BUCKLING OF CYLINDRICAL, THIN WALL, TRAILER TRUCK TANKS AND ASME SECTION XII

GA A FABRICATION OF A 35 GHz WAVEGUIDE TWT CIRCUIT USING RAPID PROTOTYPE TECHNIQUES by J.P. ANDERSON, R. OUEDRAOGO, and D.

National Accelerator Laboratory

k SLAC-PUB-7583 July 1997 Co/vF PULSE TRANSFORMER R&D FOR NLC KLYSTRON PULSE MODULATOR*

Los A LA-UR Los Alamos National Laboratory Los Alamos, New Mexico 87545

Introduction to Radar Systems. Radar Antennas. MIT Lincoln Laboratory. Radar Antennas - 1 PRH 6/18/02

cycle to cycle, so errors can be used to update the reference waveforms for future cycles. At A P S, updates are

Reducing space charge tune shift with a barrier cavity

MAPPING INDUCED POLARIZATION WITH NATURAL ELECTROMAGNETIC FIELDS FOR EXPLORATION AND RESOURCES CHARACTERIZATION BY THE MINING INDUSTRY

DISCLAIMER. Portions of this document may be illegible in electronic image products. Images are produced from the best available original document.

Quarterly Technical Progress Report. For the period starting July, , ending September 30, Xiaodi Huang and Richard Gertsch

STP-PT-054 CONCENTRATED SOLAR POWER (CSP) CODES AND STANDARDS GAP ANALYSIS

Detection of Targets in Noise and Pulse Compression Techniques

GA A23741 DATA MANAGEMENT, CODE DEPLOYMENT, AND SCIENTIFIC VISUALIZATION TO ENHANCE SCIENTIFIC DISCOVERY THROUGH ADVANCED COMPUTING

IMPACT TESTING EXEMPTION CURVES

GA A25836 PRE-IONIZATION EXPERIMENTS IN THE DIII-D TOKAMAK USING X-MODE SECOND HARMONIC ELECTRON CYCLOTRON HEATING

Performance of Image Intensifiers in Radiographic Systems

&wf-9+/ob/--21*~~ II. Ron Harper and Robert A. Hike

The Development of an Enhanced Strain Measurement Device to Support Testing of Radioactive Material Packages*

J.C. Courtney Nuclear Science Center Louisiana State University Baton Rouge, LA

GA A22712 DIII D ICRF HIGH VOLTAGE POWER SUPPLY REGULATOR UPGRADE

Controlling Changes Lessons Learned from Waste Management Facilities 8

A REGULATED POWER SUPPLY FOR THE FILAMENTS OF A HIGH POWER GYROTRON

CORRECTED RMS ERROR AND EFFECTIVE NUMBER OF BITS FOR SINEWAVE ADC TESTS

Sub-nanometer Interferometry Aspheric Mirror Fabrication

Y f OAK RIDGE Y4 2 PLANT. Lionel Levinson General Electric Company. November 24, Approved for Public Release; distribution is unlimited.

Los Alamos LA-UR-96- High Performance Computing for Domestic Petroleum Reservoir Simulation

SHADOWGRAPH ILLUMINIATION TECHNIQUES FOR FRAMING CAMERAS

DEVELOPMENT OF THE PULSE TRANSFORMER FOR NLC KLYSTRON PULSE MODULATOR* Abstract

Laser Surface Profiler

Spatial Frequency Domain Error Budget. Debbie Krulewich and Herman Hauschildt

GA A27238 MEASUREMENT OF DEUTERIUM ION TOROIDAL ROTATION AND COMPARISON TO NEOCLASSICAL THEORY IN THE DIII-D TOKAMAK

IMU integration into Sensor suite for Inspection of H-Canyon

GA A23281 EXTENDING DIII D NEUTRAL BEAM MODULATED OPERATIONS WITH A CAMAC BASED TOTAL ON TIME INTERLOCK

GA A26816 DESIGNS OF NEW COMPONENTS FOR ITER ECH&CD TRANSMISSION LINES

Assisting DOE EM 4.12, Office of Groundwater and Subsurface Closure

v-~ -8 m w Abstract Framework for Sandia Technolow Transfer Process Introduction

U.S. Air Force Phillips hboratoq, Kirtland AFB, NM 87117, 505/ , FAX:

AN ELM=RESlLlENT RF ARC DETECTION SYSTEM FOR DIII-D BASED ON ELECTROMAGNETIC AND SOUND EMISSIONS FROM THE ARC

Evaluation of Roof Bolting Requirements Based on In-Mine Roof Bolter Drilling

GA A22577 AN ELM-RESILIENT RF ARC DETECTION SYSTEM FOR DIII D BASED ON ELECTROMAGNETIC AND SOUND EMISSIONS FROM THE ARC

Measurements of edge density profile modifications during IBW on TFTR

+o GENEML ATOMfCS. RF POWER DIAGNOSTICS AND CONTROL ON THE DIII-D, 4 MW MHz FAST WAVE CURRENT DRIVE SYSTEM (FWCD)

INFRARED MEASUREMENTS OF THE SYNTHETIC DIAMOND WINDOW OF A 110 GHz HIGH POWER GYROTRON

DISCLAIMER. Portions of this document may be illegible in electronic image products. Images are produced from the best available original document.

1997 Particle Accelerator Conference, Vancouver, B.C., Canada, May 12-16, 1997 BNL

U.S. DEPARTMENT OF ENERGY. YlAMT-485 Y-I 2. Project Accomplishment Summary for Project Number 92-Y12P-013-B2 HYDROFORMING DESIGN AND PROCESS ADVISOR


Stimulated Emission from Semiconductor Microcavities

High Explosive Radio Telemetry System. Federal Manufacturing & Technologies. R. Johnson, FM&T; B. Mclaughlin, FM&T;

RECEIVED OCT % QSTI. v, J. Blasy, M. Rowland, V. Teryohin, G. Tzyg M. O'Brien, R. Wallace, M. Ehinger, L. Fishbone, T. James, J.

D&D Knowledge Management through Contributions in Wikipedia

Optical Control, Diagnostic and Power Supply System for a Solid State Induction Modulator

GA A22869 BOUNCE COATING INDUCED DOMES ON GLOW DISCHARGE POLYMER COATED SHELLS

Observation of amplification of a 1ps pulse by SRS of a 1 ns pulse in a plasma with conditions relevant to pulse compression

High Dynamic Range Imaging

MASTER --3. Gtl.- DISTRIBUTION. THiS DOCUMENT IS UNLIMITED PNL-SA Shaw Whiteman Anderson Alzheimer G. A. March 1995

AN IN-LINE POWER MONITOR FOR HE11 LOW LOSS TRANSMISSION LINES

CPC s Primary Objectives

Integration of MGDS Design into the Licensing Process' This paper presents an overview of how the Mined Geologic Disposal System (MGDS) design

Armlication For United States Patent For HOT CELL SHIELD PLUG EXTRACTION APPARATUS. Philip A. Knapp Moore, ID. and. Larry K. Manhart Pingree, ID

GA A22583 FAST WAVE ANTENNA ARRAY FEED CIRCUITS TOLERANT OF TIME-VARYING LOADING FOR DIII D

INITIAL RESULTS FROM THE MULTI-MEGAWATT 110 GHz ECH SYSTEM FOR THE DIII D TOKAMAK

Tip-Tilt Correction for Astronomical Telescopes using Adaptive Control. Jim Watson

DESIGNING MICROELECTROMECHANICAL SYSTEMS-ON-A-CHIP IN A 5-LEVEL SURF ACE MICROMACHINE TECHNOLOGY

Quarterly Technical Progress Report October 15,1998

Higher Efficiency Active Quenching Circuit for Avalanche Photodiodes

UTag: Long-range Ultra-wideband Passive Radio Frequency Tags

Introduction to Radar Systems. The Radar Equation. MIT Lincoln Laboratory _P_1Y.ppt ODonnell

DISCLAIMER. Portions of this document may be illegible in electronic image products. Images are produced from the best available original document.

GA A SOLID-STATE HIGH VOLTAGE MODULATOR WITH OUTPUT CONTROL UTILIZING SERIES-CONNECTED IGBTs by J.F. TOOKER and P. HUYNH

Parasitic Pencil Beams Caused by Lens Reflections in Laser Amplifier Chains

A Multilevel Voltage-Source Converter System with Balanced DC Voltages' Abstract

Transcription:

UCRL-JC-116695 PREPRINT Image Enhancement by Edge-Preserving Filtering Yiu-fai Wong This paper was prepared for submittal to the First IEEE International Conference on Image Processing Austin, TX November 13-26,1994 November 1994 I I Thisisapreprintofapaperintendedforpublicationina joumalorproceedings. Since changes may be made before publication, this preprint is made available with the understanding that it will not be cited or reproduced without the permission of the author. \

DISCLAIMER This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, make any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or. any agency thereof.

DISCLAIMER Portions of this document may be illegible in electronic image products. images are produced from the best available original document.

IMAGE ENHANCEMENT BY EDGE-PRESERWG FILTERING Yiu-fai Wong Institute for Scientific Computing Research, L416 Lawrence Livermore National Laboratory Livermore, CA 94550 ABSTRACT Image enhancement is useful when the details in an image are lost due to various reasons. It is common to subtract a mask from a given image to enhance the details. The trick is how to obtain a good mask. We describe here how an edge-preserving filter can be used to generate a mask which is smooth over areas with fine details, yet preserving most of the edges. Experiments with real images show that our scheme is very effective. 1. INTRODUCTION Image enhancement is needed when the dynamic range of a recording/display media is smaller than that of the signal it is exposed to. The media can be film, monitors and printers. Without any postprocessing, the visibility of important details in both the light and dark areas is greatly reduced. Such areas may appear flat or washed out (see [l]for an excellent review). This can also occur in other imagingmodalitiessuch as Magnetic Resonance Imaging. A common procedure for enhancement is to subtract a low-passed version of an image from the original, increasing the high-frequency components. But it is known that unpleasant effects, overshoots which look like bright and dark bands near the edges, can be created. Thus, sharp niasks are needed. An expensive and time-consuming way to obtain good quality pictures is to have photographers under-expose or over-expose the scene [2]. In this work, we propose and demonstrate the use of an edge-preserving filter called clusiering filter for image enhancement The filter will be briefly described below and the details can be found in [3, 41. Let xi be the coordinate of z? ~pixel2 and yi its gray Work was supported by Lawrence Livermore National Laboratory through DOE contract No. W-7405-ENG-48. We are parlly motivated by the work in [SI in w v h i d ~a saturating resistive grid was used to generate a mask. 2xi = ( i, j ) for an image with rectangular grids. We use tlris notation for simplicity. level. The output y at pixel location x is given by ~ ~is~the. scale in the input where wi = e - a ~ ~ x i - x a space and governs the effective neighborhood of data used in the filter. Given a and the data, p is computed as follows: Let g = y i w i / C i wi and u2 = xi(yj Select g ) ) " w i / Z.wi. One then chooses p = (2uY)-l. % yo = g; iterate equation (1)afew times. The converged y is the filtered output. As demonstrated in [3,4], the clustering filter is capable of 1) preserving edges, 2) removing impulsive noise and 3) providing improved smoothing of nonimpulsive noise. xi 2. ENHANCEMENT SCHEME AND RFSULTS A good mask should be smooth in areas with fine de- tails, yet has sharp edges between the larger-scale structures. Subtracting this mask from the original image and rescaling will enhance the details while avoiding the halos around the edges. This is the essence of the scheme below using the clustering filter: 1. Filter a given image I using a = 1/2 recursively k times to obtain 4. k = 5 in our case. 2. Compute the difference signal Id = I - Ii. 3. Compute the local mean M and variance V for every pixel of Id. ~ 2.5V(x) Ii(x) if IId(X) - M ( x ) < 4. I;n(x) = I ( x ) otherwise. 5. Compute I, = I - sim; s = 0.5 in our experiments. 6. Rescale Io appropriately. For example, let m and v be the mean and variance of I,. Io is stretched to the range between m - 2.571 and m + 2.511. { Let us explain each of the steps above. Step 1) is used to generate a smoothed version of the image to

serve as an initial mask. The mask is such that local details are smoothed while edges are preserved. Iiowever, smoothing means that not all useful signals can be preserved. For example, bright/dark spots can be removed completely by the filter; edges of smaller scale can be smoothed out; corners of edges can be rounded to some extent. Thus, the difference image computed in step 2) contains the signal which needs to be restored in the mask. To decide which difference signal should be added to the initial mask Ii, we compute the local means and variances of la in step 3). At corners and bright spots, the difference signal will be large a m pared to their local variances. Step 4) restores these signals to the mask by thresholding. After subtracting the mask from the original image in step 5), rescaling in step 6) brings the image to its full dynamic range. Figures la, lb and ICshow the Cameraman, the mask Imand the enhanced image Io respectively. One can see that the coat and the ground (grass) are enhanced while the brightness relations of the local areas are preserved, which is really important in image enhancement [l]. The scenes on the background also appear sharper. To demonstrate the use of our procedure for dynamic range compression, Figure 2a shows an original 8-bit MR image. Figure 2b shows the result when only 4 bits are used, which is obtained by quantizing I, after step 6). It is seen that the anatomical structures become much more visible. In both cases, the enhanced images are visibly more pleasant than those obtained by histogram equalization or adaptive histogram equalization, which are not shown here. The experiments demonstrate that the clustering filter is suitable for dynamic range compression and image enhancement. 3. SUMMARY We argue that edge-preserving filtering is essential for generating a mask such that halo effects can be avoided in image enhancement. We further demonstrate that the clustering filter can do a very good job of generating a good mask due to its nice smoothing and edgepreserving properties. Results on real images show that our scheme is very effective. 4. ACKNOWLEDGEMENT I am indebted to the late Edward Posner for his encouragements. I also acknowledge assistance from Thomas Leung who was supported by the Summer Undergraduate Research Fellowship Program at Caltech. Discussions with Andy Moore have been helpful. 5. R E F E R E N C E S [l]w.f. Schreiber, "ImageProcessing for Quality Improvement," IEEE Proc., Vol. 66, No. 12, 1978. [2] J.A.C. Yule, Principles of Color Reproduction. New York, Wiley, 1967. [3] Yiu-fai Wong, "Saddlenode Dynamics for Edgepreserving and Scale-space Filtering," Proc. IEEE Intl. Con$ Image Processing. November, 1994. [4] Yiu-fai Wong, "Nonlinear Scale-space Filtering and Multiresolution System," to appear in IEEE Bans. Image Processing. [5] A.J. Moore, Spatial Filtering in Tdne Reproduction and Vision, PhD Thesis, California Institute of Technology, 1992.

d M 2 W a