- •Getting Started
- •Product Description
- •Key Features
- •Configuration Notes
- •Related Products
- •Compilability
- •Image Import and Export
- •Introduction
- •Step 1: Read and Display an Image
- •Step 2: Check How the Image Appears in the Workspace
- •Step 3: Improve Image Contrast
- •Step 4: Write the Image to a Disk File
- •Step 5: Check the Contents of the Newly Written File
- •Image Enhancement and Analysis
- •Introduction
- •Step 1: Read Image
- •Step 2: Use Morphological Opening to Estimate the Background
- •Step 3: View the Background Approximation as a Surface
- •Step 4: Subtract the Background Image from the Original Image
- •Step 5: Increase the Image Contrast
- •Step 6: Threshold the Image
- •Step 7: Identify Objects in the Image
- •Step 8: Examine One Object
- •Step 9: View All Objects
- •Step 10: Compute Area of Each Object
- •Step 11: Compute Area-based Statistics
- •Step 12: Create Histogram of the Area
- •Getting Help
- •Product Documentation
- •Image Processing Examples
- •MATLAB Newsgroup
- •Acknowledgments
- •Introduction
- •Images in MATLAB
- •Expressing Image Locations
- •Pixel Indices
- •Spatial Coordinates
- •Intrinsic Coordinates
- •World Coordinates
- •Image Types in the Toolbox
- •Overview of Image Types
- •Binary Images
- •Indexed Images
- •Grayscale Images
- •Truecolor Images
- •Converting Between Image Types
- •Converting Between Image Classes
- •Overview of Image Class Conversions
- •Losing Information in Conversions
- •Converting Indexed Images
- •Working with Image Sequences
- •Overview of Toolbox Functions That Work with Image Sequences
- •Process Image Sequences
- •Process Multi-Frame Image Arrays
- •Image Arithmetic
- •Overview of Image Arithmetic Functions
- •Image Arithmetic Saturation Rules
- •Nesting Calls to Image Arithmetic Functions
- •Getting Information About a Graphics File
- •Reading Image Data
- •Writing Image Data to a File
- •Overview
- •Specifying Format-Specific Parameters
- •Reading and Writing Binary Images in 1-Bit Format
- •Determining the Storage Class of the Output File
- •Converting Between Graphics File Formats
- •Working with DICOM Files
- •Overview of DICOM Support
- •Reading Metadata from a DICOM File
- •Handling Private Metadata
- •Creating Your Own Copy of the DICOM Dictionary
- •Reading Image Data from a DICOM File
- •Viewing Images from DICOM Files
- •Writing Image Data or Metadata to a DICOM File
- •Writing Metadata with the Image Data
- •Understanding Explicit Versus Implicit VR Attributes
- •Removing Confidential Information from a DICOM File
- •Example: Creating a New DICOM Series
- •Working with Mayo Analyze 7.5 Files
- •Working with Interfile Files
- •Working with High Dynamic Range Images
- •Understanding Dynamic Range
- •Reading a High Dynamic Range Image
- •Creating a High Dynamic Range Image
- •Viewing a High Dynamic Range Image
- •Writing a High Dynamic Range Image to a File
- •Image Display and Exploration Overview
- •Displaying Images Using the imshow Function
- •Overview
- •Specifying the Initial Image Magnification
- •Controlling the Appearance of the Figure
- •Displaying Each Image in a Separate Figure
- •Displaying Multiple Images in the Same Figure
- •Dividing a Figure Window into Multiple Display Regions
- •Using the subimage Function to Display Multiple Images
- •Using the Image Tool to Explore Images
- •Image Tool Overview
- •Opening the Image Tool
- •Specifying the Initial Image Magnification
- •Specifying the Colormap
- •Importing Image Data from the Workspace
- •Exporting Image Data to the Workspace
- •Using the getimage Function to Export Image Data
- •Saving the Image Data Displayed in the Image Tool
- •Closing the Image Tool
- •Printing the Image in the Image Tool
- •Exploring Very Large Images
- •Overview
- •Creating an R-Set File
- •Opening an R-Set File
- •Using Image Tool Navigation Aids
- •Navigating an Image Using the Overview Tool
- •Starting the Overview Tool
- •Moving the Detail Rectangle to Change the Image View
- •Specifying the Color of the Detail Rectangle
- •Getting the Position and Size of the Detail Rectangle
- •Printing the View of the Image in the Overview Tool
- •Panning the Image Displayed in the Image Tool
- •Zooming In and Out on an Image in the Image Tool
- •Specifying the Magnification of the Image
- •Getting Information about the Pixels in an Image
- •Determining the Value of Individual Pixels
- •Saving the Pixel Value and Location Information
- •Determining the Values of a Group of Pixels
- •Selecting a Region
- •Customizing the View
- •Determining the Location of the Pixel Region Rectangle
- •Printing the View of the Image in the Pixel Region Tool
- •Determining the Display Range of an Image
- •Measuring the Distance Between Two Pixels
- •Using the Distance Tool
- •Exporting Endpoint and Distance Data
- •Customizing the Appearance of the Distance Tool
- •Adjusting Image Contrast Using the Adjust Contrast Tool
- •Understanding Contrast Adjustment
- •Starting the Adjust Contrast Tool
- •Using the Histogram Window to Adjust Image Contrast
- •Using the Window/Level Tool to Adjust Image Contrast
- •Example: Adjusting Contrast with the Window/Level Tool
- •Modifying Image Data
- •Saving the Modified Image Data
- •Cropping an Image Using the Crop Image Tool
- •Viewing Image Sequences
- •Overview
- •Viewing Image Sequences in the Movie Player
- •Example: Viewing a Sequence of MRI Images
- •Configuring the Movie Player
- •Specifying the Frame Rate
- •Specifying the Color Map
- •Getting Information about the Image Frame
- •Viewing Image Sequences as a Montage
- •Converting a Multiframe Image to a Movie
- •Displaying Different Image Types
- •Displaying Indexed Images
- •Displaying Grayscale Images
- •Displaying Grayscale Images That Have Unconventional Ranges
- •Displaying Binary Images
- •Changing the Display Colors of a Binary Image
- •Displaying Truecolor Images
- •Adding a Colorbar to a Displayed Image
- •Printing Images
- •Printing and Handle Graphics Object Properties
- •Setting Toolbox Preferences
- •Retrieving the Values of Toolbox Preferences Programmatically
- •Setting the Values of Toolbox Preferences Programmatically
- •Overview
- •Displaying the Target Image
- •Creating the Modular Tools
- •Overview
- •Associating Modular Tools with a Particular Image
- •Getting the Handle of the Target Image
- •Specifying the Parent of a Modular Tool
- •Tools With Separate Creation Functions
- •Example: Embedding the Pixel Region Tool in an Existing Figure
- •Positioning the Modular Tools in a GUI
- •Specifying the Position with a Position Vector
- •Build a Pixel Information GUI
- •Adding Navigation Aids to a GUI
- •Understanding Scroll Panels
- •Example: Building a Navigation GUI for Large Images
- •Customizing Modular Tool Interactivity
- •Overview
- •Build Image Comparison Tool
- •Creating Your Own Modular Tools
- •Overview
- •Create Angle Measurement Tool
- •Spatial Transformations
- •Resizing an Image
- •Overview
- •Specifying the Interpolation Method
- •Preventing Aliasing by Using Filters
- •Rotating an Image
- •Cropping an Image
- •Perform General 2-D Spatial Transformations
- •Spatial Transformation Procedure
- •Translate Image Using maketform and imtransform
- •Step 1: Import the Image to Be Transformed
- •Step 2: Define the Spatial Transformation
- •Step 3: Create the TFORM Structure
- •Step 4: Perform the Transformation
- •Step 5: View the Output Image
- •Defining the Transformation Data
- •Using a Transformation Matrix
- •Using Sets of Points
- •Creating TFORM Structures
- •Performing the Spatial Transformation
- •Specifying Fill Values
- •Performing N-Dimensional Spatial Transformations
- •Register Image Using XData and YData Parameters
- •Step 1: Read in Base and Unregistered Images
- •Step 2: Display the Unregistered Image
- •Step 3: Create a TFORM Structure
- •Step 4: Transform the Unregistered Image
- •Step 5: Overlay Base Image Over Registered Image
- •Step 6: Using XData and YData Input Parameters
- •Step 7: Using xdata and ydata Output Values
- •Image Registration
- •Image Registration Techniques
- •Control Point Registration
- •Using cpselect in a Script
- •Example: Registering to a Digital Orthophoto
- •Step 1: Read the Images
- •Step 2: Choose Control Points in the Images
- •Step 3: Save the Control Point Pairs to the MATLAB Workspace
- •Step 4: Fine-Tune the Control Point Pair Placement (Optional)
- •Step 6: Transform the Unregistered Image
- •Geometric Transformation Types
- •Selecting Control Points
- •Specifying Control Points Using the Control Point Selection Tool
- •Starting the Control Point Selection Tool
- •Using Navigation Tools to Explore the Images
- •Using Scroll Bars to View Other Parts of an Image
- •Using the Detail Rectangle to Change the View
- •Panning the Image Displayed in the Detail Window
- •Zooming In and Out on an Image
- •Specifying the Magnification of the Images
- •Locking the Relative Magnification of the Input and Base Images
- •Specifying Matching Control Point Pairs
- •Picking Control Point Pairs Manually
- •Using Control Point Prediction
- •Moving Control Points
- •Deleting Control Points
- •Exporting Control Points to the Workspace
- •Saving Your Control Point Selection Session
- •Using Correlation to Improve Control Points
- •Intensity-Based Automatic Image Registration
- •Registering Multimodal MRI Images
- •Step 1: Load Images
- •Step 2: Set up the Initial Registration
- •Step 3: Improve the Registration
- •Step 4: Improve the Speed of Registration
- •Step 5: Further Refinement
- •Step 6: Deciding When Enough is Enough
- •Step 7: Alternate Visualizations
- •Designing and Implementing 2-D Linear Filters for Image Data
- •Overview
- •Convolution
- •Correlation
- •Performing Linear Filtering of Images Using imfilter
- •Data Types
- •Correlation and Convolution Options
- •Boundary Padding Options
- •Multidimensional Filtering
- •Relationship to Other Filtering Functions
- •Filtering an Image with Predefined Filter Types
- •Designing Linear Filters in the Frequency Domain
- •FIR Filters
- •Frequency Transformation Method
- •Frequency Sampling Method
- •Windowing Method
- •Creating the Desired Frequency Response Matrix
- •Computing the Frequency Response of a Filter
- •Transforms
- •Fourier Transform
- •Definition of Fourier Transform
- •Visualizing the Fourier Transform
- •Discrete Fourier Transform
- •Relationship to the Fourier Transform
- •Visualizing the Discrete Fourier Transform
- •Applications of the Fourier Transform
- •Frequency Response of Linear Filters
- •Fast Convolution
- •Locating Image Features
- •Discrete Cosine Transform
- •DCT Definition
- •The DCT Transform Matrix
- •DCT and Image Compression
- •Radon Transform
- •Radon Transformation Definition
- •Plotting the Radon Transform
- •Viewing the Radon Transform as an Image
- •Detecting Lines Using the Radon Transform
- •The Inverse Radon Transformation
- •Inverse Radon Transform Definition
- •Improving the Results
- •Reconstruct Image from Parallel Projection Data
- •Fan-Beam Projection Data
- •Fan-Beam Projection Data Definition
- •Computing Fan-Beam Projection Data
- •Image Reconstruction Using Fan-Beam Projection Data
- •Reconstruct Image From Fanbeam Projections
- •Morphological Operations
- •Morphology Fundamentals: Dilation and Erosion
- •Understanding Dilation and Erosion
- •Processing Pixels at Image Borders (Padding Behavior)
- •Understanding Structuring Elements
- •The Origin of a Structuring Element
- •Creating a Structuring Element
- •Structuring Element Decomposition
- •Dilating an Image
- •Eroding an Image
- •Combining Dilation and Erosion
- •Morphological Opening
- •Skeletonization
- •Perimeter Determination
- •Morphological Reconstruction
- •Understanding Morphological Reconstruction
- •Understanding the Marker and Mask
- •Pixel Connectivity
- •Defining Connectivity in an Image
- •Choosing a Connectivity
- •Specifying Custom Connectivities
- •Flood-Fill Operations
- •Specifying Connectivity
- •Specifying the Starting Point
- •Filling Holes
- •Finding Peaks and Valleys
- •Terminology
- •Understanding the Maxima and Minima Functions
- •Finding Areas of High or Low Intensity
- •Suppressing Minima and Maxima
- •Imposing a Minimum
- •Creating a Marker Image
- •Applying the Marker Image to the Mask
- •Distance Transform
- •Labeling and Measuring Objects in a Binary Image
- •Understanding Connected-Component Labeling
- •Remarks
- •Selecting Objects in a Binary Image
- •Finding the Area of the Foreground of a Binary Image
- •Finding the Euler Number of a Binary Image
- •Lookup Table Operations
- •Creating a Lookup Table
- •Using a Lookup Table
- •Getting Image Pixel Values Using impixel
- •Creating an Intensity Profile of an Image Using improfile
- •Displaying a Contour Plot of Image Data
- •Creating an Image Histogram Using imhist
- •Getting Summary Statistics About an Image
- •Computing Properties for Image Regions
- •Analyzing Images
- •Detecting Edges Using the edge Function
- •Detecting Corners Using the corner Function
- •Tracing Object Boundaries in an Image
- •Choosing the First Step and Direction for Boundary Tracing
- •Detecting Lines Using the Hough Transform
- •Analyzing Image Homogeneity Using Quadtree Decomposition
- •Example: Performing Quadtree Decomposition
- •Analyzing the Texture of an Image
- •Understanding Texture Analysis
- •Using Texture Filter Functions
- •Understanding the Texture Filter Functions
- •Example: Using the Texture Functions
- •Gray-Level Co-Occurrence Matrix (GLCM)
- •Create a Gray-Level Co-Occurrence Matrix
- •Specifying the Offsets
- •Derive Statistics from a GLCM and Plot Correlation
- •Adjusting Pixel Intensity Values
- •Understanding Intensity Adjustment
- •Adjusting Intensity Values to a Specified Range
- •Specifying the Adjustment Limits
- •Setting the Adjustment Limits Automatically
- •Gamma Correction
- •Adjusting Intensity Values Using Histogram Equalization
- •Enhancing Color Separation Using Decorrelation Stretching
- •Simple Decorrelation Stretching
- •Adding a Linear Contrast Stretch
- •Removing Noise from Images
- •Understanding Sources of Noise in Digital Images
- •Removing Noise By Linear Filtering
- •Removing Noise By Median Filtering
- •Removing Noise By Adaptive Filtering
- •ROI-Based Processing
- •Specifying a Region of Interest (ROI)
- •Overview of ROI Processing
- •Using Binary Images as a Mask
- •Creating a Binary Mask
- •Creating an ROI Without an Associated Image
- •Creating an ROI Based on Color Values
- •Filtering an ROI
- •Overview of ROI Filtering
- •Filtering a Region in an Image
- •Specifying the Filtering Operation
- •Filling an ROI
- •Image Deblurring
- •Understanding Deblurring
- •Causes of Blurring
- •Deblurring Model
- •Importance of the PSF
- •Deblurring Functions
- •Deblurring with the Wiener Filter
- •Refining the Result
- •Deblurring with a Regularized Filter
- •Refining the Result
- •Deblurring with the Lucy-Richardson Algorithm
- •Overview
- •Reducing the Effect of Noise Amplification
- •Accounting for Nonuniform Image Quality
- •Handling Camera Read-Out Noise
- •Handling Undersampled Images
- •Example: Using the deconvlucy Function to Deblur an Image
- •Refining the Result
- •Deblurring with the Blind Deconvolution Algorithm
- •Example: Using the deconvblind Function to Deblur an Image
- •Refining the Result
- •Creating Your Own Deblurring Functions
- •Avoiding Ringing in Deblurred Images
- •Color
- •Displaying Colors
- •Reducing the Number of Colors in an Image
- •Reducing Colors Using Color Approximation
- •Quantization
- •Colormap Mapping
- •Reducing Colors Using imapprox
- •Dithering
- •Converting Color Data Between Color Spaces
- •Understanding Color Spaces and Color Space Conversion
- •Converting Between Device-Independent Color Spaces
- •Supported Conversions
- •Example: Performing a Color Space Conversion
- •Color Space Data Encodings
- •Performing Profile-Based Color Space Conversions
- •Understanding Device Profiles
- •Reading ICC Profiles
- •Writing Profile Information to a File
- •Example: Performing a Profile-Based Conversion
- •Specifying the Rendering Intent
- •Converting Between Device-Dependent Color Spaces
- •YIQ Color Space
- •YCbCr Color Space
- •HSV Color Space
- •Neighborhood or Block Processing: An Overview
- •Performing Sliding Neighborhood Operations
- •Understanding Sliding Neighborhood Processing
- •Determining the Center Pixel
- •General Algorithm of Sliding Neighborhood Operations
- •Padding Borders in Sliding Neighborhood Operations
- •Performing Distinct Block Operations
- •Understanding Distinct Block Processing
- •Implementing Block Processing Using the blockproc Function
- •Applying Padding
- •Block Size and Performance
- •TIFF Image Characteristics
- •Choosing Block Size
- •Using Parallel Block Processing on large Image Files
- •What is Parallel Block Processing?
- •When to Use Parallel Block Processing
- •How to Use Parallel Block Processing
- •Working with Data in Unsupported Formats
- •Learning More About the LAN File Format
- •Parsing the Header
- •Reading the File
- •Examining the LanAdapter Class
- •Using the LanAdapter Class with blockproc
- •Understanding Columnwise Processing
- •Using Column Processing with Sliding Neighborhood Operations
- •Using Column Processing with Distinct Block Operations
- •Restrictions
- •Code Generation for Image Processing Toolbox Functions
- •Supported Functions
- •Examples
- •Introductory Examples
- •Image Sequences
- •Image Representation and Storage
- •Image Display and Visualization
- •Zooming and Panning Images
- •Pixel Values
- •Image Measurement
- •Image Enhancement
- •Brightness and Contrast Adjustment
- •Cropping Images
- •GUI Application Development
- •Edge Detection
- •Regions of Interest (ROI)
- •Resizing Images
- •Image Registration and Alignment
- •Image Filtering
- •Fourier Transform
- •Image Transforms
- •Feature Detection
- •Discrete Cosine Transform
- •Image Compression
- •Radon Transform
- •Image Reconstruction
- •Fan-beam Transform
- •Morphological Operations
- •Binary Images
- •Image Histogram
- •Image Analysis
- •Corner Detection
- •Hough Transform
- •Image Texture
- •Image Statistics
- •Color Adjustment
- •Noise Reduction
- •Filling Images
- •Deblurring Images
- •Image Color
- •Color Space Conversion
- •Block Processing
- •Index
- •Summary of Modular Tools
- •Rules for Dilation and Erosion
- •Rules for Padding Images
- •Supported Connectivities
- •Distance Metrics
- •File Header Content
14 Color
Reducing the Number of Colors in an Image
In this section...
“Reducing Colors Using Color Approximation” on page 14-4 “Reducing Colors Using imapprox” on page 14-10 “Dithering” on page 14-11
Reducing Colors Using Color Approximation
On systems with 24-bit color displays, truecolor images can display up to 16,777,216 (i.e., 224) colors. On systems with lower screen bit depths, truecolor images are still displayed reasonably well, because MATLAB automatically uses color approximation and dithering if needed. Color approximation is the process by which the software chooses replacement colors in the event that direct matches cannot be found.
Indexed images, however, might cause problems if they have a large number of colors. In general, you should limit indexed images to 256 colors for the following reasons:
•On systems with 8-bit display, indexed images with more than 256 colors will need to be dithered or mapped and, therefore, might not display well.
•On some platforms, colormaps cannot exceed 256 entries.
•If an indexed image has more than 256 colors, MATLAB cannot store the image data in a uint8 array, but generally uses an array of class double instead, making the storage size of the image much larger (each pixel uses 64 bits).
•Most image file formats limit indexed images to 256 colors. If you write an indexed image with more than 256 colors (using imwrite) to a format that does not support more than 256 colors, you will receive an error.
To reduce the number of colors in an image, use the rgb2ind function. This function converts a truecolor image to an indexed image, reducing the number of colors in the process. rgb2ind provides the following methods for approximating the colors in the original image:
14-4
Reducing the Number of Colors in an Image
•Quantization (described in “Quantization” on page 14-5)
-Uniform quantization
-Minimum variance quantization
•Colormap mapping (described in “Colormap Mapping” on page 14-9)
The quality of the resulting image depends on the approximation method you use, the range of colors in the input image, and whether or not you use dithering. Note that different methods work better for different images. See “Dithering” on page 14-11 for a description of dithering and how to enable or disable it.
Quantization
Reducing the number of colors in an image involves quantization. The function rgb2ind uses quantization as part of its color reduction algorithm. rgb2ind supports two quantization methods: uniform quantization and minimum variance quantization.
An important term in discussions of image quantization is RGB color cube, which is used frequently throughout this section. The RGB color cube is a three-dimensional array of all of the colors that are defined for a particular data type. Since RGB images in MATLAB can be of type uint8, uint16, or double, three possible color cube definitions exist. For example, if an RGB image is of class uint8, 256 values are defined for each color plane (red, blue, and green), and, in total, there will be 224 (or 16,777,216) colors defined by the color cube. This color cube is the same for all uint8 RGB images, regardless of which colors they actually use.
The uint8, uint16, and double color cubes all have the same range of colors. In other words, the brightest red in a uint8 RGB image appears the same as the brightest red in a double RGB image. The difference is that the double RGB color cube has many more shades of red (and many more shades of all colors). The following figure shows an RGB color cube for a uint8 image.
14-5
14 Color
RGB Color Cube for uint8 Images
Quantization involves dividing the RGB color cube into a number of smaller boxes, and then mapping all colors that fall within each box to the color value at the center of that box.
Uniform quantization and minimum variance quantization differ in the approach used to divide up the RGB color cube. With uniform quantization, the color cube is cut up into equal-sized boxes (smaller cubes). With minimum variance quantization, the color cube is cut up into boxes (not necessarily cubes) of different sizes; the sizes of the boxes depend on how the colors are distributed in the image.
Uniform Quantization. To perform uniform quantization, call rgb2ind and specify a tolerance. The tolerance determines the size of the cube-shaped boxes into which the RGB color cube is divided. The allowable range for a tolerance setting is [0,1]. For example, if you specify a tolerance of 0.1, the edges of the boxes are one-tenth the length of the RGB color cube and the maximum total number of boxes is
n = (floor(1/tol)+1)^3
The commands below perform uniform quantization with a tolerance of 0.1.
RGB = imread('peppers.png'); [x,map] = rgb2ind(RGB, 0.1);
14-6
Reducing the Number of Colors in an Image
The following figure illustrates uniform quantization of a uint8 image. For clarity, the figure shows a two-dimensional slice (or color plane) from the color cube where red=0 and green and blue range from 0 to 255. The actual pixel values are denoted by the centers of the x’s.
Uniform Quantization on a Slice of the RGB Color Cube
After the color cube has been divided, all empty boxes are thrown out. Therefore, only one of the boxes is used to produce a color for the colormap. As shown earlier, the maximum length of a colormap created by uniform quantization can be predicted, but the colormap can be smaller than the prediction because rgb2ind removes any colors that do not appear in the input image.
Minimum Variance Quantization. To perform minimum variance quantization, call rgb2ind and specify the maximum number of colors in the output image’s colormap. The number you specify determines the number of boxes into which the RGB color cube is divided. These commands use minimum variance quantization to create an indexed image with 185 colors.
RGB = imread('peppers.png'); [X,map] = rgb2ind(RGB,185);
14-7
14 Color
Minimum variance quantization works by associating pixels into groups based on the variance between their pixel values. For example, a set of blue pixels might be grouped together because they have a small variance from the center pixel of the group.
In minimum variance quantization, the boxes that divide the color cube vary in size, and do not necessarily fill the color cube. If some areas of the color cube do not have pixels, there are no boxes in these areas.
While you set the number of boxes, n, to be used by rgb2ind, the placement is determined by the algorithm as it analyzes the color data in your image. Once the image is divided into n optimally located boxes, the pixels within each box are mapped to the pixel value at the center of the box, as in uniform quantization.
The resulting colormap usually has the number of entries you specify. This is because the color cube is divided so that each region contains at least one color that appears in the input image. If the input image uses fewer colors than the number you specify, the output colormap will have fewer than n colors, and the output image will contain all of the colors of the input image.
The following figure shows the same two-dimensional slice of the color cube as shown in the preceding figure (demonstrating uniform quantization). Eleven boxes have been created using minimum variance quantization.
14-8
