CN103913779A - 多能ct成像系统以及成像方法 - Google Patents
多能ct成像系统以及成像方法 Download PDFInfo
- Publication number
- CN103913779A CN103913779A CN201310008543.3A CN201310008543A CN103913779A CN 103913779 A CN103913779 A CN 103913779A CN 201310008543 A CN201310008543 A CN 201310008543A CN 103913779 A CN103913779 A CN 103913779A
- Authority
- CN
- China
- Prior art keywords
- data
- energy
- ray
- angular interval
- imaging system
- 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
- 238000013170 computed tomography imaging Methods 0.000 title claims abstract description 44
- 238000003384 imaging method Methods 0.000 title abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 52
- 238000012545 processing Methods 0.000 claims abstract description 24
- 238000001514 detection method Methods 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims abstract description 11
- 238000009826 distribution Methods 0.000 claims description 21
- 238000010494 dissociation reaction Methods 0.000 claims description 20
- 230000005593 dissociations Effects 0.000 claims description 20
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 238000000354 decomposition reaction Methods 0.000 claims description 11
- 230000008034 disappearance Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000009977 dual effect Effects 0.000 description 28
- 230000000694 effects Effects 0.000 description 13
- 230000033001 locomotion Effects 0.000 description 13
- 239000000463 material Substances 0.000 description 12
- 238000010586 diagram Methods 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000002203 pretreatment Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 210000004872 soft tissue Anatomy 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 239000004411 aluminium Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003705 background correction Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002939 conjugate gradient method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000012631 diagnostic technique Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 108010048734 sclerotin Proteins 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- 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/54—Control of apparatus or devices for radiation diagnosis
-
- 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/226—Active interrogation, i.e. by irradiating objects or goods using external radiation sources, e.g. using gamma rays or cosmic rays using tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4869—Determining body composition
-
- 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/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- 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
-
- 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/40—Imaging
- G01N2223/419—Imaging computed tomograph
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Optics & Photonics (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Theoretical Computer Science (AREA)
- Immunology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及多能CT成像系统及成像方法。多能CT成像系统具有:载物台,载置受检查物品;电压可调的X射线发生装置,发出X射线;探测器,接收从X射线发生装置发出且穿过受检查物品的X射线,输出探测数据;机架,装配X射线发生装置及探测器;数据处理及控制装置,对载物台、X射线发生装置、探测器及机架控制进行,处理探测数据,在X射线发生装置的一周扫描过程中,根据预定的能量数将X射线发生装置的圆轨道平均分成预定个数的角度区间,对每个角度区间设定不同的X射线发生装置的高压,当机架从当前的角度区间旋转进入下一个角度区间时,对X射线发生装置进行控制,将X射线发生装置切换为在下一个角度区间设定的电压。
Description
技术领域
本发明涉及断层影像(CT)成像系统以及成像方法,特别涉及一种多能CT成像系统以及成像方法。
背景技术
断层影像(CT)的图像对比度与扫描所用的X光源能谱分布有很大关系。传统CT使用具有能谱分布的一个光源进行成像,有时会出现信息模糊致使两种不同材料在CT图像上完全相同,相对于此,双能CT使用两个不同分布的能谱对物体进行成像,能够消除单能谱情况下的信息模糊。双能X光断层成像技术(CT)利用物质在不同能量下的衰减差异能够获得关于物体的多重物理特性参数的分布信息,例如,电子密度分布、等效原子序数分布、多个能量下的单能衰减图像。这样,双能X光CT能够被用来进行传统CT的射线硬化校正、临床高对比度能谱图像获取、工业和安检的特殊和危险物品检测等。与传统的X光CT成像技术相比,双能CT在其成像功能上的突破对医学诊断技术、无损检测以及安全检查等应用都具有重大意义,因此,在近年来得到越来越广泛的关注。
此外,相对于传统的单能CT成像技术,双能CT成像技术不仅能够提供被检测物体的衰减系数和几何结构信息,还能够提供物质组成信息。因此,在安全检查领域,双能CT技术能够提供被检测物的电子密度和有效原子序数信息,从而进行危险物质的识别,另外,在医疗领域,双能CT能够提供不同组织成分的图像,例如骨质、软组织、对比度造影剂等。
目前,双能CT系统主要有双源双探测器、单源双层探测器、快速能量切换这三种实现方式。关于双源双探测器,如图1所示,这种系统由两套X光源以及探测器组成,即,这种系统包括高能X光源11、高能探测器12、以及低能X光源21、低能探测器22,并且,高能X光源11以及高能探测器12与低能X光源21以及低能探测器22分别呈90°交错放置。在数据采集时,这两套X光机发出不同能量(KVp)的射线,对应的探测器独立采集数据,由此,能够获得高能投影数据和低能投影数据这两组投影数据。但是,双源双探测器系统的成本十分昂贵,并且,对于旋转机架的机械结构的稳定性和强度的设计要求非常高。此外,关于单源双层探测器,如图2所示,在这种系统设计中,在传统的单能CT的探测器(图2中的低能探测器)之后增加了一个低能滤波片和一个探测器(图2中的高能探测器),由此,组成双能探测器。在X射线穿过第一个探测器(低能探测器)和滤波片后,低能部分的X射线被过滤除去,高能部分的X射线到达第二层探测器(高能探测器)。这两个探测器同时进行工作,从而能够采集到低能投影数据以及高能投影数据这两组投影数据。但是,这样的单源双层探测器的成本也较高。
此外,对于快速能量切换的方式来说,这种系统需要选用特殊的X光机,能够快速切换高压,交替发射不同能量(KVp)的射线。图3是示出用于实现双能CT系统的能量快速切换双能CT的示意图。如图3(B)所示,通过快速切换X光机高压值的方式能够实现高低能数据的采集方式。这种CT系统在扫描时机架正常旋转,例如转速0.5秒/圈,一圈进行1000次采样,则X光机的高压在每次采样中会切换一次,探测器则会读取两次数据,第一次为高能的投影数据,第二次为低能的投影数据。此时,X光机高压在一圈即0.5秒内切换1000次。但是,在快速能量切换系统中,需要搭配新型X光机,这导致此种方式成本变高并且也不容易普及应用。
如上所述,由于以上三种双能CT的制造成本都远大于传统的单能CT,所以,难以在一般的检测需求中推广应用。此外,由于双能CT不能够准确地反映X射线与物质的真实作用过程,对某些具有特征吸收的物质的重建结果的准确性较差,相对于此,多能CT则有望部分解决这一问题。因此,多能CT成像系统受到广泛关注。
发明内容
本发明是为了解决上述课题而提出的,其目的在于提供一种多能CT成像系统以及成像方法,能够在一周扫描过程中分多次改变X光机高压,从而能够以多能进行成像。
此外,本发明是在传统的单能CT成像系统的基础上得到的,能够在不增加硬件成本的情况下解决上述课题。
此外,本发明在硬件上利用与传统的单能CT基本相同的系统设计以及相应的数据处理和图像重建方法,提供多能成像功能,因此成本低,能够广泛应用于安全检查、医疗诊断等领域。
本发明提供一种多能CT成像系统,其特征在于,具有:
载物台,用于载置受检查物品;
电压可调的X射线发生装置,用于发出照射所述受检查物品的X射线;
探测器,用于接收从所述X射线发生装置发出且穿过所述受检查物品的X射线,并且,输出探测数据;
机架,装配有所述X射线发生装置以及所述探测器;以及
数据处理及控制装置,对所述载物台、所述X射线发生装置、所述探测器、所述机架进行控制,处理来自所述探测器的所述探测数据,在所述X射线发生装置的一周扫描过程中,所述数据处理及控制装置根据预定的能量数将所述X射线发生装置的圆轨道平均分成预定个数的角度区间,对每一个角度区间设定不同的所述X射线发生装置的高压,当所述机架从当前的角度区间旋转进入下一个角度区间时,根据来自所述数据处理及控制装置的指令,对所述X射线发生装置进行控制,将所述X射线发生装置切换为在下一个角度区间设定的电压。
此外,本发明提供一种多能CT成像系统的成像方法,其特征在于,具有如下步骤:
(a)按照预定的能量数将所述多能CT成像系统的X射线发生装置的圆轨道平均分成预定个数的角度区间,对每个角度区间设定不同的所述X射线发生装置的高压;
(b)所述X射线发生装置沿着圆轨道进行旋转,在每个角度区间以在该角度区间所设定的高压进行工作,发出X射线,照射受检查物品;
(c)利用所述多能CT成像系统的探测器采集不同能量下的投影数据;
(d)利用所述多能CT成像系统的数据处理及控制装置得到所述不同能量下的投影数据的线积分数据;
(e)将所述不同能量下的线积分数据依次进行拼接,由此得到覆盖360°的完整的线积分数据;
(f)根据所述完整的线积分数据进行重建,由此得到混合重建结果;
(g)将所述混合重建结果作为先验图像,对各段不完整数据分别进行重建,得到不同能量下的衰减系数重建结果;
(h)根据所述衰减系数重建结果得到分解系数,并且,利用所述分解系数得到原子序数分布图像和电子密度分布图像。
在本发明的多能CT成像系统的成像方法中,在所述步骤(a)中将圆轨道分为N角度区间并且使在所述步骤(d)中得到的不同能量下的线积分数据为P1,P2,…,Pn,并且,在所述步骤(f)中利用普通圆轨道滤波反投影法进行重建,使得到的混合重建结果为f0,N和n是大于等于2的整数,在所述步骤(g)中,通过求解有约束条件的最小化优化问题 得到满足约束条件的衰减系数重建结果f i ,其中,G(f i ,f0)为描述先验图像与重建结果的差异的函数,H i 为第i个能量下的投影算子,ε为与测量数据噪声方差有关的参数,i=1,2,…,n。
在本发明的多能CT成像系统的成像方法中,在所述步骤(h)中,使用后处理基效应分解法得到分解系数,并且,利用所述分解系数得到原子序数分布图像和电子密度分布图像。
在本发明的多能CT成像系统的成像方法中,在所述步骤(h)中,使用后处理基材料分解法得到分解系数,并且,利用所述分解系数得到原子序数分布图像和电子密度分布图像。
在本发明的多能CT成像系统的成像方法中,代替所述步骤(h)而具有如下的步骤(j):对在所述步骤(g)中得到的衰减系数重建结果f i 进行正投影处理,从而得到补全后的分角度扫面缺失的投影数据,然后,利用前处理方法得到不同的基函数的系数,从而得到原子序数分布图像和电子密度分布图像。
根据本发明,在不改变现有的传统单能CT系统硬件设计的前提下,将X射线发生装置的圆轨道分为预定的角度区间,对每个角度区间设定预定的X射线发生装置的高压,利用本发明中的重建方法提供双能成像功能,扩展了传统CT的使用功能,并且,与现有的双能CT设备相比,成本较低。
附图说明
图1是示出用于实现双能CT系统的双源双探测器方式的示意图。
图2是示出用于实现双能CT系统的单源双层探测器方式的示意图。
图3是示出用于实现双能CT系统的能量快速切换双能CT的示意图。
图4是示出本发明的多能CT成像系统的示意图。
图5是利用本发明的双能扫描得到的正弦图。
图6是对图5进行卷积所得到的结果。
具体实施方式
以下,参照附图对本发明的实施方式进行说明。
图4是本发明的多能CT成像系统的示意图,(A)是系统结构以及多能扫描方式的示意图,(B)是机架旋转角度与X光机电压之间的关系的图。如图4所示,在本发明的多能CT成像系统中,具备:载物台,用于载置受检查物品(即,图4中的待测物体);用于发出照射上述受检查物品的X射线的电压可调的X射线发生装置(即,图4中的X光源);探测器,接收从上述X射线发生装置发出且穿过了上述受检查物品的X射线,并且,输出探测数据(例如,不同能量下的投影数据等);机架;数据处理及控制装置,对上述载物台、上述X射线发生装置、上述探测器、上述机架进行控制,处理来自上述探测器的上述探测数据。
此外,在本发明中所使用的X射线发生装置是具有可调电压功能的X光机,能够根据数据处理及控制装置的控制来切换电压,发出不同能量的X射线。此外,关于本发明中所使用的X射线发生装置,并没有特别限定,只要是能够根据数据处理及控制装置的控制进行电压切换的X光源,就能够应用于本发明。
此外,X射线发生装置以及探测器装配在环形的机架上,机架能够连续旋转,并且,通过滑环机构完成与X射线发生装置的控制信号的传递以及探测器所采集的数据的传输。此外,还能够对X射线发生装置进行供电(此处,滑环的功能与发电机电刷类似)。此外,例如,在CT成像系统的运行当中,机架始终进行匀速旋转运动,如果X射线发生装置进行圆周旋转时所在的平面为XY平面(此处例如为纸面),则受检查物品在传送皮带上沿着Z轴方向直线运动(即,垂直于纸面运动)。
此外,在X射线发生装置的一周扫描过程中,数据处理及控制装置根据预定的能量数将圆轨道平均分成预定的N(N是大于等于2的整数)个角度区间,对每一个角度区间设定不同的X射线发生装置的高压,当机架从当前的角度区间旋转进入下一个角度区间时,根据来自数据处理及控制装置的指令,对X射线发生装置进行控制,将电压切换为在下一个角度区间设定的电压,由此,发射出不同能量的X射线。
载物台、机架分别做圆轨道旋转运动和直线平移运动或者两者的组合,从而能够实现螺旋扫描轨迹或者圆周扫描轨迹。具体地说,如图4(A)所示,本发明中的载物台是放置物体的结构。当机架进行圆周旋转而载物台静止、或者机架静止而载物台进行圆周运动时,能够实圆周轨道扫描方式,这两者是等价的。此外,螺旋运动能够分解为XY平面内的圆周运动和沿Z方向的直线运动的这两个分运动,使载物台和机架分别完成其中的一个分运动,将它们组合起来就是螺旋运动,例如载物台进行直线运动而机架进行圆周运动。此外,也可以是载物台进行旋转而机架进行平移,这种方式适用于大型的物体检测,这是因为大型物体平移需要太大场地空间而旋转较容易实现。
此外,由于载物台与机架的运动属于相对运动,所以,在本发明中,按照使载物台固定并且使机架旋转的圆周扫描轨迹进行描述。
此外,数据处理及控制装置由例如计算机等构成,用于完成多能CT成像系统运行过程的控制,这些控制包括机械转动、电气控制、安全联锁控制等。在X射线发生装置的一周扫描过程中,根据所需要的能量数将X射线发生装置的圆轨道平均分成多个角度区间,例如,在双能的情况下将圆轨道分为两个180°角度区间,即,分为0°~180°的角度区间以及180°~360°的角度区间,在三能的情况下(图4中所示的情况)将圆轨道分为三个120°角度区间,即,分为0°~120°的角度区间、120°~240°的角度区间以及240°~360°的角度区间。对每一个角度区间设定不同的X光机高压,例如,在将圆轨道分为两个180°角度区间的情况下,X光机高压分别为120KV以及80KV。此外,当机架从当前的角度区间旋转进入下一个角度区间时,数据处理及控制装置发送指令,对X光机进行控制,将X光机电压切换为在下一个角度区间设定的电压,由于此单步切换过程能够快速完成,所以,并不会影响数据采集。
此外,在本发明中所使用的探测器是完整的面阵X射线探测器,但是,探测器也不限于此,只要能够完成数据探测,能够使用任意的探测器。此外,在本发明的多能CT成像系统中还具有读出电路、采集触发信号电路及数据传输电路,这些电路与面阵X射线探测器一起组成对数据进行采集的数据采集系统。
以下,对本发明的多能CT成像系统中的成像方法进行说明。
本发明的多能CT成像系统在X射线发生装置的圆轨道扫描一圈角度范围内分不同能量来采集数据。例如,预先选定N个能量(记作E1,E2,…,En),将圆轨道平均分为N个角度区间(也就是对N个角度区间的每一个设定预定的电压),N是大于等于2的整数。在此情况下,每个能量下的投影数据只覆盖圆轨道的360°/N的范围,从而可以看出这属于有限角度重建的范畴。
(1)首先,多能CT成像系统进行工作,数据处理及控制装置进行控制,使X射线发生装置沿着圆轨道进行旋转,并且,X射线发生装置在由数据处理及控制装置预先划分的每个角度区间以在该角度区间所设定的电压进行工作,发射出X射线,照射受检查物品。相应地,利用由探测器、读出电路、采集触发信号电路以及数据传输电路组成的数据采集系统对数据进行采集,由此,得到不同能量E1,E2,…,En(n是大于等于2的整数)下的投影数据,并且,将不同能量E1,E2,…,En下的投影数据进行本底校正、不一致性校正、负对数运算等操作,从而得到物体在能量E1,E2,…,En下的线积分数据P1,P2,…,Pn。
将上述的所有能量E1,E2,…,En下的线积分数据P1,P2,…,Pn依次进行拼接,由此,能够得到覆盖X射线发生装置的圆轨道360°的完整的线积分数据。接着,利用普通圆轨道滤波反投影算法等的方法进行重建,由此,能够得到混合重建结果f0。此外,关于混合重建结果f0,因为不同能量下的结构信息不同,也就是说f0是多个不完整信息混在一起的一个中间结果,所以没有严格的物理含义。此外,虽然拼接得到的线积分数据在角度方向上有不连续性,但是,在滤波方向上是连续的,因此,进行反投影所得到的结果能够准确地反映物体(受检查物品)的几何结构信息。
此外,为了便于理解,如图5所示,以二维为例进行说明。在图5所示的例子中,将圆轨道分为0°~180°和180°~360°这两个角度区间(即,双能扫描)。在此处,拼接得到的结果指的是CT中常说的正弦图。关于上述的拼接,就是不管能量上的差异,将这些数据(线积分数据P1,P2,…,Pn)按照单能CT的数据构成方式按照角度的次序放到一起,所以,在图5中可以看到,在180°时切换能量,在能量改变的地方有一个明显的跳跃边缘(即,在180°处会出现明显的不连续现象)。即,在角度这个方向(图5中的纵轴)数据是不连续的。另一方面,虽然这两个角度区间的能量不同,但是,各部分都包含被照射物体(受检查物品)的结构信息。此外,由于滤波是在探测器方向上做一维滤波,即对图5中的线积分数据一行一行独立地做一维卷积,所以,180°的阶跃不会影响卷积结果,在图6中示出卷积的结果。
(2)为了进行各个能量E1,E2,…,En下的衰减系数的有限角度范围的重建,将如上所述那样求出的混合重建结果f0作为先验图像。所谓先验图像,是指与待重建的图像之间具有一定联系并且能够反映出待重建图像的某些信息(例如,结构、边缘、特征等)的图像。在本发明中,先验图像是指由图5的线积分数据所重建出来的图像。此外,由于线积分数据是由不同能量的数据拼接得到的,所以,重建结果不完全等于任何一个能量下的等效衰减系数,但是都比较接近,并且在几何结构上能够完全反应物体信息。
接着,为了得到各个能量E1,E2,…,En下的衰减系数重建结果f i (i=1,2,…,n),构造先验约束的优化问题,使得,即,这表示寻找一个参数f i ,在满足的前提下,使得函数达到全局极小值;或者,另一种等价的描述为:在所有满足条件的变量f i 中,寻找使得为最小的那个f i 。其中,为描述先验图像与重建结果的差异的函数,此处,例如能够定义为或者等。H i 为第i个能量下的投影算子并且它只包含第i个能量下覆盖的角度范围所对应的部分,为与测量数据噪声方差有关的参数。此外,利用拉格朗日乘数法将上述问题转化为无约束优化问题,并且,利用非线性共轭梯度法求解,得到各个能量E1,E2,…,En下的衰减系数重建结果f i (i=1,2,…,n)。
由于先验图像f0具有完整的边缘和结构信息,此外,有限角度线积分数据包含了真实的衰减系数信息。因此,通过将这两者结合,从而能够同时保证几何结构与衰减系数的准确性。
(3)如以上所述那样,得到了不同能量E1,E2,…,En下的衰减系数的重建结果f1,f2,…,fn,能够使用后处理基函数分解方法得到分解系数,进而能够得到原子序数和电子密度。
通常,在双能分解重建原子序数和电子密度中,用表示线性衰减系数,可以是某个能量值E的衰减系数,也可以是X光机高压为E时所得到的等效衰减系数,能够根据不同方法来决定采用哪个定义,物质对X射线的衰减主要有光电效应和康普顿效应,这两种效应能够分解为与能量相关的光电函数和康普顿函数(Klein-Nishina公式),如下述的式(1)所示:
……(1),
其中,、为常数,为阿伏伽德罗常数,光电函数以及康普顿函数是已知函数,两个分解系数为Z3.5 ρ e 以及ρ e 。此外,在基效应方法中就是利用上述关系得到原子序数Z和电子密度ρ e 的。前处理方法是首先得到分解系数的线积分,再重建出分解系数,然后再得到和,后处理方法则是直接由衰减系数图像求出分解系数,然后再得到和。
此外,根据双能理论的假设,X射线的衰减系数也可以分解为两种已知材料的线性衰减系数的组合,如下述的式(2)所示:
……(2),
其中,、为已知材料的线性衰减系数函数,例如碳和铝、骨头和软组织。分解方法也分别前处理和后处理,在得到分解系数和之后,能够通过以下公式(3)计算和,
……(3),
其中,、、、为已知的分解材料的原子序数和电子密度。
在本发明中,以双能为例,后处理基材料分解模型如以下的式(4)所示:
……(4),
其中,为所选用的已知的两种基材料在高低能等效能量下的衰减系数,根据以上的式(4)能够得到基材料分解系数。式(4)中的基材料分解系数c 1、c 2相当于上述式(2)中的b 1、b 2,因此,能够进一步根据上述的式(3)得到原子序数Z和电子密度ρ e ,进而能够得到原子序数分布图像和电子密度分布图像。
(4)此外,在上述的(2)中已经重建出不同能量E1,E2,…,En下的衰减系数图像,对这些图像进行正投影运算,从而能够得到补全了分角度扫描缺失的投影数据。在得到了完整的投影数据之后,能够利用前处理方法得到不同基函数的系数。此处,在本发明中以双能为例,前处理基效应分解模型为如下的式(5):
……(5),
在式(5)中,为双能有限角度范围直接测量投影值与重建图像离散投影值拼接得到的完整投影数据,为光电效应系数和康普顿效应系数的线积分,为不同能量的能谱与探测器响应函数的乘积,称为等效能谱。根据上述式(5)能够得到两种效应的线积分,再利用传统FBP算法即可重建得到光电系数和康普顿系数,进而得到原子序数和电子密度。由此,能够得到原子序数分布图像和电子密度分布图像。
具体地说,在已经得到了不同能量E1,E2,…,En下的衰减系数重建结果f1,f2,…,fn的情况下,对n个能量补全各自缺失角度下的投影数据。此处,以采用两个能量将圆轨道分为0~180度(记作“第一角度区间”)和180~360度(记作“第二角度区间”)这两个角度区间为例。第一个能量下(即,第一角度区间对应的能量)的数据缺180~360度的投影数据,所以,利用所得到的f1进行投影计算180~360度的数据来进行补全,第二个能量下(即,第二角度区间对应的能量)的数据缺0~180度的投影数据,所以,利用得到的f2进行投影计算0~180度的数据来进行补全,即,在分角度扫描缺失的情况下补全投影数据,在得到了完整的投影数据之后,利用前处理方法得到不同基函数的系数。
对完整的投影数据进行基函数分解,得到基函数对应的系数,也就是式(5)中的是根据光电效应和康普顿效应取基函数的结果。实际上,也能够取某些材料随能谱变化的衰减系数函数作为基函数,即本领域内的基材料分解法。分解后所得到的数据是投影域的数据,利用传统FBP算法就能够得到系数的分布,进而根据系数的分布能够计算得到原子序数分布和电子密度分布。
此外,通常在医疗应用中优选使用骨和软组织这两种材料进行后处理分解,相对于此,在安检应用中优选使用光电效应和康普顿效应的前处理分解。
如上所述,根据本发明的多能CT成像系统以及成像方法,在将X射线发生装置的圆轨道分为多个角度区间的情况下也能够准确地得到待测物体(受检查物品)的原子序数Z和电子密度ρ e ,进而能够得到原子序数分布图像以及电子密度分布图像。此外,由于本发明的多能CT成像系统没有增加X射线发生装置或者探测器的数量,所以,成本不会增加,但是能够提供多能成像功能。此外,本发明并不象现有技术的快速能量切换系统那样需要装配新型X光机,所以,能够被广泛应用。
以下,根据本发明设计双能CT成像系统,并举出如下的实施例。
例如,将X光机的圆轨道分为0~180度和180~360度这两个角度区间,并且,在0~180度角度区间的X光机高压设定为120KV,在180~360度角度区间的X光机高压设定为80KV。
首先,将本发明的双能CT成像系统中的X光机高压设定为120KV,开机扫描。
接着,双能CT成像系统的机架连续旋转,并且,以1°为采样间隔采集投影数据。
接着,当机架旋转到180°时,数据处理及控制装置发出指令,迅速将X光机高压切换至80KV,继续扫描,并且扫描采样间隔仍为1°。
接着,当机架旋转到360°回到初始位置时,停止旋转,同时关闭光机,数据采集完成。
接着,双能CT成像系统的数据处理及控制装置对在X光机高压120KV和80KV下得到的线积分数据进行拼接,得到覆盖360度的完整的线积分数据,并且,根据该完整的线积分数据得到作为先验图像的混合重建结果f0。
接着,根据先验图像f0分别得到120KV和80KV下的衰减系数重建结果f1、f2。
接着,根据上述的重建结果f1、f2得到受检查物品的原子序数Z和电子密度ρ e ,根据原子序数Z和电子密度ρ e 显示原子序数分布图像或者电子密度分布图像,或者也可以显示各个不同能量下的衰减系数图像。
根据本申请发明,在不改变现有的传统单能CT系统硬件设计的情况下,将X射线发生装置的圆轨道分为预定的角度区间,对每个角度区间设定预定的X射线发生装置的高压,利用本发明中的重建方法提供多能成像功能(特别是双能CT功能),并且,与现有的双能CT设备相比,成本较低。
Claims (6)
1. 一种多能CT成像系统,其特征在于,具有:
载物台,用于载置受检查物品;
电压可调的X射线发生装置,用于发出照射所述受检查物品的X射线;
探测器,用于接收从所述X射线发生装置发出且穿过所述受检查物品的X射线,并且,输出探测数据;
机架,装配有所述X射线发生装置以及所述探测器;以及
数据处理及控制装置,对所述载物台、所述X射线发生装置、所述探测器、所述机架进行控制,处理来自所述探测器的所述探测数据,在所述X射线发生装置的一周扫描过程中,所述数据处理及控制装置根据预定的能量数将所述X射线发生装置的圆轨道平均分成预定个数的角度区间,对每一个角度区间设定不同的所述X射线发生装置的高压,当所述机架从当前的角度区间旋转进入下一个角度区间时,根据来自所述数据处理及控制装置的指令,对所述X射线发生装置进行控制,将所述X射线发生装置切换为在下一个角度区间设定的电压。
2. 一种多能CT成像系统的成像方法,其特征在于,具有如下步骤:
(a)按照预定的能量数将所述多能CT成像系统的X射线发生装置的圆轨道平均分成预定个数的角度区间,对每个角度区间设定不同的所述X射线发生装置的高压;
(b)所述X射线发生装置沿着圆轨道进行旋转,在每个角度区间以在该角度区间所设定的高压进行工作,发出X射线,照射受检查物品;
(c)利用所述多能CT成像系统的探测器采集不同能量下的投影数据;
(d)利用所述多能CT成像系统的数据处理及控制装置得到所述不同能量下的投影数据的线积分数据;
(e)将所述不同能量下的线积分数据依次进行拼接,由此得到覆盖360°的完整的线积分数据;
(f)根据所述完整的线积分数据进行重建,由此得到混合重建结果;
(g)将所述混合重建结果作为先验图像,得到不同能量下的衰减系数重建结果;
(h)根据所述衰减系数重建结果得到分解系数,并且,利用所述分解系数得到原子序数分布图像和电子密度分布图像。
3. 如权利要求2所述的多能CT成像系统的成像方法,其特征在于,
在所述步骤(a)中将圆轨道分为N角度区间并且使在所述步骤(d)中得到的不同能量下的线积分数据为P1,P2,…,Pn,并且,在所述步骤(f)中利用普通圆轨道滤波反投影法进行重建,使得到的混合重建结果为f0,N和n是大于等于2的整数,
在所述步骤(g)中,利用 得到满足条件的衰减系数重建结果f i ,其中,G(f i ,f0)为描述先验图像与重建结果的差异的函数,H i 为第i个能量下的投影算子,ε为与测量数据噪声方差有关的参数,i=1,2,…,n。
4. 如权利要求3所述的多能CT成像系统的成像方法,其特征在于,
在所述步骤(h)中,使用后处理基函数分解法得到分解系数,并且,利用所述分解系数得到原子序数分布图像和电子密度分布图像。
5. 如权利要求3所述的多能CT成像系统的成像方法,其特征在于,
在所述步骤(h)中,使用后处理基材料分解法得到分解系数,并且,利用所述分解系数得到原子序数分布图像和电子密度分布图像。
6. 如权利要求3所述的多能CT成像系统的成像方法,其特征在于,
代替所述步骤(h)而具有如下的步骤(j):对在所述步骤(g)中得到的衰减系数重建结果f i 进行正投影处理,从而得到补全后的分角度扫描缺失的投影数据,然后,利用前处理方法得到不同的基函数的系数,从而得到原子序数分布图像和电子密度分布图像。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310008543.3A CN103913779B (zh) | 2012-12-31 | 2013-01-10 | 多能ct成像系统以及成像方法 |
PCT/CN2013/001576 WO2014101285A1 (zh) | 2012-12-31 | 2013-12-17 | 多能ct成像系统以及成像方法 |
EP13198192.0A EP2749219B1 (en) | 2012-12-31 | 2013-12-18 | Multi-energy CT imaging system and imaging method |
US14/135,023 US9326743B2 (en) | 2012-12-31 | 2013-12-19 | Multi-energy CT imaging system and imaging method |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210589336 | 2012-12-31 | ||
CN201210589336.7 | 2012-12-31 | ||
CN2012105893367 | 2012-12-31 | ||
CN201310008543.3A CN103913779B (zh) | 2012-12-31 | 2013-01-10 | 多能ct成像系统以及成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103913779A true CN103913779A (zh) | 2014-07-09 |
CN103913779B CN103913779B (zh) | 2017-07-07 |
Family
ID=49943111
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310008543.3A Active CN103913779B (zh) | 2012-12-31 | 2013-01-10 | 多能ct成像系统以及成像方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US9326743B2 (zh) |
EP (1) | EP2749219B1 (zh) |
CN (1) | CN103913779B (zh) |
WO (1) | WO2014101285A1 (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104545960A (zh) * | 2014-12-18 | 2015-04-29 | 沈阳东软医疗系统有限公司 | 一种双能ct的扫描电压的选择方法及设备 |
CN105806858A (zh) * | 2014-12-31 | 2016-07-27 | 北京固鸿科技有限公司 | Ct检测方法和ct设备 |
CN106226333A (zh) * | 2016-08-01 | 2016-12-14 | 中国科学院深圳先进技术研究院 | 一种微聚焦x光源类同轴相衬成像自动化系统 |
CN106910157A (zh) * | 2017-02-17 | 2017-06-30 | 公安部第研究所 | 一种多级并行的图像重建方法及装置 |
CN107822652A (zh) * | 2016-09-15 | 2018-03-23 | 西门子医疗有限公司 | 用于重建光谱结果图像数据的方法 |
CN110720939A (zh) * | 2019-10-23 | 2020-01-24 | 明峰医疗系统股份有限公司 | 一种用于智能ma变换的滑环角度校准方法 |
CN110742635A (zh) * | 2019-10-08 | 2020-02-04 | 南京安科医疗科技有限公司 | 一种复合能谱ct成像方法 |
WO2020051915A1 (zh) * | 2018-09-14 | 2020-03-19 | 西安大医集团有限公司 | 一种放疗设备及其控制方法和装置 |
CN111413358A (zh) * | 2020-03-30 | 2020-07-14 | 赵磊 | 扫描组件的设置方法及其扫描组件和扫描设备 |
CN111433635A (zh) * | 2017-10-19 | 2020-07-17 | 约翰内斯堡大学 | 伽玛射线层析摄影 |
CN111552002A (zh) * | 2020-05-19 | 2020-08-18 | 重庆大学 | 一种三源摆动螺旋安检ct成像装置及方法 |
CN111812130A (zh) * | 2020-07-15 | 2020-10-23 | 深圳市金园智能科技有限公司 | 一种基于x射线的材料内部3d成像方法及装置 |
CN111896566A (zh) * | 2020-07-20 | 2020-11-06 | 上海交通大学 | 一种增加同步辐射光源成像范围的装置及方法 |
CN112690817A (zh) * | 2020-12-28 | 2021-04-23 | 明峰医疗系统股份有限公司 | 基于双能ct的数据采集方法、系统及计算机可读存储介质 |
CN112964738A (zh) * | 2021-01-29 | 2021-06-15 | 山东大学 | 一种工业ct快速扫描系统及方法 |
CN113687434A (zh) * | 2021-09-29 | 2021-11-23 | 北京航星机器制造有限公司 | 一种安检ct图像重建的物体扫描位置确定方法及装置 |
CN113885087A (zh) * | 2021-09-29 | 2022-01-04 | 北京航星机器制造有限公司 | 一种安检ct的图像重建方法 |
CN115105107A (zh) * | 2022-06-20 | 2022-09-27 | 北京朗视仪器股份有限公司 | 一种能谱成像方法及能谱成像系统 |
CN116046815A (zh) * | 2023-02-21 | 2023-05-02 | 上海福柯斯智能科技有限公司 | 一种双能ct成像方法、装置及系统 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10169909B2 (en) * | 2014-08-07 | 2019-01-01 | Pixar | Generating a volumetric projection for an object |
CN104483711B (zh) * | 2014-12-17 | 2020-02-21 | 同方威视技术股份有限公司 | 基于分布式光源的辐射成像系统 |
CN106551702B (zh) * | 2015-09-30 | 2021-05-14 | 通用电气公司 | Ct图像的射束硬化伪影校正装置及方法 |
CN105510362B (zh) * | 2015-12-07 | 2019-01-08 | 华中农业大学 | 基于微型ct的水稻分蘖性状无损测量装置及其测量方法 |
CN108267464B (zh) * | 2016-12-29 | 2020-03-10 | 同方威视技术股份有限公司 | 液体探测方法和装置及安检设备 |
EP3431007B1 (en) | 2017-07-21 | 2020-06-17 | Koninklijke Philips N.V. | Creation of electron density datasets from spectral ct datasets |
US10610173B2 (en) * | 2018-01-16 | 2020-04-07 | General Electric Company | System and method to improve spatial resolution in computed tomography |
JP7169853B2 (ja) * | 2018-11-09 | 2022-11-11 | キヤノン株式会社 | 画像処理装置、放射線撮影装置、および画像処理方法 |
CN110353716B (zh) * | 2019-07-25 | 2023-06-27 | 沈阳开普医疗影像技术有限公司 | 双能ct短程扫描协议中的数字域分解方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1502309A (zh) * | 2002-11-27 | 2004-06-09 | GEҽ��ϵͳ����������˾ | 便于减少伪影的方法和装置 |
CN1765327A (zh) * | 2004-10-25 | 2006-05-03 | 西门子公司 | 产生多重能量图像的断层造影设备及其方法 |
US20090097611A1 (en) * | 2007-10-15 | 2009-04-16 | Ge Medical Systems Global Technology Company, Llc | X-ray ct apparatus and method for processing x-ray projection data |
CN101772324A (zh) * | 2007-08-06 | 2010-07-07 | 株式会社日立医药 | X射线ct装置 |
CN203149136U (zh) * | 2012-12-31 | 2013-08-21 | 清华大学 | 多能ct成像系统 |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5661774A (en) * | 1996-06-27 | 1997-08-26 | Analogic Corporation | Dual energy power supply |
US6236709B1 (en) * | 1998-05-04 | 2001-05-22 | Ensco, Inc. | Continuous high speed tomographic imaging system and method |
US6999549B2 (en) * | 2002-11-27 | 2006-02-14 | Ge Medical Systems Global Technology, Llc | Method and apparatus for quantifying tissue fat content |
CN100435733C (zh) * | 2005-12-31 | 2008-11-26 | 清华大学 | X-ct扫描系统 |
DE102006009222B4 (de) * | 2006-02-28 | 2008-02-28 | Siemens Ag | Verfahren und Vorrichtung zur Bestimmung der Konzentration einer Substanz in einem Körpermaterial mittels Mehr-Energie-Computertomographie |
US7298812B2 (en) * | 2006-03-31 | 2007-11-20 | General Electric Company | Image-based material decomposition |
JP5107568B2 (ja) * | 2006-12-12 | 2012-12-26 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
US7945013B2 (en) * | 2007-02-14 | 2011-05-17 | Hitachi Medical Corporation | X-ray CT apparatus |
US7697657B2 (en) * | 2007-03-23 | 2010-04-13 | General Electric Company | System and method of density and effective atomic number imaging |
JP5525706B2 (ja) * | 2007-07-25 | 2014-06-18 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
US7813474B2 (en) * | 2007-10-09 | 2010-10-12 | General Electric Company | Method and apparatus for performing dual-spectrum CT with fast KV modulation at multiple-view intervals |
US20090129539A1 (en) * | 2007-11-21 | 2009-05-21 | General Electric Company | Computed tomography method and system |
US8165264B2 (en) * | 2009-01-28 | 2012-04-24 | Kabushiki Kaisha Toshiba | Method of pre-reconstruction decomposition for fast kV-switching acquisition in dual energy computed tomography (CT) |
JP5535532B2 (ja) * | 2009-06-26 | 2014-07-02 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
-
2013
- 2013-01-10 CN CN201310008543.3A patent/CN103913779B/zh active Active
- 2013-12-17 WO PCT/CN2013/001576 patent/WO2014101285A1/zh active Application Filing
- 2013-12-18 EP EP13198192.0A patent/EP2749219B1/en active Active
- 2013-12-19 US US14/135,023 patent/US9326743B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1502309A (zh) * | 2002-11-27 | 2004-06-09 | GEҽ��ϵͳ����������˾ | 便于减少伪影的方法和装置 |
CN1765327A (zh) * | 2004-10-25 | 2006-05-03 | 西门子公司 | 产生多重能量图像的断层造影设备及其方法 |
CN101772324A (zh) * | 2007-08-06 | 2010-07-07 | 株式会社日立医药 | X射线ct装置 |
US20090097611A1 (en) * | 2007-10-15 | 2009-04-16 | Ge Medical Systems Global Technology Company, Llc | X-ray ct apparatus and method for processing x-ray projection data |
CN203149136U (zh) * | 2012-12-31 | 2013-08-21 | 清华大学 | 多能ct成像系统 |
Non-Patent Citations (1)
Title |
---|
高洋 等: "结合全变差最小化的双能CT重建", 《计算机应用研究》 * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104545960A (zh) * | 2014-12-18 | 2015-04-29 | 沈阳东软医疗系统有限公司 | 一种双能ct的扫描电压的选择方法及设备 |
US10105110B2 (en) | 2014-12-18 | 2018-10-23 | Shenyang Neusoft Medical Systems Co., Ltd. | Selecting scanning voltages for dual energy CT scanning |
US10175182B2 (en) | 2014-12-31 | 2019-01-08 | Grandpect Company Limited | CT detection method and CT device |
CN105806858A (zh) * | 2014-12-31 | 2016-07-27 | 北京固鸿科技有限公司 | Ct检测方法和ct设备 |
CN105806858B (zh) * | 2014-12-31 | 2019-05-17 | 北京固鸿科技有限公司 | Ct检测方法和ct设备 |
CN106226333A (zh) * | 2016-08-01 | 2016-12-14 | 中国科学院深圳先进技术研究院 | 一种微聚焦x光源类同轴相衬成像自动化系统 |
CN107822652B (zh) * | 2016-09-15 | 2019-07-12 | 西门子医疗有限公司 | 用于重建光谱结果图像数据的方法 |
CN107822652A (zh) * | 2016-09-15 | 2018-03-23 | 西门子医疗有限公司 | 用于重建光谱结果图像数据的方法 |
CN106910157B (zh) * | 2017-02-17 | 2020-06-26 | 公安部第一研究所 | 一种多级并行的图像重建方法及装置 |
CN106910157A (zh) * | 2017-02-17 | 2017-06-30 | 公安部第研究所 | 一种多级并行的图像重建方法及装置 |
CN111433635A (zh) * | 2017-10-19 | 2020-07-17 | 约翰内斯堡大学 | 伽玛射线层析摄影 |
CN111433635B (zh) * | 2017-10-19 | 2023-09-15 | 约翰内斯堡大学 | 伽玛射线层析摄影 |
WO2020051915A1 (zh) * | 2018-09-14 | 2020-03-19 | 西安大医集团有限公司 | 一种放疗设备及其控制方法和装置 |
CN111194183A (zh) * | 2018-09-14 | 2020-05-22 | 西安大医集团有限公司 | 一种放疗设备及其控制方法和装置 |
WO2021068643A1 (zh) * | 2019-10-08 | 2021-04-15 | 南京安科医疗科技有限公司 | 一种复合能谱ct成像方法 |
CN110742635A (zh) * | 2019-10-08 | 2020-02-04 | 南京安科医疗科技有限公司 | 一种复合能谱ct成像方法 |
CN110742635B (zh) * | 2019-10-08 | 2021-10-08 | 南京安科医疗科技有限公司 | 一种复合能谱ct成像方法 |
CN110720939A (zh) * | 2019-10-23 | 2020-01-24 | 明峰医疗系统股份有限公司 | 一种用于智能ma变换的滑环角度校准方法 |
CN111413358A (zh) * | 2020-03-30 | 2020-07-14 | 赵磊 | 扫描组件的设置方法及其扫描组件和扫描设备 |
CN111552002A (zh) * | 2020-05-19 | 2020-08-18 | 重庆大学 | 一种三源摆动螺旋安检ct成像装置及方法 |
CN111812130A (zh) * | 2020-07-15 | 2020-10-23 | 深圳市金园智能科技有限公司 | 一种基于x射线的材料内部3d成像方法及装置 |
CN111896566A (zh) * | 2020-07-20 | 2020-11-06 | 上海交通大学 | 一种增加同步辐射光源成像范围的装置及方法 |
CN112690817A (zh) * | 2020-12-28 | 2021-04-23 | 明峰医疗系统股份有限公司 | 基于双能ct的数据采集方法、系统及计算机可读存储介质 |
CN112964738A (zh) * | 2021-01-29 | 2021-06-15 | 山东大学 | 一种工业ct快速扫描系统及方法 |
US11821853B2 (en) | 2021-01-29 | 2023-11-21 | Shandong University | Fast industrial CT scanning system and method |
CN113687434A (zh) * | 2021-09-29 | 2021-11-23 | 北京航星机器制造有限公司 | 一种安检ct图像重建的物体扫描位置确定方法及装置 |
CN113885087A (zh) * | 2021-09-29 | 2022-01-04 | 北京航星机器制造有限公司 | 一种安检ct的图像重建方法 |
CN113687434B (zh) * | 2021-09-29 | 2024-04-09 | 北京航星机器制造有限公司 | 一种安检ct图像重建的物体扫描位置确定方法及装置 |
CN113885087B (zh) * | 2021-09-29 | 2024-05-24 | 北京航星机器制造有限公司 | 一种安检ct的图像重建方法 |
CN115105107A (zh) * | 2022-06-20 | 2022-09-27 | 北京朗视仪器股份有限公司 | 一种能谱成像方法及能谱成像系统 |
CN116046815A (zh) * | 2023-02-21 | 2023-05-02 | 上海福柯斯智能科技有限公司 | 一种双能ct成像方法、装置及系统 |
CN116046815B (zh) * | 2023-02-21 | 2023-11-03 | 上海福柯斯智能科技有限公司 | 一种双能ct成像方法、装置及系统 |
Also Published As
Publication number | Publication date |
---|---|
US20140185741A1 (en) | 2014-07-03 |
EP2749219B1 (en) | 2021-08-25 |
US9326743B2 (en) | 2016-05-03 |
WO2014101285A1 (zh) | 2014-07-03 |
EP2749219A1 (en) | 2014-07-02 |
CN103913779B (zh) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103913779A (zh) | 多能ct成像系统以及成像方法 | |
CN203149136U (zh) | 多能ct成像系统 | |
CN100482164C (zh) | 便于减少伪影的方法和装置 | |
EP1324696B1 (en) | System and method for cone beam volume computed tomography using circle-plus-multiple-arc orbit | |
US7778383B2 (en) | Effective dual-energy x-ray attenuation measurement | |
CN103384498B (zh) | 探测装置 | |
RU2619839C1 (ru) | Способ исследования с помощью кт и устройство для кт | |
RU2595796C2 (ru) | Аппарат для обнаружения, содержащий два сцинтиллятора для обнаружения рентгеновского излучения | |
US8548118B2 (en) | Apparatus and method for spectral projection imaging with fast KV switching | |
CN103472074A (zh) | Ct成像系统和方法 | |
CN110179486B (zh) | 多能量ct成像系统及其应用 | |
CN100376211C (zh) | 运动检查对象的计算机断层造影中设置图像的方法及设备 | |
CN104422704A (zh) | 对x光能谱ct的能谱信息进行分解的方法和对应的重建方法 | |
CN1742296A (zh) | 锥形光束计算机断层成像的立体像素驱动的螺旋重建方法 | |
CN105488826A (zh) | 一种基于fbp的能谱ct迭代成像方法和成像系统 | |
CN103620393B (zh) | 成像装置 | |
JP2008529637A (ja) | 螺旋相対運動及び円錐ビーム束を用いたコンピュータ断層撮影方法 | |
CN110428478B (zh) | 交替光源扇束x射线ct采样方法及装置 | |
CN109632844A (zh) | 基于直线扫描轨迹的双能ct成像系统及方法 | |
US9173624B2 (en) | System and method for multi-energy X-ray imaging using a polychromatic source | |
CN105806860A (zh) | 一种全场x射线成像系统及成像方法 | |
JP2007508560A (ja) | コヒーレント散乱コンピュータ断層撮影法 | |
CN101453951A (zh) | 锥束ct半周期闭合螺旋轨迹 | |
CN106097411B (zh) | Ct机图像重建方法及高分辨ct扫描机 | |
CN104323790A (zh) | 同轴相衬成像方法及系统和相衬ct方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: DE Ref document number: 1199762 Country of ref document: HK |
|
GR01 | Patent grant | ||
GR01 | Patent grant | ||
REG | Reference to a national code |
Ref country code: HK Ref legal event code: GR Ref document number: 1199762 Country of ref document: HK |