[发明专利]一种确定岩石铸体薄片图片孔隙区域的方法、面孔率计算方法及其应用

权利要求书

1.一种确定岩石铸体薄片图片孔隙区域的方法，其特征在于，包括以下步骤：

步骤1，选取一岩石铸体薄片图片，框出明显的孔隙区域作为提取区域，提取所述提取区域的RGB值，分别两两进行交会，得到R与G交会图、R与B交会图和G与B交会图，手动标注以上三个交会图中条带的上限直线和下限直线，构建三个交会图中所述上限直线和下限直线的临界线方程，由此确定粗略的阈值界限；

步骤2，选取10张以上岩石铸体薄片图片，根据步骤1中获得的粗略的阈值界限提取每一岩石铸体薄片图片中的对应区域的RGB值，利用多元线性回归得到平面方程B(R，G)，选取图像中的孔隙误导区域记录其RGB值，修正所述平面方程B(R，G)，确定孔隙提取的精确阈值条件，利用该精确阈值条件确定岩石铸体薄片图片中的孔隙区域；所述步骤2中，获取所述孔隙误导区域的RGB值为B，计算L＝B-B(R，G)，统计得到最小的正L，记作L_up，统计得到最大的负L，记作L_down，所述精确阈值条件为：

B_up(R，G)＝B(R，G)+L_up

B_down(R，G)＝B(R，G)-L_down。

2.一种定量计算岩石面孔率的方法，其特征在于，包括以下步骤：

步骤1，选取一岩石铸体薄片图片，框出明显的孔隙区域作为提取区域，提取所述提取区域的RGB值，分别两两进行交会，得到R与G交会图、R与B交会图和G与B交会图，手动标注以上三个交会图中条带的上限直线和下限直线，构建三个交会图中所述上限直线和下限直线的临界线方程，由此确定粗略的阈值界限；

步骤2，选取10张以上岩石铸体薄片图片，根据步骤1中获得的粗略的阈值界限提取每一岩石铸体薄片图片中对应区域的RGB值，利用多元线性回归得到平面方程B(R，G)，选取图像中的孔隙误导区域记录其RGB值，修正所述平面方程B(R，G)，确定孔隙提取的精确阈值条件；

步骤3，输入待估面孔率的岩石铸体薄片图片，根据步骤2得到的精确阈值条件判断每个像素是否为孔隙；

步骤4，提取面孔并计算面孔率；

所述步骤2中，获取所述孔隙误导区域的RGB值为B，计算L＝B-B(R，G)，统计得到最小的正L，记作L_up，统计得到最大的负L，记作L_down，所述精确阈值条件为：

B_up(R，G)＝B(R，G)+L_up

B_down(R，G)＝B(R，G)-L_down。

3.如权利要求2所述的定量计算岩石面孔率的方法，其特征在于，

所述步骤1中的临界线方程为

其中x、y为R、G或B值，Li为三个交会图的临界线L1、L2、L3、L4、L5和L6，(x₁，y₁)、(x₂，y₂)为人为确定的临界线上的两个点；

所述阈值界限为：

f_L2(R)≤G≤f_L1(R) (2)

f_L4(R)≤B≤f_L3(R) (3)

f_L6(G)≤B≤f_L5(G) (4)。

4.如权利要求2所述的定量计算岩石面孔率的方法，其特征在于，提取RGB值时，利用铸体薄片孔隙提取软件ctsPore中孔隙初步提取模块或者利用软件Matlab中颜色阈值模块提取。