CN105427355A - 罩壳、x射线图像的散射成分计算、重建的方法及装置 - Google Patents
罩壳、x射线图像的散射成分计算、重建的方法及装置 Download PDFInfo
- Publication number
- CN105427355A CN105427355A CN201510834741.4A CN201510834741A CN105427355A CN 105427355 A CN105427355 A CN 105427355A CN 201510834741 A CN201510834741 A CN 201510834741A CN 105427355 A CN105427355 A CN 105427355A
- Authority
- CN
- China
- Prior art keywords
- image
- ray
- modulation
- ray image
- case
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000000203 mixture Substances 0.000 claims description 60
- 238000003384 imaging method Methods 0.000 claims description 24
- 238000001914 filtration Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000000903 blocking effect Effects 0.000 claims description 7
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 230000000737 periodic effect Effects 0.000 claims description 4
- 238000001959 radiotherapy Methods 0.000 description 13
- 238000011282 treatment Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 239000000463 material Substances 0.000 description 4
- 230000007423 decrease Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000000729 antidote Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000007408 cone-beam computed tomography Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000002786 image-guided radiation therapy Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- G—PHYSICS
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Optics & Photonics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种罩壳，以及X射线图形的散射成分计算、X射线图像重建的方法和装置。罩壳围绕床板，设置在床板的外部，包括：沿第一方向延伸的第一部件；所述第一部件在与所述第一方向垂直的第二方向上重复排列；所述第一部件在沿第一方向延伸的两侧边缘相互平行；所述第一部件对X射线产生遮挡。本发明的技术方案实现对X射线图像中散射成分的计算，有效去除了X射线图像的散射成分，提高了图像质量。并且罩壳的固定方式稳定，不会因为振动引起图像的运动模糊。
Description
技术领域
本发明涉及医学影像领域，特别涉及一种形成调制X射线图像的罩壳、X射线图像散射成分、X射线图像重建的方法及装置。
背景技术
在X射线图像中，散射线一直都是导致图像质量不理想的一个重要因素。特别是对于有较大探测器面积的X射线***来说，例如放射治疗设备，用到的锥束CT(conebeamcomputedtomography)采用的是电子射野影像装置(EPID，ElectronicPortalImagingDevice)作为接收X射线的探测器，在其接收到的X射线中散射线占了很大的比例。散射线在图像上引起杯状伪影，降低了图像的对比度，极大降低了重建图像的图像精度。在放射治疗过程中，锥束CT用来实时获得病人图像，对放射治疗进行图像引导，图像精度的降低导致了放射治疗的位置误差，无法达到对肿瘤区域精确放疗的目的。
目前对图像中散射线的矫正方法主要有两种，一种是硬件矫正，是指在X射线影像***中添加若干硬件工具，减少达到电子射野影像装置的散射线，从而达到散射矫正目的。常用的硬件矫正工具包括X射线准直器、防散射滤线栅等。第二种是软件矫正方法，是指对已经采集到的X射线投影图像通过对投影图像的分析和被照射物体的估计，得到一个散射分布图，由散射分布图对X射线投影图像进行散射矫正。
发明内容
本发明要解决的问题是提供一种减少X射线图像散射中散射成分，提高X射线图像的图像质量的方法和装置。
为解决上述问题，本发明提供了一种，用于形成调制X射线图像的罩壳，所述罩壳围绕床板，设置在床板的外部，包括：沿第一方向延伸的第一部件；所述第一部件在与所述第一方向垂直的第二方向上间隔重复排列；所述第一部件在沿第一方向延伸的两侧边缘相互平行；所述第一部件对X射线产生遮挡。
优选的，所述罩壳固定于病床。
优选的，所述罩壳固定于X射线成像设备的孔径处。
优选的，所述第一方向和探测器方向相同。
优选的，所述第一方向与床板的长度方向垂直。
优选的，还包括：所述第一部件同时在第一方向上间隔重复排列。
优选的，还包括：沿所述第一方向延伸的第二部件；所述第二部件在与所述第一方向垂直的第二方向上间隔重复排列，和所述第一部件交替设置；所述第二部件在沿第一方向延伸的两侧边缘相互平行；所述第一部件对X射线的遮挡作用远大于所述第二部件。
本发明还提供了一种用于形成调制X射线图像的罩壳，所述罩壳围绕床板，设置在床板的外部，包括：沿第一方向延伸的第三部件；所述第三部件的厚度分布为：沿所述第一方向上相同，在沿与所述第一方向垂直的第二方向上呈周期性分布；所述第三部件对X射线产生遮挡。
本发明还提供了一种X射线图像散射成分的计算方法，包括：接收依次穿过上述所述的罩壳、成像对象的X射线，形成调制X射线图像；对所述调制X射线图像进行低通滤波，得到所述调制X射线图像的低频部分和散射成分；对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分；对所述调制X射线图像的高频部分解调制后结合权重，计算得到所述调制X射线低频部分的估计值；将经所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
本发明还提供了一种X射线图像的重建方法，包括：采集得到成像对象各投影角度下的X射线投影图像；根据上述的计算方法计算得到所述投影图像的散射成分；将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
在上述X射线图像散射成分计算方法、X射线图像重建方法的基础上，本发明还分别提供了对应的X射线图像散射成分计算装置、X射线图像重建装置。
与现有技术相比，本发明的技术方案采用罩壳对X射线进行强度调制，实现对X射线图像中散射成分的计算，有效去除了X射线图像的散射成分，提高了图像质量。更进一步的，罩壳的固定方式稳定，不会因为振动引起图像的运动模糊。
附图说明
图1是一种放射治疗***的结构图；
图2是本发明X射线图像散射成分计算方法流程图；
图3是从a方向看本发明罩壳和病床的位置关系图；
图4是本发明罩壳的结构图；
图5是本发明调制X射线图像的示意图；
图6是本发明散射矫正流程示意图；
图7是本发明另一种实施方式中罩壳表面厚度分布示意图；
图8是本发明另一种实施方式中调制X射线图像示意图；
图9是本发明另一种实施方式中罩壳的结构图；
图10是本发明X射线图像重建的方法流程图；
图11是本发明X射线图像重建装置的结构图。
具体实施方式
为使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发明的具体实施方式做详细的说明。在以下描述中阐述了具体细节以便于充分理解本发明。但是本发明能够以多种不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广。因此本发明不受下面公开的具体实施方式的限制。
图1是一种放射治疗***的结构图，如图1所示，放射治疗***100包括固定部分101和旋转部分102，旋转部分102安装在固定部分101上，旋转部分101可以绕中心轴106进行旋转，从而实现在不同角度对患者进行放射治疗。
旋转部分102的一侧是治疗头103，治疗头103可以产生高能级的X射线(通常为兆伏级)，对在病床105上的患者进行放射治疗。病床可沿轴107方向(也是病床的长度方向)载着患者运动到治疗头103下方。对于同源双束的放射治疗***，治疗头103还可以产生低能级的X射线(通常为千伏级)，低能级X射线可以用来对患者进行成像，利用得到的患者图像对患者进行图像引导放疗(ImageGuidedRadiationTherapy)。放射治疗设备100在旋转部分102的中心还存在孔径108，病床105可以载着患者进入到孔径108内，以便更好的进行X射线成像或治疗。
在由低能级X射线进行成像时，治疗头103发出锥束X射线，旋转部分102另一侧的电子射野影像装置104接收到穿过患者的X射线，形成该角度下的投影图像(projectionimage)。治疗头103在不同角度照射，从而形成多个角度的投影图，再对投影图进行图像重建(imagereconstruction)得到患者的CT(computedtomography)图像。
图2是本发明X射线图像散射成分的计算方法流程图，参见图2，包括以下步骤：
S201，获得调制X射线图像。
结合图3至图5对调制X射线图像的获得方法进行说明。在本发明的技术方案中，提出了一种利用罩壳对X射线进行调制的方法，罩壳固定安装在病床的外部。“调制X射线图像”即为X射线穿过该罩壳，再穿过成像对象，最终在探测器上采集得到的图像。
图3为从a方向(结合图1)看病床105和罩壳300的位置关系示意图，罩壳300安装在病床105的外部并围绕病床105设置，罩壳300优选的采用中空的薄层结构。在一种优选的实施方式中，罩壳300固定于病床105，在某些情况下为了得到最佳的成像角度，病床105会发生整体的旋转或者平移，这个时候罩壳300会与病床105一起运动。在一种优选的实施方式中，罩壳300固定于X射线成像设备的孔径108处。在其他的实施方式中，罩壳300也可以直接固定在地面、固定在墙体，或者固定在其他的在成像时不会发生运动或者振动的装置处。
罩壳300的具体结构参见图4，罩壳300可以是由可对X射线进行遮挡的多个第一部件301构成，第一部件301所采用的材料可以对X射线进行有效的遮挡。这里所说的遮挡，可以是仅遮挡X射线强度的一部分，使其有一定的透射率，也可以是对X射线的完全遮挡。第一部件301在第一方向303延伸，各个第一部件301在与第一方向303垂直的第二方向(此时可以认为是罩壳的圆周方向)上间隔重复排列。第一部件301在沿第一方向延伸的两侧边缘3011、3012相互平行。第一部件301的厚度均匀，对不同的空间位置处保持对X射线的透射率相同。第一方向303此时是和病床的长度方向保持一致，在本技术方案的实际实施过程中，第一方向也可以不和病床长度方向相同，即罩壳整体的方向可以和病床方向呈一定的夹角。优选的，罩壳整体的方向和探测器104方向一致，这里探测器方向指的是一排探测器单元的排列方向。
罩壳300的形状并不局限于圆环形，其也可以是方柱形、三棱柱形或者多边棱柱形。罩壳300也可以不是完整的周向方向的一整圈，可以仅覆盖周向的一部分，例如半圈，或者小于半圈。
在本发明的一种实施方式中，第一部件301在第一方向303上可以是间隔性的延伸，即在第一方向303上保持一定间隔的重复性设置。结合在第二方向上的间隔重复排列，第一部件301整体可以是n*n的矩阵式设置方式。
罩壳300可以仅由多个第一部件301以及在第一部件301边缘起固定作用的框架组成，整体类似于一个鸟笼。也可以在各第一部件301的间隙之间另设置第二部件302，第二部件302同样沿第一方向303延伸，并且第二部件在与第一方向303垂直的第二方向上重复排列，和第一部件301交替设置，填补了第一部件301之间的空隙，整体组成一个圆环形状类似于一个纸筒。但如果要设置第二部件302，在对第二部件302材料的选择上，需要第一部件301对X射线遮挡作用远大于第二部件301对X射线的遮挡作用，这样能够使最终得到的调制X射线图像效果更好，若两者所用材料接近则效果较差。这里所说的“远大于”是指其对X射线的遮挡效果在数倍甚至数十倍以上，从而对经过第一部件301和第二部件302的X射线强度能够在图像上产生显著的差别。
由治疗头103发出的低能级X射线，经过罩壳300进行了X射线的调制，之后穿过成像对象，能够在电子射野影像装置104上形成如图5中所示的明暗相间条纹的图像500。其中，亮条纹的白色代表的是X射线没有经过罩壳300的遮挡，或者只是经过了第二部件302的遮挡在探测器104处采集得到的；而暗条纹的黑色代表的是X射线经过罩壳300的第一部件301被较大程度遮挡后，射线强度下降，在探测器104处采集得到的。最终形成的图像500可以称为调制X射线图像。
第一部件301在第二方向间隔重复排列，其具体的间隔距离可以根据调制X射线图像上明暗相间条纹的间隔距离所对应探测器上的物理间隔距离，结合射线源焦点到探测器距离、射线源焦点到罩壳距离，由对应比例计算得到。之后的步骤S202至步骤S205可以参考公开号为US7463712B2的美国专利说明书记载内容。以下结合图6，仅作简单描述说明。
对X射线图像500在计算时，可以对每个明暗条纹的中心线处进行下采样，避免图像中准直器边缘处的边缘效应对图像上X射线强度产生的影响。
S202，对调制X射线图像进行低通滤波，得到经调制X射线图像的低频部分和散射成分。
由于X射线图像中散射成分的特性，步骤S201获取的调制X射线图像中散射成分主要分布在图像数据的低频部分。当对投影图601(这里的投影图是指一定投影角度下的调制X射线图像500)进行低通滤波后，可以得到图像的低频部分和散射成分之和602的图像数据。
S203，对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分。
而对投影图601进行高通滤波，由于散射成分几乎不包含在调制X射线图像的高频部分，可以得到图像不包含散射成分的高频部分603。
S204，对所述高频部分进行解调制并结合权重，计算得到所述调制X射线低频部分的估计值。
将高频部分数据603解调制后乘以权重系数其中∝为准直器对X射线的透射系数，从而计算得到低频部分的估计值604。
S205，将所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
由步骤S202低通滤波后得到的低频部分和散射成分602之和，减去步骤S204计算得到的低频部分估计604，得到散射成分估计605。
在本发明技术方案的一种实施方式中，罩壳可以由一种完整覆盖罩壳周向的材料(这里称为第三部件)形成，第三部件在第一方向延伸，其厚度分布在第一方向上其厚度均匀，在第二方向上采用图7所示的结构，图7是沿第二方向的截面图，在第二方向上，厚度分布呈周期分布。形成的调制X射线图像情况参见图8，信号强度呈周期性分布。若采用本发明之前的实施方式(即采用第一部件或者第一部件和第二部件结合的方式)时，罩壳的厚度分布均匀，产生的调制X射线图像如图8左侧所示，在空间上的X射线数据强度是方波信号。而当采用本实施方式的罩壳(厚度呈周期性分布)时，图像如右侧所示，在空间上的X射线数据可以是余弦信号。采用余弦信号形式的调制X射线图像也同样可以执行S202至S205的步骤。这种方式下罩壳厚度的具体分布数值。可以采用预期合理的X射线分布，之后通过反投影的方式确定各部分的厚度分布。
图9是另一种罩壳300的结构示意图，此时，第一方向不再是303方向，第一方向为罩壳的圆周方向，第一部件沿罩壳的圆周方向延伸。此时与第一方向垂直的第二方向即是方向303。
本发明还提供了一种X射线图像的重建方法，如图10所示，包括：
S1001，采集得到各投影角度下的X射线投影图像。
如图1所示的放射治疗设备在进行成像时，旋转机构102旋转，治疗头103发出低能级X射线，探测器104接收多个旋转角度下成像对象的X射线投影图像，此时采集到的投影图像是常用X射线成像方法采集得到的投影图像。
S1002，由上述X射线图像散射成分的计算方法计算得到投影图像的散射成分。
由前述的X射线图像散射成分计算方法，首先采集到成像对象的调制X射线图像，之后计算得到调制X射线图像的散射成分，由于成像对象相同，散射情况也相同，所以该方法计算得到的散射成分也是步骤S1001得到的X射线投影图像的散射成分。
S1003，将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
将步骤S1001得到的投影图像减去步骤S1002计算得到的散射成分，完成对投影图像的散射去除，之后对投影图像进行图像重建，重建得到成像对象的CT图像。具体图像重建过程为本领域技术人员已知，这里不再赘述。在进行X射线图像重建时，可以采用周向覆盖角度小于或等于半圈的罩壳进行快速重建。在进行多个投影角度的X射线图像采集时，可以用电机控制罩壳和机器旋转部分的旋转速度相同，在罩壳和治疗头共同旋转(相对静止)的情况下进行不同投影角度下投影图像的采集。
与亮条纹投影区域相比，暗条纹投影区域的泊松量子噪声增大，图像信噪比较低。为了减少噪声引起的图像质量下降，在每个投影角度进行X射线图像采集时，可以采集两幅投影图像，两幅投影图像之间的明暗相间条纹互补，利用两幅图像的亮条纹组合起来进行后续处理。例如可以进行两圈的旋转扫描，第一圈旋转扫描同步骤S201相同，第一圈结束后，通过电机控制旋转罩壳，使图像的明暗相间条纹位置互换，重复步骤201采集第二圈图像，利用两次图像的亮条纹组合后续处理，进一步提高图像质量。
在上述X射线图像散射成分计算方法的基础上，本发明还提供了一种X射线图像散射成分计算装置，参见图11，X射线图像散射成分计算装置1110包括：
图像获取单元1111，用于接收依次穿过罩壳、成像对象的X射线，获得调制X射线图像；
低通滤波单元1112，用于对所述调制X射线图像进行低通滤波，得到所述调制X射线图像的低频部分和散射成分；
高通滤波单元1113，用于对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分；
低频部分估计单元1114，用于对所述调制X射线图像的高频部分解调制后结合权重，计算得到所述调制X射线低频部分的估计值；
散射计算单元1115，用于将经所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
在本发明X射线图像重建方法的基础上，本发明还提供了一种X射线图像重建装置，参见图11，X射线图像重建装置1120包括：
采集单元1121，用于采集得到成像对象各投影角度下的投影图像；
射线图像散射成分计算单元1110，用于计算所述投影图像的散射成分；
重建单元1122，用于将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
本发明X射线图像散射成分计算装置以及X射线图像重建装置的实施方式可以分别参考对应方法的实施方式，这里不再赘述。
本发明的技术方案采用罩壳对X射线进行强度调制，实现对X射线图像中散射成分的计算，有效去除了X射线图像的散射成分，提高了图像质量。更进一步的，罩壳的固定方式稳定，不会因为振动引起图像的运动模糊。
本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。
Claims (12)
1.一种用于形成调制X射线图像的罩壳，其特征在于，所述罩壳围绕床板，设置在床板的外部，包括：
沿第一方向延伸的第一部件；所述第一部件在与所述第一方向垂直的第二方向上间隔重复排列；
所述第一部件在沿第一方向延伸的两侧边缘相互平行；
所述第一部件对X射线产生遮挡。
2.根据权利要求1所述的罩壳，其特征在于，所述罩壳固定于病床。
3.根据权利要求1所述的罩壳，其特征在于，所述罩壳固定于X射线成像设备的孔径处。
4.根据权利要求1所述的罩壳，其特征在于，所述第一方向和探测器方向相同。
5.根据权利要求1所述的罩壳，其特征在于，所述第一方向与床板的长度方向垂直。
6.根据权利要求1所述的罩壳，其特征在于，还包括：所述第一部件同时在第一方向上间隔重复排列。
7.根据权利要求1所述的罩壳，其特征在于，还包括：
沿所述第一方向延伸的第二部件；
所述第二部件在与所述第一方向垂直的第二方向上间隔重复排列，和所述第一部件交替设置；
所述第二部件在沿第一方向延伸的两侧边缘相互平行；
所述第一部件对X射线的遮挡作用远大于所述第二部件。
8.一种用于形成调制X射线图像的罩壳，其特征在于，所述罩壳围绕床板，设置在床板的外部，包括：
沿第一方向延伸的第三部件；
所述第三部件的厚度分布为：沿所述第一方向上相同，在沿与所述第一方向垂直的第二方向上呈周期性分布；
所述第三部件对X射线产生遮挡。
9.一种X射线图像散射成分的计算方法，其特征在于，包括：
接收依次穿过权利要求1或8任一项所述的罩壳、成像对象的X射线，形成调制X射线图像；
对所述调制X射线图像进行低通滤波，得到所述调制X射线图像的低频部分和散射成分；
对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分；
对所述调制X射线图像的高频部分解调制后结合权重，计算得到所述调制X射线低频部分的估计值；
将经所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
10.一种X射线图像的重建方法，其特征在于，包括：
采集得到成像对象各投影角度下的X射线投影图像；
根据权利要求9所述的计算方法计算得到所述投影图像的散射成分；
将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
11.一种X射线图像散射成分计算装置，其特征在于，包括：
图像获取单元，用于接收依次穿过权利要求1或8任一项所述的罩壳、成像对象的X射线，获得调制X射线图像；
低通滤波单元，用于对所述调制X射线图像进行低通滤波，得到所述调制X射线图像的低频部分和散射成分；
高通滤波单元，用于对所述调制X射线图像进行高通滤波，得到所述调制X射线图像的高频部分；
低频部分估计单元，用于对所述调制X射线图像的高频部分解调制后结合权重，计算得到所述调制X射线低频部分的估计值；
散射计算单元，用于将经所述低通滤波后得到的调制X射线图像的低频部分和散射成分，减去所述调制X射线低频部分的估计值，计算得到散射成分。
12.一种X射线图像重建装置，其特征在于，包括：
采集单元，用于采集得到成像对象各投影角度下的投影图像；
权利要求11所述的X射线图像散射成分计算单元，用于计算所述投影图像的散射成分；
重建单元，用于将所述投影图像去掉所述散射成分后进行图像重建，得到重建图像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510834741.4A CN105427355B (zh) | 2015-11-26 | 2015-11-26 | 罩壳、x射线图像的散射成分计算、重建的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510834741.4A CN105427355B (zh) | 2015-11-26 | 2015-11-26 | 罩壳、x射线图像的散射成分计算、重建的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105427355A true CN105427355A (zh) | 2016-03-23 |
| CN105427355B CN105427355B (zh) | 2018-04-03 |
Family
ID=55505535
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510834741.4A Active CN105427355B (zh) | 2015-11-26 | 2015-11-26 | 罩壳、x射线图像的散射成分计算、重建的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105427355B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109646821A (zh) * | 2019-01-30 | 2019-04-19 | 张文山 | 一种实时图像引导的技术方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1207888A (zh) * | 1997-05-22 | 1999-02-17 | 西门子公司 | 用于生成x射线投影的计算机x射线断层造影机 |
| JP3721632B2 (ja) * | 1996-03-31 | 2005-11-30 | 株式会社島津製作所 | Ｘ線画像デジタル処理装置 |
| CN1849672A (zh) * | 2003-09-12 | 2006-10-18 | 皇家飞利浦电子股份有限公司 | 用于准直电磁辐射的装置 |
| CN102764137A (zh) * | 2012-07-27 | 2012-11-07 | 苏州生物医学工程技术研究所 | 一种静态ct扫描仪及其散射x光子校正方法 |
| CN202715140U (zh) * | 2012-08-24 | 2013-02-06 | 山东新华医疗器械股份有限公司 | 适形调强放射治疗装置 |
| CN103330570A (zh) * | 2013-04-27 | 2013-10-02 | 中国人民解放军北京军区总医院 | X射线准直器、x射线准直***及移动ct扫描仪 |
-
2015
- 2015-11-26 CN CN201510834741.4A patent/CN105427355B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3721632B2 (ja) * | 1996-03-31 | 2005-11-30 | 株式会社島津製作所 | Ｘ線画像デジタル処理装置 |
| CN1207888A (zh) * | 1997-05-22 | 1999-02-17 | 西门子公司 | 用于生成x射线投影的计算机x射线断层造影机 |
| CN1849672A (zh) * | 2003-09-12 | 2006-10-18 | 皇家飞利浦电子股份有限公司 | 用于准直电磁辐射的装置 |
| CN102764137A (zh) * | 2012-07-27 | 2012-11-07 | 苏州生物医学工程技术研究所 | 一种静态ct扫描仪及其散射x光子校正方法 |
| CN202715140U (zh) * | 2012-08-24 | 2013-02-06 | 山东新华医疗器械股份有限公司 | 适形调强放射治疗装置 |
| CN103330570A (zh) * | 2013-04-27 | 2013-10-02 | 中国人民解放军北京军区总医院 | X射线准直器、x射线准直***及移动ct扫描仪 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109646821A (zh) * | 2019-01-30 | 2019-04-19 | 张文山 | 一种实时图像引导的技术方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105427355B (zh) | 2018-04-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105326519B (zh) | X射线图像形成、散射成分计算以及重建的方法及装置 | |
| CN101987021B (zh) | Ct***的散射校正方法及ct*** | |
| CN103384498B (zh) | 探测装置 | |
| JP6066596B2 (ja) | Ｘ線撮像における散乱補正の方法及びシステム | |
| US8270561B2 (en) | Motion weighting in computed tomography (CT) with cone angle | |
| US20110176723A1 (en) | Motion Correction in Cone-Beam CT by Tracking Internal and External Markers Using Cone-Beam Projection From a kV On-Board Imager: Four-Dimensional Cone-Beam CT and Tumor Tracking Implications | |
| US8406373B2 (en) | Optimized aperture selection imaging computed tomography system and method | |
| JP2011502679A (ja) | 三次元ｘ線画像における改良された画質のための可動式くさび | |
| KR20100126658A (ko) | 치과용 엑스선 장치 및 연관된 방법 | |
| JP2005152658A (ja) | 非一様ビュー加重式トモシンセシス方法及び装置 | |
| JP5642444B2 (ja) | 放射線治療装置の作動方法および放射線治療装置制御装置 | |
| CN103415253B (zh) | 用于计算机断层摄影（ct）的患者后动态过滤器 | |
| US11350892B2 (en) | Collimator structure for an imaging system | |
| EP1478273B1 (en) | Sequential computed tomography method | |
| Shimohigashi et al. | Image quality of four-dimensional cone-beam computed tomography obtained at various gantry rotation speeds for liver stereotactic body radiation therapy with fiducial markers | |
| Mainprize et al. | Dependence of image quality on geometric factors in breast tomosynthesis | |
| US9901314B2 (en) | Adjustable bow-tie filter for achieving optimal SNR in helical computed tomography | |
| CN105427355A (zh) | 罩壳、x射线图像的散射成分计算、重建的方法及装置 | |
| US9545231B2 (en) | Optimized aperture selection imaging computed tomography system and method | |
| Du et al. | TIGRE-VarianCBCT for on-board cone-beam computed tomography, an open-source toolkit for imaging, dosimetry and clinical research | |
| Chen et al. | Super-resolution imaging in a multiple layer EPID | |
| Saha et al. | Single photon emission computed tomography | |
| JP2014230615A (ja) | 放射線断層撮影装置および投影データ補正方法並びにプログラム | |
| Stankovic et al. | Improving linac integrated cone beam computed tomography image quality using tube current modulation | |
| Cai | Artefact Reduction Methods for Iterative Reconstruction in Full-fan Cone Beam CT Radiotherapy Applications |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into 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 |
Address after: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 No. 2258 Chengbei Road, Jiading Industrial Zone, Jiading District, Shanghai. Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |