bloodvessel 这个程序代码用于提取眼底图像的血管

这个程序代码用于提取眼底图像的血管，包括去噪、边缘检测

%function J = bv（I）
% Thitiporn Chanwimaluang
% tchanwim@gmail.com
% Department of Electrical and Computer Engineering
% Oklahoma State University

% This is the main function to perform blood vessel extraction

I=imread（‘im0162.ppm‘）; % Read an image
IG = rgb2gray（I）;
[MN] = size（IG）;
G（::2） = I（::2）;
G（::1） = 0;
G（::3） = 0;
Y = rgb2gray（G）;

sa = 2.0;
rt = mim（Ysa）;

[tt1e1cmtx] = myThreshold（rt）;

ms = 45;    
mk = msk（IGms）;

rt2 = 255*ones（MN）;
for i=1:M
    for j=1:N
        if rt（ij）>=tt1 & mk（ij）==255
            rt2（ij）=0;
        end
    end
end
J = im2bw（rt2）; 

J= ~J;
[LabelNum] = bwlabel（J）;
Lmtx = zeros（Num+11）;
for i=1:M
    for j=1:N
        Lmtx（double（Label（ij））+1） = Lmtx（double（Label（ij））+1） + 1;
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1164  2020-10-21 09:18  bloodvessel\main.m
     文件        1113  2020-10-21 09:18  bloodvessel\mim.m
     文件         643  2020-10-21 09:18  bloodvessel\msk.m
     文件         816  2020-10-21 09:18  bloodvessel\myThreshold.m
     文件        1304  2020-10-21 09:18  Retinal_Blood_Vessel_Extraction\license.txt
     文件      102054  2020-10-21 09:18  Retinal_Blood_Vessel_Extraction\Retinal_Blood_Vessel_Extraction\Input.bmp
     文件         869  2020-10-21 09:18  Retinal_Blood_Vessel_Extraction\Retinal_Blood_Vessel_Extraction\Main.m
     文件        1627  2020-10-21 09:18  Retinal_Blood_Vessel_Extraction\Retinal_Blood_Vessel_Extraction\VesselExtract.m