Digital Image Processing

Size: px

Start display at page:

Download "Digital Image Processing"

Alan Hall
6 years ago
Views:

1 Digital Image Processing Using MATLAB Rafael C. Gonzalez University of Tennessee Richard E. Woods MedData Interactive Steven L. Eddins The MathWorks, Inc. Upper Saddle River, NJ 07458

2 Library of Congress Cataloging-in-Publication Data on File Vice President and Editorial Director, ECS: Marcia Horton Vice President and Director of Production and Manufacturing, ESM: David W. Riccardi Publisher: Tom Robbins Editorial Assistant: Carole Snyder Executive Managing Editor: Vince O Brien Managing Editor: David A. George Production Editor: Rose Kernan Director of Creative Services: Paul Belfanti Creative Director: Carole Anson Art Director: Jayne Conte Cover Designer: Richard E. Woods Art Editor: Xiaohong Zhu Manufacturing Manager: Trudy Pisciotti Manufacturing Buyer: Lisa McDowell Senior Marketing Manager: Holly Stark 2004 by Pearson Education, Inc. Pearson Prentice-Hall Pearson Education, Inc. Upper Saddle River, New Jersey All rights reserved. No part of this book may be reproduced or transmitted in any form or by any means, without permission in writing from the publisher. Pearson Prentice Hall is a trademark of Pearson Education, Inc. MATLAB is a registered trademark of The MathWorks, Inc., 3 Apple Hill Drive, Natick, MA The author and publisher of this book have used their best efforts in preparing this book. These efforts include the development, research, and testing of the theories and programs to determine their effectiveness. The author and publisher shall not be liable in any event for incidental or consequential damages with, or arising out of, the furnishing, performance, or use of these programs. Printed in the United States of America ISBN Pearson Education Ltd., London Pearson Education Australia Pty., Ltd., Sydney Pearson Education Singapore, Pte. Ltd. Pearson Education North Asia Ltd., Hong Kong Pearson Education Canada, Inc., Toronto Pearson Education de Mexico, S.A. de C.V. Pearson Education Japan, Tokyo Pearson Education Malaysia, Pte. Ltd. Pearson Education, Inc., Upper Saddle River, New Jersey

3 Contents Preface xi Acknowledgments About the Authors xii xiii 1 Introduction 1 Preview Background What Is Digital Image Processing? Background on MATLAB and the Image Processing Toolbox Areas of Image Processing Covered in the Book The Book Web Site Notation The MATLAB Working Environment The MATLAB Desktop Using the MATLAB Editor to Create M-files Getting Help Saving and Retrieving a Work Session How References Are Organized in the Book 11 Summary 11 2 Fundamentals 12 Preview Digital Image Representation Coordinate Conventions Images as Matrices Reading Images Displaying Images Writing Images Data Classes Image Types Intensity Images Binary Images A Note on Terminology Converting between Data Classes and Image Types Converting between Data Classes Converting between Image Classes and Types Array Indexing Vector Indexing Matrix Indexing Selecting Array Dimensions 37 v

4 vi Contents 2.9 Some Important Standard Arrays Introduction to M-Function Programming M-Files Operators Flow Control Code Optimization Interactive I/O A Brief Introduction to Cell Arrays and Structures 62 Summary 64 3 Intensity Transformations and Spatial Filtering 65 Preview Background Intensity Transformation Functions Function imadjust Logarithmic and Contrast-Stretching Transformations Some Utility M-Functions for Intensity Transformations Histogram Processing and Function Plotting Generating and Plotting Image Histograms Histogram Equalization Histogram Matching (Specification) Spatial Filtering Linear Spatial Filtering Nonlinear Spatial Filtering Image Processing Toolbox Standard Spatial Filters Linear Spatial Filters Nonlinear Spatial Filters 104 Summary Frequency Domain Processing 108 Preview The 2-D Discrete Fourier Transform Computing and Visualizing the 2-D DFT in MATLAB Filtering in the Frequency Domain Fundamental Concepts Basic Steps in DFT Filtering An M-function for Filtering in the Frequency Domain Obtaining Frequency Domain Filters from Spatial Filters Generating Filters Directly in the Frequency Domain Creating Meshgrid Arrays for Use in Implementing Filters in the Frequency Domain Lowpass Frequency Domain Filters Wireframe and Surface Plotting 132

5 Contents vii 4.6 Sharpening Frequency Domain Filters Basic Highpass Filtering High-Frequency Emphasis Filtering 138 Summary Image Restoration 141 Preview A Model of the Image Degradation/Restoration Process Noise Models Adding Noise with Function imnoise Generating Spatial Random Noise with a Specified Distribution Periodic Noise Estimating Noise Parameters Restoration in the Presence of Noise Only Spatial Filtering Spatial Noise Filters Adaptive Spatial Filters Periodic Noise Reduction by Frequency Domain Filtering Modeling the Degradation Function Direct Inverse Filtering Wiener Filtering Constrained Least Squares (Regularized) Filtering Iterative Nonlinear Restoration Using the Lucy-Richardson Algorithm Blind Deconvolution Geometric Transformations and Image Registration Geometric Spatial Transformations Applying Spatial Transformations to Images Image Registration Summary 193 Color Image Processing 194 Preview Color Image Representation in MATLAB RGB Images Indexed Images IPT Functions for Manipulating RGB and Indexed Images Converting to Other Color Spaces NTSC Color Space The YCbCr Color Space The HSV Color Space The CMY and CMYK Color Spaces The HSI Color Space The Basics of Color Image Processing Color Transformations 216

6 viii Contents 6.5 Spatial Filtering of Color Images Color Image Smoothing Color Image Sharpening Working Directly in RGB Vector Space Color Edge Detection Using the Gradient Image Segmentation in RGB Vector Space Summary 241 Wavelets 242 Preview Background The Fast Wavelet Transform FWTs Using the Wavelet Toolbox FWTs without the Wavelet Toolbox Working with Wavelet Decomposition Structures Editing Wavelet Decomposition Coefficients without the Wavelet Toolbox Displaying Wavelet Decomposition Coefficients The Inverse Fast Wavelet Transform Wavelets in Image Processing Summary 281 Image Compression 282 Preview Background Coding Redundancy Huffman Codes Huffman Encoding Huffman Decoding Interpixel Redundancy Psychovisual Redundancy JPEG Compression JPEG JPEG Summary Morphological Image Processing 334 Preview Preliminaries Some Basic Concepts from Set Theory Binary Images, Sets, and Logical Operators Dilation and Erosion Dilation Structuring Element Decomposition The strel Function Erosion 345

7 Contents ix 9.3 Combining Dilation and Erosion Opening and Closing The Hit-or-Miss Transformation Using Lookup Tables Function bwmorph Labeling Connected Components Morphological Reconstruction Opening by Reconstruction Filling Holes Clearing Border Objects Gray-Scale Morphology Dilation and Erosion Opening and Closing Reconstruction 374 Summary Image Segmentation 378 Preview Point, Line, and Edge Detection Point Detection Line Detection Edge Detection Using Function edge Line Detection Using the Hough Transform Hough Transform Peak Detection Hough Transform Line Detection and Linking Thresholding Global Thresholding Local Thresholding Region-Based Segmentation Basic Formulation Region Growing Region Splitting and Merging Segmentation Using the Watershed Transform Watershed Segmentation Using the Distance Transform Watershed Segmentation Using Gradients Marker-Controlled Watershed Segmentation 422 Summary Representation and Description 426 Preview Background Cell Arrays and Structures Some Additional MATLAB and IPT Functions Used in This Chapter Some Basic Utility M-Functions 433

8 x Contents 11.2 Representation Chain Codes Polygonal Approximations Using Minimum-Perimeter Polygons Signatures Boundary Segments Skeletons Boundary Descriptors Some Simple Descriptors Shape Numbers Fourier Descriptors Statistical Moments Regional Descriptors Function regionprops Texture Moment Invariants Using Principal Components for Description 474 Summary Object Recognition 484 Preview Background Computing Distance Measures in MATLAB Recognition Based on Decision-Theoretic Methods Forming Pattern Vectors Pattern Matching Using Minimum-Distance Classifiers Matching by Correlation Optimum Statistical Classifiers Adaptive Learning Systems Structural Recognition Working with Strings in MATLAB String Matching 508 Summary 513 Appendix A Function Summary 514 Appendix B ICE and MATLAB Graphical User Interfaces 527 Appendix C M-Functions 552 Bibliography 594 Index 597

Digital Image Processing

Digital Image Processing Third Edition Rafael C. Gonzalez University of Tennessee Richard E. Woods MedData Interactive PEARSON Prentice Hall Pearson Education International Contents Preface xv Acknowledgments