Matlab实现血管骨架提取

以冠状动脉血管为例,用MATLAB进行血管骨架提取,包括血管分割技术,去除背景

function untitled3（imgpathoutputpath）
I = imread （imgpath）;
% I=imread（‘F:\电工\dicom\IM000102_0040.jpg‘）;
I_gray=rgb2gray（I）;%转化为8位灰度BMP图像
sigma1 = 100;
gausFilter = fspecial（‘gaussian‘[1000 1000]sigma1）; 
blur=imfilter（I_graygausFilter‘replicate‘）;
imwrite（I_gray‘vessel.bmp‘）;
I=imread（‘vessel.bmp‘）;
template=[50*ones（1100） 20*ones（150） 40*ones（110）];
HI=histeq（Itemplate）;
FM=fspecial（‘motion‘30）;
FI=imfilter（HIFM‘replicate‘）;
FS=fspecial（‘unsharp‘）;
S=imfilter（FIFS‘replicate‘）;
% imshow（S）;
B=medfilt2（S）;%中值滤波
blur2=blur-B;
% imshow（blur2）
img=edge（ blur2‘sobel‘）;
imwrite（imgstrcat（outputpath）‘jpg‘）;
 end                
            

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3890  2018-11-17 19:12  test2.m

     文件      15673  2018-11-17 19:13  test3.jpg

     文件       3913  2018-11-17 19:13  test3.m

     文件      33251  2018-11-17 19:21  untitled.jpg

     文件        699  2018-11-17 19:20  Edge.m

     文件      28205  2018-11-17 19:07  test1.jpg

     文件       3960  2018-11-17 19:08  test1.m

     文件      27255  2018-11-17 19:09  test2.jpg

----------- ---------  ---------- -----  ----

               116846                    8