CN115082649B - 一种肝脏分段方法、装置、电子设备及存储介质 - Google Patents

一种肝脏分段方法、装置、电子设备及存储介质 Download PDF

Info

Publication number
CN115082649B
CN115082649B CN202210979778.6A CN202210979778A CN115082649B CN 115082649 B CN115082649 B CN 115082649B CN 202210979778 A CN202210979778 A CN 202210979778A CN 115082649 B CN115082649 B CN 115082649B
Authority
CN
China
Prior art keywords
liver
portal vein
temperature
blood supply
radius
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
CN202210979778.6A
Other languages
English (en)
Other versions
CN115082649A (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.)
Guangdong Oupu Mandi Technology Co ltd
Original Assignee
Guangdong Optomedic Technology Co Ltd
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 Guangdong Optomedic Technology Co Ltd filed Critical Guangdong Optomedic Technology Co Ltd
Priority to CN202210979778.6A priority Critical patent/CN115082649B/zh
Publication of CN115082649A publication Critical patent/CN115082649A/zh
Application granted granted Critical
Publication of CN115082649B publication Critical patent/CN115082649B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • 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/10—Segmentation; Edge detection
    • G06T7/162—Segmentation; Edge detection involving graph-based methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
    • 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/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10081—Computed x-ray tomography [CT]
    • 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/20081—Training; Learning
    • 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]
    • 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/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30056—Liver; Hepatic
    • 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/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30096—Tumor; Lesion
    • 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/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/41—Medical

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Graphics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Architecture (AREA)
  • Computer Hardware Design (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本申请提供了一种肝脏分段方法、装置、电子设备及存储介质，涉及医疗技术领域，其技术方案要点是：获取各门静脉的半径以及所述各门静脉在肝脏中的位置；根据所述各门静脉的半径以及所述各门静脉在所述肝脏中的位置确定所述各门静脉的供血区域；根据所述供血区域对所述肝脏进行分段。本申请提供的一种肝脏分段方法、装置、电子设备及存储介质具有对不同个体的肝脏分段准确的优点。

Description

一种肝脏分段方法、装置、电子设备及存储介质
技术领域
本申请涉及医疗技术领域，具体而言，涉及一种肝脏分段方法、装置、电子设备及存储介质。
背景技术
CT（Computed Tomography）扫描是临床上常用的检查手段，随着精准外科手术的发展，在肝胆外科领域，基于CT影像的术前规划手术系统已经越来越成熟。借助于三维重建方法和肝脏分段技术，医生可以直观地分析肝脏组织中不同血管之间的相对关系以及肿瘤在肝脏中的三维位置，从而作出详细的手术规划，在术中尽可能减少对健康组织的切除，减少手术创伤和术中出血。
现有的肝脏分段方法，主要为Couinaud分段法，根据门静脉的分布特点将肝脏分为1-8段，每段认为是一个独立功能区，有独立血流供应系统。但这种固定的经验分段方法未考虑个体化差异，不同患者门脉分布特点不一致和肝脏变异均会导致上述方法分段不准确。
此外，现有的个性化交互式分段方法，主要是人工定义不同段对应的门静脉分支，然后通过迭代计算每个肝脏单元与标记血管的距离进行自动分段，主要方式还是根据肝脏单元与选中静脉之间的距离进行简单的分段，将距离目标肝脏单元最近的门静脉类别赋予目标肝脏单元，未考虑门静脉血管真实形态影响，因此上述方法并不准确。
针对上述问题，本申请提出一种新的解决方案。
发明内容
本申请的目的在于提供一种肝脏分段方法、装置、电子设备及存储介质，具有对不同个体的肝脏分段准确的优点。
第一方面，本申请提供了一种肝脏分段方法，技术方案如下：
包括：
获取各门静脉的半径以及所述各门静脉在肝脏中的位置；
根据所述各门静脉的半径以及所述各门静脉在所述肝脏中的位置确定所述各门静脉的供血区域；
根据所述供血区域对所述肝脏进行分段。
通过各门静脉的半径以及各门静脉在肝脏中的位置来确定供血区域，根据供血区域来对肝脏进行分段，具有对不同个体的肝脏分段准确的有益效果。
进一步地，在本申请中，所述根据所述各门静脉的半径以及所述各门静脉在所述肝脏中的位置确定所述各门静脉的供血区域的步骤包括：
根据所述各门静脉的半径设置初始温度；
根据所述各门静脉在所述肝脏中的位置，分别以所述各门静脉作为热源对应计算所述肝脏各个位置在任意时刻的传递温度；
根据所述肝脏各个位置在任意时刻的所述传递温度确定所述各门静脉的供血区域。
进一步地，在本申请中，所述根据所述肝脏各个位置在任意时刻的所述传递温度确定所述各门静脉的供血区域的步骤包括：
根据所述肝脏各个位置在任意时刻的所述传递温度获取在经过相同热传递时间后，所述各门静脉对所述肝脏任意同一位置的所述传递温度；
根据所述肝脏任意同一位置的所述传递温度得到最大传递温度；
将所述肝脏任意同一位置作为所述最大传递温度对应的所述门静脉的供血区域；
遍历计算所述肝脏的所有位置以确定所述各门静脉的供血区域。
进一步地，在本申请中，所述分别以所述各门静脉作为热源对应计算所述肝脏各个位置在任意时刻的传递温度的公式为：
Figure 323969DEST_PATH_IMAGE001
；
Figure 100002_DEST_PATH_IMAGE002
；
Figure 858855DEST_PATH_IMAGE003
；
Figure 100002_DEST_PATH_IMAGE004
；
Figure 283889DEST_PATH_IMAGE005
；
Figure 100002_DEST_PATH_IMAGE006
表示的是所述各门静脉标号；
Figure 271437DEST_PATH_IMAGE007
表示的是以所述各门静脉为热源，所述肝脏任意位置在任意时刻下的所述传递温度的温度值；
Figure 100002_DEST_PATH_IMAGE008
表示的是所述肝脏任意位置的坐标值；
Figure 574373DEST_PATH_IMAGE009
表示的是任意时刻；
Figure 104712DEST_PATH_IMAGE010
表示的是热传递系数；
Figure 982538DEST_PATH_IMAGE011
表示的是圆周率；
Figure 863906DEST_PATH_IMAGE012
表示的是所述门静脉的半径；
Figure 740507DEST_PATH_IMAGE013
表示的是空间温度，
Figure 390932DEST_PATH_IMAGE014
表示的是空间内的坐标值；
Figure 439659DEST_PATH_IMAGE015
表示的是自然常数；
Figure 808323DEST_PATH_IMAGE016
表示的是空间中一点的温度对时间的变化率；
Figure 718642DEST_PATH_IMAGE017
表示的是拉普拉斯算子。
进一步地，在本申请中，所述将所述肝脏任意同一位置作为所述最大传递温度对应的所述门静脉的供血区域的公式为：
Figure 754731DEST_PATH_IMAGE018
；
Figure 584147DEST_PATH_IMAGE019
表示的是所述肝脏任意位置所属的所述门静脉的供血区域。
进一步地，在本申请中，根据所述各门静脉的半径设置初始温度的公式为：
Figure 564741DEST_PATH_IMAGE020
；
Figure 403384DEST_PATH_IMAGE021
表示的是所述肝脏任意位置的初始温度；
Figure 136723DEST_PATH_IMAGE022
表示的是在
Figure 402619DEST_PATH_IMAGE023
处所述门静脉的半径数值，将该数值作为初始温度。
进一步地，在本申请中，还包括：
获取CT数据；
根据所述CT数据以及对所述肝脏进行分段的分段结果进行三维重建，得到经过分段后的三维模型。
第二方面，本申请还提供一种肝脏分段装置，包括：
获取模块，用于获取各门静脉的半径以及所述各门静脉在肝脏中的位置；
计算模块，用于根据所述各门静脉的半径以及所述各门静脉在所述肝脏中的位置确定所述各门静脉的供血区域；
分段模块，用于根据所述供血区域对所述肝脏进行分段。
第三方面，本申请还提供一种电子设备，包括处理器以及存储器，所述存储器存储有计算机可读取指令，当所述计算机可读取指令由所述处理器执行时，运行如上所述方法中的步骤。
第四方面，本申请还提供一种存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时，运行如上所述方法中的步骤。
由上可知，本申请提供的一种肝脏分段方法、装置、电子设备及存储介质，通过各门静脉的半径以及各门静脉在肝脏中的位置来确定供血区域，根据供血区域来对肝脏进行分段，具有对不同个体的肝脏分段准确的有益效果。
附图说明
图1为本申请提供的一种肝脏分段方法的流程图。
图2为本申请提供的一种肝脏分段装置的结构示意图。
图3为本申请提供的一种电子设备的结构示意图。
图中：210、获取模块；220、计算模块；230、分段模块；310、处理器；320、存储器。
具体实施方式
下面将结合本申请中附图，对本申请中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本申请的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围，而是仅仅表示本申请的选定实施例。基于本申请的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本申请的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
目前，现有的肝脏分段方法最常用的是Couinaud分段法，根据门静脉的分布特点将肝脏分为1-8段，这种分段方法没有考虑个体差异，因此，对肝脏的分段并不准确，随着科学技术的发展，目前也有个性化交互式的分段方法，但是现有的个性化交互式分段方法需要通过迭代计算每个肝脏与门静脉的之间的距离进行分段，一般会获取门静脉的中心线，根据肝脏单元与门静脉的中心线的距离进行分段，这种分段方法没有考虑到血管的真实形态，因此同样存在分段不准确的问题。
对此，请参照图1，本申请提出了一种肝脏分段方法，技术方案具体包括：
S110、获取各门静脉的半径以及各门静脉在肝脏中的位置；
S120、根据各门静脉的半径以及各门静脉在肝脏中的位置确定各门静脉的供血区域；
S130、根据供血区域对肝脏进行分段。
通过上述方案，本申请根据各门静脉的半径以及各门静脉在肝脏中的位置来确定供血区域，根据供血区域来对肝脏进行分段，而不是简单的根据肝脏单元与门静脉的距离来进行分段，由于不同半径的门静脉的供血能力会有差别，即使是同一根门静脉，在不同的个体上，其半径以及在肝脏中的分布位置、走向也会存在区别，因此，本申请根据门静脉的半径以及其在肝脏中的位置来确定对应的供血区域，然后根据供血区域对肝脏进行分段，具有对不同个体的肝脏分段准确的有益效果。
本申请方案的核心在于通过门静脉的半径以及对应在肝脏中的位置来确定供血区域，然后根据供血区域来对肝脏进行分段，然而，如何根据门静脉的半径以及对应在肝脏中的位置来确定供血区域将成为一个重要的问题，在本申请公开上述方案的基础上，可以采用计算肝脏单元到门静脉外壁的距离，根据距离来进行分段，但是，这种方式的计算复杂，难以应用在实际中。
对此，本申请提出了一种全新的思路，采用热传递方程柯西问题求解方法来计算各门静脉对应的供血区域，将计算门静脉的供血区域的问题转换为热传递的计算问题，将血管作为热源，肝脏组织作为热传递介质，可以计算出门静脉在不同肝脏区域的热传递情况，通过比较不同门静脉某个特定时刻在肝脏区域内的热传递情况，就可以判断出不同门静脉的供血区域。
具体的，在其中一些可以实现的实施例中，根据各门静脉的半径以及各门静脉在肝脏中的位置确定各门静脉的供血区域的步骤包括：
根据各门静脉的半径设置初始温度；
根据各门静脉在肝脏中的位置，分别以各门静脉作为热源对应计算肝脏各个位置在任意时刻的传递温度；
根据肝脏各个位置在任意时刻的传递温度确定各门静脉的供血区域。
通过门静脉的半径来设置初始温度，因为不同半径的门静脉的供血能力不同，因此，不同半径的门静脉的热传递能力也不同，所以，根据门静脉的半径设定设置初始温度，不同的半径对应有不同的设置初始温度。
具体的，设置初始温度可以用以下公式表达：
Figure 100002_DEST_PATH_IMAGE024
；
Figure 463985DEST_PATH_IMAGE025
表示的是肝脏任意位置的初始温度。
即，以其中一根门静脉作为热源的时候，该门静脉内的初始温度标记为R，
Figure 480220DEST_PATH_IMAGE026
表示的是在
Figure 959743DEST_PATH_IMAGE023
位置处门静脉的截面半径，肝脏的其余区域的初始温度则标记为0。
然后，分别以每根门静脉作为热源，对应计算在肝脏区域内的热传递情况，具体的计算公式为：
Figure 521174DEST_PATH_IMAGE001
；
Figure 86148DEST_PATH_IMAGE002
；
Figure 532172DEST_PATH_IMAGE003
；
Figure 741568DEST_PATH_IMAGE004
；
Figure 349267DEST_PATH_IMAGE005
；
Figure 526170DEST_PATH_IMAGE027
表示的是各门静脉标号；
Figure 510307DEST_PATH_IMAGE028
表示的是以各门静脉为热源，肝脏任意位置在任意时刻下的传递温度的温度值，是时间变数
Figure DEST_PATH_IMAGE029
与空间变数
Figure 801305DEST_PATH_IMAGE023
的函数；
Figure 704539DEST_PATH_IMAGE023
表示的是肝脏任意位置的坐标值；
Figure 244105DEST_PATH_IMAGE029
表示的是任意时刻；
Figure 641720DEST_PATH_IMAGE030
表示的是热传递系数；
Figure 950341DEST_PATH_IMAGE031
表示的是圆周率；
Figure 165422DEST_PATH_IMAGE032
表示的是门静脉的半径；
Figure 51338DEST_PATH_IMAGE033
表示的是空间温度，为初始温度，
Figure 377277DEST_PATH_IMAGE034
表示的是空间内的坐标值，即
Figure 179886DEST_PATH_IMAGE033
是在坐标
Figure 300289DEST_PATH_IMAGE034
下的温度，值得注意的是，
Figure 939081DEST_PATH_IMAGE034
是计算过程中的一个积分单元；
Figure 803131DEST_PATH_IMAGE035
表示的是自然常数；
Figure 696132DEST_PATH_IMAGE036
表示的是空间中一点的温度对时间的变化率；
Figure 518595DEST_PATH_IMAGE037
表示的是拉普拉斯算子。
通过上述技术方案，将任意肝脏单元所属门静脉供血区域的问题转换成热传递的计算问题，分别以各根门静脉作为热源，计算门静脉在不同肝脏区域的热传递情况，通过比较不同的热传递情况，就可以确定肝脏单元属于哪一根门静脉的供血区域。
在一些优选方案中，在计算求解的时候，对于将不同的门静脉作为热源进行热传递，肝脏单元的热传递系数是一致的，因此，可以将热传递系数
Figure 254469DEST_PATH_IMAGE030
统一设置为1，以简化计算过程，并且，在求解的过程中，需要比较在同一时刻下，各门静脉对肝脏任意同一位置的传递温度，因此，也可以将
Figure 46845DEST_PATH_IMAGE029
统一设置为1，将上述计算在肝脏区域内的热传递情况的公式离散化，得到公式：
Figure 184565DEST_PATH_IMAGE038
其中，
Figure 286251DEST_PATH_IMAGE039
为求解范围，在这里表示的是肝脏区域。
进一步地，在其中一些可以实现的实施例中，根据肝脏各个位置在任意时刻的传递温度确定各门静脉的供血区域的步骤包括：
根据肝脏各个位置在任意时刻的传递温度获取在经过相同热传递时间后，各门静脉对肝脏任意同一位置的传递温度；
根据肝脏任意同一位置的传递温度得到最大传递温度；
将肝脏任意同一位置作为最大传递温度对应的门静脉的供血区域；
遍历计算肝脏的所有位置以确定各门静脉的供血区域。
具体的，将肝脏任意同一位置作为最大传递温度对应的门静脉的供血区域的公式为：
Figure 509422DEST_PATH_IMAGE040
；
Figure 839909DEST_PATH_IMAGE019
表示的是肝脏任意位置所属的门静脉的供血区域。
即，分别以不同的门静脉作为热源，在经过相同进行热传递的时间后，计算肝脏任意的同一位置的传递温度，可以得到与不同的门静脉对应的多个传递温度，然后比较多个传递温度，得到最大传递温度，即，该任意的同一位置则属于该最大传递温度所对应的门静脉的供血区域，通过计算热传递的方式计算出肝脏中任意位置所属的门静脉，可以完成每根门静脉的供血区域的计算。
进一步地，在其中一些可以实现的实施例中，还包括：
获取CT数据；
根据CT数据以及对肝脏进行分段的分段结果进行三维重建，得到经过分段后的三维模型。
通过上述技术方案，根据CT数据以及对肝脏进行分段的分段结果进行三维重建，得到经过分段后的三维模型，进而更好的完成术前规划。
综上，本申请提出了一种新的肝脏分段方法，可以对不同个体的门静脉进行标记，再根据门静脉的位置以及半径，采用热传递方程柯西问题求解方法来实现肝脏的分段，不仅考虑到门静脉在肝脏中的位置分布，同时还考虑到血管半径对供血分布的影响，对血管的供血区域的分区更加准确。
在实际应用过程中，首先是基于CT数据完成肝脏、门静脉等组织的三维重建，一般情况下，对于肝脏、肿瘤、下腔静脉以及门静脉组织，均是基于门脉期或者静脉期图像进行组织分割和重建。
在读入CT数据后，对于肝脏等大组织，可采用图割法完成肝脏区域图像分割，对于下腔静脉和门静脉，由于结构较为复杂，传统方法无法获取较好的分割图像，可以采用基于深度学习的nnUnet完成门静脉和下腔静脉的分割，具体流程为：预先获取足够数量的肝脏静脉期数据，例如100套数据，对目标组织进行标注，制作训练集；然后采用nnUNet对训练集进行训练，得到训练模型；最后采用训练模型对新数据进行推理，完成目标组织分割。
完成目标组织分割后，采用Marchingcube方法对分割结果进行三维重建，得到目标组织的三维模型，并采用VTK进行显示和交互操作。
在得到目标组织的三维模型后，对肝脏进行交互分段。
肝脏分段的目标是根据标记的血管，将肝脏分成不同标记的肝脏单元，每个肝脏单元主要受匹配标记血管影响，最终得到标记血管的供血区域，从而优化手术切割路径。个性化的肝脏分段充分考虑了不同个体的肝脏、门静脉解剖性差异，能够很好地根据肿瘤、血管解剖位置来确定手术路径。对于相邻两根不同标记的不同粗细血管，目前主要采用距离变换方法，根据肝脏单元与血管的欧氏距离来确定标记所属，未考虑血管半径对组织供血的情况。一般情况下，血管半径越大，所供血组织区域越大，反之越小，因此在分段过程中，应同时将血管半径、组织与血管之间的距离考虑进去。
在本申请的方案中，先通过交互方式对门静脉三维模型的各个血管分支进行标记，即标记为不同的门静脉，完成标记后，通过标记结果计算出每个肝脏单元所属标记，从而完成肝脏三维模型的分段，每一段对应一个标记。
肝脏单元所属标记的计算即为各门静脉供血区域的计算，其过程与热传递计算类似，因此，可以将不同标记的门静脉分别作为热源，肝脏组织作为热传递介质，则可以计算出不同标记的门静脉在不同肝脏区域的热传递情况，通过比较不同标记的门静脉某个特定时刻在肝脏区域内的热传递情况，则可以判断出任意肝脏单元所属的门静脉。
以门静脉为热源，肝脏组织为热传递介质的热传递过程符合热传递方程柯西问题求解情况，即：
Figure 97715DEST_PATH_IMAGE041
；
Figure 996401DEST_PATH_IMAGE042
；
其中，
Figure 316655DEST_PATH_IMAGE043
，是任意时刻下任意位置的温度值；
Figure 326200DEST_PATH_IMAGE030
为热传递系数；
Figure 563146DEST_PATH_IMAGE044
是
Figure 367154DEST_PATH_IMAGE045
时刻的空间温度，为初始温度。
上述方程的解析解为：
Figure 938818DEST_PATH_IMAGE046
；
假设标记的门静脉的标号为1，2，3…n，则通过上述方程可以计算出对应的
Figure 752053DEST_PATH_IMAGE047
、
Figure 843506DEST_PATH_IMAGE048
、
Figure 818416DEST_PATH_IMAGE049
…
Figure 378841DEST_PATH_IMAGE050
，将
Figure 995767DEST_PATH_IMAGE051
固定到某一共同时刻，则可以求出不同标记的门静脉对于肝脏任意位置的分布情况：
Figure 817093DEST_PATH_IMAGE052
，
Figure 353116DEST_PATH_IMAGE053
。对于任意的肝脏位置
Figure 525472DEST_PATH_IMAGE054
，其所属标记的计算公式为：
Figure 794691DEST_PATH_IMAGE055
，
Figure 736102DEST_PATH_IMAGE004
。
对于被标记的门静脉而言，存在初始温度，相对的，其它位置的温度可以设置为0，而被标记的门静脉的初始温度由其血管半径决定，因此，初始温度的设定为：
Figure 177448DEST_PATH_IMAGE024
；
对于不同标记的门静脉，肝脏单元的热传递系数可以认为是一致的，因此，可以将热传递系数设置为1，另外，在以不同的门静脉作为热源进行计算的时候，最好保证时间
Figure 650148DEST_PATH_IMAGE051
相等，因此，可以将
Figure 202353DEST_PATH_IMAGE051
设置为1，将上述计算公式离散化得到：
Figure 998270DEST_PATH_IMAGE056
；
对于
Figure 17042DEST_PATH_IMAGE057
的时刻，任意肝脏位置
Figure 272312DEST_PATH_IMAGE058
所属的计算公式为：
Figure 503573DEST_PATH_IMAGE059
，
Figure 278631DEST_PATH_IMAGE004
。
根据上述方法完成肝脏任意位置所属的计算后，可以对每个标记下的肝脏通过Marchingcube的方法生成相应的三维模型，每个标记对应一段，将所有标记的肝脏组织生成三维模型后，即完成了肝脏分段重建的过程。
本申请提供的一种肝脏分段方法，通过门静脉的半径以及门静脉在肝脏中的位置，基于热传递的原理求出不同门静脉的供血区域，求出的供血区域更符合实际解剖情况，有非常高的准确度。
同时，通过引入热传递方程柯西问题求解过程，使得求解区域符合热传递模型，保证求解区域的平滑性，并且，作为解析方法，无需进行迭代运算，可以并行高速的进行计算处理，提高了复杂血管肝脏分段的计算速度。
第二方面，参照图2，本申请还提供一种肝脏分段装置，包括：
获取模块210，用于获取各门静脉的半径以及各门静脉在肝脏中的位置；
计算模块220，用于根据各门静脉的半径以及各门静脉在肝脏中的位置确定各门静脉的供血区域；
分段模块230，用于根据供血区域对肝脏进行分段。
通过上述技术方案，本申请根据各门静脉的半径以及各门静脉在肝脏中的位置来确定供血区域，根据供血区域来对肝脏进行分段，而不是简单的根据肝脏单元与门静脉的距离来进行分段，由于不同半径的门静脉的供血能力会有差别，即使是同一根门静脉，在不同的个体上，其半径以及在肝脏中的分布位置、走向也会存在区别，因此，本申请根据门静脉的半径以及其在肝脏中的位置来确定对应的供血区域，然后根据供血区域对肝脏进行分段，具有对不同个体的肝脏分段准确的有益效果。
在一些优选的实施例中，本申请提供的一种肝脏分段装置可以执行上述一种肝脏分段方法中的任意一项步骤。
第三方面，参照图3，本申请还提供一种电子设备，包括处理器310以及存储器320，存储器320存储有计算机可读取指令，当计算机可读取指令由处理器310执行时，运行如上方法中的步骤。
通过上述技术方案，处理器310和存储器320通过通信总线和/或其他形式的连接机构（未标出）互连并相互通讯，存储器320存储有处理器310可执行的计算机程序，当电子设备运行时，处理器310执行该计算机程序，以执行时执行上述实施例的任一可选的实现方式中的方法，以实现以下功能：获取各门静脉的半径以及各门静脉在肝脏中的位置；根据各门静脉的半径以及各门静脉在肝脏中的位置确定各门静脉的供血区域；根据供血区域对肝脏进行分段。
第四方面，本申请还提供一种存储介质，其上存储有计算机程序，计算机程序被处理器执行时，运行如上方法中的步骤。
通过上述技术方案，计算机程序被处理器执行时，执行上述实施例的任一可选的实现方式中的方法，以实现以下功能：获取各门静脉的半径以及各门静脉在肝脏中的位置；根据各门静脉的半径以及各门静脉在肝脏中的位置确定各门静脉的供血区域；根据供血区域对肝脏进行分段。
其中，存储介质可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器（Static Random Access Memory, 简称SRAM），电可擦除可编程只读存储器（Electrically Erasable Programmable Read-Only Memory, 简称EEPROM），可擦除可编程只读存储器（Erasable Programmable Read Only Memory, 简称EPROM），可编程只读存储器（Programmable Red-Only Memory, 简称PROM），只读存储器（Read-OnlyMemory, 简称ROM），磁存储器，快闪存储器，磁盘或光盘。
在本申请所提供的实施例中，应该理解到，所揭露装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
另外，作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
再者，在本申请各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
以上所述仅为本申请的实施例而已，并不用于限制本申请的保护范围，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (9)

1.一种肝脏分段方法，其特征在于，包括：
获取各门静脉的半径以及所述各门静脉在肝脏中的位置；
根据所述各门静脉的半径以及所述各门静脉在所述肝脏中的位置确定所述各门静脉的供血区域；
根据所述供血区域对所述肝脏进行分段；
所述根据所述各门静脉的半径以及所述各门静脉在所述肝脏中的位置确定所述各门静脉的供血区域的步骤包括：
根据所述各门静脉的半径设置初始温度；
根据所述各门静脉在所述肝脏中的位置，分别以所述各门静脉作为热源对应计算所述肝脏各个位置在任意时刻的传递温度；
根据所述肝脏各个位置在任意时刻的所述传递温度确定所述各门静脉的供血区域。
2.根据权利要求1所述的一种肝脏分段方法，其特征在于，所述根据所述肝脏各个位置在任意时刻的所述传递温度确定所述各门静脉的供血区域的步骤包括：
根据所述肝脏各个位置在任意时刻的所述传递温度获取在经过相同热传递时间后，所述各门静脉对所述肝脏任意同一位置的所述传递温度；
根据所述肝脏任意同一位置的所述传递温度得到最大传递温度；
将所述肝脏任意同一位置作为所述最大传递温度对应的所述门静脉的供血区域；
遍历计算所述肝脏的所有位置以确定所述各门静脉的供血区域。
3.根据权利要求1所述的一种肝脏分段方法，其特征在于，所述分别以所述各门静脉作为热源对应计算所述肝脏各个位置在任意时刻的传递温度的公式为：
Figure DEST_PATH_IMAGE002
；
Figure DEST_PATH_IMAGE004
；
Figure DEST_PATH_IMAGE006
；
Figure DEST_PATH_IMAGE008
；
Figure DEST_PATH_IMAGE010
；
Figure DEST_PATH_IMAGE012
表示的是所述各门静脉标号；
Figure DEST_PATH_IMAGE014
表示的是以所述各门静脉为热源，所述肝脏任意位置在任意时刻下的所述传递温度的温度值；
Figure DEST_PATH_IMAGE016
表示的是所述肝脏任意位置的坐标值；
Figure DEST_PATH_IMAGE018
表示的是任意时刻；
Figure DEST_PATH_IMAGE020
表示的是热传递系数；
Figure DEST_PATH_IMAGE022
表示的是圆周率；
Figure DEST_PATH_IMAGE024
表示的是所述门静脉的半径；
Figure DEST_PATH_IMAGE026
表示的是空间温度，
Figure DEST_PATH_IMAGE028
表示的是空间内的坐标值；
Figure DEST_PATH_IMAGE030
表示的是自然常数；
Figure DEST_PATH_IMAGE032
表示的是空间中一点的温度对时间的变化率；
Figure DEST_PATH_IMAGE034
表示的是拉普拉斯算子。
4.根据权利要求2所述的一种肝脏分段方法，其特征在于，所述将所述肝脏任意同一位置作为所述最大传递温度对应的所述门静脉的供血区域的公式为：
Figure DEST_PATH_IMAGE036
；
Figure 331570DEST_PATH_IMAGE016
表示的是所述肝脏任意位置的坐标值；
Figure 920814DEST_PATH_IMAGE014
表示的是以所述各门静脉为热源，所述肝脏任意位置在任意时刻下的所述传递温度的温度值；
Figure 887502DEST_PATH_IMAGE012
表示的是所述各门静脉标号；
Figure 954815DEST_PATH_IMAGE018
表示的是任意时刻；
Figure DEST_PATH_IMAGE038
表示的是所述肝脏任意位置所属的所述门静脉的供血区域。
5.根据权利要求3所述的一种肝脏分段方法，其特征在于，根据所述各门静脉的半径设置初始温度的公式为：
Figure DEST_PATH_IMAGE040
；
Figure DEST_PATH_IMAGE042
表示的是所述肝脏任意位置的初始温度；
Figure DEST_PATH_IMAGE044
表示的是在
Figure 80772DEST_PATH_IMAGE016
处所述门静脉的半径数值，将该数值作为初始温度。
6.根据权利要求1所述的一种肝脏分段方法，其特征在于，还包括：
获取CT数据；
根据所述CT数据以及对所述肝脏进行分段的分段结果进行三维重建，得到经过分段后的三维模型。
7.一种肝脏分段装置，其特征在于，包括：
获取模块，用于获取各门静脉的半径以及所述各门静脉在肝脏中的位置；
计算模块，用于根据所述各门静脉的半径以及所述各门静脉在所述肝脏中的位置确定所述各门静脉的供血区域，具体包括：
根据所述各门静脉的半径设置初始温度；
根据所述各门静脉在所述肝脏中的位置，分别以所述各门静脉作为热源对应计算所述肝脏各个位置在任意时刻的传递温度；
根据所述肝脏各个位置在任意时刻的所述传递温度确定所述各门静脉的供血区域；
分段模块，用于根据所述供血区域对所述肝脏进行分段。
8.一种电子设备，其特征在于，包括处理器以及存储器，所述存储器存储有计算机可读取指令，当所述计算机可读取指令由所述处理器执行时，运行如权利要求1-6任一项所述方法中的步骤。
9.一种存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，运行如权利要求1-6任一项所述方法中的步骤。
CN202210979778.6A 2022-08-16 2022-08-16 一种肝脏分段方法、装置、电子设备及存储介质 Active CN115082649B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210979778.6A CN115082649B (zh) 2022-08-16 2022-08-16 一种肝脏分段方法、装置、电子设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210979778.6A CN115082649B (zh) 2022-08-16 2022-08-16 一种肝脏分段方法、装置、电子设备及存储介质

Publications (2)

Publication Number Publication Date
CN115082649A CN115082649A (zh) 2022-09-20
CN115082649B true CN115082649B (zh) 2022-12-02

Family

ID=83244011

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210979778.6A Active CN115082649B (zh) 2022-08-16 2022-08-16 一种肝脏分段方法、装置、电子设备及存储介质

Country Status (1)

Country Link
CN (1) CN115082649B (zh)

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8600719B2 (en) * 2010-02-09 2013-12-03 Fraunhofer Gesellschaft zur Förderung der angewandten Forschung e.V. Ablated object region determining apparatuses and methods
US9846765B2 (en) * 2012-11-09 2017-12-19 Siemens Healthcare Gmbh System and method for patient specific modeling of liver tumor ablation
CN103268603B (zh) * 2013-05-08 2015-10-28 同济大学 基于层级血管树划分的肝脏图像分段方法
CN103810752B (zh) * 2014-02-18 2017-01-11 海信集团有限公司 基于医学图像的肝脏分段方法及其肝脏分段系统
JP6491391B1 (ja) * 2018-11-22 2019-03-27 国立大学法人信州大学 血流動態解析システム、血流動態解析方法、およびプログラム
CN109712695B (zh) * 2018-12-25 2022-04-26 上海联影医疗科技股份有限公司 血管图像的处理方法及交互显示方法、计算机装置
CN114445424B (zh) * 2022-01-17 2022-11-18 推想医疗科技股份有限公司 模型训练方法和装置、肝段分割方法和装置

Also Published As

Publication number Publication date
CN115082649A (zh) 2022-09-20

Similar Documents

Publication Publication Date Title
JP4388958B2 (ja) 内視鏡的パスプランニングのための方法及びシステム
US7840044B2 (en) Method and system for liver lobe segmentation and pre-operative surgical planning
Seitel et al. Computer‐assisted trajectory planning for percutaneous needle insertions
US11744643B2 (en) Systems and methods facilitating pre-operative prediction of post-operative tissue function
US10176630B2 (en) Updating an electroanatomical map
US20190371474A1 (en) Methods and systems for modeling a necrotized tissue volume in an ablation procedure
US8150138B2 (en) Medical instrument and medical instrument control method
US11461895B2 (en) Automatic identification and processing of anatomical structures in an anatomical map
Ferrari et al. Value of multidetector computed tomography image segmentation for preoperative planning in general surgery
RU2635177C2 (ru) Устройство, способ и программа маркировки изображения объекта
CN102693540A (zh) 一种肝脏分段的方法及其系统
CN107818564B (zh) 一种肝脏3d医学影像分段方法
CA3138208A1 (en) Method for planning tissue ablation based on deep learning
WO2018165478A1 (en) Shell-constrained localization of vasculature
CN114332013B (zh) 基于肺动脉树分级的ct图像靶肺段识别方法
US12051156B2 (en) System and method for linking a segmentation graph to volumetric data
Haerle et al. Localized intraoperative virtual endoscopy (LIVE) for surgical guidance in 16 skull base patients
Nakao et al. Novel three-dimensional image simulation for lung segmentectomy developed with surgeons’ perspective
CN113876420B (zh) 用于规划手术路径的路径规划方法、系统、设备和介质
CN115082649B (zh) 一种肝脏分段方法、装置、电子设备及存储介质
CN112116694B (zh) 一种虚拟支气管镜辅助系统中三维模型绘制的方法及装置
US12016632B2 (en) System and method for planning surgical resection of lesions by a linear cutting stapler
Yoshida et al. Renal artery-based kidney segmentation on CT for patients with renal cell carcinoma: Feasibility of segmental artery clamping simulation
WO2021029973A1 (en) Methods and systems for modeling a necrotized tissue volume in an ablation procedure
US20220175455A1 (en) System and method for providing surgical planning and guidance with three-dimensional visualization

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
GR01 Patent grant
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 528253 Room 503, Floor 5, Building A, Jingu Zhichuang Industrial Community, No. 2, Yong'an North Road, Dawu Community, Guicheng Street, Nanhai District, Foshan City, Guangdong Province (residence declaration)

Patentee after: Guangdong Oupu Mandi Technology Co.,Ltd.

Address before: 528251 room 503, floor 5, building a, Jingu Zhichuang industrial community, No. 2, Yong'an North Road, Dawei community, Guicheng Street, Nanhai District, Foshan City, Guangdong Province

Patentee before: GUANGDONG OPTOMEDIC TECHNOLOGY CO.,LTD.