Matlab指纹细节提取

目前我对指纹验证和细节点提取很感兴趣。我在网上找到了以下代码，不知道是否有人可以解释一下？我查过质心、区域属性等，我对这些有一点了解，但下面的代码让我感到困惑！

fun=@minutie;
L = nlfilter(K,[3 3],fun);

%% Termination
LTerm=(L==1);
imshow(LTerm)
LTermLab=bwlabel(LTerm);
propTerm=regionprops(LTermLab,'Centroid');

CentroidTerm=round(cat(1,propTerm(:).Centroid));
imshow(~K)
set(gcf,'position',[1 1 600 600]);
hold on
plot(CentroidTerm(:,1),CentroidTerm(:,2),'ro')

%% Bifurcation
LBif=(L==3);
LBifLab=bwlabel(LBif);
propBif=regionprops(LBifLab,'Centroid','Image');
CentroidBif=round(cat(1,propBif(:).Centroid));
plot(CentroidBif(:,1),CentroidBif(:,2),'go')

该代码首先使用 3x3 像素的邻域过滤二值图像。nfilter是移动滤波器函数。它将遍历作为参数给出的图像中的所有像素，并根据相邻像素的值应用操作。

我不知道具体内容minutie过滤器，但从代码的其余部分来看，它可能会计算全 1 邻域中值为 1 的像素。换句话说，在线段末尾，它等于 1；当有 3 个分支（分叉）时，它等于 3。

Example:

让过滤器对邻近的过滤器进行求和，如下所示：

sum(block(1,1:3), block(3,1:3), block(2,1), block(2,3))*block(2, 2);

where block表示二值图像的每个像素周围的邻域。

在下面的左矩阵中（如果忽略边界例外），有一个位置的 1 在其 3x3 邻域中恰好有 1 个 1，在右矩阵中，有一个位置的 1 在其 3x3 邻域中恰好有 3 个 1邻里。

[0 0 0 0 0        [0 0 1 0 0
 0 0 0 0 0         0 0 1 0 0
 0 0 1 0 0         1 1 1 0 0
 0 0 1 0 0         0 0 1 0 0
 0 0 1 0 0]        0 0 1 0 0]

过滤后的输出将是：

[0 0 0 0 0        [0 0 0 0 0
 0 0 0 0 0         0 0 0 0 0
 0 0 1 0 0         0 0 3 0 0
 0 0 0 0 0         0 0 0 0 0
 0 0 0 0 0]        0 0 0 0 0]

它在左侧矩阵中发现了一个终止点，在右侧矩阵中发现了一个分叉点。

然后将滤波后的图像设置为值 1 和 3 的阈值，然后使用bwlabel and regionprops对我来说有些神秘†，因为分叉和终止点是单点，它们的位置就是它们的索引。我认为您可以使用以下方法简单地实现端点和分叉坐标的检测：

[It Jt]= find(L==1);
[Ib Jb]= find(L==3);

†我能想到的一个原因是，在matlab中图像和数组中的坐标是不同的，并且这两个函数以图像格式输出坐标，这更容易在原始图像上绘制。