%%将图1中的虹膜区域提取出来放在矩形final14中
final14=zeros(64,512);%%定义一个矩形
for i=1:64
for j=1:512
rad=j*2*pi/512;
xpos3=round((i+r11)*cos(rad)+row11);%r11是内圆半径，（row11，col11)是原的半径（内外圆同心）
ypos3=round((i+r11)*sin(rad)+col11);
final14(i,j)=final13(xpos3,ypos3);%为final13,图2为final14,从图2中看出，fianl14并没有得到final13中的虹膜部分
end
end
图1
图2