[发明专利]一种基于三维欧氏距离的岩石空间结构距离图的计算方法

权利要求书

1.一种基于三维欧氏距离的岩石空间结构距离图的计算方法，其特征在于，包括以下步骤，

S1.利用高精度扫描仪器，获得待测岩石的三维图像，并采用分水岭算法对图像进行二值分割，得到岩石三维二值图像；

S2.以x-z,y-z平面为基准，将原始三维图像均等分为多块，并将每一块沿x轴顺序编号，按编号分别将编号对应的图像读取入不同计算机，采用多台计算机并行运算；

S3.对每个块图像的所有孔隙点的像素值赋无穷大，固体点为0；

S4.按编号顺序提取一个图像块，设置距离游标及扫描方向，生成空间距离游标，记录每一个目标点距离最近的背景点的地址, 同时利用目标点到每个邻居的距离和它每个邻居到最近背景点的距离可以得到当前这个点对应的欧式距离值；

S5.判断当前编号是否为最大值，若是，进行下一步，若不是重复步骤S4；

S6.对编号相邻的块进行融合，具体表现为通过一个三维空间距离游标在接触面上滑动，如果某个面的值发生了改变，那么对应内部所有距离值重复S4-S6；

S7.对所有相邻块均进行融合后，便可以得到岩石内部空间结构的距离图；

两个块进行融合的实质是通过面的融合来进行的，即将两个接触面再加一层赋值无穷大数值的面拼接在一起后，用空间游标进行正反两次遍历，判断新融合的块中是否有距离原来块上点更近的固体点，如果有，那么更新整个块内部，如果没有，那么块内部值保持不变。

2.根据权利要求1所述的一种基于三维欧氏距离的岩石空间结构距离图的计算方法，其特征在于，所述步骤S4中，若游标中有体素数值为“0”的点，计算中心体素到所有数值为“0”体素点的距离，取其中最小值替换初始所设的无穷大。

3.根据权利要求2所述的一种基于三维欧氏距离的岩石空间结构距离图的计算方法，其特征在于，所述步骤S4中，若游标中没有体素数值为“0”的点，被测点的欧式距离=邻域中欧式距离最小值+被测点到该邻居点的距离+两倍的两距离向量之积，并更新数据体中所有像素点的欧式距离值。

4.根据权利要求1-3中任一项所述的一种基于三维欧氏距离的岩石空间结构距离图的计算方法，其特征在于，所述步骤S4中，设置与首次扫描方向相反的反方向再扫描，检验被测点的值是否为距离最近固体点的距离值，若是，保留数据，反之取距离最近固体点的距离值更新数据体，完成单次欧式距离变换。