CN107928690B - 确定主射线比例方法及图像去散射方法、系统及存储介质 - Google Patents
确定主射线比例方法及图像去散射方法、系统及存储介质 Download PDFInfo
- Publication number
- CN107928690B CN107928690B CN201711489052.XA CN201711489052A CN107928690B CN 107928690 B CN107928690 B CN 107928690B CN 201711489052 A CN201711489052 A CN 201711489052A CN 107928690 B CN107928690 B CN 107928690B
- Authority
- CN
- China
- Prior art keywords
- photon
- received
- detector crystal
- detector
- imaging condition
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000000790 scattering method Methods 0.000 title description 3
- 239000013078 crystal Substances 0.000 claims abstract description 123
- 230000005855 radiation Effects 0.000 claims abstract description 15
- 238000003384 imaging method Methods 0.000 claims description 90
- 230000015654 memory Effects 0.000 claims description 23
- 238000004590 computer program Methods 0.000 claims description 13
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000002601 radiography Methods 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 2
- 238000007408 cone-beam computed tomography Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010295 mobile communication Methods 0.000 description 2
- 229920003229 poly(methyl methacrylate) Polymers 0.000 description 2
- 239000004926 polymethyl methacrylate Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000010977 jade Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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/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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种确定探测器晶体接收到的主射线比例的方法及医学图像的去散射方法、系统及存储介质。确定探测器晶体接收到的主射线比例的方法包括：获取射线经过模体后，探测器接收到的光子信息；基于探测器晶体位置及光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量；基于所述主射光子能量及所有光子能量获得探测器晶体接收到的主射线比例。本实施例提供了可以通过软件的方式去除散射线，不需要通过滤线栅的方式来去除散射线，在一定程度上降低了设备的成本的同时，也提高了设备的普适度。并且还降低了受检者接收到的辐射剂量，避免受检者接收不必要的辐射。
Description
技术领域
本发明涉及医疗设备成像技术领域，特别是涉及一种确定探测器晶体接收到的主射线比例的方法及医学图像的去散射方法、系统及存储介质。
背景技术
数字化X射线摄影系统(Digital Radiography，DR)是一种新兴的医学成像方法，它利用平板探测器接收X射线，直接获得数字图像信号，具有图像分辨率高、动态范围宽，成像速度快，对人体辐射小等显著优点，在临床医学和科学研究中都得到了广泛的应用。
对于X射线摄影系统而言，球管发射的X射线穿过人体后被探测器所接收，由于X射线在穿过人体组织、器官以及骨骼时会发生不同程度的衰减，因此探测器的各个晶体接收到的剂量有所不同，进而输出对应的投影图像。由于X射线经过人体时，会在人体内发生散射，形成散射线，因此，被探测器接收到的光线中，除了主射线外，还包括经过人体后形成的散射线。而在生成受检者的投影图像时，所需要的只有主射线，因此散射线就成为了干扰。为了消除散射线干扰，传统技术使用滤线栅对散射线进行滤除。
然而采用滤线栅对散射线进行滤除时，滤线栅只能滤除掉大部分散射线，并不能将所有的散射线滤除，且滤线栅还会滤除掉部分的主射线，进而影响投影图像的质量。而为了提高投影图像的质量，势必会增加X射线的剂量，进而导致受检者接收过多的X射线剂量，对其造成一定的伤害。
发明内容
基于此，有必要针对现有通过滤线栅去散射可能会导致受检者接收过多剂量，提供一种确定探测器晶体接收到的主射线比例的方法及医学图像的去散射方法、系统及存储介质。
本发明实施例提供了一种确定探测器晶体接收到的主射线比例的方法包括：获取射线经过模体后，探测器接收到的光子信息；基于探测器晶体位置及所述光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量；基于所述主射光子能量及所有光子能量获得探测器晶体接收到的主射线比例。
进一步，所述光子信息包括:光子角度及能量，所述光子包括：主射光子及散射光子。
进一步，基于蒙特卡罗方法统计探测器接收到的光子信息。
本发明实施例还提供了一种医学图像的去散射方法包括：获取预设成像条件下射线经过模体后，探测器接收到的光子信息；基于所述预设成像条件下探测器晶体位置及所述光子信息确定所述预设成像条件下所述探测器晶体接收到的所有光子能量及主射光子能量；基于所述预设成像条件下所述主射光子能量及所有光子能量获得所述预设成像条件下探测器晶体接收到的主射线比例；获取医学图像及采集所述医学图像时的成像条件；基于与所述成像条件相匹配的预设成像条件下探测器晶体接收到的主射线比例对所述医学图像去散射。
进一步，所述基于与所述成像条件相匹配的预设成像条件下探测器晶体接收到的主射线比例对所述医学图像去散射包括：将所述医学图像各像素点的灰度值乘以与该像素点对应的探测器晶体接收到的主射线比例以获得去散射后的医学图像。
本发明实施例还提供了一种确定探测器晶体接收到的主射线比例的系统，包括计算机，所述计算机包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时可用于执行一种确定探测器晶体接收到的主射线比例的方法，所述方法包括：获取射线经过模体后，探测器接收到的光子信息；基于探测器晶体位置及所述光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量；基于所述主射光子能量及所有光子能量获得探测器晶体接收到的主射线比例。
本发明实施例还提供了一种医学图像的去散射系统，包括计算机，所述计算机包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时可用于执行一种医学图像的去散射方法，所述方法包括：获取预设成像条件下射线经过模体后，探测器接收到的光子信息；基于所述预设成像条件下探测器晶体位置及所述光子信息确定所述预设成像条件下所述探测器晶体接收到的所有光子能量及主射光子能量；基于所述预设成像条件下所述主射光子能量及所有光子能量获得所述预设成像条件下探测器晶体接收到的主射线比例；获取医学图像及采集所述医学图像时的成像条件；基于与所述成像条件相匹配的预设成像条件下探测器晶体接收到的主射线比例对所述医学图像去散射。
进一步，所述处理器执行所述程序时可用于执行一种医学图像的去散射方法，所述方法还包括：将所述医学图像各像素点的灰度值乘以与该像素点对应的探测器晶体的主射线比例以获得去散射后的医学图像。
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时可用于执行一种确定探测器晶体接收到的主射线比例的方法中任一项所述的方法。
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时可用于执行一种医学图像的去散射方法中任一项所述的方法。
本发明通过获取射线经过模体后，探测器接收到的光子信息，并基于探测器晶体位置及所述光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量；进而确定探测器晶体接收到的主射线的比例，提供了可以通过软件的方式去除散射线，不需要通过滤线栅的方式来去除散射线，在一定程度上降低了设备的成本的同时，也提高了设备的普适度。
通过获取预设成像条件下探测器晶体接收到的主射线比例，并对采集医学图像时的成像条件和预设成像条件进行匹配，以获得与成像条件对应的预设成像条件下探测器晶体接收到的主射线比例，进而去除医学图像的散射，较好的去除了散射线对生成医学图像时的干扰，获得了符合实际临床需求的图像的同时也降低了受检者接收到的辐射剂量，避免受检者接收不必要的辐射。
附图说明
图1为本发明实施例提供的一种确定探测器晶体接收到的主射线比例的方法的流程图；
图2为本发明实施例提供的一种医学图像的去散射方法的流程图；
图3为本发明实施例提供的一种确定探测器晶体接收到的主射线比例的系统的结构示意图；
图4为本发明实施例提供的一种医学图像的去散射系统的结构示意图。
附图标记：610、810为处理器，620、820为存储器，630、830为输入装置，640、840为输出装置。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明实施例中的医学图像可以为利用数字化X射线摄影系统或锥形束计算机断层扫描系统或计算机断层扫描系统等各种利用X射线生成的医学图像，本申请不对获取医学图像所采用的设备进行限定。
本发明实施例公开了一种确定探测器晶体接收到的主射线比例的方法，通过获取探测器接收到的所有光子信息，基于所有光子的光子信息得到探测器各个晶体接收到的光子能量及其中主射光子的能量，进一步的计算出探测器各个晶体主射线比例。
请参阅图1，图1为本发明实施例提供的一种确定探测器晶体接收到的主射线比例的方法的流程图。
如图1所示，一种确定探测器晶体接收到的主射线比例的方法可以包括以下步骤S110至S130：
步骤S110：获取射线经过模体后，探测器接收到的光子信息。
具体的，可以使用EGSnrc或Geant4等利用蒙特卡罗算法的工具模拟成像过程，例如：射线经过模体(可以为人体)后，探测器接收到射线的过程，并通过模拟获得的光子信息文件来得到光子信息；此外也可以利用实际测量工具对模体进行检测，再利用蒙特卡罗算法对检测到的数据进行统计确定。本申请不对获取光子信息的具体方式进行限定。本实施例中光子信息包括：光子角度及能量，所述光子包括：主射光子及散射光子。其中模体模拟了人体的组成，可以通过蒙特卡洛软件进行模拟，也可以为可以模拟人体组成的水模或者PMMA模体等。
步骤S120：基于探测器晶体位置及所述光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量。
具体的，根据探测器晶体(探测器包括多个晶体)的位置以及探测器晶体所接收到的所有光子的光子角度，来确定该光子是否被相应的探测器晶体捕获，或者说确定与探测器晶体相交的光子有哪些。而和探测器晶体相交的光子中包括主射光子(与模体未发生反应的光子)和散射光子。实际应用中，在采用蒙卡软件模拟成像过程时，可以将模体在蒙卡软件中作标记，以知悉光子信息文件中哪些光子是主射光子，哪些光子是散射光子；或者通过判断哪些光子与模体未发生碰撞(光子进入模体前后其入射方向和能量都不变)，此类光子为主射光子，如果与模体发生碰撞，则此类光子为散射光子。此外光子信息文件中也会给出不同光子的能量。将探测器每一个晶体接收到的所有光子的能量相加，则得到所有光子能量；将探测器每一个晶体接收到的主射光子的光子能量相加，得到主射光子能量。
步骤S130：基于所述主射光子能量及所有光子能量获得探测器晶体接收到的主射线比例。
具体的，将探测器晶体接收到的主射光子的能量与探测器晶体接收到的所有光子的能量相比得到探测器该晶体接收到的主射光子的能量占探测器接收到的所有光子的能量的比例，即主射线比例。依次通过同样的方式即可获得探测器每一个晶体接收到的主射线比例。
请参阅图2，图2为本发明实施例提供的一种医学图像的去散射方法的流程图。
如图2所示，一种医学图像的去散射方法可以包括以下步骤S210至S250：
步骤S210：获取预设成像条件下射线经过模体后，探测器接收到的光子信息。
具体的，可以使用EGSnrc或Geant4等利用蒙特卡罗算法的工具模拟成像过程，例如：射线经过模体(可以为人体)后，探测器接收到射线的过程，并通过模拟获得的光子信息文件来得到光子信息；此外也可以利用实际测量工具对模体进行检测，再利用蒙特卡罗算法对检测到的数据进行统计确定。本申请不对获取光子信息的具体方式进行限定。其中光子信息包括：光子角度及能量，所述光子包括：主射光子及散射光子。其中所述成像条件为：射野尺寸以及受检者/模体与探测器之间的距离。具体的：射野尺寸为球管呈圆锥形射出的射线，在限束器的作用下将圆锥形光束遮挡为矩形光束，所述矩形光束在探测器平面的投影尺寸即为射野尺寸。其中模体模拟了人体的组成，可以通过蒙特卡洛软件进行模拟，也可以为可以模拟人体组成的水模或者PMMA模体等。其中各种预设成像条件为当射野尺寸不同时，选择不同的模体与探测器之间的距离。不同的射野尺寸以及模体与探测器之间的距离能够满足采集不同受检者不同部位的图像时所需的成像条件。
步骤S220：基于所述预设成像条件下探测器晶体位置及所述光子信息确定所述预设成像条件下所述探测器晶体接收到的所有光子能量及主射光子能量。
具体的，根据探测器晶体(探测器包括多个晶体)的位置以及探测器晶体所接收到的所有光子的光子角度，来确定该光子是否被相应的探测器晶体捕获，或者说确定与探测器晶体相交的光子有哪些。而和探测器晶体相交的光子中包括主射光子(与模体未发生反应的光子)和散射光子。实际应用中，在采用蒙卡软件模拟成像过程时，可以将模体在蒙卡软件中作标记，以知悉光子信息文件中哪些光子是主射光子，哪些光子是散射光子；或者通过判断哪些光子与模体未发生碰撞(光子进入模体前后其入射方向和能量都不变)，此类光子为主射光子，如果与模体发生碰撞，则此类光子为散射光子。此外光子信息文件中也会给出不同光子的能量。将探测器每一个晶体接收到的所有光子的能量相加，则得到所有光子能量；将探测器每一个晶体接收到的主射光子的光子能量相加，得到主射光子能量。
步骤S230：基于所述预设成像条件下所述主射光子能量及所有光子能量获得所述预设成像条件下探测器晶体接收到的主射线比例。
具体的，将预设成像条件下探测器晶体接收到的主射光子能量与预设成像条件下探测器接收到的所有光子的所有光子能量相比得到预设成像条件下探测器晶体接收到的主射光子的能量占探测器接收到的所有光子的光子能量的比例，即主射线比例。计算预设成像条件下探测器晶体接收到的主射线比例。
步骤S240：获取医学图像及采集所述医学图像时的成像条件。
具体的，医学图像可以为利用数字化X射线摄影系统或锥形束计算机断层扫描系统或计算机断层扫描系统等各种利用X射线生成的医学图像。并且所述医学图像可以为实时检测得到的医学图像，还可以为之前检测后生成的医学图像，再对之前检测后生成的医学图像进行读取得到。基于医学图像获取采集医学图像时的成像条件，所述成像条件包括：射野尺寸及受检者者与探测器之间的距离。
步骤S250：基于与所述成像条件相匹配的预设成像条件下探测器晶体接收到的主射线比例对所述医学图像去散射。
具体的，基于采集医学图像时的成像条件匹配与采集医学图像时的成像条件相同的预设成像条件下探测器各个晶体的主射线比例。其中医学图像的各个像素点分别与预设成像条件下探测器各个晶体相对应，即医学图像中每一个像素点的灰度值是根据相对应的探测器的晶体接收的光子能量所生成。对所述医学图像去散射是将医学图像各像素点的灰度值乘以与该像素点对应的探测器晶体接收到的主射线比例以获得去散射后的医学图像。由于医学图像的各个像素点分别与预设成像条件下探测器各个晶体相对应，因此医学图像每个像素点对应一个探测器晶体，将医学图像的每一个像素点的灰度值分别乘以相应的探测器晶体的主射线比例，得到去散射后的医学图像。
请参阅图3，图3为本发明实施例提供的一种确定探测器晶体接收到的主射线比例的系统的结构示意图。
如图3所示，本发明实施例还提供了一种执行确定探测器晶体接收到的主射线比例的方法的确定探测器晶体接收到的主射线比例的系统，该系统包括计算机，其中计算机包括一个或多个处理器610以及存储器620，图3以一个处理器610为例进行说明。
执行确定探测器晶体接收到的主射线比例的方法的确定探测器晶体接收到的主射线比例的系统还可以包括：输入装置630和输出装置640。
处理器610、存储器620、输入装置630和输出装置640可以通过总线或其他方式连接，图3中以通过总线连接的方式为例。
本领域技术人员可以理解，图3中示出的系统的结构并不构成对本发明实施例的限定，它既可以是总线形结构，也可以是星型结构，还可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
处理器610可以由集成电路(IntegratedCircuit，简称IC)组成，例如可以由单颗封装的IC所组成，也可以由连接多颗相同功能或不同功能的封装IC而组成。举例来说，处理器610可以仅包括中央处理器(CentralProcessingUnit，简称CPU)，也可以是CPU、数字信号处理器(digitalsignalprocessor，简称DSP)、图形处理器(GraphicProcessingUnit，简称GPU)及各种控制芯片的组合。在本发明实施方式中，CPU可以是单运算核心，也可以包括多运算核心。
存储器620可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；在本发明实施例中，操作系统可以是Android系统、iOS系统或Windows操作系统等等。存储数据区可存储依据列表项操作的处理装置的使用所创建的数据等。此外，存储器620可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或者其他非暂态固态存储器件。在一些实施例中，存储器620可选包括相对于处理器610远程设置的存储器，这些远程存储器可以通过网络连接至列表项操作的处理装置。上述网络的实施例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置630可接收输入的数字或字符信息，以及产生与列表操作的处理装置的用户设置以及功能控制有关的键信号输入。输入装置可以包括触摸屏、键盘等，也可以包括有线接口、无线接口等。
输出装置640可包括显示屏、扬声器等设备，也可以包括有线接口、无线接口等。
所述处理器执行程序时可用于执行一种确定探测器晶体接收到的主射线比例的方法，所述方法包括：获取射线经过模体后，探测器接收到的光子信息；基于探测器晶体位置及所述光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量；基于所述主射光子能量及所有光子能量获得探测器晶体接收到的主射线比例。
请参阅图4，图4为本发明实施例提供的一种医学图像的去散射系统的结构示意图。如图4所示，本发明实施例还提供了一种执行医学图像的去散射方法的医学图像的去散射系统，该系统包括计算机，其中计算机包括一个或多个处理器810以及存储器820，图4以一个处理器810为例进行说明。
执行医学图像的去散射方法的医学图像的去散射系统还可以包括：输入装置830和输出装置840。
处理器810、存储器820、输入装置830和输出装置840可以通过总线或其他方式连接，图4中以通过总线连接的方式为例。
本领域技术人员可以理解，图4中示出的系统的结构并不构成对本发明实施例的限定，它既可以是总线形结构，也可以是星型结构，还可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件布置。
处理器810可以由集成电路(IntegratedCircuit，简称IC)组成，例如可以由单颗封装的IC所组成，也可以由连接多颗相同功能或不同功能的封装IC而组成。举例来说，处理器810可以仅包括中央处理器(CentralProcessingUnit，简称CPU)，也可以是CPU、数字信号处理器(digitalsignalprocessor，简称DSP)、图形处理器(GraphicProcessingUnit，简称GPU)及各种控制芯片的组合。在本发明实施方式中，CPU可以是单运算核心，也可以包括多运算核心。
存储器820可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；在本发明实施例中，操作系统可以是Android系统、iOS系统或Windows操作系统等等。存储数据区可存储依据列表项操作的处理装置的使用所创建的数据等。此外，存储器820可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或者其他非暂态固态存储器件。在一些实施例中，存储器820可选包括相对于处理器810远程设置的存储器，这些远程存储器可以通过网络连接至列表项操作的处理装置。上述网络的实施例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置830可接收输入的数字或字符信息，以及产生与列表操作的处理装置的用户设置以及功能控制有关的键信号输入。输入装置可以包括触摸屏、键盘等，也可以包括有线接口、无线接口等。
输出装置840可包括显示屏、扬声器等设备，也可以包括有线接口、无线接口等。
所述处理器执行程序时可用于执行一种医学图像的去散射方法，所述方法包括：获取预设成像条件下射线经过模体后，探测器接收到的光子信息；基于所述预设成像条件下探测器晶体位置及所述光子信息确定所述预设成像条件下所述探测器晶体接收到的所有光子能量及主射光子能量；基于所述预设成像条件下所述主射光子能量及所有光子能量获得所述预设成像条件下探测器晶体接收到的主射线比例；获取医学图像及采集所述医学图像时的成像条件；基于与所述成像条件相匹配的预设成像条件下探测器晶体接收到的主射线比例对所述医学图像去散射。
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时可用于执行本实施例中一种确定探测器晶体接收到的主射线比例的方法。所述方法包括：获取射线经过模体后，探测器接收到的光子信息；基于探测器晶体位置及所述光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量；基于所述主射光子能量及所有光子能量获得探测器晶体接收到的主射线比例。
本发明实施例还提供了另一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时可用于执行本实施例中一种医学图像的去散射方法。所述方法包括：获取预设成像条件下射线经过模体后，探测器接收到的光子信息；基于所述预设成像条件下探测器晶体位置及所述光子信息确定所述预设成像条件下所述探测器晶体接收到的所有光子能量及主射光子能量；基于所述预设成像条件下所述主射光子能量及所有光子能量获得所述预设成像条件下探测器晶体的接收到主射线比例；获取医学图像及采集所述医学图像时的成像条件；基于与所述成像条件相匹配的预设成像条件下探测器晶体接收到的主射线比例对所述医学图像去散射。
本申请各部分操作所需的计算机程序编码可以用任意一种或多种程序语言编写，包括面向对象编程语言如Java、Scala、Smalltalk、Eiffel、JADE、Emerald、C++、C#、VB.NET、Python等，常规程序化编程语言如C语言、Visual Basic、Fortran 2003、Perl、COBOL 2002、PHP、ABAP，动态编程语言如Python、Ruby和Groovy，或其他编程语言等。该程序编码可以完全在用户计算机上运行、或作为独立的软件包在用户计算机上运行、或部分在用户计算机上运行部分在远程计算机运行、或完全在远程计算机或服务器上运行。在后种情况下，远程计算机可以通过任何网络形式与用户计算机连接，比如局域网(LAN)或广域网(WAN)，或连接至外部计算机(例如通过因特网)，或在云计算环境中，或作为服务使用如软件即服务(SaaS)。
本发明通过获取射线经过模体后，探测器接收到的光子信息，并基于探测器晶体位置及所述光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量；进而确定探测器晶体接收到的主射线的比例，提供了可以通过软件的方式去除散射线，不需要通过滤线栅的方式来去除散射线，在一定程度上降低了设备的成本的同时，也提高了设备的普适度。
通过获取预设成像条件下探测器晶体接收到的主射线比例，并对采集医学图像时的成像条件和预设成像条件进行匹配，以获得与成像条件对应的预设成像条件下探测器晶体接收到的主射线比例，进而去除医学图像的散射，较好的去除了散射线对生成医学图像时的干扰，获得了符合实际临床需求的图像的同时也降低了受检者接收到的辐射剂量，避免受检者接收不必要的辐射。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种确定探测器晶体接收到的主射线比例的方法，其特征在于，包括：
获取射线经过模体后，探测器接收到的光子信息，其中，所述光子信息包括光子角度及能量；
基于探测器晶体位置及所述光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量；
基于所述主射光子能量及所有光子能量获得探测器晶体接收到的主射线比例。
2.如权利要求1所述的确定探测器晶体接收到的主射线比例的方法，其特征在于，所述光子包括：主射光子及散射光子。
3.如权利要求1所述的确定探测器晶体接收到的主射线比例的方法，其特征在于，基于蒙特卡罗方法统计探测器接收到的光子信息。
4.一种医学图像的去散射方法，其特征在于，包括：
获取预设成像条件下射线经过模体后，探测器接收到的光子信息，其中，所述光子信息包括光子角度及能量；
基于所述预设成像条件下探测器晶体位置及所述光子信息确定所述预设成像条件下所述探测器晶体接收到的所有光子能量及主射光子能量；
基于所述预设成像条件下所述主射光子能量及所有光子能量获得所述预设成像条件下探测器晶体接收到的主射线比例；
获取医学图像及采集所述医学图像时的成像条件；
基于与所述成像条件相匹配的预设成像条件下探测器晶体接收到的主射线比例对所述医学图像去散射。
5.如权利要求4所述的医学图像的去散射方法，其特征在于，所述基于与所述成像条件相匹配的预设成像条件下探测器晶体接收到的主射线比例对所述医学图像去散射包括：将所述医学图像各像素点的灰度值乘以与该像素点对应的探测器晶体接收到的主射线比例以获得去散射后的医学图像。
6.一种确定探测器晶体接收到的主射线比例的系统，包括计算机，所述计算机包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时可用于执行一种确定探测器晶体接收到的主射线比例的方法，所述方法包括：
获取射线经过模体后，探测器接收到的光子信息，其中，所述光子信息包括光子角度及能量；
基于探测器晶体位置及所述光子信息确定所述探测器晶体接收到的所有光子能量及主射光子能量；
基于所述主射光子能量及所有光子能量获得探测器晶体接收到的主射线比例。
7.一种医学图像的去散射系统，包括计算机，所述计算机包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时可用于执行一种医学图像的去散射方法，所述方法包括：
获取预设成像条件下射线经过模体后，探测器接收到的光子信息，其中，所述光子信息包括光子角度及能量；
基于所述预设成像条件下探测器晶体位置及所述光子信息确定所述预设成像条件下所述探测器晶体接收到的所有光子能量及主射光子能量；
基于所述预设成像条件下所述主射光子能量及所有光子能量获得所述预设成像条件下探测器晶体接收到的主射线比例；
获取医学图像及采集所述医学图像时的成像条件；
基于与所述成像条件相匹配的预设成像条件下探测器晶体接收到的主射线比例对所述医学图像去散射。
8.如权利要求7所述的去散射系统，其特征在于，所述处理器执行所述程序时可用于执行一种医学图像的去散射方法，所述方法还包括：将所述医学图像各像素点的灰度值乘以与该像素点对应的探测器晶体的主射线比例以获得去散射后的医学图像。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时可用于执行权利要求1-3任一项所述的方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时可用于执行权利要求4-5任一项所述的方法。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711489052.XA CN107928690B (zh) | 2017-12-29 | 2017-12-29 | 确定主射线比例方法及图像去散射方法、系统及存储介质 |
| PCT/CN2018/121049 WO2019128731A1 (en) | 2017-12-29 | 2018-12-14 | Systems and methods for scatter correction of image |
| EP18896349.0A EP3731759A4 (en) | 2017-12-29 | 2018-12-14 | SYSTEMS AND METHODS FOR SCRAP CORRECTION OF AN IMAGE |
| CN201880084682.5A CN111526796B (zh) | 2017-12-29 | 2018-12-14 | 用于图像散射校正的系统和方法 |
| US16/914,538 US11538160B2 (en) | 2017-12-29 | 2020-06-29 | Systems and methods for scatter correction of image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711489052.XA CN107928690B (zh) | 2017-12-29 | 2017-12-29 | 确定主射线比例方法及图像去散射方法、系统及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107928690A CN107928690A (zh) | 2018-04-20 |
| CN107928690B true CN107928690B (zh) | 2021-06-22 |
Family
ID=61937160
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711489052.XA Active CN107928690B (zh) | 2017-12-29 | 2017-12-29 | 确定主射线比例方法及图像去散射方法、系统及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107928690B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019128731A1 (en) * | 2017-12-29 | 2019-07-04 | Shenzhen United Imaging Healthcare Co., Ltd. | Systems and methods for scatter correction of image |
| CN108577872B (zh) * | 2017-12-29 | 2021-07-16 | 上海联影医疗科技股份有限公司 | 医学图像的去散射方法、系统及存储介质 |
| CN109655473B (zh) * | 2018-12-18 | 2021-10-22 | 北京应用物理与计算数学研究所 | 点探测器计数的闪光照相图像接收装置的模拟方法及系统 |
| CN111783292B (zh) * | 2020-06-23 | 2023-12-12 | 上海联影医疗科技股份有限公司 | X射线成像设备的建模方法、装置、设备及存储介质 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6163589A (en) * | 1998-06-13 | 2000-12-19 | General Electric Company | Monte Carlo scatter correction method for computed tomography of general object geometries |
| JP3919724B2 (ja) * | 2003-09-19 | 2007-05-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線計算断層画像装置および断層像データ生成方法 |
| WO2006082557A2 (en) * | 2005-02-01 | 2006-08-10 | Koninklijke Philips Electronics N.V. | Apparatus and method for correction or extension of x-ray projections |
| US9226722B2 (en) * | 2005-09-19 | 2016-01-05 | Yongwang Li | Medical imaging system and method with separate primary and scattered components |
| US8488902B2 (en) * | 2009-01-29 | 2013-07-16 | Koninklijke Philips Electronics N.V. | Detection values correction apparatus |
| JP5815048B2 (ja) * | 2011-12-12 | 2015-11-17 | 株式会社日立メディコ | Ｘ線ｃｔ装置 |
| CN104840211B (zh) * | 2015-05-18 | 2018-12-25 | 上海联影医疗科技有限公司 | 一种投影图像的散射矫正方法及装置 |
-
2017
- 2017-12-29 CN CN201711489052.XA patent/CN107928690B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107928690A (zh) | 2018-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107928690B (zh) | 确定主射线比例方法及图像去散射方法、系统及存储介质 | |
| CN109621228B (zh) | 放射剂量的计算设备、装置和存储介质 | |
| CN107516330B (zh) | 模型生成方法、图像处理方法及医学成像设备 | |
| US11295158B2 (en) | Image processing apparatus, image processing method, and storage medium for extracting an irradiation field of a radiograph | |
| US20190385279A1 (en) | System and method for image reconstruction | |
| US10271811B2 (en) | Scatter simulation with a radiative transfer equation using direct integral spherical harmonics method for computed tomography | |
| US7262417B2 (en) | Method and system for improved image reconstruction and data collection for compton cameras | |
| JP6746676B2 (ja) | 画像処理装置、画像処理方法、及びプログラム | |
| US10314558B2 (en) | Image processing apparatus, image processing system, image processing method, and recording medium | |
| CN107802280B (zh) | 校正曲线生成方法、投影图像的校正方法、系统及存储介质 | |
| CN108577872B (zh) | 医学图像的去散射方法、系统及存储介质 | |
| US11635392B2 (en) | Radiation imaging apparatus, radiation imaging method, and non-transitory computer-readable storage medium | |
| CN111001097B (zh) | 一种放疗剂量评估系统、装置及存储介质 | |
| US20160349195A1 (en) | Image processing devices, image processing system, image processing method, and non-transitory recording medium | |
| WO2021062885A1 (zh) | Ct成像方法、装置、存储介质及医学成像系统 | |
| JP2017189392A (ja) | 画像処理装置、画像処理方法、画像処理システム及びプログラム。 | |
| US20190192090A1 (en) | Decomposition method and apparatus based on basis material combination | |
| US10231688B2 (en) | Radiation image analysis device, method, and program | |
| JP6744749B2 (ja) | 放射線撮影装置、放射線撮影システム、放射線撮影方法、及びプログラム | |
| US10417795B2 (en) | Iterative reconstruction with system optics modeling using filters | |
| KR20150117417A (ko) | 컴퓨터 단층 촬영 영상 처리 장치 및 이를 이용한 3차원 복원 영상 생성 방법 | |
| JP6803768B2 (ja) | 材料分析装置、材料分析方法及びｘ線ｃｔ装置 | |
| KR101958099B1 (ko) | 양전자 방출 단층 촬영기의 비행시간 정보를 이용한 감쇠보정 방법 | |
| JP2021122674A (ja) | 画像処理装置および方法、プログラム | |
| Jo et al. | Scatter correction method with primary modulator for dual energy digital radiography: a preliminary study |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |