CN102768699A - 基于ct图像精确重构异质材料微观有限元网格模型的方法 - Google Patents
基于ct图像精确重构异质材料微观有限元网格模型的方法 Download PDFInfo
- Publication number
- CN102768699A CN102768699A CN2012101964058A CN201210196405A CN102768699A CN 102768699 A CN102768699 A CN 102768699A CN 2012101964058 A CN2012101964058 A CN 2012101964058A CN 201210196405 A CN201210196405 A CN 201210196405A CN 102768699 A CN102768699 A CN 102768699A
- Authority
- CN
- China
- Prior art keywords
- image
- pixel
- model
- finite element
- sequence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000000463 material Substances 0.000 title claims abstract description 104
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000002591 computed tomography Methods 0.000 title abstract 7
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 230000011218 segmentation Effects 0.000 claims abstract description 11
- 238000004458 analytical method Methods 0.000 claims abstract description 10
- 238000005520 cutting process Methods 0.000 claims description 20
- 238000013507 mapping Methods 0.000 claims description 14
- 230000001186 cumulative effect Effects 0.000 claims description 12
- 238000005315 distribution function Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 7
- 230000002950 deficient Effects 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 5
- 238000013519 translation Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 claims description 4
- 230000008034 disappearance Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 5
- 238000009825 accumulation Methods 0.000 abstract description 2
- 230000003044 adaptive effect Effects 0.000 abstract 1
- 238000005457 optimization Methods 0.000 abstract 1
- 238000005192 partition Methods 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 14
- 238000003672 processing method Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000003709 image segmentation Methods 0.000 description 2
- 240000007762 Ficus drupacea Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 210000003918 fraction a Anatomy 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 239000005445 natural material Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000000059 patterning Methods 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 229920002994 synthetic fiber Polymers 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明提出了一种利基于CT图像精确重构异质材料微观有限元网格模型的方法,该方法通过工业CT采集序列断层图像,并在数字化和阈值分割的基础上将断层图像中的微结构信息映射到重构的有限元网格模型中,使得异质材料中的任何细节结构信息都可以在重构模型中得到重现。本发明通过对比度受限自适应直方均衡化、中值滤波和像素点插值提高重构精度,通过图像裁剪和像素点合并提高重构效率。本发明直接重构出有限元分析精度较高的矩形(二维)和长方体(三维)单元网格模型,避免了现有几何重构方法在重构和网格划分等环节中的误差积累,提高了重构精度和效率。本发明可广泛应用于异质材料的性能预测及优化设计等领域。
Description
技术领域
本发明属于有限元网格模型重构技术领域,涉及异质材料有限元网格模型重构方法,尤其是一种基于CT图像精确重构异质材料微观有限元网格模型的方法。
背景技术
由多种组分构成的异质材料(如复合材料、多孔材料等)是天然材料及合成材料中最为常见的材料之一,这些材料的宏观性能(如刚度、强度和韧度等)主要由其微观结构决定。因此,透彻地研究异质材料的微观结构对其宏观性能的影响对于设计和开发新的高性能异质材料具有重要意义。有限元方法是研究材料微观结构与宏观性能之间关系的最有效的方法之一,该方法需要首先建立能够反映材料真实微观结构的模型。然而,由于异质材料微观结构的复杂性,建立其微结构模型是材料设计的一大难题。
在异质材料性能的研究中,人们提出了多种克服建模困难的方法,其中,比较传统的是单胞模型方法。该方法基于异质材料在宏观上表现出的均质特性,将其微观结构理想化为具有简单几何特征的模型。这种方法能够有效地预测异质材料的组分性能、体积分数、形状和分布对宏观性能的影响,但是该方法忽略了异质材料真实的微观几何构造,难以用于对材料进行精确的分析。
随着数字图像处理技术的不断发展,研究人员提出了利用异质材料断层扫描图像重构其微观几何模型的方法。该方法采集断面图像的基础上,通过将图像分割为不同的组分相并进而提取各相边界信息的方法建立材料的几何表面或实体模型。该几何重构方法考虑了异质材料固有的形态学、聚集态和分布特征,并且尽可能少地引入微结构假设,这些特点使得该方法被广泛地用于分析异质材料的宏观行为。
尽管近年来几何重构方法趋于完善,对异质材料进行宏观性能的有限元预测仍然需要发展新的微观结构重构方法。因为采用现有几何重构方法建立的模型需要进一步划分网格,而基于复杂表面或实体划分的网格多为三角形或四面体单元,有限元分析精度不高;另外,异质材料微观结构的复杂性限制了网格的质量,且网格剖分过程中容易丢失材料微观结构中相对细小的信息。这些缺陷使得几何重构方法难以更加广泛地应用于材料的研究设计中,因而发展能够克服这些缺陷的新型重构方法显得非常重要。
发明内容
本发明的目的在于克服上述现有技术的缺点,提供一种基于CT图像精确重构异质材料微观有限元网格模型的方法,该方法通过工业CT采集材料的序列断层图像,并将断层图像中的微结构信息映射到重构的有限元网格模型中,使得异质材料中的任何细节结构信息都可以在模型中得到重现;为了提高重构的精确度,引入了对比度受限自适应直方均衡化、中值滤波和像素点插值等图像处理方法;另外,提出了图像裁剪和像素点合并等方法,用于提高后期有限元分析效率。
本发明的目的是通过以下技术方案来解决的:
这种基于CT图像精确重构异质材料微观有限元网格模型的方法,包括以下步骤:
1)利用工业CT采集异质材料中各平行截面的序列断层图像,并将图像数字化;
2)对图像进行阈值分割,并通过对比度受限自适应直方均衡化、中值滤波和像素点插值提高图像反映材料真实微观结构的精度;
3)建立与序列断层图像形成映射关系的有限元网格拓扑模型,利用图像中的组分相灰度信息确定有限元模型中各个单元的材料属性,并通过图像裁剪和像素点合并来减小有限元网格模型的单元规模。
上述步骤1)具体按照以下方法进行:
通过工业CT机在试件表面选取一个宽为w,高为h的能够代表异质材料微观结构特征的代表性矩形区域I1,采集该区域图像;然后将CT机的扫描位置沿试件的轴向平移微小距离dn,其中n=1~D-1;并采集新的图像In+1,得到总层数为D的一系列断层图像;断层图像用像素点矩形表示,且矩阵沿i和j方向的大小分别记为W和H,故单个像素点的尺寸为s=w/W=h/H;像素点(i,j)的灰度用离散函数g(i,j)表示,故单张断层图像可以表示为灰度矩阵[g(i,j)](i=1~W,j=1~H);离散函数g(i,j,k)用于表示序列断层图像中像素点(i,j,k)的灰度,故序列断层图像表示为三维数组{g(i,j,k)}(k=1~D)。
上述步骤2)具体按照以下方法进行:
(1)图像阈值分割;引入统计图像中给定灰度值的像素点数目的灰度直方图h(g):
对于具有M(M≥2)相的异质材料,灰度直方图具有M-1个波谷,将这些波谷对应的灰度值记为灰度阈值Tn(n=1,2,…,M-1),则异质材料中各组分相由这一系列灰度阈值来分割,即图像中灰度值低于阈值T1的像素点表示密度最低的组分相,灰度值在区间(T1,T2]内的像素点表示密度第二低的组分相,依此类推便可将断层图像完全地分割为M个组分相;各组分相的体积分数由以下累积分布函数f(g)来确定:
(2)对比度受限自适应直方均衡化;采用对比度受限自适应直方均衡化降低环境及仪器误差对断层图像的影响,从而提高图像反映异质材料真实微观结构的精度;
(3)中值滤波;对断层图像进行中值滤波处理,以降低图像噪声;
(4)像素点插值;为了能够精确地重构材料沿z轴方向的微结构特征,序列图像的层间距要尽可能地接近图像像素点的尺寸,即层间距需要满足下式:
E(d)+σ(d)≤ξs其中,
参数ξ依照精度要求而定,推荐该值的取值范围为1.0≤ξ≤2.0;如果E(d)+σ(d)超过了上限ξs,序列图像便需要通过像素点插值的方法来提高重构精度;原始图像In和In+1之间插值的图像层数Dn以及新的序列图像之间的层间距由下式确定:
将插值形成的新图像和原始图像按照他们的z坐标进行排序,则插值图像各像素点的灰度值由下式确定:
g(i,j,k)=P(Zk,g(i,j,n)),k∈Ni,n∈Na
(11)
式中,Ni和Na是分别指向插值新图像和原始图像序号的集合,Zk为新序列图像中第k张图像的z坐标,P(Z,g)是以(Zn,g(i,j,n)),n∈Na为插值点的插值函数。
上述步骤3)具体按照以下方法进行:
(1)二维重构;建立与单层图像各像素点形成映射的二维有限元网格拓扑模型,该有限元模型由四节点的矩形单元构成,并且单元尺寸与图像像素点的尺寸s相等,模型中第n个节点的坐标由下式确定:
第n个单元的材料属性编号为:
(2)三维重构;建立与序列断层图像像素点形成映射关系的三维有限元网格模型,该模型由八节点的长方体单元构成,单元沿x和y方向的尺寸与图像像素点的尺寸s相等,沿z方向的尺寸与序列图像的层间距相等;模型中第n个节点的坐标为:
模型中第n个单元沿逆时针方向排列的八个节点由下式确定:
模型中第n个单元的材料属性编号由下式确定:
(3)图像裁剪及像素点合并;重构得到的模型中包含的单元数目非常大,使得后续有限元分析非常耗时;由于断层图像中靠近图像边缘的像素点受到的噪声干扰最多,故采用裁剪图像周界像素点的方式来减小重构模型的单元规模;然而,裁剪的像素点数目越多,图像中包含的材料聚集团簇的数目也就越少,使得该方法的应用具有一定的局限性;这个缺陷可以通过像素点合并来予以克服,即通过下式来合并图像中相邻的像素点:
本发明具有以下有益效果:
(1)本发明直接精确地重构出反映异质材料真实微观结构的有限元网格模型,避免了现有几何方法在重构和划分网格等环节中的误差积累,节约了重构时间,提高了重构效率;
(2)本发明采用了能够使断层图像更加精确地反映材料真实微观结构的图像处理方法,并且重构出的网格分别为有限元分析精度较高的矩形(二维)和长方体(三维)单元,相对于现有几何重构方法,大大提高了有限元分析的精度;
(3)本发明通过采用能够减小有限元网格模型单元规模的图像处理方法,有效地平衡了有限元分析精度与计算时间之间的矛盾,在保证分析精度的前提下提高了计算效率;
(4)本发明重构的有限元网格模型可以方便地用于分析异质材料中各组分相的形态、聚集特征以及分布等微观结构与宏观性能之间的关系,从而为材料的开发设计提供依据。
附图说明
图1为本发明的流程图;
图2为序列断层图像的采集;
图3为图像的数字化;
图4为与单张断层图像形成映射关系的二维有限元网格模型的建立;
图5为二维有限元网格模型的节点和单元编号;
图6为与序列断层图像形成映射关系的三维有限元网格模型的建立;
图7为三维有限元网格模型的节点和单元编号;
图8为某两相异质材料的单张断层图像;
图9为图8中图像的灰度直方图和累积分布函数;
图10为经过CLAHE和MF处理后的单张断层图像;
图11为图10中图像的灰度直方图和累积分布函数;
图12为对某两相异质材料重构出的二维有限元网格模型;
图13为对图10中图像进行的裁剪处理及其阈值分割结果;
图14为对图10中图像进行的像素点合并处理及其阈值分割结果;
图15为经图像裁剪和像素点合并处理之后重构得到的二维有限元网格模型;
图16为某两相异质材料的序列断层图像;
图17为图16中序列断层图像的灰度直方图和累积分布函数;
图18为对某两相异质材料重构出的三维有限元网格模型;
图19为经像素点合并处理之后重构得到的三维有限元网格模型。
具体实施方式
本发明首先通过工业CT采集异质材料的序列断层图像,并将采集的图像数字化;然后通过灰度直方图和累积分布函数对图像进行阈值分割,并通过对比度受限自适应直方均衡化、中值滤波和像素点插值提高图像反映材料真实微观结构的精度;进而在建立有限元网格拓扑模型的基础上,通过拓扑模型与断层图像之间的映射关系确定模型的单元材料属性,并采用图像裁剪和像素点合并减小有限元网格模型中的单元规模。整个有限元网格模型重构的具体实施流程如图1所示,下面根据该流程详细地描述具体技术问题。
1.序列图像采集及数字化
本发明通过工业CT采集异质材料的序列断层图像,采集过程如图2所示。首先将打磨抛光后的异质材料试件固定在工业CT机的转台上,调节转台及平移***,在试件表面选取一个宽为w,高为h的能够代表异质材料微观结构特征的代表性矩形区域I1。该区域的选择不仅要使采集的图像中包含足够多的相聚集簇,同时又要限制图像中相聚集簇的数目以保证较高的图像分辨率和较快的重构速度。确定代表性矩形区域I1后,用CT机采集该区域的图像,然后将CT机的扫描位置沿试件的轴向平移微小距离d1(注意,在扫描位置的平移过程中必须保证在x和y方向不产生位移),采集新的矩形区域I2的图像,重复以上过程便可采集总层数为D的一系列断层图像。上述过程中,相邻两层图像之间的距离dn(n=1~D-1)要尽可能接近图像像素点的尺寸,以使序列断层图像能够精确地反映材料在z方向的微结构特征。另外,层数D的选择需要保证由序列图像重构出的模型在z方向包含足够数目的组分相聚集簇。
每一张断层图像都是以矩形排列的一个个像素点构成的,因此这些图像可以在图3所示的像素点坐标中以矩阵的形式予以表示。如图,图像矩阵沿i和j方向的大小分别记为W和H,每个像素点在宽度和高度方向的尺寸相等,并记为s=w/W=h/H。在每个像素点中,存储着反映图像视觉效果的信息,这些信息由名为灰度的整数表示,因此像素点(i,j)的灰度可以用离散函数g(i,j)予以表示,单张断层图像可以表示为灰度矩阵[g(i,j)](i=1~W,j=1~H)。对于序列断层图像,我们建立像素点左手坐标系下的k坐标,用以记录这一系列图像的序号,则离散函数g(i,j,k)可用于表示像素点(i,j,k)的灰度,三维数组{g(i,j,k)}(k=1~D)可用于表示序列断层图像。
2.图像阈值分割及图像处理
(1)图像阈值分割
工业CT采集的断层图像中各像素点的灰度值主要由异质材料各组分相的密度决定,即组分相密度越高,其吸收的X射线越多,该相占据的图像区域中的灰度值也相应地越高。根据这一规律,图像中的不同组分相可以由图像的灰度值予以区分。为了分割图像中不同的组分相,引入灰度直方图h(g),该直方图通过统计图像中给定灰度值的像素点数目来确定,即:
对所有的灰度值进行遍历统计便得到图像的灰度直方图。
对于具有M(M≥2)相的异质材料,其不同的组分相具有不同的密度,从而在CT图像中由不同的灰度值表征。灰度直方图具有M-1个波谷,这些波谷便是各相材料的分界点。将M-1个波谷处对应的灰度值记为灰度阈值Tn(n=1,2,…,M-1),则异质材料中各组分相由这一系列灰度阈值来分割,即图像中灰度值低于阈值T1的像素点表示密度最低的组分相,灰度值在区间(T1,T2]内的像素点表示密度第二低的组分相,依此类推便可将断层图像完全地分割为M个组分相。
将断层图像分割之后,各组分相的体积分数便可以由累积分布函数f(g)来确定。累积分布函数由以下经过像素点总数目归一化处理的公式计算,即:
(2)对比度受限自适应直方均衡化
断层图像的采集容易受到环境及仪器误差等一系列因素的影响,因此采集到的图像不能精确地反映异质材料真实的微观结构,而这一缺陷可以通过对比度增强和降噪等图像处理方法得以改善。
对比度增强有助于重现异质材料的真实微观结构,因为该方法能够增强图像中组分材料之间的对比度。本发明采用对比度增强系列方法中应用最为广泛的对比度受限自适应直方均衡化(CLAHE),该方法的数学描述详见Reza的综述文章(A.M.Reza,J VLSI Sig Proc Syst,38(2004)35-44)。
(3)中值滤波
除了对比度增强,降噪也能提高图像反映材料真实微观结构的精度。图像降噪方法很多,包括均值滤波、中值滤波和自适应滤波等,本发明选用其中降噪效果较好的中值滤波(MF)。对于特定形式的随机噪声,中值滤波能够在尽可能少地模糊边界和其他细节的情况下提供优异的降噪性能。中值滤波的理论和实现方法详见Gonzalez关于图像处理的书籍(R.C.Gonzalez,R.E.Woods,Digital Image Processing,Second Ed.,Prentice Hall,New Jersey,2001)。
(4)像素点插值
对于序列断层图像,为了能够精确地重构材料沿z轴方向的微结构特征,相邻图像之间的距离要尽可能地接近图像像素点的尺寸,也就是说,层间距需要满足下式:
E(d)+σ(d)≤ξs (20)
式中,E(d)和σ(d)分别为层间距的数学期望和标准差,即
参数ξ依照具体的精度要求而定,但是本发明推荐该值的取值范围为1.0≤ξ≤2.0。
公式(22)表明只要大多数的层间距不超过规定的上限ξs,序列图像便可直接用于重构异质材料的微观模型。相反,如果E(d)+σ(d)超过了上限ξs,序列图像便需要通过像素点插值的方法来提高重构精度。在像素点插值的过程中,原始图像之间将产生一系列新的图像,而这些图像用于重现原始图像之间的微结构细节。原始图像In和In+1之间插值的图像层数Dn以及新的序列图像之间的层间距由下式确定:
插值形成的新图像和原始图像需要按照他们的z坐标进行重新排序,并且新序列图像的记号与原始图像保持一致。插值图像各像素点的灰度值由以下插值公式确定:
g(i,j,k)=P(Zk,g(i,j,n)),k∈Ni,n∈Na (24)
式中,Ni和Na是分别指向插值新图像和原始图像序号的集合,Zk为新序列图像中第k张图像的z坐标,P(Z,g)是以(Zn,g(i,j,n)),n∈Na为插值点的插值函数。
3.有限元网格模型重构
(1)二维重构
本发明可以根据序列断层图像中的任何一张图像In来重构异质材料的二维真实微观结构。建立与单层图像各像素点形成映射的二维有限元网格拓扑模型,如图4所示。该有限元模型由四节点的矩形单元构成,并且单元沿x和y方向的尺寸与图像像素点的尺寸s相等。为了便于建立有限元模型与图像像素点之间简单的数学映射关系,对有限元模型节点和单元编号如图5所示,则模型第n个节点的坐标由下式确定:
式中,xn和yn分别为节点n沿x和y方向的坐标,“%”为整数除法取余,取小于运算对象的最大整数。模型第n个单元逆时针方向排列的四个节点由下式确定:
在有限元数据文件中,不同的单元材料属性由其属性编号确定,且属性编号用从1开始的自然数予以表示,编号越小,代表单元的材料密度越小。从而根据阈值分割后的断层图像和映射关系可确定第n个单元的材料属性编号为:
其中,i和j与n之间的关系为:
(2)三维重构
建立与序列断层图像像素点形成映射关系的三维有限元网格模型,如图6所示。该有限元网格模型由八节点的长方体单元构成,单元沿x和y方向的尺寸与图像像素点的尺寸s相等,沿z方向的尺寸与序列图像的层间距相等。对有限元网格模型的节点和单元编号如图7所示,则模型第n个节点的坐标为:
模型中第n个单元沿逆时针方向排列的八个节点由下式确定:
与二维重构类似,模型中第n个单元的材料属性编号可由下式确定:
其中,i、j和k与n之间的关系为
(3)图像裁剪及像素点合并
根据以上二维和三维重构方法建立的有限元网格模型能够保留异质材料详细而真实的微结构特征,但是,重构的模型中包含的单元数目非常大,而且随着图像分辨率的提高,模型中将包含更多的单元,这将使得后续有限元分析非常耗时。由于断层图像中靠近图像边缘的像素点受到的噪声干扰最多,故可以通过裁剪图像周界像素点的方式来减小重构模型的单元规模。
然而,裁剪的像素点数目越多,图像中包含的材料聚集团簇的数目也就越少,使得该方法的应用具有一定的局限性。这个缺陷可以通过像素点合并来予以克服,即通过下式来合并图像中相邻的像素点:
4.实施例
(1)二维重构
首先利用工业CT机采集如图8所示的某种两相异质材料的单张断层图像,其基体相材料的密度低于增强相材料的密度,故图像中趋近黑色的部分为基体相材料,趋近白色的为增强相材料。该图像i和j方向的像素点数目分别为W=606和H=597,x和y方向的实际尺寸分别为w=16.12μm和h=15.88μm,从而每个像素点的尺寸为s=0.0266μm。图9给出了该图像的灰度直方图和累积分布函数,由图可知,灰度阈值T1=154,基体相的体积分数为实验测得材料基体相的体积分数为55.00%,即通过图像分割得到的结果具有2.64%的相对误差。为了减小建模误差,对断层图像实施以下四种图像处理:
(a)CLAHE处理,处理后T1=140,v1=54.92%,相对误差为0.08%;
(b)MF处理,处理后T1=154,v1=53.63%,相对误差为2.49%;
(c)先CLAHE后MF处理,处理后T1=140,v1=54.98%,相对误差为0.04%;
(d)先MF后CLAHE处理,处理后T1=140,v1=54.54%,相对误差为0.84%。
由以上这一系列图像处理结果来看,先对断层图像进行CLAHE处理,然后进行MF处理时,图像能够最为精确地反映异质材料的真实微观结构。经过这种处理之后的断层图像如图10所示,其灰度直方图和累积分布函数如图11所示。利用本发明提出的二维重构方法对处理后的图像进行重构,得到如图12所示的二维有限元网格模型。
由于重构的有限元模型中单元数目过多,对图像进行进一步的图像裁剪和图像合并处理。首先,对图像实施如图13所示的一系列图像裁剪处理,由结果可知,处理后的图像的相对误差均不超过2.79%。图像裁剪的程度可以在平衡重构精度和计算效率的基础上予以确定,如果要求相对误差不超过1.00%,那么较好平衡计算效率的处理方法是裁减掉图像28个像素点的周界,在这种情况下,图像的相对误差为0.68%。
图像裁剪之后单元数目仍然过多,即单一的图像裁剪无法适当地减小模型单元规模。尝试对图像开展如图14所示的一系列像素点合并处理,由结果可知,处理后的图像相对误差均不超过1.94%,仍然可以通过衡量重构精度和计算效率来确定像素点合并的数目。
经过一系列试凑发现,首先裁剪掉图像28个像素点的周界,然后对图像进行5个像素点的合并处理能够较好地平衡重构精度和计算效率。经过处理后W=110,H=108,T1=140,v1=54.81%,相对误差为0.352%。由该图像重构得到的二维有限元网格模型如图15所示。
(2)三维重构
利用工业CT机采集如图16所示的上述两相异质材料的序列断层图像,其W=H=256,w=h=6.8096μm,s=0.0266μm,D=11,dn=0.158μm。为了提高各层图像重现材料微观结构的精度,首先对序列图像进行CLAHE处理,然后对其进行MF处理。处理之后的图像的灰度直方图和累积分布函数如图17所示,由图可知,T1=139,v1=54.66%,相对误差为0.62%。
取参数ξ为1.0,则E(d)+σ(d)超过了上限ξs=0.0266μm,故先对序列图像进行自然样条插值,得到新的序列图像中包含61层图像,且新序列图像相邻两层之间的距离为0.0263μm,这足以重构材料沿z方向的微结构细节信息。由新序列图像重构得到的三维有限元网格模型如图18所示,重构时T1=139,v1=55.21%%,相对误差为0.38%。
重构得到的三维模型中包含过多的单元,需要对其进行处理以减小单元规模。由于该序列图像沿i和j方向的像素点数目不多,这里只对其进行像素点合并处理,且处理过程中沿两个方向只合并两个像素点。经过图像处理之后重构出的三维有限元网格模型如图19所示,其相对误差为0.49%。
上述建模重构实例是针对两相异质材料而进行的,但是本发明所提出的技术并不仅限于两相异质材料,由前文述及的技术原理可知,本发明可以容易地对任意多相的异质材料进行有限元网格重构。
Claims (4)
1.一种基于CT图像精确重构异质材料微观有限元网格模型的方法,其特征在于,包括以下步骤:
1)利用工业CT采集异质材料中各平行截面的序列断层图像,并将图像数字化;
2)对图像进行阈值分割,并通过对比度受限自适应直方均衡化、中值滤波和像素点插值提高图像反映材料真实微观结构的精度;
3)建立与序列断层图像形成映射关系的有限元网格拓扑模型,利用图像中的组分相灰度信息确定有限元模型中各个单元的材料属性,并通过图像裁剪和像素点合并来减小有限元网格模型的单元规模。
2.根据权利要求1所述的基于CT图像精确重构异质材料微观有限元网格模型的方法,其特征在于,步骤1)具体按照以下方法进行:
通过工业CT机在试件表面选取一个宽为w,高为h的能够代表异质材料微观结构特征的代表性矩形区域I1,采集该区域图像;然后将CT机的扫描位置沿试件的轴向平移微小距离dn,其中n=1~D-1;并采集新的图像In+1,得到总层数为D的一系列断层图像;断层图像用像素点矩形表示,且矩阵沿i和j方向的大小分别记为W和H,故单个像素点的尺寸为s=w/W=h/H;像素点(i,j)的灰度用离散函数g(i,j)表示,故单张断层图像可以表示为灰度矩阵[g(i,j)](i=1~W,j=1~H);离散函数g(i,j,k)用于表示序列断层图像中像素点(i,j,k)的灰度,故序列断层图像表示为三维数组{g(i,j,k)}(k=1~D)。
3.根据权利要求1所述的基于CT图像精确重构异质材料微观有限元网格模型的方法,其特征在于,步骤2)按照以下方法进行:
(1)图像阈值分割;引入统计图像中给定灰度值的像素点数目的灰度直方图h(g):
对于具有M(M≥2)相的异质材料,灰度直方图具有M-1个波谷,将这些波谷对应的灰度值记为灰度阈值Tn(n=1,2,…,M-1),则异质材料中各组分相由这一系列灰度阈值来分割,即图像中灰度值低于阈值T1的像素点表示密度最低的组分相,灰度值在区间(T1,T2]内的像素点表示密度第二低的组分相,依此类推便可将断层图像完全地分割为M个组分相;各组分相的体积分数由以下累积分布函数f(g)来确定:
(2)对比度受限自适应直方均衡化;采用对比度受限自适应直方均衡化降低环境及仪器误差对断层图像的影响,从而提高图像反映异质材料真实微观结构的精度;
(3)中值滤波;对断层图像进行中值滤波处理,以降低图像噪声;
(4)像素点插值;为了能够精确地重构材料沿z轴方向的微结构特征,序列图像的层间距要尽可能地接近图像像素点的尺寸,即层间距需要满足下式:
E(d)+σ(d)≤ξs其中,
将插值形成的新图像和原始图像按照他们的z坐标进行排序,则插值图像各像素点的灰度值由下式确定:
g(i,j,k)=P(Zk,g(i,j,n)),k∈Ni,n∈Na (1)
式中,Ni和Na是分别指向插值新图像和原始图像序号的集合,Zk为新序列图像中第k张图像的z坐标,P(Z,g)是以(Zn,g(i,j,n)),n∈Na为插值点的插值函数。
4.根据权利要求1所述的基于CT图像精确重构异质材料微观有限元网格模型的方法,其特征在于,步骤3)按照以下方法进行:
(1)二维重构;建立与单层图像各像素点形成映射的二维有限元网格拓扑模型,该有限元模型由四节点的矩形单元构成,并且单元尺寸与图像像素点的尺寸s相等,模型中第n个节点的坐标由下式确定:
第n个单元的材料属性编号为:
(2)三维重构;建立与序列断层图像像素点形成映射关系的三维有限元网格模型,该模型由八节点的长方体单元构成,单元沿x和y方向的尺寸与图像像素点的尺寸s相等,沿z方向的尺寸与序列图像的层间距相等;模型中第n个节点的坐标为:
模型中第n个单元沿逆时针方向排列的八个节点由下式确定:
模型中第n个单元的材料属性编号由下式确定:
(3)图像裁剪及像素点合并;重构得到的模型中包含的单元数目非常大,使得后续有限元分析非常耗时;由于断层图像中靠近图像边缘的像素点受到的噪声干扰最多,故采用裁剪图像周界像素点的方式来减小重构模型的单元规模;然而,裁剪的像素点数目越多,图像中包含的材料聚集团簇的数目也就越少,使得该方法的应用具有一定的局限性;这个缺陷可以通过像素点合并来予以克服,即通过下式来合并图像中相邻的像素点:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210196405.8A CN102768699B (zh) | 2012-06-14 | 2012-06-14 | 基于ct图像精确重构异质材料微观有限元网格模型的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210196405.8A CN102768699B (zh) | 2012-06-14 | 2012-06-14 | 基于ct图像精确重构异质材料微观有限元网格模型的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102768699A true CN102768699A (zh) | 2012-11-07 |
CN102768699B CN102768699B (zh) | 2016-08-24 |
Family
ID=47096100
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210196405.8A Expired - Fee Related CN102768699B (zh) | 2012-06-14 | 2012-06-14 | 基于ct图像精确重构异质材料微观有限元网格模型的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102768699B (zh) |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903303A (zh) * | 2012-12-27 | 2014-07-02 | 清华大学 | 三维模型创建方法和设备 |
CN104063902A (zh) * | 2014-06-22 | 2014-09-24 | 湘潭大学 | 一种基于材料真实微观组织结构的有限元建模方法 |
CN104184957A (zh) * | 2014-09-17 | 2014-12-03 | 中国科学院光电技术研究所 | 一种适用于空间探测成像的基于自适应预期图像平均亮度的自动曝光控制方法 |
ES2530424A1 (es) * | 2013-09-02 | 2015-03-02 | Quibim S L | Método de segmentación por umbral adaptativo variable para la obtención de valores de referencia del aire corte a corte en estudios de imagen por tomografía computarizada |
CN104952108A (zh) * | 2015-05-20 | 2015-09-30 | 中国矿业大学(北京) | 一种ct逆向化建模技术的网格模型优化方法 |
CN106469236A (zh) * | 2015-08-21 | 2017-03-01 | 波音公司 | 利用映射到结构上的不一致而建模的结构的分析 |
CN106682366A (zh) * | 2017-02-04 | 2017-05-17 | 哈尔滨工大泰铭科技有限公司 | 一种基于图像处理的有限元节点坐标提取方法 |
CN107392992A (zh) * | 2017-06-26 | 2017-11-24 | 北京农业信息技术研究中心 | 玉米茎节维管束三维建模方法 |
CN107391855A (zh) * | 2017-07-26 | 2017-11-24 | 华中科技大学 | 一种面向多种微观结构的材料结构一体化构建方法 |
CN107766640A (zh) * | 2017-10-16 | 2018-03-06 | 北京理工大学 | 考虑微观组织界面的颗粒增强复合材料有限元建模方法 |
CN108021770A (zh) * | 2017-12-04 | 2018-05-11 | 北京理工大学 | 基于ct扫描的叶片寿命评价方法 |
CN108053478A (zh) * | 2017-11-10 | 2018-05-18 | 北京理工大学 | 一种基于像素理论的颗粒增强复合材料有限元建模方法 |
CN108088864A (zh) * | 2017-12-15 | 2018-05-29 | 浙江隆劲电池科技有限公司 | 一种材料三维微观结构重构方法及*** |
CN108387524A (zh) * | 2018-01-22 | 2018-08-10 | 青岛理工大学 | 一种多相材料连续剖面图像获取方法 |
CN109087396A (zh) * | 2018-07-26 | 2018-12-25 | 东南大学 | 一种基于混凝土ct图像像素特征的细观结构重构方法 |
CN109596618A (zh) * | 2018-12-26 | 2019-04-09 | 太原理工大学 | 一种基于断面轮廓序列的微观多相结构三维建模测量方法 |
CN110245626A (zh) * | 2019-06-19 | 2019-09-17 | 北京万里红科技股份有限公司 | 一种在虹膜图像中精确检测眼睫毛图像的方法 |
CN110796693A (zh) * | 2019-09-11 | 2020-02-14 | 重庆大学 | 一种工业ct切片图像直接生成二维有限元模型的方法 |
CN111047687A (zh) * | 2019-12-18 | 2020-04-21 | 浙江大学 | 一种基于三维t样条的异质材料实体建模方法 |
CN111508075A (zh) * | 2020-04-09 | 2020-08-07 | 北京理工大学 | 一种闭孔泡沫铝的三维真实有限元模型建模方法 |
CN112613208A (zh) * | 2020-12-16 | 2021-04-06 | 武汉科技大学 | 基于三维随机重构的无机非金属材料性能分析方法 |
CN113570544A (zh) * | 2021-05-31 | 2021-10-29 | 清华大学 | 松质骨建模方法、装置、存储介质及电子设备 |
CN114693660A (zh) * | 2022-04-11 | 2022-07-01 | 中国人民解放军空军工程大学 | 一种基于ict的固体火箭发动机装药计算网格生成方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007045000A2 (en) * | 2005-10-14 | 2007-04-19 | Vantus Technology Corporation | Personal fit medical implants and orthopedic surgical instruments and methods for making |
CN102236737A (zh) * | 2011-07-14 | 2011-11-09 | 西安交通大学 | 一种基于序列图像的多相材料细观构造有限元重构方法 |
-
2012
- 2012-06-14 CN CN201210196405.8A patent/CN102768699B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007045000A2 (en) * | 2005-10-14 | 2007-04-19 | Vantus Technology Corporation | Personal fit medical implants and orthopedic surgical instruments and methods for making |
CN102236737A (zh) * | 2011-07-14 | 2011-11-09 | 西安交通大学 | 一种基于序列图像的多相材料细观构造有限元重构方法 |
Non-Patent Citations (1)
Title |
---|
MING HUANG等: "X-ray tomography image-based reconstruction of microstructural finite element mesh models for heterogeneous materials", 《COMPUTATIONAL MATERIALS SCIENCE》, vol. 67, 28 February 2013 (2013-02-28), pages 63 - 72 * |
Cited By (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903303A (zh) * | 2012-12-27 | 2014-07-02 | 清华大学 | 三维模型创建方法和设备 |
ES2530424A1 (es) * | 2013-09-02 | 2015-03-02 | Quibim S L | Método de segmentación por umbral adaptativo variable para la obtención de valores de referencia del aire corte a corte en estudios de imagen por tomografía computarizada |
CN104063902A (zh) * | 2014-06-22 | 2014-09-24 | 湘潭大学 | 一种基于材料真实微观组织结构的有限元建模方法 |
CN104184957A (zh) * | 2014-09-17 | 2014-12-03 | 中国科学院光电技术研究所 | 一种适用于空间探测成像的基于自适应预期图像平均亮度的自动曝光控制方法 |
CN104952108A (zh) * | 2015-05-20 | 2015-09-30 | 中国矿业大学(北京) | 一种ct逆向化建模技术的网格模型优化方法 |
CN106469236B (zh) * | 2015-08-21 | 2021-07-20 | 波音公司 | 利用映射到结构上的不一致而建模的结构的分析 |
CN106469236A (zh) * | 2015-08-21 | 2017-03-01 | 波音公司 | 利用映射到结构上的不一致而建模的结构的分析 |
CN106682366A (zh) * | 2017-02-04 | 2017-05-17 | 哈尔滨工大泰铭科技有限公司 | 一种基于图像处理的有限元节点坐标提取方法 |
CN106682366B (zh) * | 2017-02-04 | 2020-04-28 | 哈尔滨工大泰铭科技有限公司 | 一种基于图像处理的有限元节点坐标提取方法 |
CN107392992A (zh) * | 2017-06-26 | 2017-11-24 | 北京农业信息技术研究中心 | 玉米茎节维管束三维建模方法 |
CN107391855B (zh) * | 2017-07-26 | 2018-03-09 | 华中科技大学 | 一种面向多种微观结构的材料结构一体化构建方法 |
CN107391855A (zh) * | 2017-07-26 | 2017-11-24 | 华中科技大学 | 一种面向多种微观结构的材料结构一体化构建方法 |
CN107766640A (zh) * | 2017-10-16 | 2018-03-06 | 北京理工大学 | 考虑微观组织界面的颗粒增强复合材料有限元建模方法 |
CN107766640B (zh) * | 2017-10-16 | 2020-06-30 | 北京理工大学 | 考虑微观组织界面的颗粒增强复合材料有限元建模方法 |
CN108053478A (zh) * | 2017-11-10 | 2018-05-18 | 北京理工大学 | 一种基于像素理论的颗粒增强复合材料有限元建模方法 |
CN108053478B (zh) * | 2017-11-10 | 2021-02-09 | 北京理工大学 | 一种基于像素理论的颗粒增强复合材料有限元建模方法 |
CN108021770A (zh) * | 2017-12-04 | 2018-05-11 | 北京理工大学 | 基于ct扫描的叶片寿命评价方法 |
CN108088864A (zh) * | 2017-12-15 | 2018-05-29 | 浙江隆劲电池科技有限公司 | 一种材料三维微观结构重构方法及*** |
CN108387524A (zh) * | 2018-01-22 | 2018-08-10 | 青岛理工大学 | 一种多相材料连续剖面图像获取方法 |
CN108387524B (zh) * | 2018-01-22 | 2020-06-05 | 青岛理工大学 | 一种多相材料连续剖面图像获取方法 |
CN109087396A (zh) * | 2018-07-26 | 2018-12-25 | 东南大学 | 一种基于混凝土ct图像像素特征的细观结构重构方法 |
CN109087396B (zh) * | 2018-07-26 | 2023-09-29 | 东南大学 | 一种基于混凝土ct图像像素特征的细观结构重构方法 |
CN109596618A (zh) * | 2018-12-26 | 2019-04-09 | 太原理工大学 | 一种基于断面轮廓序列的微观多相结构三维建模测量方法 |
CN109596618B (zh) * | 2018-12-26 | 2021-02-26 | 太原理工大学 | 一种基于断面轮廓序列的微观多相结构三维建模测量方法 |
CN110245626B (zh) * | 2019-06-19 | 2021-06-22 | 北京万里红科技股份有限公司 | 一种在虹膜图像中精确检测眼睫毛图像的方法 |
CN110245626A (zh) * | 2019-06-19 | 2019-09-17 | 北京万里红科技股份有限公司 | 一种在虹膜图像中精确检测眼睫毛图像的方法 |
CN110796693A (zh) * | 2019-09-11 | 2020-02-14 | 重庆大学 | 一种工业ct切片图像直接生成二维有限元模型的方法 |
CN110796693B (zh) * | 2019-09-11 | 2023-03-21 | 重庆大学 | 一种工业ct切片图像直接生成二维有限元模型的方法 |
CN111047687A (zh) * | 2019-12-18 | 2020-04-21 | 浙江大学 | 一种基于三维t样条的异质材料实体建模方法 |
CN111047687B (zh) * | 2019-12-18 | 2021-10-22 | 浙江大学 | 一种基于三维t样条的异质材料实体建模方法 |
CN111508075A (zh) * | 2020-04-09 | 2020-08-07 | 北京理工大学 | 一种闭孔泡沫铝的三维真实有限元模型建模方法 |
CN111508075B (zh) * | 2020-04-09 | 2022-08-19 | 北京理工大学 | 一种闭孔泡沫铝的三维真实有限元模型建模方法 |
CN112613208A (zh) * | 2020-12-16 | 2021-04-06 | 武汉科技大学 | 基于三维随机重构的无机非金属材料性能分析方法 |
CN112613208B (zh) * | 2020-12-16 | 2022-09-13 | 武汉科技大学 | 基于三维随机重构的无机非金属材料性能分析方法 |
CN113570544A (zh) * | 2021-05-31 | 2021-10-29 | 清华大学 | 松质骨建模方法、装置、存储介质及电子设备 |
CN114693660A (zh) * | 2022-04-11 | 2022-07-01 | 中国人民解放军空军工程大学 | 一种基于ict的固体火箭发动机装药计算网格生成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102768699B (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102768699A (zh) | 基于ct图像精确重构异质材料微观有限元网格模型的方法 | |
CN103426169B (zh) | 一种医学图像的分割方法 | |
Géraud et al. | A quasi-linear algorithm to compute the tree of shapes of n D images | |
CN102236737B (zh) | 一种基于序列图像的多相材料细观构造有限元重构方法 | |
Halley et al. | Uses and abuses of fractal methodology in ecology | |
Goodchild et al. | Development and test of an error model for categorical data | |
CN112669463B (zh) | 三维点云的曲面重建方法、计算机设备和计算机可读存储介质 | |
Huang et al. | X-ray tomography image-based reconstruction of microstructural finite element mesh models for heterogeneous materials | |
CN106600683A (zh) | 一种面向骨骼ct序列图像的网格模型自适应重建方法 | |
CN103136535A (zh) | 一种用于点云简化的k近邻搜索方法 | |
CN110442994A (zh) | 一种基于图论的3d打印切片处理方法 | |
Barequet et al. | Straight-skeleton based contour interpolation | |
CN106126918A (zh) | 一种基于相互作用力的地理空间异常聚集区扫描统计方法 | |
Henderson et al. | Using geometry to uncover relationships between isotropy, homogeneity, and modularity in cortical connectivity | |
Zhou et al. | A streaming framework for seamless building reconstruction from large-scale aerial lidar data | |
CN107610028A (zh) | 一种基于无线云传感网的大气污染物在线监测*** | |
Liu et al. | Developing a scene-based triangulated irregular network (TIN) technique for individual tree crown reconstruction with LiDAR data | |
Lecigne et al. | Annual shoot segmentation and physiological age classification from TLS data in trees with acrotonic growth | |
CN108287940B (zh) | 样点布设的方法、装置、电子设备和存储介质 | |
CN110415339A (zh) | 计算输入三维形体间的匹配关系的方法和装置 | |
CN110967738A (zh) | 基于深度学习的地震数据插值处理方法及*** | |
CN110222742B (zh) | 基于分层多回波的点云分割方法、装置、存储介质及设备 | |
Kim et al. | Digital elevation modeling via curvature interpolation for LiDAR data | |
CN104766091A (zh) | 空间与光谱协同结构及遥感影像的多尺度理解方法、*** | |
Schultze et al. | Space carving and filtered back-projection as preconditioners for proton computed tomography reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160824 |
|
CF01 | Termination of patent right due to non-payment of annual fee |