CN107592935A - X射线ct的自动***校准方法 - Google Patents
X射线ct的自动***校准方法 Download PDFInfo
- Publication number
- CN107592935A CN107592935A CN201680026497.1A CN201680026497A CN107592935A CN 107592935 A CN107592935 A CN 107592935A CN 201680026497 A CN201680026497 A CN 201680026497A CN 107592935 A CN107592935 A CN 107592935A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- projection
- parameter
- vector
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 158
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 31
- 238000003325 tomography Methods 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims description 68
- 238000005070 sampling Methods 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 17
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 claims description 16
- 230000005855 radiation Effects 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 2
- 238000012804 iterative process Methods 0.000 abstract description 4
- 230000033001 locomotion Effects 0.000 description 27
- 238000002591 computed tomography Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 10
- 230000009467 reduction Effects 0.000 description 9
- 210000001015 abdomen Anatomy 0.000 description 8
- 230000001186 cumulative effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000002059 diagnostic imaging Methods 0.000 description 5
- 230000006872 improvement Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000005291 magnetic effect Effects 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000007408 cone-beam computed tomography Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000010603 microCT Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000002247 constant time method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005294 ferromagnetic effect Effects 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000013519 translation Methods 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
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- 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
-
- 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
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- 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/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
提供了使用迭代重建算法在计算机断层(CT)扫描中的几何校准和图像重建的***和方法。迭代重建算法能够用来重建经改善的图像,然后所述经改善的图像能够用来通过使用局部线性嵌入(LLE)方法调整不准确的参数。进而能够使用经调整的参数来重建新图像,然后能够用该新图像来进一步调整参数。能够重复该迭代过程的步骤,直到满足质量阈值。
Description
相关申请的交叉引用
本申请要求于2015年3月20日提交的序列号为62/135861的美国临时申请的权益,故通过引用的方式将包括任何图形、表格和附图在内的上述申请的全部内容并入本文。
政府支持
本发明是在美国国立卫生研究院授予的第R01EBO 16977和U01 EB017140号拨款的政府支持下进行的。政府对本发明有一定的权利。
背景技术
在医学成像领域,X射线计算机断层扫描(CT)提供了关键的诊断信息。近年来,已经开发了CT技术,来以低辐射剂量提供高质量的图像。患者运动或不准确的机器能够导致X射线源与被扫描对象的中心之间不准确的投影角度或距离,从而导致不准确的结果。这一领域的另一个挑战性问题是几何校准,例如C型臂CT和超高分辨率CT的几何校准。这个问题也与刚性患者运动补偿有关,因为运动在成像组件与患者身体之间是相对的。
为了进行几何校准和运动校正,最近已经提出了许多方法。已经提出了具有校准模体的分析方法和具有或不具有校准模体的迭代方法。分析方法广泛应用于工业CT,并且能够基于对锥束几何中的椭圆参数的识别。一些校准方法是迭代的,例如针对锥束CT的基于优化的校准、针对锥束CT的自校准和针对锥束微CT的自校准。在有关几何校准和运动减少的文献中有共同之处。一些运动减少方法甚至使用多源检测器***或者在避免受运动影响的数据时运用快速扫描,而其他方法估计患者运动并补偿其影响。然而,这些方法中的每一种都有局限性。
发明内容
本发明提供用于使用一个或更多个迭代重建算法在计算机断层(CT)扫描(例如X射线CT扫描)中进行几何校准和图像重建的新颖且有利的***和方法。能够使用迭代重建算法(例如基于全变差(Total Variation))以重建经改善的图像,然后能够使用经改善的图像,例如通过使用局部线性嵌入(LLE)方法,来调整不准确的参数。进而能够使用经调整的参数重建新图像,然后能够用该新图像来进一步调整参数。能够重复该迭代过程的步骤,直到满足质量阈值。这能够导致针对X射线CT的自动***校准,而不是现有技术中的基于准确但昂贵的机器的参数调整方法和***。
在实施方案中,重建CT图像的方法能够包括:获得初始CT图像;对初始CT图像执行重建算法以获得初始CT图像的重建图像;使用重建图像(例如使用LLE方法)来调整与图像相关联的一个或更多个参数;使用经调整的一个或更多个参数对重建图像执行重建算法以获得更新的重建图像;迭代地重复参数调整和图像重建的更新(每次使用最新更新的重建图像,并且每次更新相同的一个或更多个参数),直到满足预定特征的阈值。
在另一实施方案中,CT***能够包括:辐射源;用于检测来自辐射源的辐射的检测器;以及包括计算机可读介质的计算机,该计算机可读介质具有存储在其上的、用于执行本发明的方法的计算机可执行指令。
附图说明
图1示出根据本发明的实施方案的基于局部线性嵌入(LLE)的校准的流程图。
图2A示出在进行任何补偿之前的Shepp-Logan模体的图像。
图2B示出根据本发明的实施方案在25次校准迭代之后的图2A的模体的图像。
图3示出根据本发明的实施方案的均方根误差(RMSE)相对于校准的迭代的曲线图。
图4示出根据本发明的实施方案的平均角度误差相对于校准迭代的曲线图。
图5A示出角度误差相对于角度数的曲线图。
图5B示出根据本发明的实施方案的校准的角度误差相对于角度数的曲线图。
图6示出根据本发明的实施方案的平均原始源对象距离(OS)误差相对于校准迭代的曲线图。
图7A示出OS误差相对于角度数的曲线图。
图7B示出根据本发明的实施方案的校准的OS误差相对于角度数的曲线图。
图8示出说明光谱计算机断层(CT)的示意图。
图9示出在任何补偿之前胸部模体的原始图像。
图10A示出在任何补偿之前胸部模体的图像。
图10B示出根据本发明的实施方案在6次校准迭代之后的图10A的模体图像。
图11示出根据本发明的实施方案的校准的角度误差相对于角度数的曲线图。
图12示出根据本发明的实施例的校准的OS误差相对于角度数的曲线图。
图13示出说明X射线CT的扇束几何的示意图。
图14示出说明扇束X射线CT中的对象或患者移动的示意图。
图15示出来自实际CT扫描的腹部模体图像。
图16A示出图15所示的腹部模体在任何补偿之前的重建图像。
图16B示出根据本发明的实施方案在1次校准迭代之后的图15所示的腹部模体的重建图像。
图16C示出根据本发明的实施方案在2次校准迭代之后的图15所示的腹部模体的重建图像。
图16D示出根据本发明的实施方案在4次校准迭代之后的图15所示的腹部模体的重建图像。
图17A示出根据本发明的实施方案的重建图像的通用质量指数(UQI)相对于校准迭代的曲线图。
图17B示出根据本发明的实施方案的重建图像的计算成本(以小时计)相对于校准迭代的曲线图。
图17C示出根据本发明的实施方案的在重建图像角度60处的累积成本相对于校准迭代的曲线图。
图18A示出根据本发明的实施方案的重建图像的平均角度误差(以度计)相对于校准迭代的曲线图。
图18B示出根据本发明的实施方案的重建图像的源-对象距离(SOD)(以cm计)相对于校准迭代的曲线图。
图18C示出根据本发明的实施方案的重建图像的对象-检测器距离(ODD)(以cm计)相对于校准迭代的曲线图。
图18D示出根据本发明的实施方案的重建图像的检测器倾斜误差(以度计)相对于校准迭代的曲线图。
图18E示出根据本发明的实施方案的重建图像的检测器偏移误差(以cm计)相对于校准迭代的曲线图。
图19A示出在校准“迭代0”处(即在校准开始之前)的投影角度误差(以度计)相对于投影角度数的曲线图。
图19B示出根据本发明的实施方案在1次校准迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图。
图19C示出根据本发明的实施方案在2次校准迭代之后的投影角误差(以度计)相对于投影角度数的曲线图。
图19D示出根据本发明的实施方案在4次校准迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图。
图20示出说明具有未对齐旋转台的扇束X射线CT几何结构的示意图。
图21A示出在CT扫描期间未对准的对象在任何补偿之前的重建图像。
图21B示出根据本发明的实施方案在1次校准迭代之后的图21A中描绘的对象的重建图像。
图21C示出根据本发明的实施方案在2次校准迭代之后的图21A中描绘的对象的重建图像。
图21D示出根据本发明的实施方案在5次校准迭代之后的图21A中描绘的对象的重建图像。
图22A示出根据本发明的实施方案的重建图像的计算成本(以小时计)相对于校准迭代的曲线图。
图22B示出根据本发明的实施方案的重建图像在角度60处的累积成本相对于校准迭代的曲线图。
图22C示出根据本发明的实施方案的重建图像的校准旋转中心偏移(以cm计)相对于校准迭代的曲线图。
图23A示出使用CT扫描器扫描的对象在任何补偿之前的重建图像。
图23B示出根据本发明的实施方案在1次校准迭代之后图23A中描绘的对象的重建图像。
图23C示出根据本发明的实施方案在10次校准迭代之后图23A中描绘的对象的重建图像。
图23D示出根据本发明的实施方案在20次校准迭代之后图23A中描绘的对象的重建图像。
图24A示出根据本发明的实施方案的UQI相对于校准迭代的曲线图。
图24B示出根据本发明的实施方案的平均角度误差(以度计)相对于校准迭代的曲线图。
图24C示出根据本发明的实施方案的平均x坐标偏移误差(以cm计)相对于校准迭代的曲线图。
图24D示出根据本发明的实施方案的平均y坐标偏移误差(以cm计)相对于校准迭代的曲线图。
图24E示出根据本发明的实施方案的重建图像的计算成本(以小时计)相对于校准迭代的曲线图。
图24F示出根据本发明的实施方案的重建图像在角度60处的累积成本相对于校准迭代的曲线图。
图25A示出在任何校准之前(迭代0)的x坐标偏移误差(以cm计)相对于投影角度数的曲线图。
图25B示出根据本发明的实施方案在1次校准迭代之后的x坐标偏移误差(以cm计)相对于投影角度数的曲线图。
图25C示出根据本发明的实施方案在10次校准迭代之后的x坐标偏移误差(以cm计)相对于投影角度数的曲线图。
图25D示出根据本发明的实施方案在20次校准迭代之后的x坐标偏移误差(以cm计)相对于投影角度数的曲线图。
图26A示出在任何校准之前(迭代0)的y坐标偏移误差(以cm计)相对于投影角度数的曲线图。
图26B示出根据本发明的实施方案在1次校准迭代之后的y坐标偏移误差(以cm计)相对于投影角度数的曲线图。
图26C示出根据本发明的实施方案在10次校准迭代之后的y坐标偏移误差(以cm计)相对于投影角度数的曲线图。
图26D示出根据本发明的实施方案在20次校准之后的y坐标偏移误差(cm)相对于投影角数的曲线图。
图27A示出在任何校准前(迭代0)的投影角度误差(以度计)相对于投影角度数的曲线图。
图27B示出根据本发明的实施方案在1次校准迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图。
图27C示出根据本发明的实施方案在10次校准迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图。
图27D示出根据本发明的实施方案在20次校准迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图。
具体实施方式
本发明提供使用一个或更多个迭代重建算法在计算机断层(CT)扫描(例如X射线CT扫描)中进行几何校准和图像重建的新颖且有利的***和方法。能够使用迭代重建算法(例如基于全变差)以重建经改善的图像,然后能够使用经改善的图像,例如通过使用局部线性嵌入(LLE)方法,来调整不准确的参数。进而能够使用经调整的参数重建新图像,然后能够用该新图像来进一步调整参数。能够重复该迭代过程的步骤,直到满足质量阈值。这能够导致针对X射线CT的自动***校准,而不是现有技术中的基于准确但昂贵的机器的参数调整方法和***。因此,本发明的方法和***能够显著降低X射线CT的成本。此外,***几何参数能够在实际条件下直接从投影数据中提取。
本发明的实施方案能够使用数学方法调整CT扫描中的不准确参数,从而有助于重建质量更好的CT图像。因为该方法是一种迭代方法,在实施方案中,能够使用并行计算(例如GPU方法),并且能够减少计算时间。
在本发明的许多实施方案中,能够通过结合LLE和重投影对CT扫描(例如X射线CT扫描)执行几何校准(例如迭代更新的参数)。初始参数化能够用于重建CT图像,并且能够计算对参数范围进行密集地采样的经重投影的投影向量。利用经重投影的投影向量和原始投影向量,能够计算LLE的权重系数和近邻,并且能够更新参数估计。能够迭代地重建图像,直到达到令人满意的质量(例如,直到满足感兴趣的一个或更多个特征中的一个或更多个阈值)。在某些实施方案中,LLE和重投影能够用于扇束成像几何CT扫描。
LLE是2000年提出一种用于降维的无监督流形学习算法(Roweis等人,“Nonlineardimensionality reduction by locally linear embedding”,Science290.5500(2000):2323-2326;Roweis等人,“An introduction to locally linear embedding”,http://www.cs.toronto.edu/~roweis/lle/publications.html(2000);两者全部内容均通过引用并入本文)。与经典降维方法例如主成分分析(PCA)和多维尺度分析(MDS)相关,LLE计算特征向量,但是在本地计算以保持嵌入在高维空间中的本征几何。LLE易于实现,并且在数据足够密集时具有出色的性能。LLE能够包括三个步骤:步骤1能够包括根据欧几里德距离找到每个高维点的K个最近邻;步骤2能够包括用每个点的近邻线性地表示每个点并计算其近邻的权重系数;步骤3能够包括将具有权重系数的高维数据映射到本征流形上的低维表示。
在本发明的许多实施方案中,利用数据向量bi,LLE的第一步骤能够是根据欧几里德距离找到其数据集向量中的K个最近邻:
利用K个最近邻,第二步骤能够是用其近邻向量线性地表示原始数据向量:
其中是K个最近邻,wik是相应的权重系数,并且能够通过最小化以下误差求解权重矩阵:
使用约束求解方程式(3)等同于求解线性组:
其中C=(cjk)是局部协方差矩阵,其被计算为:
利用权重系数W=(wik),第三步骤能够包括通过求解方程:
来计算全局内坐标Y=(yi)。
通过Rayleitz-Ritz定理,方程式(6)的解由广义特征值问题的底部d+1特征向量给出:
MY=λY, (7)
其中λ是特征值,并且M=(I-W)T(I-W)。
图13示出说明用于X射线CT的扇束几何的示意图,该X射线CT广泛用于商业螺旋CT***。参考图13,在扇束几何中,笛卡尔坐标系的原点能够是待扫描和待重建的对象的标称中心,并且X射线源旨在处于其中心为原点的圆形轨迹上。源至对象的距离(SOD)是X射线源与***原点之间的距离,即扫描圆的半径。X射线源发射覆盖对象的扇束X射线。扇束的中心射线能够垂直于长度为L的线性检测器阵列,其在被对象衰减之后接收X射线,并且感兴趣的射线与中心射线之间的角度能够表示为γ。对象至检测器的距离(ODD)是***原点与检测器中点之间的距离。符号θ表示X射线扫描角度,并且β=θ-γ。
由于X射线源是在圆形轨迹上,所以在极坐标系中能够方便地表示扇束几何。各种实际因素能够导致视角依赖型方式中不准确的几何参数。具体地,不准确的投影角度和其他参数能够包括:
其中θi是准确的投影角度;pj是其他准确的几何参数,包括SOD、ODD、检测器偏移和检测器倾斜角度,其不依赖于θi。N是视角的数量,M是其他参数的数量,δi和εi分别是角度误差和参数误差。通常,更准确的几何参数会产生重建更好CT图像的能力。
在扇束几何CT扫描期间,几何校准中能够存在刚性患者运动问题。图14示出说明扇束X射线CT中的对象或患者移动的示意图。参考图14,能够用对象x坐标偏移xi、对象y坐标偏移yi和投影角度误差来表示患者运动,而投影角度仍然能够由方程式(8)表示。
在本发明的许多实施方案中,能够执行使用LLE的校准。能够对迭代CT重建方法进行建模,以求解以下线性方程组:
Au=bi, (10)
其中u=(u1,u2,...,uJ)是表示为J维向量的图像,J是像素数,bi=(b1,b2,...,bL)是数据,L是向量元素的数量,A=(ajk)是与几何参数相关的投影矩阵。如果***矩阵和投影数据均是已知的,则能够使用迭代算法通过减少原始数据和重投影数据之间的差异而重建具有固定投影矩阵的CT图像。在重建算法中,重要步骤是计算投影矩阵,其受到投影角度和SOD等几何参数的影响;也就是说,
A=A(P), (11)
其中P是估计参数向量,其包括SOD、ODD和检测器偏移距离、检测器倾斜角度、和/或投影角度等。
在实施方案中,能够使用距离驱动法计算投影,例如Siddon所公开的距离驱动法(Fast calculation of the exact radiological path for a three-dimensional CTarray,Medical Physics 12(2)(1985),252-255),其全部内容通过引用并入本文。基于迭代重建中的重投影方法,能够制定一种新的迭代方法来通过最小化投影数据和经重投影的投影数据之间的均方误差而估计几何参数,其能够表示为:
其中bi是从沿着不同投影视角的测量而获得的投影向量,是来自具有采样的参数的重建图像的相应的经重投影的投影向量,P是更新的参数向量。能够在密集采样的参数范围:
内,由方程式(10)和(11)计算经重投影的投影向量。
如果参数采样间隔足够小,则真实参数向量接近相邻的采样的参数向量,并且所测量的投影向量能够由与采样的参数向量相关联的K个最近的经重投影的投影向量线性地表示。也就是说,
其中是与参数的相应的K个向量相关联的K个最近重投影向量,并且wik是权重系数。关键的关系是方程式(14)和(15)的权重系数是相同的。因此,能够通过搜索K个最近的经重投影的投影向量并用权重系数和相应的采样的参数更新参数向量来改进真实参数估计。因此,几何校准问题能够通过LLE的降维而解决。
利用密集采样的参数域和相应的经重投影的投影向量,能够找到原始投影向量的K个最近的经重投影的投影向量,并且能够通过使用方程式(1)、(4)和(5)计算权重系数。利用采样的参数和相应的权重系数,能够根据方程式(15)执行参数更新。
图1示出根据本发明的实施方案的基于LLE的校准的流程图,如前文段落所讨论的。
在某些实施方案中,能够使用通用质量指数(UQI)来评估重建图像和/或作为阈值特性以确定应该何时停止迭代重建过程。UQI在Wang等人的文献(A universal imagequality index,Signal Processing Letters,IEEE 9(3)(2002),81–84)中被详细地描述,其全部内容通过引用并入本文。在某些实施方案中,能够使用平均角误差参数误差(PError)和/或平均对象坐标偏移和量化几何校准结果。这些能够与UQI一起使用或与一个或更多个其他阈值特性一起使用。
UQI通过整合三个因素来评估图像,包括相关失真、亮度失真和对比度失真。UQI的范围为-1至1。UQI越接近1,重建图像会越好。给定重建图像u*ij和S×T的本地真实图像uij,UQI定义为:
其中
平均角度误差为:
其中是校准后的投影角度,θi是原始投影角度。参数误差是:
PError=|P-POri|, (23)
其中P是校准后的参数向量,POri是真实参数向量。平均对象x坐标偏移和y坐标偏移是:
其中xi和yi分别是原始对象x坐标偏移和y坐标偏移,而和分别为经校准的偏移。
在许多实施方案中,能够使用有序子集同时代数重建技术(OS-SART)来执行图像重建。在Wang等人的文献(Ordered-subset simultaneous algebraic reconstructiontechniques,Journal of X-ray Science and Technology,12(3)(2004),169-177)中论述了OS-SART,其全部内容通过引用并入本文。此外,能够使用全变差(TV)正则化。在Sidky等人的文献(Accurate image reconstruction from few-views and limited-angle datain divergent beam CT,Journal of X-ray Science and Technology,14(2)(2006),119–139)中详细论述了TV,其全部内容通过引用并入本文。
在实施方案中,CT***能够包括辐射源(例如X射线源)、用于从辐射源检测辐射(例如X射线)的检测器、以及计算机***和/或一个或更多个计算机可读介质,其上存储有用于执行本申请所公开的方法的计算机可执行指令。例如,该计算机可执行指令能够执行用于CT(例如X射线CT)的迭代重建算法和/或几何校准。
本发明的实施方案提供了视角式(view-wise)随机几何参数的CT***(例如X射线CT***,如二维X射线CT***)的校准。LLE能够作为评估具有固有的低维度一致性的几何参数的重要步骤,这已证明被能够提供显著的改进(参见实施例)。LLE能够使用相应的采样的参数找到K个最近的经重投影的投影向量,并且基于具有经由LLE计算的权重系数和采样的参数的线性组合来更新几何参数。
本申请公开的统一图像重建和参数估计方案能够包括通过重投影和LLE迭代地更新投影矩阵和底层图像。这能够应用于CT***(例如X射线CT***)校准和刚性患者运动补偿。虽然已经具体讨论了扇束几何和刚性患者运动,但是这仅仅是为了示例性的目的,不应理解为限制性的。本发明的实施方案能够与其他类型的几何结构(例如锥束几何)和其他类型的运动(例如非刚性患者运动)一起使用。
本发明的实施方案能够使用一个或更多个迭代重建算法(例如基于全变差)来重建更好的图像,然后使用经重建的图像通过LLE方法调整不准确的参数。进而能够使用经调整的参数重建新图像,然后能够使用新图像调整不准确的参数。能够重复该迭代过程(例如直到满足一个或更多个特征(例如UQI)的一个或更多个阈值)以获得期望的参数。这能够产生CT(例如X射线CT)的自动***校准,从而通过消除对某些昂贵机器的需求而允许经济上的节省。不准确的参数能够用数学方法进行调整,从而能够重建质量更好的CT图像。
对于(X射线)CT,几何校准和患者运动补偿(例如刚性患者运动补偿)是优化图像重建质量的相关联问题。在CT扫描期间,未校准的***几何和患者运动会导致重建图像中的条纹状模糊和其他伪影。本申请所公开的LLE方法解决了该挑战,并且能够在刚性二维对象假设下进行,由此以比任何现有的相关技术方法更通用的方式来解决几何校准和患者运动补偿的挑战。投影能够通过LLE向上采样的近邻线性表示,并且能够从投影数据迭代地估计CT***参数。
本发明的方法和***能够用于当前的CT***中,特别是在具有不准确的机器(因此需要参数调整)的低成本CT***中。在现有相关技术的CT***中,大量的资金花在准确的机器上,导致小型医院和诊所可能无法负担高成本的CT***。本发明的数学参数调整方法和***不需要高成本和高精度机器,能够仅使用计算机设备(例如,如计算机和/或一个或更多个计算机可读介质的处理设备)来调整参数。通过允许降低CT***的生产成本,本发明的方法和***对于CT制造商也是有价值的,这能够导致开发低成本***和/或便携式CT***,能够增加总体CT***的销售。
本申请所描述的方法和过程能够被实现为代码和/或数据。本申请所描述的软件代码和数据能够存储在一个或更多个计算机可读介质上,其可以包括能够存储供计算机***使用的代码和/或数据的任何设备或介质。当计算机***读取并执行存储在计算机可读介质上的代码和/或数据时,计算机***执行体现为存储在计算机可读存储介质内的数据结构和代码的方法和过程。
本领域技术人员应当理解,计算机可读介质包括能够用于存储诸如计算机可读指令、数据结构、程序模块和计算***/环境所使用的其它数据的信息的可移除和不可移除的结构/设备。计算机可读介质包括但不限于:易失性存储器,例如随机存取存储器(RAM、DRAM、SRAM);非易失性存储器,例如闪速存储器、各种只读存储器(ROM、PROM、EPROM、EEPROM)、磁性和铁磁/铁电存储器(MRAM、FeRAM)、以及磁性和光存储设备(硬盘驱动器、磁带、CD、DVD);网络设备;或现在已知或未来开发的能够存储计算机可读信息/数据的介质。计算机可读介质不应该被理解或解释为包括任何传播信号。本发明的计算机可读介质能够是例如光盘(CD)、数字视频盘(DVD)、闪速存储器设备、易失性存储器、或硬盘驱动器(HDD),如外部HDD或计算设备的HDD,但是实施方案不限于此。计算设备能够是例如膝上型计算机、台式计算机、服务器、手机、或平板电脑,但是实施方案不限于此。
本发明包括但不限于以下示例性实施方案。
实施方案1.一种重建计算机断层扫描(CT)图像的方法,所述方法包括:
i)获得初始CT图像;
ii)对所述初始CT图像执行重建算法以获得所述初始CT图像的重建图像;
iii)使用重建图像来调整与所述图像(例如所述初始CT图像或先前的重建图像)或与从其获得初始CT图像的CT***相关联的一个或更多个参数;以及
iv)使用经调整的一个或更多个参数对重建图像执行重建算法以获得更新的重建图像;
v)重复步骤iii)至iv)(在每个重复的步骤iii)中使用最近更新的重建图像,并且在每个重复的步骤iii))中更新相同的参数,直到达到预定特征的阈值。
实施方案2.根据实施方案1所述的方法,其中,所述重建算法是基于全变差(TV)的。
实施方案3.根据实施方案1至2中任一项所述的方法,其中,使用所述重建图像来调整与所述图像相关联的一个或更多个参数包括对所述一个或更多个参数中的每个参数使用局部线性嵌入(LLE)方法来完成。
实施方案4.根据实施方案3所述的方法,其中所述LLE方法包括:
步骤1)根据欧几里德距离找出每个高维点的K个最近邻;
步骤2)用每个点的近邻线性地表示每个点,并计算其近邻的权重系数;
步骤3)将具有所述权重系数的高维数据映射到本征流形上的低维表示。
实施方案5.根据实施方案3至4中任一项所述的方法,其中,所述LLE方法包括:
根据所述欧几里德距离:
在数据向量bi的数据集向量中找到其K个最近邻。
实施方案6.根据实施方案3至5中任一项所述的方法,其中,所述LLE方法(还)包括:
使用K个最近向量,用其近邻向量:
来线性地表示原始数据向量,
其中是K个最近邻,wik是各自的权重系数,并且能够通过最小化以下误差来求解权重矩阵:
实施方案7.根据实施方案6所述的方法,其中,使用约束来最小化所述误差等同于求解线性组:
其中C=(cjk)是局部协方差矩阵,其被计算为:
实施方案8.根据实施方案3至7中任一项所述的方法,其中,所述LLE方法(还)包括利用权重系数W=(wik),通过求解方程:
来计算全局内坐标Y=(yi)。
实施方案9.根据实施方案8所述的方法,其中,能够被求解以计算Y=(yi)的方程式的解由广义特征值问题:
MY=λY
的底部的d+1特征向量给出。
其中λ是特征值,并且M=(I-W)T(I-W)。
实施方案10.根据实施方案3至9中任一项所述的方法,其中在所述LLE方法中使用的最近邻的数量K=2。
实施方案11.根据实施方案3至9中任一项所述的方法,其中在所述LLE方法中使用的最近邻的数量至少为2。
实施方案12.根据实施方案1至11中任一项所述的方法,其中,必须满足其阈值才能停止重复步骤iii)至iv)的预定特征是通用质量指数(UQI)。
实施方案13.根据实施方案12所述的方法,其中,所述阈值为0.6。
实施方案14.根据实施方案12所述的方法,其中,所述阈值至少为0.6。
实施方案15.根据实施方案12所述的方法,其中,所述阈值为0.7。
实施方案16.根据实施方案12所述的方法,其中,所述阈值至少为0.7。
实施方案17.根据实施方案12所述的方法,其中,所述阈值为0.8。
实施方案18.根据实施方案12所述的方法,其中,所述阈值至少为0.8。
实施方案19.根据实施方案12所述的方法,其中,所述阈值为0.9。
实施方案20.根据实施方案12所述的方法,其中,所述阈值至少为0.9。
实施方案21.根据实施方案12所述的方法,其中,所述阈值为0.95。
实施方案22.根据实施方案12所述的方法,其中,所述阈值至少为0.95。
实施方案23.根据实施方案12所述的方法,其中,所述阈值为0.5。
实施方案24.根据实施方案12所述的方法,其中,所述阈值至少为0.5。
实施方案25.根据实施方案1至24中任一项所述的方法,其中,所述一个或更多个参数包括以下各项中的至少一个:投影角度、源到对象的距离(SOD)、对象到检测器的距离(ODD)、检测器偏移、检测器倾斜角度、对象x坐标偏移、对象y坐标偏移、和投影角度误差。
实施方案26.根据实施方案1至24中任一项所述的方法,其中,所述一个或更多个参数包括以下各项中的至少两个:投影角度、SOD、ODD、检测器偏移、检测器倾斜角度、对象x坐标偏移、对象y坐标偏移、和投影角度误差。
实施方案27.根据实施方案1至24中任一项所述的方法,其中,所述一个或更多个参数包括以下各项中的至少三个:投影角度、SOD、ODD、检测器偏移、检测器倾斜角度、对象x坐标偏移、对象y坐标偏移、和投影角度误差。
实施方案28.根据实施方案1至24中任一项所述的方法,其中,所述一个或更多个参数包括以下各项中的至少四个:投影角度、SOD、ODD、检测器偏移、检测器倾斜角度、对象x坐标偏移、对象y坐标偏移、和投影角度误差。
实施方案29.根据实施方案1至24中任一项所述的方法,其中,所述一个或更多个参数包括投影角度、SOD、ODD、检测器偏移、和检测器倾斜角度。
实施方案30.根据实施方案1至24中任一项所述的方法,其中,所述一个或更多个参数包括对象x坐标偏移、对象y坐标偏移、和投影角误差。
实施方案31.根据实施方案3至30中任一项所述的方法,其中,所述LLE方法包括求解以下线性方程组:
Au=bi,
其中u=(u1,u2,...,uJ)是表示为J维向量的图像,J是像素数,bi=(b1,b2,...,bL)是数据,J是向量元素的数量,A=(ajk)是与所述一个或更多个(几何)参数相关的投影矩阵。
实施方案32.根据实施方案31所述的方法,其中,执行重建算法包括计算受到一个或更多个(几何)参数影响的投影矩阵:
A=A(P),
其中P是包括所述一个或更多个参数的估计参数向量。
实施方案33.根据实施方案32所述的方法,其中,使用距离驱动法计算所述投影矩阵的投影。
实施方案34.根据实施方案32至33中任一项所述的方法,其中,使用所述重建图像来调整一个或更多个参数包括:最小化所述投影数据和经重投影的投影数据之间的均方误差,其能够表示为:
其中bi是从沿着不同投影视角的测量获得的投影向量,是来自具有采样的参数的重建图像的相应的经重投影的投影向量,P是更新的参数向量。
实施方案35.根据实施方案34所述的方法,其中能够通过在密集采样的参数范围:
内,用实施方案31和32中提供的方程式来计算经重投影的投影向量。
实施方案36.根据实施方案35所述的方法,其中,真实参数向量接近相邻的采样的参数向量,并且所测量的投影向量能够由与采样的参数向量相关联的K个最近的经重投影的投影向量来线性地表示,使得:
和
其中是与参数的相应的K个向量相关联的K个最近的重投影向量,而wik是权重系数。
实施方案37.根据实施方案1至36中任一项所述的方法,其中,获得所述初始CT图像包括获得初始参数向量,所述初始参数向量包括与所述图像相关联的所述一个或更多个参数的初始值(例如所述初始CT图像或所述先前的重建图像)或从其获得初始CT图像的CT***。
实施方案38.根据实施方案37所述的方法,执行重建算法包括基于(最近)重建图像,用一个或更多个参数的更新值来更新参数向量。
实施方案39.根据实施方案1至38中任一项所述的方法,其中,所述重建算法包括使用有序子集同时代数重建技术(OS-SART)。
实施方案40.根据实施方案1至39中任一项所述的方法,其中,所述重建算法包括使用TV进行正则化。
实施方案41.根据实施方案1至40中任一项所述的方法,其中,步骤ii)由处理器执行。
实施方案42.根据实施方案1至41中任一项所述的方法,其中,步骤iii)由处理器执行。
实施方案43.根据实施方案1至42中任一项所述的方法,其中,步骤iv)由处理器执行。
实施方案44.根据实施方案1至43中任一项所述的方法,其中,步骤v)由处理器执行。
实施方案45.一种计算机可读介质,其具有用于执行根据实施方案1至44中任一项所述的方法的计算机可执行指令(存储在其上)。
实施方案46.一种CT***,其包括:
辐射源(例如X射线源);
用于从辐射源检测辐射(例如X射线)的检测器;
具有根据实施方案45所述的计算机可读介质的计算机。
从以下通过说明的方式给出的实施例能够更好地理解本发明及其许多优点。以下实施例说明了本发明的一些方法、应用、实施方案和变体。当然,它们不认为是限制本发明。能够对本发明进行许多改变和修改。
示例1
使用来自诊所的腹部图像模体评估基于LLE的校准方法,如图15所示。腹部模体大小为512×512。使用该模体,生成具有角度和其他几何误差的投影数据。然后,将本文描述的算法应用于具有几何参数误差的投影数据。该算法在刚性患者运动问题中进行了测试。计算环境是在具有Intel Core 2Duo [email protected]、4.00GB RAM和64位操作***的计算机上的Matlab 2010a。
为了评估重建图像,使用了UQI。此外,几何校准结果用平均角度误差参数误差PError和平均对象坐标偏移量和进行量化。
图像重建过程使用具有TV正则化的OS-SART进行。作为重建和校准迭代停止标准,在图像和投影域中均监测到令人满意的UQI值。
评估了校准方法用于减少投影角度和其他几何误差的效用。投影角度数N=360。投影角度如下随机扰动:
θi=360×(i-1)/N+δi,i=1,2,...,N。 (26)
所使用的其他参数如表1所示,包括实际值、初始值和最终估计值。探测器阵列的长度为80cm。使用几何参数,假设入射X射线光子数为105,产生具有泊松噪声的投影数据。
为了有效地校准参数,首先将几何参数逐一校准,而不是将所有参数一起校准。校准过程按以下顺序进行:检测器偏移;SOD;ODD;检测器倾斜;投影角度。每个参数的角度采样步长分别为0.04cm、0.1cm、0.1cm、0.01°、0.02°。采样范围分别为[0cm,0.4cm]、[46cm,54cm]、[46cm,54cm]、[0°,2°]和[-1°,1°]。最近邻的数量为2。重建的停止UQI阈值为0.9,校准的停止UQI阈值为0.6。
图16A示出图15在任何补偿之前的重建图像;图16B至16D分别示出图15在几何校准的1次迭代、2次迭代和4次迭代之后的重建图像。校准参数如表1所示。
表格1.实际几何、初始、和校准参数
图17A示出重建图像的UQI相对于校准迭代的曲线图;图17B示出重建图像相对于校准迭代的计算成本(以小时计)的曲线图;图17C示出重建图像相对于校准迭代的在角度60°处的累积成本的曲线图。图18A示出重建图像相对于校准迭代的平均角度误差(以度计)的曲线图;图18B示出重建图像相对于校准迭代的SOD(以cm计)的曲线图;图18C示出重建图像与校准迭代的ODD(以cm计)的曲线图;图18D示出重建图像相对于校准迭代的检测器倾斜误差(以度计)的曲线图;图18E示出重建图像相对于校准迭代的检测器偏移误差(以cm计)的曲线图。图19A示出在校准“迭代0”处(即在校准开始之前)的投影角度误差(以度计)相对于投影角度数的曲线图;图19B示出在1次校准迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图;图19C示出在2次校准迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图;图19D示出在4次校准迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图。
参考图17A至17C、18A至18E和19A至19D以及表1,能够看出,重建图像的所有参数和特征都通过校准过程得到显著改善。对于许多参数,只需要少量迭代就看到剧烈的改善。
示例2
在具有几何参数误差的实际投影数据的情况下,应用了本文所讨论的校准算法。标称SOD和ODD值均为38cm,检测器元件数为1024,检测器长度为13.0048cm。在角度范围[0°,360°]中的投影角度数为900。如图20所示,当旋转中心不在X射线源和检测器中心之间的线上时,对象在旋转台上,其中红点(阴影圆内包含的斑点)是旋转中心,其不在***源上。因此,该算法用于校准旋转中心偏移。采样范围为[-1cm,1cm],采样步长为0.05cm。最近邻的数量为2。重建的UQI停止阈值为0.9,校准的UQI停止阈值为0.6。
图21A示出在任何补偿之前的重建图像。图21B示出在1次校准迭代之后的图21A所示的重建图像;图21C示出在2次迭代之后的重建图像;图21D示出在5次迭代之后的重建图像。另外,图22A示出重建图像相对于校准迭代的计算成本(以小时计)的曲线图;图22B示出重建图像相对于校准迭代在角度60°处的累积成本的曲线图;图22C示出重建图像相对于校准迭代的校准旋转中心偏移(以cm计)的曲线图。校准旋转中心偏移为-0.3847cm。
参考图21A至21D和22A至22C,能够看出,重建图像的所有参数和特征都通过校准过程得到显著改善。对于许多参数,只需要少量迭代就看到剧烈的改善。
示例3
针对涉及刚性患者运动的情况,测试了本文所讨论的校准。如讨论的,刚性患者运动能够表示为对象x坐标偏移误差、对象y坐标偏移误差和投影角度误差。用腹部模体模拟刚性患者运动校准。投影角度等参数与实施例1相同。对象x坐标偏移误差和对象y坐标偏移误差均在[-1cm,1cm]的范围内随机选择,投影角度误差在[-1°,1°]的范围内随机选择。
校准按以下顺序:对象x坐标偏移;对象y坐标偏移;投影角度。参数的采样率分别为0.05cm、0.05cm和0.1°。采样范围分别为[-1cm,1cm]、[-1cm,1cm]和[-1°,1°]。最近邻的数量K=2。重建的UQI停止值为0.9,校准的UQI停止值为0.8。
图23A示出在任何补偿之前的重建图像。图23B示出在1次校准迭代之后的图23A的重建图像;图23C示出在10次迭代之后的重建图像;图23D示出在20次迭代之后的重建图像。图24A示出根据本发明的实施方案的UQI相对于校准迭代的曲线图;图24B示出平均角度误差(以度计)相对于校准迭代的曲线图;图24C示出平均x坐标偏移误差(以cm计)相对于校准迭代的曲线图;图24D示出平均y坐标偏移误差(以cm计)相对于校准迭代的曲线图;图24E示出重建图像相对于校准迭代的计算成本(以小时计)的曲线图;图24F示出重建图像相对于校准迭代在角度60°处的累积成本的曲线图。
图25A示出x坐标偏移误差(以cm计)相对于在任何校准之前(迭代0)的投影角度数的曲线图;图25B示出在1次校准迭代之后x坐标偏移误差(以cm计)相对于投影角度数的曲线图;图25C示出在10次校准迭代之后的x坐标偏移误差(以cm计)相对于投影角度数的曲线图;图25D示出在20次迭代之后x坐标偏移误差(以cm计)相对于投影角度数的曲线图。图26A示出了在任何校准之前(迭代0)的y坐标偏移误差(以cm计)相对于投影角度数的曲线图;图26B示出在1次校准迭代之后的y坐标偏移误差(以cm计)相对于投影角度数的曲线图;图26C示出在10次校准迭代之后的y坐标偏移误差(以cm计)相对于投影角度数的曲线图;图26D示出在20次校准迭代之后的y坐标偏移误差(以cm计)相对于投影角度数的曲线图。图27A示出在任何校准之前(迭代0)的投影角度误差(以度计)相对于投影角度数的曲线图;图27B示出在1次校准迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图;图27C示出在10次迭代之后的投影角度误差(以度计)相对于投影角度数的曲线图;图27D示出在20次迭代之后的投影角度误差(以度计)相对于投影角度数的关系曲线。
参考图23A至23D、24A至24F、25A至25D、26A至26D和27A至27D,能够看出,重建图像的所有参数和特性都通过校准过程得到显著改善。对于许多参数,只需要少量迭代就看到剧烈的改善。
几何参数对重建质量具有不同的敏感度。在几何校准中,检测器偏移显著地影响重建。检测器倾斜和投影角度误差比SOD误差和ODD误差对重建质量具有更强的影响。这有助于解释检测器偏移、检测器倾斜和投影角度的校准结果优于SOD和ODD的校准结果。在患者运动校准中,对象x和y坐标偏移对于重建图像的质量比投影角度误差更有效。
密切相关的问题是重投影的参数抽样率。采样率越高,校准结果越准确,但计算成本越高。选择采样率能够考虑到校准结果与计算开销之间的平衡。
使用TV的重要问题是其重量。如果太小,TV会无法减少伪影和噪音。如果太大,TV会使CT图像过度平滑。TV参数取决于伪影和噪音的水平。在实施例1至3中,几何校准后TV参数经验性地设定为0.1。在噪音较弱的情况下,能够使用较小的TV参数。
示例4
使用Shepp-Logan模体评估基于LLE的校准方法。图2A示出在任何校准之前的模体,图2B示出在25次校准迭代之后的重建图像。图3示出均方根误差(RMSE)相对于校准迭代的曲线图,图4示出平均角度误差相对于校准迭代的曲线图。图5A示出该校准的角度误差相对于角度数的曲线图,图5B示出校准的角度误差相对于角度数的曲线图。图6示出平均原始SOD(OS)误差相对于校准迭代的曲线图;图7A示出OS误差相对于角度数的曲线图;图7B示出校准的OS误差相对于角度数的曲线图。
参考图2A至2B、3、4、5A至5B、6和7A至7B,能够看出,重建图像的所有参数和特性都通过校准过程得到显著改善。对于许多参数,只需要少量迭代就看到剧烈的改善。
示例5
使用具有105个光子的泊松噪声的胸部模体评估基于LLE的校准方法。图9示出了在该实施例中使用的胸部模体图像。图10A示出在任何补偿之前的胸部模体图像,图10B示出在6次校准迭代之后的模体图像。图11示出校准的角度误差相对于角度数的曲线图,图12示出校准的OS误差相对于角度数的曲线图。
参考图10A至10B、11和12,能够看出,校准过程改善了重建图像。
应当理解,本文所述的实施例和实施方案仅用于说明目的,并且本领域技术人员将提出各种修改或改变,其包括在本申请的精神和范围内。
本文中参考或引用的所有专利、专利申请、临时申请和出版物(包括“参考”部分中的那些)均以其全部内容通过引用并入,包括所有图形和表格,直到其与本说明书的明确教导不一致的程度。
参考文献
-G.Wang,H.Yu,and B.D.Man."An outlook on x-ray CT research anddevelopment."Medical physics 35.3(2008):1051-1064.
-G.Wang,S.Zhao,and H.Yu,et al.,"Design,analysis and simulation fordevelopment of the first clinical micro-CT scanner."Academic Radiol.12.4(2005):511–525.
-C.Rock,U.Linsenmaier,R.Brandl,et al."Introduction of a new mobile C-arm/CT combination equipment(ISO-C-3D).Initial results of 3-D sectionalimaging."Der Unfallchirurg 104.9(2001):827-833.
-G.Rajiv,M.Grasruck,C.Suess,et al."Ultra-high resolution flat-panelvolume CT:fundamental principles,design architecture,and systemcharacterization."European radiology 16.6(2006):1191-1205.
-F.Noo,R.Clackdoyle,C.Mennessier,et al."Analytic method based onidentification of ellipse parameters for scanner calibration in cone-beamtomography."Phys.Med.Biol.45(2000)3489–3508.
-D.Panetta,N.Belcari,A.D.Guerra,et al."An optimization-based methodfor geometrical calibration in cone-beam CT without dedicated phantoms."Med.Biol.53(2008)3841–3861.
-W.Wein,A.Ladikos,A.Baumgartner."Self-calibration of geometric andradiometric parameters for cone-beam computed tomography."in:Fully 3D,Potsdam,(2011):327–330.
-V.Patel,R.N.Chityala,K.R.Hoffmann,et al."Self-calibration of a cone-beam micro-CT system."Med.Phys.36(2009)48–58.
-B.D.Man,S.Basu,D.Bequé,et al."Multi-source inverse geometry CT:a newsystem concept for X-ray computed tomography."Medical Imaging.InternationalSociety for Optics and Photonics(2007).
-Y.Liu,H.Liu,Y.Wang,et al."Half-scan cone-beam CT fluoroscopy withmultiple sources."Med.Phys.28.7(2001):1466–1471.
-J.Hsieh."Computed Tomography,Principles,Design,Artifacts,and RecentAdvances.Bellingham."WA:SPIE,2003.
-C.R.Crawford and K.F.King."Computed-tomography scanning withsimultaneous patient translation."Med.Phys.17.6(1990):967–982.
-G.Wang and M.W.Vannier."Preliminary study on helical CT algorithmsfor patient motion estimation and compensation."IEEE Trans.Med.Imag.14.2(1995):205–211.
-W.G.Lu and T.R.Mackie."Tomographic motion detection and correctiondirectly in sinogram space."Phys.,Med.,Biol.47.8(2002):1267–1284.
-H.Yu,and G.Wang."Data consistency based rigid motion artifactreduction in fan-beam CT."Medical Imaging,IEEE Transactions on 26.2(2007):249-260.
-H.Yu,Y.Wei,J.Hsieh,et al."Data consistency based translationalmotion artifact reduction in fan-beam CT."Medical Imaging,IEEE Transactionson 25.6(2006):792-803.
-S.Leng,B.Nett,M.Speidel,et al."Motion artifact reduction in fan-beamand cone-beam computed tomography via the fan-beam data consistency condition(FDCC)."Medical Imaging.International Society for Optics and Photonics,2007.
-S.T.Roweis,and S.K.Lawrence."Nonlinear dimensionality reduction bylocally linear embedding."Science 290.5500(2000):2323-2326.
-S.T.Roweis,and S.K.Lawrence."An introduction to locally linearembedding."unpublished.Available at:http://www.cs.toronto.edu/~roweis/lle/publications.html(2000).
-R.A.Horn and C.R.Johnson."Matrix Analysis."Cambridge UniversityPress,Cambridge,1990
-R.L.Siddon."Fast calculation of the exact radiological path for athree-dimensional CT array."Medical physics 12.2(1985):252-255.
-Z.Wang,and A.C.Bovik."A universal image quality index."SignalProcessing Letters,IEEE 9.3(2002):81-84.
-G.Wang,and M.Jiang."Ordered-subset simultaneous algebraicreconstruction techniques(OS-SART)."Journal of X-ray Science and Technology12.3(2004):169-177.
-Y.E.Sidky,C.M.Kao,and X.Pan."Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT."Journal of X-ray Scienceand Technology 14.2(2006):119-139.
Claims (40)
1.一种重建计算机断层(CT)图像的方法,所述方法包括:
i)获得初始CT图像;
ii)对所述初始CT图像执行重建算法以获得所述初始CT图像的重建图像;
iii)使用所述重建图像来调整与所述图像相关联的一个或更多个参数;
iv)使用经调整的所述一个或更多个参数对所述重建图像执行重建算法以获得更新的重建图像;以及
v)使用在每个重复的步骤iii)中最近更新的重建图像来重复步骤iii)至iv),并且在每个重复的步骤iii)中更新相同的一个或更多个参数,直到满足预定特征的阈值。
2.根据权利要求1所述的方法,其中,所述重建算法是基于全变差(TV)的。
3.根据权利要求1至2中任一项所述的方法,其中,使用所述重建图像来调整与所述图像相关联的一个或更多个参数包括对所述一个或更多个参数中的每个参数使用局部线性嵌入(LLE)方法来完成。
4.根据权利要求3所述的方法,其中,所述LLE方法包括:
a)根据欧几里德距离找出每个高维点的K个最近邻;
b)用每个点的近邻线性地表示每个点,并计算其近邻的权重系数;
c)将具有所述权重系数的高维数据映射到本征流形上的低维表示。
5.根据权利要求3至4中任一项所述的方法,其中,所述LLE方法包括:
根据所述欧几里德距离:
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>b</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>,</mo>
</mrow>
在数据向量bi的数据集向量中找到其K个最近邻。
6.根据权利要求3至5中任一项所述的方法,其中,所述LLE方法(还)包括:
使用K个最近向量,用其近邻向量:
<mrow>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>b</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
来线性地表示原始数据向量,
其中是K个最近邻,wik是各自的权重系数,并且能够通过最小化以下误差来求解权重矩阵:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>min</mi>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mi>k</mi>
</munder>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>b</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>j</mi>
</munder>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>.</mo>
</mrow>
7.根据权利要求6所述的方法,其中,使用约束来最小化所述误差等同于求解线性组:
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>j</mi>
</munder>
<msub>
<mi>c</mi>
<mrow>
<mi>j</mi>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
</mrow>
其中C=(cjk)是局部协方差矩阵,其被计算为:
<mrow>
<mi>C</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>b</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>b</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
8.根据权利要求3至7中任一项所述的方法,其中,所述LLE方法(还)包括,利用权重系数W=(wik),通过求解方程:
<mrow>
<munder>
<mi>min</mi>
<mi>Y</mi>
</munder>
<munder>
<mo>&Sigma;</mo>
<mi>i</mi>
</munder>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<munder>
<mo>&Sigma;</mo>
<mi>k</mi>
</munder>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>,</mo>
</mrow>
来计算全局内部坐标Y=(yi)。
9.根据权利要求8所述的方法,其中,能够求解以计算Y=(yi)的方程式的解由广义特征值问题:
MY=λY
的底部的d+1特征向量给出,
其中λ是特征值,并且M=(I-W)T(I-W)。
10.根据权利要求3至9中任一项所述的方法,其中,在所述LLE方法中使用的最近邻的数量K=2。
11.根据权利要求3至9中任一项所述的方法,其中,在所述LLE方法中使用的最近邻的数量至少为2。
12.根据权利要求3至11中任一项所述的方法,其中,所述LLE方法包括求解以下线性方程组:
Au=bi,
其中u=(u1,u2,...,uJ)是表示为J维向量的图像,J是像素数,bi=(b1,b2,...,bL)是数据,L是向量元素的数量,A=(ajk)是与所述一个或更多个(几何)参数相关的投影矩阵。
13.根据权利要求12所述的方法,其中,执行所述重建算法包括计算受到所述一个或更多个(几何)参数影响的投影矩阵:
A=A(P),
其中P是估计参数向量,其包括所述一个或多个更多个参数。
14.根据权利要求13所述的方法,其中,使用距离驱动法计算所述投影矩阵的投影。
15.根据权利要求13至14中任一项所述的方法,其中,使用所述重建图像来调整一个或更多个参数包括:最小化所述投影数据和经重投影的投影数据之间的均方误差,其能够表示为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>P</mi>
<mo>=</mo>
<mi>arg</mi>
<mi> </mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>b</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>|</mo>
<msubsup>
<mo>|</mo>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>P</mi>
<mo>)</mo>
</mrow>
<mi>u</mi>
<mo>=</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中bi是从沿着不同投影视角的测量获得的投影向量,是来自具有采样的参数的重建图像的相应的经重投影的投影向量,P是更新的参数向量。
16.根据权利要求15所述的方法,其中,能够通过在密集采样的参数范围:
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>~</mo>
</mover>
<mi>j</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
内,通过权利要求12和13中提供的方程式来计算所述经重投影的投影向量。
17.根据权利要求16所述的方法,其中,真实参数向量接近相邻的采样的参数向量,并且所测量的投影向量能够由与采样的参数向量相关联的K个最近的经重投影的投影向量来线性地表示,使得:
和
<mrow>
<mi>P</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>w</mi>
<mrow>
<mi>i</mi>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mover>
<mi>P</mi>
<mo>~</mo>
</mover>
<mi>k</mi>
</msub>
<mo>,</mo>
</mrow>
其中是与参数的相应的K个向量相关联的K个最近重投影向量,并且wik是权重系数。
18.根据权利要求1至17中任一项所述的方法,其中,必须满足其阈值才能停止重复步骤iii)至iv)的所述预定特征是通用质量指数(UQI)。
19.根据权利要求18所述的方法,其中,所述阈值至少为0.6。
20.根据权利要求12所述的方法,其中,所述阈值至少为0.7。
21.根据权利要求18所述的方法,其中,所述阈值至少为0.8。
22.根据权利要求18所述的方法,其中,所述阈值为至少0.9。
23.根据权利要求18所述的方法,其中,所述阈值至少为0.95。
24.根据权利要求1至23中任一项所述的方法,其中,所述一个或更多个参数包括以下各项中的至少一个:投影角度、源到对象的距离(SOD)、对象到检测器的距离(ODD)、检测器偏移、检测器倾斜角度、对象x坐标偏移、对象y坐标偏移、和投影角度误差。
25.根据权利要求1-23中任一项所述的方法,其中,所述一个或更多个参数包括以下各项中的至少两个:投影角度、SOD、ODD、检测器偏移、检测器倾斜角度、对象x坐标偏移、对象y坐标偏移、和投影角度误差。
26.根据权利要求1至23中任一项所述的方法,其中,所述一个或更多个参数包括以下各项中的至少三个:投影角度、SOD、ODD、检测器偏移、检测器倾斜角度、对象x坐标偏移、对象y坐标偏移、和投影角度误差。
27.根据权利要求1至23中任一项所述的方法,其中,所述一个或更多个参数包括以下各项中的至少四个:投影角度、SOD、ODD、检测器偏移、检测器倾斜角度、对象x坐标偏移、对象y坐标偏移、和投影角度误差。
28.根据权利要求1至23中任一项所述的方法,其中,所述一个或更多个参数包括:投影角、SOD、ODD、检测器偏移、和检测器倾斜角。
29.根据权利要求1至23中任一项所述的方法,其中,所述一个或更多个参数包括:对象x坐标偏移、对象y坐标偏移、和投影角误差。
30.根据权利要求1至29中任一项所述的方法,其中,获得所述初始CT图像包括:获得包含与所述图像相关联的所述一个或更多个参数的初始值的初始参数向量。
31.根据权利要求30所述的方法,执行所述重建算法包括:基于最近的重建图像,用所述一个或多个参数的更新值来更新所述参数向量。
32.根据权利要求1至31中任一项所述的方法,其中,所述重建算法包括使用有序子集同时代数重建技术(OS-SART)。
33.根据权利要求1至32中任一项所述的方法,其中,所述重建算法包括使用TV进行正则化。
34.根据权利要求1至33中任一项所述的方法,其中,步骤ii)由处理器执行。
35.根据权利要求1至34中任一项所述的方法,其中步骤iii)由处理器执行。
36.根据权利要求1至35中任一项所述的方法,其中步骤iv)由处理器执行。
37.根据权利要求1至36中任一项所述的方法,其中步骤v)由处理器执行。
38.一种计算机可读介质,其具有用于执行根据权利要求1至37中任一项的方法的计算机可执行指令。
39.一种CT***,其包括:
辐射源,
检测器,其用于检测来自所述辐射源的辐射;以及
计算机,其具有根据权利要求38所述的计算机可读介质。
40.根据权利要求39所述的CT***,其中所述辐射源是X射线源。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562135861P | 2015-03-20 | 2015-03-20 | |
US62/135,861 | 2015-03-20 | ||
PCT/US2016/023460 WO2016154136A1 (en) | 2015-03-20 | 2016-03-21 | Automatic system calibration method of x-ray ct |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107592935A true CN107592935A (zh) | 2018-01-16 |
Family
ID=56977665
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201680026497.1A Pending CN107592935A (zh) | 2015-03-20 | 2016-03-21 | X射线ct的自动***校准方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10643354B2 (zh) |
CN (1) | CN107592935A (zh) |
WO (1) | WO2016154136A1 (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108903961A (zh) * | 2018-07-19 | 2018-11-30 | 深圳市倍康美医疗电子商务有限公司 | 一种cbct成像方法、存储介质及*** |
CN110084855A (zh) * | 2019-04-19 | 2019-08-02 | 合肥中科离子医学技术装备有限公司 | 一种改进cbct几何参数标定算法 |
CN110755099A (zh) * | 2019-10-10 | 2020-02-07 | 东软医疗***股份有限公司 | 一种偏转角度的检测方法、校正方法、装置及终端设备 |
CN110974280A (zh) * | 2019-12-23 | 2020-04-10 | 常州博恩中鼎医疗科技有限公司 | 一种适用于锥束ct的几何参数精确自动化校正方法 |
CN111449670A (zh) * | 2020-05-07 | 2020-07-28 | 南京安科医疗科技有限公司 | 一种移动ct***的步进成像方法 |
CN116921269A (zh) * | 2023-09-18 | 2023-10-24 | 无锡日联科技股份有限公司 | 一种半导体x射线3d自动检测***与检测方法 |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10539689B1 (en) | 2015-01-23 | 2020-01-21 | Rensselaer Polytechnic Institute | Spectral X-ray detectors with dynamic electronic control and computational methods |
WO2016154136A1 (en) | 2015-03-20 | 2016-09-29 | Rensselaer Polytechnic Institute | Automatic system calibration method of x-ray ct |
WO2016197127A1 (en) | 2015-06-04 | 2016-12-08 | Rensselaer Polytechnic Institute | Attenuation map reconstruction from tof pet data |
US10410343B2 (en) * | 2015-07-08 | 2019-09-10 | Hitachi, Ltd. | Image computing device, image computing method, and tomograph |
US10729397B2 (en) | 2015-07-20 | 2020-08-04 | Rensselaer Polutechnic Institute | X-ray phase contrast and dark-field information extraction with electric fringe scanning and/or active pixel processing |
WO2017019782A1 (en) | 2015-07-27 | 2017-02-02 | Rensselaer Polytechnic Institute | Combination of an x-ray tube and a source grating with electron beam manipulation |
US10360697B2 (en) * | 2015-09-02 | 2019-07-23 | Siemens Healthcare Gmbh | Fast sparse computed tomography image reconstruction from few views |
US11278250B2 (en) | 2015-11-13 | 2022-03-22 | Rensselaer Polytechnic Institute | Simultaneous interior MRI X-ray imaging system (MRX) |
US10127691B2 (en) * | 2016-11-21 | 2018-11-13 | The Boeing Company | Geometry correction for computed tomography |
US10127658B2 (en) | 2016-11-21 | 2018-11-13 | The Boeing Company | Geometry correction for computed tomography |
WO2019118387A1 (en) * | 2017-12-11 | 2019-06-20 | Dentsply Sirona Inc. | Methods, systems, apparatuses, and computer program products for extending the field of view of a sensor and obtaining a synthetic radiagraph |
KR102591672B1 (ko) | 2018-08-29 | 2023-10-20 | 한국전자통신연구원 | 영상 생성 장치, 영상 생성 장치를 포함하는 이미징 시스템 및 이미징 시스템의 동작 방법 |
CN110349234B (zh) * | 2019-06-27 | 2023-06-20 | 深圳安科高技术股份有限公司 | 一种ct锥形束图像重建方法及其*** |
CN110717959B (zh) * | 2019-10-16 | 2022-07-08 | 首都师范大学 | 基于曲率约束的x射线有限角ct图像重建方法和装置 |
CN112381904B (zh) * | 2020-11-26 | 2024-05-24 | 南京医科大学 | 一种基于DTw-SART-TV迭代过程的有限角度CT图像重建方法 |
JP2022191670A (ja) * | 2021-06-16 | 2022-12-28 | トヨタ自動車株式会社 | X線ct装置及び画像生成方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103514615A (zh) * | 2012-06-22 | 2014-01-15 | 通用电气公司 | 用于迭代重建的方法和设备 |
CN104103086A (zh) * | 2014-06-06 | 2014-10-15 | 华南理工大学 | 一种稀疏采样角度下基于变分不等式的ct图像重建方法 |
Family Cites Families (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6332030B1 (en) * | 1998-01-15 | 2001-12-18 | The Regents Of The University Of California | Method for embedding and extracting digital data in images and video |
WO2005050533A2 (en) * | 2003-11-13 | 2005-06-02 | Honda Motor Co., Ltd. | Image clustering with metric, local linear structure, and affine symmetry |
DE602005014214D1 (zh) | 2004-03-04 | 2009-06-10 | Philips Intellectual Property | |
US7526060B2 (en) * | 2004-03-10 | 2009-04-28 | Koninklijke Philips Electronics N.V. | Artifact correction |
US8005184B2 (en) * | 2007-04-24 | 2011-08-23 | Wisconsin Alumni Research Foundation | Ultra low radiation dose X-ray CT scanner |
TWI382354B (zh) * | 2008-12-02 | 2013-01-11 | Nat Univ Tsing Hua | 臉部辨識方法 |
US8189735B2 (en) | 2010-07-22 | 2012-05-29 | General Electric Company | System and method for reconstruction of X-ray images |
US8885898B2 (en) * | 2010-10-07 | 2014-11-11 | Siemens Medical Solutions Usa, Inc. | Matching of regions of interest across multiple views |
JP5258915B2 (ja) * | 2011-02-28 | 2013-08-07 | 株式会社デンソーアイティーラボラトリ | 特徴変換装置、それを含む類似情報探索装置、コード化パラメータ生成方法、及びコンピュータ・プログラム |
WO2013116709A1 (en) | 2012-02-01 | 2013-08-08 | The Research Foundation of States University of New York | Computerized image reconstruction method and apparatus |
US10285659B2 (en) | 2013-12-06 | 2019-05-14 | Rensselaer Polytechnic Institute | Stored luminescence computed tomography |
US9730657B2 (en) | 2013-12-17 | 2017-08-15 | Rensselaer Polytechnic Institute | Computed tomography based on linear scanning |
US9600910B2 (en) | 2014-01-08 | 2017-03-21 | Rensselaer Polytechnic Institute | Attenuation map reconstruction from TOF PET data |
US20170043041A1 (en) | 2014-04-21 | 2017-02-16 | Rensselaer Polytechnic Institute | Nanoparticle-enabled x-ray magnetic resonance imaging (nxmri) |
US10092257B2 (en) | 2014-10-27 | 2018-10-09 | Rensselaer Polytechnic Institute | Dynamic bowtie filter for cone-beam/multi-slice CT |
US9801600B2 (en) | 2014-11-17 | 2017-10-31 | Rensselaer Polytechnic Institute | X-ray phase-contrast imaging |
US10737111B2 (en) | 2014-12-16 | 2020-08-11 | Rensselaer Polytechnic Institute | X-optogenetics / U-optogenetics |
WO2016106348A1 (en) | 2014-12-22 | 2016-06-30 | Rensselaer Polytechinic Institute | Novel spectral ct systems and methods |
US10539689B1 (en) | 2015-01-23 | 2020-01-21 | Rensselaer Polytechnic Institute | Spectral X-ray detectors with dynamic electronic control and computational methods |
WO2016154136A1 (en) | 2015-03-20 | 2016-09-29 | Rensselaer Polytechnic Institute | Automatic system calibration method of x-ray ct |
WO2016197127A1 (en) | 2015-06-04 | 2016-12-08 | Rensselaer Polytechnic Institute | Attenuation map reconstruction from tof pet data |
US10729397B2 (en) | 2015-07-20 | 2020-08-04 | Rensselaer Polutechnic Institute | X-ray phase contrast and dark-field information extraction with electric fringe scanning and/or active pixel processing |
WO2017019782A1 (en) | 2015-07-27 | 2017-02-02 | Rensselaer Polytechnic Institute | Combination of an x-ray tube and a source grating with electron beam manipulation |
US11020077B2 (en) | 2015-09-14 | 2021-06-01 | Rensselaer Polytechnic Institute | Simultaneous CT-MRI image reconstruction |
US11278250B2 (en) | 2015-11-13 | 2022-03-22 | Rensselaer Polytechnic Institute | Simultaneous interior MRI X-ray imaging system (MRX) |
US11266363B2 (en) | 2016-02-17 | 2022-03-08 | Rensselaer Polytechnic Institute | Energy-sensitive multi-contrast cost-effective CT system |
US11337663B2 (en) | 2016-04-08 | 2022-05-24 | Rensselaer Polytechnic Institute | Rapid filtration methods for dual-energy X-ray CT |
-
2016
- 2016-03-21 WO PCT/US2016/023460 patent/WO2016154136A1/en active Application Filing
- 2016-03-21 CN CN201680026497.1A patent/CN107592935A/zh active Pending
- 2016-03-21 US US15/558,053 patent/US10643354B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103514615A (zh) * | 2012-06-22 | 2014-01-15 | 通用电气公司 | 用于迭代重建的方法和设备 |
CN104103086A (zh) * | 2014-06-06 | 2014-10-15 | 华南理工大学 | 一种稀疏采样角度下基于变分不等式的ct图像重建方法 |
Non-Patent Citations (3)
Title |
---|
GE WANG等: "Ordered-subset simultaneous algebraic reconstruction techniques (OS-SART)", 《JOURNAL OF X-RAY SCIENCE AND TECHNOLOGY》 * |
张俊等: "采用高频能量的CT几何参数自标定方法", 《红外与激光工程》 * |
张兴福等: "局部线性嵌入算法综述", 《2011 INTERNATIONAL CONFERENCE ON FUTURE COMPUTER SCIENCE AND APPLICATION(FCSA 2011)》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108903961A (zh) * | 2018-07-19 | 2018-11-30 | 深圳市倍康美医疗电子商务有限公司 | 一种cbct成像方法、存储介质及*** |
CN110084855A (zh) * | 2019-04-19 | 2019-08-02 | 合肥中科离子医学技术装备有限公司 | 一种改进cbct几何参数标定算法 |
CN110084855B (zh) * | 2019-04-19 | 2020-12-15 | 合肥中科离子医学技术装备有限公司 | 一种改进cbct几何参数标定方法 |
CN110755099A (zh) * | 2019-10-10 | 2020-02-07 | 东软医疗***股份有限公司 | 一种偏转角度的检测方法、校正方法、装置及终端设备 |
CN110755099B (zh) * | 2019-10-10 | 2023-10-31 | 东软医疗***股份有限公司 | 一种偏转角度的检测方法、校正方法、装置及终端设备 |
CN110974280A (zh) * | 2019-12-23 | 2020-04-10 | 常州博恩中鼎医疗科技有限公司 | 一种适用于锥束ct的几何参数精确自动化校正方法 |
CN111449670A (zh) * | 2020-05-07 | 2020-07-28 | 南京安科医疗科技有限公司 | 一种移动ct***的步进成像方法 |
CN111449670B (zh) * | 2020-05-07 | 2022-12-30 | 南京安科医疗科技有限公司 | 一种移动ct***的步进成像方法 |
CN116921269A (zh) * | 2023-09-18 | 2023-10-24 | 无锡日联科技股份有限公司 | 一种半导体x射线3d自动检测***与检测方法 |
CN116921269B (zh) * | 2023-09-18 | 2024-01-23 | 无锡日联科技股份有限公司 | 一种半导体x射线3d自动检测***与检测方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2016154136A1 (en) | 2016-09-29 |
US20180068467A1 (en) | 2018-03-08 |
US10643354B2 (en) | 2020-05-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107592935A (zh) | X射线ct的自动***校准方法 | |
CN104107065B (zh) | 3d图像集在不同空间之间的最佳变换 | |
CN102346924B (zh) | 用于x射线图像的重建的***和方法 | |
Bieniosek et al. | Characterization of custom 3D printed multimodality imaging phantoms | |
CN1957847B (zh) | 再现对象的断层造影图像的方法和断层造影设备 | |
CN103153192B (zh) | X射线ct装置以及图像再构成方法 | |
Aguiar et al. | Geometrical and Monte Carlo projectors in 3D PET reconstruction | |
Qi et al. | Extraction of tumor motion trajectories using PICCS‐4DCBCT: a validation study | |
Grangeat | Tomography | |
CN114067013A (zh) | 用于经由单应重采样变换进行重投影和反投影的***和方法 | |
Chen et al. | X-ray CT geometrical calibration via locally linear embedding | |
Chen et al. | Development and validation of an open data format for CT projection data | |
Chang et al. | A geometric calibration method for the digital chest tomosynthesis with dual-axis scanning geometry | |
US9495770B2 (en) | Practical model based CT construction | |
Bajpai et al. | High resolution 3d image reconstruction using the algebraic method for cone-beam geometry over circular and helical trajectories | |
US20070230762A1 (en) | Fast iterative 3D pet image reconstruction using a set of 2D linogram transformations | |
Gao et al. | Principal component reconstruction (PCR) for cine CBCT with motion learning from 2D fluoroscopy | |
Nakano et al. | Cone-beam CT reconstruction for non-periodic organ motion using time-ordered chain graph model | |
Miao | Comparative studies of different system models for iterative CT image reconstruction | |
Chen et al. | A novel extension of the parallel‐beam projection‐slice theorem to divergent fan‐beam and cone‐beam projections | |
Grimmer et al. | Cone‐beam CT image reconstruction with extended range | |
Zhong et al. | A dual‐view digital tomosynthesis imaging technique for improved chest imaging | |
Rodríguez-Alvarez et al. | Expectation maximization (EM) algorithms using polar symmetries for computed tomography (CT) image reconstruction | |
US9508164B2 (en) | Fast iterative image reconstruction method for 3D computed tomography | |
CN107251095A (zh) | 图像重建***、方法和计算机程序 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180116 |
|
WD01 | Invention patent application deemed withdrawn after publication |