matlab手掌静脉处理

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