CN106725567A - 电子计算机x射线断层扫描仪 - Google Patents
电子计算机x射线断层扫描仪 Download PDFInfo
- Publication number
- CN106725567A CN106725567A CN201611243857.1A CN201611243857A CN106725567A CN 106725567 A CN106725567 A CN 106725567A CN 201611243857 A CN201611243857 A CN 201611243857A CN 106725567 A CN106725567 A CN 106725567A
- Authority
- CN
- China
- Prior art keywords
- ray
- detector
- intensity
- instrument
- scanner
- 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
- 230000005855 radiation Effects 0.000 claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000008569 process Effects 0.000 claims abstract description 10
- 230000003068 static effect Effects 0.000 claims abstract description 6
- 230000000007 visual effect Effects 0.000 claims description 13
- 230000007246 mechanism Effects 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 3
- 230000003287 optical effect Effects 0.000 claims description 3
- ZOKXTWBITQBERF-UHFFFAOYSA-N Molybdenum Chemical compound [Mo] ZOKXTWBITQBERF-UHFFFAOYSA-N 0.000 claims description 2
- 229910052750 molybdenum Inorganic materials 0.000 claims description 2
- 239000011733 molybdenum Substances 0.000 claims description 2
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 claims description 2
- 229910052721 tungsten Inorganic materials 0.000 claims description 2
- 239000010937 tungsten Substances 0.000 claims description 2
- 230000000873 masking effect Effects 0.000 claims 1
- 229910052751 metal Inorganic materials 0.000 claims 1
- 239000002184 metal Substances 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 4
- 238000005728 strengthening Methods 0.000 abstract 1
- 238000000691 measurement method Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000000053 physical method Methods 0.000 description 2
- 238000004846 x-ray emission Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 239000011133 lead Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000002250 progressing effect Effects 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4021—Arrangements for generating radiation specially adapted for radiation diagnosis involving movement of the focal spot
-
- 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/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4291—Arrangements for detecting radiation specially adapted for radiation diagnosis the detector being combined with a grid or grating
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Pulmonology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供了一种电子计算机X射线断层扫描仪。本发明涉及一种CT扫描仪,目的是提供一种实际散焦强度的测量和在此基础上的散焦校正。扫描仪主控台被配置为执行一旋转扫描曝光,以及以预定的采样率从检测器阵列的至少部分检测器中获取辐射强度值;其中在旋转扫描曝光过程中,X射线球管及检测器阵列围绕扫描仪的机架旋转中心旋转,而在旋转的每一圈中,处于静止状态的模体会逐渐先遮蔽后退出检测器阵列中各检测器的X射线最大辐射区,从而使各检测器所接收的X射线辐射强度相应地先减弱后增强;扫描仪的建像机被配置为依据模体逐渐进入并退出至少检测器的X射线最大辐射区期间相邻采样点的检测器的辐射强度值变化,分别计算对应位置的X射线球管散焦强度,得到X射线球管的散焦强度分布。
Description
本申请是申请日为2012年11月05日、申请号为“201210437063.4”、发明名称为“CT扫描仪及其散焦强度测量方法及散焦校正方法”的发明专利申请的分案申请。
技术领域
本发明涉及电子计算机X射线断层扫描仪(简称CT扫描仪),尤其是涉及CT扫描仪的散焦强度测量方法及散焦校正方法。
背景技术
CT扫描仪是利用计算机技术对被测物体断层扫描图像进行重建获得三维断层图像的设备。CT扫描仪使用X射线球管发射X射线到焦点区,从而聚焦到被照射物体(如人体器官)。散焦辐射是X射线球管内从焦点以外区域辐射出X射线的现象,它是由于二次电子和场发射电子轰击球管靶区焦点以外区域产生X射线造成的。散焦现象会导致被照射物体边缘部分对比度退化或产生阴影,影响甚至误导医生根据图像进行诊断。因此需要在CT扫描仪包含对散焦现象的校正。
美国专利US 6628744 B1提出了一种CT扫描仪中散焦辐射的校正方法。该方法使用依据理论模型计算的散焦校正系数在数据域进行散焦校正。理论模型计算的缺陷在于,它可能与实际产品的散焦辐射情况有出入,从而导致所计算的散焦校正系数不准确。
为了解决上述问题,本发明提供一种CT扫描仪的散焦强度测量方法以及根据该测量方法获得的数据进行散焦校正的方法。
发明内容
本发明所要解决的技术问题是提供一种CT扫描仪的散焦强度测量方法及散焦校正方法,它可以实际测量每一CT扫描仪的散焦强度,并作为散焦校正的依据。
本发明为解决上述技术问题而采用的技术方案是提出一种电子计算机X射线断层扫描仪的散焦强度的测量方法,包括以下步骤:将一能够遮蔽X射线的模体置于该扫描仪的机架孔径内;使该扫描仪的X射线球管及检测器阵列围绕扫描仪的机架旋转中心旋转以执行旋转曝光扫描,在旋转的每一圈中,处于静止状态的该模体会逐渐先遮蔽后退出该检测器阵列中各检测器的X射线最大辐射区,从而使各检测器所接收的X射线辐射强度相应地先减弱后增强;以预定的采样率从该检测器阵列的至少部分检测器中获取辐射强度值;以及依据该模体逐渐进入并退出该至少部分检测器的X射线最大辐射区期间相邻采样点的检测器的辐射强度值变化,分别计算对应位置的X射线球管散焦强度,从而得到该X射线球管的散焦强度分布。
在本发明的一实施例中,该模体垂直于该扫描仪的机架光平面。
在本发明的一实施例中,上述方法包括将该散焦强度分布作为该扫描仪出厂前的预存设置。
在本发明的一实施例中,上述方法包括在该扫描仪出厂后定期地更新该散焦强度分布。
本发明还提出一种电子计算机X射线断层扫描仪的散焦校正方法,包括以下步骤:提供该扫描仪的X射线球管的散焦强度分布;获取未进行散焦校正的基础图像;对基础图像进行正投影,得到原始投影值;将原始投影值转化为原始强度值;根据该原始强度值以及X射线球管的散焦强度分布,计算散焦引起的误差强度;根据该误差强度计算误差投影;将误差投影进行图像重建,得到误差图像;以及从该基础图像中减去误差图像,得到最终校正后的图像。
在本发明的一实施例中,对基础图像进行正投影,得到原始投影值之前还包括:判断该基础图像是否包含它所对应断面的整个被扫描物体,若是,则执行该正投影步骤,否则以更大的视野重建该基础图像。
在本发明的一实施例中,所述散焦校正是在该扫描仪的建像机中执行。
在本发明的一实施例中,所述误差强度通过如下方式获得:其中OffR(i)是该散焦强度分布和一焦点强度的比例系数,N是对散焦强度测量时的采样点数,原始投影(i)是对应于第i条散焦射线经过的被扫描物体的投影值。
本发明还提出一种电子计算机X射线断层扫描仪,包括一机架、一能够遮蔽X射线的模体、主控台以及一建像机。该机架具有一孔径,该机架内部具有一旋转机构,该旋转机构包含设于该孔径一侧的X射线球管和设于该孔径另一侧的检测器阵列。该模体适于放置于该扫描仪的机架孔径内。该主控台被配置为执行一旋转扫描曝光,以及以预定的采样率从该检测器阵列的至少部分检测器中获取辐射强度值;其中在旋转扫描曝光过程中,该X射线球管及该检测器阵列围绕扫描仪的机架旋转中心旋转,而在旋转的每一圈中,处于静止状态的该模体会逐渐先遮蔽后退出该检测器阵列中各检测器的X射线最大辐射区,从而使各检测器所接收的X射线辐射强度相应地先减弱后增强。该建像机,被配置为依据该模体逐渐进入并退出至少检测器的X射线最大辐射区期间相邻采样点的检测器的辐射强度值变化,分别计算对应位置的X射线球管散焦强度,从而得到该X射线球管的散焦强度分布。
在本发明的一实施例中,在上述建像机中执行如下步骤:获取未进行散焦校正的基础图像;对基础图像进行正投影,得到原始投影值;将原始投影值转化为原始强度值;根据该原始强度值以及X射线球管的散焦强度分布,计算散焦引起的误差强度;根据该误差强度计算误差投影;将误差投影进行图像重建,得到误差图像;以及从该基础图像中减去误差图像,得到最终校正后的图像。
本发明由于采用以上技术方案,使之与现有技术相比,具有如下显著优点:
1、本发明测量散焦强度的方式是基于实测的方法,相比理论模型的计算方法,具有准确、操作方便、结果实用性强、适应性强的特点。
2、本发明根据测量的结果在图像域进行散焦校正,对于薄切片扫描重建厚图像的情况,基于图像的校正方式可以减少运算量,节约校正时间;另外,由于默认用户仅对所选视野内组织感兴趣,仅针对图像视野内区域进行散焦校正,降低参与运算的人体体积,也即降低数据量,节约校正时间。
附图说明
为让本发明的上述目的、特征和优点能更明显易懂,以下结合附图对本发明的具体实施方式作详细说明,其中:
图1示出本发明的CT扫描仪成像系统的示意图。
图2示出本发明物理测量散焦强度的示意图,其中模体未遮挡X射线球管的光子辐射区。
图3示出本发明的物理测量散焦强度的示意图,其中模体逐步遮挡X射线球管的光子辐射区。
图4示出本发明散焦强度测量方法的流程示意图。
图5示出根据本发明散焦强度测量方法测得的散焦强度分布进行散焦校正的建像流程示意图。
图6示出本发明一实施例的散焦校正示意图。
具体实施方式
图1示出本发明一实施例的CT扫描仪成像系统的示意图。参照图1所示,CT扫描仪100包括机架110,该机架110包括一旋转机构,该旋转机构具有一孔径111。在孔径111的一侧设有X射线球管112。X射线球管112产生的X射线主要地由焦点O射出,然后射向位于孔径111内的被照射物体(如人体)。孔径111的另一侧设有检测器阵列114,用于检测穿过被照射物体后的X射线强度。当X射线球管112和检测器阵列114设置在旋转机构上,当该旋转机构旋转时,通过X射线球管112连续地照射及检测器阵列114连续地检测,就能得到被照射物体各个角度的辐射强度。
在本实施例中,CT扫描仪100的孔径111内设置有一模体113,用于部分地遮蔽X射线球管112射出的X射线。为达到此目的,模体113使用可屏蔽X射线的材料,例如钼、钨、或铅等强衰减物质。在本实施例中,模体113为边缘平滑、厚度均匀的金属板。
在实施散焦强度测量时,模体113置于孔径111内,其垂直于机架110的光平面,并在机架轴向方向上能够覆盖所有X光视野。在本发明的实施例中,可以使用机架外罩或其他装置来支撑模体113,使其稳定地固定,并不随旋转机构旋转。在正式使用时,模体113可自机架上取下。
检测器阵列114所获得的数据将通过数据传输链路120传输到建像机130。在建像机130内将根据检测器阵列114所获得的数据完成影像的重建。所重建的影像可在图形显示设备150中显示。主控台140用于对CT扫描仪100的控制。例如,当执行旋转曝光扫描时,主控台140可控制旋转机构旋转,并以一定的采样率从检测器阵列114获得输出辐射强度数据。
图2和图3示出本发明散焦校正系数的测量示意图。参照图2,未示出的X射线球管112发射的X射线可以视为被很大范围低强度X光发射源(散焦区)环绕的高强度X光发射源(焦点区)。X射线的辐射区以标记201表示,该辐射区201包括焦点区域O和散焦区OFF。
在图2中示出检测器阵列114中的单个检测器114a,它能接收的X射线最大辐射区为W。
在旋转曝光扫描时,X射线球管112发射X射线,检测器114a检测其最大辐射区W内的辐射,可以获得一圈多个角度下的扇束投影,对每个角度下的投影称为一个检测视野(view)。这个过程可以等效视为机架110的旋转机构静止,而模体113沿孔径111绕机架110的旋转中心旋转一周,旋转方向为A,而模体113未遮挡X射线球管112的X光子区域时模体113相对于机架110的位置为P。
对每个检测器114a而言,当模体113未‘切割’该检测器114a所能接收的最大辐射区时(如图2的位置P),该检测器114a接收到全部最大辐射区W的能量;而当模体113进入最大辐射区W,模体113会逐步遮挡(或相反过程,释放)该检测器114a的最大辐射区W。
例如图3所示,在旋转曝光扫描过程中的时刻i-1,模体113处在位置Pi-1,此时假设模体113未‘切割’该检测器114a所能接收的球管X射线辐射区时,该检测器114a接收到全部最大辐射区W的能量,其辐射强度以Si-1表示。在时刻i,模体113刚进入最大辐射区W,遮挡一部分的该检测器114a的最大辐射区W,检测器114a只能接收到最大辐射区W的其余部分的能量,其辐射强度以Si表示。辐射强度Si-1和Si之间的差异在图3中直观地示出。
模体113离开最大辐射区W的过程与图3所示恰好相反,在此不再描述。
从上述实例可以看出,当模体113‘切割’最大辐射区W时,可以视为,相邻扇形束投影的强度差异(如图中第i-1和第i个检测视野的强度差异),就是该检测器在该几何角度下球管散焦强度值(ΔS)。
对每个检测器,在一圈的多个角度扇束投影中,都可以找到一组连续的检测视野,检测器输出强度值在这一组检测视野之间的差异,即可视为将球管X射线辐射区离散化的强度分布。提供足够的采样率,可使该强度分布达到所需的精度。
这样,对所有检测器通道进行此类运算,就可以获得所有检测器通道所接收的球管X射线辐射区的强度分布(例如图2所示的X射线的辐射区201所示出的强度分布)。
基于上述描述,归纳本发明一实施例的散焦强度测量方法流程如图4所示,其过程描述如下:
在步骤401,将一能够遮蔽X射线的模体置于CT扫描仪的机架上;
在步骤402,令CT扫描仪执行旋转曝光扫描,使X射线球管及检测器阵列围绕CT扫描仪的机架旋转中心旋转。在旋转的每一圈中,处于静止状态的模体会逐渐先遮蔽后退出检测器阵列中各检测器的X射线最大辐射区,从而使各检测器所接收的X射线辐射强度相应地先减弱后增强;
在步骤403,以预定的采样率从检测器阵列的各检测器中获取辐射强度值;
在步骤404,依据模体逐渐进入并退出各检测器的X射线最大辐射区期间,相邻采样点的检测器的辐射强度值变化,计算对应位置的X射线球管散焦强度。所有位置的X射线球管散焦强度组合成X射线球管的散焦强度分布。
上述的散焦强度测量方法可以在各个CT扫描仪组装完成后、出厂前实施,所获得的散焦强度分布可以预存于CT扫描仪中,以便在使用时进行校正。另外,考虑到X射线球管在使用过程中的变化,上述的散焦强度测量方法也可以在各个CT扫描仪的使用期间定期地实施,以获得更新的散焦强度分布。
在获得上述散焦强度分布图后,可执行基于图像的散焦校正,这一过程可在CT扫描仪100的建像机130中执行,其方法流程如下:
步骤501,提供X射线球管的散焦强度分布。它可以按照上述实施例的方法获得,并预存于CT扫描仪中。
步骤502,CT扫描仪获取未进行散焦校正的图像,作为基础图像521。例如,CT扫描仪按照正常的流程对进入机架孔径中被扫描物体进行旋转曝光扫描,依据所获得的X射线强度重建有关该被扫描物体的基础图像。
基础图像需要包含它所对应断面的整个被扫描物体,或至少大于被重建视野一定范围,以保障获得用于校正的所有正投影数据。因此在步骤503,判断重建视野是否包含断面的整个被扫描物体,若是,则继续进入到步骤504,否则在步骤505以更大的视野重建。扩大视野的程度取决于步骤501执行过程所获得的球管X射线辐射区的宽度。
在步骤506,对基础图像521进行正投影,得到原始投影值523。
在步骤507,将原始投影值523转化为原始强度值。
在步骤508,根据对原始强度采样值,以及步骤501得到的球管散焦强度的分布,计算散焦对检测器的接收强度的影响(以误差强度表征)。
在步骤509,根据该误差强度计算出投影域的误差,称为误差投影。
误差投影=-log(原始强度-误差强度)-原始投影。
在步骤510,将误差投影进行图像重建,得到误差图像524。
在步骤511,从原始未校正的基础图像521中减去误差图像524,得到最终校正后的图像525。
在上述的步骤508中,
其中OffR(i)是在步骤501中得到的散焦强度分布和焦点强度的比例系数,N是对散焦强度测量时的采样点数。原始投影(i)是对应于第i条散焦射线经过的被扫描物体的投影值。
以图6为例,以下描述原始投影的计算方式:
设CT扫描仪机架的旋转中心距离某原始投影路径的距离是D,则该投影所属通道编号是:
通道号=中心通道号+arcsin(D/焦心距)/ΔchR
其中中心通道号是焦点经过旋转中心到达检测器通道的位置编号,焦心距是焦点到达旋转中心的距离,ΔchR是每个通道对应的扇角。
原始投影所属检测视野(view)和当前view的距离是:
其中α是原始投影位置和“焦点——旋转中心”连线的夹角,ΔviewR是每个view采样划过的角度。
所以原始投影(i)中的每一项均可从相邻view某通道的正投影值获得。
在本发明的实施例中,π+扇形束角度的正投影范围获得的数据即可用于校正。
在本发明的实施例中,由于散焦本身的低频特性,可以采用较低的通道数和检测视野数进行正投影。
在本发明的实施例中,根据投影值域值判断,若投影区域全是空气,则不进行误差计算。
虽然本发明已以较佳实施例揭示如上,然其并非用以限定本发明,任何本领域技术人员,在不脱离本发明的精神和范围内,当可作些许的修改和完善,因此本发明的保护范围当以权利要求书所界定的为准。
Claims (9)
1.一种电子计算机X射线断层扫描仪,包括:
一机架,具有一孔径,该机架内部具有一旋转机构,该旋转机构包含设于该孔径一侧的X射线球管和设于该孔径另一侧的检测器阵列;
一能够遮蔽X射线的模体,适于放置于该扫描仪的机架孔径内,其中所述模体为金属板;
一主控台,被配置为执行一旋转扫描曝光,以及以预定的采样率从该检测器阵列的至少部分检测器中获取辐射强度值;其中在旋转扫描曝光过程中,该X射线球管及该检测器阵列围绕扫描仪的机架旋转中心旋转,而在旋转的每一圈中,处于静止状态的该模体会逐渐先遮蔽后退出该检测器阵列中各检测器的X射线最大辐射区,从而使各检测器所接收的X射线辐射强度相应地先减弱后增强;以及
一建像机,被配置为依据该模体逐渐进入并退出至少检测器的X射线最大辐射区期间相邻采样点的检测器的辐射强度值变化,分别计算对应位置的X射线球管散焦强度,从而得到该X射线球管的散焦强度分布。
2.如权利要求1所述的电子计算机X射线断层扫描仪,其特征在于:该模体垂直于该扫描仪的机架光平面。
3.如权利要求1所述的电子计算机X射线断层扫描仪,其特征在于,该扫描仪将该散焦强度分布作为出厂前的预存设置。
4.如权利要求1所述的电子计算机X射线断层扫描仪,其特征在于,该扫描仪出厂后定期地更新该散焦强度分布。
5.如权利要求1所述的电子计算机X射线断层扫描仪,其特征在于,所述金属板由钼、钨或铅中至少一种物质组成。
6.如权利要求1所述的电子计算机X射线断层扫描仪,其特征在于,所述模体由所述扫描仪的机架外罩支撑。
7.如权利要求1所述的电子计算机X射线断层扫描仪,其特征在于:在所述建像机中执行如下步骤:
获取未进行散焦校正的基础图像;
对基础图像进行正投影,得到原始投影值;
将原始投影值转化为原始强度值;
根据该原始强度值以及X射线球管的散焦强度分布,计算散焦引起的误差强度;
根据该误差强度计算误差投影;
将误差投影进行图像重建,得到误差图像;以及
从该基础图像中减去误差图像,得到最终校正后的图像。
8.如权利要求7所述的电子计算机X射线断层扫描仪,其特征在于,对基础图像进行正投影,得到原始投影值之前还包括:判断该基础图像是否包含它所对应断面的整个被扫描物体,若是,则执行该正投影步骤,否则以更大的视野重建该基础图像。
9.如权利要求7所述的电子计算机X射线断层扫描仪,其特征在于,所述误差强度通过如下方式获得:其中OffR(i)是该散焦强度分布和一焦点强度的比例系数,N是对散焦强度测量时的采样点数,原始投影(i)是对应于第i条散焦射线经过的被扫描物体的投影值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611243857.1A CN106725567B (zh) | 2012-11-05 | 2012-11-05 | 电子计算机x射线断层扫描仪 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210437063.4A CN103800025B (zh) | 2012-11-05 | 2012-11-05 | Ct扫描仪及其散焦强度测量方法及散焦校正方法 |
CN201611243857.1A CN106725567B (zh) | 2012-11-05 | 2012-11-05 | 电子计算机x射线断层扫描仪 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210437063.4A Division CN103800025B (zh) | 2012-11-05 | 2012-11-05 | Ct扫描仪及其散焦强度测量方法及散焦校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106725567A true CN106725567A (zh) | 2017-05-31 |
CN106725567B CN106725567B (zh) | 2020-02-07 |
Family
ID=50697759
Family Applications (4)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210437063.4A Active CN103800025B (zh) | 2012-11-05 | 2012-11-05 | Ct扫描仪及其散焦强度测量方法及散焦校正方法 |
CN201611243499.4A Active CN106667515B (zh) | 2012-11-05 | 2012-11-05 | Ct扫描仪散焦校正方法 |
CN201611244684.5A Active CN106725568B (zh) | 2012-11-05 | 2012-11-05 | Ct扫描仪散焦强度测量方法 |
CN201611243857.1A Active CN106725567B (zh) | 2012-11-05 | 2012-11-05 | 电子计算机x射线断层扫描仪 |
Family Applications Before (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210437063.4A Active CN103800025B (zh) | 2012-11-05 | 2012-11-05 | Ct扫描仪及其散焦强度测量方法及散焦校正方法 |
CN201611243499.4A Active CN106667515B (zh) | 2012-11-05 | 2012-11-05 | Ct扫描仪散焦校正方法 |
CN201611244684.5A Active CN106725568B (zh) | 2012-11-05 | 2012-11-05 | Ct扫描仪散焦强度测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (4) | CN103800025B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110559001A (zh) * | 2019-09-03 | 2019-12-13 | 明峰医疗系统股份有限公司 | 一种ct扫描仪散焦辐射强度分布测量方法 |
CN110755101A (zh) * | 2019-10-23 | 2020-02-07 | 东软医疗系统股份有限公司 | 一种ct扫描方法和装置 |
CN117064427A (zh) * | 2023-10-07 | 2023-11-17 | 北京富通康影科技有限公司 | 一种三源ct增强扫描中造影剂浓度监控及跟踪方法与系统 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106361367B (zh) * | 2016-12-01 | 2019-10-08 | 上海联影医疗科技有限公司 | 一种检测器的校正方法和使用该校正方法的装置及设备 |
CN108169787B (zh) * | 2017-12-27 | 2020-10-13 | 上海联影医疗科技有限公司 | Pet探测器、pet系统及pet探测方法 |
JP6988732B2 (ja) * | 2018-08-02 | 2022-01-05 | 株式会社島津製作所 | 放射線画像処理装置および放射線画像処理方法 |
CN110108230B (zh) * | 2019-05-06 | 2021-04-16 | 南京理工大学 | 基于图像差分与lm迭代的二值光栅投影离焦程度评估方法 |
CN111789624B (zh) * | 2020-06-29 | 2023-08-08 | 上海联影医疗科技股份有限公司 | 散焦辐射测量方法、装置、计算机设备和可读存储介质 |
CN113017673B (zh) * | 2021-03-03 | 2022-08-26 | 明峰医疗系统股份有限公司 | 一种模体的定位装置 |
CN114903507B (zh) * | 2022-05-16 | 2023-06-09 | 张海光 | 一种医学图像数据处理系统及方法 |
CN115105108B (zh) * | 2022-06-30 | 2023-06-09 | 赛诺威盛科技(北京)股份有限公司 | 散焦校正方法、装置、散焦校正模体及电子装置 |
CN117679060B (zh) * | 2024-02-01 | 2024-06-21 | 赛诺威盛科技(北京)股份有限公司 | 基于薄圆盘模体的球管机械位置校正方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5757951A (en) * | 1994-12-20 | 1998-05-26 | Picker International, Inc. | Correction of off-focal radiation |
US6628744B1 (en) * | 2002-09-26 | 2003-09-30 | Koninklijke Philips Electronics N.V. | Off-focal radiation correction in CT |
US20100246918A1 (en) * | 2009-03-26 | 2010-09-30 | Steffen Kappler | Iterative extra-focal radiation correction in the reconstruction of ct images |
CN102138804A (zh) * | 2009-10-01 | 2011-08-03 | 西门子公司 | 在计算开支最小化条件下降低焦外辐射效应的ct数据处理 |
CN102727234A (zh) * | 2011-03-31 | 2012-10-17 | 西门子公司 | 产生检查对象的图像数据的方法、处理装置和x射线系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4809172A (en) * | 1987-03-11 | 1989-02-28 | Kabushiki Kaisha Toshiba | Method for obtaining image data with a tomographic apparatus |
JP2006102299A (ja) * | 2004-10-07 | 2006-04-20 | Ge Medical Systems Global Technology Co Llc | X線線量補正方法およびx線ct装置 |
-
2012
- 2012-11-05 CN CN201210437063.4A patent/CN103800025B/zh active Active
- 2012-11-05 CN CN201611243499.4A patent/CN106667515B/zh active Active
- 2012-11-05 CN CN201611244684.5A patent/CN106725568B/zh active Active
- 2012-11-05 CN CN201611243857.1A patent/CN106725567B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5757951A (en) * | 1994-12-20 | 1998-05-26 | Picker International, Inc. | Correction of off-focal radiation |
US6628744B1 (en) * | 2002-09-26 | 2003-09-30 | Koninklijke Philips Electronics N.V. | Off-focal radiation correction in CT |
US20100246918A1 (en) * | 2009-03-26 | 2010-09-30 | Steffen Kappler | Iterative extra-focal radiation correction in the reconstruction of ct images |
CN102138804A (zh) * | 2009-10-01 | 2011-08-03 | 西门子公司 | 在计算开支最小化条件下降低焦外辐射效应的ct数据处理 |
CN102727234A (zh) * | 2011-03-31 | 2012-10-17 | 西门子公司 | 产生检查对象的图像数据的方法、处理装置和x射线系统 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110559001A (zh) * | 2019-09-03 | 2019-12-13 | 明峰医疗系统股份有限公司 | 一种ct扫描仪散焦辐射强度分布测量方法 |
CN110755101A (zh) * | 2019-10-23 | 2020-02-07 | 东软医疗系统股份有限公司 | 一种ct扫描方法和装置 |
CN110755101B (zh) * | 2019-10-23 | 2023-07-28 | 东软医疗系统股份有限公司 | 一种ct扫描方法和装置 |
CN117064427A (zh) * | 2023-10-07 | 2023-11-17 | 北京富通康影科技有限公司 | 一种三源ct增强扫描中造影剂浓度监控及跟踪方法与系统 |
CN117064427B (zh) * | 2023-10-07 | 2024-05-14 | 北京富通康影科技有限公司 | 一种三源ct增强扫描中造影剂浓度监控及跟踪方法与系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106667515A (zh) | 2017-05-17 |
CN106725567B (zh) | 2020-02-07 |
CN106725568A (zh) | 2017-05-31 |
CN103800025A (zh) | 2014-05-21 |
CN103800025B (zh) | 2016-12-21 |
CN106725568B (zh) | 2019-12-20 |
CN106667515B (zh) | 2019-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103800025B (zh) | Ct扫描仪及其散焦强度测量方法及散焦校正方法 | |
US10022564B2 (en) | Systems, methods, and devices for radiation beam alignment and radiation beam measurements using electronic portal imaging devices | |
US20070274457A1 (en) | Method and apparatus to control radiation tube focal spot size | |
US6980623B2 (en) | Method and apparatus for z-axis tracking and collimation | |
US6370218B1 (en) | Methods and systems for determining x-ray beam position in multi-slice computed tomography scanners | |
US8077826B2 (en) | CT scanner with scatter radiation correction and method of using same | |
CN110559001B (zh) | 一种ct扫描仪散焦辐射强度分布测量方法 | |
JP2000516834A (ja) | 調整可能なコンピュータ断層撮影装置 | |
EP2638558A1 (en) | Determining changes in the x-ray emission yield of an x-ray source | |
CN106687045B (zh) | 数据处理装置、x射线ct装置以及参考校正方法 | |
US7372936B2 (en) | Radiation computed tomography apparatus and tomographic image data generating method | |
US5579359A (en) | Methods and apparatus for calibrating detector cell output signals | |
US10542947B2 (en) | Photon-counting CT apparatus | |
US20040022364A1 (en) | Method for determining correction coefficients for detector channels for a computed tomograph | |
JPH08510090A (ja) | 診断用x線管に印加される電圧の間接測定 | |
EP0982000B1 (en) | Apparatus for dose verification in an imaging system | |
US20150190108A1 (en) | 3D Image Generation Method and Device for G-arm X-ray Machine and G-arm X-ray Machine | |
JP5931394B2 (ja) | X線診断装置及び線量分布データ生成方法 | |
US4260895A (en) | Radiation diagnostic apparatus for generating tomographic images | |
Moore et al. | Cone beam CT with zonal filters for simultaneous dose reduction, improved target contrast and automated set-up in radiotherapy | |
US20050226366A1 (en) | Method, system and storage medium for reference normalization for blocked reference channels | |
CN106955118B (zh) | Ct形状过滤器缺陷的单能矫正方法 | |
US7101078B1 (en) | Methods and systems for imaging system radiation source alignment | |
JP6777556B2 (ja) | X線ct装置 | |
JPS62284250A (ja) | 産業用ctスキヤナ |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |