一种光弹试验力链信息自动提取方法

1.本发明涉及数字图像处理技术领域，具体涉及一种光弹试验力链信息自动提取方法。


背景技术：

2.颗粒物质是由大量离散固体颗粒相互作用而形成的复杂体系。岩体工程中，因岩石颗粒介质流动而造成的山体滑坡、尾矿坝溃坝、泥石流、隧道塌方、溜井堵塞等工程地质灾害，严重威胁施工人员与设备安全。揭示岩石颗粒介质体系间的宏观力学行为和微观力学机理是开展上述灾害预测与防护的依据。颗粒介质的宏、微观力学性质及接触力传递路径常通过光弹试验的力链参量进行表征，光弹力链信息的量化与提取是其研究的关键环节，具有重要的科学与工程意义。
3.光弹试验是一种基于光学的应力测量方法。早在1957年，研究人员第一次利用光弹试验技术，清晰的观测到了力在颗粒介质中的传递光带。该力传递光带于1995年被定义为力链。此后20多年来，众多学者利用光弹试验对力链展开了大量的研究。如有文献表明，有学者借助光弹试验量化了顶煤开采过程中颗粒间接触力的大小，阐述了工作面矿压的形成机理；还有学者研究了颗粒物质的土拱效应，探讨了力链在内拱和外拱中的传递；此外，有人分析了颗粒物质在封堵地层裂缝中的作用，讨论了强弱力链的演化规律。而申请号为cn202210609208.8，名称为“基于力链特性量化的沥青混合料传递荷载行为评价方法”的中国发明专利申请，通过数值模拟的数据，得到了沥青混合料间接触信息，完成了从力链角度对沥青混合料中力传递行为的评价，但该发明方法却不适用于实际物理光弹试验中的力链信息提取。
4.上述研究成果对于促进颗粒物质力学的发展有着重要意义。但通过分析发现，关于光弹试验中力链信息的提取与分析，现有研究多针对力链宏观分布特征、接触力的量化以及强弱力链的区分开展研究。而且对于每条力链分支点的精确识别、单条力链信息的定量提取方法及力链动态的演化规律研究仍较少。


技术实现要素：

5.本发明的目的是提供一种光弹试验力链信息自动提取方法，以克服光弹试验中力链分支点的精确判别问题，实现力链信息的逐条量化，促进颗粒物质力学的发展。
6.为解决上述问题，本发明采用的技术方案是：
7.一种光弹试验力链信息自动提取方法，包括以下步骤：
8.(1)通过有偏振镜及无偏振镜分别获取同一次试验状态下的颗粒图像，并对图像进行预处理；
9.(2)提取颗粒几何信息：对预处理后的无偏振镜的光弹图像，进行霍夫变换与欧式距离计算，提取出颗粒圆心坐标、颗粒半径、颗粒间接触点位置和颗粒间接触角度；
10.(3)提取颗粒接触力信息：基于识别出的颗粒几何信息，通过在有偏振镜的图像上
进行掩膜运算，得到颗粒的彩色梯度均方值；再结合接触力与彩色梯度均方值的标定曲线，提取出颗粒的接触力信息；
11.(4)提取力链网络结构信息：结合颗粒的几何信息、接触力信息及力链的构成判据，将角度阈值转换为距离阈值，进而判别出强、弱力链类型，提取出力链宏观分布状态，识别出力链长度、力链数量和力链方位信息。
12.所述无偏振镜的颗粒图像预处理按以下方法进行：将无偏振镜的颗粒图像进行二值化、腐蚀、膨胀、边缘检测并进行孔洞填充。
13.所述有偏振镜的颗粒图像预处理是将有偏振镜的颗粒图像进行滤波。
14.所述力链的构成判据包括强弱力链判据、力链长度判据以及距离阈值判据，所述强弱力链判据为：当颗粒接触力f大于或等于颗粒体系内的平均接触力时，判定该颗粒位于强力链上；反之，判定该颗粒位于弱力链上或处于不受外力作用的状态。
15.所述力链长度判据为：一条力链上的颗粒个数即为力链长度，其值须大于等于3，即接触点个数须大于等于2。
16.所述距离阈值判据为：两相邻接触点之间的距离d，需大于等于距离阈值dc，距离阈值的计算公式如下：
17.dc＝(2-0.0002358*θ
c2.04
)*r
18.其中，dc为距离阈值；θc为角度阈值，取值越小，力链越直；r为颗粒半径。
19.所述步骤(4)的具体过程如下：
20.4-1：输入已识别出的颗粒几何信息和颗粒接触力信息，并计算出颗粒体系的平均接触力
[0021][0022]
4-2：生成空矩阵c
all
,c1,c2,c3,c4,c
fc
，用于后续数据的存储；
[0023]
4-3：选取出的强接触颗粒信息存储于c
all
中，并按接触点坐标y值由大到小排序；
[0024]
4-4：判别c
all
中的接触点数量是否大于等于2：如果是，则选取一个起始接触点a，并对a周边的接触点ai进行搜索，将接触点a和搜索结果ai存储于c1中；如果否，则结束步骤(4)；
[0025]
4-5：判别c1中的接触点数量是否大于等于2：如果是，则将c1拷贝一份为c2；如果否，则从c
all
中删除与c1的交集，清空c1，跳转至步骤4-4；
[0026]
4-6：从c1中删除与起始接触点a构成4个颗粒的接触点ai；
[0027]
4-7：从c1中删除与c4的交集；
[0028]
4-8：从c1中选取距起始接触点a最远的接触点ai，并计算两接触点间的距离，赋值为d；
[0029]
4-9：判别距离d是否大于等于距离阈值dc：如果是，则将接触点a和ai存储于c3，另新的起始接触点为ai，c2拷贝为c4，清空c1和c2，并跳转至步骤4-4；如果否，则执行如下步骤：
[0030]
4-9-1：判别c3中接触点数量是否大于等于2：如果是，则将c3拷贝为c
fc
，并转至步骤4-9-2；如果否，则直接转至步骤4-9-2；
[0031]
4-9-2：从c
all
中删除接触点a，从c
all
中删除与c3的交集；
[0032]
4-9-3：清空c1,c2,c3,c4，并跳转至步骤4-4，开始新一条力链信息的检索。
[0033]
采用上述技术方案所产生的有益效果是：本发明提供的一种光弹试验力链信息自动提取方法，能够准确地从物理光弹试验图像中提取出力链网络，动态刻画出颗粒物质的力链特征，在颗粒几何信息提取、力链分支点判别、力链检索、力链分布及其方位角表征等方面存在显著优势。研究成果为颗粒物质宏观力学行为和微观力学机理的研究提供了科学依据。
附图说明
[0034]
图1为本发明的整体流程图。
[0035]
图2为本发明的基于成链判据的力链网络结构信息自动提取流程图。
[0036]
图3为本发明实施例提供的光弹试验物理模型图。
[0037]
图4为本发明实施例提供的光弹试验物理图像局部放大图，其中，a为有偏振镜下的光弹试验图像，b为无偏振镜下的光弹试验图像。
[0038]
图5为本发明实施例提供的光弹试验物理图像预处理结果图。
[0039]
图6为本发明实施例提供的光弹试验物理图像几何信息提取结果图，其中，a为颗粒识别结果，b为接触网络识别结果。
[0040]
图7为本发明实施例提供的接触力f与彩色梯度均方值的标定曲线图。
[0041]
图8为本发明实施例提供的接触点筛选与力链生成示意图。
[0042]
图9为本发明实施例提供的强力链信息提取结果局部图。
[0043]
图10为本发明实施例提供的接触力信息演化规律图。
[0044]
图11为本发明实施例提供的第10次放矿下光弹试验物理图像与力链信息提取结果图，其中，a为光弹试验物理图像整体图，b为力链网络信息提取结果图，c为强力链方位角提取结果，d为弱力链方位角提取结果。
具体实施方式
[0045]
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
[0046]
本实施例以大量放矿同步充填无顶柱留矿采矿方法下的物理光弹试验为例，采用本发明方法进行了光弹试验的力链信息自动提取，整体流程如图1所示，具体步骤如下：
[0047]
步骤1：通过有偏振镜及无偏振镜分别获取同一次试验状态下的颗粒图像，并对图像进行预处理。输入的同一次试验状态下有偏振镜和无偏振镜的两张光弹试验颗粒图像如图4所示，其中图4a为有偏振镜下的光弹试验图像，图4b为无偏振镜下的光弹试验图像。在对图4b的形态学预处理过程中，二值化阈值为0.125，删除小面积斑点面积为1300像素，边缘检测所用算子为canny算子，预处理结果如图5所示。
[0048]
步骤2：提取颗粒几何信息：对预处理后的无偏振镜的光弹图像，进行霍夫变换与欧式距离计算，提取出颗粒圆心坐标、颗粒半径、颗粒间接触点位置和颗粒间接触角度，其中距离判别时的修正系数为1.05，颗粒间接触角度范围为-90
°
～90
°
，提取结果如图6所示。
[0049]
步骤3：提取颗粒接触力信息：基于识别出的颗粒几何信息，通过在有偏振镜的图像上进行掩膜运算，得到颗粒的彩色梯度均方值；再结合接触力与彩色梯度均方值的标定
曲线，提取出颗粒的接触力信息。接触力与彩色梯度均方值的标定曲线如图7所示。
[0050]
步骤4：提取力链网络结构信息：结合颗粒的几何信息、接触力信息及力链的构成判据，将角度阈值转换为距离阈值，进而判别出强、弱力链类型，提取出力链宏观分布状态，识别出力链长度、力链数量和力链方位信息。
[0051]
力链的构成判据包括：
[0052]
(1)强弱力链判据：当颗粒接触力f大于或等于颗粒体系内的平均接触力f时，判定该颗粒位于强力链上；反之，判定该颗粒位于弱力链上或处于不受外力作用的状态。
[0053]
(2)力链长度判据：一条力链上的颗粒个数即为力链长度，其值须大于等于3，即接触点个数须大于等于2。
[0054]
(3)距离阈值判据：两相邻接触点之间的距离d，需大于等于距离阈值dc。距离阈值的计算公式如下：
[0055]
dc＝(2-0.0002358*θ
c2.04
)*r
[0056]
其中，dc为距离阈值；θc为角度阈值，取值越小，力链越直；r为颗粒半径。此实施例中，距离阈值为1.66r。
[0057]
步骤4中，基于成链判据的力链网络结构信息自动提取流程，如图2所示，具体步骤如下(以图8为例)：
[0058]
步骤4-1：输入已识别出的颗粒几何信息和颗粒接触力信息，并计算出颗粒体系的平均接触力f。
[0059]
步骤4-2：生成空矩阵c
all
,c1,c2,c3,c4,c
fc
，用于后续数据的存储。
[0060]
步骤4-3：选取出的强接触颗粒信息存储于c
all
中，并按接触点坐标y值由大到小排序。如需选取出弱接触颗粒信息，只需将判别条件改为即可。此时c
all
中的接触点序号有1号、2号、3号、4号、5号、6号、7号、8号、9号、10号、11号、12号、13号。
[0061]
步骤4-4：判别c
all
中的接触点数量是否大于等于2。如果是，则选取一个起始接触点a，并对a周边的接触点ai进行搜索，将接触点a和搜索结果ai存储于c1中。如果否，则结束步骤4。此实施例中，接触点搜索距离为1.1r。
[0062]
因图8a的计算结果较为简单，现默认已完成一次检索，现以起始接触点a为3号进行说明，相关结果如图8b所示。此时起始检索点a为3号，搜索到的周边接触点ai为1号、2号、4号、5号、6号、7号、8号、9号、10号。c1中的接触点序号有1号、2号、3号、4号、5号、6号、7号、8号、9号、10号。c4中的接触点序号有1号、2号、3号、4号。
[0063]
步骤4-5：判别c1中的接触点数量是否大于等于2。如果是，则将c1拷贝一份为c2。如果否，则从c
all
中删除与c1的交集，清空c1，跳转至步骤4-4。此时c2中的接触点序号有1号、2号、3号、4号、5号、6号、7号、8号、9号、10号。
[0064]
步骤4-6：从c1中删除与起始接触点a构成4个颗粒的接触点ai，即删除5号和10号。此时c1中的接触点序号有1号、2号、3号、4号、6号、7号、8号、9号。
[0065]
步骤4-7：从c1中删除与c4的交集，即删除1号、2号、3号、4号。此时c1中的接触点序号有6号、7号、8号、9号。
[0066]
步骤4-8：从c1中选取距起始接触点a最远的接触点ai，并计算两接触点间的距离，赋值为d。此时起始接触点a为3号，最远的接触点ai为7号，距离d为2r。
[0067]
步骤4-9：判别距离d是否大于等于距离阈值dc。如果是，则将接触点a和ai存储于
c3，另新的起始接触点为ai，c2拷贝为c4，清空c1和c2，并跳转至步骤4-4。此实施例中d＝2r＞dc＝1.66r。此时c3中接触点序号有1号、3号、7号。新的起始接触点为7号。c4中的接触点序号有1号、2号、3号、4号、5号、6号、7号、8号、9号、10号。
[0068]
因每次循环过程判断一致，在此不过多赘述，相关循环检索结果如图8c和图8d所示。最终识别结果如图9至图11所示。
[0069]
通过本实施例说明，本发明提供的一种光弹试验力链信息自动提取方法，能够准确地从物理光弹试验图像中提取出力链网络，动态刻画出颗粒物质的力链特征，在颗粒几何信息提取、力链分支点判别、力链检索、力链分布及其方位角表征等方面存在显著优势。研究成果为颗粒物质宏观力学行为和微观力学机理的研究提供了科学依据。
[0070]
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。