MATLAB血管骨架提取方法详解，包括血管分割和背景去除技巧

function untitled3（imgpathoutpath） 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=[50ones（1100） 20ones（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=m