一种用于从HRTEM图像定量评估硫化物催化剂活性相形态的晶格条纹几何聚类方法

《Applied Sciences》：A Geometric Clustering Method of Lattice Fringes for Quantitative Assessment of the Active Phase Morphology of Sulfide Catalysts from HRTEM Images

　　

本文提出了一种用于高分辨透射电子显微镜（HRTEM）图像中晶格条纹的几何聚类方法，以直接从显微照片实现硫化物催化剂活性相的定量形态学分析。研究人员开发了一个复合距离函数，用于晶格条纹的折线表示，该函数结合了催化剂的结构和形态特征，并具有直接的几何解释。所提出的函数满足半度量公理，确保其与不需要质心表示的基于密度和图的聚类算法（如DBSCAN）的兼容性。对于手动标注的晶格条纹分组，采用DBSCAN算法，最小簇大小为1。在合成数据和真实数据上的实验验证以及消融研究，证实了该方法的高精度和物理可解释性，使其成为硫化物催化剂高通量形态学分析和催化系统理性设计的有前景工具。

论文解读文章

研究背景与问题：硫化物催化剂（如MoS₂、WS₂及其促进形式NiMoS、CoMoS等）在生物质、木质素和热解油的催化处理中广泛用于脱除杂原子化合物。催化活性直接由活性相颗粒的形态决定，特别是平均片层长度和堆叠数。高分辨透射电子显微镜（HRTEM）是研究这些催化剂纳米颗粒形态的主要技术，能清晰显示对应于MoS₂/WS₂层状结构的晶格条纹。然而，手动分析HRTEM图像极其耗时且主观性强。现有自动或人工智能辅助方法多针对其他对象（如缺陷识别、相位重建、层数计数等），没有专门针对硫化物催化剂中将晶格条纹分组为堆叠并提取形态参数（片层长度和堆叠数）的物理驱动聚类框架。因此，研究人员开展了本研究，旨在提出一种基于折线间距离函数的几何聚类方法，用于自动化分组晶格条纹，实现高通量形态分析。

研究人员开展的研究及结论：研究人员提出了一种复合距离函数，用于衡量表示晶格条纹的折线（polyline）之间的相似性。该函数整合了硫化物纳米颗粒的物理化学约束：固定层间距（0.62–0.65 nm）、层间平行性、长度统一性以及对横向偏移的容忍性。函数由四个分量组成：基于归一化平均距离的指数分量、长度比、距离均匀性（控制平行性）、对角线比（控制偏移），并具有半度量性质，适用于基于密度和图的聚类算法。采用DBSCAN算法（最小簇大小设为1）进行聚类，将折线分组为对应纳米颗粒堆叠的簇，从而得到每个颗粒的堆叠数和平均片层长度。在合成数据上通过五个受控实验验证了距离函数对中心距、长度比、角度偏差、平行偏移和均匀缩放等几何参数的敏感性，证实了函数的单调性、可区分性和相对尺度不变性。在真实数据上，使用两个HRTEM图像片段（C0_0023和C0_0025），由专家手动标注了273条折线（194个颗粒）和196条折线（114个颗粒）。通过与专家标注对比，评估聚类质量，采用Kullback–Leibler散度（KL距离）和推土机距离（EMD）作为指标。优化后ε=4.33时，组合数据集的KL距离约为7×10^?4，EMD约为1×10^?2，且每个图像的聚类错误不超过3处。消融研究进一步证明，所有四个距离函数分量对实现高精度聚类都是必要的，完整变体在两个图像上的归一化互信息（NMI）错误率分别降低至0.17和0.18，优于仅使用平均欧氏距离的基线。论文发表在《Applied Sciences》。该研究建立了首个针对硫化物催化剂HRTEM图像中晶格条纹分组的物理驱动聚类框架，为高通量形态分析和理性设计催化系统提供了有前景的工具。

主要关键的技术方法（不超过250字）：研究人员开发了一个复合距离函数，包含四个几何解释的分量：（1）基于归一化平均距离的指数分量，用于控制层间距；（2）长度比分量，防止长度差异过大的层合并；（3）距离均匀性分量，用于评估平行性；（4）对角线比分量，用于检测不允许的横向偏移。该函数满足半度量公理。采用基于密度的DBSCAN聚类算法（最小簇大小设为1），输入为所有折线对的距离矩阵。在合成数据上通过五个受控实验系统验证了距离函数对中心距、长度比、角度偏差、平行偏移和均匀缩放的敏感性。在真实数据上，使用两个来自HRTEM图像的片段（C0_0023和C0_0025），由专家手动标注折线及对应颗粒分组作为金标准。聚类质量通过KL距离、EMD和NMI评估，并进行了消融研究以验证各分量的必要性。

研究结果：

6. 合成数据实验结果
6.1 数据模型：研究人员定义了五个关键几何参数（中心距、长度比、角度、平行偏移、均匀缩放），使用一对直线折线作为初始数据。
6.2 实验1：中心距依赖性。通过改变垂直偏移（1–32像素），验证了距离函数随中心距单调指数增长，主要由指数分量驱动，有效分离空间远离的折线。
6.3 实验2：长度比依赖性。固定第一折线长度77像素，变化第二折线长度（1–77像素），发现长度比分量呈双曲线型增长，当长度差异显著时，距离均匀性分量也增大，确保只合并长度相近的层。
6.4 实验3：角度依赖性。改变角度（?10°到10°），距离函数对称增长，在0°（平行）时最小，偏离时距离均匀性和对角线分量迅速增大，有效惩罚非平行性。
6.5 实验4：平行偏移依赖性。固定垂直偏移4像素，变化水平偏移（?40到40像素），距离函数对称增长，当偏移超过折线半长时加速，区分明显的横向错位。
6.6 实验5：均匀缩放依赖性。变化两折线长度（比例保持1.5），发现距离函数主要受相对长度差影响，绝对尺度变化影响小，确保尺度不变性。
6.7 合成数据实验讨论：五个实验证实距离函数对中心距最敏感，长度比、角度、平行偏移均有效响应，且对均匀缩放鲁棒，符合物理约束。

7. 真实数据实验结果
7.1 数据集：专家手动标注了两个HRTEM图像片段，C0_0023含273条折线（194个颗粒），C0_0025含196条折线（114个颗粒）。所有标注经第二位专家独立验证以保证一致性。距离函数中参考层间距设为16像素（对应0.63 nm）。
7.2 形态参数计算与质量评估：通过DBSCAN聚类得到每个颗粒的堆叠数和平均片层长度。与专家标注的堆叠数分布直方图对比，采用KL距离和EMD评估。优化DBSCAN邻域半径ε，对每个图像单独优化（分别为4.125和4.368），然后对组合数据集选择ε=4.33，得到KL≈7×10^?4，EMD≈1×10^?2。
7.3 聚类结果定性分析：每个图像仅出现3处聚类错误（C0_0023倾向于过度聚合，C0_0025倾向于将大颗粒拆分），表明算法与专家标注高度一致。

8. 消融研究
研究人员使用NMI（转换为1?NMI作为误差指标）评估不同距离函数变体的聚类质量。基线为仅用平均欧氏距离；三个部分变体分别移除长度比、距离均匀性、对角线分量；完整变体为公式(1)。结果显示完整变体在两个图像上均达到最低误差（C0_0023为0.17，C0_0025为0.18），优于所有部分变体，证明每个分量均对准确聚类必要。

总结讨论与结论：研究人员指出，自动分析HRTEM显微照片可降低劳动强度和解释主观性。所提出的距离函数首次形式化了硫化物催化剂纹理的模型表示，满足半度量性质，适用于密度和图聚类算法。在合成和真实数据上的实验均验证了其准确性和物理有效性。未来将结合自动线检测方法形成完整分析流程。
结论翻译：总体而言，所开发的方法展现出高精度和物理有效性，使其成为硫化物催化剂材料科学研究以及高活性催化系统理性设计的有前景工具。