CN111462838B - 一种直接将图像像素转换成有限元单元的方法 - Google Patents

Abstract

本发明公开了一种直接将图像像素转换成有限元单元的方法，包括步骤：采用扫描电子显微镜获取研究对象的微观图像，并对图像进行图像处理以及进行不同相的划分；将图像的像素点矩阵坐标值转换成二维空间坐标系的点坐标值；将每个像素点定义成一个有限单元并计算其相应的节点坐标，对每个单元与节点坐标分别进行编号；将网格模型中每个单元的材料属性根据其图像像素颜色进行区分，并分别定义好相对应的单元属性；建立特殊相以及各个边界的单元和节点集，便于边界条件的施加与定义。本发明的方法适用于基于图像重建模型的有限元模拟分析，为图像模型有限元模拟领域提供一种直接将像素转换成单元与节点的方法。

Description

一种直接将图像像素转换成有限元单元的方法

技术领域

本发明涉及计算机模拟技术领域，尤其涉及一种直接将图像像素转换成有限元单元的方法。

背景技术

一直以来实验是学术研究者进行材料表征、结构性能测试等科学研究采用的最佳手段。实验带来的数据是直观、客观且可靠的，但是受实验环境、实验边界条件、材料均匀性、样品精度等各方面的因素影响，实验往往存在局限性，实验过程总是不可避免地出现实验误差，并且一些周期过长、重复次数多的实验需要耗费大量时间精力。于是研究者们采用计算机模拟计算来代替一些难以实现的实验以及验证一些难以评判可靠性的实验数据。其中有限元模拟是计算机数值计算的一类有效方法，有限元模拟需要先建立精确的结构模型，并确立符合实际的材料属性，定义精准的边界条件，才能获得有价值的仿真结果。

近年来，随着电子显微镜的发展，材料的微观结构可以由扫描电镜获取清晰的图像，而随着有限元技术的广泛应用，基于图像的微观结构重建模型及其有限元模拟研究正在日益增长。在许多学者的工作中，微观结构模型的重建与有限元模拟往往是区分为两个步骤，其中需要使用三维绘图软件，并且需要借助商用有限元软件进行有限单元模型的网格划分，对整个计算过程带来极大的不便利。于是本发明提供一种直接将图像像素转换成有限元单元的方法。

发明内容

本发明的目的在于克服现有技术的不足，提供一种直接将图像像素转换成有限元单元的方法。本发明可以直接将图像像素转换成有限单元模型，无需多余的软件手动操作，能够高效便捷地得到有限元模型以及完成计算所需要的参数定义，使其应用于有限元模拟计算。

本发明的目的能够通过以下技术方案实现：

一种直接将图像像素转换成有限元单元的方法，包括步骤：

采用扫描电子显微镜(SEM)获取研究对象的微观图像，并对图像进行图像处理以及进行不同相的划分；

将图像的像素点矩阵坐标值转换成二维空间坐标系的点坐标值；

将每个像素点定义成一个有限单元并计算其相应的节点坐标，对每个单元与节点坐标分别进行编号；

将网格模型中每个单元的材料属性根据其图像像素颜色进行区分，并分别定义好相对应的单元属性；

建立特殊相以及各个边界的单元和节点集，便于边界条件的施加与定义。

具体地，所述并对图像进行图像处理以及进行不同相的划分的步骤中，包括：

通过Matlab编程平台将SEM图像读取为RGB三维矩阵；

将SEM图像从彩色RGB图像格式转换成灰度图；

采用默认阈值将灰度图转换成二值图，即分为2个相，前景值为1，背景值为0；若模型包括三个及以上的相，则根据不同相的灰度值情况进行多阈值分割；

采用形态学处理对二值图进行优化，所述优化包括膨胀、开运算和闭运算等处理方法。

具体地，所述将图像的像素点矩阵坐标值转换成二维空间坐标系的点坐标值的步骤中，坐标转换公式为：

其中，(x₀,y₀)为大小为m×n的像素矩阵中某一个像素点(i,j)转换到平面坐标系中的坐标值。

具体地，所述将每个像素点定义成一个有限单元并计算其相应的节点坐标，对每个单元与节点坐标分别进行编号的步骤中，包括：

将每个像素点(i,j)转换成一个平面矩形四节点单元，并且像素点坐标(x₀,y₀)直接作为单元左下角节点坐标；

每个单元包括4个节点，从左下角开始进行逆时针编号为1～4，以1个像素点单位为节点间隔，则4个节点的坐标值表示为：

分别对单元与节点进行编号，并界定好每个单元对应的节点编号，即完成了图像矩阵像素点坐标转换成有限元单元节点坐标。

更进一步地，由于相邻的像素点在转换成相邻单元时会出现节点共享的情况，且节点定义中需要避免重复定义，因此在计算某单元的节点坐标以及对单元与节点的归属关系进行界定时，需要优先判断该单元相邻单元的节点是否已完成定义。

由于像素矩阵的扫描顺序是从左上角像素点开始，并按照从左往右，从上往下的顺序进行逐行扫描，因此节点定义的原则为：

基于扫描转换顺序，不管该单元的上邻与左邻是否已存在单元，右下角节点都未生成，优先根据像素点坐标值定义2号节点坐标值；

如果该单元上邻已有定义单元，即该单元的3号和4号节点坐标可以分别直接提取上邻单元的2号和1号节点坐标；

如果该单元左邻也有已定义单元，则该单元的1号节点坐标直接提取左邻单元的2号节点坐标；否则，则需要根据像素点坐标值定义1号节点的坐标值；

如果该单元上邻没有已定义单元，则需要根据像素点坐标值定义3号节点的坐标值；

如果左邻有已定义单元，则该单元的1号和4号节点坐标可以分别直接提取上邻单元的2号和3号节点坐标；否则，则需要根据像素点坐标值定义1、4号节点的坐标值。

具体地，所述将网格模型中每个单元的材料属性根据其图像像素颜色进行区分，并分别定义好相对应的单元属性的步骤中，包括：

在对SEM图像的灰度图进行图像处理时，对不同材料的区域进行相的区分；

给不同材料的区域赋予不同颜色，分别获取对应的二值图，建立不同相的单元模型；

根据图像中不同颜色对应的材料相对有限元网格模型进行不同材料属性的定义。

本发明相较于现有技术，具有以下的有益效果：

1、本发明提供的像素转换成有限元单元及节点的方法简单且有效，保证了网格连续性。

2、本发明中的有限元模型构建方法可以直接将图像转换成模型，并且不同材料相有对应的材料属性，可应用于多材料体系模型。

3、本发明最终得到的有限元模型可以直接应用于有限元模拟计算，转换步骤清晰、简单，无需第三方商用建模软件与有限元软件的参与，大大提高了计算智能化，有效简化了基于图像重建模型的有限元模拟计算方法。

附图说明

图1是本发明中一种直接将图像像素转换成有限元单元的方法的流程图；

图2是本发明实施例中图像像素转换成单元模型的示意图。

图3是本发明实施例中单元对应的节点编号示意图。

具体实施方式

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

实施例

在本实施例中，如图1所示为一种直接将图像像素转换成有限元单元的方法的流程图，所述方法包括步骤：

(1)采用扫描电子显微镜(SEM)获取研究对象的微观图像，并对图像进行图像处理以及进行不同相的划分，包括步骤：

(1-1)通过Matlab编程平台将SEM图像读取为RGB三维矩阵；

(1-2)将SEM图像从彩色RGB图像格式转换成灰度图SEM-Gray；

(1-3)采用默认阈值将灰度图SEM-Gray转换成二值图，前景值为1，背景值为0；

(1-4)采用形态学处理对二值图进行优化，所述优化包括膨胀、开运算和闭运算等处理方法。

如图2所示为图像像素转换成单元模型的示意图，从图2中可以看出图像在计算机中由像素矩阵组成，放大来看一个面是由点组成的。每个像素点可以转换成一个单元，相对应的每个单元由一定的节点组成。

(2)将图像的像素点矩阵坐标值转换成二维空间坐标系的点坐标值，具体转换公式如下：

其中，(x₀,y₀)为大小为m×n的像素矩阵中某一个像素点(i,j)转换到平面坐标系中的坐标值。

(3)将每个像素点定义成一个有限单元并计算其相应的节点坐标，对每个单元与节点坐标分别进行编号，包括步骤：

(3-1)将像素点Pn(i,j)转换成一个平面矩形四节点单元n，并且像素点坐标(x₀,y₀)直接作为单元左下角节点坐标；

(3-2)单元n包括4个节点，从左下角开始进行逆时针编号1～4，具体编号如图3所示，以1个单位为节点间隔，即4个节点的坐标值如下：

(3-3)分别对单元与节点进行编号，并界定好每个单元对应的节点编号，即完成了图像矩阵像素点坐标转换成有限元单元节点坐标。

更进一步地，由于相邻的像素点在转换成相邻单元时会出现节点共享的情况，且节点定义中需要避免重复定义，因此在计算某单元的节点坐标以及对单元与节点的归属关系进行界定时，需要优先判断该单元相邻单元的节点是否已完成定义。

由于像素矩阵的扫描顺序是从左上角像素点开始，并按照从左往右，从上往下的顺序进行逐行扫描，因此节点定义的原则为：

基于扫描转换顺序，不管该单元的上邻与左邻是否已存在单元，右下角节点都未生成，优先根据像素点坐标值定义2号节点坐标值；

如果该单元上邻已有定义单元，即该单元的3号和4号节点坐标可以分别直接提取上邻单元的2号和1号节点坐标；

如果该单元左邻也有已定义单元，则该单元的1号节点坐标直接提取左邻单元的2号节点坐标；否则，则需要根据像素点坐标值定义1号节点的坐标值；

如果该单元上邻没有已定义单元，则需要根据像素点坐标值定义3号节点的坐标值；

如果左邻有已定义单元，则该单元的1号和4号节点坐标可以分别直接提取上邻单元的2号和3号节点坐标；否则，则需要根据像素点坐标值定义1、4号节点的坐标值。

(4)网格模型中每个单元的材料属性根据其图像像素颜色进行区分，分别定义好相对应的单元属性，包括步骤：

(4-1)在对SEM图像的灰度图进行图像处理时，对不同材料的区域进行相的区分；

(4-2)给不同材料的区域赋予不同颜色，分别获取对应的二值图，建立不同相的单元模型；

(4-3)根据图像中不同颜色对应的材料相对有限元网格模型进行不同材料属性的定义。

(5)建立特殊相以及各个边界的单元、节点集，便于边界条件的施加与定义。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (3)

1.一种直接将图像像素转换成有限元单元的方法，其特征在于，包括步骤：

采用扫描电子显微镜获取研究对象的微观图像，并对图像进行图像处理以及进行不同相的划分；

所述并对图像进行图像处理以及进行不同相的划分的步骤中，包括：

通过Matlab编程平台将SEM图像读取为RGB三维矩阵；

将SEM图像从彩色RGB图像格式转换成灰度图；

采用默认阈值将灰度图转换成二值图，前景值为1，背景值为0；

采用形态学处理对二值图进行优化；

将图像的像素点矩阵坐标值转换成二维空间坐标系的点坐标值；

所述将图像的像素点矩阵坐标值转换成二维空间坐标系的点坐标值的步骤中，坐标转换公式为：

其中，(x₀,y₀)为大小为m×n的像素矩阵中某一个像素点(i,j)转换到平面坐标系中的坐标值；

将每个像素点定义成一个有限单元并计算其相应的节点坐标，对每个单元与节点坐标分别进行编号；

所述将每个像素点定义成一个有限单元并计算其相应的节点坐标，对每个单元与节点坐标分别进行编号的步骤中，包括：

将每个像素点(i,j)转换成一个平面矩形四节点单元，并且像素点坐标(x₀,y₀)直接作为单元左下角节点坐标；

每个单元包括4个节点，从左下角开始进行逆时针编号为1～4，以1个单位为节点间隔，则4个节点的坐标值表示为：

分别对单元与节点进行编号，并界定好每个单元对应的节点编号，即完成了图像矩阵像素点坐标转换成有限元单元节点坐标；

将有限元网格模型中每个单元的材料属性根据其图像像素颜色进行区分，并分别定义好相对应的单元属性；

所述将网格模型中每个单元的材料属性根据其图像像素颜色进行区分，并分别定义好相对应的单元属性的步骤中，包括：

在对SEM图像的灰度图进行图像处理时，对不同材料的区域进行相的区分；

给不同材料的区域赋予不同颜色，分别获取对应的二值图，建立不同相的单元模型；

根据图像中不同颜色对应的材料相对有限元网格模型进行不同材料属性的定义；

建立特殊相以及各个边界的单元和节点集。

2.根据权利要求1所述的方法，其特征在于，由于相邻的像素点在转换成相邻单元时会出现节点共享的情况，且节点定义中需要避免重复定义，因此在计算某单元的节点坐标以及对单元与节点的归属关系进行界定时，需要优先判断该单元相邻单元的节点是否已完成定义。

3.根据权利要求2所述的方法，其特征在于，由于像素矩阵的扫描顺序是从左上角像素点开始，并按照从左往右，从上往下的顺序进行逐行扫描，因此节点定义的原则为：

基于扫描转换顺序，不管该单元的上邻与左邻是否已存在单元，右下角节点都未生成，优先根据像素点坐标值定义2号节点坐标值；

如果该单元上邻已有定义单元，即该单元的3号和4号节点坐标可以分别直接提取上邻单元的2号和1号节点坐标；

如果该单元左邻也有已定义单元，则该单元的1号节点坐标直接提取左邻单元的2号节点坐标；否则，则需要根据像素点坐标值定义1号节点的坐标值；

如果该单元上邻没有已定义单元，则需要根据像素点坐标值定义3号节点的坐标值；

如果左邻有已定义单元，则该单元的1号和4号节点坐标可以分别直接提取上邻单元的2号和3号节点坐标；否则，则需要根据像素点坐标值定义1、4号节点的坐标值。

Images