hello everyone
I really need help as soon as possible
I want a solution to calculate the diameter of each white line in the image(Vessel Diameter Measurement)
MATLAB: How to calculate the diameter of a vessel
calculatediameterimage processing
Related Solutions
Michael, try this:
clc; % Clear the command window.
close all; % Close all figures (except those of imtool.)
clear; % Erase all existing variables. Or clearvars if you want.
workspace; % Make sure the workspace panel is showing.
format long g;format compact;fontSize = 20;% Read in mat file.
storedStructure = load('binaryimage.mat');binaryImage = storedStructure.BWfinal;% Display the image.
subplot(2, 2, 1);imshow(binaryImage, []);title('Original Binary Image', 'FontSize', fontSize, 'Interpreter', 'None');% Set up figure properties:
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);% Get rid of tool bar and pulldown menus that are along top of figure.
set(gcf, 'Toolbar', 'none', 'Menu', 'none');% Give a name to the title bar.
set(gcf, 'Name', 'Demo by ImageAnalyst', 'NumberTitle', 'Off') drawnow;% Fill the image.
binaryImage(1,:) = true;binaryImage(end,:) = true;binaryImage = imfill(binaryImage, 'holes');% Display the image.subplot(2, 2, 2);imshow(binaryImage, []);title('Original Binary Image', 'FontSize', fontSize, 'Interpreter', 'None');drawnow;% Get the Euclidean Distance Transform.
binaryImage(1,:) = false;binaryImage(end,:) = false;edtImage = bwdist(~binaryImage);% Display the image.subplot(2, 2, 3);imshow(edtImage, []);title('Distance Transform Image', 'FontSize', fontSize, 'Interpreter', 'None');drawnow;skelImage = bwmorph(binaryImage, 'skel', inf);skelImage = bwmorph(skelImage, 'spur', 40);% There should be just one now. Let's check
[labeledImage, numLines] = bwlabel(skelImage);fprintf('Found %d lines\n', numLines);% Display the image.subplot(2, 2, 4);imshow(skelImage, []);title('Skeleton Image', 'FontSize', fontSize, 'Interpreter', 'None');% Measure the radius be looking along the skeleton of the distance transform.
meanRadius = mean(edtImage(skelImage))meanDiameter = 2 * meanRadiusmessage = sprintf('Mean Radius = %.1f pixels.\nMean Diameter = %.1f pixels',... meanRadius, meanDiameter);uiwait(helpdlg(message));
You forgot to attach your image. Usually when people ask for image processing advice, they attach an image.
I would compute the Euclidean Distance Transform for the binary image with bwdist().
edtImage = bwdist(binaryImage);Then ask regionprops for the pixel values of the EDT image.
labeledImage = bwlabel(binaryImage);props = regionprops(labeledImage, edtImage, 'PixelValues');If the max pixel value for that ring is acceptable, store that index.
numRings = length(props)keepers = false(1, numRings)for k = 1 : numRings if max(props(k).PixelValues) > someValue % This ring is thick enough. So keep it.
keepers(k) = true; endendThen use ismember(mask, indexes) to extract just the rings that meet your criteria.
filteredBinaryImage = ismember(labeledImage, find(keepers));The code is just off the top of my head - untested, so you may have to adapt it some. If you want the average thickness of the ring, then you'll have to call bwmorph() to skeletonize it and then multiply the skeleton image by the EDT image to get the thickness just along the centerlines. Then get the mean of all those pixels.
If you need more help, post your image and ring thickness criteria for keeping and excluding rings.
Best Answer