CN111494813B - Modeling method, verification method, device, equipment and storage medium - Google Patents
Modeling method, verification method, device, equipment and storage medium Download PDFInfo
- Publication number
- CN111494813B CN111494813B CN202010316472.3A CN202010316472A CN111494813B CN 111494813 B CN111494813 B CN 111494813B CN 202010316472 A CN202010316472 A CN 202010316472A CN 111494813 B CN111494813 B CN 111494813B
- Authority
- CN
- China
- Prior art keywords
- detector
- particle
- modeling
- photon
- response
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
- A61N5/1071—Monitoring, verifying, controlling systems and methods for verifying the dose delivered by the treatment plan
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/103—Treatment planning systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N5/1048—Monitoring, verifying, controlling systems and methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N2005/1085—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy characterised by the type of particles applied to the patient
- A61N2005/1089—Electrons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N5/00—Radiation therapy
- A61N5/10—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
- A61N2005/1085—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy characterised by the type of particles applied to the patient
- A61N2005/1091—Kilovoltage or orthovoltage range photons
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
Description
技术领域technical field
本发明实施例涉及生物医学信号处理技术领域,尤其涉及一种建模方法、验证方法、装置、设备及存储介质。Embodiments of the present invention relate to the technical field of biomedical signal processing, and in particular, to a modeling method, verification method, device, equipment, and storage medium.
背景技术Background technique
电子射野影像装置(Electronic Portal Imaging Device,EPID)是安装在直线加速器上的平板探测器系统,其可作为基于图像引导的放射治疗(Image-guided radiationtherapy,IGRT)技术中的图像引导装置,用于获取患者图像以辅助医生判断患者摆位是否准确、肿瘤位置或形状是否发生变化等,以减少正常组织接受照射的可能,由此提高了放射治疗的精准性和高效性。Electronic Portal Imaging Device (EPID) is a flat-panel detector system installed on a linear accelerator, which can be used as an image-guided device in image-guided radiation therapy (IGRT) technology. It is used to obtain patient images to assist doctors in judging whether the patient's positioning is accurate, whether the tumor location or shape has changed, etc., so as to reduce the possibility of normal tissue receiving irradiation, thereby improving the accuracy and efficiency of radiation therapy.
近年来,随着探测器技术的快速发展,基于EPID的剂量仪的研究逐渐成为热点。EPID凭借其灵活、方便、快速和高分辨等特点,已成为常用的治疗质量保证工具,其本质上是将图像信息转换为剂量信息,以监测或是重建患者在放射治疗过程中接受的真实三维剂量,而该监测或是重建过程实现的一个重要前提就是EPID响应的精确建模。但是,现有的EPID响应的建模方案多是从单一维度出发,建模精度有待提高。In recent years, with the rapid development of detector technology, research on EPID-based dosimeters has gradually become a hot spot. EPID has become a commonly used treatment quality assurance tool because of its flexibility, convenience, speed and high resolution. It essentially converts image information into dose information to monitor or reconstruct the real three-dimensional images received by patients during radiotherapy. Dose, and an important prerequisite for the monitoring or reconstruction process to be realized is the accurate modeling of the EPID response. However, most of the existing EPID response modeling schemes start from a single dimension, and the modeling accuracy needs to be improved.
发明内容Contents of the invention
本发明实施例提供了一种建模方法、验证方法、装置、设备及存储介质,以实现探测器响应的精确建模的效果。The embodiment of the present invention provides a modeling method, a verification method, a device, a device and a storage medium, so as to realize the effect of accurate modeling of the detector response.
第一方面,本发明实施例提供了一种探测器响应的建模方法,可以包括:In a first aspect, an embodiment of the present invention provides a method for modeling a detector response, which may include:
获取待建模的探测器的结构模型;Obtain the structural model of the detector to be modeled;
获取探测器上的第一粒子的第一粒子信息,其中,第一粒子信息包括类型和入射角度中的至少一个,以及能量,类型包括光子或电子;acquiring first particle information of the first particle on the detector, wherein the first particle information includes at least one of type and incident angle, and energy, and the type includes photon or electron;
根据结构模型和第一粒子信息,模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果。According to the structure model and the information of the first particle, the response value of the first particle on the detector is simulated, and the modeling result is determined according to the response value.
可选的,第一粒子包括第一光子和第一电子,根据结构模型和第一粒子信息,模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果,包括:Optionally, the first particle includes the first photon and the first electron, and the response value of the first particle on the detector is simulated according to the structure model and the information of the first particle, and the modeling result is determined according to the response value, including:
根据结构模型和第一光子的第一光子信息,模拟出第一光子在探测器上的光子响应值,并根据光子响应值确定光子建模结果;Simulating the photon response value of the first photon on the detector according to the structure model and the first photon information of the first photon, and determining the photon modeling result according to the photon response value;
根据结构模型和第一电子的第一电子信息,模拟出第一电子在探测器上的电子响应值,并根据电子响应值确定电子建模结果。According to the structure model and the first electronic information of the first electron, the electronic response value of the first electron on the detector is simulated, and the electronic modeling result is determined according to the electronic response value.
可选的,模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果,可以包括:Optionally, simulating the response value of the first particle on the detector, and determining the modeling result according to the response value may include:
分别模拟出与每个第一粒子信息对应的第一粒子在探测器上的能量沉积,根据各能量沉积以及与各能量沉积分别对应的第一粒子信息,确定建模结果。The energy deposition of the first particles corresponding to each first particle information on the detector is respectively simulated, and the modeling results are determined according to each energy deposition and the first particle information respectively corresponding to each energy deposition.
可选的,探测器可以包括电子射野影像装置,和/或,建模结果可以包括探测器响应曲面。Optionally, the detector may include an electronic portal imaging device, and/or the modeling result may include a detector response surface.
第二方面,本发明实施例还提供了一种放射治疗剂量的验证方法,可包括:In the second aspect, the embodiment of the present invention also provides a method for verifying radiation therapy dose, which may include:
获取探测器上的第二粒子的第二粒子信息以及按照本发明任意实施例所述的探测器响应的建模方法确定的建模结果,根据第二粒子信息和建模结果,预测出第二粒子照射在探测器上的预测图像;Acquiring the second particle information of the second particle on the detector and the modeling result determined according to the modeling method of the detector response described in any embodiment of the present invention, and predicting the second particle information and the modeling result according to the second particle information and the modeling result Predicted images of particles falling on the detector;
获取第二粒子照射在探测器上的实测图像,并对比所述实测图像和所述预测图像间的相似性;acquiring the measured image of the second particle irradiating on the detector, and comparing the similarity between the measured image and the predicted image;
其中,预测图像表征第二粒子在探测器上的预测剂量分布,实测图像表征第二粒子在探测器上的真实剂量分布。Wherein, the predicted image represents the predicted dose distribution of the second particle on the detector, and the measured image represents the real dose distribution of the second particle on the detector.
可选的,建模结果包括光子建模结果和电子建模结果,根据第二粒子信息和建模结果,预测出第二粒子照射在探测器上的预测图像,可以包括:Optionally, the modeling results include photon modeling results and electronic modeling results. According to the second particle information and modeling results, the predicted image of the second particle irradiating on the detector is predicted, which may include:
根据第二粒子信息和光子建模结果,预测出第二粒子照射在探测器上的光子预测图像,并根据第二粒子信息和电子建模结果,预测出第二粒子照射在探测器上的电子预测图像;According to the second particle information and photon modeling results, predict the photon prediction image of the second particle irradiating on the detector, and predict the electrons of the second particle irradiating on the detector according to the second particle information and electron modeling results predicted image;
根据光子预测图像和电子预测图像,生成预测图像。A predicted image is generated based on the photon predicted image and the electron predicted image.
第三方面,本发明实施例还提供了一种探测器响应的建模装置,可包括:In the third aspect, the embodiment of the present invention also provides a detector response modeling device, which may include:
结构模型获取模块,用于获取待建模的探测器的结构模型;A structural model obtaining module, configured to obtain the structural model of the detector to be modeled;
第一粒子信息获取模块,用于获取探测器上的第一粒子的第一粒子信息,其中,第一粒子信息可以包括类型和入射角度中的至少一个,以及能量,类型可以包括光子或电子;A first particle information acquisition module, configured to acquire first particle information of the first particle on the detector, wherein the first particle information may include at least one of type and incident angle, and energy, and the type may include photons or electrons;
建模模块,用于根据结构模型和第一粒子信息,模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果。The modeling module is used to simulate the response value of the first particle on the detector according to the structure model and the information of the first particle, and determine the modeling result according to the response value.
第四方面,本发明实施例还提供了一种放射治疗剂量的验证装置,可包括:In the fourth aspect, the embodiment of the present invention also provides a radiotherapy dose verification device, which may include:
预测图像预测模块,用于获取探测器上的第二粒子的第二粒子信息以及按照本发明任意实施例所述的探测器响应的建模方法确定的建模结果,根据第二粒子信息和建模结果,预测出第二粒子照射在探测器上的预测图像;A predicted image prediction module, configured to obtain the second particle information of the second particle on the detector and the modeling result determined according to the modeling method of the detector response described in any embodiment of the present invention, according to the second particle information and the model predict the predicted image of the second particle irradiated on the detector;
图像对比模块,用于获取第二粒子照射在探测器上的实测图像,并对比实测图像和预测图像间的相似性;The image comparison module is used to obtain the measured image of the second particle irradiated on the detector, and compare the similarity between the measured image and the predicted image;
其中,预测图像表征第二粒子在探测器上的预测剂量分布,实测图像表征第二粒子在探测器上的真实剂量分布。Wherein, the predicted image represents the predicted dose distribution of the second particle on the detector, and the measured image represents the real dose distribution of the second particle on the detector.
第五方面,本发明实施例还提供了一种设备,该设备可以包括:In the fifth aspect, the embodiment of the present invention also provides a device, which may include:
一个或多个处理器;one or more processors;
存储器,用于存储一个或多个程序;memory for storing one or more programs;
当一个或多个程序被一个或多个处理器执行,使得一个或多个处理器实现本发明任意实施例提供的探测器响应的建模方法或是放射治疗剂量的验证方法。When one or more programs are executed by one or more processors, the one or more processors implement the detector response modeling method or the radiotherapy dose verification method provided by any embodiment of the present invention.
第六方面,本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现本发明任意实施例所提供的探测器响应的建模方法或是放射治疗剂量的验证方法。In the sixth aspect, an embodiment of the present invention also provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the method for modeling the detector response provided by any embodiment of the present invention is implemented. Or verification methods for radiation therapy doses.
本发明实施例的技术方案,通过获取待建模的探测器的结构模型以及模拟出的照射在探测器上的第一粒子的第一粒子信息,可根据结构模型和第一粒子信息模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果。上述技术方案,同时考虑到了第一粒子的类型和/或入射角度、以及能量,由此,在探测器响应的建模过程中加入了电子响应建模和/或入射角度建模,实现了探测器响应的多粒子类型多维度的建模效果,提高了探测器响应的建模精度。According to the technical solution of the embodiment of the present invention, by obtaining the structural model of the detector to be modeled and the simulated first particle information of the first particle irradiated on the detector, the first particle information can be simulated according to the structural model and the first particle information. The response value of a particle on the detector, and the modeling result is determined according to the response value. The above technical solution takes into account the type and/or incident angle and energy of the first particle at the same time, thus, electronic response modeling and/or incident angle modeling are added to the modeling process of the detector response to realize the detection The multi-particle type and multi-dimensional modeling effect of the detector response improves the modeling accuracy of the detector response.
附图说明Description of drawings
图1是本发明实施例一中的一种探测器响应的建模方法的流程图;Fig. 1 is a flowchart of a modeling method of a detector response in
图2是本发明实施例一中的一种探测器响应的建模方法中第一粒子的第一粒子信息的应用示意图;2 is a schematic diagram of the application of the first particle information of the first particle in a detector response modeling method in
图3是本发明实施例一中的一种探测器响应的建模方法的第一工作流图;Fig. 3 is a first workflow diagram of a modeling method of a detector response in
图4是本发明实施例一中的一种探测器响应的建模方法的第二工作流图;Fig. 4 is a second workflow diagram of a modeling method of a detector response in
图5是本发明实施例一中的一种探测器响应的建模方法中光子的探测器响应曲面的示意图和电子的探测器响应曲面的示意图;5 is a schematic diagram of a photon detector response surface and a schematic diagram of an electron detector response surface in a detector response modeling method in
图6是本发明实施例二中的一种放射治疗剂量的验证方法的流程图;Fig. 6 is a flow chart of a method for verifying radiation therapy dose in Embodiment 2 of the present invention;
图7是本发明实施例二中的一种放射治疗剂量的验证方法的工作流图;FIG. 7 is a workflow diagram of a method for verifying radiation therapy doses in Embodiment 2 of the present invention;
图8a是本发明实施例二中的一种放射治疗剂量的验证方法的结果示意图;Fig. 8a is a schematic diagram of the results of a verification method of radiation therapy dose in Example 2 of the present invention;
图8b是现有技术中的一种放射治疗剂量的验证方法的结果示意图;Fig. 8b is a schematic diagram of the results of a verification method of radiotherapy dose in the prior art;
图9是本发明实施例三中的一种探测器响应的建模装置的结构框图;9 is a structural block diagram of a detector response modeling device in
图10是本发明实施例四中的一种放射治疗剂量的验证装置的结构框图;Fig. 10 is a structural block diagram of a radiotherapy dose verification device in Embodiment 4 of the present invention;
图11是本发明实施例五中的一种设备的结构示意图。Fig. 11 is a schematic structural diagram of a device in Embodiment 5 of the present invention.
具体实施方式Detailed ways
下面结合附图和实施例对本发明作进一步详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。The present invention will be described in further detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention, but not to limit the present invention. In addition, it should be noted that, for the convenience of description, only some structures related to the present invention are shown in the drawings but not all structures.
实施例一Embodiment one
图1是本发明实施例一中提供的一种探测器响应的建模方法的流程图。本实施例可适用于对探测器的响应进行建模的情况,尤其适用于从多维度出发,对探测器的响应进行精确建模的情况。该方法可以由本发明实施例提供的探测器响应的建模装置来执行,该装置可以由软件和/或硬件的方式实现,该装置可以集成在各种用户终端或服务器上。FIG. 1 is a flow chart of a modeling method for a detector response provided in
参见图1,本发明实施例的方法具体包括如下步骤:Referring to Fig. 1, the method of the embodiment of the present invention specifically includes the following steps:
S110、获取待建模的探测器的结构模型。S110. Acquire a structural model of the detector to be modeled.
其中,探测器可以是任意的二维平板探测器,EPID就是一种常用的二维平板探测器。探测器的结构模型也可称为探测器的物理模型,其可根据探测器的制作材料和几何尺寸预先建立,其是探测器响应的建模过程的重要因素。Wherein, the detector can be any two-dimensional flat panel detector, and EPID is a commonly used two-dimensional flat panel detector. The structural model of the detector can also be referred to as the physical model of the detector, which can be pre-established according to the material and geometric dimensions of the detector, and is an important factor in the modeling process of the detector response.
S120、获取探测器上的第一粒子的第一粒子信息,其中,第一粒子信息包括类型和入射角度中的至少一个,以及能量,类型包括光子或电子。S120. Acquire first particle information of the first particle on the detector, where the first particle information includes at least one of type and incident angle, and energy, and the type includes photon or electron.
其中,第一粒子是模拟出的照射在探测器上的粒子,或是说模拟出的探测器接受到的粒子,该第一粒子可通过笔形束的方式照射在探测器上,该笔形束是大量的能量相同且方向一致的同一类型的第一粒子的集合。Wherein, the first particle is the simulated particle irradiated on the detector, or the simulated particle received by the detector, and the first particle can be irradiated on the detector by means of a pencil beam, and the pencil beam is A collection of a large number of first particles of the same type with the same energy and direction.
第一粒子的第一粒子信息可以是第一粒子的能量和类型,可以是第一粒子的能量和入射角度,还可以是第一粒子的能量、类型和入射角度,等等,在此未做具体限定。第一粒子的类型可以是光子或是电子,也就是说,第一粒子可能是光子也可能是电子,在此未做具体限定。示例性的,以探测器是EPID为例,获取到EPID上的第一粒子的第一粒子信息可以如图2所示,第一粒子信息可以包括由各光子构成的光子束或是由各电子构成的电子束,入射角度α,及能量沉积层中的能量。The first particle information of the first particle can be the energy and type of the first particle, it can be the energy and angle of incidence of the first particle, it can also be the energy, type and angle of incidence of the first particle, etc., which are not done here Specific limits. The type of the first particle may be a photon or an electron, that is, the first particle may be a photon or an electron, which is not specifically limited here. Exemplarily, taking the detector as an EPID as an example, the first particle information of the first particle obtained on the EPID can be shown in Figure 2, and the first particle information can include a photon beam composed of photons or a The formed electron beam, the incident angle α, and the energy in the energy deposition layer.
需要说明的是,这样设置第一粒子信息的原因在于,首先,考虑到直线加速器发射出的第一粒子可能是光子也可能是电子,特别地,当由大量的第一粒子构成的放射治疗射束中的模体厚度较大时,根据模拟计算结果可知,在探测器的响应图像中,电子响应占比可达10%以上。由此,相较于常规的探测器响应的建模方案中只考虑光子响应建模,电子响应建模的加入可在较大程度上提高探测器响应的建模精度。It should be noted that the reason for setting the first particle information in this way is that, first of all, considering that the first particles emitted by the linear accelerator may be photons or electrons, especially when the radiotherapy radiation composed of a large number of first particles When the thickness of the phantom in the beam is relatively large, according to the simulation calculation results, in the response image of the detector, the electronic response accounts for more than 10%. Therefore, compared with the conventional modeling scheme of detector response, which only considers photon response modeling, the addition of electronic response modeling can greatly improve the modeling accuracy of detector response.
其次,常规的探测器响应的建模方案并未区分各第一粒子的入射角度,即假设具有不同的入射角度的第一粒子具有相同的能量效率。但是,根据实验验证可知,当放射治疗射束中存在大厚度和/或复杂的模体时,入射角度建模的加入可明显提高探测器响应的建模精度,即基于能量和入射角度共同计算出的响应值更加贴近真实的响应值。Secondly, the conventional modeling scheme of the detector response does not distinguish the incident angles of the first particles, that is, it is assumed that the first particles with different incident angles have the same energy efficiency. However, according to experimental verification, when there are large thickness and/or complex phantoms in the radiotherapy beam, the addition of incident angle modeling can significantly improve the modeling accuracy of the detector response, that is, based on the joint calculation of energy and incident angle The resulting response value is closer to the real response value.
S130、根据结构模型和第一粒子信息,模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果。S130. Simulate the response value of the first particle on the detector according to the structure model and the information of the first particle, and determine a modeling result according to the response value.
其中,根据已获取的结构模型和第一粒子信息,可以模拟出该第一粒子信息对应的第一粒子在具有该结构模型的探测器上的响应值,进而根据响应值确定出建模结果。实际上,该建模结果可以认为是通过事先模拟第一粒子在探测器内的输运得到的,其可以通过多种形式呈现,如探测器响应曲面等等。Wherein, according to the obtained structure model and first particle information, the response value of the first particle corresponding to the first particle information on the detector with the structure model can be simulated, and then the modeling result can be determined according to the response value. In fact, the modeling result can be considered to be obtained by simulating the transport of the first particle in the detector in advance, which can be presented in various forms, such as the detector response surface and so on.
在此基础上,可选的,该响应值可以是能量沉积,由此,可分别模拟出与每个第一粒子信息对应的第一粒子在探测器上的能量沉积,根据各能量沉积以及与各能量沉积分别对应的第一粒子信息,确定建模结果。由此,一种可选方案,可采用开源的EGSnrc软件中的DOSEXYZnrc模块模拟出探测器的结构模型以得到第一粒子在探测器表面的能量沉积,具体的,可将第一粒子信息作为DOSEXYZnrc模块的输入信息,以使DOSEXYZnrc模块根据输入信息和结构模型计算出与输入信息对应的第一粒子在探测器表面的能量沉积。On this basis, optionally, the response value can be energy deposition, thus, the energy deposition on the detector of the first particle corresponding to each first particle information can be simulated respectively, according to each energy deposition and The first particle information corresponding to each energy deposition determines the modeling result. Therefore, as an alternative, the DOSEXYZnrc module in the open source EGSnrc software can be used to simulate the structural model of the detector to obtain the energy deposition of the first particle on the surface of the detector. Specifically, the information of the first particle can be used as DOSEXYZnrc The input information of the module, so that the DOSEXYZnrc module can calculate the energy deposition of the first particle corresponding to the input information on the detector surface according to the input information and the structure model.
在实际应用中,为提高探测响应的建模结果的鲁棒性,在建模过程中可涉及到各种第一粒子信息,该各种第一粒子信息可包括各种类型和/或各种入射角度,以及各种能量,在此未做具体限定。示例性的,如图3所示,以第一粒子信息包括能量和入射角度为例,在获取到已建立的探测器的结构模型之后,可设置第一粒子的能量和入射角度,并根据结构模型、能量和入射角度计算出第一粒子在探测器上的响应值;进而,判断是否各种能量的第一粒子都计算完毕,若否则更新能量以计算其余能量的第一粒子在探测器上的响应值,若是则判断是否各种入射角度的第一粒子都计算完毕;若否则更新入射角度以计算其余入射角度的第一粒子在探测器上的响应值,若是则根据各种能量和各种入射角度,以及与各种能量和各种入射角度分别对应的响应值,生成探测器响应曲面。In practical applications, in order to improve the robustness of the modeling results of the detection response, various first particle information may be involved in the modeling process, and the various first particle information may include various types and/or various The incident angle and various energies are not specifically limited here. Exemplarily, as shown in Figure 3, taking the first particle information including energy and incident angle as an example, after obtaining the established detector structure model, the energy and incident angle of the first particle can be set, and according to the structure Calculate the response value of the first particle on the detector based on the model, energy and incident angle; then, judge whether the first particle of various energies has been calculated, if not, update the energy to calculate the first particle of the remaining energy on the detector If so, judge whether the first particles at various incident angles have been calculated; if not, update the incident angles to calculate the response values of the first particles at other incident angles on the detector, and if so, according to various energies and Various incident angles, and the response values corresponding to various energies and various incident angles respectively, generate a detector response surface.
在此基础上,可选的,若第一粒子包括第一光子和第一电子,则根据结构模型和第一光子的第一光子信息,可以模拟出第一光子在探测器上的光子响应值,并根据光子响应值确定光子建模结果;且根据结构模型和第一电子的第一电子信息,可以模拟出第一电子在探测器上的电子响应值,并根据电子响应值确定电子建模结果。由此,若第一光子信息包括第一光子的入射角度和能量,第一电子信息包括第一电子的入射角度和能量,则上述建模过程可如图4所示。另外,以光子建模结果为光子的探测器响应曲面且电子建模结果为电子的探测器响应曲面为例,示例性的,如图5所示,根据第一光子的各种入射角度和各种能量可绘制出光子的探测器响应曲面,以及根据第一电子的各种入射角度和各种能量可绘制出电子的探测器响应曲面。On this basis, optionally, if the first particle includes the first photon and the first electron, the photon response value of the first photon on the detector can be simulated according to the structural model and the first photon information of the first photon , and determine the photon modeling result according to the photon response value; and according to the structure model and the first electronic information of the first electron, the electronic response value of the first electron on the detector can be simulated, and the electronic modeling can be determined according to the electronic response value result. Thus, if the first photon information includes the incident angle and energy of the first photon, and the first electronic information includes the incident angle and energy of the first electron, the above modeling process can be shown in FIG. 4 . In addition, taking the photon modeling result as a photon detector response surface and the electronic modeling result as an electron detector response surface as an example, as shown in FIG. 5 , according to various incident angles and various A detector response surface for photons can be drawn for each energy, and a detector response surface for electrons can be drawn for various incident angles and various energies of the first electron.
本发明实施例的技术方案,通过获取待建模的探测器的结构模型以及模拟出的照射在探测器上的第一粒子的第一粒子信息,可根据结构模型和第一粒子信息模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果。上述技术方案,同时考虑到了第一粒子的类型和/或入射角度、以及能量,由此,在探测器响应的建模过程中加入了电子响应建模和/或入射角度建模,实现了探测器响应的多粒子类型多维度的建模效果,提高了探测器响应的建模精度。According to the technical solution of the embodiment of the present invention, by obtaining the structural model of the detector to be modeled and the simulated first particle information of the first particle irradiated on the detector, the first particle information can be simulated according to the structural model and the first particle information. The response value of a particle on the detector, and the modeling result is determined according to the response value. The above technical solution takes into account the type and/or incident angle and energy of the first particle at the same time, thus, electronic response modeling and/or incident angle modeling are added to the modeling process of the detector response to realize the detection The multi-particle type and multi-dimensional modeling effect of the detector response improves the modeling accuracy of the detector response.
需要说明的是,首先,上述技术方案可同时应用于穿透模式的EPID剂量仪和非穿透式的EPID剂量仪。特别地,在穿透式模式中,放射治疗射束中存在患者和/或模体,根据模拟计算结果可知,此时的照射到EPID表面的第一粒子中,电子占比有所提高且入射角度较大,因此,同时考虑到电子响应建模和入射角度建模的技术方案,对穿透模型的精度提升尤为明显。It should be noted that, first of all, the above technical solution can be applied to both the penetrating mode EPID dosimeter and the non-penetrating EPID dosimeter. In particular, in the penetrating mode, there are patients and/or phantoms in the radiotherapy beam. According to the simulation calculation results, the proportion of electrons in the first particles irradiated on the surface of the EPID at this time is increased and the incident The angle is relatively large. Therefore, considering the technical solutions of electronic response modeling and incident angle modeling at the same time, the accuracy improvement of the penetration model is particularly obvious.
其次,上述技术方案可同时应用于前向模式(forward method)的EPID剂量仪和反向模式(backward method)的EPID剂量仪中,且可提高二者的计算精度。具体的,在前向模式中,需要在放射治疗前预测出探测器生成的预测图像,因此,探测器响应的精确建模至关重要。在反向模式中,需要通过探测器实际生成的实测图像反推直线加速器的出射通量,而探测器响应的精确建模是精确重建直线加速器的出射通量的基础。另外,本发明实施例所述的探测器响应的建模方法可应用于多种不同的物理算法,如笔形束算法、卷积算法、蒙特卡罗算法等等,它们都可以应用于前向模式和反向模式的预测图像中。Secondly, the above technical solution can be applied to the EPID dosimeter of the forward mode (forward method) and the EPID dosimeter of the reverse mode (backward method) at the same time, and can improve the calculation accuracy of both. Specifically, in the forward mode, the predicted image generated by the detector needs to be predicted before the radiotherapy, so accurate modeling of the detector response is crucial. In the reverse mode, the output flux of the linac needs to be deduced from the measured image actually generated by the detector, and the accurate modeling of the detector response is the basis for accurately reconstructing the output flux of the linac. In addition, the modeling method of the detector response described in the embodiment of the present invention can be applied to a variety of different physical algorithms, such as pencil beam algorithm, convolution algorithm, Monte Carlo algorithm, etc., and they can all be applied to the forward mode and in the predicted image for the reverse mode.
实施例二Embodiment two
在介绍本发明实施例二之前,先对本发明实施例二的应用场景进行示例性说明:为了实现精准放射治疗,IGRT技术广泛应用于临床中,其可以在治疗前精确定位肿瘤的位置,由此减少了正常组织接受照射的可能性,降低了放射治疗的副作用,且提高了放射治疗的效率。在放射治疗中,为了确保一个放射治疗计划准确地照射到患者身上,使用治疗计划系统(Treatment Plan System,TPS)计算出来的剂量分布通常需要得到实验验证。本发明实施例二详细阐述了对剂量分布进行快速准确验证的具体实现过程,该剂量分布是基于待验证的放射治疗计划模拟出的照射在探测器上的放射治疗射束的剂量分布,由此实现了放射治疗计划的合理性的高效验证。Before introducing the second embodiment of the present invention, the application scenario of the second embodiment of the present invention is exemplified: in order to achieve precise radiotherapy, IGRT technology is widely used in clinics, which can accurately locate the tumor location before treatment, thus It reduces the possibility of normal tissue being irradiated, reduces the side effects of radiation therapy, and improves the efficiency of radiation therapy. In radiotherapy, in order to ensure that a radiotherapy plan is irradiated to patients accurately, the dose distribution calculated using a treatment plan system (Treatment Plan System, TPS) usually needs to be verified experimentally. Embodiment 2 of the present invention elaborates in detail the specific implementation process of quickly and accurately verifying the dose distribution, which is the dose distribution of the radiotherapy beam irradiated on the detector simulated based on the radiotherapy plan to be verified, thus The efficient verification of the rationality of the radiotherapy plan is realized.
图6是本发明实施例二中提供的一种放射治疗剂量的验证方法的流程图。本实施例可适用于对基于待验证的放射治疗计划模拟出的照射在探测器上的放射治疗射束的剂量分布进行验证的情况。该方法可以由本发明实施例提供的放射治疗剂量的验证装置来执行,该装置可以由软件和/或硬件的方式实现,该装置可以集成在各种用户终端或服务器上。Fig. 6 is a flow chart of a radiation therapy dose verification method provided in Embodiment 2 of the present invention. This embodiment is applicable to the case of verifying the dose distribution of the radiotherapy beam irradiated on the detector simulated based on the radiotherapy plan to be verified. The method can be executed by the radiotherapy dose verification device provided by the embodiment of the present invention, the device can be realized by software and/or hardware, and the device can be integrated on various user terminals or servers.
参见图6,本发明实施例的方法具体包括如下步骤:Referring to Figure 6, the method of the embodiment of the present invention specifically includes the following steps:
S210、获取探测器上的第二粒子的第二粒子信息以及按照本发明实施例一所述的探测器响应的建模方法确定的建模结果,根据第二粒子信息和建模结果,预测出第二粒子照射在探测器上的预测图像,其中,预测图像表征第二粒子在探测器上的预测剂量分布。S210. Acquire the second particle information of the second particle on the detector and the modeling result determined according to the detector response modeling method described in
其中,第二粒子是根据待验证的放射治疗计划模拟出的照射在探测器上的粒子,该第二粒子的第二粒子信息与本发明实施例一中所述的第一粒子信息相一致,比如,若第一粒子信息包括能量、类型和入射角度,则第二粒子信息也包括能量、类型和入射角度。由此,根据第二粒子信息和已生成的建模结果,则可模拟出第二粒子照射在探测器上的预测图像。Wherein, the second particle is the particle irradiated on the detector simulated according to the radiotherapy plan to be verified, and the second particle information of the second particle is consistent with the first particle information described in
具体的,可选的,可采用蒙特卡罗算法计算出虚拟的放射治疗射束在探测器上的能注量分布图,每个放射治疗射束是由大量的能量相同且方向一致的同一类型的第二粒子构成的:根据能注量分布图和探测器的建模结果计算出预测图像,该预测图像可以表征第二粒子在探测器上的预测剂量分布,即为虚拟的放射治疗射束在探测器上的剂量分布。Specifically, optionally, the Monte Carlo algorithm can be used to calculate the energy fluence distribution map of the virtual radiotherapy beam on the detector. Each radiotherapy beam is composed of a large number of the same type Composed of the second particles: the predicted image is calculated according to the energy fluence distribution map and the modeling results of the detector, and the predicted image can represent the predicted dose distribution of the second particle on the detector, which is a virtual radiotherapy beam Dose distribution on the detector.
在此基础上,若建模结果包括光子建模结果和电子建模结果,则根据第二粒子信息和光子建模结果,可以预测出第二粒子照射在探测器上的光子预测图像;根据第二粒子信息和电子建模结果,可以预测出第二粒子照射在探测器上的电子预测图像;进而,根据光子预测图像和电子预测图像,可以生成预测图像。也就是说,预测图像是光子建模结果和电子建模结果的综合预测结果。On this basis, if the modeling results include photon modeling results and electronic modeling results, then according to the second particle information and photon modeling results, the photon prediction image of the second particle irradiating on the detector can be predicted; according to the The information of the two particles and the result of electron modeling can predict the electron prediction image of the second particle irradiated on the detector; furthermore, according to the photon prediction image and the electron prediction image, the prediction image can be generated. That is, the predicted image is a comprehensive predicted result of photon modeling results and electronic modeling results.
S220、获取第二粒子照射在探测器上的实测图像,并对比实测图像和预测图像间的差异性,实测图像表征第二粒子在探测器上的真实剂量分布。S220. Acquire a measured image of the second particle irradiating on the detector, and compare the difference between the measured image and the predicted image. The measured image represents a real dose distribution of the second particle on the detector.
其中,直线加速器根据待验证的放射治疗计划发射出具有一定剂量分布的放射治疗射束,每个放射治疗射束由大量的能量相同且方向一致的同一类型的第二粒子构成,该放射治疗射束照射在探测器的表面时,探测器根据其内部的感光元件对照射到其上的放射治疗射束中的第二粒子的感应信号进行成像,得到探测器实际检测到的实测图像,因此,该实测图像可以表征基于待验证的放射治疗计划发射出的放射治疗射束在探测器上的真实剂量分布,即为实际发射出的第二粒子在探测器上的真实剂量分布。Among them, the linear accelerator emits radiotherapy beams with a certain dose distribution according to the radiotherapy plan to be verified, and each radiotherapy beam is composed of a large number of second particles of the same type with the same energy and the same direction. When the beam is irradiated on the surface of the detector, the detector will image the induction signal of the second particle in the radiotherapy beam irradiated on it according to its internal photosensitive element, and obtain the measured image actually detected by the detector. Therefore, The measured image can represent the real dose distribution of the radiation therapy beams emitted on the detector based on the radiation therapy plan to be verified, that is, the real dose distribution of the second particles actually emitted on the detector.
需要说明的是,通常情况下,因探测器处于高速采集状态,针对放射治疗计划中的每个放射治疗射束,探测器均可检测到多幅实测图像,由此,可以将该多幅实测图像相叠加,如将各幅实测图像相对应位置像素点的灰度值相叠加,并将叠加结果作为一幅实测图像。进一步,若放射治疗计划中具有多个放射治疗射束,则可对每个放射治疗射束下产生的多幅实测图像分别进行处理,由此,每个放射治疗射束可对应一幅实测图像。It should be noted that, usually, because the detector is in a high-speed acquisition state, the detector can detect multiple actual images for each radiotherapy beam in the radiotherapy plan. Therefore, the multiple actual images can be The images are superimposed, such as superimposing the gray values of the corresponding pixel points of each measured image, and using the superimposed result as a measured image. Further, if there are multiple radiation therapy beams in the radiation therapy plan, multiple measured images generated under each radiation therapy beam can be processed separately, so that each radiation therapy beam can correspond to a measured image .
在获取到第二粒子在探测器上模拟出的预测图像和真实采集的实测图像后,可对二者的差异性进行对比,该差异性可以表征实测图像和预测图像间的吻合程度,根据该吻合程度可验证相应的放射治疗计划是否合理,从而辅助医生对放射治疗计划进行调整。在此基础上,可选的,这一对比过程可通过伽马通过评价算法实现。具体的,通过伽马函数计算图像伽马值,该伽马函数如下所示:After obtaining the predicted image simulated by the second particle on the detector and the measured image actually collected, the difference between the two can be compared, and the difference can represent the degree of agreement between the measured image and the predicted image. According to the The degree of coincidence can verify whether the corresponding radiation therapy plan is reasonable, thereby assisting doctors to adjust the radiation therapy plan. On this basis, optionally, this comparison process can be realized by a gamma pass evaluation algorithm. Specifically, the gamma value of the image is calculated through the gamma function, and the gamma function is as follows:
其中,分别为计算点(即,预测图像中的像素点)和参考点(即,实测图像中的像素点),为计算点和参考点的剂量差,为计算点和参考点的距离,ΔD和Δd分别为预先设置的剂量容差和距离容差。由此,根据伽马函数可以计算参考点的伽马值:in, are calculation points (that is, pixels in the predicted image) and reference points (that is, pixels in the measured image), respectively, is the dose difference between the calculation point and the reference point, For the distance between the calculation point and the reference point, ΔD and Δd are the preset dose tolerance and distance tolerance, respectively. Thus, the gamma value of the reference point can be calculated according to the gamma function:
相应的伽马评价准则是,若则通过,若则未通过。这样一来,计算通过点占总计算点百分比即可获得实测图像的伽马通过率,根据伽马通过率即可确定实测图像和预测图像间的差异性。The corresponding gamma evaluation criterion is, if pass if is not passed. In this way, the gamma pass rate of the measured image can be obtained by calculating the percentage of pass points in the total calculation points, and the difference between the measured image and the predicted image can be determined according to the gamma pass rate.
本发明实施例的技术方案,通过获取到的探测器上的第二粒子的第二粒子信息以及该探测器的建模结果,可以模拟出第二粒子照射在探测器上的预测图像;进而,通过获取到的第二粒子照射在探测器上的实测图像,可以对比出实测图像和预测图像间的相似性,该相似性可用于验证与第二粒子对应的放射治疗计划的合理性。上述技术方案,因利用了已生成的探测器的建模结果,减小了预测图像的计算耗时且提高了预测图像的计算精度,由此实现了放射治疗计划中的放射治疗剂量的快速且精确验证的效果。According to the technical solution of the embodiment of the present invention, by obtaining the second particle information of the second particle on the detector and the modeling result of the detector, a predicted image of the second particle irradiating on the detector can be simulated; furthermore, Through the obtained measured image of the second particle irradiating on the detector, the similarity between the measured image and the predicted image can be compared, and the similarity can be used to verify the rationality of the radiotherapy plan corresponding to the second particle. The above technical solution reduces the time-consuming calculation of the predicted image and improves the calculation accuracy of the predicted image due to the use of the generated detector modeling results, thus realizing the rapid and accurate calculation of the radiation treatment dose in the radiation treatment plan. The effect of precise verification.
为了更好地理解上述步骤的具体实现过程,下面结合具体示例,对本实施例的放射治疗剂量的验证方法进行示例性的说明。In order to better understand the specific implementation process of the above steps, the method for verifying the radiotherapy dose in this embodiment will be exemplarily described below in conjunction with specific examples.
示例性的,如图7所示,以一个头颈部动态调强计划(dIMRT)为例,该放射治疗计划共有9个放射治疗射束,其中放置固体水(一种模体),该固体水的尺寸为30cm×30cm×15cm,固体水中心位于机器等中心。根据光源模型和模体信息,利用蒙特卡罗算法可计算出探测器表面的第二粒子的第二粒子信息,该第二粒子信息包括类型、能量和入射角度。需要说明的是,光源模型中的光源可以认为是大量的第二粒子的集合,从光源模型中提取出大量的第二粒子在照射到模体后,模体会改变这些第二粒子的第二粒子信息,因此,根据光源模型和模体信息可以计算出探测器表面上的第二粒子的第二粒子信息。进一步,根据第二粒子信息和已生成的探测器响应的建模结果计算出预测图像,并将预测图像和实测图像进行对比,得到实测图像的伽马通过率。Exemplarily, as shown in Figure 7, taking a head and neck dynamic intensity-modulated plan (dIMRT) as an example, the radiotherapy plan has 9 radiotherapy beams in total, in which solid water (a phantom) is placed, and the solid The dimensions of the water are 30cm x 30cm x 15cm with the center of solid water at the isocenter of the machine. According to the light source model and the phantom information, the second particle information of the second particle on the detector surface can be calculated by using the Monte Carlo algorithm, and the second particle information includes type, energy and incident angle. It should be noted that the light source in the light source model can be considered as a collection of a large number of second particles. After a large number of second particles are extracted from the light source model and irradiated to the phantom, the phantom will change the second particles of these second particles Therefore, the second particle information of the second particle on the detector surface can be calculated according to the light source model and phantom information. Further, the predicted image is calculated according to the second particle information and the generated modeling results of the detector response, and the predicted image is compared with the measured image to obtain the gamma pass rate of the measured image.
经试验验证,若探测器响应的建模过程中同时考虑到光子的能量和入射角度以及电子的能量和入射角度,则第一个放射治疗射束的预测图像、实测图像和伽马值图像如图8a所示,其中,(a)是预测图像,(b)是实测图像,(c)是伽马值图像;若探测器响应的建模过程中只考虑到光子的能量,则第一个放射治疗射束的预测图像、实测图像和伽马值图像如图8b所示,其中,(a)是预测图像,(b)是实测图像,(c)是伽马值图像,对比图8a和图8b可知,在图8b中,相较于实测图像,预测图像的数值偏低,伽马通过率明显下降。It has been verified by experiments that if the energy and incident angle of photons and the energy and incident angle of electrons are considered in the modeling process of the detector response, the predicted image, measured image and gamma value image of the first radiotherapy beam are as follows: As shown in Figure 8a, (a) is the predicted image, (b) is the measured image, and (c) is the gamma value image; if only the energy of the photon is considered in the modeling process of the detector response, the first The predicted image, measured image and gamma value image of the radiotherapy beam are shown in Figure 8b, where (a) is the predicted image, (b) is the measured image, and (c) is the gamma value image, compared with Figure 8a and It can be seen from Figure 8b that in Figure 8b, compared with the measured image, the value of the predicted image is low, and the gamma pass rate drops significantly.
上述技术方案是一种放射治疗前的治疗质量保证方案,其采用前向模式,在治疗患者前实施放射治疗计划,并通过探测器记录实测图像。需要说明的是,该技术方案是用于验证放射治疗剂量,并非用于治疗,因此在获取实测图像时,放射治疗射束中无任何模体。而且,基于蒙特卡罗算法,TPS可以精确预测出与实测图像对应的预测图像。进而,通过伽马通过率可定量比较实测图像和预测图像间的差异性,由此验证出直线加速器的出束状态是否达到预期。The technical solution above is a treatment quality assurance solution before radiotherapy, which adopts a forward mode, implements a radiotherapy plan before treating patients, and records actual measured images through a detector. It should be noted that this technical solution is used to verify the dose of radiation therapy, not for treatment, so when acquiring the measured image, there is no phantom in the radiation therapy beam. Moreover, based on the Monte Carlo algorithm, TPS can accurately predict the predicted image corresponding to the measured image. Furthermore, the gamma pass rate can be used to quantitatively compare the difference between the measured image and the predicted image, thereby verifying whether the beam output state of the linac meets expectations.
另外,探测器响应的建模结果还可应用于探测器成像相关领域,通过考虑第二粒子信息完成图像校正。例如,在EPID成像相关领域中,辅助完成图像引导放疗也是EPID的一项重要功能,而要获取高质量的锥形束投照电子计算机断层扫描(Cone Beam ComputerTomography,CBCT)图像,通常需要对EPID图像实施一系列的图像校正,如能谱硬化校正、散射校正等等。这一图像校正方案的通常做法是先建立探测器响应的建模结果,通过模拟和定量计算能谱硬化和散射等因素对图像质量的影响,再进一步对实测图像完成校正。In addition, the modeling results of the detector response can also be applied to the related fields of detector imaging, and the image correction can be completed by considering the second particle information. For example, in the field of EPID imaging, assisting in the completion of image-guided radiotherapy is also an important function of EPID, and to obtain high-quality cone beam computer tomography (CBCT) images, it usually requires EPID The image implements a series of image corrections, such as spectral hardening correction, scatter correction and so on. The usual method of this image correction scheme is to establish the modeling results of the detector response first, and then to further complete the correction of the measured image by simulating and quantitatively calculating the influence of factors such as energy spectrum hardening and scattering on the image quality.
实施例三Embodiment Three
图9为本发明实施例三提供的探测器响应的建模装置的结构框图,该装置用于执行上述任意实施例所提供的探测器响应的建模方法。该装置与上述各实施例的探测器响应的建模方法属于同一个发明构思,在探测器响应的建模装置的实施例中未详尽描述的细节内容,可以参考上述探测器响应的建模方法的实施例。参见图9,该装置具体可包括:结构模型获取模块310、第一粒子信息获取模块320和建模模块330。FIG. 9 is a structural block diagram of a detector response modeling device provided in
其中,结构模型获取模块310,用于获取待建模的探测器的结构模型;Wherein, the structure
第一粒子信息获取模块320,用于获取探测器上的第一粒子的第一粒子信息,其中,第一粒子信息可以包括类型和入射角度中的至少一个,以及能量,类型可以包括光子或电子;The first particle
建模模块330,用于根据结构模型和第一粒子信息,模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果。The
可选的,建模模块330,具体可以包括:Optionally, the
光子建模单元,用于第一粒子包括第一光子,根据结构模型和第一光子的第一光子信息,模拟出第一光子在探测器上的光子响应值,并根据光子响应值确定光子建模结果;The photon modeling unit is used for the first particle including the first photon, simulates the photon response value of the first photon on the detector according to the structure model and the first photon information of the first photon, and determines the photon modeling unit according to the photon response value model result;
电子建模单元,用于第一粒子还包括第一电子,根据结构模型和第一电子的第一电子信息,模拟出第一电子在探测器上的电子响应值,并根据电子响应值确定电子建模结果。The electronic modeling unit is used for the first particle and also includes the first electron, and simulates the electronic response value of the first electron on the detector according to the structure model and the first electronic information of the first electron, and determines the electron response value according to the electronic response value. modeling results.
可选的,建模模块330,具体可以包括:Optionally, the
能量沉积模拟单元,用于分别模拟出与每个第一粒子信息对应的第一粒子在探测器上的能量沉积,根据各能量沉积以及与各能量沉积分别对应的第一粒子信息,确定建模结果。The energy deposition simulation unit is used to separately simulate the energy deposition of the first particles corresponding to each first particle information on the detector, and determine the modeling according to each energy deposition and the first particle information corresponding to each energy deposition result.
可选的,探测器可以包括电子射野影像装置,和/或,建模结果可以包括探测器响应曲面。Optionally, the detector may include an electronic portal imaging device, and/or the modeling result may include a detector response surface.
本发明实施例三提供的探测器响应的建模装置,通过结构模型获取模块、第一粒子信息获取模块和建模模块相互配合,获取待建模的探测器的结构模型以及模拟出的照射在探测器上的第一粒子的第一粒子信息,可根据结构模型和第一粒子信息模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果。上述装置,同时考虑到了第一粒子的类型和/或入射角度、以及能量,由此,在探测器响应的建模过程中加入了电子响应建模和/或入射角度建模,实现了探测器响应的多粒子类型多维度的建模效果,提高了探测器响应的建模精度。The detector response modeling device provided in
本发明实施例所提供的探测器响应的建模装置可执行本发明任意实施例所提供的探测器响应的建模方法,具备执行方法相应的功能模块和有益效果。The detector response modeling device provided in the embodiments of the present invention can execute the detector response modeling method provided in any embodiment of the present invention, and has corresponding functional modules and beneficial effects for executing the method.
值得注意的是,上述探测器响应的建模装置的实施例中,所包括的各个单元和模块只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。It is worth noting that, in the above-mentioned embodiments of the detector response modeling device, the included units and modules are only divided according to functional logic, but are not limited to the above-mentioned divisions, as long as the corresponding functions can be realized ; In addition, the specific names of each functional unit are only for the convenience of distinguishing each other, and are not used to limit the protection scope of the present invention.
实施例四Embodiment Four
图10为本发明实施例四提供的放射治疗剂量的验证装置的结构框图,该装置用于执行上述任意实施例所提供的放射治疗剂量的验证方法。该装置与上述各实施例的放射治疗剂量的验证方法属于同一个发明构思,在放射治疗剂量的验证装置的实施例中未详尽描述的细节内容,可以参考上述放射治疗剂量的验证方法的实施例。参见图10,该装置具体可包括:预测图像预测模块410和图像对比模块420。Fig. 10 is a structural block diagram of a radiotherapy dose verification device provided in Embodiment 4 of the present invention, and the device is used to implement the radiotherapy dose verification method provided in any of the above-mentioned embodiments. This device belongs to the same inventive concept as the verification method of radiotherapy dose in the above-mentioned embodiments. For the details not described in detail in the embodiment of the verification device of radiotherapy dose, you can refer to the above-mentioned embodiment of the verification method of radiotherapy dose . Referring to FIG. 10 , the device may specifically include: a predicted
其中,预测图像预测模块410,用于获取探测器上的第二粒子的第二粒子信息以及按照本发明任意实施例所述的探测器响应的建模方法确定的建模结果,根据第二粒子信息和建模结果,预测出第二粒子照射在探测器上的预测图像;Among them, the predicted
图像对比模块420,用于获取第二粒子照射在探测器上的实测图像,并对比实测图像和预测图像间的相似性;An
其中,预测图像表征第二粒子在探测器上的预测剂量分布,实测图像表征第二粒子在探测器上的真实剂量分布。Wherein, the predicted image represents the predicted dose distribution of the second particle on the detector, and the measured image represents the real dose distribution of the second particle on the detector.
可选的,预测图像预测模块410,具体可以包括:Optionally, the predicted
光子预测图像预测单元,用于建模结果包括光子建模结果,根据第二粒子信息和光子建模结果,预测出第二粒子照射在探测器上的光子预测图像;The photon prediction image prediction unit is used to predict the photon prediction image of the second particle irradiated on the detector according to the second particle information and the photon modeling result for the modeling result including the photon modeling result;
电子预测图像预测单元,用于建模结果还包括电子建模结果,根据第二粒子信息和电子建模结果,预测出第二粒子照射在探测器上的电子预测图像;The electron predictive image prediction unit is used for modeling results and also includes electronic modeling results, and predicts the electron predictive image of the second particle irradiating on the detector according to the second particle information and the electronic modeling results;
预测图像生成单元,用于根据光子预测图像和电子预测图像生成预测图像。A predictive image generation unit is used to generate a predictive image according to the photon predictive image and the electron predictive image.
本发明实施例四提供的放射治疗剂量的验证装置,通过预测图像预测模块获取到的探测器上的第二粒子的第二粒子信息以及该探测器的建模结果,可以模拟出第二粒子照射在探测器上的预测图像;进而,图像对比模块获取到的第二粒子照射在探测器上的实测图像,可以对比出实测图像和预测图像间的相似性,该相似性可用于验证与第二粒子对应的放射治疗计划的合理性。上述装置,因利用了已生成的探测器的建模结果,减小了预测图像的计算耗时且提高了预测图像的计算精度,由此实现了放射治疗计划中的放射治疗剂量的快速且精确验证的效果。The radiotherapy dose verification device provided in Embodiment 4 of the present invention can simulate the second particle irradiation by predicting the second particle information of the second particle on the detector acquired by the image prediction module and the modeling result of the detector. The predicted image on the detector; furthermore, the measured image of the second particle irradiated on the detector obtained by the image comparison module can compare the similarity between the measured image and the predicted image, and the similarity can be used to verify the similarity with the second particle. Particles correspond to the rationality of the radiotherapy plan. The above device reduces the time-consuming calculation of the predicted image and improves the calculation accuracy of the predicted image due to the use of the generated detector modeling results, thus realizing the rapid and accurate calculation of the radiation treatment dose in the radiation treatment plan Verify the effect.
本发明实施例所提供的放射治疗剂量的验证装置可执行本发明任意实施例所提供的放射治疗剂量的验证方法,具备执行方法相应的功能模块和有益效果。The radiotherapy dose verification device provided in the embodiments of the present invention can execute the radiotherapy dose verification method provided in any embodiment of the present invention, and has corresponding functional modules and beneficial effects for executing the method.
值得注意的是,上述放射治疗剂量的验证装置的实施例中,所包括的各个单元和模块只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。It is worth noting that, in the embodiment of the radiation therapy dose verification device described above, the units and modules included are only divided according to functional logic, but are not limited to the above division, as long as the corresponding functions can be realized; In addition, the specific names of the functional units are only for the convenience of distinguishing each other, and are not used to limit the protection scope of the present invention.
实施例五Embodiment five
图11为本发明实施例五提供的一种设备的结构示意图,如图11所示,该设备包括存储器510、处理器520、输入装置530和输出装置540。设备中的处理器520的数量可以是一个或多个,图11中以一个处理器520为例;设备中的存储器510、处理器520、输入装置530和输出装置540可以通过总线或其它方式连接,图11中以通过总线550连接为例。FIG. 11 is a schematic structural diagram of a device provided by Embodiment 5 of the present invention. As shown in FIG. 11 , the device includes a
存储器510作为一种计算机可读存储介质,可用于存储软件程序、计算机可执行程序以及模块,如本发明实施例中的探测器响应的建模方法对应的程序指令/模块(例如,探测器响应的建模装置中的结构模型获取模块310、第一粒子信息获取模块320和建模模块330),再如本发明实施例中的放射治疗剂量的验证方法对应的程序指令/模块(例如,放射治疗剂量的验证装置中的预测图像预测模块410和图像对比模块420)。处理器520通过运行存储在存储器510中的软件程序、指令以及模块,从而执行设备的各种功能应用以及数据处理,即实现上述的探测器响应的建模方法或是上述的放射治疗剂量的验证方法。The
存储器510可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序;存储数据区可存储根据设备的使用所创建的数据等。此外,存储器510可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中,存储器510可进一步包括相对于处理器520远程设置的存储器,这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。The
输入装置530可用于接收输入的数字或字符信息,以及产生与装置的用户设置以及功能控制有关的键信号输入。输出装置540可包括显示屏等显示设备。The
实施例六Embodiment six
本发明实施例六提供了一种包含计算机可执行指令的存储介质,所述计算机可执行指令在由计算机处理器执行时用于执行一种探测器响应的建模方法,该方法包括:Embodiment 6 of the present invention provides a storage medium containing computer-executable instructions, and the computer-executable instructions are used to perform a modeling method of a detector response when executed by a computer processor. The method includes:
获取待建模的探测器的结构模型;Obtain the structural model of the detector to be modeled;
获取探测器上的第一粒子的第一粒子信息,其中,第一粒子信息包括类型和入射角度中的至少一个,以及能量,类型包括光子或电子;acquiring first particle information of the first particle on the detector, wherein the first particle information includes at least one of type and incident angle, and energy, and the type includes photon or electron;
根据结构模型和第一粒子信息,模拟出第一粒子在探测器上的响应值,并根据响应值确定建模结果。According to the structure model and the information of the first particle, the response value of the first particle on the detector is simulated, and the modeling result is determined according to the response value.
当然,本发明实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的方法操作,还可以执行本发明任意实施例所提供的探测器响应的建模方法中的相关操作。Certainly, a storage medium containing computer-executable instructions provided by an embodiment of the present invention, the computer-executable instructions are not limited to the method operations described above, and may also implement the detector response suggestions provided by any embodiment of the present invention. Related operations in the modulo method.
实施例七Embodiment seven
本发明实施例七提供一种包含计算机可执行指令的存储介质,所述计算机可执行指令在由计算机处理器执行时用于执行一种放射治疗剂量的验证方法,该方法包括:Embodiment 7 of the present invention provides a storage medium containing computer-executable instructions, and the computer-executable instructions are used to perform a radiation therapy dose verification method when executed by a computer processor. The method includes:
获取探测器上的第二粒子的第二粒子信息以及按照本发明任意实施例所述的探测器响应的建模方法确定的建模结果,根据第二粒子信息和建模结果,预测出第二粒子照射在探测器上的预测图像;Acquiring the second particle information of the second particle on the detector and the modeling result determined according to the modeling method of the detector response described in any embodiment of the present invention, and predicting the second particle information and the modeling result according to the second particle information and the modeling result Predicted images of particles falling on the detector;
获取第二粒子照射在探测器上的实测图像,并对比所述实测图像和所述预测图像间的相似性;acquiring the measured image of the second particle irradiating on the detector, and comparing the similarity between the measured image and the predicted image;
其中,预测图像表征第二粒子在探测器上的预测剂量分布,实测图像表征第二粒子在探测器上的真实剂量分布。Wherein, the predicted image represents the predicted dose distribution of the second particle on the detector, and the measured image represents the real dose distribution of the second particle on the detector.
当然,本发明实施例所提供的一种包含计算机可执行指令的存储介质,其计算机可执行指令不限于如上所述的方法操作,还可以执行本发明任意实施例所提供的放射治疗剂量的验证方法中的相关操作。Certainly, a storage medium containing computer-executable instructions provided by an embodiment of the present invention, the computer-executable instructions are not limited to the method operations described above, and can also perform the verification of radiation therapy dose provided by any embodiment of the present invention Related operations in the method.
通过以上关于实施方式的描述,所属领域的技术人员可以清楚地了解到,本发明可借助软件及必需的通用硬件来实现,当然也可以通过硬件实现,但很多情况下前者是更佳的实施方式。依据这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。Through the above description about the implementation mode, those skilled in the art can clearly understand that the present invention can be realized by means of software and necessary general-purpose hardware, and of course it can also be realized by hardware, but in many cases the former is a better implementation mode . Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as a floppy disk of a computer , read-only memory (Read-Only Memory, ROM), random access memory (Random Access Memory, RAM), flash memory (FLASH), hard disk or optical disc, etc., including several instructions to make a computer device (which can be a personal computer, A server, or a network device, etc.) executes the methods described in various embodiments of the present invention.
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。Note that the above are only preferred embodiments of the present invention and applied technical principles. Those skilled in the art will understand that the present invention is not limited to the specific embodiments described herein, and that various obvious changes, readjustments and substitutions can be made by those skilled in the art without departing from the protection scope of the present invention. Therefore, although the present invention has been described in detail through the above embodiments, the present invention is not limited to the above embodiments, and can also include more other equivalent embodiments without departing from the concept of the present invention, and the present invention The scope is determined by the scope of the appended claims.
Claims (9)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010316472.3A CN111494813B (en) | 2020-04-21 | 2020-04-21 | Modeling method, verification method, device, equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010316472.3A CN111494813B (en) | 2020-04-21 | 2020-04-21 | Modeling method, verification method, device, equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111494813A CN111494813A (en) | 2020-08-07 |
| CN111494813B true CN111494813B (en) | 2022-11-08 |
Family
ID=71848334
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010316472.3A Active CN111494813B (en) | 2020-04-21 | 2020-04-21 | Modeling method, verification method, device, equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111494813B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022143333A1 (en) * | 2020-12-31 | 2022-07-07 | 中硼(厦门)医疗器械有限公司 | Radiation irradiation system and control method therefor |
| CN112581475B (en) * | 2021-02-25 | 2021-05-25 | 四川大学华西医院 | A method for predicting gamma pass rate for radiotherapy planning |
| CN113570705B (en) * | 2021-07-28 | 2024-04-30 | 广州瑞多思医疗科技有限公司 | Three-dimensional dose reconstruction method, device, computer equipment and storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5821541A (en) * | 1996-02-02 | 1998-10-13 | Tuemer; Tuemay O. | Method and apparatus for radiation detection |
| EP1917999A1 (en) * | 2006-11-03 | 2008-05-07 | Ion Beam Applications S.A. | Method and device for IMRT verification |
| EP2116277A1 (en) * | 2008-05-06 | 2009-11-11 | Ion Beam Applications S.A. | Device and method for particle therapy monitoring and verification |
Family Cites Families (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2356956A1 (en) * | 1976-02-26 | 1978-01-27 | Drukier Andnej | 3-D image reconstruction of physical magnitudes inside bodies - creating position-electron pairs and reconstructing annihilation photon emission point from detection coordinates |
| JPS62278480A (en) * | 1986-05-28 | 1987-12-03 | Toshiba Corp | Radiation detecting circuit |
| US6345114B1 (en) * | 1995-06-14 | 2002-02-05 | Wisconsin Alumni Research Foundation | Method and apparatus for calibration of radiation therapy equipment and verification of radiation treatment |
| US6448560B1 (en) * | 1996-02-02 | 2002-09-10 | Tumay O. Tumer | Method and apparatus for gamma ray detection |
| US6236050B1 (en) * | 1996-02-02 | 2001-05-22 | TüMER TüMAY O. | Method and apparatus for radiation detection |
| US6563906B2 (en) * | 2000-08-28 | 2003-05-13 | University Of New Brunswick | X-ray compton scattering density measurement at a point within an object |
| CN1481756A (en) * | 2002-08-14 | 2004-03-17 | ��ʽ���綫֥ | Concentrated radiation therapy equipment |
| CN101523171A (en) * | 2005-07-11 | 2009-09-02 | 瑞沃瑞公司 | Method and system for non-destructive distribution profiling of an element in a film |
| DE102009002273A1 (en) * | 2009-04-07 | 2010-10-14 | Gsi Helmholtzzentrum Für Schwerionenforschung Gmbh | Method for operating a system for irradiating a particle beam onto a body and device |
| CN103135125A (en) * | 2011-12-01 | 2013-06-05 | 中国辐射防护研究院 | Energy spectrum analytical method of X-ray filtering reference radiation |
| US9057788B2 (en) * | 2012-02-03 | 2015-06-16 | Analogic Corporatiom | Photon counting-based virtual detector |
| US9020092B2 (en) * | 2013-02-19 | 2015-04-28 | Kabushiki Kaisha Toshiba | Apparatus and method for angular response calibration of photon-counting detectors in sparse spectral computed tomography imaging |
| CN105324684B (en) * | 2013-06-13 | 2019-11-05 | 皇家飞利浦有限公司 | Detectors for Radiation Therapy Disposition Guidance and Validation |
| CN104436450B (en) * | 2013-09-25 | 2018-07-31 | 苏州雷泰医疗科技有限公司 | Radiotherapy unit light limiting means and radiotherapy unit |
| CN105468883B (en) * | 2014-08-20 | 2018-11-06 | 上海联影医疗科技有限公司 | The source model method for building up and device of electron linear accelerator |
| US10114129B2 (en) * | 2016-01-28 | 2018-10-30 | The Research Foundation For The State University Of New York | Semiconductor detector for x-ray single-photon detection |
| CZ29250U1 (en) * | 2016-01-29 | 2016-03-08 | Advacam S.R.O. | Ionizing radiation laminated pixel detector |
| CN105866151B (en) * | 2016-04-22 | 2018-07-24 | 西北核技术研究所 | A kind of crystal rocking curve measurement method based on energy-resolved detector |
| CN107688194B (en) * | 2016-08-03 | 2020-12-29 | 清华大学 | Data processing method and data processing device |
| CN106388841B (en) * | 2016-09-23 | 2019-05-21 | 东软医疗系统股份有限公司 | A kind of image reconstruction method, apparatus and equipment |
| CN108926356B (en) * | 2018-05-31 | 2022-03-25 | 上海联影医疗科技股份有限公司 | Detector photon arrival information determination method, device and medical device |
| CN108615250B (en) * | 2018-05-31 | 2022-04-26 | 上海联影医疗科技股份有限公司 | Image reconstruction method, device, system and computer readable storage medium |
| CN109683189A (en) * | 2019-01-15 | 2019-04-26 | 南昌大学 | A kind of cherenkov radiation detector and detection method |
| CN110025329B (en) * | 2019-04-24 | 2022-12-27 | 上海联影医疗科技股份有限公司 | Coincidence counting chord graph generation method and device, computer equipment and storage medium |
-
2020
- 2020-04-21 CN CN202010316472.3A patent/CN111494813B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5821541A (en) * | 1996-02-02 | 1998-10-13 | Tuemer; Tuemay O. | Method and apparatus for radiation detection |
| EP1917999A1 (en) * | 2006-11-03 | 2008-05-07 | Ion Beam Applications S.A. | Method and device for IMRT verification |
| EP2116277A1 (en) * | 2008-05-06 | 2009-11-11 | Ion Beam Applications S.A. | Device and method for particle therapy monitoring and verification |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111494813A (en) | 2020-08-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105854191B (en) | 3-dimensional dose verification system and verification method in a kind of radiotherapy | |
| CN111408072B (en) | Device for evaluating radiation dose delivery and system for verifying radiation dose delivery | |
| US8515011B2 (en) | System and method for dose verification radiotherapy | |
| CN104548372A (en) | Radiotherapy planning method and device, radiotherapy dose determining method and device and radiotherapy quality guaranteeing method and device | |
| CN103127623B (en) | Method of online authentication of accelerator out-beam accuracy in radiation therapy | |
| JP6892649B2 (en) | A device with a radiation therapy planning function and how to operate the device | |
| CN111494813B (en) | Modeling method, verification method, device, equipment and storage medium | |
| US9849306B2 (en) | Systems and methods for generating beam-specific planning target volume design outputs | |
| CN111001097B (en) | Radiotherapy dose evaluation system, device and storage medium | |
| WO2017054316A1 (en) | Method of online monitoring and verifying dose in radiotherapy | |
| Olaciregui-Ruiz et al. | Virtual patient 3D dose reconstruction using in air EPID measurements and a back-projection algorithm for IMRT and VMAT treatments | |
| Bresciani et al. | Comparison of two different EPID-based solutions performing pretreatment quality assurance: 2D portal dosimetry versus 3D forward projection method | |
| Olaciregui‐Ruiz et al. | Site‐specific alert criteria to detect patient‐related errors with 3D EPID transit dosimetry | |
| CN110302475B (en) | Cloud Monte Carlo dose verification analysis method, equipment and storage medium | |
| CN110237445B (en) | In vivo three-dimensional dose monitoring and verification method based on EPID | |
| CN104027128A (en) | Offline dose verification method based on improved CBCT (cone beam computed tomography) images | |
| Amit et al. | Automated planning of breast radiotherapy using cone beam CT imaging | |
| Zhu et al. | Fast 3D dosimetric verifications based on an electronic portal imaging device using a GPU calculation engine | |
| Ecclestone et al. | Experimental validation of the van Herk margin formula for lung radiation therapy | |
| Helmbrecht et al. | Systematic analysis on the achievable accuracy of PT-PET through automated evaluation techniques | |
| Grogg et al. | Feasibility of using distal endpoints for in-room PET range verification of proton therapy | |
| Zhong et al. | Generation of a novel phase‐space‐based cylindrical dose kernel for IMRT optimization | |
| Liu et al. | A method to reconstruct and apply 3D primary fluence for treatment delivery verification | |
| Aland | Quality assurance of complex radiotherapy treatments using high-resolution 2D dosimeters | |
| CN117727423A (en) | Particle beam range real-time verification method and system |
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 |