【血管检测】 mom方法结合Hessian和曲线拟合血管直径测量【含Matlab源码 1970期】.zip

% 正常血管，检测单根比较准确 clc;clear all;close all; %% 图像预处理 Im = imread('..\masks\IMG_0062_mask.png'); % Im = imread('1.png'); % IM = Im; im = imread('..\seg_result\IMG_0062_out.png'); % im = imread('otsu_1.png'); IM1(:,:,1) = im2double(Im); IM1(:,:,2) = im2double(Im); IM1(:,:,3) = im2double(Im); IM2(:,:,1) = im2double(im); IM2(:,:,2) = im2double(im); IM2(:,:,3) = im2double(im); % se = strel('disk',1); % im = imerode(im,se); % 分割的结果太宽了，腐蚀一下，去留随意 S = IM1; BW = bwmorph(im,'thin',inf); % 细化得到骨架 [m,n] = find(BW==1); for i = 1:length(m) S(m(i),n(i),1) = 255; S(m(i),n(i),2) = 0; S(m(i),n(i),3) = 0; end hf = figure('NumberTitle', 'off', 'Name', '骨架'); imshow(S) %% % 手动输入 % BW(:,1:48) = 0; % 截取 % BW(1:187,:) = 0; % BW(:,83:end) = 0; % BW(219:end,:) = 0; % 交互式 rect = getrect(hf); BW(:,1:round(rect(1))) = 0; % 截取 BW(1:round(rect(2)),:) = 0; BW(:,round(rect(1)+rect(3)):end) = 0; BW(round(rect(2)+rect(4)):end,:) = 0; %% 中心线坐标 [y_l,x_l] = find(BW==1); sp=spap2(1,4,x_l,y_l); % 拟合曲线 [x,idx] = unique(x_l,'stable'); % 中点纵坐标 k_s1 = fnval(fnder(sp), x); % 血管流向斜率 k_s2 = -1./k_s1; % 法一：拟合曲线求法线斜率 y = fnval(sp, x); % 中点横坐标 %% 画截面 && 算真实半径 real_R = zeros(1,length(x)); % 多个中心点的半径 r = 10; for j = 1:length(x) % ang = HessianAng(m, x0, y0)+ pi/2; % 法二：Hessian求特征向量(效果差) ang = atan(k_s2(j)); x0 = round(y(j)); y0 = x(j); m = 2; % rr = MoM(x0, y0, k_s1(j), k_s2(j), r, m, Im2); % mom评价方法 rr = MoMforSeg2(x0, y0, r ,Im); for i=-rr:rr xx = x0 + i * sin(ang); yy = y0 + i * cos(ang); IM1(round(xx), round(yy),3) = 255; IM1(round(xx), round(yy),1) = 0; IM1(round(xx), round(yy),2) = 0; end real_R(j) = rr; end figure('NumberTitle', 'off', 'Name', '真实图像的直径'); imshow(IM1) %% 求分割后的直径 R = zeros(1,length(x)); % 多个中心点的半径 for j = 1:length(x) x0 = round(y(j)); y0 = x(j); ang = atan(k_s2(j)); r = 11; rr = MoMforSeg1(x0, y0, r ,im); % %画截面 for i=-rr:rr xx = x0 + i * sin(ang); yy = y0 + i * cos(ang); IM2(round(xx), round(yy),1) = 255; IM2(round(xx), round(yy),2) = 0; IM2(round(xx), round(yy),3) = 0; end R(j) = rr; end figure('NumberTitle', 'off', 'Name', '分割图像的直径'); imshow(IM2) figure('NumberTitle', 'off', 'Name', '直径对比'); plot(R,'r^'); hold on plot(real_R,'b-'); hold off xlabel('沿血管路径/像素','FontSize',11); ylabel('血管直径/像素','FontSize',11); legend({'测量值';'真实值'}); axis([0 length(real_R)+1 0 8]); %% 误差计算 e1 = abs(R - real_R); % 绝对误差 e2 = abs(R - real_R)./real_R; % 相对误差 disp('血管直径为：'); fprintf('%c', 8); % 去除行末回车 disp(real_R); disp('绝对误差为：'); fprintf('%c', 8); disp(e1); disp('相对误差为：'); fprintf('%c', 8); disp(e2); %% 新增部分 RMSE1 = sqrt(mean(e1.^2)); RMSE2 = sqrt(mean(e2.^2)); disp('RMSE1 = '); fprintf('%c', 8); disp(RMSE1); disp('RMSE2 = '); fprintf('%c', 8); disp(RMSE2);