CN110246586A - 一种肝肿瘤微波消融术前仿真的有限元模型建立方法 - Google Patents

一种肝肿瘤微波消融术前仿真的有限元模型建立方法 Download PDF

Info

Publication number
CN110246586A
CN110246586A CN201910450725.3A CN201910450725A CN110246586A CN 110246586 A CN110246586 A CN 110246586A CN 201910450725 A CN201910450725 A CN 201910450725A CN 110246586 A CN110246586 A CN 110246586A
Authority
CN
China
Prior art keywords
liver
tumour
vascular pattern
finite element
emulation
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.)
Pending
Application number
CN201910450725.3A
Other languages
English (en)
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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201910450725.3A priority Critical patent/CN110246586A/zh
Publication of CN110246586A publication Critical patent/CN110246586A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • G06T2207/10012Stereo images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30056Liver; Hepatic

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Geometry (AREA)
  • Medical Informatics (AREA)
  • Data Mining & Analysis (AREA)
  • Epidemiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Pathology (AREA)
  • Computer Graphics (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明公开了一种肝肿瘤微波消融术前仿真的有限元模型建立方法,属于医学图像建模领域。该方法包括以下步骤:步骤1,消融术前对病人腹部进行CT扫描,得到肝脏周围断层图像;步骤2,在医学图像三维可视化模块中,对每层断层图像使用阈值分割、边界提取等方法结合手动勾画,分割出肝脏,肿瘤以及血管模型;步骤3,在逆向建模模块中将分割出的肝脏,肿瘤以及血管模型表面进行光滑;步骤4,在多物理场仿真模块中将平滑后的肝脏,肿瘤以及血管模型划分网格,建立有限元模型。本发明可以解决对于临近血管及危及器官类肿瘤,难以制定合适的消融热剂量以及调节微波消融针进针方式的问题。

Description

一种肝肿瘤微波消融术前仿真的有限元模型建立方法
技术领域
本发明涉及一种肝肿瘤微波消融术前仿真的有限元模型建立方法,属于医学图像建模领域。
背景技术
微波消融作为一种针对肝癌的有效治疗方式已进入NCCN(美国国立综合癌症网络)治疗指南中,然而对于其手术规划仍存在部分待解决的难题。现有的肝肿瘤微波消融术前规划在热剂量设定上是基于理想的二维轴对称模型计算得出的仿真温度场,而病人病情的差异使得医生需要对热剂量设定以及消融针穿刺位置做出调整。
肝脏临近多种重要器官且内部血管结构丰富,在临近大血管以及危及器官的肝肿瘤消融术中,缺少更加精准的术前规划,以明确治疗效果。
常用的对断层图像实现三维可视化的医学图像软件多是实现观测功能,而缺乏计算功能,生物组织在经过医学图像三维重建后的模型是极其不规则的且网格顺序混乱,用于有限元计算则需要对模型更进一步的修复。
发明内容
本发明提出了一种肝肿瘤微波消融术前仿真的有限元模型建立方法,利用病人术前CT断层图像数据,建立三维有限元模型,用于微波消融术前手术方案规划。
本发明为解决其技术问题采用如下技术方案:
一种肝肿瘤微波消融术前仿真的有限元模型建立方法,包括以下步骤:
步骤1,消融术前对病人腹部进行CT扫描,得到肝脏周围断层图像;
步骤2,在医学图像三维可视化模块中,对每层断层图像使用阈值分割、边界提取方法结合手动勾画,分割出肝脏,肿瘤以及血管模型;
步骤3,在逆向建模模块中将分割出的肝脏,肿瘤以及血管模型表面进行光滑处理;
步骤4,在多物理场仿真模块中将平滑后的肝脏,肿瘤以及血管模型划分网格,建立有限元模型。
步骤2的具体过程如下:
1) 导入病人断层图像之后,在医学图像三维可视化模块中建立三个分割对象,分割对象内容分别为肝脏,肿瘤以及血管;
2) 分割肝脏模型时,首先采用阈值分割的方法,选择合适的阈值范围,将肝脏组织CT值包含在内;
3) 阈值分割后的肝脏模型中包含其他组织,先将大块多余组织剪除后再逐层细化擦除;
4) 分割肿瘤以及血管模型时,首先采取边界识别的分割方法再逐层细化擦除;
5) 对上述细化擦除过的肝脏,肿瘤以及血管模型进行滤波处理。
步骤3的具体过程如下:
1) 对肝脏,肿瘤以及血管模型分别进行处理,初步进行表面光滑处理,减少非流行边,自相交,钉状物特征;
2) 对减少特征后的肝脏,肿瘤以及血管模型进行网格简化,减少网格数并且松弛模型以减少单独多边形之间的角度;
3) 局部细化平滑肝脏,肿瘤以及血管模型。
步骤4中的具体过程如下:
1) 光滑后的肝脏,肿瘤以及血管模型导入后生成实体,分别对其简化处理并修复缺陷;
2) 对导入的肝脏,肿瘤以及血管模型划分网格,肝脏,肿瘤以及血管模型使用不同的网格精度,生成有限元模型。
所述肝脏,肿瘤以及血管模型使用的网格精度如下:
对肝脏划分四边形网格单元最大尺寸为8.5 mm,最小尺寸1.53 mm,对肿瘤以及血管模型划分四边形网格单元最大尺寸为6.8 mm,最小尺寸0.85 mm。
本发明的有益效果如下:
本发明能解决对不同病人构建个性精准化三维微波消融术前仿真有限元模型的问题,将病人实际肿瘤位置、大小,临近器官及大血管的相对位置信息纳入有限元模型中,在后期的仿真中,可以综合调节微波消融针进针位置及消融热剂量,以达到最优治疗效果。
附图说明
图1是本发明有限元模型建立流程图。
图2是本发明使用医学图像三维可视化软件3D Slicer对影像数据阈值分割重建的组织图。
图3(a) 是本发明使用医学图像三维可视化软件3D Slicer建立的肝脏模型图;图3(b) 是本发明使用医学图像三维可视化软件3D Slicer建立的肿瘤模型图;图3(c) 是本发明使用医学图像三维可视化软件3D Slicer建立的血管模型图。
图4是本发明使用逆向建模软件Geomagic Studio对导入的肿瘤模型检测的非流行边,自相交,钉状物等特征图。
图5是本发明使用逆向建模软件Geomagic Studio对导入的肿瘤模型消除特征后的图。
图6(a) 是本发明使用逆向建模软件Geomagic Studio对初步修复后具有较为复杂的网格单元的肿瘤模型图;图6(b) 是本发明使用逆向建模软件Geomagic Studio对导入的肿瘤模型简化网格单元图。
图7(a) 是本发明使用逆向建模软件Geomagic Studio对导入的肿瘤模型突出表面孔洞化处理图;图7(b) 是本发明使用逆向建模软件Geomagic Studio对导入的肿瘤模型表面填补孔洞图。
图8(a) 是本发明使用逆向建模软件Geomagic Studio对导入的肝脏平滑后的模型图;图8(b) 是本发明使用逆向建模软件Geomagic Studio对导入的肿瘤平滑后的模型图;图8(c) 是本发明使用逆向建模软件Geomagic Studio对导入的血管修剪远离肿瘤处,保留临近血管并平滑后的模型图。
图9是本发明使用多物理场仿真软件Comsol Mutiphysics对平滑后的模型划分不同精度网格,建立的有限元模型图。
具体实施方式
下面结合附图对本发明进一步说明。
实施例:以腹部CT,临近血管处肝肿瘤模型为例。
图1表示本发明流程,一种肝肿瘤微波消融术前仿真的有限元模型建立方法,包括以下步骤:
(一)消融术前对病人腹部进行CT扫描,得到肝脏周围断层图像;
(二)在医学图像三维可视化模块中,对每层断层图像使用阈值分割、边界提取等方法结合手动勾画,分割出肝脏,肿瘤以及血管模型;
1) 导入病人断层图像之后,在医学图像三维可视化软件3D Slicer中建立三个分割对象,分割对象内容分别为肝脏,肿瘤以及血管;
2) 分割肝脏模型时,首先采用阈值分割的方法,选择合适的阈值范围,将肝脏组织CT值包含在内。本例设定-18 - 120的阈值,阈值分割结果如图2 所示;
3) 阈值分割后的肝脏模型包含其他组织,将大块多余组织剪除后再逐层细化擦除;
4) 分割肿瘤以及血管模型时,首先采取边界识别的分割方法再逐层细化擦除;
5) 对上述细化擦出过的肝脏,肿瘤以及血管模型进行滤波处理。对肝脏进行内核尺寸为3 mm的开运算和闭运算,对肿瘤和血管进行内核尺寸为1 mm的开运算和闭运算。图3(a),(b),(c) 所示分别为使用3D Slicer医学图像软件建立的肝脏、肿瘤和血管模型图;;
(三)在逆向建模模块中将分割出的肝脏,肿瘤以及血管模型表面进行光滑;
1) 在逆向建模软件Geomagic Studio中对肝脏,肿瘤以及血管模型分别处理,初步进行表面光滑,减少非流行边,自相交,钉状物等特征。图4所示为肿瘤包含的非流行边、自相交等多项特征的模型图,图5所示是减少此类特征后的模型图;
2) 对减少特征后的肝脏,肿瘤以及血管模型进行网格简化,减少网格数并且松弛模型以减少单独多边形之间的角度。图6(a),(b) 所示的分别是肿瘤网格简化前以及网格简化后的模型图;
3) 局部细化平滑肝脏,肿瘤以及血管模型。对表面极不平整处孔洞化处理后再用曲率填充孔洞。图7(a) 所示的是肿瘤局部孔洞化处理的模型图,图7(b) 所示的是肿瘤局部孔洞化处理后曲率填充孔洞的模型图;
4) 删除血管模型中距离肿瘤较远的部分以减少后期的计算量;图8(a),(b),(c) 所示的分别是平滑处理后的肝脏、肿瘤以及血管模型。
(四)在多物理场仿真模块中将平滑后的肝脏,肿瘤以及血管模型划分网格,建立有限元模型;
1) 光滑后的肝脏,肿瘤以及血管模型导入多物理场仿真软件Comsol Multiphysics后,分别对其简化处理并修复缺陷;
2) 对导入的模型划分网格,肝脏,肿瘤以及血管模型使用不同的网格精度,生成有限元模型。本例对肝脏划分四边形网格单元最大尺寸为8.5 mm,最小尺寸1.53 mm,对肿瘤以及血管模型划分四边形网格单元最大尺寸为6.8 mm,最小尺寸0.85 mm。生成的有限元模型如图9所示。

Claims (5)

1.一种肝肿瘤微波消融术前仿真的有限元模型建立方法,其特征在于,包括以下步骤:
步骤1,消融术前对病人腹部进行CT扫描,得到肝脏周围断层图像;
步骤2,在医学图像三维可视化模块中,对每层断层图像使用阈值分割、边界提取方法结合手动勾画,分割出肝脏,肿瘤以及血管模型;
步骤3,在逆向建模模块中将分割出的肝脏,肿瘤以及血管模型表面进行光滑处理;
步骤4,在多物理场仿真模块中将平滑后的肝脏,肿瘤以及血管模型划分网格,建立有限元模型。
2.根据权利要求1所述的一种肝肿瘤微波消融术前仿真的有限元模型建立方法,其特征在于:步骤2的具体过程如下:
1) 导入病人断层图像之后,在医学图像三维可视化模块中建立三个分割对象,分割对象内容分别为肝脏,肿瘤以及血管;
2) 分割肝脏模型时,首先采用阈值分割的方法,选择合适的阈值范围,将肝脏组织CT值包含在内;
3) 阈值分割后的肝脏模型中包含其他组织,先将大块多余组织剪除后再逐层细化擦除;
4) 分割肿瘤以及血管模型时,首先采取边界识别的分割方法再逐层细化擦除;
5) 对上述细化擦除过的肝脏,肿瘤以及血管模型进行滤波处理。
3.根据权利要求1所述的一种肝肿瘤微波消融术前仿真的有限元模型建立方法,其特征在于:步骤3的具体过程如下:
1) 对肝脏,肿瘤以及血管模型分别进行处理,初步进行表面光滑处理,减少非流行边,自相交,钉状物特征;
2) 对减少特征后的肝脏,肿瘤以及血管模型进行网格简化,减少网格数并且松弛模型以减少单独多边形之间的角度;
3) 局部细化平滑肝脏,肿瘤以及血管模型。
4.根据权利要求1所述的一种肝肿瘤微波消融术前仿真的有限元模型建立方法,其特征在于:步骤4中的具体过程如下:
1) 光滑后的肝脏,肿瘤以及血管模型导入后生成实体,分别对其简化处理并修复缺陷;
2) 对导入的肝脏,肿瘤以及血管模型划分网格,肝脏,肿瘤以及血管模型使用不同的网格精度,生成有限元模型。
5.根据权利要求4所述的一种肝肿瘤微波消融术前仿真的有限元模型建立方法,其特征在于:所述肝脏,肿瘤以及血管模型使用的网格精度如下:
对肝脏划分四边形网格单元最大尺寸为8.5 mm,最小尺寸1.53 mm,对肿瘤以及血管模型划分四边形网格单元最大尺寸为6.8 mm,最小尺寸0.85 mm。
CN201910450725.3A 2019-05-28 2019-05-28 一种肝肿瘤微波消融术前仿真的有限元模型建立方法 Pending CN110246586A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910450725.3A CN110246586A (zh) 2019-05-28 2019-05-28 一种肝肿瘤微波消融术前仿真的有限元模型建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910450725.3A CN110246586A (zh) 2019-05-28 2019-05-28 一种肝肿瘤微波消融术前仿真的有限元模型建立方法

Publications (1)

Publication Number Publication Date
CN110246586A true CN110246586A (zh) 2019-09-17

Family

ID=67885295

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910450725.3A Pending CN110246586A (zh) 2019-05-28 2019-05-28 一种肝肿瘤微波消融术前仿真的有限元模型建立方法

Country Status (1)

Country Link
CN (1) CN110246586A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111402321A (zh) * 2020-03-22 2020-07-10 华北理工大学 一种固体颗粒在高温熔池中目标识别及体积估计的方法
CN111685870A (zh) * 2020-06-23 2020-09-22 南京诺源医疗器械有限公司 一种用于肝脏的微波消融碳化调控方法
CN113378879A (zh) * 2021-05-06 2021-09-10 上海美杰医疗科技有限公司 术后肿瘤评估方法、装置及计算机存储介质
CN113468827A (zh) * 2021-06-17 2021-10-01 江苏大学 一种研究化疗药物抑制肿瘤生长的等几何仿真系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102117378A (zh) * 2009-12-31 2011-07-06 苏州瑞派宁科技有限公司 一种基于三维多模影像的肝肿瘤综合手术规划模拟方法及其系统
US20140136174A1 (en) * 2012-11-09 2014-05-15 Institut National De Recherche En Informatique Et En Automatique System and Method for Patient Specific Modeling of Liver Tumor Ablation
CN105208957A (zh) * 2013-02-26 2015-12-30 西门子公司 用于射频消融疗法的交互式患者特定仿真的系统和方法
CN105912874A (zh) * 2016-04-29 2016-08-31 青岛大学附属医院 基于dicom医学图像构建的肝脏三维数据库系统
CN107049475A (zh) * 2017-04-19 2017-08-18 纪建松 肝癌局部消融方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102117378A (zh) * 2009-12-31 2011-07-06 苏州瑞派宁科技有限公司 一种基于三维多模影像的肝肿瘤综合手术规划模拟方法及其系统
US20140136174A1 (en) * 2012-11-09 2014-05-15 Institut National De Recherche En Informatique Et En Automatique System and Method for Patient Specific Modeling of Liver Tumor Ablation
CN105208957A (zh) * 2013-02-26 2015-12-30 西门子公司 用于射频消融疗法的交互式患者特定仿真的系统和方法
CN105912874A (zh) * 2016-04-29 2016-08-31 青岛大学附属医院 基于dicom医学图像构建的肝脏三维数据库系统
CN107049475A (zh) * 2017-04-19 2017-08-18 纪建松 肝癌局部消融方法及系统

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111402321A (zh) * 2020-03-22 2020-07-10 华北理工大学 一种固体颗粒在高温熔池中目标识别及体积估计的方法
CN111685870A (zh) * 2020-06-23 2020-09-22 南京诺源医疗器械有限公司 一种用于肝脏的微波消融碳化调控方法
CN113378879A (zh) * 2021-05-06 2021-09-10 上海美杰医疗科技有限公司 术后肿瘤评估方法、装置及计算机存储介质
CN113378879B (zh) * 2021-05-06 2023-06-30 上海美杰医疗科技有限公司 术后肿瘤评估方法、装置及计算机存储介质
CN113468827A (zh) * 2021-06-17 2021-10-01 江苏大学 一种研究化疗药物抑制肿瘤生长的等几何仿真系统

Similar Documents

Publication Publication Date Title
CN110246586A (zh) 一种肝肿瘤微波消融术前仿真的有限元模型建立方法
US10878639B2 (en) Interactive voxel manipulation in volumetric medical imaging for virtual motion, deformable tissue, and virtual radiological dissection
US10950338B2 (en) Method and apparatus for generating an artificial intelligence 3D dataset and performing interactive manipulation and rendering of the dataset
Fuessinger et al. Planning of skull reconstruction based on a statistical shape model combined with geometric morphometrics
CN104462650B (zh) 一种可实现内外结构的实体化心脏3d模型制作方法
US10748438B2 (en) System and method for interactive patient specific simulation of radiofrequency ablation therapy
Reed et al. Automatic segmentation of whole breast using atlas approach and deformable image registration
CN106355582B (zh) 一种基于形状先验的肝脏医学图像分割的方法
CN108038862B (zh) 一种交互式医学图像智能分割建模方法
Whitfield et al. Automated delineation of radiotherapy volumes: are we going in the right direction?
CN107067398A (zh) 用于三维医学模型中缺失血管的补全方法及装置
CN104091347A (zh) 一种基于3d打印技术的颅内肿瘤手术规划模拟方法
Gerber et al. A multiscale imaging and modelling dataset of the human inner ear
EP1593087A2 (en) A method and system for image processing and contour assessment
Zachow et al. 3D reconstruction of individual anatomy from medical image data: Segmentation and geometry processing
CN102693540A (zh) 一种肝脏分段的方法及其系统
CN102254317A (zh) 牙种植导航中牙弓曲面的自动提取方法
Gall et al. Computer-aided planning and reconstruction of cranial 3D implants
CN108597017A (zh) 一种基于测量参数的骨骼变形模板构建方法
CN106725846A (zh) 一种基于人体器官3d模型的手术仿真系统及方法
CN111724389B (zh) 髋关节ct图像分割方法、装置、存储介质和计算机设备
Lv et al. Automatic segmentation of temporal bone structures from clinical conventional CT using a CNN approach
Stojkovic et al. Reverse modelling of human long bones using T-splines-case of tibia
Galstyan et al. Applications of 3D printing in breast cancer management
Wust et al. Evaluation of segmentation algorithms for generation of patient models in radiofrequency hyperthermia

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190917