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

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

Info

Publication number
CN112304987B
CN112304987B CN202011120012.XA CN202011120012A CN112304987B CN 112304987 B CN112304987 B CN 112304987B CN 202011120012 A CN202011120012 A CN 202011120012A CN 112304987 B CN112304987 B CN 112304987B
Authority
CN
China
Prior art keywords
atomic number
equivalent atomic
photon counting
energy spectrum
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.)
Active
Application number
CN202011120012.XA
Other languages
English (en)
Other versions
CN112304987A (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

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公式,其等效原子序数可表示为
其中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中,利用标定物的等效原子序数及其相对比值,拟合等效原子序数与相对比值的关系曲线为其中Zk,eff是步骤1中获取的标定物等效原子序数,R为相应标定物在SR(x,y)中的相对比值,a、b、c为关系曲线参数。
在上述步骤5中,根据步骤4的关系曲线R,可以得到等效原子序数图像计算公式再根据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为标定物序号,根据公式:
可计算得到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)中存在噪声,分别取各标定物的多个像素的灰度平均值。利用标定物的等效原子序数和相对比值拟合关系曲线其中R为相应标定物在SR(x,y)中的相对比值,a、b、c为关系曲线参数,拟合可得关系曲线参数为a=6.375,b=23390,c=125100。
步骤5:根据步骤4的关系曲线R,可以得到等效原子序数图像计算公式再根据SR(x,y)即可计算得到待测样品和标定物的等效原子序数图像Zeff(x,y),其中a、b、c的取值与步骤4中相同。
本发明提供的上述方法,所得到的等效原子序数图像计算公式,在同一光子计数能谱CT和相同扫描参数的情况下可以重复利用,即此时只需对待测样品(无标定物)执行步骤2、步骤3和步骤5,即可得到待测样品的等效原子序数图像。

Claims (2)

1.一种基于光子计数能谱CT的含能材料等效原子序数测量方法,其特征在于包括下述步骤:
步骤1:获取3个材质组成已知的标定物的等效原子序数Zk,eff,k为标定物序号,设第k个标定物由原子序数为Z1,Z2,…,Zn的原子组成,其等效原子序数表示为
其中ai为标定物中第i种原子的质量百分比,Ai为第i种原子的原子质量,ni为第i种原子的原子数目,Zi为第i种原子的原子序数;
步骤2:利用光子计数能谱CT,将待测样品和3个标定物同时进行CT扫描,并重建低能区CT图像SL(x,y)和高能区CT图像SH(x,y);
步骤3:将SL(x,y)除以SH(x,y),得到待测样品和标定物的相对比值图像SR(x,y);
步骤4:利用标定物的等效原子序数Zk,eff及其相对比值,拟合等效原子序数与相对比值的关系曲线R为相应标定物在SR(x,y)中的相对比值,a、b、c为关系曲线参数;
步骤5:利用SR(x,y)和步骤4的关系曲线R,计算得到待测样品和标定物的等效原子序数图像其中a、b、c的取值与步骤4中相同。
2.根据权利要求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 CN112304987A (zh) 2021-02-02
CN112304987B true 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)

Families Citing this family (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有效原子序数估计方法
CN114813798B (zh) * 2022-05-18 2023-07-07 中国工程物理研究院化工材料研究所 用于表征材料内部结构及成分的ct检测装置和成像方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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射线成像系统及方法
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成像中自适应能谱优化方法及其应用

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7190757B2 (en) * 2004-05-21 2007-03-13 Analogic Corporation Method of and system for computing effective atomic number images in multi-energy computed tomography
US10451568B2 (en) * 2014-08-22 2019-10-22 Canon Medical Systems Corporation Photon counting X-ray CT apparatus

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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射线成像系统及方法
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
"Precise material identification method based on a photon counting technique with correction of the beam hardening effect in X-ray spectra";N. Kimoto et al.;《Applied Radiation and Isotopes》(第124期);16-26 *
"基于光子计数能谱CT的含能材料等效原子序数测量方法";杨亚飞等;《光谱学与光谱分析》;第42卷(第5期);1400-1406 *
"多能谱X 射线成像技术及其在CT 中的应用";郝佳等;《CT 理论与应用研究》;第20卷(第1期);141-150 *

Also Published As

Publication number Publication date
CN112304987A (zh) 2021-02-02

Similar Documents

Publication Publication Date Title
US9687207B2 (en) Pre-reconstruction calibration, data correction, and material decomposition method and apparatus for photon-counting spectrally-resolving X-ray detectors and X-ray imaging
US9025815B2 (en) System and method for multi-material correction of image data
Beyer et al. The use of X-ray CT for attenuation correction of PET data
EP2377467A1 (en) System and method for determining the composition of an object
CN112304987B (zh) 基于光子计数能谱ct的含能材料等效原子序数测量方法
JP2016147055A (ja) X線ct装置及び補正方法
JP2016127907A (ja) 医用情報処理装置、x線ct装置、医用情報処理方法、及び医用情報処理プログラム
WO2009102990A2 (en) System and method for quantitative imaging of chemical composition to decompose multiple materials
CN110084864A (zh) 一种基于能谱ct的电子密度图像重建方法
CN108603849B (zh) 用于确定x射线管的有效能谱的装置
US11879856B2 (en) Two-step material decomposition calibration method for a full size photon counting computed tomography system
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
Chi et al. Experimental feasibility of dual-energy computed tomography based on the Thomson scattering X-ray source
JP2022113115A (ja) ビームハードニング補正方法、x線ct装置及びビームハードニング補正プログラム
Jumanazarov et al. Significance of the spectral correction of photon counting detector response in material classification from spectral x-ray CT
Baek et al. Diagnostic X-ray spectra detection by Monte Carlo simulation
Bisogni et al. Experimental study of beam hardening artifacts in photon counting breast computed tomography
CN106373100B (zh) 一种ct图像的硬化伪影校正方法
Rosso et al. Preliminary study of the advantages of X-ray energy selection in CT imaging
US11300696B2 (en) Radiation imaging device and photon counting type detector calibration method
EP4400872A1 (en) Method for calibrating detector and apparatus for calibrating detector
US20240230934A1 (en) Pcct energy calibration from x-ray tube spectra using a neural network
US12102471B2 (en) Iterative method for material decomposition calibration in a full size photon counting computed tomography system

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