CN112304987A - 基于光子计数能谱ct的含能材料等效原子序数测量方法 - Google Patents

基于光子计数能谱ct的含能材料等效原子序数测量方法 Download PDF

Info

Publication number
CN112304987A
CN112304987A CN202011120012.XA CN202011120012A CN112304987A CN 112304987 A CN112304987 A CN 112304987A CN 202011120012 A CN202011120012 A CN 202011120012A CN 112304987 A CN112304987 A CN 112304987A
Authority
CN
China
Prior art keywords
atomic number
equivalent atomic
photon counting
calibration object
sample
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
Application number
CN202011120012.XA
Other languages
English (en)
Other versions
CN112304987B (zh
Inventor
黄魁东
杨亚飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202011120012.XA priority Critical patent/CN112304987B/zh
Publication of CN112304987A publication Critical patent/CN112304987A/zh
Application granted granted Critical
Publication of CN112304987B publication Critical patent/CN112304987B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating 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/02Investigating 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/04Investigating 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/046Investigating 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]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/36Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/03Investigating materials by wave or particle radiation by transmission
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/10Different kinds of radiation or particles
    • G01N2223/101Different kinds of radiation or particles electromagnetic radiation
    • G01N2223/1016X-ray

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pulmonology (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明提供一种基于光子计数能谱CT的含能材料等效原子序数测量方法。该方法只需利用光子计数能谱CT对待测样品和3种已知材料的标定物进行CT扫描重建,进而由低、高能量区间CT图像的相对比值图像拟合出等效原子序数与相对比值的关系曲线,利用该关系曲线即可进一步计算得到待测样品和标定物的等效原子序数图像。该方法不依赖于光子计数能谱CT的专业知识,实施简便,具有很好的鲁棒性和通用性,可以大大降低现有方法的设备要求和算法复杂度。

Description

基于光子计数能谱CT的含能材料等效原子序数测量方法
技术领域
本发明涉及一种基于光子计数能谱CT的含能材料等效原子序数测量方法,属于工业CT无损检测技术领域。
背景技术
计算机断层成像(Computed Tomography,CT)在医学诊断及工业无损检测中有着广泛的应用。CT图像可看作是材料的线性衰减系数分布,该系数取决于X射线能量和材料特性(质量密度和等效原子序数)。在实际应用中,X射线管产生的光子能谱很宽(这与传统CT理论中X射线源只能发出特定能量光子的假设相违背),积分型探测器对整个能谱的光子进行积分,导致质量密度和等效原子序数均不相同的材料在给定的能谱下可能表现为相同的线性衰减系数,这最终导致传统CT无法有效完成材质的识别。以含能材料检测为例,在实际生产中可以经常发现含能材料样品中含有高线性衰减系数的杂质,其衰减值远高于正常含能材料,甚至高于金属铝。从CT理论分析,高线性衰减系数杂质可能是由于含能材料高聚集造成的(等效原子序数与其他接近,密度远高于正常含能材料),或是由于实际生产过程中混入了高原子序数的杂质。由于传统CT无法有效获得的质量密度和等效原子序数,因此很难判断该类杂质的组成并针对杂质成分改进生产工艺。
作为传统CT的补充,双能CT和能谱CT利用不同能谱的X射线对物体进行CT扫描。结合双能CT和多能CT分解算法,可以获得被扫描物体的有效原子序数和电子密度的分布。但实际应用中上述方法存在能谱测量困难,参考材料标定过程复杂,且可能放大噪声等缺点。据此基于同步辐射以氧气为标准元素的迭代方法被提出,但是该技术设备要求高,难以满足实际需要。后续基于光子计数探测器的标定方法相继被提出,但上述方法必须设置很窄的能量区间(0.5keV或1keV),严重降低了光子计数探测器有效能量区间内接收的光子数,导致CT扫描十分耗时。此外,待测材料的原子序数精度仍可能受重建误差、探测器响应准确性、散射效应等因素影响。
现有的技术缺点主要包括:
(1)双能CT和能谱CT存在能谱测量困难、算法复杂度高、参考材料标定过程复杂,且可能放大噪声等问题。
(2)基于迭代的同步辐射测量方法设备要求高,难以广泛应用。
(3)基于光子计数探测器的标定的方法扫描效率低,难以满足工业需要,且精度易受重建误差、探测器响应准确性、散射效应等因素影响。
综上所述,现有的技术方法存在算法复杂、设备要求高、扫描效率低等问题,难以满足实际含能材料生成中的高效检测需求。
发明内容
为提高实际含能材料生成中等效原子序数的识别精度,并降低设备要求和算法复杂度,本发明提供一种基于光子计数能谱CT的等效原子序数测量方法。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1:获取3个材质组成已知的标定物的等效原子序数Zk,eff,k为标定物序号;
步骤2:利用光子计数能谱CT,将待测样品和3个标定物同时进行CT扫描,并重建低能区CT图像SL(x,y)和高能区CT图像SH(x,y);
步骤3:将SL(x,y)除以SH(x,y),得到待测样品和标定物的相对比值图像SR(x,y);
步骤4:利用标定物的等效原子序数及其相对比值,拟合等效原子序数与相对比值的关系曲线;
步骤5:利用SR(x,y)和步骤4的关系曲线,计算得到待测样品和标定物的等效原子序数图像Zeff(x,y)。
在上述步骤1中,设第k个标定物由原子序数为Z1,Z2,…,Zn的原子组成,根据Glasser公式,其等效原子序数可表示为
Figure BDA0002731666340000021
其中ai为标定物中第i种原子的质量百分比,Ai为第i种原子的原子质量,ni为第i种原子的原子数目,Zi为第i种原子的原子序数。
在上述步骤3中,当实际图像中的噪声较大时,需要先对SL(x,y)和SH(x,y)进行降噪。为避免异常点,对于SH(x,y)中灰度值小于T的点不进行计算,直接将SR(x,y)中的对应点置为0,其中T为灰度阈值。
在上述步骤4中,利用标定物的等效原子序数及其相对比值,拟合等效原子序数与相对比值的关系曲线为
Figure BDA0002731666340000022
其中Zk,eff是步骤1中获取的标定物等效原子序数,R为相应标定物在SR(x,y)中的相对比值,a、b、c为关系曲线参数。
在上述步骤5中,根据步骤4的关系曲线R,可以得到等效原子序数图像计算公式
Figure BDA0002731666340000031
再根据SR(x,y)即可计算得到待测样品和标定物的等效原子序数图像Zeff(x,y),其中a、b、c的取值与步骤4中相同。
本发明提供的上述方法,所得到的等效原子序数图像计算公式,在同一光子计数能谱CT和相同扫描参数的情况下可以重复利用,即此时只需对待测样品(无标定物)执行步骤2、步骤3和步骤5,即可得到待测样品的等效原子序数图像。
本发明的有益效果是:该方法在实际应用中不依赖于光子计数能谱CT的专业知识,实施简便,具有很好的鲁棒性和通用性,可以大大降低现有方法的设备要求和算法复杂度,并且所得到的等效原子序数图像计算公式,在同一光子计数能谱CT和相同扫描参数的情况下可以重复利用。
下面结合附图和实施例对本发明做进一步说明。
附图说明
图1为本发明算法示意图。
具体实施方式
利用微焦点光子计数能谱CT(光子计数探测器为XCounter Hydra FX50)应用本发明,执行以下步骤:
步骤1:选择3个材质组成已知的标定物,依次为尼龙6、聚四氟乙烯、铝,获取其等效原子序数Zk,eff,k为标定物序号,根据公式:
Figure BDA0002731666340000032
可计算得到3种标定物的等效原子序数分别为6.6008、8.6163和13.0000。
步骤2:利用微焦点光子计数能谱CT,选择射线源管电压100kV和管电流70μA,同时扫描标定物和待测样品,扫描几何参数为:射线源到探测器距离420mm,射线源到旋转中心距离175mm,探测器像素尺寸为0.10mm,探测器采集速度为20帧/s,20帧图像求和得到一幅投影;低能量阈值为40keV,高能量阈值为60keV;单幅投影大小为1024×64,投影幅数为900幅,重建矩阵大小为1024×1024×64。利用FDK算法分别重建低能区CT图像SL(x,y)和高能区CT图像SH(x,y)。
步骤3:将SL(x,y)除以SH(x,y),得到待测样品和标定物的相对比值图像SR(x,y)。考虑到实际图像中噪声较大,先对SL(x,y)和SH(x,y)利用BM3D算法进行降噪。为避免异常点,对于SH(x,y)中灰度值小于T的点不进行计算,直接将SR(x,y)中的对应点置为0,其中T为灰度阈值,本实例中T取0.01。
步骤4:根据SR(x,y)获得尼龙6、聚四氟乙烯、铝3种标定物的相对比值分别为1.0912、1.1684和1.4854。获取3种标定物的相对比值时,考虑到SR(x,y)中存在噪声,分别取各标定物的多个像素的灰度平均值。利用标定物的等效原子序数和相对比值拟合关系曲线
Figure BDA0002731666340000041
其中R为相应标定物在SR(x,y)中的相对比值,a、b、c为关系曲线参数,拟合可得关系曲线参数为a=6.375,b=23390,c=125100。
步骤5:根据步骤4的关系曲线R,可以得到等效原子序数图像计算公式
Figure BDA0002731666340000042
再根据SR(x,y)即可计算得到待测样品和标定物的等效原子序数图像Zeff(x,y),其中a、b、c的取值与步骤4中相同。
本发明提供的上述方法,所得到的等效原子序数图像计算公式,在同一光子计数能谱CT和相同扫描参数的情况下可以重复利用,即此时只需对待测样品(无标定物)执行步骤2、步骤3和步骤5,即可得到待测样品的等效原子序数图像。

Claims (4)

1.一种基于光子计数能谱CT的含能材料等效原子序数测量方法,其特征在于包括下述步骤:
步骤1:获取3个材质组成已知的标定物的等效原子序数Zk,eff,k为标定物序号;
步骤2:利用光子计数能谱CT,将待测样品和3个标定物同时进行CT扫描,并重建低能区CT图像SL(x,y)和高能区CT图像SH(x,y);
步骤3:将SL(x,y)除以SH(x,y),得到待测样品和标定物的相对比值图像SR(x,y);
步骤4:利用标定物的等效原子序数及其相对比值,拟合等效原子序数与相对比值的关系曲线;
步骤5:利用SR(x,y)和步骤4的关系曲线,计算得到待测样品和标定物的等效原子序数图像Zeff(x,y)。
2.根据权利要求1所述的一种基于光子计数能谱CT的含能材料等效原子序数测量方法,其特征在于:在所述步骤4中,利用标定物的等效原子序数及其相对比值,拟合等效原子序数与相对比值的关系曲线为
Figure FDA0002731666330000011
R为相应标定物在SR(x,y)中的相对比值,a、b、c为关系曲线参数。
3.根据权利要求1所述的一种基于光子计数能谱CT的含能材料等效原子序数测量方法,其特征在于:在所述步骤5中,根据步骤4的关系曲线R,可以得到等效原子序数图像计算公式
Figure FDA0002731666330000012
再根据SR(x,y)即可计算得到待测样品和标定物的等效原子序数图像Zeff(x,y),其中a、b、c的取值与步骤4中相同。
4.根据权利要求1所述的一种基于光子计数能谱CT的含能材料等效原子序数测量方法,其特征在于:该方法所得到的等效原子序数图像计算公式,在同一光子计数能谱CT和相同扫描参数的情况下可以重复利用,即此时只需对待测样品(无标定物)执行步骤2、步骤3和步骤5,即可得到待测样品的等效原子序数图像。
CN202011120012.XA 2020-10-19 2020-10-19 基于光子计数能谱ct的含能材料等效原子序数测量方法 Active CN112304987B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011120012.XA CN112304987B (zh) 2020-10-19 2020-10-19 基于光子计数能谱ct的含能材料等效原子序数测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011120012.XA CN112304987B (zh) 2020-10-19 2020-10-19 基于光子计数能谱ct的含能材料等效原子序数测量方法

Publications (2)

Publication Number Publication Date
CN112304987A true CN112304987A (zh) 2021-02-02
CN112304987B CN112304987B (zh) 2024-04-12

Family

ID=74328283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011120012.XA Active CN112304987B (zh) 2020-10-19 2020-10-19 基于光子计数能谱ct的含能材料等效原子序数测量方法

Country Status (1)

Country Link
CN (1) CN112304987B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114403912A (zh) * 2022-01-17 2022-04-29 明峰医疗***股份有限公司 一种基于深度学习的能谱ct有效原子序数估计方法
CN114813798A (zh) * 2022-05-18 2022-07-29 中国工程物理研究院化工材料研究所 用于表征材料内部结构及成分的ct检测装置和成像方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050259781A1 (en) * 2004-05-21 2005-11-24 Zhengrong Ying Method of and system for computing effective atomic number images in multi-energy computed tomography
AU2014268157A1 (en) * 2011-07-29 2014-12-11 The Johns Hopkins University Differential phase contrast X-ray imaging and X-ray Illumination Methods
CN104346820A (zh) * 2013-07-26 2015-02-11 清华大学 一种x光双能ct重建方法
CN104970815A (zh) * 2014-04-04 2015-10-14 曹红光 基于光栅相位衬度和光子计数的x射线成像***及方法
US20160054453A1 (en) * 2014-08-22 2016-02-25 Kabushiki Kaisha Toshiba Photon counting x-ray ct apparatus
CN105559813A (zh) * 2014-11-04 2016-05-11 株式会社东芝 医用图像诊断装置以及医用图像处理装置
CN106353828A (zh) * 2015-07-22 2017-01-25 清华大学 在安检***中估算被检查物体重量的方法和装置
CN107356615A (zh) * 2016-05-10 2017-11-17 清华大学 一种用于双能x射线ct的方法和***
CN111707688A (zh) * 2020-06-18 2020-09-25 清华大学 光子计数能谱ct成像中自适应能谱优化方法及其应用

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050259781A1 (en) * 2004-05-21 2005-11-24 Zhengrong Ying Method of and system for computing effective atomic number images in multi-energy computed tomography
AU2014268157A1 (en) * 2011-07-29 2014-12-11 The Johns Hopkins University Differential phase contrast X-ray imaging and X-ray Illumination Methods
CN104346820A (zh) * 2013-07-26 2015-02-11 清华大学 一种x光双能ct重建方法
CN104970815A (zh) * 2014-04-04 2015-10-14 曹红光 基于光栅相位衬度和光子计数的x射线成像***及方法
US20160054453A1 (en) * 2014-08-22 2016-02-25 Kabushiki Kaisha Toshiba Photon counting x-ray ct apparatus
CN105559813A (zh) * 2014-11-04 2016-05-11 株式会社东芝 医用图像诊断装置以及医用图像处理装置
CN106353828A (zh) * 2015-07-22 2017-01-25 清华大学 在安检***中估算被检查物体重量的方法和装置
CN107356615A (zh) * 2016-05-10 2017-11-17 清华大学 一种用于双能x射线ct的方法和***
CN111707688A (zh) * 2020-06-18 2020-09-25 清华大学 光子计数能谱ct成像中自适应能谱优化方法及其应用

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
N. KIMOTO ET AL.: ""Precise material identification method based on a photon counting technique with correction of the beam hardening effect in X-ray spectra"", 《APPLIED RADIATION AND ISOTOPES》, no. 124, pages 16 - 26, XP029979981, DOI: 10.1016/j.apradiso.2017.01.049 *
杨亚飞等: ""基于光子计数能谱CT的含能材料等效原子序数测量方法"", 《光谱学与光谱分析》, vol. 42, no. 5, pages 1400 - 1406 *
郝佳等: ""多能谱X 射线成像技术及其在CT 中的应用"", 《CT 理论与应用研究》, vol. 20, no. 1, pages 141 - 150 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114403912A (zh) * 2022-01-17 2022-04-29 明峰医疗***股份有限公司 一种基于深度学习的能谱ct有效原子序数估计方法
CN114813798A (zh) * 2022-05-18 2022-07-29 中国工程物理研究院化工材料研究所 用于表征材料内部结构及成分的ct检测装置和成像方法

Also Published As

Publication number Publication date
CN112304987B (zh) 2024-04-12

Similar Documents

Publication Publication Date Title
GB2551029B (en) Method for performing material decoposition using a dual-energy X-Ray CT and corresponding Dual-energy X-Ray CT Apparatus
Le et al. Least squares parameter estimation methods for material decomposition with energy discriminating detectors
EP2377467A1 (en) System and method for determining the composition of an object
US7298812B2 (en) Image-based material decomposition
US7190757B2 (en) Method of and system for computing effective atomic number images in multi-energy computed tomography
US7197172B1 (en) Decomposition of multi-energy scan projections using multi-step fitting
Zhao et al. Dual-dictionary learning-based iterative image reconstruction for spectral computed tomography application
US10507005B2 (en) Spectral calibration of spectral computed tomography (CT)
CN112304987B (zh) 基于光子计数能谱ct的含能材料等效原子序数测量方法
Brambilla et al. Basis material decomposition method for material discrimination with a new spectrometric X-ray imaging detector
US20220028127A1 (en) Energy weighting of photon counts for conventional imaging
Clements et al. Multi-contrast CT imaging using a high energy resolution CdTe detector and a CZT photon-counting detector
US7050530B2 (en) Method and apparatus for the spatially-resolved determination of the element concentrations in objects to be examined
US10799192B2 (en) Method and apparatus for partial volume identification from photon-counting macro-pixel measurements
Jumanazarov et al. Significance of the spectral correction of photon counting detector response in material classification from spectral x-ray CT
Chi et al. Experimental feasibility of dual-energy computed tomography based on the Thomson scattering X-ray source
US11169286B2 (en) Methods of calibrating semiconductor radiation detectors using K-edge filters
Fu et al. Methods determining the angular increment of a continuous scan cone-beam CT system
CN116183647A (zh) 一种物质识别方法
Ogawa et al. Identification of a material with a photon counting x-ray CT system
Kumar et al. Behaviour of reconstructed attenuation values with X-ray tube voltage in an experimental third-generation industrial CT system using Xscan linear detector array
Haase et al. Estimation of statistical weights for model-based iterative CT reconstruction
CN111973205B (zh) 一种多层x射线探测器图像处理方法及终端
Martín et al. Processing of the signal obtained from X-ray fluorescence emissions to estimate localization in biomedical applications
Honkanen et al. Monochromatic computed tomography using laboratory-scale setup: proof-of-concept

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