CN110297003A - 一种基于轴、径环状分层的ct扫描后处理方法 - Google Patents
一种基于轴、径环状分层的ct扫描后处理方法 Download PDFInfo
- Publication number
- CN110297003A CN110297003A CN201910586551.3A CN201910586551A CN110297003A CN 110297003 A CN110297003 A CN 110297003A CN 201910586551 A CN201910586551 A CN 201910586551A CN 110297003 A CN110297003 A CN 110297003A
- Authority
- CN
- China
- Prior art keywords
- short column
- sample image
- section
- sample
- type
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/044—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using laminography or tomosynthesis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/03—Investigating materials by wave or particle radiation by transmission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/646—Specific applications or type of materials flaws, defects
- G01N2223/6466—Specific applications or type of materials flaws, defects flaws comparing to predetermined standards
Landscapes
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Theoretical Computer Science (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Pulmonology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开一种基于轴、径环状分层的CT扫描后处理方法,它包括以下步骤:步骤1):对标准试样进行CT扫描获取其CT扫描图像,并得到CT值;步骤2):将标准试样图像均分为多段短柱试样图像;步骤3):将每段短柱试样图像均分为多层环状短柱试样图像;步骤4):根据步骤1)中得到的CT扫描图像来获取经过步骤2)和步骤3)划分后的环状短柱试样图像CT值,然后计算试样的CT轴、径平均值;本发明能够精确地分析试样内部的不均匀性并准确确定受损范围,改善现有的CT扫描图像分析无法精细化的问题。
Description
技术领域
本发明涉及土木工程领域,具体地指一种基于轴、径环状分层的CT扫描后处理方法。
背景技术
近年来,国内外运用CT扫描对土木、水利等领域的结构进行分析的研究越来越重视。CT扫描主要运用于如下几方面:1、分析混凝土、岩石等试样中裂隙的分布;2、分析多场耦合作用后混凝土、岩石等试样的损伤规律;3、分析冻融循环后混凝土、岩石等试样的损伤机理。大量的工程实例和试验研究表明,裂隙的存在影响了结构的长期变形和工程的稳定性,多场耦合作用后所造成的结果原因探析不明,冻融循环之后试样内部损伤分析困难。故而CT扫描技术在土木、水利等领域的运用日益增多。
常规CT扫描主要有三种形式。一、在对试样扫描层面位置确定后,分析各层面的CT扫描图,对比扫描图形的亮度,图形越亮,表示试样致密性越好,中部的少数黑点表示有局部空隙,利用此方法虽然可以直观的观测出岩石内部的孔隙对比,大致得出岩石的损伤规律,却无法对其进行一个量化的分析,无法精确确定受损范围。二、层面上选取有限个圈取面内的CT平均值对整体进行分析,但在圈取面的选择时,往往无法划出具有代表性的范围,所以得出的结论与实际差异较大。三、取整个扫描层面的CT平均值作为扫描层面致密性的判断标准。但在该区域内,致密性往往不均匀,而以该CT值作为判断依据的话,无法反映出扫描层面的不均匀性。
为了能更精确的分析试样内部的不均匀性,进一步确定试样的重点受损范围,提出一种基于轴-径环状分层的CT扫描后处理方法。利用CT值与轴-径环状分层后形成的柱壳体积之间的关系得到每一个柱壳区域内的体积CT值,分析比较每一个柱壳的CT值变化以达到对试样进行损伤定量的动态分析。
发明内容
本发明的目的在于克服上述不足,提供一种基于轴、径环状分层的CT扫描后处理方法,以精确地分析试样内部的不均匀性并准确确定受损范围,改善现有的CT扫描图像分析无法精细化的问题。
本发明为解决上述技术问题,所采用的技术方案是:一种基于轴、径环状分层的CT扫描后处理方法,它包括以下步骤:
步骤1):对标准试样进行CT扫描获取其CT扫描图像,并得到CT值;
步骤2):将标准试样图像均分为多段短柱试样图像;
步骤3):将每段短柱试样图像均分为多层环状短柱试样图像;
步骤4):根据步骤1)中得到的CT扫描图像来获取经过步骤2)和步骤3)划分后的环状短柱试样图像CT值,然后计算试样的CT轴、径平均值。
优选地,所述步骤2)中,将标准试样图像均分为多段短柱试样图像的具体分段方法如下:
对于高度为L,直径为D的圆柱体试样,可径向将其均分为I段短柱试样,并且初步令每段短柱试样高度如果n不是整数,则处理过程如下:
n*=[n],其中[n]表示不超过n的最大整数,并令I=n*;
这时最终确定将每段短柱岩样高度变为h=L/I。
优选地,所述步骤3)中,将每段短柱试样图像均分为多层环状短柱试样图像的具体分层方法如下:
对于高度为L,直径为D的圆柱体试样,可将每段短柱试样的上表面或下表面均分为J个环状,初步令每个环状宽度为d,J=m*=[m],其中[m]表示不超过m的最大整数,其中r为试样划分时最内侧圆形半径;
这时最终确定将每层环状宽度变为
优选地,所述步骤2)中,将标准试样图像均分为多段短柱试样图像的具体分段方法如下:
对于高度为L,横截面边长均为D的长方体试样,可径向将其均分为I段短柱试样,并且初步令每段短柱试样高度如果n不是整数,则处理过程如下:
n*=[n],其中[n]表示不超过n的最大整数,并令I=n*;
这时最终确定将每段短柱岩样高度变为h=L/I。
优选地,所述步骤3)中,将每段短柱试样图像均分为多层环状短柱试样图像的具体分层方法如下:
对于高度为L,横截面边长均为D的长方体试样,可将每段短柱试样的上表面或下表面均分为J个环状,初步令每个环状宽度为d,J=m*=[m],其中[m]表示不超过m的最大整数,其中r为试样划分时最内侧正方形边长的一半;
这时最终确定将每层环状宽度变为
优选地,所述步骤5)中,计算试样的CT轴、径平均值具体方法如下:
分别获取每个划分后的环状短柱试样图像上、下表面圆环的CT值,CTSij表示第i段、第j层环状短柱试样图像的上表面圆环外侧CT值与内侧CT值之差,CTXij表示第i段、第j层环状短柱试样图像的下表面圆环外侧CT值与内侧CT值之差;Vij表示第i段、第j层环状短柱试样图像体积,则对于其中任意一个试样存在以下计算公式:
本发明的有益效果:
1、本发明首次提出了一种基于轴、径环状分层的CT扫描后处理方法,将微分的思想运用于CT图像处理,可以精确分析试样在某种状态下的物理特性的空间分布规律,同时不再局限于通过肉眼识别图像明暗度以及平面内特定区域的平均CT值进行分析,借助数学公式将CT值转化为单元柱壳(即环状短柱试样)的体积CT值,实现了精确量化分析。
2、本发明可对试样进行轴、径单元划分,得到若干单元柱壳,通过柱壳的体积CT值计算,得到各柱壳单元的体积CT值,再通过对比各柱壳单元的体积CT值,得到各柱壳单元物理特性变化程度,进一步缩小试样的特性变化范围。
3、本方法不仅可适用于岩石试样,同时还可适用于混凝土试样以及其他工程试验中的试样,而且对试样的形状规格要求不高,常见的圆柱状、立方体等形态试样均可适用,故而本发明具有高度的普适性。
4、本发明中对于试样CT图像的划分可有效避免任意划分带来的离散性,提高了测验结果的准确性
5、本发明相比于常规的CT扫描图像分析方法,其操作简便,分析结果可信度高。
附图说明
图1为本发明所涉及的圆柱体分段和分层示意图;
图2为图1的俯视图;
图3为图1中涉及的每个单元环状短柱试样图像;
图4为本发明所涉及的长方体分段和分层示意图;
图5为图4的俯视图;
图6为图4中涉及的每个单元环状短柱试样图像;
图7为本发明所涉及的常规CT扫描实图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细描述。
如图1所示,一种基于轴、径环状分层的CT扫描后处理方法,它包括以下步骤:
步骤1):对标准试样进行CT扫描获取其CT扫描图像,并得到CT值;
步骤2):将标准试样图像均分为多段短柱试样图像;
步骤3):将每段短柱试样图像均分为多层环状短柱试样图像;
步骤4):根据步骤1)中得到的CT扫描图像来获取经过步骤2)和步骤3)划分后的环状短柱试样图像CT值,然后计算试样的CT轴、径平均值。
实施例1:如图1、2和3所示,对于高度为L,直径为D的圆柱体试样,所述步骤2)中,将标准试样图像均分为多段短柱试样图像的具体分段方法如下:
可径向将其均分为I段短柱试样(本实施例中I=5),并且初步令每段短柱试样高度如果n不是整数,则处理过程如下:
n*=[n],其中[n]表示不超过n的最大整数,并令I=n*;
这时最终确定将每段短柱岩样高度变为h=L/I。
将每段短柱试样图像均分为多层环状短柱试样图像的具体分层方法如下:
可将每段短柱试样的上表面或下表面均分为J个环状(本实施例中J=2),初步令每个环状宽度为d,J=m*=[m],其中[m]表示不超过m的最大整数,其中r为试样划分时最内侧圆形半径;
这时最终确定将每层环状宽度变为
所述步骤5)中,计算试样的CT轴、径平均值具体方法如下:
分别获取每个划分后的环状短柱试样图像上、下表面圆环的CT值,CTSij表示第i段、第j层环状短柱试样图像的上表面圆环外侧CT值与内侧CT值之差,CTXij表示第i段、第j层环状短柱试样图像的下表面圆环外侧CT值与内侧CT值之差;Vij表示第i段、第j层环状短柱试样图像体积,则对于其中任意一个试样存在以下计算公式:
实施例2:如图4、5和6所示,对于高度为L,横截面边长均为D的长方体试样,所述步骤2)中,将标准试样图像均分为多段短柱试样图像的具体分段方法如下:
可径向将其均分为I段短柱试样(本实施例中I=5),并且初步令每段短柱试样高度如果n不是整数,则处理过程如下:
n*=[n],其中[n]表示不超过n的最大整数,并令I=n*;
这时最终确定将每段短柱岩样高度变为h=L/I。
将每段短柱试样图像均分为多层环状短柱试样图像的具体分层方法如下:
可将每段短柱试样的上表面或下表面均分为J个环状(本实施例中J=4),初步令每个环状宽度为d,J=m*=[m],其中[m]表示不超过m的最大整数,其中r为试样划分时最内侧正方形边长的一半;
这时最终确定将每层环状宽度变为
所述步骤5)中,计算试样的CT轴、径平均值具体方法如下:
分别获取每个划分后的环状短柱试样图像上、下表面圆环的CT值,CTSij表示第i段、第j层环状短柱试样图像的上表面圆环外侧CT值与内侧CT值之差,CTXij表示第i段、第j层环状短柱试样图像的下表面圆环外侧CT值与内侧CT值之差;Vij表示第i段、第j层环状短柱试样图像体积,则对于其中任意一个试样存在以下计算公式:
上述公式中在计算第i段、第j层环状短柱试样图像的体积Vij时,可以通过第j层外侧圆柱体或长方体体积减去内侧的圆柱体或长方体体积即得到环形体积,以高度为L,直径为D的圆柱体试样为例,其Vij的计算公式如下:
vij=h{π(r+jd)2-π[(j-1)d+r]2}
另外需要说明的是,在对试样CT图像划分后,其最里面的一个圆柱体或长方体的CT值计算方法并不适用上述公式,可直接将其上、下表面的CT值相加除以其体积即可,不必取其内外侧CT值之差。
上述的实施例仅为本发明的优选技术方案,而不应视为对于本发明的限制,本申请中的实施例及实施例中的特征在不冲突的情况下,可以相互任意组合。本发明的保护范围应以权利要求记载的技术方案,包括权利要求记载的技术方案中技术特征的等同替换方案为保护范围。即在此范围内的等同替换改进,也在本发明的保护范围之内。
Claims (6)
1.一种基于轴、径环状分层的CT扫描后处理方法,其特征在于:它包括以下步骤:
步骤1):对标准试样进行CT扫描获取其CT扫描图像,并得到CT值;
步骤2):将标准试样图像均分为多段短柱试样图像;
步骤3):将每段短柱试样图像均分为多层环状短柱试样图像;
步骤4):根据步骤1)中得到的CT扫描图像来获取经过步骤2)和步骤3)划分后的环状短柱试样图像CT值,然后计算试样的CT轴、径平均值。
2.根据权利要求1所述的一种基于轴、径环状分层的CT扫描后处理方法,其特征在于:所述步骤2)中,将标准试样图像均分为多段短柱试样图像的具体分段方法如下:
对于高度为L,直径为D的圆柱体试样,可径向将其均分为I段短柱试样,并且初步令每段短柱试样高度如果n不是整数,则处理过程如下:
n*=[n],其中[n]表示不超过n的最大整数,并令I=n*;
这时最终确定将每段短柱岩样高度变为h=L/I。
3.根据权利要求1所述的一种基于轴、径环状分层的CT扫描后处理方法,其特征在于:所述步骤3)中,将每段短柱试样图像均分为多层环状短柱试样图像的具体分层方法如下:
对于高度为L,直径为D的圆柱体试样,可将每段短柱试样的上表面或下表面均分为J个环状,初步令每个环状宽度为d,J=m*=[m],其中[m]表示不超过m的最大整数,其中r为试样划分时最内侧圆形半径;
这时最终确定将每层环状宽度变为
4.根据权利要求1所述的一种基于轴、径环状分层的CT扫描后处理方法,其特征在于:所述步骤2)中,将标准试样图像均分为多段短柱试样图像的具体分段方法如下:
对于高度为L,横截面边长均为D的长方体试样,可径向将其均分为I段短柱试样,并且初步令每段短柱试样高度如果n不是整数,则处理过程如下:
n*=[n],其中[n]表示不超过n的最大整数,并令I=n*;
这时最终确定将每段短柱岩样高度变为h=L/I。
5.根据权利要求1所述的一种基于轴、径环状分层的CT扫描后处理方法,其特征在于:所述步骤3)中,将每段短柱试样图像均分为多层环状短柱试样图像的具体分层方法如下:
对于高度为L,横截面边长均为D的长方体试样,可将每段短柱试样的上表面或下表面均分为J个环状,初步令每个环状宽度为d,J=m*=[m],其中[m]表示不超过m的最大整数,其中r为试样划分时最内侧正方形边长的一半;
这时最终确定将每层环状宽度变为
6.根据权利要求1所述的一种基于轴、径环状分层的CT扫描后处理方法,其特征在于:所述步骤5)中,计算试样的CT轴、径平均值具体方法如下:
分别获取每个划分后的环状短柱试样图像上、下表面圆环的CT值,CTSij表示第i段、第j层环状短柱试样图像的上表面圆环外侧CT值与内侧CT值之差,CTXij表示第i段、第j层环状短柱试样图像的下表面圆环外侧CT值与内侧CT值之差;Vij表示第i段、第j层环状短柱试样图像体积,则对于其中任意一个试样存在以下计算公式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910586551.3A CN110297003A (zh) | 2019-07-01 | 2019-07-01 | 一种基于轴、径环状分层的ct扫描后处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910586551.3A CN110297003A (zh) | 2019-07-01 | 2019-07-01 | 一种基于轴、径环状分层的ct扫描后处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110297003A true CN110297003A (zh) | 2019-10-01 |
Family
ID=68029784
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910586551.3A Pending CN110297003A (zh) | 2019-07-01 | 2019-07-01 | 一种基于轴、径环状分层的ct扫描后处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110297003A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117848921A (zh) * | 2023-11-30 | 2024-04-09 | 中国人民解放军陆军工程大学 | 一种基于细观模型构建的岩土介质含气量、饱和度测量及孕灾、减灾评价方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101871861A (zh) * | 2010-06-07 | 2010-10-27 | 中国水电顾问集团华东勘测设计研究院 | 深部岩样初始损伤范围及程度测定方法 |
CN104266862A (zh) * | 2014-10-24 | 2015-01-07 | 武汉大学 | 一种基于ct扫描的深部岩石取样方法 |
CN106228558A (zh) * | 2016-07-26 | 2016-12-14 | 三峡大学 | 一种对岩石进行层进式损伤分析的方法 |
CN109816655A (zh) * | 2019-02-01 | 2019-05-28 | 华院数据技术(上海)有限公司 | 基于ct图像的肺结节图像特征检测方法 |
-
2019
- 2019-07-01 CN CN201910586551.3A patent/CN110297003A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101871861A (zh) * | 2010-06-07 | 2010-10-27 | 中国水电顾问集团华东勘测设计研究院 | 深部岩样初始损伤范围及程度测定方法 |
CN104266862A (zh) * | 2014-10-24 | 2015-01-07 | 武汉大学 | 一种基于ct扫描的深部岩石取样方法 |
CN106228558A (zh) * | 2016-07-26 | 2016-12-14 | 三峡大学 | 一种对岩石进行层进式损伤分析的方法 |
CN109816655A (zh) * | 2019-02-01 | 2019-05-28 | 华院数据技术(上海)有限公司 | 基于ct图像的肺结节图像特征检测方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117848921A (zh) * | 2023-11-30 | 2024-04-09 | 中国人民解放军陆军工程大学 | 一种基于细观模型构建的岩土介质含气量、饱和度测量及孕灾、减灾评价方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101539502B (zh) | 一种固体表面液滴接触角的测量方法及装置 | |
CN111337992B (zh) | 一种基于位场数据向下延拓的场源深度获得方法 | |
CN105806861A (zh) | 一种基于虚拟单轴贯入试验的沥青混合料均匀性评价方法 | |
CN105631876A (zh) | 一种基于全局二值化的ct图像分辨率自动测试方法 | |
JPWO2010041423A1 (ja) | 病理組織診断支援システム、病理組織診断支援プログラム、病理組織診断支援方法 | |
CN113158490B (zh) | 致密砂岩储层渗透率计算模型的建立方法 | |
CN110297003A (zh) | 一种基于轴、径环状分层的ct扫描后处理方法 | |
CN105114067A (zh) | 一种岩性测井相的方法 | |
CN111027117A (zh) | 胶凝砂砾石配合比关键指标预测抗压强度的bp神经网络分析方法 | |
CN106228558B (zh) | 一种对岩石进行层进式损伤分析的方法 | |
CN107923830A (zh) | 材料的断裂表面的开裂机制的表征方法 | |
CN102621221B (zh) | 一种基于相控超声波的缺陷分类方法 | |
CN106127227B (zh) | 采用拉曼光谱对物质进行聚类分析识别的方法 | |
Liu et al. | Computer Vision‐Based Structural Displacement Monitoring and Modal Identification with Subpixel Localization Refinement | |
Tovar-Valencia et al. | Short-term setup of jacked piles in a calibration chamber | |
CN109738955B (zh) | 一种基于成分-结构分类下的变质岩岩性综合判别方法 | |
CN110895704B (zh) | 微生物丘滩复合体储集层类型识别方法、装置及存储介质 | |
CN110108800A (zh) | 基于超声相控阵图像碳纤维复合材料缺陷分类方法和装置 | |
CN114609163B (zh) | 基于色差识别的岩石微观结构射线成像探测法 | |
CN106548474A (zh) | 一种微结构表面检测方法 | |
CN115128244A (zh) | 一种考虑岩石流体作用力的注二氧化碳封存潜力分析方法 | |
CN108872391A (zh) | 用于评价岩体稳定状态的物探分析方法 | |
Wang et al. | An image-based soil type classification method considering the impact of image acquisition distance factor | |
CN112070110B (zh) | 一种致密储层微观孔喉结构分级模式的预测方法 | |
CN113253354A (zh) | 一种核磁t2谱截止值的确定方法、装置、介质及电子设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191001 |