基于MATLAB的数字图像处理（指纹增强）

川长思鸟来 2024-04-06 09:22 77阅读 0赞

## 数字图像处理 ##

### 需求分析 ###

指纹在生活中的各个领域都有丰富的应用，如在刑侦时利用指纹追踪嫌犯，但很多情况下获取的指纹不够清晰，或纹线之间有断裂、粘连，因此指纹增强是很有必要的。通过滤波、形态学变换、时域频域转换等操作，使指纹图像在图片背景中变得清晰可见，是指纹增强所追求的。

在本次作业中，需要对如图 1 所示三张指纹图像进行增强，为表述简便，下文分别称为图片 1、图片 2 和图片 3。

![在这里插入图片描述][9586b389dce14a588cc2aaaafc652812.png]  
![在这里插入图片描述][1bb0310d8c99491a8e152500c475b86c.png]  
![在这里插入图片描述][1dfddfca186e4cff8eaf0bce0be4a16a.png]

图 1 需增强指纹

### 算法设计 ###

指纹增强主要有以下几个步骤——前景分割、局部脊线方向估计、局部脊线频率估计、脊线增强，其中前三个步骤之间并无先后关系，在实际操作中相辅相成，现对各步骤分别进行介绍。

### 前景分割 ###

前景分割即从原图像中找出指纹区域，主要通过检测指纹区域的频域特征和均、方差特征来实现。首先对原图进行切分，将原图分成多个相邻的 8*8 小块，为保证切割完整性，将图像边界向外扩充，将长和宽变为 8 的倍数，为了傅里叶变换的效果更好，对每个 8*8 小块向外扩充至 32\*32。

通常指纹纹线可以近似看作二位正弦函数，因此可以通过二位正弦函数的特征来判断指纹区域。使用2(2())函数进行 2D DFT，取绝对值。首先使用2()找到\*32 区域中灰度值最大的像素点，一般应是中心直流点，去掉最亮点之后，再次使用 2()找到次亮点，若次两点的个数为 2 且关于中心点对称，则认为可能是指纹区域。计算两个两点到中心点的平均距离，作为局部脊线频率的测度，同时可以通过距离筛选掉一些不属于指纹区域的背景部分。

上述步骤完成后，对于图片 1 和图片 2 已经基本可以分出指纹区域，而对于背景较为杂乱的图片 3，上述操作只能找到一些离散的块，对这些块做扩展，认定其周围 48\*48 的区域均为指纹区域，此时得到离散块中，有一些白色部分较多的或方差较大的，可以通过均值筛掉，这样可以获得一个如图 2 所示的 mask，mask 为 1 的部分是认定的指纹区域。之后再对图片 3 的原图进行处理，进行一系列腐蚀、膨胀的形态学运算，通过灰度值、方差等条件大致去除背景中其他部分，与 mask 结合后效果如图 3.

![在这里插入图片描述][3f7c60d692a14426b3dd92120aefd22c.png]

图 2 图片 3 指纹区域 mask

![在这里插入图片描述][51cd4ac773874ec193a0f5f4376ac816.png]  
![在这里插入图片描述][7ff48298ef3f4a139274439f6ca704d8.png]

图 3 图片 3 前景分割效果

局部脊线方向估计

已知对于二维正弦函数，其 2D DFT 幅度谱两个分量的连线与正弦波方向垂直。在 2.1 的操作中，已获得指纹区域 32\*32 小块的 2D DFT 幅度谱，找到两个次亮点后，通过两个亮点的连线与竖直方向的夹角计算指纹脊线方向，此过程中对于不同方向的部分要分类讨论，将最终获得的脊线方向角存入矩阵。以图片 1 为例，使用函数展示方向图如图 4

，可以看到该方向图能反映出指纹的大致走势，但仍有些较为杂乱的部分。

接下来对方向图进行平滑，将方向矩阵每个元素乘以 2，计算其正弦和余弦，分别对正弦图和余弦图进行平滑，再对正余弦矩阵调用 atan2 函数，除以 2，得到平滑后的方向矩阵，再使用函数展示方向图如图 4（b），可以看到此时方向图已较为平滑。

对于图片 3，由于在 2.1 的操作中，去除了穿指纹图像而过的两条黑线，生成的方向图也有明显的断裂，对于此现象做方向的手动调整和两次平滑，最终处理效果如图 5(b)，可以看到中间断裂处已被基本连接。

![在这里插入图片描述][9807d499619f4d60b36221247cee55c9.png]  
![在这里插入图片描述][72e496b89a414d8caa1f552e72bd8a2f.png]

初始 (b)平滑后图 4 图片 1 方向图

![在这里插入图片描述][780990ab3e794a92b18584de7075b824.png]  
![在这里插入图片描述][d91a9bcf46924042899fc08b5f642fb2.png]

初始 (b)处理后图 5 图片 3 方向图

### 局部脊线频率估计 ###

已知对于二维正弦函数，其 2D DFT 幅度谱两个分量连线距离反映了该点频率大小。在 2.1 的操作中，用两个次亮点到中心点的平均距离作为指纹区域的判断条件，判断为是指纹区域后，将存入频率矩阵，获得频率矩阵后，对其进行空域平滑滤波，以图片 1 为例，将距离大小转换为灰度值大小显示如图 6.

![在这里插入图片描述][af18e36e6119410ca740218f25f7d9ac.png]  
![在这里插入图片描述][2a036832ed8941dcbe328834a02e54e0.png]

图 6 图片 1 频率图

### 脊线增强 ###

最后一步即脊线增强，本次作业中分别使用了notch滤波和Gabor滤波。其中，notch滤波是在频域上进行的，通过之前获得的指纹区域方向图和频率图，计算滤波器两个峰的位置，即仅允许这两个位置处的频率分量通过，滤除其他频率的分量，从而获得增强的指纹图像；而Gabor滤波是在空域上实现的，使用函数(, , )，其中是滤波器两个峰之间的距离，是方向角，其转移函数如图7所示。  
![在这里插入图片描述][2f9b796a5ccc476fac71eceadbdb1888.png]

图7 Gabor滤波器转移函数

对三幅图片分别使用Notch和Gabor滤波，再进行二值化处理，对比如图8所示，可以看到Notch滤波对指纹信息如指纹方向的保留较好，但指纹纹线之间会有断裂；而Gabor滤波增强后的指纹纹线较为完整、平滑，但纹线方向有一定扭曲和毛刺，另外，在运行中也发现

### Gabor滤波的用时长于Notch滤波。 ###

在指纹增强时，不妨将两种滤波结果结合起来看，能较好的得到指纹信息。

![在这里插入图片描述][7a1574f6cf0d4901a4faa4daf0555225.png]

![在这里插入图片描述][6ce2a5a23d644cfc81ddc541dc23fecb.png]

Notch (b) Gabor 图 8 图片 1 增强结果

实验结果每个步骤的实验结果以图片形式展示，顺序为图片 1、图片 2、图片 3.

前景分割

![在这里插入图片描述][8ceb84f84b1d4a05a26f082ee113e63c.png]  
![在这里插入图片描述][e6b1b402a84440dcb87154ef64345ab2.png]  
![在这里插入图片描述][36b7d9dae9f34f21bd6d513a1d334d2d.png]

方向图

初始

![在这里插入图片描述][453b8b7aa8244800ab67eb754d11ceda.png]  
![在这里插入图片描述][4b16a8d331474273b98372308771242c.png]  
![在这里插入图片描述][c2dcbb42be984756b6f05ac78def0077.png]

平滑后

![在这里插入图片描述][649215ff11fa437c8aca7108ebc30dba.png]  
![在这里插入图片描述][f1d4002648f94e058b39540d69b22b33.png]  
![在这里插入图片描述][955371af7bdd4ecea3f3e83409e960b9.png]

频率图

![在这里插入图片描述][d4389d61cf0643909e76ec53b7b6446b.png]

![在这里插入图片描述][1b7b0cf4b1394b1c96ea1cf48b6e9a8d.png]  
![在这里插入图片描述][08faa366572d4220bb3ff7e785d91868.png]

增强图

Notch

![在这里插入图片描述][44127fead5c241cabda4cdf6b97d9ed8.png]  
![在这里插入图片描述][e2cd7e4acd544288a23728d4ff2548c7.png]  
![在这里插入图片描述][bf758ba3f8cd4f3ab330a5aefb9f5f55.png]

Gabor

![在这里插入图片描述][32934c619276444e8527e7c4f4a881a2.png]  
![在这里插入图片描述][e47e494321894cccab990113bea6ae2c.png]  
![在这里插入图片描述][7ce196f8fb8c4b7ca1c9cd24779c97c9.png]

[9586b389dce14a588cc2aaaafc652812.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/de2f06f277f64a9483b95e37fefc406c.png
[1bb0310d8c99491a8e152500c475b86c.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/85022dc59c454e2fb116de8debe491ca.png
[1dfddfca186e4cff8eaf0bce0be4a16a.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/2677d575767742eda184abc385533f80.png
[3f7c60d692a14426b3dd92120aefd22c.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/8929e07684544ca9a227635f58cfbbc0.png
[51cd4ac773874ec193a0f5f4376ac816.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/0a67aa5d2038449783039ee663b0d4d6.png
[7ff48298ef3f4a139274439f6ca704d8.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/a39f48794acf4d8284307e38eeb06ef1.png
[9807d499619f4d60b36221247cee55c9.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/b4180feaea564b56aea66af77b946f33.png
[72e496b89a414d8caa1f552e72bd8a2f.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/610815a87d214cca8cc739252915c619.png
[780990ab3e794a92b18584de7075b824.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/7b7c98b2ac4241339ab0ae2e259f027d.png
[d91a9bcf46924042899fc08b5f642fb2.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/85fa8fb165dd4901a91334c605b2ae5e.png
[af18e36e6119410ca740218f25f7d9ac.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/973193f5fa4a46399abc0c990539ad25.png
[2a036832ed8941dcbe328834a02e54e0.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/e8d945938854445e9dab08fa55298868.png
[2f9b796a5ccc476fac71eceadbdb1888.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/a8bac3ebeb264e278a15f7c79778c1fa.png
[7a1574f6cf0d4901a4faa4daf0555225.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/51d62c2aada74ce09bf52b6ba3933960.png
[6ce2a5a23d644cfc81ddc541dc23fecb.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/b140ada33d454d0a81ee5c4b0050036d.png
[8ceb84f84b1d4a05a26f082ee113e63c.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/8b37e2deb5424d2b8a76504626f7145f.png
[e6b1b402a84440dcb87154ef64345ab2.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/177f071dc67a422a94e27a1e61daba89.png
[36b7d9dae9f34f21bd6d513a1d334d2d.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/aee927c2046a4c87801c7f988cb8ca40.png
[453b8b7aa8244800ab67eb754d11ceda.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/f7c62c76b35741aab3428f0758565353.png
[4b16a8d331474273b98372308771242c.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/58cafd9f40b0482c9b0312fe43dbdaa2.png
[c2dcbb42be984756b6f05ac78def0077.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/dbf6b1c7e2224d43b586e3ee31ad30fb.png
[649215ff11fa437c8aca7108ebc30dba.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/243852a0e24647adb70de060cd46b1a5.png
[f1d4002648f94e058b39540d69b22b33.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/df91609942944aefa7e9a361f7823767.png
[955371af7bdd4ecea3f3e83409e960b9.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/d8cc99341bda41938912a44156812757.png
[d4389d61cf0643909e76ec53b7b6446b.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/3d3ee8bb873a4dfead02ad05994c8546.png
[1b7b0cf4b1394b1c96ea1cf48b6e9a8d.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/d6736f51426e4beda5bb60180072339b.png
[08faa366572d4220bb3ff7e785d91868.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/ea52e5ad80cd450d8b857de228c31216.png
[44127fead5c241cabda4cdf6b97d9ed8.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/d8e120fa0dd7487bad786739e42b77ba.png
[e2cd7e4acd544288a23728d4ff2548c7.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/702f055313f54e6abff760bbf0a30fc6.png
[bf758ba3f8cd4f3ab330a5aefb9f5f55.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/7e9da488cabf4969a5b20fafcf6ea559.png
[32934c619276444e8527e7c4f4a881a2.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/ca3d1247c926452987069bdd7a0233f3.png
[e47e494321894cccab990113bea6ae2c.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/0903f5c8356d41199aaa6600f2d4b86f.png
[7ce196f8fb8c4b7ca1c9cd24779c97c9.png]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/06/1aba3dc73d2343a3a106d779ae8a4e58.png