[发明专利]一种自动分析扫描式劳厄衍射图谱中衍射峰峰形的方法

权利要求书

1.一种自动分析扫描式劳厄衍射图谱中衍射峰峰形的方法，包括以下步骤：

S100：从来自同一实验的所有扫描式劳厄衍射图谱中任选一张衍射图谱I进行二维中值滤波，获得X射线荧光分布图谱I_f并计算所述X射线荧光分布图谱I_f上所有像素点的强度平均值I_fave；

S200：读取所述实验中的一张衍射图谱，计算所述衍射图谱I上各像素点I_ij的强度平均值I_ave和所述衍射图谱I去除荧光背景后的各像素点I_ij的强度值(I_re)_ij，获得去除荧光背景后的衍射图谱I_re；

S300：寻找所述去除荧光背景后的衍射图谱I_re上的衍射峰p，具体包括以下步骤：

S301：计算所述衍射图谱I的二值化图谱I_b，方法如下：

其中，I_reave表示去除荧光背景后的衍射图谱I_re上所有像素点的强度平均值；(I_b)_ij表示二值化图谱I_b中各点的强度值；α表示二值化阈值；

S302：使用二值图像连通域标记算法在二值图像上标记所述衍射峰p并通过提取目标正外接矩形的算法取得标记后的所述衍射峰p的正外接矩形以框取所述二值化图谱I_b中所有衍射峰p，对框取的衍射峰p使用二维高斯拟合法拟合出衍射峰p的中心位置；

S400：分析步骤S300中寻找的衍射峰p的峰形，具体包括以下步骤：

S401：通过坐标变换将被框取的衍射峰p中所有像素点转移至所述衍射峰p的倒易空间q；

S402：在所述倒易空间q内以衍射峰p的中心位置为坐标原点定义坐标系，计算衍射峰p在所述坐标系x轴方向上的积分拉长值S_p，方法如下：

其中，n表示转化到倒易空间中所有像素点的个数，I_i表示第i个像素点的强度值，r_i表示第i个像素点到x轴的距离；

S403：不断旋转步骤S402中所定义的坐标系，并计算S_p值，直到S_p值取到最小；

S404：通过高斯拟合求得衍射峰p在倒易空间q中x轴方向和y轴方向上的半高宽，并计算广义峰宽，方法如下：

其中，W_long表示衍射峰p在x轴上的半高宽，W_short表示衍射峰p在y轴上的半高宽，h_long表示高斯拟合求得的衍射峰p在x轴方向上的峰高，h_short表示高斯拟合求得的衍射峰p在y轴方向上的峰高，I_long表示x轴上各像素点的强度，I_short表示y轴上各像素点的强度；S500：判断所有需要分析的衍射图谱是否已经分析完成，若完成，则分析结束，否则重复执行步骤S200至步骤S400。

2.根据权利要求1所述的方法，其特征在于，步骤S100中，对所述衍射图谱I进行二维中值滤波时，确定取中值的邻域的大小。

3.根据权利要求1所述的方法，其特征在于，步骤S200中，所述衍射图谱I去除荧光背景后的各像素点I_ij的强度值(I_re)_ij通过下式获得：