l = dir('asl/*.jpg');
s = strcat('asl/',l(10).name)

s =

asl/ASL-A.jpg

I = imread(s);
Ir = I(:,:,1);


BWr = im2bw(Ir,graythresh(Ir));
figure, imshow(BWr);
help bwlabel
 BWLABEL Label connected components in binary image.
    L = BWLABEL(BW,N) returns a matrix L, of the same size as BW, containing
    labels for the connected components in BW. N can have a value of either
    4 or 8, where 4 specifies 4-connected objects and 8 specifies
    8-connected objects; if the argument is omitted, it defaults to 8.
 
    The elements of L are integer values greater than or equal to 0.  The
    pixels labeled 0 are the background.  The pixels labeled 1 make up one
    object, the pixels labeled 2 make up a second object, and so on.
 
    [L,NUM] = BWLABEL(BW,N) returns in NUM the number of connected objects
    found in BW.
 
    Note: Comparing BWLABEL and BWLABELN
    ------------------------------------
    BWLABEL supports 2-D inputs only, whereas BWLABELN support any 
    input dimension.  In some cases you might prefer to use BWLABELN even
    for 2-D problems because it can be faster.  If you have a 2-D input
    whose objects are relatively "thick" in the vertical direction,
    BWLABEL will probably be faster; otherwise BWLABELN will probably be
    faster.
 
    Class Support
    -------------
    BW can be logical or numeric, and it must be real, 2-D, and
    nonsparse.  L is double.
 
    Example
    -------
        BW = logical([1 1 1 0 0 0 0 0
                      1 1 1 0 1 1 0 0
                      1 1 1 0 1 1 0 0
                      1 1 1 0 0 0 1 0
                      1 1 1 0 0 0 1 0
                      1 1 1 0 0 0 1 0
                      1 1 1 0 0 1 1 0
                      1 1 1 0 0 0 0 0]);
        L = bwlabel(BW,4);
        [r,c] = find(L == 2);
 
    See also <a href="matlab:help bweuler">bweuler</a>, <a href="matlab:help bwlabeln">bwlabeln</a>, <a href="matlab:help bwselect">bwselect</a>.

    Reference page in Help browser
       <a href="matlab:doc bwlabel">doc bwlabel</a>




L = bwlabel(BWr,4);

whos 
  Name      Size                           Bytes  Class

  BWr     240x320                          76800  logical array
  I       240x320x3                       230400  uint8 array
  Ir      240x320                          76800  uint8 array
  L       240x320                         614400  double array
  l        28x1                             8852  struct array
  s         1x13                              26  char array

Grand total is 461793 elements using 1007278 bytes

max(max(L))

ans =

     8

for i=1:8,
ind = find(L == i);
length(i)
end

ans =

     1


ans =

     1


ans =

     1


ans =

     1


ans =

     1


ans =

     1


ans =

     1


ans =

     1

for i=1:8,
ind = find(L == i);
length(ind)
end

ans =

     1


ans =

       17136


ans =

     1


ans =

     1


ans =

     1


ans =

     1


ans =

     1


ans =

     1

for i=1:8,
ind = find(L == i);
if (length(ind) < 10)
   L(ind) = 0;
emd
end
end
??? Undefined function or variable 'emd'.

for i=1:8,
ind = find(L == i);
if (length(ind) < 10)
L(ind) = 0;
end
end


max(max(L))

ans =

     2

ind0 = find(L == 0);
whos
  Name       Size                           Bytes  Class

  BWr      240x320                          76800  logical array
  I        240x320x3                       230400  uint8 array
  Ir       240x320                          76800  uint8 array
  L        240x320                         614400  double array
  ans        1x1                                8  double array
  i          1x1                                8  double array
  ind        1x1                                8  double array
  ind0   59664x1                           477312  double array
  l         28x1                             8852  struct array
  s          1x13                              26  char array

Grand total is 521460 elements using 1484614 bytes

ind = find(L == 2);
whos
  Name       Size                           Bytes  Class

  BWr      240x320                          76800  logical array
  I        240x320x3                       230400  uint8 array
  Ir       240x320                          76800  uint8 array
  L        240x320                         614400  double array
  ans        1x1                                8  double array
  i          1x1                                8  double array
  ind    17136x1                           137088  double array
  ind0   59664x1                           477312  double array
  l         28x1                             8852  struct array
  s          1x13                              26  char array

Grand total is 538595 elements using 1621694 bytes

BWr(ind0) = 0;
figure, imshow(BWr);
help bwlabel
 BWLABEL Label connected components in binary image.
    L = BWLABEL(BW,N) returns a matrix L, of the same size as BW, containing
    labels for the connected components in BW. N can have a value of either
    4 or 8, where 4 specifies 4-connected objects and 8 specifies
    8-connected objects; if the argument is omitted, it defaults to 8.
 
    The elements of L are integer values greater than or equal to 0.  The
    pixels labeled 0 are the background.  The pixels labeled 1 make up one
    object, the pixels labeled 2 make up a second object, and so on.
 
    [L,NUM] = BWLABEL(BW,N) returns in NUM the number of connected objects
    found in BW.
 
    Note: Comparing BWLABEL and BWLABELN
    ------------------------------------
    BWLABEL supports 2-D inputs only, whereas BWLABELN support any 
    input dimension.  In some cases you might prefer to use BWLABELN even
    for 2-D problems because it can be faster.  If you have a 2-D input
    whose objects are relatively "thick" in the vertical direction,
    BWLABEL will probably be faster; otherwise BWLABELN will probably be
    faster.
 
    Class Support
    -------------
    BW can be logical or numeric, and it must be real, 2-D, and
    nonsparse.  L is double.
 
    Example
    -------
        BW = logical([1 1 1 0 0 0 0 0
                      1 1 1 0 1 1 0 0
                      1 1 1 0 1 1 0 0
                      1 1 1 0 0 0 1 0
                      1 1 1 0 0 0 1 0
                      1 1 1 0 0 0 1 0
                      1 1 1 0 0 1 1 0
                      1 1 1 0 0 0 0 0]);
        L = bwlabel(BW,4);
        [r,c] = find(L == 2);
 
    See also <a href="matlab:help bweuler">bweuler</a>, <a href="matlab:help bwlabeln">bwlabeln</a>, <a href="matlab:help bwselect">bwselect</a>.

    Reference page in Help browser
       <a href="matlab:doc bwlabel">doc bwlabel</a>



help bwlabeln
 BWLABELN Label connected components in N-D binary image.
    L = BWLABELN(BW) returns a label matrix, L, containing labels for the
    connected components in BW.  BW can have any dimension; L is the same
    size as BW.  The elements of L are integer values greater than or equal
    to 0.  The pixels labeled 0 are the background.  The pixels labeled 1
    make up one object, the pixels labeled 2 make up a second object, and so
    on.  The default connectivity is 8 for two dimensions, 26 for three
    dimensions, and CONNDEF(NDIMS(BW),'maximal') for higher dimensions.
 
    [L,NUM] = BWLABELN(BW) returns the number of connected objects found
    in BW.
 
    [L,NUM] = BWLABELN(BW,CONN) specifies the desired connectivity.  CONN
    may have the following scalar values:
 
        4     two-dimensional four-connected neighborhood
        8     two-dimensional eight-connected neighborhood
        6     three-dimensional six-connected neighborhood
        18    three-dimensional 18-connected neighborhood
        26    three-dimensional 26-connected neighborhood
 
    Connectivity may be defined in a more general way for any dimension by
    using for CONN a 3-by-3-by- ... -by-3 matrix of 0s and 1s.  The 1-valued
    elements define neighborhood locations relative to the center element of
    CONN.  CONN must be symmetric about its center element.
 
    Note: Comparing BWLABEL and BWLABELN
    ------------------------------------
    BWLABEL supports 2-D inputs only, whereas BWLABELN support any 
    input dimension.  In some cases you might prefer to use BWLABELN even
    for 2-D problems because it can be faster.  If you have a 2-D input
    whose objects are relatively "thick" in the vertical direction,
    BWLABEL will probably be faster; otherwise BWLABELN will probably be
    faster.
 
    Class Support
    -------------
    BW can be numeric or logical, and it must be real and nonsparse.  L
    is double.
 
    Example
    -------
        BW = cat(3,[1 1 0; 0 0 0; 1 0 0],...
                   [0 1 0; 0 0 0; 0 1 0],...
                   [0 1 1; 0 0 0; 0 0 1])
        bwlabeln(BW)
 
    See also <a href="matlab:help bwlabel">bwlabel</a>.

    Reference page in Help browser
       <a href="matlab:doc bwlabeln">doc bwlabeln</a>



help bwlabel
 BWLABEL Label connected components in binary image.
    L = BWLABEL(BW,N) returns a matrix L, of the same size as BW, containing
    labels for the connected components in BW. N can have a value of either
    4 or 8, where 4 specifies 4-connected objects and 8 specifies
    8-connected objects; if the argument is omitted, it defaults to 8.
 
    The elements of L are integer values greater than or equal to 0.  The
    pixels labeled 0 are the background.  The pixels labeled 1 make up one
    object, the pixels labeled 2 make up a second object, and so on.
 
    [L,NUM] = BWLABEL(BW,N) returns in NUM the number of connected objects
    found in BW.
 
    Note: Comparing BWLABEL and BWLABELN
    ------------------------------------
    BWLABEL supports 2-D inputs only, whereas BWLABELN support any 
    input dimension.  In some cases you might prefer to use BWLABELN even
    for 2-D problems because it can be faster.  If you have a 2-D input
    whose objects are relatively "thick" in the vertical direction,
    BWLABEL will probably be faster; otherwise BWLABELN will probably be
    faster.
 
    Class Support
    -------------
    BW can be logical or numeric, and it must be real, 2-D, and
    nonsparse.  L is double.
 
    Example
    -------
        BW = logical([1 1 1 0 0 0 0 0
                      1 1 1 0 1 1 0 0
                      1 1 1 0 1 1 0 0
                      1 1 1 0 0 0 1 0
                      1 1 1 0 0 0 1 0
                      1 1 1 0 0 0 1 0
                      1 1 1 0 0 1 1 0
                      1 1 1 0 0 0 0 0]);
        L = bwlabel(BW,4);
        [r,c] = find(L == 2);
 
    See also <a href="matlab:help bweuler">bweuler</a>, <a href="matlab:help bwlabeln">bwlabeln</a>, <a href="matlab:help bwselect">bwselect</a>.

    Reference page in Help browser
       <a href="matlab:doc bwlabel">doc bwlabel</a>



help bwselect
 BWSELECT Select objects in binary image.
    BW2 = BWSELECT(BW1,C,R,N) returns a binary image containing
    the objects that overlap the pixel (R,C). R and C can be
    scalars or equal-length vectors. If R and C are vectors, BW2
    contains the set of objects overlapping with any of the
    pixels (R(k),C(k)). N can have a value of either 4 or 8 (the
    default), where 4 specifies 4-connected objects and 8
    specifies 8-connected objects. Objects are connected sets of
    "on" pixels (i.e., having value of 1).
 
    BW2 = BWSELECT(BW1,N) displays the image BW1 on the screen
    and lets you select the (R,C) coordinates using the mouse. If
    you omit BW1, BWSELECT operates on the image in the current
    axes. Use normal button clicks to add points. Pressing
    <BACKSPACE> or <DELETE> removes the previously selected
    point. A shift-click, right-click, or double-click selects
    the final point; pressing <RETURN> finishes the selection
    without adding a point.
 
    [BW2,IDX] = BWSELECT(...) returns the linear indices of the 
    pixels belonging to the selected objects.
 
    BW2 = BWSELECT(X,Y,BW1,Xi,Yi,N) uses the vectors X and Y to
    establish a nondefault spatial coordinate system for BW1. Xi
    and Yi are scalars or equal-length vectors that specify
    locations in this coordinate system.
 
    [X,Y,BW2,IDX,Xi,Yi] = BWSELECT(...) returns the XData and
    YData in X and Y; the output image in BW2; linear indices of
    the pixels belonging to the selected objects in IDX; and the 
    specified spatial coordinates Xi and Yi.
 
    If bwselect is called with no output arguments, the resulting
    image is displayed in a new figure.
 
    Class Support
    ------------- 
    The input image BW1 can be logical or any numeric type and 
    must be 2-D and nonsparse.  The output image BW2 is logical.
 
    Example
    -------
        BW1 = imread('text.png');
        c = [126 187 11];
        r = [34 172 20];
        BW2 = bwselect(BW1,c,r,4);
        imview(BW1)
        imview(BW2)
 
    See also <a href="matlab:help imfill">imfill</a>, <a href="matlab:help bwlabel">bwlabel</a>, <a href="matlab:help roipoly">roipoly</a>, <a href="matlab:help roifill">roifill</a>.

    Reference page in Help browser
       <a href="matlab:doc bwselect">doc bwselect</a>




help images
  Image Processing Toolbox
  Version 4.2 (R14) 05-May-2004
 
  Image display.
    colorbar       - Display colorbar (MATLAB Toolbox).
    <a href="matlab:help getimage">getimage</a>       - Get image data from axes.
    image          - Create and display image object (MATLAB Toolbox).
    imagesc        - Scale data and display as image (MATLAB Toolbox).
    <a href="matlab:help immovie">immovie</a>        - Make movie from multiframe image.
    <a href="matlab:help imshow">imshow</a>         - Display image in Handle Graphics figure.
    <a href="matlab:help imview">imview</a>         - Display image in the image viewer.
    <a href="matlab:help montage">montage</a>        - Display multiple image frames as rectangular montage.
    movie          - Play recorded movie frames (MATLAB Toolbox).
    <a href="matlab:help subimage">subimage</a>       - Display multiple images in single figure.
    <a href="matlab:help truesize">truesize</a>       - Adjust display size of image.
    <a href="matlab:help warp">warp</a>           - Display image as texture-mapped surface.
 
  Image file I/O.
    <a href="matlab:help dicominfo">dicominfo</a>      - Read metadata from a DICOM message.
    <a href="matlab:help dicomread">dicomread</a>      - Read a DICOM image.
    <a href="matlab:help dicomuid">dicomuid</a>       - Generate DICOM unique identifier.
    <a href="matlab:help dicomwrite">dicomwrite</a>     - Write a DICOM image.
    <a href="matlab:help dicom-dict.txt">dicom-dict.txt</a> - Text file containing DICOM data dictionary.
    imfinfo        - Return information about image file (MATLAB Toolbox).
    imread         - Read image file (MATLAB Toolbox).
    imwrite        - Write image file (MATLAB Toolbox).
 
  Image arithmetic.
    <a href="matlab:help imabsdiff">imabsdiff</a>      - Compute absolute difference of two images.
    <a href="matlab:help imadd">imadd</a>          - Add two images, or add constant to image.
    <a href="matlab:help imcomplement">imcomplement</a>   - Complement image.
    <a href="matlab:help imdivide">imdivide</a>       - Divide two images, or divide image by constant.
    <a href="matlab:help imlincomb">imlincomb</a>      - Compute linear combination of images.
    <a href="matlab:help immultiply">immultiply</a>     - Multiply two images, or multiply image by constant.
    <a href="matlab:help imsubtract">imsubtract</a>     - Subtract two images, or subtract constant from image.
    <a href="matlab:help ippl">ippl</a>           - Check for presence of IPPL.
 
  Geometric transformations.
    <a href="matlab:help checkerboard">checkerboard</a>   - Create checkerboard image.
    <a href="matlab:help findbounds">findbounds</a>     - Find output bounds for geometric transformation.
    <a href="matlab:help fliptform">fliptform</a>      - Flip the input and output roles of a TFORM struct.
    <a href="matlab:help imcrop">imcrop</a>         - Crop image.
    <a href="matlab:help imresize">imresize</a>       - Resize image.
    <a href="matlab:help imrotate">imrotate</a>       - Rotate image.
    <a href="matlab:help imtransform">imtransform</a>    - Apply geometric transformation to image.
    <a href="matlab:help makeresampler">makeresampler</a>  - Create resampler structure.
    <a href="matlab:help maketform">maketform</a>      - Create geometric transformation structure (TFORM).
    <a href="matlab:help tformarray">tformarray</a>     - Apply geometric transformation to N-D array.
    <a href="matlab:help tformfwd">tformfwd</a>       - Apply forward geometric transformation.
    <a href="matlab:help tforminv">tforminv</a>       - Apply inverse geometric transformation.
 
  Image registration.
    <a href="matlab:help cpstruct2pairs">cpstruct2pairs</a> - Convert CPSTRUCT to valid pairs of control points.
    <a href="matlab:help cp2tform">cp2tform</a>       - Infer geometric transformation from control point pairs.
    <a href="matlab:help cpcorr">cpcorr</a>         - Tune control point locations using cross-correlation. 
    <a href="matlab:help cpselect">cpselect</a>       - Control point selection tool.
    <a href="matlab:help normxcorr2">normxcorr2</a>     - Normalized two-dimensional cross-correlation.
 
  Pixel values and statistics.
    <a href="matlab:help corr2">corr2</a>          - Compute 2-D correlation coefficient.
    <a href="matlab:help imcontour">imcontour</a>      - Create contour plot of image data.
    <a href="matlab:help imhist">imhist</a>         - Display histogram of image data.
    <a href="matlab:help impixel">impixel</a>        - Determine pixel color values.
    <a href="matlab:help improfile">improfile</a>      - Compute pixel-value cross-sections along line segments.
    <a href="matlab:help mean2">mean2</a>          - Compute mean of matrix elements.
    <a href="matlab:help pixval">pixval</a>         - Display information about image pixels.
    <a href="matlab:help regionprops">regionprops</a>    - Measure properties of image regions.
    <a href="matlab:help std2">std2</a>           - Compute standard deviation of matrix elements.
 
  Image analysis.
    <a href="matlab:help bwboundaries">bwboundaries</a>    - Trace region boundaries in binary image.
    <a href="matlab:help bwtraceboundary">bwtraceboundary</a> - Trace object in binary image.
    <a href="matlab:help edge">edge</a>            - Find edges in intensity image.
    <a href="matlab:help qtdecomp">qtdecomp</a>        - Perform quadtree decomposition.
    <a href="matlab:help qtgetblk">qtgetblk</a>        - Get block values in quadtree decomposition.
    <a href="matlab:help qtsetblk">qtsetblk</a>        - Set block values in quadtree decomposition.
 
  Image enhancement.
    <a href="matlab:help adapthisteq">adapthisteq</a>    - Perform adaptive histogram equalization using CLAHE.
    <a href="matlab:help decorrstretch">decorrstretch</a>  - Apply a decorrelation stretch to a multichannel image.
    <a href="matlab:help histeq">histeq</a>         - Enhance contrast using histogram equalization.
    <a href="matlab:help imadjust">imadjust</a>       - Adjust image intensity values or colormap.
    <a href="matlab:help imnoise">imnoise</a>        - Add noise to an image.
    <a href="matlab:help medfilt2">medfilt2</a>       - Perform 2-D median filtering.
    <a href="matlab:help ordfilt2">ordfilt2</a>       - Perform 2-D order-statistic filtering.
    <a href="matlab:help stretchlim">stretchlim</a>     - Find limits to contrast stretch an image.
    <a href="matlab:help uintlut">uintlut</a>        - Computes new array values based on lookup table.
    <a href="matlab:help wiener2">wiener2</a>        - Perform 2-D adaptive noise-removal filtering.
 
  Linear filtering.
    <a href="matlab:help convmtx2">convmtx2</a>       - Compute 2-D convolution matrix.
    <a href="matlab:help fspecial">fspecial</a>       - Create predefined filters.
    <a href="matlab:help imfilter">imfilter</a>       - Filter 2-D and N-D images.
 
  Linear 2-D filter design.
    freqspace      - Determine 2-D frequency response spacing (MATLAB Toolbox).
    <a href="matlab:help freqz2">freqz2</a>         - Compute 2-D frequency response.
    <a href="matlab:help fsamp2">fsamp2</a>         - Design 2-D FIR filter using frequency sampling.
    <a href="matlab:help ftrans2">ftrans2</a>        - Design 2-D FIR filter using frequency transformation.
    <a href="matlab:help fwind1">fwind1</a>         - Design 2-D FIR filter using 1-D window method.
    <a href="matlab:help fwind2">fwind2</a>         - Design 2-D FIR filter using 2-D window method.
 
  Image deblurring.
    <a href="matlab:help deconvblind">deconvblind</a>    - Deblur image using blind deconvolution.
    <a href="matlab:help deconvlucy">deconvlucy</a>     - Deblur image using Lucy-Richardson method.
    <a href="matlab:help deconvreg">deconvreg</a>      - Deblur image using regularized filter.
    <a href="matlab:help deconvwnr">deconvwnr</a>      - Deblur image using Wiener filter.
    <a href="matlab:help edgetaper">edgetaper</a>      - Taper edges using point-spread function.
    <a href="matlab:help otf2psf">otf2psf</a>        - Optical transfer function to point-spread function.
    <a href="matlab:help psf2otf">psf2otf</a>        - Point-spread function to optical transfer function.
 
  Image transforms.
    <a href="matlab:help dct2">dct2</a>           - 2-D discrete cosine transform.
    <a href="matlab:help dctmtx">dctmtx</a>         - Discrete cosine transform matrix.
    <a href="matlab:help fan2para">fan2para</a>       - Convert fan-beam projections to parallel-beam.
    <a href="matlab:help fanbeam">fanbeam</a>        - Compute fan-beam transform.
    fft2           - 2-D fast Fourier transform (MATLAB Toolbox).
    fftn           - N-D fast Fourier transform (MATLAB Toolbox).
    fftshift       - Reverse quadrants of output of FFT (MATLAB Toolbox).
    <a href="matlab:help idct2">idct2</a>          - 2-D inverse discrete cosine transform.
    ifft2          - 2-D inverse fast Fourier transform (MATLAB Toolbox).
    ifftn          - N-D inverse fast Fourier transform (MATLAB Toolbox).
    <a href="matlab:help ifanbeam">ifanbeam</a>       - Compute inverse fan-beam transform.
    <a href="matlab:help iradon">iradon</a>         - Compute inverse Radon transform.
    <a href="matlab:help para2fan">para2fan</a>       - Convert parallel-beam projections to fan-beam.
    <a href="matlab:help phantom">phantom</a>        - Generate a head phantom image.
    <a href="matlab:help radon">radon</a>          - Compute Radon transform.
 
  Neighborhood and block processing.
    <a href="matlab:help bestblk">bestblk</a>        - Choose block size for block processing.
    <a href="matlab:help blkproc">blkproc</a>        - Implement distinct block processing for image.
    <a href="matlab:help col2im">col2im</a>         - Rearrange matrix columns into blocks.
    <a href="matlab:help colfilt">colfilt</a>        - Columnwise neighborhood operations.
    <a href="matlab:help im2col">im2col</a>         - Rearrange image blocks into columns.
    <a href="matlab:help nlfilter">nlfilter</a>       - Perform general sliding-neighborhood operations.
 
  Morphological operations (intensity and binary images).
    <a href="matlab:help conndef">conndef</a>        - Default connectivity.
    <a href="matlab:help imbothat">imbothat</a>       - Perform bottom-hat filtering.
    <a href="matlab:help imclearborder">imclearborder</a>  - Suppress light structures connected to image border.
    <a href="matlab:help imclose">imclose</a>        - Close image.
    <a href="matlab:help imdilate">imdilate</a>       - Dilate image.
    <a href="matlab:help imerode">imerode</a>        - Erode image.
    <a href="matlab:help imextendedmax">imextendedmax</a>  - Extended-maxima transform.
    <a href="matlab:help imextendedmin">imextendedmin</a>  - Extended-minima transform.
    <a href="matlab:help imfill">imfill</a>         - Fill image regions and holes.
    <a href="matlab:help imhmax">imhmax</a>         - H-maxima transform.
    <a href="matlab:help imhmin">imhmin</a>         - H-minima transform.
    <a href="matlab:help imimposemin">imimposemin</a>    - Impose minima.
    <a href="matlab:help imopen">imopen</a>         - Open image.
    <a href="matlab:help imreconstruct">imreconstruct</a>  - Morphological reconstruction.
    <a href="matlab:help imregionalmax">imregionalmax</a>  - Regional maxima.
    <a href="matlab:help imregionalmin">imregionalmin</a>  - Regional minima.
    <a href="matlab:help imtophat">imtophat</a>       - Perform tophat filtering.
    <a href="matlab:help watershed">watershed</a>      - Watershed transform.
 
  Morphological operations (binary images).
    <a href="matlab:help applylut">applylut</a>       - Perform neighborhood operations using lookup tables.
    <a href="matlab:help bwarea">bwarea</a>         - Compute area of objects in binary image.
    <a href="matlab:help bwareaopen">bwareaopen</a>     - Binary area open (remove small objects).
    <a href="matlab:help bwdist">bwdist</a>         - Compute distance transform of binary image.
    <a href="matlab:help bweuler">bweuler</a>        - Compute Euler number of binary image.
    <a href="matlab:help bwhitmiss">bwhitmiss</a>      - Binary hit-miss operation.
    <a href="matlab:help bwlabel">bwlabel</a>        - Label connected components in 2-D binary image.
    <a href="matlab:help bwlabeln">bwlabeln</a>       - Label connected components in N-D binary image.
    <a href="matlab:help bwmorph">bwmorph</a>        - Perform morphological operations on binary image.
    <a href="matlab:help bwpack">bwpack</a>         - Pack binary image.
    <a href="matlab:help bwperim">bwperim</a>        - Determine perimeter of objects in binary image.
    <a href="matlab:help bwselect">bwselect</a>       - Select objects in binary image.
    <a href="matlab:help bwulterode">bwulterode</a>     - Ultimate erosion.
    <a href="matlab:help bwunpack">bwunpack</a>       - Unpack binary image.
    <a href="matlab:help makelut">makelut</a>        - Construct lookup table for use with applylut.
 
  Structuring element (STREL) creation and manipulation.
    <a href="matlab:help getheight">getheight</a>      - Get strel height.
    <a href="matlab:help getneighbors">getneighbors</a>   - Get offset location and height of strel neighbors
    <a href="matlab:help getnhood">getnhood</a>       - Get strel neighborhood.
    <a href="matlab:help getsequence">getsequence</a>    - Get sequence of decomposed strels.
    <a href="matlab:help isflat">isflat</a>         - Return true for flat strels.
    <a href="matlab:help reflect">reflect</a>        - Reflect strel about its center.
    <a href="matlab:help strel">strel</a>          - Create morphological structuring element.
    <a href="matlab:help translate">translate</a>      - Translate strel.
 
  Region-based processing.
    <a href="matlab:help poly2mask">poly2mask</a>      - Convert region-of-interest polygon to mask.
    <a href="matlab:help roicolor">roicolor</a>       - Select region of interest, based on color.
    <a href="matlab:help roifill">roifill</a>        - Smoothly interpolate within arbitrary region.
    <a href="matlab:help roifilt2">roifilt2</a>       - Filter a region of interest.
    <a href="matlab:help roipoly">roipoly</a>        - Select polygonal region of interest.
 
  Colormap manipulation.
    brighten       - Brighten or darken colormap (MATLAB Toolbox).
    <a href="matlab:help cmpermute">cmpermute</a>      - Rearrange colors in colormap.
    <a href="matlab:help cmunique">cmunique</a>       - Find unique colormap colors and corresponding image.
    colormap       - Set or get color lookup table (MATLAB Toolbox).
    <a href="matlab:help imapprox">imapprox</a>       - Approximate indexed image by one with fewer colors.
    rgbplot        - Plot RGB colormap components (MATLAB Toolbox).
 
  Color space conversions.
    <a href="matlab:help applycform">applycform</a>     - Apply device-independent color space transformation.
    hsv2rgb        - Convert HSV values to RGB color space (MATLAB Toolbox).
    <a href="matlab:help iccread">iccread</a>        - Read ICC color profile.
    <a href="matlab:help lab2double">lab2double</a>     - Convert Lab color values to double.
    <a href="matlab:help lab2uint16">lab2uint16</a>     - Convert Lab color values to uint16.
    <a href="matlab:help lab2uint8">lab2uint8</a>      - Convert Lab color values to uint8.
    <a href="matlab:help makecform">makecform</a>      - Create device-independent color space transform structure.
    <a href="matlab:help ntsc2rgb">ntsc2rgb</a>       - Convert NTSC values to RGB color space.
    rgb2hsv        - Convert RGB values to HSV color space (MATLAB Toolbox).
    <a href="matlab:help rgb2ntsc">rgb2ntsc</a>       - Convert RGB values to NTSC color space.
    <a href="matlab:help rgb2ycbcr">rgb2ycbcr</a>      - Convert RGB values to YCBCR color space.
    <a href="matlab:help whitepoint">whitepoint</a>     - Returns XYZ values of standard illuminants.
    <a href="matlab:help xyz2double">xyz2double</a>     - Convert XYZ color values to double.
    <a href="matlab:help xyz2uint16">xyz2uint16</a>     - Convert XYZ color values to uint16.
    <a href="matlab:help ycbcr2rgb">ycbcr2rgb</a>      - Convert YCBCR values to RGB color space.
 
  Interactive mouse utility functions.
    <a href="matlab:help getline">getline</a>        - Select polyline with mouse.
    <a href="matlab:help getpts">getpts</a>         - Select points with mouse.
    <a href="matlab:help getrect">getrect</a>        - Select rectangle with mouse.
 
  ICC color profiles.
    <a href="matlab:help lab8.icm">lab8.icm</a>       - 8-bit Lab profile.
    <a href="matlab:help monitor.icm">monitor.icm</a>    - Typical monitor profile.
                     Sequel Imaging, Inc., used with permission.
    <a href="matlab:help sRGB.icm">sRGB.icm</a>       - sRGB profile.
                     Hewlett-Packard, used with permission.
    <a href="matlab:help swopcmyk.icm">swopcmyk.icm</a>   - CMYK input profile.
                     Eastman Kodak, used with permission.
 
  Array operations.
    circshift      - Shift array circularly (MATLAB Toolbox).
    <a href="matlab:help padarray">padarray</a>       - Pad array.
 
  Image types and type conversions.
    <a href="matlab:help dither">dither</a>         - Convert image using dithering.
    <a href="matlab:help gray2ind">gray2ind</a>       - Convert intensity image to indexed image.
    <a href="matlab:help grayslice">grayslice</a>      - Create indexed image from intensity image by thresholding.
    <a href="matlab:help graythresh">graythresh</a>     - Compute global image threshold using Otsu's method.
    <a href="matlab:help im2bw">im2bw</a>          - Convert image to binary image by thresholding.
    <a href="matlab:help im2double">im2double</a>      - Convert image array to double precision.
    im2java        - Convert image to Java image (MATLAB Toolbox).
    <a href="matlab:help im2java2d">im2java2d</a>      - Convert image to Java BufferedImage.
    <a href="matlab:help im2uint8">im2uint8</a>       - Convert image array to 8-bit unsigned integers.
    <a href="matlab:help im2uint16">im2uint16</a>      - Convert image array to 16-bit unsigned integers.
    <a href="matlab:help ind2gray">ind2gray</a>       - Convert indexed image to intensity image.
    ind2rgb        - Convert indexed image to RGB image (MATLAB Toolbox).
    <a href="matlab:help label2rgb">label2rgb</a>      - Convert label matrix to RGB image.
    <a href="matlab:help mat2gray">mat2gray</a>       - Convert matrix to intensity image.
    <a href="matlab:help rgb2gray">rgb2gray</a>       - Convert RGB image or colormap to grayscale.
    <a href="matlab:help rgb2ind">rgb2ind</a>        - Convert RGB image to indexed image.
 
  Toolbox preferences.
    <a href="matlab:help iptgetpref">iptgetpref</a>     - Get value of Image Processing Toolbox preference.
    <a href="matlab:help iptsetpref">iptsetpref</a>     - Set value of Image Processing Toolbox preference.
 
  Demos.
    <a href="matlab:help iptdemos">iptdemos</a>       - Index of Image Processing Toolbox demos.




whos
  Name       Size                           Bytes  Class

  BWr      240x320                          76800  logical array
  I        240x320x3                       230400  uint8 array
  Ir       240x320                          76800  uint8 array
  L        240x320                         614400  double array
  ans        1x1                                8  double array
  i          1x1                                8  double array
  ind    17136x1                           137088  double array
  ind0   59664x1                           477312  double array
  l         28x1                             8852  struct array
  s          1x13                              26  char array

Grand total is 538595 elements using 1621694 bytes

Ir1 = uint8(BWr);
max(max(Ir1))

ans =

    1

ind = find(Ir1 == 1);
Ir1(ind) = 255;
max(mmax(Ir1))
??? Undefined command/function 'mmax'.

max(max(Ir1))

ans =

  255

Ig1 = uint8(zeros(size(Ir1)));


whos
  Name       Size                           Bytes  Class

  BWr      240x320                          76800  logical array
  I        240x320x3                       230400  uint8 array
  Ig1      240x320                          76800  uint8 array
  Ir       240x320                          76800  uint8 array
  Ir1      240x320                          76800  uint8 array
  L        240x320                         614400  double array
  ans        1x1                                1  uint8 array
  i          1x1                                8  double array
  ind    17136x1                           137088  double array
  ind0   59664x1                           477312  double array
  l         28x1                             8852  struct array
  s          1x13                              26  char array

Grand total is 692195 elements using 1775287 bytes

Ib1 = uint8(zeros(size(Ir1)));
Ic = [Ir1; Ig1; Ib1];
whos
  Name       Size                           Bytes  Class

  BWr      240x320                          76800  logical array
  I        240x320x3                       230400  uint8 array
  Ib1      240x320                          76800  uint8 array
  Ic       720x320                         230400  uint8 array
  Ig1      240x320                          76800  uint8 array
  Ir       240x320                          76800  uint8 array
  Ir1      240x320                          76800  uint8 array
  L        240x320                         614400  double array
  ans        1x1                                1  uint8 array
  i          1x1                                8  double array
  ind    17136x1                           137088  double array
  ind0   59664x1                           477312  double array
  l         28x1                             8852  struct array
  s          1x13                              26  char array

Grand total is 999395 elements using 2082487 bytes

Ic = [Ir1 Ig1 Ib1];
whos
  Name       Size                           Bytes  Class

  BWr      240x320                          76800  logical array
  I        240x320x3                       230400  uint8 array
  Ib1      240x320                          76800  uint8 array
  Ic       240x960                         230400  uint8 array
  Ig1      240x320                          76800  uint8 array
  Ir       240x320                          76800  uint8 array
  Ir1      240x320                          76800  uint8 array
  L        240x320                         614400  double array
  ans        1x1                                1  uint8 array
  i          1x1                                8  double array
  ind    17136x1                           137088  double array
  ind0   59664x1                           477312  double array
  l         28x1                             8852  struct array
  s          1x13                              26  char array

Grand total is 999395 elements using 2082487 bytes

Ic(:,:,1) = Ir1;
??? Subscripted assignment dimension mismatch.

clear Ic
Ic(:,:,1) = Ir1;
Ic(:,:,2) = Ig1;
Ic(:,:,3) = Ib1;
imshow(Ic); 
Ic(:,:,2) = Ir1;
imshow(Ic);
Ic(:,:,2) = Ib1;
Ic(:,:,3) = Ir1;
imshow(Ic);
Ic(:,:,1) = Ib1;
Ic(:,:,2) = Ir1;
imshow(Ic);
diary;