{"id":337,"date":"2010-07-30T12:47:39","date_gmt":"2010-07-30T16:47:39","guid":{"rendered":"https:\/\/blogs.mathworks.com\/steve\/2010\/07\/30\/visualizing-regionprops-ellipse-measurements\/"},"modified":"2019-10-29T13:32:16","modified_gmt":"2019-10-29T17:32:16","slug":"visualizing-regionprops-ellipse-measurements","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/steve\/2010\/07\/30\/visualizing-regionprops-ellipse-measurements\/","title":{"rendered":"Visualizing regionprops ellipse measurements"},"content":{"rendered":"
\r\n

The Image Processing Toolbox function regionprops<\/tt> gives you measurements of shape-based measurements of image regions. It's pretty useful, and I have shown many examples of\r\n it on this blog. Today I want to show you how to visualize the ellipse-based measurements produced by regionprops<\/tt>.\r\n <\/p>\r\n

There are several supported measurements that are based on approximating regions by ellipses:<\/p>\r\n

\r\n
    \r\n
  • Orientation<\/li>\r\n
  • MajorAxisLength<\/li>\r\n
  • MinorAxisLength<\/li>\r\n
  • Eccentricity<\/li>\r\n <\/ul>\r\n <\/div>\r\n

    Here's an example.<\/p>

    url = 'https:\/\/blogs.mathworks.com\/images\/steve\/2010\/rice_binary.png'<\/span>;\r\nbw = imread(url);\r\nimshow(bw)<\/pre> 
    Call regionprops<\/tt> with the list of measurements you desire:\r\n   <\/p>
    s = regionprops(bw, 'Orientation'<\/span>, 'MajorAxisLength'<\/span>, ...<\/span>\r\n    'MinorAxisLength'<\/span>, 'Eccentricity'<\/span>)<\/pre>
    \r\ns = \r\n\r\n99x1 struct array with fields:\r\n    MajorAxisLength\r\n    MinorAxisLength\r\n    Eccentricity\r\n    Orientation\r\n\r\n<\/pre>
    From the output you can infer that regionprops<\/tt> found 99 objects in the binary image. It returned the desired measurements as a 99-element struct array. Here are the measurements\r\n      for the 10th object:\r\n   <\/p>
    s(10)<\/pre>
    \r\nans = \r\n\r\n    MajorAxisLength: 28.7693\r\n    MinorAxisLength: 9.4320\r\n       Eccentricity: 0.9447\r\n        Orientation: -27.0162\r\n\r\n<\/pre>
    How can you do a visual sanity check of the output? We should probably include some regionprops<\/tt> visualization utilities in the toolbox, but in the meantime I'll show you some code you can use to superimpose the fitted\r\n      ellipses over the original image.\r\n   <\/p>\r\n   
    First, we'll one additional measurement: the centroid.<\/p>
    s = regionprops(bw, 'Orientation'<\/span>, 'MajorAxisLength'<\/span>, ...<\/span>\r\n    'MinorAxisLength'<\/span>, 'Eccentricity'<\/span>, 'Centroid'<\/span>);<\/pre>
    For each ellipse, we'll use a parametric form of the ellipse equation to plot the outline of the ellipse over the image.<\/p>
    imshow(bw)\r\nhold on<\/span>\r\n\r\nphi = linspace(0,2*pi,50);\r\ncosphi = cos(phi);\r\nsinphi = sin(phi);\r\n\r\nfor<\/span> k = 1:length(s)\r\n    xbar = s(k).Centroid(1);\r\n    ybar = s(k).Centroid(2);\r\n\r\n    a = s(k).MajorAxisLength\/2;\r\n    b = s(k).MinorAxisLength\/2;\r\n\r\n    theta = pi*s(k).Orientation\/180;\r\n    R = [ cos(theta)   sin(theta)\r\n         -sin(theta)   cos(theta)];\r\n\r\n    xy = [a*cosphi; b*sinphi];\r\n    xy = R*xy;\r\n\r\n    x = xy(1,:) + xbar;\r\n    y = xy(2,:) + ybar;\r\n\r\n    plot(x,y,'r'<\/span>,'LineWidth'<\/span>,2);\r\nend<\/span>\r\nhold off<\/span><\/pre> 
    You can zoom in to see exactly why some of the ellipses don't seem to match objects in the image very well. Normally you'd\r\n      zoom interactively using the zoom controls in the figure toolbar, but of course I can't do that in this blog post, so I'll\r\n      use the axis<\/tt> function.\r\n   <\/p>
    axis([140 225 85 170])<\/pre> 
    Now you can see that the two objects in the center of the display are actually touching at one point. That caused regionprops<\/tt> to consider them to be a single object and to fit an ellipse to both of them together.\r\n   <\/p>