WO2023178527A1 - Generation method and generation apparatus for tumor radiotherapy region - Google Patents

Generation method and generation apparatus for tumor radiotherapy region Download PDF

Info

Publication number
WO2023178527A1
WO2023178527A1 PCT/CN2022/082331 CN2022082331W WO2023178527A1 WO 2023178527 A1 WO2023178527 A1 WO 2023178527A1 CN 2022082331 W CN2022082331 W CN 2022082331W WO 2023178527 A1 WO2023178527 A1 WO 2023178527A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
network model
control area
dense network
parameter
Prior art date
Application number
PCT/CN2022/082331
Other languages
French (fr)
Chinese (zh)
Inventor
梁晓坤
产银萍
谢耀钦
张楚龙
何文丰
Original Assignee
中国科学院深圳先进技术研究院
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 中国科学院深圳先进技术研究院 filed Critical 中国科学院深圳先进技术研究院
Priority to PCT/CN2022/082331 priority Critical patent/WO2023178527A1/en
Publication of WO2023178527A1 publication Critical patent/WO2023178527A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration

Definitions

  • the present application relates to the field of medical imaging, and in particular to a method and device for generating a tumor radiotherapy area.
  • CT computed tomography
  • This application discloses a method and device for generating a tumor radiotherapy area to solve the problem that the positioning of a tumor patient's radiotherapy area requires manual intervention for correction.
  • embodiments of the present application provide a method for generating a tumor radiotherapy area, including: acquiring a first CT image actually obtained by a target object during treatment and a second CT image initially obtained before treatment; in response to External input instructions mark the position of the first control area on the second CT image to generate a third CT image; by inputting the first CT image and the third CT image into a dense network model, obtain the first translation parameters and a first rotation parameter; according to the first translation parameter and the first rotation parameter, the first CT image is transformed to generate a fourth CT image, wherein the fourth CT image is a Image of the target control area for radiotherapy treatment of the target subject.
  • embodiments of the present application provide a device for generating a tumor radiotherapy area, including: a first acquisition module for acquiring the first CT image actually obtained by the target object during the treatment and the initial CT image obtained before the treatment. a second CT image; a marking module configured to mark the position of the first control area on the second CT image in response to an external input instruction to generate a third CT image; a second acquisition module configured to generate a third CT image by A CT image and the third CT image are input into a dense network model to obtain a first translation parameter and a first rotation parameter; a generation module is used to generate the first translation parameter and the first rotation parameter according to the first translation parameter and the first rotation parameter.
  • a CT image is transformed to generate a fourth CT image, wherein the fourth CT image is an image with a target control area for radiotherapy of the target object.
  • inventions of the present application provide an electronic device.
  • the electronic device includes a processor and a memory.
  • the memory stores programs or instructions that can be run on the processor.
  • the programs or instructions are processed by the processor.
  • the processor is executed, the steps of the method described in the first aspect are implemented.
  • embodiments of the present application provide a readable storage medium.
  • Programs or instructions are stored on the readable storage medium.
  • the steps of the method described in the first aspect are implemented. .
  • Embodiments of the present application provide a method for generating a tumor radiotherapy area by acquiring the first CT image actually obtained by the target object during the treatment and the second CT image initially obtained before the treatment, and then in response to an external input instruction, The position of the first control area is marked on the second CT image to generate a third CT image, and then the first CT image and the third CT image are input into the dense network model to obtain the first translation parameter and the first rotation parameter, and then according to the A translation parameter and a first rotation parameter are used to transform the first CT image to generate a fourth CT image with a target control area for radiotherapy of the target object, thereby achieving positioning of the tumor radiotherapy area.
  • registration is achieved through the artifact-free first control area manually selected by the doctor, which can avoid the impact of uneven image information in global registration, thereby reducing potential uncertainty errors.
  • Figure 1 is a schematic flow chart of a method for generating a tumor radiotherapy area disclosed in an embodiment of the present application
  • Figure 2 is a schematic structural diagram of a device for generating a tumor radiotherapy area disclosed in an embodiment of the present application
  • FIG. 3 is a schematic structural diagram of an electronic device disclosed in an embodiment of the present application.
  • first, second, etc. in the description and claims of this application are used to distinguish similar objects and are not used to describe a specific order or sequence. It is to be understood that the figures so used are interchangeable under appropriate circumstances so that the embodiments of the present application can be practiced in orders other than those illustrated or described herein, and that "first,” “second,” etc. are distinguished Objects are usually of one type, and the number of objects is not limited. For example, the first object can be one or multiple.
  • “and/or” in the description and claims indicates at least one of the connected objects, and the character “/" generally indicates that the related objects are in an "or” relationship.
  • Figure 1 is a schematic flowchart of a method for generating a tumor radiotherapy area disclosed in an embodiment of the present application.
  • the method can be executed by an electronic device.
  • the method can be executed by software or hardware installed on the electronic device, as shown in Figure 1
  • the method includes the following steps.
  • the doctor manually selects and marks the position of the first control area on the second CT image to generate a third CT image, where the position of the first control area can be an artifact-free area that is easily identifiable by the doctor.
  • the first translation parameter and the first rotation parameter of the treatment bed are obtained.
  • the fourth CT image and the second CT image can be consistent in a spatial and anatomical sense, That is, the fourth CT image and the second CT image are most similar.
  • Embodiments of the present application provide a method for generating a tumor radiotherapy area by acquiring the first CT image actually obtained by the target object during the treatment and the second CT image initially obtained before the treatment, and then in response to an external input instruction, The position of the first control area is marked on the second CT image to generate a third CT image, and then the first CT image and the third CT image are input into the dense network model to obtain the first translation parameter and the first rotation parameter, and then according to the A translation parameter and a first rotation parameter are used to transform the first CT image to generate a fourth CT image with a target control area for radiotherapy of the target object, thereby achieving positioning of the tumor radiotherapy area.
  • registration is achieved through the artifact-free first control area manually selected by the doctor, which can avoid the impact of uneven image information in global registration, thereby reducing potential uncertainty errors.
  • the skeletal anatomy is not a simple rigid object. Therefore, for robustness, the positions of multiple first control areas can usually be marked on the second CT image.
  • the first translation parameter and the first rotation parameter are obtained, It may include: inputting the first CT image and the third CT image into a dense network model to obtain a plurality of first translation parameters and a plurality of first rotation parameters; and calculating respectively according to the plurality of first translation parameters and the plurality of first rotation parameters. The average value of the first translation parameter and the average value of the first rotation parameter.
  • multiple first translation parameters can be obtained by inputting the first CT image and the third CT image into the dense network model. and a plurality of first rotation parameters, by calculating the average of the first translation parameter and the first rotation parameter respectively, to determine the conversion parameter, thereby improving the robustness of the model.
  • the first CT image and the third CT image before inputting the first CT image and the third CT image into the dense network model, it may also include: acquiring multiple sets of data sets of the first object, wherein the data set includes the data sets actually obtained during the treatment process.
  • the position of the second control region is a randomly generated position in the sixth CT image; according to multiple sets of A training sample set is generated from a data set of an object; the training sample set is input into the dense network model to be trained for iterative training processing to obtain a trained dense network model.
  • inputting the training sample set into the dense network model to be trained for iterative training processing may include: marking the position of the second control area on the sixth CT image to generate a seventh CT image;
  • the CT image and the seventh CT image are input into the dense network model to obtain the second translation parameter and the second rotation parameter; according to the second translation parameter and the second rotation parameter, the fifth CT image is transformed to generate the eighth CT image; according to The position of the second control area, respectively extract the first target control area in the sixth CT image and the second target control area in the eighth CT image; determine the first target control area and the second target control area through normalized cross-correlation Local similarity between regions; determine the global similarity between the sixth CT image and the eighth CT image through normalized cross-correlation; determine the sixth CT image and the eighth CT image based on the local similarity and global similarity. loss function between.
  • the fifth CT image is transformed to generate the eighth CT image according to the second translation parameter and the second rotation parameter, according to the second control marked on the sixth CT image
  • the position of the area respectively extract the first target control area in the sixth CT image and the second target control area in the eighth CT image, and determine the relationship between the first target control area and the second target control area through normalized cross-correlation local similarity, determine the global similarity between the sixth CT image and the eighth CT image through normalized cross-correlation, and then determine the relationship between the sixth CT image and the eighth CT image based on the local similarity and global similarity. loss function.
  • the local similarity between the first target control area and the second target control area is used as the main part of the loss function.
  • the loss function between the sixth CT image and the eighth CT image is minimum, the image similarity between the sixth CT image and the eighth CT image is maximum.
  • LG represents the global similarity between the sixth CT image and the eighth CT image
  • L CV represents the first target control area and the second target control
  • represents the constant parameter of the relative weight of global similarity and local similarity
  • represents the learning parameter of the dense network model
  • v represents the position of the second control region in the sixth CT image
  • I d represents the third
  • I p represents the sixth CT image
  • I p,CV represents the first target control area in the sixth CT image
  • I′ d,CV represents the second target control area in the eighth CT image.
  • the local similarity between the first target control area and the second target control area, and the global similarity between the sixth CT image and the eighth CT image are calculated using negative normalized cross-correlation.
  • the specific calculation is The formula is as follows: Among them, I 1 and I 2 respectively represent two images or two target control areas, p is the index of the voxel, and It represents the average gray level of the voxels in I 1 and I 2 .
  • the loss function between the sixth CT image and the eighth CT image may also include: performing iterative training processing through the Adam (Adaptive Moment Estimation) algorithm, and saving the dense corresponding to the minimum loss function. Learning parameters of the network model.
  • the Adam algorithm can be used to continuously iteratively train the dense network model to minimize the loss function, and save the weight parameters corresponding to the minimum loss function and the learning parameters of the dense network model, so that when the dense network model is used later, By directly substituting the data, the corresponding results can be obtained without iterating again, which can speed up the process.
  • the learning rate in the Adam algorithm can be initially set to 10 -2 and then gradually reduced to 10 -6 , using a total of 10 5 iterations to train the dense network model.
  • a large amount of initial data can be generated by randomly translating and rotating the sixth CT image initially obtained before treatment, and by transforming the sixth CT image, and the dense network model to be trained can be processed through a large amount of initial data. Training can make the dense network model obtained by training more robust and avoid over-fitting.
  • the dense network model includes a convolution layer and three dense layers. Each dense layer is followed by a transition layer. The last dense layer is followed by a pooling layer and a linear layer.
  • the dense layer consists of a convolutional layer, a batch normalization layer and a sequence of rectified linear units, and the transition layer consists of a convolutional layer and a pooling layer.
  • the convolutional layer of the dense network model is used to extract different features of the input image.
  • the first convolutional layer extracts some low-level features such as edges, lines, corners, etc., and the multi-layer convolutional layer obtains deeper features and outputs Becomes a feature map; the pooling layer is used to compress the input feature map.
  • the batch normalization layer is used to make the distribution of input data relatively stable and accelerate the model learning speed;
  • the corrected linear unit sequence is used to provide the nonlinear expression modeling capability of the network;
  • the transition layer is used to reduce the feature mapping between adjacent dense blocks The dimensionality of the layer;
  • the linear layer is used to implement a linear combination or linear transformation of the previous layer, that is, convert their input features into output features.
  • the CT images actually obtained during the treatment process were collected.
  • four control regions are usually selected at the spinous process of the second cervical vertebra, the spinous process of the sixth cervical vertebra, the mandible and the skull.
  • the final transformation parameters are obtained by applying the 4 translation parameters and 4 rotation parameters output by the dense network model respectively.
  • a calculated average was determined using eight known and easily identifiable anatomical points marked by an experienced radiation oncologist as a reference.
  • the translation parameters and rotation parameters of the treatment bed using bone landmark alignment can be determined by the following formula:
  • R represents the rotation matrix of the treatment bed in the yaw, pitch and roll directions
  • T represents the translation matrix of the treatment bed in the front and rear, left and right and up and down directions
  • the value of R corresponding to when the formula obtains the minimum value is determined as the rotation parameter
  • the value of T corresponding to when the formula obtains the minimum value is determined as the translation parameter.
  • this application gives more weight to the artifact-free control area in the loss function (that is, taking local similarity as the main part of the loss function), compared with the entire cone beam CT image when using the traditional registration method
  • This application eliminates this misalignment by moving it rearward relative to its intended position.
  • the error of anatomical mapping in the control area is reduced, it still cannot be matched well due to metal image artifacts.
  • the anatomical structure can be matched very well. .
  • Table 1 is a summary of the translation and rotation errors of three different methods: the traditional registration method, the deep learning (DL) model without control volumes (CVs), and the DL model with CVs.
  • the execution subject may be a device for generating a tumor radiation treatment area.
  • a method for generating a tumor radiotherapy area performed by a device for generating a radiation treatment area for tumors is used as an example to illustrate the device for generating a radiation treatment area for tumors provided by embodiments of the application.
  • Figure 2 is a schematic structural diagram of a device for generating a tumor radiation treatment area disclosed in an embodiment of the present application.
  • the device 200 for generating a tumor radiation treatment area includes: a first acquisition module 210 , a marking module 220 , a second acquisition module 230 and a generation module 240 .
  • the first acquisition module 210 is used to acquire the first CT image actually obtained by the target object during the treatment and the second CT image initially obtained before the treatment;
  • the marking module 220 is used to respond to external input instructions. , mark the position of the first control area on the second CT image, and generate a third CT image;
  • the second acquisition module 230 is used to input the first CT image and the third CT image into a dense network model , obtain the first translation parameter and the first rotation parameter;
  • the generation module 240 is configured to transform the first CT image according to the first translation parameter and the first rotation parameter to generate a fourth CT image, wherein,
  • the fourth CT image is an image with a target control area for radiation therapy of the target subject.
  • the second acquisition module 230 when there are multiple positions of the first control areas marked on the second CT image, the second acquisition module 230 combines the first CT image and the The third CT image is input into the dense network model, and a first translation parameter and a first rotation parameter are obtained, including: inputting the first CT image and the third CT image into the dense network model, and obtaining a plurality of first translation parameters. and a plurality of first rotation parameters; according to a plurality of the first translation parameters and a plurality of the first rotation parameters, the average value of the first translation parameters and the average value of the first rotation parameters are respectively calculated.
  • the second acquisition module 230 is further configured to: acquire multiple sets of data sets of the first object before inputting the first CT image and the third CT image into the dense network model.
  • the data set includes the fifth CT image actually obtained during the treatment, the sixth CT image initially obtained before the treatment, and the position of the second control area in the sixth CT image, and the second control area
  • the position of the region is a randomly generated position in the sixth CT image;
  • a training sample set is generated according to the data sets of the multiple groups of first objects; and the training sample set is input into the dense network model to be trained for iteration Training process to obtain the dense network model that has been trained.
  • the second acquisition module 230 inputs the training sample set into the dense network model to be trained for iterative training processing, including: marking the second control on the sixth CT image The position of the area, generate a seventh CT image; by inputting the fifth CT image and the seventh CT image into the dense network model, obtain the second translation parameter and the second rotation parameter; according to the second translation parameter and the second rotation parameter, transform the fifth CT image to generate an eighth CT image; according to the position of the second control area, extract the first target control area and the first target control area in the sixth CT image respectively.
  • the second target control area in the eighth CT image determines the local similarity between the first target control area and the second target control area through normalized cross-correlation; determine the local similarity between the first target control area and the second target control area through normalized cross-correlation
  • the global similarity between the sixth CT image and the eighth CT image determine the loss between the sixth CT image and the eighth CT image according to the local similarity and the global similarity. function.
  • the second acquisition module 230 is further configured to: after determining the loss function between the sixth CT image and the eighth CT image, perform iterative training processing through the Adam algorithm. , save the learning parameters of the dense network model corresponding to the minimum loss function.
  • the device for generating a tumor radiation treatment area provided by the embodiments of the present application can implement various processes implemented by the method embodiments for generating a tumor radiation treatment area. To avoid duplication, they will not be described again here.
  • this embodiment of the present application also provides an electronic device 300, including a processor 301 and a memory 302.
  • the memory 302 stores programs or instructions that can be run on the processor 301.
  • the program or instruction is executed by the processor 301, each step of the above embodiment of the method for generating a tumor radiation treatment area is implemented, and the same technical effect can be achieved. To avoid duplication, the details will not be described here.
  • Embodiments of the present application also provide a readable storage medium, with a program or instructions stored on the readable storage medium.
  • the program or instructions are executed by a processor, each process of the above embodiment of the method for generating a tumor radiation treatment area is implemented. And can achieve the same technical effect. To avoid repetition, they will not be described again here.
  • the processor is the processor in the electronic device described in the above embodiment.
  • the readable storage medium includes computer readable storage media, such as computer read-only memory ROM, random access memory RAM, magnetic disk or optical disk, etc.
  • the methods of the above embodiments can be implemented by means of software plus the necessary general hardware platform. Of course, it can also be implemented by hardware, but in many cases the former is better. implementation.
  • the technical solution of the present application can be embodied in the form of a computer software product that is essentially or contributes to the existing technology.
  • the computer software product is stored in a storage medium (such as ROM/RAM, disk , optical disk), including several instructions to cause a terminal (which can be a mobile phone, computer, server, or network device, etc.) to execute the methods described in various embodiments of this application.

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

The present application discloses a generation method and generation apparatus for a tumor radiotherapy region. The generation method for a tumor radiotherapy region comprises: acquiring a first CT image actually obtained during a treatment process of a target object and a second CT image initially obtained before the treatment; in response to an external input instruction, marking the position of a first control region on the second CT image to generate a third CT image; obtaining a first translation parameter and a first rotation parameter by means of inputting the first CT image and the third CT image into a dense network model; transforming the first CT image according to the first translation parameter and the first rotation parameter to generate a fourth CT image, wherein the fourth CT image is an image provided with a target control region for the radiotherapy of the target object.

Description

肿瘤放射治疗区域的生成方法及生成装置Method and device for generating tumor radiotherapy area 技术领域Technical field
本申请涉及医学影像领域,尤其涉及一种肿瘤放射治疗区域的生成方法及生成装置。The present application relates to the field of medical imaging, and in particular to a method and device for generating a tumor radiotherapy area.
背景技术Background technique
对于解剖结构复杂的肿瘤(例如,头颈癌)的放射治疗来说,精准定位至关重要,而图像引导是实现高精度定位的关键步骤。目前,计算机断层扫描成像(computed tomography,CT)是图像引导中的常用工具,一般是通过求取初始得到的CT图像和在治疗过程中实际得到的CT图像间相似性最大的最优变换参数进行配准,以实现精准定位。For radiotherapy of tumors with complex anatomy (e.g., head and neck cancer), precise positioning is crucial, and image guidance is a key step to achieve high-precision positioning. At present, computed tomography (CT) is a commonly used tool in image guidance, which is generally performed by obtaining the optimal transformation parameters that maximize the similarity between the initially obtained CT image and the actual CT image obtained during the treatment process. Registration for precise positioning.
在相关技术中,虽然已经实现了将自动配准算法结合到图像引导放射治疗(image-guided radio therapy,IGRT)系统中,但是,由于图像中伪影的存在和在治疗过程中解剖结构的变异,可能会导致初始得到的CT图像和实际得到的CT图像之间的全局图像相似性度量出现偏差等因素的影响,导致肿瘤患者放射治疗区域的定位仍然需要人工干预进行校正。In the related art, although automatic registration algorithms have been implemented into image-guided radiotherapy (IGRT) systems, due to the presence of artifacts in images and variations in anatomical structures during treatment , may be affected by factors such as deviations in the global image similarity measure between the initially obtained CT image and the actually obtained CT image, resulting in the positioning of the radiotherapy area for tumor patients still requiring manual intervention for correction.
发明内容Contents of the invention
本申请公开一种肿瘤放射治疗区域的生成方法及生成装置,以解决肿瘤患者放射治疗区域的定位需要人工干预进行校正的问题。This application discloses a method and device for generating a tumor radiotherapy area to solve the problem that the positioning of a tumor patient's radiotherapy area requires manual intervention for correction.
为了解决上述问题,本申请采用下述技术方案:In order to solve the above problems, this application adopts the following technical solutions:
第一方面,本申请实施例提供了一种肿瘤放射治疗区域的生成方法,包 括:获取目标对象在治疗过程中实际得到的第一CT图像和在治疗前初始得到的第二CT图像;响应于外部输入指令,在所述第二CT图像上标记第一控制区域的位置,生成第三CT图像;通过将所述第一CT图像和所述第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数;根据所述第一平移参数和所述第一旋转参数,对所述第一CT图像进行变换生成第四CT图像,其中,所述第四CT图像为带有用于所述目标对象放射治疗的目标控制区域的图像。In a first aspect, embodiments of the present application provide a method for generating a tumor radiotherapy area, including: acquiring a first CT image actually obtained by a target object during treatment and a second CT image initially obtained before treatment; in response to External input instructions mark the position of the first control area on the second CT image to generate a third CT image; by inputting the first CT image and the third CT image into a dense network model, obtain the first translation parameters and a first rotation parameter; according to the first translation parameter and the first rotation parameter, the first CT image is transformed to generate a fourth CT image, wherein the fourth CT image is a Image of the target control area for radiotherapy treatment of the target subject.
第二方面,本申请实施例提供了一种肿瘤放射治疗区域的生成装置,包括:第一获取模块,用于获取目标对象在治疗过程中实际得到的第一CT图像和在治疗前初始得到的第二CT图像;标记模块,用于响应于外部输入指令,在所述第二CT图像上标记第一控制区域的位置,生成第三CT图像;第二获取模块,用于通过将所述第一CT图像和所述第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数;生成模块,用于根据所述第一平移参数和所述第一旋转参数,对所述第一CT图像进行变换生成第四CT图像,其中,所述第四CT图像为带有用于所述目标对象放射治疗的目标控制区域的图像。In a second aspect, embodiments of the present application provide a device for generating a tumor radiotherapy area, including: a first acquisition module for acquiring the first CT image actually obtained by the target object during the treatment and the initial CT image obtained before the treatment. a second CT image; a marking module configured to mark the position of the first control area on the second CT image in response to an external input instruction to generate a third CT image; a second acquisition module configured to generate a third CT image by A CT image and the third CT image are input into a dense network model to obtain a first translation parameter and a first rotation parameter; a generation module is used to generate the first translation parameter and the first rotation parameter according to the first translation parameter and the first rotation parameter. A CT image is transformed to generate a fourth CT image, wherein the fourth CT image is an image with a target control area for radiotherapy of the target object.
第三方面,本申请实施例提供了一种电子设备,该电子设备包括处理器和存储器,所述存储器存储可在所述处理器上运行的程序或指令,所述程序或指令被所述处理器执行时实现如第一方面所述的方法的步骤。In a third aspect, embodiments of the present application provide an electronic device. The electronic device includes a processor and a memory. The memory stores programs or instructions that can be run on the processor. The programs or instructions are processed by the processor. When the processor is executed, the steps of the method described in the first aspect are implemented.
第四方面,本申请实施例提供了一种可读存储介质,所述可读存储介质上存储程序或指令,所述程序或指令被处理器执行时实现如第一方面所述的方法的步骤。In a fourth aspect, embodiments of the present application provide a readable storage medium. Programs or instructions are stored on the readable storage medium. When the programs or instructions are executed by a processor, the steps of the method described in the first aspect are implemented. .
本申请实施例提供一种肿瘤放射治疗区域的生成方法,通过获取目标对象在治疗过程中实际得到的第一CT图像和在治疗前初始得到的第二CT图像,然后响应于外部输入指令,在第二CT图像上标记第一控制区域的位置,生成第三CT图像,再通过将第一CT图像和第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数,然后根据第一平移参数和第一旋转参数,对第一CT图像进行变换生成带有用于目标对象放射治疗的目标控制区域的第四CT图像,从而实现对肿瘤放射治疗区域的定位。另外,通过医生手动选择的无伪影的第一控制区域实现配准,能够避免全局配准中由于图像信息不均匀带来的影响,进而能够减少潜在的不确定误差。Embodiments of the present application provide a method for generating a tumor radiotherapy area by acquiring the first CT image actually obtained by the target object during the treatment and the second CT image initially obtained before the treatment, and then in response to an external input instruction, The position of the first control area is marked on the second CT image to generate a third CT image, and then the first CT image and the third CT image are input into the dense network model to obtain the first translation parameter and the first rotation parameter, and then according to the A translation parameter and a first rotation parameter are used to transform the first CT image to generate a fourth CT image with a target control area for radiotherapy of the target object, thereby achieving positioning of the tumor radiotherapy area. In addition, registration is achieved through the artifact-free first control area manually selected by the doctor, which can avoid the impact of uneven image information in global registration, thereby reducing potential uncertainty errors.
附图说明Description of the drawings
图1为本申请实施例公开的一种肿瘤放射治疗区域的生成方法的流程示意图;Figure 1 is a schematic flow chart of a method for generating a tumor radiotherapy area disclosed in an embodiment of the present application;
图2为本申请实施例公开的一种肿瘤放射治疗区域的生成装置的结构示意图;Figure 2 is a schematic structural diagram of a device for generating a tumor radiotherapy area disclosed in an embodiment of the present application;
图3为本申请实施例公开的电子设备的一种结构示意图。FIG. 3 is a schematic structural diagram of an electronic device disclosed in an embodiment of the present application.
具体实施方式Detailed ways
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员获得的所有其他实施例,都属于本申请保护的范围。The technical solutions in the embodiments of the present application will be clearly described below with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, but not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by those of ordinary skill in the art fall within the scope of protection of this application.
本申请的说明书和权利要求书中的术语“第一”、“第二”等是用于区别类似的对象,而不用于描述特定的顺序或先后次序。应该理解这样使用的 数据在适当情况下可以互换,以便本申请的实施例能够以除了在这里图示或描述的那些以外的顺序实施,且“第一”、“第二”等所区分的对象通常为一类,并不限定对象的个数,例如第一对象可以是一个,也可以是多个。此外,说明书以及权利要求中“和/或”表示所连接对象的至少其中之一,字符“/”,一般表示前后关联对象是一种“或”的关系。The terms "first", "second", etc. in the description and claims of this application are used to distinguish similar objects and are not used to describe a specific order or sequence. It is to be understood that the figures so used are interchangeable under appropriate circumstances so that the embodiments of the present application can be practiced in orders other than those illustrated or described herein, and that "first," "second," etc. are distinguished Objects are usually of one type, and the number of objects is not limited. For example, the first object can be one or multiple. In addition, "and/or" in the description and claims indicates at least one of the connected objects, and the character "/" generally indicates that the related objects are in an "or" relationship.
下面结合附图,通过具体的实施例及其应用场景对本申请实施例提供的肿瘤放射治疗区域的生成方法进行详细地说明。The method for generating a tumor radiation treatment area provided by embodiments of the present application will be described in detail below with reference to the accompanying drawings through specific embodiments and application scenarios.
图1为本申请实施例公开的一种肿瘤放射治疗区域的生成方法的流程示意图,该方法可以由电子设备执行,换言之,该方法可以由安装在电子设备的软件或硬件来执行,如图1所示,该方法包括如下步骤。Figure 1 is a schematic flowchart of a method for generating a tumor radiotherapy area disclosed in an embodiment of the present application. The method can be executed by an electronic device. In other words, the method can be executed by software or hardware installed on the electronic device, as shown in Figure 1 As shown, the method includes the following steps.
S120、获取目标对象在治疗过程中实际得到的第一CT图像和在治疗前初始得到的第二CT图像。S120. Obtain the first CT image actually obtained by the target object during the treatment and the second CT image initially obtained before the treatment.
S140、响应于外部输入指令,在第二CT图像上标记第一控制区域的位置,生成第三CT图像。S140. In response to the external input instruction, mark the position of the first control area on the second CT image to generate a third CT image.
也就是说,通过医生手动选择,在第二CT图像上标记第一控制区域的位置,生成第三CT图像,其中,第一控制区域的位置可以为医生容易识别的、无伪影的区域。That is to say, the doctor manually selects and marks the position of the first control area on the second CT image to generate a third CT image, where the position of the first control area can be an artifact-free area that is easily identifiable by the doctor.
S160、通过将第一CT图像和第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数。S160. Obtain the first translation parameter and the first rotation parameter by inputting the first CT image and the third CT image into the dense network model.
即通过将第一CT图像和标记有第一控制区域位置的第二CT图像输入密集网络模型,获取治疗床的第一平移参数和第一旋转参数。That is, by inputting the first CT image and the second CT image marked with the position of the first control area into the dense network model, the first translation parameter and the first rotation parameter of the treatment bed are obtained.
S180、根据第一平移参数和第一旋转参数,对第一CT图像进行变换生成第四CT图像,其中,第四CT图像为带有用于目标对象放射治疗的目标控制区域的图像。S180. According to the first translation parameter and the first rotation parameter, transform the first CT image to generate a fourth CT image, where the fourth CT image is an image with a target control area for radiotherapy of the target object.
通过根据第一平移参数和第一旋转参数,对第一CT图像进行变换生成第四CT图像,从而实现配准,能够使得第四CT图像与第二CT图像达到空间及解剖意义上的一致,也即第四CT图像与第二CT图像之间最相似。By transforming the first CT image to generate a fourth CT image according to the first translation parameter and the first rotation parameter, registration is achieved, and the fourth CT image and the second CT image can be consistent in a spatial and anatomical sense, That is, the fourth CT image and the second CT image are most similar.
本申请实施例提供一种肿瘤放射治疗区域的生成方法,通过获取目标对象在治疗过程中实际得到的第一CT图像和在治疗前初始得到的第二CT图像,然后响应于外部输入指令,在第二CT图像上标记第一控制区域的位置,生成第三CT图像,再通过将第一CT图像和第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数,然后根据第一平移参数和第一旋转参数,对第一CT图像进行变换生成带有用于目标对象放射治疗的目标控制区域的第四CT图像,从而实现对肿瘤放射治疗区域的定位。另外,通过医生手动选择的无伪影的第一控制区域实现配准,能够避免全局配准中由于图像信息不均匀带来的影响,进而能够减少潜在的不确定误差。Embodiments of the present application provide a method for generating a tumor radiotherapy area by acquiring the first CT image actually obtained by the target object during the treatment and the second CT image initially obtained before the treatment, and then in response to an external input instruction, The position of the first control area is marked on the second CT image to generate a third CT image, and then the first CT image and the third CT image are input into the dense network model to obtain the first translation parameter and the first rotation parameter, and then according to the A translation parameter and a first rotation parameter are used to transform the first CT image to generate a fourth CT image with a target control area for radiotherapy of the target object, thereby achieving positioning of the tumor radiotherapy area. In addition, registration is achieved through the artifact-free first control area manually selected by the doctor, which can avoid the impact of uneven image information in global registration, thereby reducing potential uncertainty errors.
以头颈癌(head and neck cancer,HNC)为例,在HNC区域,骨骼解剖并不是一个简单的刚性物体,因此,为了稳健,通常可以在第二CT图像上标记多个第一控制区域的位置,而在第二CT图像上标记的第一控制区域的位置为多个的情况下,通过将第一CT图像和第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数,可以包括:将第一CT图像和第三CT图像输入密集网络模型,获取多个第一平移参数和多个第一旋转参数;根据多个第一平移参数和多个第一旋转参数,分别计算第一平移参数的平均值和第一旋转参数的平均值。也就是说,在第二CT图像上标记的第一控制区域 的位置为多个的情况下,通过将第一CT图像和第三CT图像输入密集网络模型,可以获取到多个第一平移参数和多个第一旋转参数,通过分别计算第一平移参数和第一旋转参数的平均值,确定转换参数,进而提升模型的稳健性。Taking head and neck cancer (HNC) as an example, in the HNC area, the skeletal anatomy is not a simple rigid object. Therefore, for robustness, the positions of multiple first control areas can usually be marked on the second CT image. , and when there are multiple positions of the first control areas marked on the second CT image, by inputting the first CT image and the third CT image into the dense network model, the first translation parameter and the first rotation parameter are obtained, It may include: inputting the first CT image and the third CT image into a dense network model to obtain a plurality of first translation parameters and a plurality of first rotation parameters; and calculating respectively according to the plurality of first translation parameters and the plurality of first rotation parameters. The average value of the first translation parameter and the average value of the first rotation parameter. That is to say, when there are multiple positions of the first control areas marked on the second CT image, multiple first translation parameters can be obtained by inputting the first CT image and the third CT image into the dense network model. and a plurality of first rotation parameters, by calculating the average of the first translation parameter and the first rotation parameter respectively, to determine the conversion parameter, thereby improving the robustness of the model.
在本申请实施例中,在将第一CT图像和第三CT图像输入密集网络模型之前,还可以包括:获取多组第一对象的数据集,其中,数据集包括在治疗过程中实际得到的第五CT图像、在治疗前初始得到的第六CT图像和第六CT图像中第二控制区域的位置,第二控制区域的位置为在第六CT图像中随机生成的位置;根据多组第一对象的数据集生成训练样本集;将训练样本集输入待训练的密集网络模型进行迭代训练处理,得到训练完成的密集网络模型。本申请通过采用在第六CT图像中随机生成的第二控制区域的位置,无需手动标注数据,能够在不依赖于标注数据集,极大地节省人力的情况下,完成对密集网络模型的训练。In the embodiment of the present application, before inputting the first CT image and the third CT image into the dense network model, it may also include: acquiring multiple sets of data sets of the first object, wherein the data set includes the data sets actually obtained during the treatment process. The fifth CT image, the sixth CT image initially obtained before treatment, and the position of the second control region in the sixth CT image. The position of the second control region is a randomly generated position in the sixth CT image; according to multiple sets of A training sample set is generated from a data set of an object; the training sample set is input into the dense network model to be trained for iterative training processing to obtain a trained dense network model. By using the position of the second control area randomly generated in the sixth CT image, this application eliminates the need for manual annotation of data and can complete the training of dense network models without relying on annotation data sets, greatly saving manpower.
在一种实现方式中,将训练样本集输入待训练的密集网络模型进行迭代训练处理,可以包括:在第六CT图像上标记第二控制区域的位置,生成第七CT图像;通过将第五CT图像和第七CT图像输入所述密集网络模型,获取第二平移参数和第二旋转参数;根据第二平移参数和第二旋转参数,对第五CT图像进行变换生成第八CT图像;根据第二控制区域的位置,分别提取第六CT图像中的第一目标控制区域和第八CT图像中的第二目标控制区域;通过归一化互相关确定第一目标控制区域和第二目标控制区域之间的局部相似度;通过归一化互相关确定第六CT图像和第八CT图像之间的全局相似度;根据局部相似度和全局相似度,确定第六CT图像和第八CT图像之间的损失函数。In one implementation, inputting the training sample set into the dense network model to be trained for iterative training processing may include: marking the position of the second control area on the sixth CT image to generate a seventh CT image; The CT image and the seventh CT image are input into the dense network model to obtain the second translation parameter and the second rotation parameter; according to the second translation parameter and the second rotation parameter, the fifth CT image is transformed to generate the eighth CT image; according to The position of the second control area, respectively extract the first target control area in the sixth CT image and the second target control area in the eighth CT image; determine the first target control area and the second target control area through normalized cross-correlation Local similarity between regions; determine the global similarity between the sixth CT image and the eighth CT image through normalized cross-correlation; determine the sixth CT image and the eighth CT image based on the local similarity and global similarity. loss function between.
也就是说,在训练密集网络模型的过程中,在根据第二平移参数和第二旋转参数,对第五CT图像进行变换生成第八CT图像之后,根据第六CT图像上标记的第二控制区域的位置,分别提取第六CT图像中的第一目标控制区域和第八CT图像中的第二目标控制区域,通过归一化互相关确定第一目标控制区域和第二目标控制区域之间的局部相似度,通过归一化互相关确定第六CT图像和第八CT图像之间的全局相似度,然后根据局部相似度和全局相似度,确定第六CT图像和第八CT图像之间的损失函数。并且为了避免局部伪影或解剖差异对全局的影响,提高医生感兴趣区域的配准精度,将第一目标控制区域和第二目标控制区域之间的局部相似度作为损失函数的主要部分。在第六CT图像和第八CT图像之间的损失函数最小的情况下,第六CT图像和第八CT图像之间的图像相似性最大。That is to say, in the process of training the dense network model, after the fifth CT image is transformed to generate the eighth CT image according to the second translation parameter and the second rotation parameter, according to the second control marked on the sixth CT image The position of the area, respectively extract the first target control area in the sixth CT image and the second target control area in the eighth CT image, and determine the relationship between the first target control area and the second target control area through normalized cross-correlation local similarity, determine the global similarity between the sixth CT image and the eighth CT image through normalized cross-correlation, and then determine the relationship between the sixth CT image and the eighth CT image based on the local similarity and global similarity. loss function. And in order to avoid the global impact of local artifacts or anatomical differences and improve the registration accuracy of the doctor's area of interest, the local similarity between the first target control area and the second target control area is used as the main part of the loss function. When the loss function between the sixth CT image and the eighth CT image is minimum, the image similarity between the sixth CT image and the eighth CT image is maximum.
其中,第六CT图像和第八CT图像之间的损失函数可以通过下述公式进行表示:L(θ;I d,I p,v)=L G(θ;I p,I d)+λL CV(θ;I p,CV;I′ d,CV),其中,L G表示第六CT图像和第八CT图像之间的全局相似度,L CV表示第一目标控制区域和第二目标控制区域之间的局部相似度,λ表示全局相似度和局部相似度相对权重的常数参数,θ表示密集网络模型的学习参数,v表示第六CT图像中第二控制区域的位置,I d表示第五CT图像,I p表示第六CT图像,I p,CV表示第六CT图像中的第一目标控制区域,I′ d,CV表示第八CT图像中的第二目标控制区域。 Among them, the loss function between the sixth CT image and the eighth CT image can be expressed by the following formula: L (θ; I d , I p , v) = L G (θ; I p , I d ) + λL CV (θ; I p, CV ; I′ d, CV ), where LG represents the global similarity between the sixth CT image and the eighth CT image, and L CV represents the first target control area and the second target control The local similarity between regions, λ represents the constant parameter of the relative weight of global similarity and local similarity, θ represents the learning parameter of the dense network model, v represents the position of the second control region in the sixth CT image, I d represents the third In the fifth CT image, I p represents the sixth CT image, I p,CV represents the first target control area in the sixth CT image, and I′ d,CV represents the second target control area in the eighth CT image.
而第一目标控制区域和第二目标控制区域之间的局部相似度,以及第六CT图像和第八CT图像之间的全局相似度,均使用负的归一化互相关进行计算,具体计算公式如下:
Figure PCTCN2022082331-appb-000001
其中,I 1和 I 2分别表示两幅图像或者两个目标控制区域,p是其中的体素的索引,
Figure PCTCN2022082331-appb-000002
Figure PCTCN2022082331-appb-000003
则表示I 1和I 2中体素的灰度平均值。
The local similarity between the first target control area and the second target control area, and the global similarity between the sixth CT image and the eighth CT image are calculated using negative normalized cross-correlation. The specific calculation is The formula is as follows:
Figure PCTCN2022082331-appb-000001
Among them, I 1 and I 2 respectively represent two images or two target control areas, p is the index of the voxel,
Figure PCTCN2022082331-appb-000002
and
Figure PCTCN2022082331-appb-000003
It represents the average gray level of the voxels in I 1 and I 2 .
在本申请实施例中,在确定第六CT图像和第八CT图像之间的损失函数之后,还可以包括:通过Adam(Adaptive Moment Estimation)算法进行迭代训练处理,保存损失函数最小时对应的密集网络模型的学习参数。也就是说,可以通过Adam算法不断对密集网络模型进行迭代训练,使得损失函数最小,并且保存损失函数最小时对应的权重参数和密集网络模型的学习参数,以便于在后续使用密集网络模型时,通过直接代入数据,即可得到对应结果,而无需经过再次迭代,能够加快速度。示例性的,为了减小损失函数,Adam算法中的学习率最初可以设定为10 -2,然后逐渐降低到10 -6,总共使用10 5次迭代来训练密集网络模型。此外,在进行迭代训练处理时,可以通过随机平移和旋转治疗前初始得到的第六CT图像,通过对第六CT图像进行变换,生成大量初始数据,通过大量初始数据对待训练的密集网络模型进行训练,能够使得训练得到的密集网络模型的鲁棒性更好,同时能够避免过渡拟合。 In the embodiment of the present application, after determining the loss function between the sixth CT image and the eighth CT image, it may also include: performing iterative training processing through the Adam (Adaptive Moment Estimation) algorithm, and saving the dense corresponding to the minimum loss function. Learning parameters of the network model. In other words, the Adam algorithm can be used to continuously iteratively train the dense network model to minimize the loss function, and save the weight parameters corresponding to the minimum loss function and the learning parameters of the dense network model, so that when the dense network model is used later, By directly substituting the data, the corresponding results can be obtained without iterating again, which can speed up the process. For example, in order to reduce the loss function, the learning rate in the Adam algorithm can be initially set to 10 -2 and then gradually reduced to 10 -6 , using a total of 10 5 iterations to train the dense network model. In addition, when performing iterative training processing, a large amount of initial data can be generated by randomly translating and rotating the sixth CT image initially obtained before treatment, and by transforming the sixth CT image, and the dense network model to be trained can be processed through a large amount of initial data. Training can make the dense network model obtained by training more robust and avoid over-fitting.
在本申请实施例中,密集网络模型中包括卷积层和三个密集层,每一个密集层后面都跟着一个过渡层,最后一个密集层后面紧跟着一个池化层和线性层,而密集层由卷积层、批量归一化层和校正线性单元序列组成,过渡层由卷积层和池化层组成。其中,密集网络模型的卷积层用于提取输入的图像的不同特征,第一层卷积层提取一些低级的特征如边缘、线条和角等,多层卷积层得到更深层次的特征,输出成为特征图;池化层用于对输入的特征图进行压缩,一方面使特征图变小,简化网络计算复杂度,一方面进行特征压缩,提取主要特征,还可以用于防止过拟合;批量归一化层用于使输入数据的分布相对稳定,加速模型学习速度;校正线性单元序列用于提供网络的非 线性表达建模能力;过渡层用于降低相邻密集块之间的特征映射的维数;线性层用于实现对前一层的线性组合或线性变换,即将它们的输入特征转换为输出特征。In the embodiment of this application, the dense network model includes a convolution layer and three dense layers. Each dense layer is followed by a transition layer. The last dense layer is followed by a pooling layer and a linear layer. The dense layer The layer consists of a convolutional layer, a batch normalization layer and a sequence of rectified linear units, and the transition layer consists of a convolutional layer and a pooling layer. Among them, the convolutional layer of the dense network model is used to extract different features of the input image. The first convolutional layer extracts some low-level features such as edges, lines, corners, etc., and the multi-layer convolutional layer obtains deeper features and outputs Becomes a feature map; the pooling layer is used to compress the input feature map. On the one hand, it makes the feature map smaller and simplifies the network calculation complexity. On the other hand, it compresses the features and extracts the main features. It can also be used to prevent overfitting; The batch normalization layer is used to make the distribution of input data relatively stable and accelerate the model learning speed; the corrected linear unit sequence is used to provide the nonlinear expression modeling capability of the network; the transition layer is used to reduce the feature mapping between adjacent dense blocks The dimensionality of the layer; the linear layer is used to implement a linear combination or linear transformation of the previous layer, that is, convert their input features into output features.
以HNC为例,在实验验证时,对于每个患者,均采集了治疗过程中实际得到的CT图像、在治疗前初始得到的CT图像和在初始得到的CT图像中手动选择的控制区域。为了稳健,通常选择第二颈椎棘突、第六颈椎棘突、下颌骨和颅骨处的四个控制区域,最终的转换参数是通过对密集网络模型输出的4个平移参数和4个旋转参数分别计算平均值确定,并以一名经验丰富的放射肿瘤学家标记的八个已知且易于识别的解剖点作为参考。Taking HNC as an example, during the experimental verification, for each patient, the CT images actually obtained during the treatment process, the CT images initially obtained before treatment, and the manually selected control area in the initially obtained CT images were collected. For robustness, four control regions are usually selected at the spinous process of the second cervical vertebra, the spinous process of the sixth cervical vertebra, the mandible and the skull. The final transformation parameters are obtained by applying the 4 translation parameters and 4 rotation parameters output by the dense network model respectively. A calculated average was determined using eight known and easily identifiable anatomical points marked by an experienced radiation oncologist as a reference.
在本申请中,使用骨标志对齐的治疗床的平移参数和旋转参数可以通过下述公式确定:
Figure PCTCN2022082331-appb-000004
其中,R表示治疗床在偏航、俯仰和侧滚方向上的旋转矩阵,T表示治疗床在前后、左右以及上下方向的平移矩阵,
Figure PCTCN2022082331-appb-000005
表示治疗过程中实际得到的CT图像中的控制区域,
Figure PCTCN2022082331-appb-000006
表示在治疗前初始得到的CT图像中的控制区域,将公式取得最小值时对应的R的值确定为旋转参数,将公式取得最小值时对应的T的值确定为平移参数。
In this application, the translation parameters and rotation parameters of the treatment bed using bone landmark alignment can be determined by the following formula:
Figure PCTCN2022082331-appb-000004
Among them, R represents the rotation matrix of the treatment bed in the yaw, pitch and roll directions, T represents the translation matrix of the treatment bed in the front and rear, left and right and up and down directions,
Figure PCTCN2022082331-appb-000005
Represents the control area in the CT image actually obtained during the treatment process,
Figure PCTCN2022082331-appb-000006
Represents the control area in the CT image initially obtained before treatment. The value of R corresponding to when the formula obtains the minimum value is determined as the rotation parameter, and the value of T corresponding to when the formula obtains the minimum value is determined as the translation parameter.
由于本申请在损失函数中对无伪影的控制区域赋予了更多的权重(即将局部相似度作为损失函数的主要部分),因此,相较于使用传统配准方法时整个锥形束CT图像相对于其预期位置向后移动来说,本申请消除了这种未对齐。另外,通过传统方法配准后,控制区域中的解剖映射的误差虽然减小,但是由于金属图像伪影,仍然不能很好地匹配,而通过本申请提供的方法,解剖结构能非常好地匹配。Since this application gives more weight to the artifact-free control area in the loss function (that is, taking local similarity as the main part of the loss function), compared with the entire cone beam CT image when using the traditional registration method This application eliminates this misalignment by moving it rearward relative to its intended position. In addition, after registration through traditional methods, although the error of anatomical mapping in the control area is reduced, it still cannot be matched well due to metal image artifacts. However, through the method provided by this application, the anatomical structure can be matched very well. .
表1是对传统配准方法、无控制区域(control volumes,CVs)的深度学习(DeepLearning,DL)模型、有CVs的DL模型三种不同方法的平移和旋转误差的总结。Table 1 is a summary of the translation and rotation errors of three different methods: the traditional registration method, the deep learning (DL) model without control volumes (CVs), and the DL model with CVs.
表1 不同方法的平移和旋转误差总结Table 1 Summary of translation and rotation errors of different methods
Figure PCTCN2022082331-appb-000007
Figure PCTCN2022082331-appb-000007
通过表1,可以看到,在平移和旋转中,本方法(有CVs的DL模型)测量和参考之间的系统/随机定位误差分别小于0.47/1.13mm和0.17/0.29°。并且,在预先定义的临床可接受的公差(2.0mm/1.0°)内,所提出的方法的测试分数在平移和旋转中的比例分别为87.88%和100.00%。HNC病例在临床可 接受的范围内的比例由传统方法的66.67%提高到本方法的90.91%。并且,通过基于深度学习的配准结果的定量比较,在有CVs时,我们发现90.91%的病例在临床上是可接受的,这显著提高了63.64%的没有CVs的可接受比例。Through Table 1, it can be seen that in translation and rotation, the systematic/random positioning error between the measurement and reference of this method (DL model with CVs) is less than 0.47/1.13mm and 0.17/0.29° respectively. And, within the predefined clinically acceptable tolerance (2.0mm/1.0°), the test scores of the proposed method were 87.88% and 100.00% in translation and rotation, respectively. The proportion of HNC cases within the clinically acceptable range increased from 66.67% with the traditional method to 90.91% with this method. And, through quantitative comparison of deep learning-based registration results, we found that 90.91% of cases were clinically acceptable with CVs, which significantly improved the acceptable ratio of 63.64% without CVs.
以上数据均可表明本申请提出的方法能够高效的实现高精度的患者定位。The above data can show that the method proposed in this application can effectively achieve high-precision patient positioning.
本申请实施例提供的肿瘤放射治疗区域的生成方法,执行主体可以为肿瘤放射治疗区域的生成装置。本申请实施例中以肿瘤放射治疗区域的生成装置执行肿瘤放射治疗区域的生成方法为例,说明本申请实施例提供的肿瘤放射治疗区域的生成装置。For the method for generating a tumor radiation treatment area provided by an embodiment of the present application, the execution subject may be a device for generating a tumor radiation treatment area. In the embodiments of the present application, a method for generating a tumor radiotherapy area performed by a device for generating a radiation treatment area for tumors is used as an example to illustrate the device for generating a radiation treatment area for tumors provided by embodiments of the application.
图2为本申请实施例公开的一种肿瘤放射治疗区域的生成装置的结构示意图。如图2所示,肿瘤放射治疗区域的生成装置200包括:第一获取模块210、标记模块220、第二获取模块230和生成模块240。Figure 2 is a schematic structural diagram of a device for generating a tumor radiation treatment area disclosed in an embodiment of the present application. As shown in FIG. 2 , the device 200 for generating a tumor radiation treatment area includes: a first acquisition module 210 , a marking module 220 , a second acquisition module 230 and a generation module 240 .
在本申请中,第一获取模块210,用于获取目标对象在治疗过程中实际得到的第一CT图像和在治疗前初始得到的第二CT图像;标记模块220,用于响应于外部输入指令,在所述第二CT图像上标记第一控制区域的位置,生成第三CT图像;第二获取模块230,用于通过将所述第一CT图像和所述第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数;生成模块240,用于根据所述第一平移参数和所述第一旋转参数,对所述第一CT图像进行变换生成第四CT图像,其中,所述第四CT图像为带有用于所述目标对象放射治疗的目标控制区域的图像。In this application, the first acquisition module 210 is used to acquire the first CT image actually obtained by the target object during the treatment and the second CT image initially obtained before the treatment; the marking module 220 is used to respond to external input instructions. , mark the position of the first control area on the second CT image, and generate a third CT image; the second acquisition module 230 is used to input the first CT image and the third CT image into a dense network model , obtain the first translation parameter and the first rotation parameter; the generation module 240 is configured to transform the first CT image according to the first translation parameter and the first rotation parameter to generate a fourth CT image, wherein, The fourth CT image is an image with a target control area for radiation therapy of the target subject.
在一种实现方式中,在所述第二CT图像上标记的所述第一控制区域的位置为多个的情况下,所述第二获取模块230通过将所述第一CT图像和所述第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数,包 括:将所述第一CT图像和所述第三CT图像输入所述密集网络模型,获取多个第一平移参数和多个第一旋转参数;根据多个所述第一平移参数和多个所述第一旋转参数,分别计算所述第一平移参数的平均值和所述第一旋转参数的平均值。In one implementation, when there are multiple positions of the first control areas marked on the second CT image, the second acquisition module 230 combines the first CT image and the The third CT image is input into the dense network model, and a first translation parameter and a first rotation parameter are obtained, including: inputting the first CT image and the third CT image into the dense network model, and obtaining a plurality of first translation parameters. and a plurality of first rotation parameters; according to a plurality of the first translation parameters and a plurality of the first rotation parameters, the average value of the first translation parameters and the average value of the first rotation parameters are respectively calculated.
在一种实现方式中,所述第二获取模块230还用于:在所述将所述第一CT图像和所述第三CT图像输入密集网络模型之前,获取多组第一对象的数据集,其中,所述数据集包括在治疗过程中实际得到的第五CT图像、在治疗前初始得到的第六CT图像和所述第六CT图像中第二控制区域的位置,所述第二控制区域的位置为在所述第六CT图像中随机生成的位置;根据所述多组第一对象的数据集生成训练样本集;将所述训练样本集输入待训练的所述密集网络模型进行迭代训练处理,得到训练完成的所述密集网络模型。In one implementation, the second acquisition module 230 is further configured to: acquire multiple sets of data sets of the first object before inputting the first CT image and the third CT image into the dense network model. , wherein the data set includes the fifth CT image actually obtained during the treatment, the sixth CT image initially obtained before the treatment, and the position of the second control area in the sixth CT image, and the second control area The position of the region is a randomly generated position in the sixth CT image; a training sample set is generated according to the data sets of the multiple groups of first objects; and the training sample set is input into the dense network model to be trained for iteration Training process to obtain the dense network model that has been trained.
在一种实现方式中,所述第二获取模块230将所述训练样本集输入待训练的所述密集网络模型进行迭代训练处理,包括:在所述第六CT图像上标记所述第二控制区域的位置,生成第七CT图像;通过将所述第五CT图像和所述第七CT图像输入所述密集网络模型,获取第二平移参数和第二旋转参数;根据所述第二平移参数和所述第二旋转参数,对所述第五CT图像进行变换生成第八CT图像;根据所述第二控制区域的位置,分别提取所述第六CT图像中的第一目标控制区域和所述第八CT图像中的第二目标控制区域;通过归一化互相关确定所述第一目标控制区域和所述第二目标控制区域之间的局部相似度;通过归一化互相关确定所述第六CT图像和所述第八CT图像之间的全局相似度;根据所述局部相似度和所述全局相似度,确定所述第六CT图像和所述第八CT图像之间的损失函数。In one implementation, the second acquisition module 230 inputs the training sample set into the dense network model to be trained for iterative training processing, including: marking the second control on the sixth CT image The position of the area, generate a seventh CT image; by inputting the fifth CT image and the seventh CT image into the dense network model, obtain the second translation parameter and the second rotation parameter; according to the second translation parameter and the second rotation parameter, transform the fifth CT image to generate an eighth CT image; according to the position of the second control area, extract the first target control area and the first target control area in the sixth CT image respectively. the second target control area in the eighth CT image; determine the local similarity between the first target control area and the second target control area through normalized cross-correlation; determine the local similarity between the first target control area and the second target control area through normalized cross-correlation The global similarity between the sixth CT image and the eighth CT image; determine the loss between the sixth CT image and the eighth CT image according to the local similarity and the global similarity. function.
在一种实现方式中,所述第二获取模块230,还用于:在所述确定所述第六CT图像和所述第八CT图像之间的损失函数之后,通过Adam算法进行迭代训练处理,保存所述损失函数最小时对应的所述密集网络模型的学习参数。In one implementation, the second acquisition module 230 is further configured to: after determining the loss function between the sixth CT image and the eighth CT image, perform iterative training processing through the Adam algorithm. , save the learning parameters of the dense network model corresponding to the minimum loss function.
本申请实施例提供的肿瘤放射治疗区域的生成装置能够实现肿瘤放射治疗区域的生成方法实施例实现的各个过程,为避免重复,这里不再赘述。The device for generating a tumor radiation treatment area provided by the embodiments of the present application can implement various processes implemented by the method embodiments for generating a tumor radiation treatment area. To avoid duplication, they will not be described again here.
可选地,如图3所示,本申请实施例还提供一种电子设备300,包括处理器301和存储器302,存储器302上存储有可在所述处理器301上运行的程序或指令,该程序或指令被处理器301执行时实现上述肿瘤放射治疗区域的生成方法实施例的各个步骤,且能达到相同的技术效果,为避免重复,这里不再赘述。Optionally, as shown in Figure 3, this embodiment of the present application also provides an electronic device 300, including a processor 301 and a memory 302. The memory 302 stores programs or instructions that can be run on the processor 301. When the program or instruction is executed by the processor 301, each step of the above embodiment of the method for generating a tumor radiation treatment area is implemented, and the same technical effect can be achieved. To avoid duplication, the details will not be described here.
本申请实施例还提供一种可读存储介质,所述可读存储介质上存储有程序或指令,该程序或指令被处理器执行时实现上述肿瘤放射治疗区域的生成方法实施例的各个过程,且能达到相同的技术效果,为避免重复,这里不再赘述。Embodiments of the present application also provide a readable storage medium, with a program or instructions stored on the readable storage medium. When the program or instructions are executed by a processor, each process of the above embodiment of the method for generating a tumor radiation treatment area is implemented. And can achieve the same technical effect. To avoid repetition, they will not be described again here.
其中,所述处理器为上述实施例中所述的电子设备中的处理器。所述可读存储介质,包括计算机可读存储介质,如计算机只读存储器ROM、随机存取存储器RAM、磁碟或者光盘等。Wherein, the processor is the processor in the electronic device described in the above embodiment. The readable storage medium includes computer readable storage media, such as computer read-only memory ROM, random access memory RAM, magnetic disk or optical disk, etc.
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过 程、方法、物品或者装置中还存在另外的相同要素。此外,需要指出的是,本申请实施方式中的方法和装置的范围不限按示出或讨论的顺序来执行功能,还可包括根据所涉及的功能按基本同时的方式或按相反的顺序来执行功能,例如,可以按不同于所描述的次序来执行所描述的方法,并且还可以添加、省去、或组合各种步骤。另外,参照某些示例所描述的特征可在其他示例中被组合。It should be noted that, in this document, the terms "comprising", "comprising" or any other variations thereof are intended to cover a non-exclusive inclusion, such that a process, method, article or device that includes a series of elements not only includes those elements, It also includes other elements not expressly listed or inherent in the process, method, article or apparatus. Without further limitation, an element defined by the statement "comprises a..." does not exclude the presence of other identical elements in the process, method, article or device that includes the element. In addition, it should be pointed out that the scope of the methods and devices in the embodiments of the present application is not limited to performing functions in the order shown or discussed, but may also include performing functions in a substantially simultaneous manner or in reverse order according to the functions involved. Functions may be performed, for example, the methods described may be performed in an order different from that described, and various steps may be added, omitted, or combined. Additionally, features described with reference to certain examples may be combined in other examples.
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以计算机软件产品的形式体现出来,该计算机软件产品存储在一个存储介质(如ROM/RAM、磁碟、光盘)中,包括若干指令用以使得一台终端(可以是手机,计算机,服务器,或者网络设备等)执行本申请各个实施例所述的方法。Through the above description of the embodiments, those skilled in the art can clearly understand that the methods of the above embodiments can be implemented by means of software plus the necessary general hardware platform. Of course, it can also be implemented by hardware, but in many cases the former is better. implementation. Based on this understanding, the technical solution of the present application can be embodied in the form of a computer software product that is essentially or contributes to the existing technology. The computer software product is stored in a storage medium (such as ROM/RAM, disk , optical disk), including several instructions to cause a terminal (which can be a mobile phone, computer, server, or network device, etc.) to execute the methods described in various embodiments of this application.
上面结合附图对本申请的实施例进行了描述,但是本申请并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本申请的启示下,在不脱离本申请宗旨和权利要求所保护的范围情况下,还可做出很多形式,均属于本申请的保护之内。The embodiments of the present application have been described above in conjunction with the accompanying drawings. However, the present application is not limited to the above-mentioned specific implementations. The above-mentioned specific implementations are only illustrative and not restrictive. Those of ordinary skill in the art will Inspired by this application, many forms can be made without departing from the purpose of this application and the scope protected by the claims, all of which fall within the protection of this application.

Claims (10)

  1. 一种肿瘤放射治疗区域的生成方法,包括:A method for generating a tumor radiotherapy area, including:
    获取目标对象在治疗过程中实际得到的第一CT图像和在治疗前初始得到的第二CT图像;Obtaining the first CT image actually obtained by the target object during the treatment and the second CT image initially obtained before the treatment;
    响应于外部输入指令,在所述第二CT图像上标记第一控制区域的位置,生成第三CT图像;In response to an external input instruction, mark the position of the first control area on the second CT image to generate a third CT image;
    通过将所述第一CT图像和所述第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数;By inputting the first CT image and the third CT image into a dense network model, a first translation parameter and a first rotation parameter are obtained;
    根据所述第一平移参数和所述第一旋转参数,对所述第一CT图像进行变换生成第四CT图像,其中,所述第四CT图像为带有用于所述目标对象放射治疗的目标控制区域的图像。According to the first translation parameter and the first rotation parameter, the first CT image is transformed to generate a fourth CT image, wherein the fourth CT image is with a target for radiotherapy of the target object. Image of the control area.
  2. 根据权利要求1所述的生成方法,其中,在所述第二CT图像上标记的所述第一控制区域的位置为多个的情况下,所述通过将所述第一CT图像和所述第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数,包括:The generation method according to claim 1, wherein when there are multiple positions of the first control areas marked on the second CT image, the first CT image and the The third CT image is input into the dense network model to obtain the first translation parameter and the first rotation parameter, including:
    将所述第一CT图像和所述第三CT图像输入所述密集网络模型,获取多个第一平移参数和多个第一旋转参数;Input the first CT image and the third CT image into the dense network model to obtain a plurality of first translation parameters and a plurality of first rotation parameters;
    根据多个所述第一平移参数和多个所述第一旋转参数,分别计算所述第一平移参数的平均值和所述第一旋转参数的平均值。According to a plurality of the first translation parameters and a plurality of the first rotation parameters, an average value of the first translation parameters and an average value of the first rotation parameters are respectively calculated.
  3. 根据权利要求1所述的生成方法,其中,在所述将所述第一CT图像和所述第三CT图像输入密集网络模型之前,还包括:The generation method according to claim 1, wherein before inputting the first CT image and the third CT image into the dense network model, it further includes:
    获取多组第一对象的数据集,其中,所述数据集包括在治疗过程中实际 得到的第五CT图像、在治疗前初始得到的第六CT图像和所述第六CT图像中第二控制区域的位置,所述第二控制区域的位置为在所述第六CT图像中随机生成的位置;Acquire multiple sets of data sets of the first object, wherein the data set includes a fifth CT image actually obtained during the treatment, a sixth CT image initially obtained before treatment, and a second control in the sixth CT image. The position of the region, the position of the second control region is a randomly generated position in the sixth CT image;
    根据所述多组第一对象的数据集生成训练样本集;Generate a training sample set according to the data sets of the plurality of first objects;
    将所述训练样本集输入待训练的所述密集网络模型进行迭代训练处理,得到训练完成的所述密集网络模型。The training sample set is input into the dense network model to be trained for iterative training processing to obtain the trained dense network model.
  4. 根据权利要求3所述的生成方法,其中,所述将所述训练样本集输入待训练的所述密集网络模型进行迭代训练处理,包括:The generation method according to claim 3, wherein said inputting the training sample set into the dense network model to be trained for iterative training processing includes:
    在所述第六CT图像上标记所述第二控制区域的位置,生成第七CT图像;Mark the position of the second control area on the sixth CT image to generate a seventh CT image;
    通过将所述第五CT图像和所述第七CT图像输入所述密集网络模型,获取第二平移参数和第二旋转参数;By inputting the fifth CT image and the seventh CT image into the dense network model, a second translation parameter and a second rotation parameter are obtained;
    根据所述第二平移参数和所述第二旋转参数,对所述第五CT图像进行变换生成第八CT图像;Transform the fifth CT image to generate an eighth CT image according to the second translation parameter and the second rotation parameter;
    根据所述第二控制区域的位置,分别提取所述第六CT图像中的第一目标控制区域和所述第八CT图像中的第二目标控制区域;According to the position of the second control area, respectively extract the first target control area in the sixth CT image and the second target control area in the eighth CT image;
    通过归一化互相关确定所述第一目标控制区域和所述第二目标控制区域之间的局部相似度;Determining the local similarity between the first target control area and the second target control area by normalized cross-correlation;
    通过归一化互相关确定所述第六CT图像和所述第八CT图像之间的全局相似度;Determine the global similarity between the sixth CT image and the eighth CT image by normalized cross-correlation;
    根据所述局部相似度和所述全局相似度,确定所述第六CT图像和所述第八CT图像之间的损失函数。According to the local similarity and the global similarity, a loss function between the sixth CT image and the eighth CT image is determined.
  5. 根据权利要求4所述的生成方法,其中,在所述确定所述第六CT图像和所述第八CT图像之间的损失函数之后,还包括:The generation method according to claim 4, wherein, after determining the loss function between the sixth CT image and the eighth CT image, further comprising:
    通过Adam算法进行迭代训练处理,保存所述损失函数最小时对应的所述密集网络模型的学习参数。Iterative training is performed through the Adam algorithm, and the learning parameters of the dense network model corresponding to the minimum loss function are saved.
  6. 一种肿瘤放射治疗区域的生成装置,包括:A device for generating a tumor radiotherapy area, including:
    第一获取模块,用于获取目标对象在治疗过程中实际得到的第一CT图像和在治疗前初始得到的第二CT图像;The first acquisition module is used to acquire the first CT image actually obtained by the target object during the treatment and the second CT image initially obtained before treatment;
    标记模块,用于响应于外部输入指令,在所述第二CT图像上标记第一控制区域的位置,生成第三CT图像;A marking module, configured to mark the position of the first control area on the second CT image in response to an external input instruction, and generate a third CT image;
    第二获取模块,用于通过将所述第一CT图像和所述第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数;a second acquisition module, configured to acquire the first translation parameter and the first rotation parameter by inputting the first CT image and the third CT image into a dense network model;
    生成模块,用于根据所述第一平移参数和所述第一旋转参数,对所述第一CT图像进行变换生成第四CT图像,其中,所述第四CT图像为带有用于所述目标对象放射治疗的目标控制区域的图像。Generating module, configured to transform the first CT image to generate a fourth CT image according to the first translation parameter and the first rotation parameter, wherein the fourth CT image is a Image of the target control area of a subject receiving radiation therapy.
  7. 根据权利要求6所述的生成装置,其中,在所述第二CT图像上标记的所述第一控制区域的位置为多个的情况下,所述第二获取模块通过将所述第一CT图像和所述第三CT图像输入密集网络模型,获取第一平移参数和第一旋转参数,包括:The generating device according to claim 6, wherein when there are multiple positions of the first control area marked on the second CT image, the second acquisition module The image and the third CT image are input into the dense network model to obtain the first translation parameter and the first rotation parameter, including:
    将所述第一CT图像和所述第三CT图像输入所述密集网络模型,获取多个第一平移参数和多个第一旋转参数;Input the first CT image and the third CT image into the dense network model to obtain a plurality of first translation parameters and a plurality of first rotation parameters;
    根据多个所述第一平移参数和多个所述第一旋转参数,分别计算所述第一平移参数的平均值和所述第一旋转参数的平均值。According to a plurality of the first translation parameters and a plurality of the first rotation parameters, an average value of the first translation parameters and an average value of the first rotation parameters are respectively calculated.
  8. 根据权利要求6所述的生成装置,其中,所述第二获取模块还用于:The generating device according to claim 6, wherein the second acquisition module is also used for:
    在所述将所述第一CT图像和所述第三CT图像输入密集网络模型之前,获取多组第一对象的数据集,其中,所述数据集包括在治疗过程中实际得到 的第五CT图像、在治疗前初始得到的第六CT图像和所述第六CT图像中第二控制区域的位置,所述第二控制区域的位置为在所述第六CT图像中随机生成的位置;Before inputting the first CT image and the third CT image into the dense network model, a plurality of data sets of the first object are obtained, wherein the data set includes a fifth CT actually obtained during the treatment process. image, the sixth CT image initially obtained before treatment and the position of the second control region in the sixth CT image, where the position of the second control region is a randomly generated position in the sixth CT image;
    根据所述多组第一对象的数据集生成训练样本集;Generate a training sample set according to the data sets of the plurality of first objects;
    将所述训练样本集输入待训练的所述密集网络模型进行迭代训练处理,得到训练完成的所述密集网络模型。The training sample set is input into the dense network model to be trained for iterative training processing to obtain the trained dense network model.
  9. 根据权利要求8所述的生成装置,其中,所述第二获取模块将所述训练样本集输入待训练的所述密集网络模型进行迭代训练处理,包括:The generating device according to claim 8, wherein the second acquisition module inputs the training sample set into the dense network model to be trained for iterative training processing, including:
    在所述第六CT图像上标记所述第二控制区域的位置,生成第七CT图像;Mark the position of the second control area on the sixth CT image to generate a seventh CT image;
    通过将所述第五CT图像和所述第七CT图像输入所述密集网络模型,获取第二平移参数和第二旋转参数;By inputting the fifth CT image and the seventh CT image into the dense network model, a second translation parameter and a second rotation parameter are obtained;
    根据所述第二平移参数和所述第二旋转参数,对所述第五CT图像进行变换生成第八CT图像;Transform the fifth CT image to generate an eighth CT image according to the second translation parameter and the second rotation parameter;
    根据所述第二控制区域的位置,分别提取所述第六CT图像中的第一目标控制区域和所述第八CT图像中的第二目标控制区域;According to the position of the second control area, respectively extract the first target control area in the sixth CT image and the second target control area in the eighth CT image;
    通过归一化互相关确定所述第一目标控制区域和所述第二目标控制区域之间的局部相似度;Determining the local similarity between the first target control area and the second target control area by normalized cross-correlation;
    通过归一化互相关确定所述第六CT图像和所述第八CT图像之间的全局相似度;Determine the global similarity between the sixth CT image and the eighth CT image by normalized cross-correlation;
    根据所述局部相似度和所述全局相似度,确定所述第六CT图像和所述第八CT图像之间的损失函数。According to the local similarity and the global similarity, a loss function between the sixth CT image and the eighth CT image is determined.
  10. 根据权利要求9所述的生成装置,其中,所述第二获取模块,还用于:The generating device according to claim 9, wherein the second acquisition module is also used for:
    在所述确定所述第六CT图像和所述第八CT图像之间的损失函数之后,通过Adam算法进行迭代训练处理,保存所述损失函数最小时对应的所述密集网络模型的学习参数。After the loss function between the sixth CT image and the eighth CT image is determined, an iterative training process is performed through the Adam algorithm, and the learning parameters of the dense network model corresponding to the minimum loss function are saved.
PCT/CN2022/082331 2022-03-22 2022-03-22 Generation method and generation apparatus for tumor radiotherapy region WO2023178527A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/CN2022/082331 WO2023178527A1 (en) 2022-03-22 2022-03-22 Generation method and generation apparatus for tumor radiotherapy region

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2022/082331 WO2023178527A1 (en) 2022-03-22 2022-03-22 Generation method and generation apparatus for tumor radiotherapy region

Publications (1)

Publication Number Publication Date
WO2023178527A1 true WO2023178527A1 (en) 2023-09-28

Family

ID=88099614

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/082331 WO2023178527A1 (en) 2022-03-22 2022-03-22 Generation method and generation apparatus for tumor radiotherapy region

Country Status (1)

Country Link
WO (1) WO2023178527A1 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105469366A (en) * 2015-11-23 2016-04-06 山东科技大学 Analytic method for eliminating metal artifact of CT image
US20200034948A1 (en) * 2018-07-27 2020-01-30 Washington University Ml-based methods for pseudo-ct and hr mr image estimation
US20210233645A1 (en) * 2020-01-23 2021-07-29 GE Precision Healthcare LLC Methods and systems for characterizing anatomical features in medical images
US20210268313A1 (en) * 2020-02-27 2021-09-02 Elekta ltd Image Registration in Treatment Planning

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105469366A (en) * 2015-11-23 2016-04-06 山东科技大学 Analytic method for eliminating metal artifact of CT image
US20200034948A1 (en) * 2018-07-27 2020-01-30 Washington University Ml-based methods for pseudo-ct and hr mr image estimation
US20210233645A1 (en) * 2020-01-23 2021-07-29 GE Precision Healthcare LLC Methods and systems for characterizing anatomical features in medical images
US20210268313A1 (en) * 2020-02-27 2021-09-02 Elekta ltd Image Registration in Treatment Planning

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
LIANG XIAOKUN; BASSENNE MAXIME; HRISTOV DIMITRE H.; ISLAM MD TAUHIDUL; ZHAO WEI; JIA MENGYU; ZHANG ZHICHENG; GENSHEIMER MICHAEL; B: "Human-level comparable control volume mapping with a deep unsupervised-learning model for image-guided radiation therapy", COMPUTERS IN BIOLOGY AND MEDICINE, NEW YORK, NY, US, vol. 141, 17 December 2021 (2021-12-17), US , XP086945813, ISSN: 0010-4825, DOI: 10.1016/j.compbiomed.2021.105139 *

Similar Documents

Publication Publication Date Title
CN106920234B (en) Combined automatic radiotherapy planning method
CN106485695B (en) Medical image Graph Cut dividing method based on statistical shape model
CN106934821B (en) Conical beam CT and CT image registration method based on ICP algorithm and B spline
US8068650B2 (en) Lesion quantification and tracking using multiple modalities
CN109509193B (en) Liver CT atlas segmentation method and system based on high-precision registration
CN110378881B (en) Tumor positioning system based on deep learning
CN107240128B (en) X-ray and color photo registration method based on contour features
CN110246580B (en) Cranial image analysis method and system based on neural network and random forest
CN110751187B (en) Training method of abnormal area image generation network and related product
CN113421226B (en) CT-DR multi-modal esophageal image registration method and system based on mutual information
WO2012069965A1 (en) Interactive deformation map corrections
CN108898578B (en) Medical image processing method and device and computer storage medium
CN104599268A (en) Local area accurate deformation registration algorithm combining point registration
CN114792326A (en) Surgical navigation point cloud segmentation and registration method based on structured light
Dai et al. Locating anatomical landmarks on 2D lateral cephalograms through adversarial encoder-decoder networks
Nakao et al. Image-to-graph convolutional network for deformable shape reconstruction from a single projection image
CN111261296A (en) Tumor clinical target area automatic delineation method and system based on conditional random vector field
CN111462018B (en) Image alignment method in image, computer equipment and storage medium
CN111192268A (en) Medical image segmentation model construction method and CBCT image bone segmentation method
WO2020114332A1 (en) Segmentation-network-based ct lung tumor segmentation method, apparatus and device, and medium
CN110378910B (en) Abdominal cavity multi-organ segmentation method and device based on map fusion
WO2023178527A1 (en) Generation method and generation apparatus for tumor radiotherapy region
CN114820740A (en) Image registration method, storage medium and computer device
Sun et al. Imaging study of pseudo-CT images of superposed ultrasound deformation fields acquired in radiotherapy based on step-by-step local registration
CN116168097A (en) Method, device, equipment and medium for constructing CBCT sketching model and sketching CBCT image

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22932596

Country of ref document: EP

Kind code of ref document: A1