COGS 119/219 MATLAB for Experimental Research. Fall 2017 Image Processing in Matlab

Transcription

1 COGS 119/219 MATLAB for Experimental Research Fall 2017 Image Processing in Matlab

2 What is an image? An image is an array, or a matrix of square pixels (picture elements) arranged in rows and columns.

3 8-bit grey scale image In a (8-bit) grey scale image, each picture element has an assigned intensity that ranges from 0 to 255. A grey scale image includes many shades of grey. Each pixel has a value from 0 (black) to 255 (white). The possible range of the pixel values depend on the color depth of the image, here 8-bit = 256 tones or grey scales.

4 True-color A true-color image has 24 bit color depth ~= 16 million colors. A true color image assembled from three grey scale images colored red, green and blue.

5 16 bit grey scale image Some grey scale images have more grey scales, e.g. 16 bits = grey scales. In principle, three grey scale images can be combined to form an image with 281,474,976,710,656 grey scales ( = x x 65536)

6 Image formats There are two general groups of images: vector graphics (or line art) and bitmaps (pixel-based). Some of the most common file formats are:

7 Colors - RGB The RGB color model relates very closely to the way we perceive color with the r, g and b receptors in or retinas. RGB uses additive color mixing. RGB is the basic color model used in television, computers, for web graphics, and in any medium that projects color with light. The secondary colors of RGB cyan, magenta and yellow are formed by mixing the two of the primary colors (red, green or blue) and excluding the third color.

8 RGB Red and green combine to make yellow, green and blue to make cyan, blue and red form magenta. The combination of red, green and blue in full intensity makes white.

9 CMYK The 4-color CMYK model used in printing lays down overlapping layers of varying percentages of transparent cyan (C), magenta (M), and yellow (Y) inks. In addition, a layer of black (B) ink can be added. CMYK color model uses subtractive color model.

10 Figures and Images >> help figure figure(): creates a figure window figure(h): makes H the current figure, forces it to become visible, and raises it above all other figures on the screen. If Figure H does not exist, and H is an integer, a new figure is created with handle H. >> help image image(x): displays matrix X as an image. Each element of X specifies the color of a rectilinear patch in the image.

11 Figures and Images >> x = [ ]; >> figure(1) >> image(x) >> x = [ ; ]; >> figure(2) >> image(x) How do the colors refer to numbers?

12 Indexed images in Matlab An indexed image consists of a data matrix X and a colormap matrix, which is N x 3 array of class double containing values in the range [0,1]. Each row of the map specifies red, blue and green components of a single color. The color of each image pixel is determined by using the corresponding value of X as an index into the map. The value 1 points to the first row in the map, 2 points to the second row, and so on. Values of X must be integers.

13 Colormap You can use defaults or make your own. >> help colormap colormap(map) sets the current figure's colormap to MAP. Some built-in colormaps are: Jet: Default colormap

14 Colormap >> figure(3); >> image(x); >> colormap( autumn ) >> figure(4); >> image(x); >> colormap( spring )

15 Colormap A colormap is also a matrix. For color images, it is an N x 3 matrix where the three columns corresponds to the RGB (red, green, blue) indices.

16 Colormap >> c = colormap( spring ); >> c2 = colormap( autumn );

17 Colormap The reason we don t see much of a color range is, the color map has many possible colors (e.g. 64 colors) so we are using only a subset of the colors. But we can easily scale the values of the data matrix to use the full color map. >> help imagesc Scale data and display as image. >> figure(5) >> imagesc(x) >> colormap( spring )

18 Colormap Alternatively, we can create a new image matrix that uses a large range of values (between 0 and 64) in the first place. >> img = round(rand(100,100)*63)+1; >> figure(6) >> image(img) >> colormap( spring );

19 Colormap >> figure(3) >> colormap( hot ) >> colormap(3, winter ) will change the figure whose handle is 3

20 Change axis elements >> help axis >> figure(3); >> axis off; >> axis square;

21 User-defined colormaps We can also make our own colormaps. >> mycmap1 = [0 0 0 ; ; ; ; 1 1 1]; >> mycmap1 = reshape(mycmap1, 3,5 )'; >> figure(7); >> axis off; >> axis square; >> image(img); >> colormap(mycmap1);

22 User-defined colormaps We can also make our own colormaps. >> mycmap1 = [0 0 0 ; ; ; ; 1 1 1]; >> mycmap1 = reshape(mycmap1, 3,5 )'; >> figure(7); >> axis off; >> axis square; >> image(img); >> colormap(mycmap1); Why aren t there more colors?

23 User-defined colormaps We can also make our own colormaps. >> mycmap1 = [0 0 0 ; ; ; ; 1 1 1]; >> mycmap1 = reshape(mycmap1, 3,5 )'; >> figure(7); >> axis off; >> axis square; >> image(img); >> colormap(mycmap1); Now the issue is my colormap has 5 colors but my image has 64. Since I have haven t used imagesc, any number greater than or equal to 5in this image will be displayed as the 5 th color, and in this example that is a medium-light blue.try changing the 5 th row of mycmap1 and draw figure again to check

24 User-defined colormaps We can also make our own colormaps. >> mycmap1 = [0 0 0 ; ; ; ; 1 1 1]; >> mycmap1 = reshape(mycmap1, 3,5 )'; >> figure(7); >> axis off; >> axis square; >> image(img); >> colormap(mycmap1); Try imagesc instead of image

25 User-defined colormaps We can also make our own colormaps. >> mycmap1 = [0 0 0 ; ; ; ; 1 1 1]; >> mycmap1 = reshape(mycmap1, 3,5 )'; >> figure(7); >> axis off; >> axis square; >> imagesc(img); >> colormap(mycmap1); Notice colors in the image are the same but now it s interpolating such that the range of numbers 1:64 are being assigned more uniformly to the range of colors available

26 Let s play with colormaps New.m file called mycolormaps.m clear all; close all; colormap(gray(256)); myimg = reshape(1:256,16,16); image(myimg); axis square; axis off; pause figure(1); for i = 1:200 cmaps = rand(256,3); colormap(cmaps); drawnow end What does the code do?