CN104091329A - 一种ct图像的标定方法、装置和一种ct系统 - Google Patents
一种ct图像的标定方法、装置和一种ct系统 Download PDFInfo
- Publication number
- CN104091329A CN104091329A CN201410295062.XA CN201410295062A CN104091329A CN 104091329 A CN104091329 A CN 104091329A CN 201410295062 A CN201410295062 A CN 201410295062A CN 104091329 A CN104091329 A CN 104091329A
- Authority
- CN
- China
- Prior art keywords
- image
- calibration element
- fixed calibration
- theoretical value
- actual reconstruction
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000013507 mapping Methods 0.000 claims abstract description 60
- 238000002591 computed tomography Methods 0.000 claims abstract description 41
- 239000000463 material Substances 0.000 claims description 13
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 8
- 229910052782 aluminium Inorganic materials 0.000 claims description 8
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 229910001868 water Inorganic materials 0.000 claims description 8
- 239000004411 aluminium Substances 0.000 claims description 6
- 229910002804 graphite Inorganic materials 0.000 claims description 6
- 239000010439 graphite Substances 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 6
- 238000012423 maintenance Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000009977 dual effect Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000032683 aging Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000012491 analyte Substances 0.000 description 1
- 238000007630 basic procedure Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- 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
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
- G01V5/20—Detecting prohibited goods, e.g. weapons, explosives, hazardous substances, contraband or smuggled objects
- G01V5/22—Active interrogation, i.e. by irradiating objects or goods using external radiation sources, e.g. using gamma rays or cosmic rays
- G01V5/224—Multiple energy techniques using one type of radiation, e.g. X-rays of different energies
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
- G01V5/20—Detecting prohibited goods, e.g. weapons, explosives, hazardous substances, contraband or smuggled objects
- G01V5/22—Active interrogation, i.e. by irradiating objects or goods using external radiation sources, e.g. using gamma rays or cosmic rays
- G01V5/226—Active interrogation, i.e. by irradiating objects or goods using external radiation sources, e.g. using gamma rays or cosmic rays using tomography
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Pulmonology (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Spectroscopy & Molecular Physics (AREA)
Abstract
本发明提供一种CT图像的标定方法、装置和一种CT系统,包括:在CT扫描装置的通道区域外,且在其最大重建区域内设置固定式标定件,存储所述固定式标定件的理论值;采集所述固定式标定件的投影数据,得到所述固定式标定件的实际重建图像;将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数。本发明可以有效提高标定质量,提升图像标定效果,提高了CT扫描装置的可靠性,节约了维护成本,具有很高的实际应用价值。
Description
技术领域
本发明涉及辐射成像检测技术领域,尤其涉及一种CT图像的标定方法、装置和一种CT系统。
背景技术
CT技术由于能够消除物体重叠的影响,在安全检查等场合中发挥了重要作用。传统的CT采用滑环装置通过X光源和探测器的旋转来获取不同角度上的投影数据,通过重建方法来获取断层图像,从而获得被检测行李物品的内部信息。多能谱分析是指利用某种材料对不同能量的X射线吸收性能上的差异来分辨该材料。CT技术配合多能谱分析技术,使得目前的行李物品检查设备可以对被检物质内部任意位置的原子序数和电子密度进行重建,实现物质种类的识别,在爆炸物、毒品等检测中起到了较好的效果。
然而,不论是传统CT还是多能谱CT,由于光机和探测器的不一致性,同一型号的CT设备的重建图像在数值上总会存在差异,而随着工作时间的增加、光机、探测器的老化,即使是同一台CT设备,重建图像在数值上也会存在漂移,偏离出厂时的情况。为了解决上述问题,CT设备在出厂和定期维护时,都会进行标定工作,通过扫描特定材料的工件,将重建图像数值和工件真实值进行比较,确定后标定参数,实现CT图像的后标定。
但是,上述传统的后标定方法有着如下的不足:一方面,CT设备的老化是个持续的过程,而标定流程只能定期进行,如每年1~2次,这意味着后标定参数并不能在所有时间内都是最优的,并且由于标定流程技术性较强,必须有售后工程师来实现,维护成本增加,而且有可能占用设备正常工作时间;另一方面,受探测器、皮带、通道形状、重建算法等因素的影响,重建数值和理想值之间的差异在断层图像上并非各个区域都相同,而传统的标定工件通常放置于一个中等体积或者较小体积的工具箱中,放于CT传送带上送入扫描通道,实现标定,因此这种方法下往往只有皮带附近区域的数值后标定的较好,而其他区域并不能得到最佳的校正效果。。
发明内容
(一)要解决的技术问题
本发明提供一种CT图像的标定方法、装置和一种CT系统,以解决现有技术CT图像的后标定技术中标定不够及时,准确性不高的技术问题。
(二)技术方案
为解决上述技术问题,本发明提供一种CT图像的标定方法,包括:
在CT扫描装置的通道区域外,且在其最大重建区域内设置固定式标定件,存储所述固定式标定件的理论值;
采集所述固定式标定件的投影数据,得到所述固定式标定件的实际重建图像;
将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数。
进一步地,
所述固定式标定件为一组或多组,每组标定件为一个或多个。
进一步地,当所述固定式标定件为一组且为多个时,所述将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数还包括:
将组内多个固定式标定件互相作为参考,建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的映射函数。
进一步地,当所述固定式标定件为多组时,所述将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数还包括:
在同一组内建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的分组映射函数;
根据每一组的分组映射函数和各组之间的空间位置关系,使用线性和/或多项式插值手段,得到CT扫描装置的最大重建区域内不同空间位置的映射函数。
进一步地,
所述固定式标定件的实际重建图像包括:固定式标定件的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种重建图像;
相应地,所述固定式标定件的理论值包括:固定式标定件的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种理论值。
进一步地,所述方法还包括:
对CT扫描装置的通道区域内的扫描对象进行投影数据采集和图像重建,得到所述扫描对象的扫描重建图像;
使用所述映射函数对所述扫描重建图像进行校正,得到校正重建图像。
另一方面,本发明还提供一种CT图像的标定装置,包括:固定式标定件、存储单元、采集单元和映射单元,其中:
所述固定式标定件设置在CT扫描装置的通道区域外,最大重建区域内;
所述存储单元用于存储所述固定式标定件的理论值;
所述采集单元用于采集所述固定式标定件的投影数据,得到所述固定式标定件的实际重建图像;
所述映射单元用于将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数。
进一步地,
所述固定式标定件为一组或多组,每组标定件为一个或多个。
进一步地,所述映射单元还用于:
当所述固定式标定件为一组且为多个时,将组内多个固定式标定件互相作为参考,建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的映射函数。
进一步地,所述映射单元还用于:
当所述固定式标定件为多组时,在同一组内建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的分组映射函数;
根据每一组的分组映射函数和各组之间的空间位置关系,使用线性和/或多项式插值手段,得到CT扫描装置的最大重建区域内不同空间位置的映射函数。
进一步地,
所述固定式标定件的实际重建图像包括:固定式标定件的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种重建图像;
相应地,所述固定式标定件的理论值包括:固定式标定件的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种理论值。
进一步地,所述固定式标定件的材料为:
水、石墨、铝中的一种或多种。
进一步地,
所述固定式标定件的设置在CT扫描装置的通道侧面、通道顶部、皮带下方中的一个或多个位置。
进一步地,
所述采集单元还用于对CT扫描装置的通道区域内的扫描对象进行投影数据采集和图像重建,得到所述扫描对象的扫描重建图像;
所述装置还包括:校正单元,与所述采集单元和映射单元分别相连,用于使用所述映射函数对所述扫描重建图像进行校正,得到校正重建图像。
再一方面,本发明还提供一种CT系统,包括如上任一项所述的CT图像的标定装置。
(三)有益效果
可见,在本发明所提供的CT图像的标定方法、装置和一种CT系统中,能够利用固定式标定件对CT扫描装置进行实时标定,免去了CT扫描装置的传统标定方法中需要定期进行人工后标定的环节,可以有效提高标定质量,提升图像标定效果,提高了CT扫描装置的可靠性,节约了维护成本,具有很高的实际应用价值;并且,通过本发明的标定方法可以得到CT扫描装置最大重建区域内各个不同部分的校正值,而不是仅仅局限于通道区域附近,能够实现更为精细的图像校正。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例CT图像的标定方法基本流程示意图;
图2是本发明实施例CT图像的标定装置基本结构示意图;
图3是本发明实施例CT图像的标定装置的实际应用流程示意图;
图4是CT扫描装置的通道区域和最大重建区域示意图;
图5是本发明实施例CT图像的标定装置中固定式标定件放置位置示意图;
图6是本发明实施例CT图像的标定装置中线性插值的标定方法示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例首先提供一种CT图像的标定方法,参见图1,包括:
步骤101:在CT扫描装置的通道区域外,且在其最大重建区域内设置固定式标定件,存储所述固定式标定件的理论值;
步骤102:采集所述固定式标定件的投影数据,得到所述固定式标定件的实际重建图像;
步骤103:将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数。
可见,在本发明实施例所提供的CT图像的标定方法中,能够利用固定式标定件对CT扫描装置进行实时标定,免去了CT扫描装置的传统标定方法中需要定期进行人工后标定的环节,可以有效提高标定质量,提升图像标定效果,提高了CT扫描装置的可靠性,节约了维护成本,具有很高的实际应用价值。
优选地,固定式标定件可以为一组或多组,每组标定件也可以为一个或多个,取决于所需测量的重建图像类型。
优选地,当固定式标定件为一组且为多个时,可以对实际重建图像进行如下处理:将组内多个固定式标定件互相作为参考,建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的映射函数。
优选地,当固定式标定件为多组时,可以首先对其进行分组处理,然后再考虑各组之间的空间位置关系进行整合,具体地,可以包括:
在同一组内建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的分组映射函数;
根据每一组的分组映射函数和各组之间的空间位置关系,使用线性和/或多项式插值手段,得到CT扫描装置的最大重建区域内不同空间位置的映射函数。
本发明实施例的方法可以对多种CT重建图像进行标定,具体地,可以包括:固定式标定件的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布等各项数据,此处的特定元素可以为碳、水或铝等元素。
优选地,可以根据本发明实施例的标定方法对CT扫描装置进行扫描重建图像校正,具体地,方法还可以包括:
对CT扫描装置的通道区域内的扫描对象进行投影数据采集和图像重建,得到所述扫描对象的扫描重建图像;
使用所述映射函数对所述扫描重建图像进行校正,得到校正重建图像。
本发明实施例还提供一种CT图像的标定装置,参见图2,包括:
固定式标定件201,设置在CT扫描装置的通道区域1外,最大重建区域2内;
存储单元202,用于存储所述固定式标定件201的理论值;
采集单元203,用于采集所述固定式标定件201的投影数据,得到所述固定式标定件201的实际重建图像;
映射单元204,用于将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数。
图2中,3为射线源,4为传送带。
优选地,固定式标定件201可以为一组或多组,每组标定件也可以为一个或多个。
优选地,映射单元204还可以用于:当固定式标定件201为一组且为多个时,将组内多个固定式标定件201互相作为参考,建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的映射函数。
优选地,映射单元204还可以用于:
当固定式标定件201为多组时,在同一组内建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的分组映射函数;
根据每一组的分组映射函数和各组之间的空间位置关系,使用线性和/或多项式插值手段,得到CT扫描装置的最大重建区域2内不同空间位置的映射函数。
优选地,固定式标定件201的实际重建图像可以包括:固定式标定件201的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种重建图像;相应地,固定式标定件201的理论值也可以包括:固定式标定件201的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种理论值,此处的特定元素可以为碳、水或铝等元素。
优选地,固定式标定件201的材料可以为:水、石墨、铝等标定材料。
优选地,固定式标定件201的位置可以在CT扫描装置的通道侧面、通道顶部、皮带下方等,以得到CT扫描装置最大重建区域2内各个不同部分的校正值,而不是仅仅局限于通道区域1附近,能够实现更为精细的图像校正。
优选地,采集单元203还可以用于对CT扫描装置的通道区域1内的扫描对象进行投影数据采集和图像重建,得到所述扫描对象的扫描重建图像;相应地,装置还可以包括:校正单元,与采集单元203和映射单元204分别相连,用于使用映射函数对所述扫描重建图像进行校正,得到校正重建图像。
下面以在双能安检CT中对本发明实施例的CT图像的标定装置的实际应用流程为例,来详细说明本发明实施例的具体实现过程,参加图3:
步骤301:将固定式标定件置于双能安检CT的通道区域外,最大重建区域内。
通常扫描通道可以放置物体的有效区域都会小于扫描装置理论上最大可以重建的区域,如图4所示,因此通道外区域有一部分是可以进行重建的,但是在实际工作时通道外不可能有物体,因此这一部分一般情况下并没有重建出来。
本发明实施例中,即可利用上述区域,将固定式标定件置于双能安检CT的通道区域外,最大重建区域内,其中,本实施例中,共使用了6组固定式标定件,每组固定式标定件均包括纯水柱5、石墨柱6和铝柱7三种,见图5,并存储上述固定式标定件的各项理论值。在双能安检CT中,利用双能CT重建原理,可以重建高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数以及电子密度四种重建图像,由于原子序数以及电子密度重建图像对于材料的识别有着重要的作用,因此本实施例中,将选择原子序数和电子密度重建图像进行标定。上述三种材料的原子序数和电子密度的理论值均为已知,本实施例中,可以利用存储单元将水、石墨和铝材料的原子序数和电子密度理论值进行存储。
步骤302:采集6组固定式标定件的投影数据,得到固定式标定件的实际重建图像。
本实施例中,通过采集单元采集上述6组不同材料的固定式标定件的投影数据,得到其实际重建图像。
步骤303:将实际重建图像与所存储的相应的理论值比较,建立将实际重建图像校正为理论值的映射函数。
本步骤中,通过映射单元来将采集单元采集到的实际重建图像和存储单元所存储的水、石墨和铝材料的相应理论值分别进行比较,具体地:
可以首先在同一组内,将各固定式标定件互相作为参考,建立实际重建图像和相应的理论值之间的多项式映射关系,得到每一组的分组映射函数;然后通过线性插值方法由6个不同位置的分组映射函数建立整个重建范围内各位置的映射函数。
或者直接根据6组不同位置的固定式标定件的空间位置直接建立整个重建范围内各位置的映射函数。图6中给出了一种通过建立三角形的方法进行线性插值的实现手段,对于给定任一位置D,首先找到包含D点的最小三角形,由A、B、C三点构成,之后由A、B、C三点三种材料块的实际重建图像,估计D点的三种材料块的实际重建图像,这一过程可以通过双线性插值实现,之后根据推算出来的D点的三种材料的实际重建图像,得到D点的映射函数,分别实现D点的原子序数以及电子密度的校正。
步骤304:对CT扫描装置的通道区域内的扫描对象进行投影数据采集和图像重建,得到所述扫描对象的扫描重建图像;使用所述映射函数对所述扫描重建图像进行校正,得到校正重建图像。
本步骤中,首先通过采集单元对CT扫描装置的通道区域内的扫描对象进行投影数据采集和图像重建,得到扫描对象的实际的扫描重建图像,然后再根据之前建立的映射函数对扫描重建图像进行校正,将扫描重建图像映射回理论值,得到校正重建图像。
至此,则完成了本发明实施例CT图像的标定装置的实际应用流程的全过程。
本发明实施例装置中的存储单元、采集单元、映射单元和校正单元可以位于CT扫描装置的数据采集和图像重建装置内,也可以与CT扫描装置相连,可以成为一个整体,也可以由多个分体构成,只要具有上述功能单元即可。本领域技术人员可以采用硬件、软件、固件或三者的结合来实现上述功能单元。
本发明实施例还提供一种CT系统,包括如上任一项所述的CT图像的标定装置。
可见,本发明实施例至少具有如下有益效果:
在本发明实施例所提供的CT图像的标定方法、装置和一种CT系统中,能够利用固定式标定件对CT扫描装置进行实时标定,免去了CT扫描装置的传统标定方法中需要定期进行人工后标定的环节,可以有效提高标定质量,提升图像标定效果,提高了CT扫描装置的可靠性,节约了维护成本,具有很高的实际应用价值;并且,通过本发明实施例的标定方法可以得到CT扫描装置最大重建区域内各个不同部分的校正值,而不是仅仅局限于通道区域附近,能够实现更为精细的图像校正。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (15)
1.一种CT图像的标定方法,其特征在于,包括:
在CT扫描装置的通道区域外,且在其最大重建区域内设置固定式标定件,存储所述固定式标定件的理论值;
采集所述固定式标定件的投影数据,得到所述固定式标定件的实际重建图像;
将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数。
2.根据权利要求1所述的CT图像的标定方法,其特征在于:
所述固定式标定件为一组或多组,每组标定件为一个或多个。
3.根据权利要求2所述的CT图像的标定方法,其特征在于,当所述固定式标定件为一组且为多个时,所述将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数还包括:
将组内多个固定式标定件互相作为参考,建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的映射函数。
4.根据权利要求2所述的CT图像的标定方法,其特征在于,当所述固定式标定件为多组时,所述将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数还包括:
在同一组内建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的分组映射函数;
根据每一组的分组映射函数和各组之间的空间位置关系,使用线性和/或多项式插值手段,得到CT扫描装置的最大重建区域内不同空间位置的映射函数。
5.根据权利要求1所述的CT图像的标定方法,其特征在于:
所述固定式标定件的实际重建图像包括:固定式标定件的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种重建图像;
相应地,所述固定式标定件的理论值包括:固定式标定件的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种理论值。
6.根据权利要求1至5中任一项所述的CT图像的标定方法,其特征在于,所述方法还包括:
对CT扫描装置的通道区域内的扫描对象进行投影数据采集和图像重建,得到所述扫描对象的扫描重建图像;
使用所述映射函数对所述扫描重建图像进行校正,得到校正重建图像。
7.一种CT图像的标定装置,其特征在于,包括:固定式标定件、存储单元、采集单元和映射单元,其中:
所述固定式标定件设置在CT扫描装置的通道区域外,最大重建区域内;
所述存储单元用于存储所述固定式标定件的理论值;
所述采集单元用于采集所述固定式标定件的投影数据,得到所述固定式标定件的实际重建图像;
所述映射单元用于将所述实际重建图像与所存储的相应的理论值比较,建立将所述实际重建图像校正为理论值的映射函数。
8.根据权利要求7所述的CT图像的标定装置,其特征在于:
所述固定式标定件为一组或多组,每组标定件为一个或多个。
9.根据权利要求8所述的CT图像的标定装置,其特征在于,所述映射单元还用于:
当所述固定式标定件为一组且为多个时,将组内多个固定式标定件互相作为参考,建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的映射函数。
10.根据权利要求8所述的CT图像的标定装置,其特征在于,所述映射单元还用于:
当所述固定式标定件为多组时,在同一组内建立线性、多项式、指数中的一种或多种形式的实际重建图像和理论值之间的映射关系,根据所述映射关系建立将所述实际重建图像校正为理论值的分组映射函数;
根据每一组的分组映射函数和各组之间的空间位置关系,使用线性和/或多项式插值手段,得到CT扫描装置的最大重建区域内不同空间位置的映射函数。
11.根据权利要求7所述的CT图像的标定装置,其特征在于:
所述固定式标定件的实际重建图像包括:固定式标定件的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种重建图像;
相应地,所述固定式标定件的理论值包括:固定式标定件的高能衰减系数、低能衰减系数、某一特定X射线能量下的衰减系数、原子序数、电子密度、以特定元素为基材料的密度分布中的一种或多种理论值。
12.根据权利要求7所述的CT图像的标定装置,其特征在于,所述固定式标定件的材料为:
水、石墨、铝中的一种或多种。
13.根据权利要求7所述的CT图像的标定装置,其特征在于:
所述固定式标定件的设置在CT扫描装置的通道侧面、通道顶部、皮带下方中的一个或多个位置。
14.根据权利要求7所述的CT图像的标定装置,其特征在于:
所述采集单元还用于对CT扫描装置的通道区域内的扫描对象进行投影数据采集和图像重建,得到所述扫描对象的扫描重建图像;
所述装置还包括:校正单元,与所述采集单元和映射单元分别相连,用于使用所述映射函数对所述扫描重建图像进行校正,得到校正重建图像。
15.一种CT系统,其特征在于,包括如权利要求7至14中任一项所述的CT图像的标定装置。
Priority Applications (8)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410295062.XA CN104091329B (zh) | 2014-06-25 | 2014-06-25 | 一种ct图像的标定方法、装置和一种ct系统 |
PCT/CN2015/077098 WO2015196852A1 (zh) | 2014-06-25 | 2015-04-21 | 一种ct图像的标定方法、装置和一种ct系统 |
ES15173858.0T ES2652544T3 (es) | 2014-06-25 | 2015-06-25 | Método y dispositivo de calibración de imágenes por CT y sistema de CT |
EP15173858.0A EP2959836B1 (en) | 2014-06-25 | 2015-06-25 | Ct image calibration method and device and ct system |
US14/751,044 US9773326B2 (en) | 2014-06-25 | 2015-06-25 | Method and device for calibrating computerized tomography (CT) image and CT system including the same |
PL15173858T PL2959836T3 (pl) | 2014-06-25 | 2015-06-25 | Sposób kalibracji obrazu tomografii komputerowej oraz urządzenie i system tomografii komputerowej |
JP2015128026A JP6047633B2 (ja) | 2014-06-25 | 2015-06-25 | Ct画像のキャリブレーション方法、ct画像のキャリブレーション装置及びctシステム |
RU2015125284/28A RU2602750C1 (ru) | 2014-06-25 | 2015-06-25 | Способ калибровки компьютерно-томографического изображения, устройство и система компьютерной томографии |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410295062.XA CN104091329B (zh) | 2014-06-25 | 2014-06-25 | 一种ct图像的标定方法、装置和一种ct系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104091329A true CN104091329A (zh) | 2014-10-08 |
CN104091329B CN104091329B (zh) | 2017-02-15 |
Family
ID=51639044
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410295062.XA Active CN104091329B (zh) | 2014-06-25 | 2014-06-25 | 一种ct图像的标定方法、装置和一种ct系统 |
Country Status (8)
Country | Link |
---|---|
US (1) | US9773326B2 (zh) |
EP (1) | EP2959836B1 (zh) |
JP (1) | JP6047633B2 (zh) |
CN (1) | CN104091329B (zh) |
ES (1) | ES2652544T3 (zh) |
PL (1) | PL2959836T3 (zh) |
RU (1) | RU2602750C1 (zh) |
WO (1) | WO2015196852A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015196852A1 (zh) * | 2014-06-25 | 2015-12-30 | 清华大学 | 一种ct图像的标定方法、装置和一种ct系统 |
CN110444276A (zh) * | 2018-05-04 | 2019-11-12 | 西门子医疗有限公司 | 产生图像数据的方法、计算机断层扫描设备、程序产品和数据载体 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPWO2018062308A1 (ja) * | 2016-09-28 | 2019-07-11 | 国立大学法人群馬大学 | X線ct装置、電子密度及び実効原子番号の測定方法、ct検査方法、検査方法 |
US10621756B2 (en) * | 2017-01-13 | 2020-04-14 | Canon Medical Systems Corporation | Apparatus and method for correcting bias in low-count computed tomography projection data |
RU195421U1 (ru) * | 2019-10-17 | 2020-01-28 | Федеральное государственное унитарное предприятие "ВСЕРОССИЙСКИЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ ОПТИКО-ФИЗИЧЕСКИХ ИЗМЕРЕНИЙ" (ФГУП "ВНИИОФИ") | Калибровочное устройство для денситометрии |
CN112697821B (zh) * | 2020-12-02 | 2022-12-02 | 赛诺威盛科技(北京)股份有限公司 | 多能谱ct扫描方法、装置、电子设备和ct设备 |
CN113256629B (zh) * | 2021-07-05 | 2022-06-24 | 之江实验室 | 一种图像标定错误检测方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101266216A (zh) * | 2007-03-14 | 2008-09-17 | 清华大学 | 标定双能ct系统的方法和图像重建方法 |
CN101848678A (zh) * | 2007-11-06 | 2010-09-29 | 皇家飞利浦电子股份有限公司 | 用于ct体积中的新生脉管系统量化的系统 |
CN103325143A (zh) * | 2013-06-13 | 2013-09-25 | 华南理工大学 | 基于模型匹配的标记点自动注册方法 |
WO2013185011A1 (en) * | 2012-06-07 | 2013-12-12 | The Johns Hopkins University | Integration of quantitative calibration systems in computed tomography scanners |
Family Cites Families (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS54121089A (en) * | 1978-03-14 | 1979-09-19 | Toshiba Corp | Test device using radiant rays |
US4571491A (en) * | 1983-12-29 | 1986-02-18 | Shell Oil Company | Method of imaging the atomic number of a sample |
JPH08266530A (ja) * | 1995-03-31 | 1996-10-15 | Hitachi Medical Corp | X線ct装置 |
US6801594B1 (en) * | 1997-11-26 | 2004-10-05 | General Electric Company | Computed tomography fluoroscopy system |
DE69929026T2 (de) * | 1998-05-05 | 2006-08-24 | Koninklijke Philips Electronics N.V. | Bildgebungsverfahren für bildgesteuerte chirurgie |
DE10305105A1 (de) * | 2003-02-07 | 2004-08-26 | Siemens Ag | Eichung der Transformation spektraler Röntgenschwächungswerte in Dichte- und Ordnungszahlinformation |
DE10333101B4 (de) * | 2003-07-21 | 2008-05-21 | Qimonda Ag | Kalibrierungseinrichtung für die Kalibrierung eines Testerkanals einer Testereinrichtung, Testersystem und Verfahren zum Kalibrieren eines Testerkanals |
US7697738B2 (en) * | 2003-08-25 | 2010-04-13 | Koninklijke Philips Electronics N.V. | Calibration image alignment in a PET-CT system |
US7379527B2 (en) * | 2005-12-22 | 2008-05-27 | General Electric Company | Methods and apparatus for CT calibration |
US7459689B2 (en) * | 2006-06-20 | 2008-12-02 | Siemens Medical Solutions Usa, Inc. | Detector head position calibration and correction for SPECT imaging apparatus using virtual CT |
US7415093B2 (en) * | 2006-10-30 | 2008-08-19 | General Electric Company | Method and apparatus of CT cardiac diagnostic imaging using motion a priori information from 3D ultrasound and ECG gating |
US8011454B2 (en) * | 2007-09-25 | 2011-09-06 | Baker Hughes Incorporated | Apparatus and methods for continuous tomography of cores |
RU2359614C1 (ru) * | 2007-10-31 | 2009-06-27 | Закрытое Акционерное Общество "Импульс" | Способ калибровки цифрового рентгеновского аппарата (варианты) |
JP2009142397A (ja) * | 2007-12-12 | 2009-07-02 | Toshiba Corp | X線ct装置およびx線ct装置における検出器補正データの収集方法 |
JP5661624B2 (ja) * | 2008-08-13 | 2015-01-28 | コーニンクレッカ フィリップス エヌ ヴェ | 三次元回転型x線スキャナシステムの機械的アラインメントに起因するリング・アーチファクトの除去 |
DE102008052792B4 (de) * | 2008-10-22 | 2010-11-04 | Airbus Deutschland Gmbh | Luftführungselement und Luftführungselementanordnung für eine Flugzeugklimaanlage |
CN102335002B (zh) * | 2010-07-16 | 2015-01-14 | Ge医疗系统环球技术有限公司 | Ct机x光发生器及探测器位置校准工具,校准工具的校准方法和ct系统校准方法 |
US8712679B1 (en) * | 2010-10-29 | 2014-04-29 | Stc.Unm | System and methods for obstacle mapping and navigation |
JP2012189517A (ja) * | 2011-03-13 | 2012-10-04 | National Institute Of Advanced Industrial & Technology | X線ct装置の校正及び評価用の標準ゲージ、並びに該x線ct装置の校正及び評価用の標準ゲージを用いたx線ct装置の校正方法及び評価方法 |
JP5938860B2 (ja) * | 2011-08-29 | 2016-06-22 | 株式会社島津製作所 | 放射線撮影装置および断層画像補正方法 |
US9192303B2 (en) * | 2011-12-19 | 2015-11-24 | The Regents Of The University Of California | Temperature-modulated fluorescence tomography |
JP5953590B2 (ja) | 2012-03-21 | 2016-07-20 | 矢崎総業株式会社 | 圧着端子付き電線、および圧着端子を電線に圧着する方法 |
CN103226113B (zh) * | 2013-03-29 | 2015-10-21 | 中国计量科学研究院 | 锥束3d-ct扫描系统重建体素尺寸的自动标定方法 |
CN204008508U (zh) * | 2014-06-25 | 2014-12-10 | 清华大学 | 一种ct图像的标定装置和一种ct系统 |
CN104091329B (zh) * | 2014-06-25 | 2017-02-15 | 清华大学 | 一种ct图像的标定方法、装置和一种ct系统 |
-
2014
- 2014-06-25 CN CN201410295062.XA patent/CN104091329B/zh active Active
-
2015
- 2015-04-21 WO PCT/CN2015/077098 patent/WO2015196852A1/zh active Application Filing
- 2015-06-25 RU RU2015125284/28A patent/RU2602750C1/ru active
- 2015-06-25 PL PL15173858T patent/PL2959836T3/pl unknown
- 2015-06-25 JP JP2015128026A patent/JP6047633B2/ja active Active
- 2015-06-25 US US14/751,044 patent/US9773326B2/en active Active
- 2015-06-25 EP EP15173858.0A patent/EP2959836B1/en active Active
- 2015-06-25 ES ES15173858.0T patent/ES2652544T3/es active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101266216A (zh) * | 2007-03-14 | 2008-09-17 | 清华大学 | 标定双能ct系统的方法和图像重建方法 |
CN101848678A (zh) * | 2007-11-06 | 2010-09-29 | 皇家飞利浦电子股份有限公司 | 用于ct体积中的新生脉管系统量化的系统 |
WO2013185011A1 (en) * | 2012-06-07 | 2013-12-12 | The Johns Hopkins University | Integration of quantitative calibration systems in computed tomography scanners |
CN103325143A (zh) * | 2013-06-13 | 2013-09-25 | 华南理工大学 | 基于模型匹配的标记点自动注册方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015196852A1 (zh) * | 2014-06-25 | 2015-12-30 | 清华大学 | 一种ct图像的标定方法、装置和一种ct系统 |
US9773326B2 (en) | 2014-06-25 | 2017-09-26 | Tsinghua University | Method and device for calibrating computerized tomography (CT) image and CT system including the same |
CN110444276A (zh) * | 2018-05-04 | 2019-11-12 | 西门子医疗有限公司 | 产生图像数据的方法、计算机断层扫描设备、程序产品和数据载体 |
Also Published As
Publication number | Publication date |
---|---|
RU2602750C1 (ru) | 2016-11-20 |
CN104091329B (zh) | 2017-02-15 |
US20150379745A1 (en) | 2015-12-31 |
US9773326B2 (en) | 2017-09-26 |
JP2016007546A (ja) | 2016-01-18 |
PL2959836T3 (pl) | 2018-03-30 |
WO2015196852A1 (zh) | 2015-12-30 |
EP2959836B1 (en) | 2017-10-04 |
JP6047633B2 (ja) | 2016-12-21 |
ES2652544T3 (es) | 2018-02-05 |
EP2959836A1 (en) | 2015-12-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104091329A (zh) | 一种ct图像的标定方法、装置和一种ct系统 | |
CN103559699B (zh) | 一种基于投影估计的多能谱ct图像重建方法 | |
US8466420B2 (en) | Charge loss correction | |
CN103393434B (zh) | 获取正电子发射断层扫描系统响应模型与图像重建的方法 | |
CN204008508U (zh) | 一种ct图像的标定装置和一种ct系统 | |
CN102715914B (zh) | G型臂x光机的三维图像生成方法及装置与g型臂x光机 | |
EP3084720A1 (en) | Inspection system and method | |
CN105574845A (zh) | 一种多相机阵列烟标叠层数量测量方法及装置 | |
EP2961320B1 (en) | Gain calibration and correction in radiation system | |
CN103759679A (zh) | 一种锥束ct系统角度偏差测量方法 | |
US9763640B2 (en) | 3D image generation method and device for G-arm X-ray machine and G-arm X-ray machine | |
CN104596449B (zh) | 基于ct图像的小间隙精确测量方法 | |
CN104878679A (zh) | 一种车载车辙实时检测系统及检测方法 | |
CA2821773C (en) | Method for generating a 3d representation of an object | |
CN103961122A (zh) | 用于非等γ角CT系统中数据转换的方法和装置 | |
CN203720109U (zh) | 基于x射线源阵列的实时在线工业ct检测系统 | |
CN101603810B (zh) | 锥束ct系统快速定位方法及其专用测具 | |
CN103282938B (zh) | 用于通过x射线计算机断层扫描装置确定位于待检查物体中的结构位置的方法和评估设备 | |
US20140175298A1 (en) | 3d mapping with two orthogonal imaging views | |
CN105222730B (zh) | 一种基于图像复原的工业ct几何尺寸测量方法 | |
CN102590239B (zh) | 在液态物品检查系统中确定ct扫描位置的方法和设备 | |
Hung et al. | Development of correction schemes for a small field of view gamma camera | |
CN104306013B (zh) | 一种ct系统中线阵探测器位置的调整方法 | |
US9351697B2 (en) | Method and system for dimensional analysis of an object | |
CN102376085B (zh) | 用于乳房成像的核医学影像设备的图像衰减校正方法 |
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 |