CN116152437B - 施源器重建方法、装置、电子设备和计算机可读存储介质 - Google Patents

施源器重建方法、装置、电子设备和计算机可读存储介质 Download PDF

Info

Publication number
CN116152437B
CN116152437B CN202310103177.3A CN202310103177A CN116152437B CN 116152437 B CN116152437 B CN 116152437B CN 202310103177 A CN202310103177 A CN 202310103177A CN 116152437 B CN116152437 B CN 116152437B
Authority
CN
China
Prior art keywords
plane
dimensional
mask image
coordinate axis
applicator
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
Application number
CN202310103177.3A
Other languages
English (en)
Other versions
CN116152437A (zh
Inventor
邱杰
于浪
白璐
李森林
陈颀
王少彬
陈宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking Union Medical College Hospital Chinese Academy of Medical Sciences
Original Assignee
Peking Union Medical College Hospital Chinese Academy of Medical Sciences
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 Peking Union Medical College Hospital Chinese Academy of Medical Sciences filed Critical Peking Union Medical College Hospital Chinese Academy of Medical Sciences
Priority to CN202310103177.3A priority Critical patent/CN116152437B/zh
Publication of CN116152437A publication Critical patent/CN116152437A/zh
Application granted granted Critical
Publication of CN116152437B publication Critical patent/CN116152437B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • 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/1001—X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy using radiation sources introduced into or applied onto the body; brachytherapy
    • 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/1049—Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00—Manipulating 3D models or images for computer graphics
    • G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G06T7/11—Region-based segmentation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20084—Artificial neural networks [ANN]

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Biomedical Technology (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Software Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Graphics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Architecture (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本申请的实施例提供了一种施源器重建方法、装置、电子设备和计算机可读存储介质，该方法包括：基于施源器管体的轴位掩膜图像，生成施源器管体的三维掩膜图像，施源器管体至少包括直线部分和弯曲部分；在预设三维坐标系下，对三维掩膜图像进行切割、旋转，使得直线部分与第一坐标轴重合、弯曲部分位于第一平面或第二平面上；生成施源器管体的冠状位切面和矢状位切面。本申请能够更精确地获取得到施源器的驻留位置，从而能够确定放射源可停留的治疗位置、可能治疗位置相对于靶区以及危及器官的距离。

Description

施源器重建方法、装置、电子设备和计算机可读存储介质
技术领域
本申请的实施例涉及图像处理技术领域，尤其涉及一种施源器重建方法、装置、电子设备和计算机可读存储介质。
背景技术
对于局部晚期宫颈癌，放疗是其重要的治疗手段，体外照射加腔内后装近距离放射治疗是其经典治疗模式，随着近年来强调放疗技术的发展以及图像引导下的三维后装技术的应用和推广，宫颈癌放射治疗的疗效得到了更大的提高。然而在施源器重建时会不可避免地出现偏差，导致获取得到的施源器的驻留位置与施源器的实际驻留位置不一致。
发明内容
为了提高施源器重建后驻留位置的准确性，本申请的实施例提供了施源器重建方法、装置、电子设备和计算机可读存储介质。
在本申请的第一方面，提供了一种施源器重建方法，包括：
基于施源器管体的轴位掩膜图像，生成所述施源器管体的三维掩膜图像，所述施源器管体至少包括直线部分和弯曲部分；
在预设三维坐标系下，对所述三维掩膜图像进行切割、旋转，使得所述直线部分与第一坐标轴重合、所述弯曲部分位于第一平面或第二平面上；
生成所述施源器管体的冠状位切面和矢状位切面；
其中，第一平面由第一坐标轴和第二坐标轴形成，第二平面由第一坐标轴和第三坐标轴形成。
在一种可能的实现方式中，所述在预设三维坐标系下，对所述三维掩膜图像进行切割、旋转，包括：
若所述施源器管体为非中间管体，
在预设三维坐标系下，在p1点处沿第一平面方向切割所述三维掩膜图像；
旋转第一平面，使得所述直线部分的中心线与第一坐标轴重合，并更新预设三维坐标系；
在更新后的预设三维坐标系下，在p1点处沿第二平面方向切割所述三维掩膜图像；
旋转第二平面，使得所述直线部分的中心线与第一坐标轴重合，并再次更新预设三维坐标系；
将所述施源器管体投影至再次更新后的预设三维坐标系的第三平面上；
旋转第三平面，使得投影部分的中心线与第二坐标轴或第三坐标轴重合；
其中，p1点为所述直线部分下端任一点，第三平面由第二坐标轴和第三坐标轴形成。
在一种可能的实现方式中，所述在预设三维坐标系下，对所述三维掩膜图像进行切割、旋转，包括：
若所述施源器管体为中间管体，
在预设三维坐标系下，在p1点处沿第一平面方向切割所述三维掩膜图像；
旋转第一平面，使得所述直线部分的中心线与第一坐标轴重合，并更新预设三维坐标系；
在更新后的预设三维坐标系下，在p1点处沿第二平面方向切割所述三维掩膜图像；
旋转第二平面，使得所述直线部分的中心线与第一坐标轴重合，并再次更新预设三维坐标系；
在再次更新后的预设三维坐标系下，确定所述直线部分的中心点q1和所述弯曲部分的中心点q2；
旋转第三平面，使得中心点q1和中心点q2的横坐标位于第三坐标轴上，或者使得中心点q1和中心点q2的纵坐标位于第二坐标轴上；
其中，p1点为所述直线部分下端任一点，第三平面由第二坐标轴和第三坐标轴形成。
在一种可能的实现方式中，确定所述直线部分的中心点q1和所述弯曲部分的中心点q2，包括：
在p1点处沿第三平面方向切割所述三维掩膜图像，得到所述施源器管体的第一轴位图像，基于所述第一轴位图像确定中心点q1；
在p2点处沿第三平面方向切割所述三维掩膜图像，得到所述施源器管体的第二轴位图像，基于所述第二轴位图像确定中心点q2；
其中，p2点为所述弯曲部分的顶点。
在一种可能的实现方式中，基于施源器管体的轴位掩膜图像，生成所述施源器管体的三维掩膜图像，包括：
将所述施源器管体的所述轴位掩膜图像进行叠放，生成所述施源器管体的所述三维掩膜图像。
在一种可能的实现方式中，所述基于施源器管体的轴位掩膜图像，生成所述施源器管体的三维掩膜图像之前，还包括：
获取多个角度的医疗影像；
根据所述医疗影像，通过U-Net网络模型，得到所述施源器管体的所述轴位掩膜图像。
在一种可能的实现方式中，所述方法还包括：
基于骨架提取算法，提取所述施源器管体的初始中心线；
延长所述初始中心线的两端，形成所述施源器管体的中心线。
在本申请的第二方面，提供了一种施源器重建装置，包括：
第一生成模块，用于基于施源器管体的轴位掩膜图像，生成所述施源器管体的三维掩膜图像，所述施源器管体至少包括直线部分和弯曲部分；
管体旋转模块，用于在预设三维坐标系下，对所述三维掩膜图像进行切割、旋转，使得所述直线部分与第一坐标轴重合、所述弯曲部分位于第一平面或第二平面上；
第二生成模块，用于生成所述施源器管体的冠状位切面和矢状位切面；
其中，第一平面由第一坐标轴和第二坐标轴形成，第二平面由第一坐标轴和第三坐标轴形成。
在本申请的第三方面，提供了一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述处理器执行所述计算机程序时实现如第一方面中任一项所述的方法。
在本申请的第四方面，提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如第一方面中任一项所述的方法。
在本申请实施例提供的施源器重建方法、装置、电子设备和计算机可读存储介质中，基于施源器管体的轴位掩膜图像生成施源器管体的三维掩膜图像，在预设三维坐标系下，对三维掩膜图像进行切割、旋转，使得施源器管体的直线部分与第一坐标轴重合、施源器管体的弯曲部分位于第一平面或第二平面上，生成施源器管体的冠状位切面和矢状位切面，通过上述方法对施源器进行空间位置重建，能够更精确地获取得到施源器的驻留位置，从而能够确定放射源可停留的治疗位置、可能治疗位置相对于靶区以及危及器官的距离。
应当理解，发明内容部分中所描述的内容并非旨在限定本申请的实施例的关键或重要特征，亦非用于限制本申请的范围。本申请的其它特征将通过以下的描述变得容易理解。
附图说明
图1示出了根据本申请的实施例的施源器重建方法的流程图。
图2示出了根据本申请的实施例的施源器mpr切面的示意图；
图3示出了根据本申请的实施例的施源器的中心线的示意图；
图4示出了根据本申请的实施例的施源器重建装置的方框图；
图5示出了适于用来实现本申请实施例的电子设备的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述。
本申请实施例提供的施源器重建方法，可以由电子设备执行，该电子设备可以为服务器也可以为终端设备，其中，该服务器可以是独立的物理服务器，也可以是多个物理服务器构成的服务器集群或者分布式系统，还可以是提供云计算服务的云服务器。终端设备可以是智能手机、平板电脑、笔记本电脑、台式计算机等，但并不局限于此。需要说明的是，若该电子设备为服务器时，需要外接显示设备来显示相关界面。
图1示出了根据本申请的实施例的施源器重建方法的流程图。在一些实施例中，图1所示施源器重建方法可以由上述电子设备执行。如图1所示，该方法包括以下步骤：
步骤101，基于施源器管体的轴位掩膜图像，生成施源器管体的三维掩膜图像。
在本申请实施例中，施源器可以包括至少一根施源器管体，每根施源器管体均包括直线部分和弯曲部分。需要说明的是，下文中是以施源器包括三根施源器管体为例来进行解释说明，并不构成对本申请的限制。
在本申请实施例中，可以通过医疗影像数据(反映施源器在病灶区域的位置)，根据训练完成的分割模型分割出每根施源器管体的轴位掩膜图像。需要说明的是，医疗影像数据例如可以使用CT扫描图像，也可以使用核磁影响，或者是其他类型的医疗影像。分割模型可以采用神经网络模型，例如U-Net深度学习网络模型。
在一些实施例中，可以将施源器管体的轴位掩膜图像进行叠放来生成施源器管体的三维掩膜图像。
步骤102，在预设三维坐标系下，对三维掩膜图像进行切割、旋转，使得直线部分与第一坐标轴重合、弯曲部分位于第一平面或第二平面上，其中，第一平面由第一坐标轴和第二坐标轴形成，第二平面由第一坐标轴和第三坐标轴形成。
需要说明的是，为了使得表述更加清楚，下文中将“第一坐标轴”称作“z轴”、“第二坐标轴”称作“y轴”、“第三坐标轴”称作“x轴”、“第一平面”称作“yz面”、“第二平面”称作“xz面”、“第三平面”称作“xy面”。
在本申请实施例中，在对三维掩膜图像进行切割、旋转时，需要确定该三维掩膜图像所对应的施源器管体是否为中间管体，其确定方式将在后文中进行介绍。下面详细介绍对三维掩膜图像进行切割、旋转。
施源器管体为非中间管体
在预设三维坐标系下，在p1点处沿yz面方向切割三维掩膜图像；
旋转yz面，使得直线部分的中心线与z轴重合，并更新预设三维坐标系；
在更新后的预设三维坐标系下，在p1点处沿xz面方向切割三维掩膜图像；
旋转xz面，使得直线部分的中心线与z轴重合，并再次更新预设三维坐标系；
将施源器管体投影至再次更新后的预设三维坐标系的xy面上；
旋转xy面，使得投影部分的中心线与y轴或x轴重合；
施源器管体为中间管体
在预设三维坐标系下，在p1点处沿yz面方向切割所述三维掩膜图像；
旋转yz面，使得直线部分的中心线与z轴重合，并更新预设三维坐标系；
在更新后的预设三维坐标系下，在p1点处沿xz面方向切割所述三维掩膜图像；
旋转xz面，使得直线部分的中心线与z轴重合，并再次更新预设三维坐标系；
在再次更新后的预设三维坐标系下，确定直线部分的中心点q1和弯曲部分的中心点q2；
旋转xy面，使得中心点q1和中心点q2的横坐标位于x轴上，或者使得中心点q1和中心点q2的纵坐标位于y轴上；
在本申请实施例中，p1点为直线部分下端任一点，p2点为所述弯曲部分的顶点，旋转yz面的旋转角度可由直线部分的中心点和z轴之间的夹角来确定，旋转xz面的旋转角度可由直线部分的中心线和z轴之间的夹角来确定，旋转xy面的旋转角度可由投影部分的中心线与y轴或x轴之间的夹角来确定。
下面来介绍确定直线部分的中心点q1和弯曲部分的中心点q2。
在p1点处沿xy面方向切割三维掩膜图像，得到施源器管体的第一轴位图像，基于第一轴位图像确定中心点q1。
在p2点处沿xy面方向切割三维掩膜图像，得到施源器管体的第二轴位图像，基于第二轴位图像确定中心点q2。
需要说明的是，可以通过人工的方式来判断施源器管体是否为中间管体，还可以根据施源器管体的三维掩膜图像的直线部分中心线与z轴重合后，其投影至xy面上的面积来进行判断，例如施源器包括三根施源器管体时，中间管体的投影小于其他两根管体的投影。
还需要说明的是，上述仅是以施源器包括三根施源器管体为例来进行说明，若施源器包括一根管体或偶数跟管体，无法判断某根管体是否为中间管体，则可以在施源器管体的直线部分与z轴重合后，通过其投影至xy面上的投影来判断，例如投影面积小于预设面积，则判断其为“中间管体”，采用中间管体的重建方式来对其进行重建。
步骤103，生成施源器管体的冠状位切面和矢状位切面。
参见图2，图2中间部分为施源器的轴位掩膜图像，右上角部分为施源器的冠状位切面，右下角部分为施源器的矢状位切面。
在一些实施例中，在生成冠状位切面和矢状位切面之后，还可以基于骨架提取算法，提取施源器管体的初始中心线，由于骨架提取算法未提取到施源器两端末尾的中心线，所以需要延长初始中心线的两端，形成施源器管体的中心线。参见图3，图3示出了施源器的中心线的示意图。
根据本申请实施例，基于施源器管体的轴位掩膜图像生成施源器管体的三维掩膜图像，在预设三维坐标系下，对三维掩膜图像进行切割、旋转，使得施源器管体的直线部分与第一坐标轴重合、施源器管体的弯曲部分位于第一平面或第二平面上，生成施源器管体的冠状位切面和矢状位切面，通过上述方法对施源器进行空间位置重建，能够更精确地获取得到施源器的驻留位置，从而能够确定放射源可停留的治疗位置、可能治疗位置相对于靶区以及危及器官的距离。
在一些实施例中，为了使得得到的施源器管体的轴位掩膜图像更加准确，可以获取多个角度的医疗影像，根据医疗影像，通过U-Net网络模型，得到施源器管体的轴位掩膜图像。
需要说明的是，对于前述的各方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本申请并不受所描述的动作顺序的限制，因为依据本申请，某些步骤可以采用其他顺序或者同时进行。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于可选实施例，所涉及的动作和模块并不一定是本申请所必须的。
以上是关于方法实施例的介绍，以下通过装置实施例，对本申请所述方案进行进一步说明。
图4示出了根据本申请的实施例的施源器中间装置的方框图。如图4所示，该装置包括：
第一生成模块401，用于基于施源器管体的轴位掩膜图像，生成所述施源器管体的三维掩膜图像，所述施源器管体至少包括直线部分和弯曲部分.
管体旋转模块402，用于在预设三维坐标系下，对所述三维掩膜图像进行切割、旋转，使得所述直线部分与第一坐标轴重合、所述弯曲部分位于第一平面或第二平面上。
第二生成模块403，用于生成所述施源器管体的冠状位切面和矢状位切面。其中，第一平面由第一坐标轴和第二坐标轴形成，第二平面由第一坐标轴和第三坐标轴形成。
在一些实施例中，管体旋转模块402具体用于：
若所述施源器管体为非中间管体，
在预设三维坐标系下，在p1点处沿第一平面方向切割所述三维掩膜图像；
旋转第一平面，使得所述直线部分的中心线与第一坐标轴重合，并更新预设三维坐标系；
在更新后的预设三维坐标系下，在p1点处沿第二平面方向切割所述三维掩膜图像；
旋转第二平面，使得所述直线部分的中心线与第一坐标轴重合，并再次更新预设三维坐标系；
将所述施源器管体投影至再次更新后的预设三维坐标系的第三平面上；
旋转第三平面，使得投影部分的中心线与第二坐标轴或第三坐标轴重合；
其中，p1点为所述直线部分下端任一点，第三平面由第二坐标轴和第三坐标轴形成。
在一些实施例中，管体旋转模块402具体还用于：
若所述施源器管体为中间管体，
在预设三维坐标系下，在p1点处沿第一平面方向切割所述三维掩膜图像；
旋转第一平面，使得所述直线部分的中心线与第一坐标轴重合，并更新预设三维坐标系；
在更新后的预设三维坐标系下，在p1点处沿第二平面方向切割所述三维掩膜图像；
旋转第二平面，使得所述直线部分的中心线与第一坐标轴重合，并再次更新预设三维坐标系；
在再次更新后的预设三维坐标系下，确定所述直线部分的中心点q1和所述弯曲部分的中心点q2；
旋转第三平面，使得中心点q1和中心点q2的横坐标位于第三坐标轴上，或者使得中心点q1和中心点q2的纵坐标位于第二坐标轴上；
其中，p1点为所述直线部分下端任一点，第三平面由第二坐标轴和第三坐标轴形成。
在一些实施例中，管体旋转模块402具体还用于：
在p1点处沿第三平面方向切割所述三维掩膜图像，得到所述施源器管体的第一轴位图像，基于所述第一轴位图像确定中心点q1；
在p2点处沿第三平面方向切割所述三维掩膜图像，得到所述施源器管体的第二轴位图像，基于所述第二轴位图像确定中心点q2；
其中，p2点为所述弯曲部分的顶点。
在一些实施例中，第一生成模块401具体用于将所述施源器管体的所述轴位掩膜图像进行叠放，生成所述施源器管体的所述三维掩膜图像。
在一些实施例中，该装置还包括第三生成模块，第三生成模块用于获取多个角度的医疗影像；根据所述医疗影像，通过U-Net网络模型，得到所述施源器管体的所述轴位掩膜图像。
在一些实施例中，该装置还包括提取模块，提取模块用于基于骨架提取算法，提取所述施源器管体的初始中心线；延长所述初始中心线的两端，形成所述施源器管体的中心线。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，所述描述的模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
本申请实施例中提供了一种电子设备，如图5所示，图5所示的电子设备500包括：处理器501和存储器503。其中，处理器501和存储器503相连，如通过总线502相连。可选地，电子设备500还可以包括收发器504。需要说明的是，实际应用中收发器504不限于一个，该电子设备500的结构并不构成对本申请实施例的限定。
处理器501可以是CPU(CentralProcessingUnit，中央处理器)，通用处理器，DSP(DigitalSignalProcessor，数据信号处理器)，ASIC(ApplicationSpecificIntegratedCircuit，专用集成电路)，FPGA(FieldProgrammableGateArray，现场可编程门阵列)或者其他可编程逻辑器件、晶体管逻辑器件、硬件部件或者其任意组合。其可以实现或执行结合本申请公开内容所描述的各种示例性的逻辑方框，模块和电路。处理器501也可以是实现计算功能的组合，例如包含一个或多个微处理器组合，DSP和微处理器的组合等。
总线502可包括一通路，在上述组件之间传送信息。总线502可以是PCI(PeripheralComponentInterconnect，外设部件互连标准)总线或EISA(ExtendedIndustryStandardArchitecture，扩展工业标准结构)总线等。总线502可以分为地址总线、数据总线、控制总线等。为便于表示，图5中仅用一条粗线表示，但并不表示仅有一根总线或一种类型的总线。
存储器503可以是ROM(ReadOnlyMemory，只读存储器)或可存储静态信息和指令的其他类型的静态存储设备，RAM(RandomAccessMemory，随机存取存储器)或者可存储信息和指令的其他类型的动态存储设备，也可以是EEPROM(ElectricallyErasableProgrammableReadOnlyMemory，电可擦可编程只读存储器)、CD-ROM(CompactDiscReadOnlyMemory，只读光盘)或其他光盘存储、光碟存储(包括压缩光碟、激光碟、光碟、数字通用光碟、蓝光光碟等)、磁盘存储介质或者其他磁存储设备、或者能够用于携带或存储具有指令或数据结构形式的期望的程序代码并能够由计算机存取的任何其他介质，但不限于此。
存储器503用于存储执行本申请方案的应用程序代码，并由处理器501来控制执行。处理器501用于执行存储器503中存储的应用程序代码，以实现前述方法实施例所示的内容。
其中，电子设备包括但不限于：移动电话、笔记本电脑、数字广播接收器、PDA(个人数字助理)、PAD(平板电脑)、PMP(便携式多媒体播放器)、车载终端(例如车载导航终端)等等的移动终端以及诸如数字TV、台式计算机等等的固定终端。图5示出的电子设备仅仅是一个示例，不应对本申请实施例的功能和使用范围带来任何限制。
本申请实施例提供了一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，当其在计算机上运行时，使得计算机可以执行前述方法实施例中相应内容。
应该理解的是，虽然附图的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，其可以以其他的顺序执行。而且，附图的流程图中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，其执行顺序也不必然是依次进行，而是可以与其他步骤或者其他步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
以上所述仅是本申请的部分实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本申请原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本申请的保护范围。

Claims (9)

1.一种施源器重建方法，其特征在于，包括：
基于施源器管体的轴位掩膜图像，生成所述施源器管体的三维掩膜图像，所述施源器管体至少包括直线部分和弯曲部分；
在预设三维坐标系下，对所述三维掩膜图像进行切割、旋转，使得所述直线部分与第一坐标轴重合、所述弯曲部分位于第一平面或第二平面上；
生成所述施源器管体的冠状位切面和矢状位切面；
其中，第一平面由第一坐标轴和第二坐标轴形成，第二平面由第一坐标轴和第三坐标轴形成；
所述在预设三维坐标系下，对所述三维掩膜图像进行切割、旋转，包括：
若所述施源器管体为中间管体，在预设三维坐标系下，在 p1 点处沿第一平面方向切割所述三维掩膜图像；
旋转第一平面，使得所述直线部分的中心线与第一坐标轴重合，并更新预设三维坐标系；
在更新后的预设三维坐标系下，在 p1 点处沿第二平面方向切割所述三维掩膜图像；
旋转第二平面，使得所述直线部分的中心线与第一坐标轴重合，并再次更新预设三维坐标系；
在再次更新后的预设三维坐标系下，确定所述直线部分的中心点 q1 和所述弯曲部分的中心点 q2；
旋转第三平面，使得中心点 q1 和中心点 q2 的横坐标位于第三坐标轴上，或者使得中心点 q1 和中心点 q2 的纵坐标位于第二坐标轴上；
其中，p1 点为所述直线部分下端任一点，第三平面由第二坐标轴和第三坐标轴形成。
2.根据权利要求 1 所述的方法，其特征在于，所述在预设三维坐标系，对所述三维掩膜图像进行切割、旋转，包括：
若所述施源器管体为非中间管体，在预设三维坐标系下，在 p1 点处沿第一平面方向切割所述三维掩膜图像；
旋转第一平面，使得所述直线部分的中心线与第一坐标轴重合，并更新预设三维坐标系；
在更新后的预设三维坐标系下，在 p1 点处沿第二平面方向切割所述三维掩膜图像；
旋转第二平面，使得所述直线部分的中心线与第一坐标轴重合，并再次更新预设三维坐标系；
将所述施源器管体投影至再次更新后的预设三维坐标系的第三平面上；
旋转第三平面，使得投影部分的中心线与第二坐标轴或第三坐标轴重合；其中，p1 点为所述直线部分下端任一点，第三平面由第二坐标轴和第三坐标轴形成。
3.根据权利要求1所述的方法，其特征在于，确定所述直线部分的中心点 q1 和所述弯曲部分的中心点 q2，包括：
在 p1 点处沿第三平面方向切割所述三维掩膜图像，得到所述施源器管体的第一轴位图像，基于所述第一轴位图像确定中心点 q1；
在 p2 点处沿第三平面方向切割所述三维掩膜图像，得到所述施源器管体的第二轴位图像，基于所述第二轴位图像确定中心点 q2；
其中，p2 点为所述弯曲部分的顶点。
4.根据权利要求 1 所述的方法，其特征在于，基于施源器管体的轴位掩膜图像，生成所述施源器管体的三维掩膜图像，包括：
将所述施源器管体的所述轴位掩膜图像进行叠放，生成所述施源器管体的所述三维掩膜图像。
5.根据权利要求 1 所述的方法，其特征在于，所述基于施源器管体的轴位掩膜图像，生成所述施源器管体的三维掩膜图像之前，还包括：
获取多个角度的医疗影像；
根据所述医疗影像，通过 U-Net 网络模型，得到所述施源器管体的所述轴位掩膜图像。
6.根据权利要求 1 所述的方法，其特征在于，还包括：
基于骨架提取算法，提取所述施源器管体的初始中心线；
延长所述初始中心线的两端，形成所述施源器管体的中心线。
7.一种施源器重建装置，其特征在于，包括：
第一生成模块，用于基于施源器管体的轴位掩膜图像，生成所述施源器管体的三维掩膜图像，所述施源器管体至少包括直线部分和弯曲部分；
管体旋转模块，用于在预设三维坐标系下，对所述三维掩膜图像进行切 割、旋转，使得所述直线部分与第一坐标轴重合、所述弯曲部分位于第一平面或第二平面上；
第二生成模块，用于生成所述施源器管体的冠状位切面和矢状位切面；
其中，第一平面由第一坐标轴和第二坐标轴形成，第二平面由第一坐标轴和第三坐标轴形成；
管体旋转模块，还用于若所述施源器管体为中间管体，在预设三维坐标系下，在 p1 点处沿第一平面方向切割所述三维掩膜图像；旋转第一平面，使得所述直线部分的中心线与第一坐标轴重合，并更新预设三维坐标系；在更新后的预设三维坐标系下，在 p1 点处沿第二平面方向切割所述三维掩膜图像；旋转第二平面，使得所述直线部分的中心线与第一坐标轴重合，并再次更新预设三维坐标系；在再次更新后的预设三维坐标系下，确定所述直线部分的中心点 q1 和所述弯曲部分的中心点 q2；旋转第三平面，使得中心点 q1 和中心点q2 的横坐标位于第三坐标轴上，或者使得中心点 q1 和中心点 q2 的纵坐标位于第二坐标轴上；其中，p1 点为所述直线部分下端任一点，第三平面由第二坐标轴和第三坐标轴形成。
8.一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求 1 至 6中任一项所述的方法。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求 1 至 6 中任一项所述的方法。
CN202310103177.3A 2023-02-13 2023-02-13 施源器重建方法、装置、电子设备和计算机可读存储介质 Active CN116152437B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310103177.3A CN116152437B (zh) 2023-02-13 2023-02-13 施源器重建方法、装置、电子设备和计算机可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310103177.3A CN116152437B (zh) 2023-02-13 2023-02-13 施源器重建方法、装置、电子设备和计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN116152437A CN116152437A (zh) 2023-05-23
CN116152437B true CN116152437B (zh) 2024-02-02

Family

ID=86357837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310103177.3A Active CN116152437B (zh) 2023-02-13 2023-02-13 施源器重建方法、装置、电子设备和计算机可读存储介质

Country Status (1)

Country Link
CN (1) CN116152437B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117274833B (zh) * 2023-11-20 2024-02-27 浙江国遥地理信息技术有限公司 一种建筑轮廓处理方法、装置、设备和存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018058291A1 (zh) * 2016-09-27 2018-04-05 深圳先进技术研究院 基于磁共振成像的施源位置定位方法和施源器
CN110675481A (zh) * 2019-09-26 2020-01-10 慧影医疗科技（北京）有限公司 影像交互联动方法、装置及可读存储介质
CN111353339A (zh) * 2018-12-21 2020-06-30 厦门歌乐电子企业有限公司 一种物体识别装置及方法
CN113822923A (zh) * 2021-09-22 2021-12-21 柏意慧心(杭州)网络科技有限公司 用于获取血管的目标切面图像的方法、设备和介质
CN115359143A (zh) * 2022-09-07 2022-11-18 浙江大学医学院附属妇产科医院 一种施源器自动重建方法和装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015181632A1 (en) * 2014-05-28 2015-12-03 Nucletron Operations B.V. Methods and systems for brachytherapy planning based on imaging data

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018058291A1 (zh) * 2016-09-27 2018-04-05 深圳先进技术研究院 基于磁共振成像的施源位置定位方法和施源器
CN111353339A (zh) * 2018-12-21 2020-06-30 厦门歌乐电子企业有限公司 一种物体识别装置及方法
CN110675481A (zh) * 2019-09-26 2020-01-10 慧影医疗科技（北京）有限公司 影像交互联动方法、装置及可读存储介质
CN113822923A (zh) * 2021-09-22 2021-12-21 柏意慧心(杭州)网络科技有限公司 用于获取血管的目标切面图像的方法、设备和介质
CN115359143A (zh) * 2022-09-07 2022-11-18 浙江大学医学院附属妇产科医院 一种施源器自动重建方法和装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Comparing multichannel cylinder and 3D-printed applicators for vaginal cuff brachytherapy with preliminary exploration of post-hysterectomy vaginal morphology;Junfang Yan等;Journal of Contemporary Brachytherapy;第13卷(第6期);641–648 *
基于CT图像的三维后装腔内治疗计划系统的研制;余辉;徐海荣;张书旭;石玉生;钱剑杨;;中国医疗器械杂志(第02期);93-96+107 *
多通道柱状施源器端到端测试;于浪等;中国医疗设备;第36卷(第4期);104-107+126 *

Also Published As

Publication number Publication date
CN116152437A (zh) 2023-05-23

Similar Documents

Publication Publication Date Title
CN108509915B (zh) 人脸识别模型的生成方法和装置
US11657509B2 (en) Method for precisely and automatically positioning reference line for integrated images
WO2021189843A1 (zh) Ct图像椎骨定位方法、装置、设备及介质
CN116152437B (zh) 施源器重建方法、装置、电子设备和计算机可读存储介质
US20210027526A1 (en) Lighting estimation
US20150348229A1 (en) Image annotation using a haptic plane
CN111149128B (zh) 标记数据的获取方法及装置、训练方法及装置、医疗设备
CN115359063B (zh) 基于目标器官三维图像的刚性配准方法及相关装置
CN111445550B (zh) Pet图像的迭代重建方法、装置和计算机可读存储介质
CN113888566B (zh) 目标轮廓曲线确定方法、装置、电子设备以及存储介质
US20210390792A1 (en) Method and device for customizing facial expressions of user
EP3843038B1 (en) Image processing method and system
US12020351B2 (en) Method, device and equipment for reconstructing PET images
Zhang et al. Generalizable and promptable artificial intelligence model to augment clinical delineation in radiation oncology
US11232612B2 (en) Highly accurate and efficient forward and back projection methods for computed tomography
CN113262048B (zh) 空间配准方法、装置、终端设备及术中导航系统
CN115909035A (zh) 手术器械的清点方法、系统、设备及介质
Eulig et al. Deep learning-aided CBCT image reconstruction of interventional material from four x-ray projections
Datteri et al. Applying the algorithm" assessing quality using image registration circuits"(aquirc) to multi-atlas segmentation
Chegeni et al. Introduction of a simple algorithm to create synthetic-Computed tomography of the head from magnetic resonance imaging
CN111368832B (zh) 感兴趣区域的标记方法、装置、设备和存储介质
Kim et al. Clinical feasibility of deep learning-based synthetic CT images from T2-weighted MR images for cervical cancer patients compared to MRCAT
CN116309924B (zh) 一种模型训练的方法、图像展示的方法及装置
EP4036865A1 (en) Method, device and system for visualizing colon
Das et al. Co-ordinate-based positional embedding that captures resolution to enhance transformer’s performance in medical image analysis

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230711

Address after: 100010 Shuaifuyuan No. 1, Wangfujing, Dongcheng District, Beijing

Applicant after: PEKING UNION MEDICAL COLLEGE Hospital

Address before: Room 02 and 03, 5 / F, block a, Tiangong building, Keda, No. 30 Xueyuan Road, Haidian District, Beijing 100089

Applicant before: Beijing yizhiying Technology Co.,Ltd.

GR01 Patent grant
GR01 Patent grant