CN114813798A - 用于表征材料内部结构及成分的ct检测装置和成像方法 - Google Patents
用于表征材料内部结构及成分的ct检测装置和成像方法 Download PDFInfo
- Publication number
- CN114813798A CN114813798A CN202210541988.7A CN202210541988A CN114813798A CN 114813798 A CN114813798 A CN 114813798A CN 202210541988 A CN202210541988 A CN 202210541988A CN 114813798 A CN114813798 A CN 114813798A
- Authority
- CN
- China
- Prior art keywords
- axis
- flat
- photon
- shaft
- detector
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 67
- 238000001514 detection method Methods 0.000 title claims abstract description 44
- 239000000463 material Substances 0.000 title claims abstract description 20
- 238000012512 characterization method Methods 0.000 claims abstract description 14
- 238000013507 mapping Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 4
- 239000003086 colorant Substances 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000000007 visual effect Effects 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 3
- 239000000523 sample Substances 0.000 description 17
- 238000005516 engineering process Methods 0.000 description 11
- 230000004927 fusion Effects 0.000 description 11
- 238000000034 method Methods 0.000 description 11
- 238000001228 spectrum Methods 0.000 description 7
- 238000013170 computed tomography imaging Methods 0.000 description 4
- 239000013078 crystal Substances 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000001066 destructive effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 229910018072 Al 2 O 3 Inorganic materials 0.000 description 1
- 229910004298 SiO 2 Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009658 destructive testing Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000004255 ion exchange chromatography Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004811 liquid chromatography Methods 0.000 description 1
- 238000001819 mass spectrum Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000000918 plasma mass spectrometry Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
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/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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20092—Interactive image processing based on input by user
- G06T2207/20104—Interactive definition of region of interest [ROI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Pulmonology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种用于表征材料内部结构及成分的CT检测装置和成像方法,底座中央设有旋转轴,旋转轴上设有被测物端X轴,被测物端X轴上设有被测物端Y轴,被测物端Y轴上设有被测物端Z轴,被测物端Z轴上设有样品台;旋转轴一侧设有射线源端X轴,射线源端X轴上设有射线源端Y轴,射线源端Y轴上设有射线源端Z轴,射线源端Z轴上设有射线源;旋转轴另一侧设有探测器端X轴,探测器端X轴上设有探测器端Y轴,探测器端Y轴上设有探测器端Z轴和平板探测器,探测器端Z轴上设有光子计数探测器;本发明利用光子计数探测器能量分辨成像的优势,结合图像映射模型,弥补CT检测装置对材料成分表征的不足,实现等效原子序数成像。
Description
技术领域
本发明涉及测试技术领域,尤其涉及一种用于表征材料内部结构及成分的CT检测装置和成像方法。
背景技术
目前主要采用液相色谱、离子色谱、等离子体质谱、红外、XPS等技术对材料成分、纯度、杂质残留进行表征分析,而这些技术均为微量、有损型测试,无法满足宏观状态下产品的无损成分分析需求,因此,针对目前材料无损测试分析的科学技术问题,探索宏观、无损的成分表征方法具有重要意义。
X射线CT成像技术由于对材料进行宏观、无损的三维结构表征的优势,被广泛用于医疗、安检、工业无损检测等领域,但目前CT技术由于射线多色性以及成像伪影的原因,对材料内部的成分及密度测量还存在困难。近年来,随着X射线光子计数探测器技术的问世,促使新型X射线能谱技术应运而生,其能将每个入射光子按能量所属区间分别计数,通过设置电子学阈值消除暗电流噪声,提高图像信噪比,减少图像伪影,同时基于能量分辨能力,利用多种X射线能量成像,经处理可获得被扫描物体的原子序数等信息,实现高精度密度均匀性检测。将能谱技术与常规CT技术的优势相结合,可实现材料宏观、无损的内部结构与成分融合成像。
当前能谱CT技术的研究主要集中于扫描成像、数据采集、图像重建等方面,没有形成能谱CT技术的系统研究,因此,需要对能谱CT融合成像系统进行详细的装置设计、成像应用方法研究。
发明内容
本发明的目的就在于为了解决上述问题而提供一种用于表征材料内部结构及成分的CT检测装置和成像方法,本发明的CT检测装置包括光子计数探测器、平板探测器、射线源、机械运动机构,能够同时对检测样品进行无损结构及成分表征分析。
本发明通过以下技术方案来实现上述目的:
一种用于表征材料内部结构及成分的CT检测装置,包括底座,所述底座中央设有旋转轴,所述旋转轴上设有被测物端X轴,所述被测物端X轴上设有被测物端Y轴,所述被测物端Y轴上设有被测物端Z轴,所述被测物端Z轴上设有样品台;
所述旋转轴一侧设有射线源端X轴,所述射线源端X轴上设有射线源端Y轴,所述射线源端Y轴上设有射线源端Z轴,所述射线源端Z轴上设有射线源;
所述旋转轴另一侧设有探测器端X轴,所述探测器端X轴上设有探测器端Y轴,所述探测器端Y轴上设有探测器端Z轴和平板探测器,所述探测器端Z轴上设有光子计数探测器;
所述被测物端X轴、射线源端X轴、探测器端X轴在所述底座上的投影在一条直线上。
进一步方案为,所述平板探测器用于检测样品内部结构成像,所述光子探测器用于检测样品内部成分成像。
本发明另一方面还提供了上述装置的成像方法,包括多能量段成像、切片重建、感兴趣区域扫描、等效原子序数表征,对被测物进行内部结构及成分表征,获得被测物内部结构和成分的融合成像结果。
进一步方案为,所述多能量段成像依靠光子计数探测器实现,对被测物的衰减投影信号进行能量阈值分段成像,将不易被平板探测器区分的信息在多能量段投影结合运算,实现直观分析评判。
进一步方案为,所述感兴趣区域扫描利用平板探测器对检测样品进行全局DR成像,通过DR图像感兴趣区域选取,结合成像坐标系与实际坐标系转化,利用机械运动实现感兴趣区域的光子计数探测成像;
所述成像坐标系与实际坐标系转化公式为:
(ΔXobj,ΔYobj)=μflat*(Xflat-roi-c-Xflat-c,Yflat-roi-c-Yflat-c)/λflat
(ΔXobj,ΔYobj)=μphoton*(Xphoton-roi-c-Xphoton-c,Yphoton-roi-c-Yphoton-c)/λphoton
式中(ΔXobj,ΔYobj)为实际被测物坐标系中感兴趣区域中心与被测物中心的偏移量,λflat和λphoton分别为平板探测器和光子计数探测器的成像几何放大比,μflat和μphoton分别为平板探测器和光子计数探测器的单像素实际尺寸,(Xflat-c,Yflat-c)为平板DR投影图像中物体的中心坐标,(Xflat-roi-c,Yflat-roi-c)为平板DR投影图像中选取的感兴趣区域中心坐标,同理,(Xphoton-c,Yphoton-c)和(Xphoton-roi-c,Yphoton-roi-c)分别为光子计数DR投影图像的中心坐标和感兴趣区域中心坐标;
在系统几何参数自动校准后,(Xflat-c,Yflat-c)与(Xphoton-c,Yphoton-c)的实际空间位置一致,依靠探测器端的X轴进行放大比调节,保证λflat=λphoton,光子计数探测平面位置利用探测器端的Y轴与Z轴调整,调整量为μflat*(Xflat-roi-c-Xflat-c,Yflat-roi-c-Yflat-c);若由于机械行程或探测条件导致λflat≠λphoton,则调整量为μflat*(Xflat-roi-c-Xflat-c,Yflat-roi-c-Yflat-c)*λphoton/λflat。
进一步方案为,所述等效原子序数表征利用光子对检测样品进行多能量段投影采集,通过切片重建将多能量段投影转为多能量段切片,结合等效原子序数计算公式进行映射,获得检测样品等效原子序数切片图像。
所述等效原子序数计算公式为:
其中代表相除切片图像,代表低能段切片图像,为高能段切片图像,高、低能衰减系数可分解为光电效应与康普顿散射,fp(E)和fc(E)分别表示光电效应、康普顿散射与入射能量的关系,n取值3.94,在能量确定的情况下,将相除图像和等效原子序数Z按三个参数的指数模型进行关系建立。
进一步方案为,将平板或光子计数探测器单能量段投影获取内部结构切片图像,并获取内部等效原子序数切片图像,在内部结构切片图像上将等效原子序数切片图像以不同颜色进行绘制,图像中心像素点为对齐基准,获得内部结构和成分的融合成像结果;
对于λflat≠λphoton的情况,等效原子序数切片图像需要进行λphoton/λflat倍缩放,在平板探测器获得的内部结构切片图像上,以图像中心像素点为对齐基准,绘制缩放后的等效原子序数切片图像,舍弃超过内部结构切片图像矩阵索引范围的像素点。
本发明的有益效果在于:
本发明的一种用于表征材料内部结构及成分的CT检测装置和成像方法,通过CT扫描进行样品的成分表征,弥补液谱、质谱、红外、XPS等微量、有损型测试技术的不足,实现材料宏观、无损的内部结构与成分融合成像。检测装置主要包括光子计数探测器、平板探测器、射线源、机械运动机构,平板探测器对检测样品内部结构进行三维成像,光子计数探测器对检测样品关键位置成分状态进行表征;该CT检测装置的成像方法包括多能量段成像、切片重建、感兴趣区域扫描、等效原子序数表征;利用光子计数探测器能量分辨成像的优势,结合图像映射模型,弥补CT检测装置对材料成分表征的不足,实现等效原子序数成像,通过融合成像方法,最终实现被测物的内部结构与成分自动CT检测。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要实用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明结构图。
1:射线源;2:射线源端Z轴;3:射线源端Y轴;4:射线源端X轴;5:底座;6:样品台;7:被测物端Z轴;8:被测物端Y轴;9:被测物端X轴;10:旋转轴;11:探测器端X轴;12:探测器端Y轴;13:探测器端Z轴;14:平板探测器;15:光子计数探测器。
图2为本发明的探测中心调节流程图;
图3为本发明的系统成像几何参数自动调节流程图;
图4为本发明的感兴趣区域探测图;
图5为本发明实施例获取的不同能量段CT切片;a)高能量段CT切片,b)低能量段CT切片,c)多能量段融合CT切片;
图6为本发明实施例获取的单晶件能谱相除切片;
1)高纯Ti、2)Al、3)LiF、4)SiO2、5)Si、6)Al2O3
图7为本发明的内部结构和成分融合成像示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将对本发明的技术方案进行详细的描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施方式,都属于本发明所保护的范围。
在任一实施例中,如图1所示,本发明的一种用于表征材料内部结构及成分的CT检测装置,包括:包括底座5,所述底座5中央设有旋转轴10,所述旋转轴10上设有被测物端X轴9,所述被测物端X轴9上设有被测物端Y轴8,所述被测物端Y轴8上设有被测物端Z轴7,所述被测物端Z轴7上设有样品台6;
所述旋转轴10一侧设有射线源端X轴4,所述射线源端X轴4上设有射线源端Y轴3,所述射线源端Y轴3上设有射线源端Z轴2,所述射线源端Z轴2上设有射线源1;
所述旋转轴10另一侧设有探测器端X轴11,所述探测器端X轴11上设有探测器端Y轴12,所述探测器端Y轴12上设有探测器端Z轴13和平板探测器14,所述探测器端Z轴13上设有光子计数探测器15;
所述被测物端X轴9、射线源端X轴4、探测器端X轴11在所述底座5上的投影在一条直线上。
在一个具体的实施例中,本发明的CT检测装置,可通过被测物端X轴和Y轴的调节,保证被测物探测中心位于旋转轴线,调节方法如图2,系统首先进行几何参数校准,确定固定扫描条件下的几何放大比λ;之后判断被测物是否为回转体,如果被测物为回转体,则旋转轴转至90°,标记被测物成像最右侧边界坐标(X1,Y1),再将旋转轴转至-90°,标记被测物成像最右侧边界坐标(X1,Y2),调节被测物X轴使得被测物成像最右侧边界位于坐标(X1,(Y2+Y1)/2),旋转轴转至0°,调节被测物Y轴使得被测物成像最右侧边界位于坐标(X1,(Y2+Y1)/2);如果被测物不为回转体,旋转轴转至90°,调节被测物端X轴使被测物成像中心(Xoc1,Yoc1)位于图像中心(Xc,Yoc1),被测物端X轴的移动距离为(Xc-Xoc1)*探测器像元尺寸/λ,之后旋转轴转至0°,调节被测物端Y轴使被测物成像中心(Xoc2,Yoc1)位于图像中心(Xc,Yoc1),被测物端Y轴的移动距离为(Xc-Xo2)c*探测器像元尺寸λ/;调节完成后,被测物探测中心已位于旋转轴线,最后被测物端的Z轴进行被测物高度调节,可将被测物投影移至探测器成像中心,便于被测物成像。
在一个具体的实施例中,本发明的CT检测装置,利用射线源端的Y轴与Z轴调节射线束的中心位置,利用探测器端的Y轴调节探测器端的中心位置;利用射线源端和探测器端的X轴进行成像几何放大比调节,检测装置可实现CT系统成像几何参数自动调节,自动调节流程如图3,同时可对被测物感兴趣区域进行聚焦成像,保证检测全过程自动化。
在任一实施例中,本发明提供一种用于表征材料内部结构及成分的CT成像方法,所述CT成像方法包括多能量段成像、切片重建、感兴趣区域扫描、等效原子序数表征,能够对被测物进行内部结构及成分表征,获得被测物内部结构和成分的融合成像结果。
在一个具体的实施例中,本发明的CT检测方法,其中所述多能量段成像依靠光子计数探测器实现,对被测物的衰减投影信号进行能量阈值分段成像,将不易被平板探测器区分的信息在多能量段结合运算(权重加、减法)中凸显,实现直观分析评判。如图5,分别为高、低两个能量段下的CT切片图像,显示窗口为0~0.1,高能段的金属伪影弱于低能段(1号物体区域),但低能段的对比度明显强于高能段(2号物体区域),通过结合运算,可获得高对比度、低伪影噪声的多能量段融合图像结果。
在一个具体的实施例中,本发明的CT检测方法,其中所述感兴趣区域扫描如图4,利用平板探测器对检测样品进行全局DR成像,通过DR图像感兴趣区域选取,结合成像坐标系与实际坐标系转化,利用机械运动实现感兴趣区域的光子计数探测成像。
所述成像坐标系与实际坐标系转化公式为:
(ΔXobj,ΔYobj)=μflat*(Xflat-roi-c-Xflat-c,Yflat-roi-c-Yflat-c)/λflat
(ΔXobj,ΔYobj)=μphoton*(Xphoton-roi-c-Xphoton-c,Yphoton-roi-c-Yphoton-c)/λphoton
式中(ΔXobj,ΔYobj)为实际被测物坐标系中感兴趣区域中心与被测物中心的偏移量,λflat和λphoton分别为平板探测器和光子计数探测器的成像几何放大比,μflat和μphoton分别为平板探测器和光子计数探测器的单像素实际尺寸,(Xflat-c,Yflat-c)为平板DR投影图像中物体的中心坐标,(Xflat-roi-c,Yflat-roi-c)为平板DR投影图像中选取的感兴趣区域中心坐标,同理,(Xphoton-c,Yphoton-c)和(Xphoton-roi-c,Yphoton-roi-c)分别为光子计数DR投影图像的中心坐标和感兴趣区域中心坐标。
在系统几何参数自动校准后,(Xflat-c,Yflat-c)与(Xphoton-c,Yphoton-c)的实际空间位置一致,依靠探测器端的X轴进行放大比调节,保证λflat=λphoton,光子计数探测平面位置利用探测器端的Y轴与Z轴调整,调整量为μflat*(Xflat-roi-c-Xflat-c,Yflat-roi-c-Yflat-c);若由于机械行程或探测条件导致λflat≠λphoton,则调整量为μflat*(Xflat-roi-c-Xflat-c,Yflat-roi-c-Yflat-c)*λphoton/λflat。
在一个具体的实施例中,本发明的CT检测方法,其中所述等效原子序数表征利用光子对检测样品进行多能量段投影采集,通过切片重建将多能量段投影转为多能量段切片,结合等效原子序数计算公式进行映射,获得检测样品等效原子序数切片图像。
所述等效原子序数计算公式为:
其中代表相除切片图像,代表低能段切片图像,为高能段切片图像,高、低能衰减系数可分解为光电效应与康普顿散射,fp(E)和fc(E)分别表示光电效应、康普顿散射与入射能量的关系,n取值3.94,在能量确定的情况下,可将相除图像和等效原子序数Z按三个参数的指数模型进行关系建立。
通过三种以上标准物进行模型参数确定,相除切片图像依据关系模型可映射到等效原子序数切片图,所述标准物采用纯度不低于99.99%的单晶成型件。为进一步提升材料识别成像量化分析精度,检测样品应与标准物同时检测,并且检测样品的等效原子序数应在标准物的原子序数极值范围内。图6为本发明实施例上获取的单晶标准在相除切片图像,表1为依据关系模型映射得到的等效原子序数结果,结果显示标准件的等效原子序数识别误差小于3%。
表1能谱CT处理结果表
根据本公开的一些实施例的CT成像方法,其中所述内部结构和成分的融合成像过程如图7,利用平板或光子计数探测器单能量段投影获取内部结构切片图像,依据等效原子序数计算公式获取内部等效原子序数切片图像,在内部结构切片图像上将等效原子序数切片图像以不同颜色进行绘制,图像中心像素点为对齐基准,获得内部结构和成分的融合成像结果。
对于λflat≠λphoton的情况,等效原子序数切片图像需要进行λphoton/λflat倍缩放,在平板探测器获得的内部结构切片图像上,以图像中心像素点为对齐基准,绘制缩放后的等效原子序数切片图像,舍弃超过内部结构切片图像矩阵索引范围的像素点。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合,为了避免不必要的重复,本发明对各种可能的组合方式不再另行说明。此外,本发明的各种不同的实施方式之间也可以进行任意组合,只要其不违背本发明的思想,其同样应当视为本发明所公开的内容。
Claims (7)
1.一种用于表征材料内部结构及成分的CT检测装置,其特征在于,包括底座,所述底座中央设有旋转轴,所述旋转轴上设有被测物端X轴,所述被测物端X轴上设有被测物端Y轴,所述被测物端Y轴上设有被测物端Z轴,所述被测物端Z轴上设有样品台;
所述旋转轴一侧设有射线源端X轴,所述射线源端X轴上设有射线源端Y轴,所述射线源端Y轴上设有射线源端Z轴,所述射线源端Z轴上设有射线源;
所述旋转轴另一侧设有探测器端X轴,所述探测器端X轴上设有探测器端Y轴,所述探测器端Y轴上设有探测器端Z轴和平板探测器,所述探测器端Z轴上设有光子计数探测器;
所述被测物端X轴、射线源端X轴、探测器端X轴在所述底座上的投影在一条直线上。
2.如权利要求1所述的CT检测装置,其特征在于,所述平板探测器用于检测样品内部结构成像,所述光子探测器用于检测样品内部成分成像。
3.如权利要求1-2任一项所述的用于表征材料内部结构及成分的CT检测装置的成像方法,其特征在于,包括多能量段成像、切片重建、感兴趣区域扫描、等效原子序数表征,对被测物进行内部结构及成分表征,获得被测物内部结构和成分的融合成像结果。
4.如权利要求3所述的成像方法,其特征在于,所述多能量段成像依靠光子计数探测器实现,对被测物的衰减投影信号进行能量阈值分段成像,将不易被平板探测器区分的信息在多能量段投影结合运算,实现直观分析评判。
5.如权利要求3所述的成像方法,其特征在于,所述感兴趣区域扫描利用平板探测器对检测样品进行全局DR成像,通过DR图像感兴趣区域选取,结合成像坐标系与实际坐标系转化,利用机械运动实现感兴趣区域的光子计数探测成像;
所述成像坐标系与实际坐标系转化公式为:
(ΔXobj,ΔYobj)=μflat*(Xflat-roi-c-Xflat-c,Yflat-roi-c-Yflat-c)/λflat
(ΔXobj,ΔYobj)=μphoton*(Xphoton-roi-c-Xphoton-c,Yphoton-roi-c-Yphoton-c)/λphoton
式中(ΔXobj,ΔYobj)为实际被测物坐标系中感兴趣区域中心与被测物中心的偏移量,λflat和λphoton分别为平板探测器和光子计数探测器的成像几何放大比,μflat和μphoton分别为平板探测器和光子计数探测器的单像素实际尺寸,(Xflat-c,Yflat-c)为平板DR投影图像中物体的中心坐标,(Xflat-roi-c,Yflat-roi-c)为平板DR投影图像中选取的感兴趣区域中心坐标,同理,(Xphoton-c,Yphoton-c)和(Xphoton-roi-c,Yphoton-roi-c)分别为光子计数DR投影图像的中心坐标和感兴趣区域中心坐标;
在系统几何参数自动校准后,(Xflat-c,Yflat-c)与(Xphoton-c,Yphoton-c)的实际空间位置一致,依靠探测器端的X轴进行放大比调节,保证λflat=λphoton,光子计数探测平面位置利用探测器端的Y轴与Z轴调整,调整量为μflat*(Xflat-roi-c-Xflat-c,Yflat-roi-c-Yflat-c);若由于机械行程或探测条件导致λflat≠λphoton,则调整量为μflat*(Xflat-roi-c-Xflat-c,Yflat-roi-c-Yflat-c)*λphoton/λflat。
7.如权利要求3所述的成像方法,其特征在于,将平板或光子计数探测器单能量段投影获取内部结构切片图像,并获取内部等效原子序数切片图像,在内部结构切片图像上将等效原子序数切片图像以不同颜色进行绘制,图像中心像素点为对齐基准,获得内部结构和成分的融合成像结果;
对于λflat≠λphoton的情况,等效原子序数切片图像需要进行λphoton/λflat倍缩放,在平板探测器获得的内部结构切片图像上,以图像中心像素点为对齐基准,绘制缩放后的等效原子序数切片图像,舍弃超过内部结构切片图像矩阵索引范围的像素点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210541988.7A CN114813798B (zh) | 2022-05-18 | 2022-05-18 | 用于表征材料内部结构及成分的ct检测装置和成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210541988.7A CN114813798B (zh) | 2022-05-18 | 2022-05-18 | 用于表征材料内部结构及成分的ct检测装置和成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114813798A true CN114813798A (zh) | 2022-07-29 |
CN114813798B CN114813798B (zh) | 2023-07-07 |
Family
ID=82514992
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210541988.7A Active CN114813798B (zh) | 2022-05-18 | 2022-05-18 | 用于表征材料内部结构及成分的ct检测装置和成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114813798B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117152077A (zh) * | 2023-08-30 | 2023-12-01 | 苏州谱睿源电子有限公司 | 基于层析融合技术的无损探测方法、装置及存储介质 |
Citations (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001198118A (ja) * | 2000-01-24 | 2001-07-24 | Hitachi Medical Corp | X線ct装置 |
US20060204972A1 (en) * | 2005-03-11 | 2006-09-14 | Roche Molecular Systems, Inc. | Systems and methods for determining real-time PCR cycle thresholds using a rotation transformation |
JP2007044275A (ja) * | 2005-08-10 | 2007-02-22 | Hitachi Medical Corp | マルチエナジーx線ct装置 |
CN101308102A (zh) * | 2008-07-16 | 2008-11-19 | 中北大学 | 一种计算机断层扫描成像装置和方法 |
US20110116594A1 (en) * | 2008-09-30 | 2011-05-19 | Hitachi Medical Corporation | X-ray ct apparatus |
JP2013116213A (ja) * | 2011-12-02 | 2013-06-13 | Hitachi Medical Corp | X線ct装置 |
CN103327901A (zh) * | 2011-09-28 | 2013-09-25 | 株式会社东芝 | X射线ct装置、物质确定方法及图像处理装置 |
CN103810712A (zh) * | 2014-03-05 | 2014-05-21 | 南方医科大学 | 一种能谱ct图像质量的评价方法 |
CN103913472A (zh) * | 2012-12-31 | 2014-07-09 | 同方威视技术股份有限公司 | Ct成像系统和方法 |
US20140362970A1 (en) * | 2013-06-10 | 2014-12-11 | General Electric Company | System and method of correcting banding artifacts in cardiac ct |
CN104422704A (zh) * | 2013-08-21 | 2015-03-18 | 同方威视技术股份有限公司 | 对x光能谱ct的能谱信息进行分解的方法和对应的重建方法 |
US20150085981A1 (en) * | 2013-09-24 | 2015-03-26 | Siemens Aktiengesellschaft | Method of image registration in a multi-source/single detector radiographic imaging system, and image acquisition apparatus |
CN104769422A (zh) * | 2012-09-07 | 2015-07-08 | 卡尔蔡司X射线显微镜公司 | 组合共焦x射线荧光和x射线计算机断层扫描的系统和方法 |
CN105092617A (zh) * | 2015-09-18 | 2015-11-25 | 重庆大学 | 一种基于x射线能谱ct和x射线荧光ct技术的双模态分子成像系统 |
CN106384342A (zh) * | 2015-07-24 | 2017-02-08 | 富士通株式会社 | 用于图像处理的计算机可读存储介质、设备和方法 |
US20170215969A1 (en) * | 2014-06-11 | 2017-08-03 | Tsinghua University | Human organ movement monitoring method, surgical navigation system and computer readable medium |
CN107610196A (zh) * | 2017-09-04 | 2018-01-19 | 沈阳东软医疗系统有限公司 | 双能量ct图像处理方法、装置以及图像处理设备 |
CN108169257A (zh) * | 2017-12-08 | 2018-06-15 | 中国科学院高能物理研究所 | 具有三种定量成像机制的纳米ct三维成像方法 |
CN108168996A (zh) * | 2017-12-13 | 2018-06-15 | 西安交通大学 | 一种基于声发射信号的pbx损伤演化ct原位分析方法 |
US20190162679A1 (en) * | 2017-05-16 | 2019-05-30 | Job Corporation, | Apparatus and method of processing data acquired in x-ray examination, and x-ray examination system equipped with the apparatus |
JP2019158447A (ja) * | 2018-03-09 | 2019-09-19 | 東芝Itコントロールシステム株式会社 | Ct撮影装置 |
CN110321943A (zh) * | 2019-06-25 | 2019-10-11 | 中国科学院自动化研究所 | 基于半监督深度学习的ct图像分类方法、系统、装置 |
CN110702706A (zh) * | 2019-09-20 | 2020-01-17 | 天津大学 | 一种能谱ct系统输出数据的模拟方法 |
US20200051229A1 (en) * | 2013-02-04 | 2020-02-13 | Ge Sensing & Inspection Technologies Gmbh | Method for the non-destructive testing of the volume of a test object and testing device configured for carrying out such a method |
EP3657528A1 (en) * | 2018-11-26 | 2020-05-27 | FEI Company | Method of imaging a sample using an electron microscope |
CN111388001A (zh) * | 2020-03-31 | 2020-07-10 | 王振 | 一种非旋转式计算机断层成像装置及方法 |
US20200302660A1 (en) * | 2019-03-21 | 2020-09-24 | Ziehm Imaging Gmbh | X-ray system for the iterative determination of an optimal coordinate transformation between overlapping volumes that have been reconstructed from volume data sets of discretely scanned object areas |
CN111879798A (zh) * | 2020-06-19 | 2020-11-03 | 中国人民解放军战略支援部队信息工程大学 | 基于采集序列细分的纳米ct投影位置漂移校正方法及装置 |
EP3739324A1 (en) * | 2019-05-16 | 2020-11-18 | Shimadzu Corporation | Reconstruction device |
CN112304987A (zh) * | 2020-10-19 | 2021-02-02 | 西北工业大学 | 基于光子计数能谱ct的含能材料等效原子序数测量方法 |
CN112365492A (zh) * | 2020-11-27 | 2021-02-12 | 上海联影医疗科技股份有限公司 | 图像扫描方法、装置、电子设备和存储介质 |
CN213239979U (zh) * | 2020-11-03 | 2021-05-18 | 扬州哈工博浩智能科技有限公司 | 一种工业ct检测设备 |
EP3831577A1 (fr) * | 2019-12-06 | 2021-06-09 | Commissariat à l'Energie Atomique et aux Energies Alternatives | Methode et systeme d'inspection en 3d d'une piece en cours de fabrication par un processus de type additif |
US20210170198A1 (en) * | 2019-12-10 | 2021-06-10 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for imaging of moving subjects |
US20210244374A1 (en) * | 2018-07-30 | 2021-08-12 | Xenselab Llc | Apparatus and methods for x-ray imaging |
CN114259298A (zh) * | 2022-02-25 | 2022-04-01 | 博志生物科技(深圳)有限公司 | 一种规划椎弓根后路内固定置钉轨迹的数字可视化方法 |
-
2022
- 2022-05-18 CN CN202210541988.7A patent/CN114813798B/zh active Active
Patent Citations (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001198118A (ja) * | 2000-01-24 | 2001-07-24 | Hitachi Medical Corp | X線ct装置 |
US20060204972A1 (en) * | 2005-03-11 | 2006-09-14 | Roche Molecular Systems, Inc. | Systems and methods for determining real-time PCR cycle thresholds using a rotation transformation |
JP2007044275A (ja) * | 2005-08-10 | 2007-02-22 | Hitachi Medical Corp | マルチエナジーx線ct装置 |
CN101308102A (zh) * | 2008-07-16 | 2008-11-19 | 中北大学 | 一种计算机断层扫描成像装置和方法 |
US20110116594A1 (en) * | 2008-09-30 | 2011-05-19 | Hitachi Medical Corporation | X-ray ct apparatus |
CN103327901A (zh) * | 2011-09-28 | 2013-09-25 | 株式会社东芝 | X射线ct装置、物质确定方法及图像处理装置 |
JP2013116213A (ja) * | 2011-12-02 | 2013-06-13 | Hitachi Medical Corp | X線ct装置 |
CN104769422A (zh) * | 2012-09-07 | 2015-07-08 | 卡尔蔡司X射线显微镜公司 | 组合共焦x射线荧光和x射线计算机断层扫描的系统和方法 |
CN103913472A (zh) * | 2012-12-31 | 2014-07-09 | 同方威视技术股份有限公司 | Ct成像系统和方法 |
US20200051229A1 (en) * | 2013-02-04 | 2020-02-13 | Ge Sensing & Inspection Technologies Gmbh | Method for the non-destructive testing of the volume of a test object and testing device configured for carrying out such a method |
US20140362970A1 (en) * | 2013-06-10 | 2014-12-11 | General Electric Company | System and method of correcting banding artifacts in cardiac ct |
CN104422704A (zh) * | 2013-08-21 | 2015-03-18 | 同方威视技术股份有限公司 | 对x光能谱ct的能谱信息进行分解的方法和对应的重建方法 |
US20150085981A1 (en) * | 2013-09-24 | 2015-03-26 | Siemens Aktiengesellschaft | Method of image registration in a multi-source/single detector radiographic imaging system, and image acquisition apparatus |
CN103810712A (zh) * | 2014-03-05 | 2014-05-21 | 南方医科大学 | 一种能谱ct图像质量的评价方法 |
US20170215969A1 (en) * | 2014-06-11 | 2017-08-03 | Tsinghua University | Human organ movement monitoring method, surgical navigation system and computer readable medium |
CN106384342A (zh) * | 2015-07-24 | 2017-02-08 | 富士通株式会社 | 用于图像处理的计算机可读存储介质、设备和方法 |
CN105092617A (zh) * | 2015-09-18 | 2015-11-25 | 重庆大学 | 一种基于x射线能谱ct和x射线荧光ct技术的双模态分子成像系统 |
US20190162679A1 (en) * | 2017-05-16 | 2019-05-30 | Job Corporation, | Apparatus and method of processing data acquired in x-ray examination, and x-ray examination system equipped with the apparatus |
CN107610196A (zh) * | 2017-09-04 | 2018-01-19 | 沈阳东软医疗系统有限公司 | 双能量ct图像处理方法、装置以及图像处理设备 |
CN108169257A (zh) * | 2017-12-08 | 2018-06-15 | 中国科学院高能物理研究所 | 具有三种定量成像机制的纳米ct三维成像方法 |
CN108168996A (zh) * | 2017-12-13 | 2018-06-15 | 西安交通大学 | 一种基于声发射信号的pbx损伤演化ct原位分析方法 |
JP2019158447A (ja) * | 2018-03-09 | 2019-09-19 | 東芝Itコントロールシステム株式会社 | Ct撮影装置 |
US20210244374A1 (en) * | 2018-07-30 | 2021-08-12 | Xenselab Llc | Apparatus and methods for x-ray imaging |
EP3657528A1 (en) * | 2018-11-26 | 2020-05-27 | FEI Company | Method of imaging a sample using an electron microscope |
US20200302660A1 (en) * | 2019-03-21 | 2020-09-24 | Ziehm Imaging Gmbh | X-ray system for the iterative determination of an optimal coordinate transformation between overlapping volumes that have been reconstructed from volume data sets of discretely scanned object areas |
EP3739324A1 (en) * | 2019-05-16 | 2020-11-18 | Shimadzu Corporation | Reconstruction device |
CN110321943A (zh) * | 2019-06-25 | 2019-10-11 | 中国科学院自动化研究所 | 基于半监督深度学习的ct图像分类方法、系统、装置 |
CN110702706A (zh) * | 2019-09-20 | 2020-01-17 | 天津大学 | 一种能谱ct系统输出数据的模拟方法 |
EP3831577A1 (fr) * | 2019-12-06 | 2021-06-09 | Commissariat à l'Energie Atomique et aux Energies Alternatives | Methode et systeme d'inspection en 3d d'une piece en cours de fabrication par un processus de type additif |
US20210170198A1 (en) * | 2019-12-10 | 2021-06-10 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for imaging of moving subjects |
CN111388001A (zh) * | 2020-03-31 | 2020-07-10 | 王振 | 一种非旋转式计算机断层成像装置及方法 |
CN111879798A (zh) * | 2020-06-19 | 2020-11-03 | 中国人民解放军战略支援部队信息工程大学 | 基于采集序列细分的纳米ct投影位置漂移校正方法及装置 |
CN112304987A (zh) * | 2020-10-19 | 2021-02-02 | 西北工业大学 | 基于光子计数能谱ct的含能材料等效原子序数测量方法 |
CN213239979U (zh) * | 2020-11-03 | 2021-05-18 | 扬州哈工博浩智能科技有限公司 | 一种工业ct检测设备 |
CN112365492A (zh) * | 2020-11-27 | 2021-02-12 | 上海联影医疗科技股份有限公司 | 图像扫描方法、装置、电子设备和存储介质 |
CN114259298A (zh) * | 2022-02-25 | 2022-04-01 | 博志生物科技(深圳)有限公司 | 一种规划椎弓根后路内固定置钉轨迹的数字可视化方法 |
Non-Patent Citations (8)
Title |
---|
TAO SHENGZHEN 等: "Multi-energy CT imaging for large patients using dual-source photon-counting detector CT", 《PHYSICS IN MEDICINE AND BIOLOGY》, vol. 65, no. 17, pages 1 - 10, XP020356666, DOI: 10.1088/1361-6560/ab99e4 * |
WANG MIAOSHI 等: "An adaptive reconstruction algorithm for spectral CT regularized by a reference image", 《PHYSICS IN MEDICINE AND BIOLOGY》, vol. 61, no. 24, pages 1 - 12 * |
ZENG-NIAN YUAN 等: "In-Situ X-ray Tomography Observation of Structure Evolution in 1,3,5-Triamino-2,4,6-Trinitrobenzene Based Polymer Bonded Explosive (TATB-PBX) under Thermo-Mechanical Loading", 《MATERIALS》, vol. 11, no. 5, pages 1 - 14 * |
曹波 等: "不同ROI的选择方法对胃癌能谱CT参数值影响的研究", 《中国医疗设备》, vol. 36, no. 10, pages 161 - 164 * |
杨亚飞 等: "基于光子计数能谱 CT 的含能材料等效原子序数测量方法", 《光谱学与光谱分析》, vol. 42, no. 5, pages 1400 - 1406 * |
盖立平 等: "能谱CT成像中测量位置与校正系数的关系", 《中国组织工程研究》, vol. 20, no. 31, pages 4677 - 4686 * |
钟晶 等: "双能量CT增强扫描在肺部血管容积再现中的价值", 《放射学实践》, vol. 31, no. 7, pages 701 - 704 * |
陈立新: "核素内照射三维剂量计算-<\'153>Sm治疗骨转移病例的应用", 《万方学术期刊数据库》, pages 1 - 97 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117152077A (zh) * | 2023-08-30 | 2023-12-01 | 苏州谱睿源电子有限公司 | 基于层析融合技术的无损探测方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN114813798B (zh) | 2023-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10521936B2 (en) | Device and method for image reconstruction at different X-ray energies, and device and method for X-ray three-dimensional measurement | |
JP2000321221A (ja) | Ctシステムで投影データを作成する方法およびシステム | |
JP2010032485A (ja) | X線分析装置及びx線分析方法 | |
CN114813798B (zh) | 用于表征材料内部结构及成分的ct检测装置和成像方法 | |
Villarraga-Gómez et al. | Effect of geometric magnification on dimensional measurements with a metrology-grade X-ray computed tomography system | |
CN108362716A (zh) | 一种基于医用ct的文物材质检测判定方法及检测装置 | |
CN104316546A (zh) | 一种复合材料天线罩质量无损跟踪与评价方法 | |
CN106934856A (zh) | 基于x射线检测技术的三维断层重建及切片显示方法 | |
CN111322970B (zh) | X射线胶路测量装置和方法 | |
CN109870471B (zh) | 一种单光栅侦测的锥束ct角度序列散射获取方法 | |
KR20200015375A (ko) | X선 투과 검사 장치 및 x선 투과 검사 방법 | |
JP2009175065A (ja) | 中性子即発ガンマ線分析による複数元素の同時3次元分布・可視化観察・計測方法及びその装置 | |
Xia et al. | Dual-energy micro-CT using GAGG: Ce and LYSO: Ce scintillators | |
CN105321206B (zh) | 一种适用于中子层析成像系统的旋转轴线偏摆角的误差补偿方法 | |
JP2008089325A (ja) | 電子プローブマイクロアナライザによる地質年代測定法 | |
Kato et al. | 3D spectroscopic computed tomography imaging using terahertz waves | |
Susanto et al. | Development of low-cost industrial x-ray computed tomography system based on digital fluoroscopy | |
JP2013079825A (ja) | X線ct画像再構成方法およびx線ct装置 | |
Seifert et al. | Resolving power of 3D x-ray microtomography systems | |
CN116413290B (zh) | 一种工业ct散射校正方法 | |
Scanavini et al. | Computed tomography | |
Jia et al. | The study of X-ray full-field imaging method for low atomic number elements based on Compton scattering | |
CN116211330A (zh) | 基于投影迹线的锥束cl系统几何参数敏感度分析方法及装置 | |
Uher et al. | RadalyX–Robotic X-ray scanner: applications in aerospace and material industry | |
CN116878430A (zh) | 一种基于x光射线微聚焦检测设备的珍珠层厚度测量方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |