function?R3=veinp（I）;
%veinp?is?used?to?obtain?vein?texture
%input:is?the?ROI?imagemust?be?rgb?format.
sigma=3;
X=rgb2gray（I）;
gausFilter?=?fspecial（‘gaussian‘[3?3]sigma）;
Y=imfilter（IgausFilter‘replicate‘）;
[MNnDims]=size（Y）;
Y=im2double（Y）;%?獲取圖像的尺寸和波段數
Z=Y;
for?i=1:nDims??%?對每個波段依次進行灰度拉伸
????Sp=Y（::i）;
????MaxDN=max（max（Sp））;
????MinDN=min（min（Sp））;
????Sp=（Sp-MinDN）/（MaxDN-MinDN）;??%?灰度拉伸公式
????Z（::i）=Sp;
end
%?將灰度拉伸結果保存在Z中
%?Z=histeq（Z）;
Z=rgb2gray（Z）;
tf=medfilt2（Z[99]）;
O=segNiBlack（tf150.08）;
O=~O;
O1=medfilt2（O[77]）;
O1=medfilt2（O1[55]）;
B=[0?1?01?1?10?1?0];
P1=imdilat