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