CN111415404A - 术中预设区域的定位方法、装置、存储介质及电子设备 - Google Patents
术中预设区域的定位方法、装置、存储介质及电子设备 Download PDFInfo
- Publication number
- CN111415404A CN111415404A CN202010180471.0A CN202010180471A CN111415404A CN 111415404 A CN111415404 A CN 111415404A CN 202010180471 A CN202010180471 A CN 202010180471A CN 111415404 A CN111415404 A CN 111415404A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- roi
- region
- mask
- medical
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000013507 mapping Methods 0.000 claims abstract description 23
- 238000013135 deep learning Methods 0.000 claims description 30
- 238000002591 computed tomography Methods 0.000 claims description 29
- 238000000605 extraction Methods 0.000 claims description 26
- 210000004872 soft tissue Anatomy 0.000 claims description 20
- 210000000988 bone and bone Anatomy 0.000 claims description 10
- 238000012952 Resampling Methods 0.000 claims description 8
- 238000003759 clinical diagnosis Methods 0.000 claims description 8
- 210000003491 skin Anatomy 0.000 claims description 7
- 238000010219 correlation analysis Methods 0.000 claims description 6
- 238000002513 implantation Methods 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000004807 localization Effects 0.000 claims 1
- 239000000126 substance Substances 0.000 claims 1
- 238000003384 imaging method Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 210000001519 tissue Anatomy 0.000 description 7
- 238000013528 artificial neural network Methods 0.000 description 6
- 230000003902 lesion Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000010882 preoperative diagnosis Methods 0.000 description 3
- 238000002583 angiography Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 239000007943 implant Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000002560 therapeutic procedure Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000002324 minimally invasive surgery Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/08—Projecting images onto non-planar surfaces, e.g. geodetic screens
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
-
- 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/13—Edge detection
-
- 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/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
-
- 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
- G06T7/73—Determining position or orientation of objects or cameras using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Multimedia (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开了一种术中预设区域的定位方法、装置、存储介质及电子设备,方法包括:获取预设区域的医学三维影像,对医学三维影像进行感兴趣区域ROI的提取;在手术中获取ROI的多张不同角度的二维投影平片;根据多张不同角度的二维投影平片、医学三维影像以及深度学习网络模型,生成ROI与每一二维投影平片对应的区域掩模;对ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射,生成术中ROI的三维定位信息,以进行预设区域的定位。本发明提供的术中预设区域的定位方法、装置、存储介质及电子设备,可以实现精确的三维定位,解决了现有技术中常规手术定位的剂量过高的问题,实现手术中精确、便捷、高效的定位需求。
Description
技术领域
本发明是涉及图像处理技术领域,特别是涉及一种术中预设区域的定位方法、装置、存储介质及电子设备。
背景技术
穿刺、介入和微创等精细化手术目前在临床检查和治疗上得到广泛应用。对于无明确边界、无法直接观测性质的内在脏器病灶(如结节,肿瘤等感兴趣区域),目前主要通过电子计算机X射线断层扫描技术(英文全称:computed tomography,英文缩写:CT)或磁共振图像的三维图像(英文全称:Magnetic Resonance,英文缩写:MR)来进行术前诊断。在诊断后的穿刺、介入和微创等精细化手术过程中,需要精确定位和精细治疗,因此需要精确的医学成像技术引导医生进行辅助定位,获取精准的三维空间位置,从而进行准确的手术治疗。
目前用于术中引导的医学成像技术,主要包括数字减影血管造影(英文全称:digital subtraction angiography,英文缩写:DSA)、超声、CT等。DSA应用于血管造影和介入治疗中,但对末端的血管和普通的软组织无法成像和精确定位;超声对无骨性组织反射的浅表性脏器结节定位效果较好,但对医生的成像手法和适应部位存在较大的限制。CT是目前应用较为广泛的术中医学成像技术,定位准确、安全性高,速度快。
基于此,本申请的发明人发现,虽然术中CT定位准确、安全性高,速度快,但是使用剂量高,需多次定位的手术场景会导致患者受照剂量的累加;设备使用场景受限,无法安装至手术室等无屏蔽的空间,对于穿梭于CT屏蔽室和无菌的手术环境有较大限制;设备价格昂贵,无法普及于常规的手术科室。无法满足手术中精确、便捷、高效的定位需求。
公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。
发明内容
为解决上述问题,本发明实施例的目的在于提供一种术中预设区域的定位方法、装置、存储介质及电子设备。
为实现上述目的,本发明提供了一种术中预设区域的定位方法,包括:获取预设区域的医学三维影像,并对所述医学三维影像进行感兴趣区域ROI的提取,所述医学三维影像为CT或MR影像,所述ROI为人体内具备临床诊断意义的特征区域或外部植入器械所在区域;在手术中获取所述ROI的多张不同角度的二维投影平片;根据所述多张不同角度的二维投影平片、所述医学三维影像以及深度学习网络模型,生成所述ROI与每一二维投影平片对应的区域掩模,其中,所述区域掩膜包括ROI在二维投影平片中的二维轮廓定位信息;对所述ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射,生成术中所述ROI的三维定位信息;根据所述ROI的三维定位信息进行预设区域的定位。
在一优选的实施方式中,在对所述医学三维影像进行感兴趣区域ROI的提取之后,还包括:获得所述ROI的三维空间掩膜。
在一优选的实施方式中,所述根据多张不同角度的二维投影平片、所述医学三维影像以及深度学习网络模型,生成所述ROI与每一二维投影平片对应的区域掩模包括:根据所述多张不同角度的二维投影平片生成术中的三维体位模型,所述三维体位模型包括皮肤、软组织、骨骼的分布信息和外轮廓信息;根据所述三维体位模型以及所述深度学习网络模型,对所述医学三维影像和三维空间掩膜进行形变,生成模拟的术中医学三维影像和模拟的术中三维空间掩膜;对所述模拟的术中医学三维影像和所述模拟的术中三维空间掩膜对应ROI进行二维空间映射,生成所述ROI与每一二维投影平片对应的区域掩模以及模拟的二维影像。
在一优选的实施方式中,所述根据所述多张不同角度的二维投影平片生成术中的三维体位模型包括:对所述多张不同角度的二维投影平片进行预处理;对预处理后的二维投影平片进行重采样,根据重采样后的二维投影平片生成术中的三维体位模型。
在一优选的实施方式中,所述生成所述ROI与每一二维投影平片对应的区域掩模之后,对ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射之前,还包括:根据所述区域掩膜以及与所述区域掩膜对应的二维投影平片,确定特定区域;通过软组织增强和边缘信息提取算法,对所述特定区域进行图像结构特征提取;根据特征提取后的数据,对所述模拟的二维影像进行相关性分析,确定精确的二维定位掩模。
本发明还提供了一种术中预设区域的定位装置,包括:提取模块,用于获取预设区域的医学三维影像,并对所述医学三维影像进行感兴趣区域ROI的提取,所述医学三维影像为CT或MR影像,所述ROI为人体内具备临床诊断意义的特征区域或外部植入器械所在区域;获取模块,用于在手术中获取所述ROI的多张不同角度的二维投影平片;深度学习模块,用于根据所述多张不同角度的二维投影平片、所述医学三维影像以及深度学习网络模型,生成所述ROI与每一二维投影平片对应的区域掩模,其中,所述区域掩膜包括ROI在二维投影平片中的二维轮廓定位信息;映射模块,用于对所述ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射,生成术中所述ROI的三维定位信息;定位模块,用于根据所述ROI的三维定位信息进行预设区域的定位。
在一优选的实施方式中,所述提取模块用于在对所述医学三维影像进行感兴趣区域ROI的提取之后,获得所述ROI的三维空间掩膜;所述深度学习模块还用于:根据所述多张不同角度的二维投影平片生成术中的三维体位模型,所述三维体位模型包括皮肤、软组织、骨骼的分布信息和外轮廓信息;根据所述三维体位模型以及所述深度学习网络模型,对所述医学三维影像和三维空间掩膜进行形变,生成模拟的术中医学三维影像和模拟的术中三维空间掩膜;对所述模拟的术中医学三维影像和所述模拟的术中三维空间掩膜对应ROI进行二维空间映射,生成所述ROI与每一二维投影平片对应的区域掩模以及模拟的二维影像。
在一优选的实施方式中,还包括精确模块,用于:根据所述区域掩膜以及与所述区域掩膜对应的二维投影平片,确定特定区域;通过软组织增强和边缘信息提取算法,对所述特定区域进行图像结构特征提取;根据特征提取后的数据,对所述模拟的二维影像进行相关性分析,确定精确的二维定位掩模。
本实施例提供的术中预设区域的定位方法、装置、存储介质及电子设备,通过DR采集有限角度数张二维平片,生成ROI与每一二维投影平片对应的区域掩模,对所述ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射,生成术中所述ROI的三维定位信息,解决了现有技术中常规手术CT定位的剂量过高、手术环境有较大限制、无法普及于常规的手术科室的问题,实现手术中精确、便捷、高效的定位需求。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出了本发明一实施例所提供的术中预设区域的定位方法的流程图;
图2示出了本发明另一实施例所提供的术中预设区域的定位方法的流程图;
图3示出了本发明一实施例所提供的术中预设区域的定位装置的结构示意图;
图4示出了本发明实施例所提供的用于执行术中预设区域的定位方法的电子设备的结构示意图。
具体实施方式
下面结合附图,对本发明的具体实施方式进行详细描述,但应当理解本发明的保护范围并不受具体实施方式的限制。
数字化X线摄影(DR)是目前应用最为广泛的X射线成像设备,对于手术定位场景,一般采用移动DR或C型臂此类便携式DR应用骨科或者介入治疗中,其具有空间分辨高、成像速度快、价格低、辐射少等优势,但其仅能提供二维的组织堆叠图像,往往无法直接观测组织密度差异小的微小病变。另一方面,由于呼吸运动,摆位误差和姿势差异,术中DR摄影和术前诊断成像的体位往往存在较大差异,医学三维影像到二维投影的空间映射关系无法直接映射,因此,多角度的DR摄影无法直接对细微的软组织病灶进行精确定位。
本发明以深度学习技术为基础,拟在术中利用DR图像与医学三维影像信息相结合,进行特征提取和三维空间变换,对感兴趣区域进行更为精确的三维定位,为手术诊断和治疗提供精确、便捷、高效的三维定位方法,同时可显著降低患者和医护人员的受照剂量。
本发明实施例提供一种术中预设区域的定位方法,该方法用于实现手术中的病灶定位,参见图1所示,该方法包括:
步骤1,获取预设区域的医学三维影像,并对所述医学三维影像进行感兴趣区域ROI的提取,所述ROI为具备临床诊断意义的特征区域或外部植入器械所在区域,所述医学三维影像为CT或MR影像。
其中,预设区域可以为病灶区域,或者术前的诊断区域。预设区域包括一个或多个感兴趣区域。该医学三维影像包含正常人体组织和ROI的密度(水成分)差异和空间信息,可为临床医生提供诊断信息。所述的ROI包括肿块、结节、狭窄区域等可疑病灶和其他具备临床诊断意义的特征区域;感兴趣区域除了上述的临床诊断区域外,还可包括临床手术器械、植入支架等外部植入物所在区域。主要是应用在手术穿刺等场景中,例如达到患者身体哪里,由此判断病人的安全、是否触及病灶等。
在一种实现方式中,可以在术前获取病灶的医学三维影像,即术前的诊断CT或MR影像,以降低对术中设备的要求。可以通过图像处理技术实现所述医学三维影像进行感兴趣区域ROI的提取,所述的图像处理技术包括图像增强、特征识别、深度学习等通用分割技术。
该CT或MR医学三维影像包含正常人体组织和ROI的密度(水成分)差异和空间信息,可为临床医生提供诊断信息。
步骤2,在手术中获取所述ROI的多张不同角度的二维投影平片DR。
由于操作者在术中进行的拍摄时,X射线源到人体距离、人体到探测器距离,射线源和探测器的偏转角的改变均可能存在差异,因此需要进行多张不同角度的二维投影平片的拍摄。
步骤3,根据所述多张不同角度的二维投影平片(DR)、所述医学三维影像以及深度学习网络模型,生成所述ROI与每一二维投影平片对应的区域掩模,其中,所述区域掩膜包括ROI在二维投影平片的二维轮廓定位信息。
具体的,多张不同角度的二维投影平片可以为2张或更多角度进行拍摄的DR影片,一个病人可能有N个ROI、每张平片上每个ROI对应一个区域掩模。
在一种实现方式中,区域掩膜还可以包括ROI的形状信息,ROI的形状信息具有边缘信息。
步骤4,对所述ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射,生成术中所述ROI的三维定位信息;
步骤5,根据所述ROI的三维定位信息进行预设区域的定位。
本实施例提供的术中预设区域的定位方法,通过DR采集有限角度数张二维平片,生成ROI与每一二维投影平片对应的区域掩模,对所述ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射,生成术中所述ROI的三维定位信息,解决了现有技术中常规手术CT定位的剂量过高、手术环境有较大限制、无法普及于常规的手术科室的问题,实现手术中精确、便捷、高效的定位需求。
在一种实现方式中,在对所述医学三维影像进行感兴趣区域ROI的提取之后,还包括:获得所述ROI的三维空间掩膜。
在一种实施方式中,本实施例对步骤3进行了进一步限定,参见图2所示,具体包括:
步骤301,根据所述多张不同角度的二维投影平片生成术中的三维体位模型,所述三维体位模型包括皮肤、软组织、骨骼的分布信息和外轮廓信息;
具体的,可以首先对所述多张不同角度的二维投影平片进行预处理;对预处理后的二维投影平片进行重采样,然后根据重采样后的二维投影平片生成术中的三维体位模型。
其中,预处理,包括但不限于灰度归一化、标准化、边缘增强、对比度增强等。所述的重采样,可仅为DR二维图像的重采样,也可为DR和三维图像的重采样;重采样的插值方法包括但不限于双线性插值、样条插值等。例如DR的图像分辨率为M*N,CT的分辨率为O*O*P,可将DR重新采样为O*P,也可将两者图像重采样成X*X*Y。以满足不同分辨率图像输入至神经网络的要求。
在一个实施例中,术中三维体位模型的生成,可通过生成对抗网络(GAN)进行生成。
首先将同一体位下的多张DR影像作为神经网络的输入,通过该体位的三维图像分割获取基本组织的三维体位模型(皮肤,软组织,骨骼等三维容积下的定位信息)作为标签。模型的损失函数可以为MSE、DICE等相似性指标,对GAN网络进行预训练,其中,之后,将多张DR图像输入至训练的GAN网络,通过网络L=GAN(I1,…,In)即可生成术中定位的三维体位模型。
在另一个实施例中,也可通过术中体位下的DR摄影条件建立成像系统矩阵,进行二维三维空间信息的转换。通过三维空间关系和投影原理,利用三维图像重建方法,将多张DR进行反向投影至三维空间目标中,利用图像保真项进行重建图像的约束,并以图像正则项修复稀疏角度造成的噪声问题,同时通过术前医学三维影像获取的三维体位模型作为图像先验项,最终重建出术中的三维体位模型。
步骤302,根据所述三维体位模型以及所述深度学习网络模型,对所述医学三维影像和三维空间掩膜进行形变,生成模拟的术中医学三维影像和模拟的术中三维空间掩膜。
具体的,可通过深度学习网络进行术中医学三维影像的生成。将术前医学三维影像进行利用通用预分割(如阈值分割),生成包含组织分布(如软组织、骨骼、空腔等分布)的术前三维体位模型,并和术中三维体位模型,一同输入至变形神经网络中。利用三维体位模型轮廓特征的相似性作为神经网络的损失函数,三维空间形变场特征值(包括体素对应形变场的方向和形变距离),通过三维空间变换层对术前三维体位模型进行三维形变,并以损失函数计算术前和术中的三维体位模型的轮廓相似性,最后反向传播至形变特征值进行迭代,最终获得满足loss的三维空间形变场。通过该三维空间形变场,对术前医学三维影像和三维空间掩模进行三维形变,最终获得模拟的术中医学三维影像和三维空间掩模。
步骤303,对所述模拟的术中医学三维影像和所述模拟的术中三维空间掩膜对应ROI进行二维空间映射,生成所述ROI与每一二维投影平片对应的区域掩模以及模拟的二维影像。
具体的,根据术中DR的探测器参数、射线源参数、三维空间位置、病人体位等信息,建立投影系统,实现二维空间映射,最终将三维空间掩模投影至二维空间中,以获得所述ROI与每一二维投影平片对应的区域掩模以及模拟的二维影像。
在步骤33之后,步骤4之前,还包括确定精确的二维定位掩模的步骤,在一种实现方式中,确定精确的二维定位掩模的步骤包括:
步骤304,根据所述区域掩膜以及与所述区域掩膜对应的二维投影平片,确定特定区域;
其中,特定区域为根据区域掩膜以及与该区域掩膜对应的二维投影平片确定的病灶所在位置。
步骤305,通过软组织增强和边缘信息提取算法,对所述特定区域进行图像结构特征提取;
对该特定区域的图像进行如软组织增强,并利用sobel或拉普拉斯算子等边缘强化处理方式对特定区域进行边缘结构特征提取。
步骤306,根据所述特征提取后的数据,对所述模拟的二维影像进行相关性分析,确定精确的二维定位掩模。
通过特征提取后的数据与二维影像进行相关性作为目标函数,进行ROI区域搜索和更新,最终定位出精确的术中二维定位掩模。
在另一种实现方式中,确定精确的二维定位掩模的步骤包括:可通过人工对DR或CT的数据进行人工或其他通用方法标注骨骼获得对应骨骼标签,然后对深度学习网络进行预先训练,实现DR平片骨骼分离;将上述特定区域的DR图像输入至深度学习网络,进行骨骼抑制、软组织增强等方法,实现ROI纹理信息的增强和边缘提取,最终获得定位掩模。
在又一种实现方式中,确定精确的二维定位掩模的步骤包括:可以通过传统的图像后处理增强算法,将骨性组织进行分割和抑制后,利用边缘提取算法进行软组织高频分量增强。具体的,可以通过以下公式获得:
En(x,y)=∑i,j∈MI(x,y)*h(i,j);
其中,En(x,y)为图像后处理增强算法之后的图像x行y列的像素的数值,I(x,y)原始图像对应x行y列的像素的数值,i,j分别为滤波算子的长及宽,h(i,j)为边缘增强滤波算子;再输入至预先训练神经网络(该网络使用经过处理的软组织图像和人工标记的ROI位置进行训练),进而获取ROI的定位掩模。再通过深度学习网络进行提取和分割。
在一种实现方式中,步骤4可以通过以下集中方式实现,所述ROI在多张二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射
在一个实施例中,可将多张二维投影平片的区域掩模和多张二维投影平片输入至生成对抗网络(GAN),直接生成ROI的三维定位信息。
在又一个实施例中,还可以通过GAN和成像该系统矩阵两个方法结合,利用GAN进行初步定位信息提取,再利用图像先验项和成像矩阵投影,进一步迭代优化,最终生成精确的三维定位信息。
由此,本实施例提供的定位方法,能通过数张DR平片实现精确的ROI三维定位,保证定位精度情况下,解决常规手术CT定位的剂量过高的问题。能够通过深度学习网络生成三维术中体位,直接生成模拟术中医学三维影像,无需多次CT照射即可实现ROI的定位。使用深度学习神经网络,融合三维二维的图像特征利用术前CT分割结果,直接定位术中DR无法直接观测的ROI信息,为精确术中定位提供新的技术思路。
参见图3所示,其为本实施例的术中预设区域的定位装置的结构示意图,包括:
提取模块1,用于获取预设区域的医学三维影像,并对所述医学三维影像进行感兴趣区域ROI的提取,所述医学三维影像为CT或MR影像,所述ROI为人体内具备临床诊断意义的特征区域或外部植入器械;
获取模块2,用于在手术中获取所述ROI的多张不同角度的二维投影平片;
深度学习模块3,用于根据所述多张不同角度的二维投影平片、所述医学三维影像以及深度学习网络模型,生成所述ROI与每一二维投影平片对应的区域掩模,其中,所述区域掩膜包括ROI在二维投影平片中的二维轮廓定位信息;
映射模块4,用于对所述ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射,生成术中所述ROI的三维定位信息;
定位模块5,用于根据所述ROI的三维定位信息进行预设区域的定位。
在一种实施方式中,所述提取模块1用于在对所述医学三维影像进行感兴趣区域ROI的提取之后,获得所述ROI的三维空间掩膜。
所述深度学习模块3还用于:根据所述多张不同角度的二维投影平片生成术中的三维体位模型,所述三维体位模型包括皮肤、软组织、骨骼的分布信息和外轮廓信息;根据所述三维体位模型以及所述深度学习网络模型,对所述医学三维影像和三维空间掩膜进行形变,生成模拟的术中医学三维影像和模拟的术中三维空间掩膜;对所述模拟的术中医学三维影像和所述模拟的术中三维空间掩膜对应ROI进行二维空间映射,生成所述ROI与每一二维投影平片对应的区域掩模以及模拟的二维影像。
深度学习模块3进一步用于,对所述多张不同角度的二维投影平片进行预处理;对预处理后的二维投影平片进行重采样,根据重采样后的二维投影平片生成术中的三维体位模型。
本实施例的定位装置还包括精确模块6,精确模块6用于:根据所述区域掩膜以及与所述区域掩膜对应的二维投影平片,确定特定区域;通过软组织增强和边缘信息提取算法,对所述特定区域进行图像结构特征提取;根据特征提取后的数据,对所述模拟的二维影像进行相关性分析,确定精确的二维定位掩模。
本实施例提供的术中预设区域的定位装置,通过DR采集有限角度数张二维平片,并结合深度学习技术,提取术中体位形变信息,结合术前三维图像生成术中DR ROI的基本二维定位信息,最终生成精确的三维定位信息。解决了现有技术中常规手术CT定位的剂量过高、手术环境有较大限制、无法普及于常规的手术科室的问题,实现手术中精确、便捷、高效的定位需求。
本发明实施例还提供了一种存储介质,所述存储介质存储有计算机可执行指令,其包含用于执行上述术中预设区域的定位方法的程序,该计算机可执行指令可执行上述任意方法实施例中的方法。
其中,所述存储介质可以是计算机能够存取的任何可用介质或数据存储设备,包括但不限于磁性存储器(例如软盘、硬盘、磁带、磁光盘(MO)等)、光学存储器(例如CD、DVD、BD、HVD等)、以及半导体存储器(例如ROM、EPROM、EEPROM、非易失性存储器(NAND FLASH)、固态硬盘(SSD))等。
图4示出了本发明的另一个实施例的一种电子设备的结构框图。所述电子设备1100可以是具备计算能力的主机服务器、个人计算机PC、或者可携带的便携式计算机或终端等。本发明具体实施例并不对电子设备的具体实现做限定。
该电子设备1100包括至少一个处理器(processor)1110、通信接口(Communications Interface)1120、存储器(memory array)1130和总线1140。其中,处理器1110、通信接口1120、以及存储器1130通过总线1140完成相互间的通信。
通信接口1120用于与上位机通信,下载数据等。
处理器1110用于执行程序。处理器1110可能是一个中央处理器CPU,或者是专用集成电路ASIC(Application Specific Integrated Circuit),或者是被配置成实施本发明实施例的一个或多个集成电路。
存储器1130用于可执行的指令。存储器1130可能包含高速RAM存储器,也可能还包括非易失性存储器(non-volatile memory),例如至少一个磁盘存储器。存储器1130也可以是存储器阵列。存储器1130还可能被分块,并且所述块可按一定的规则组合成虚拟卷。存储器1130存储的指令可被处理器1110执行,以使处理器1110能够执行上述任意方法实施例中的术中预设区域的定位方法。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种术中预设区域的定位方法,其特征在于,包括:
获取预设区域的医学三维影像,并对所述医学三维影像进行感兴趣区域ROI的提取,所述医学三维影像为CT或MR影像,所述ROI为人体内具备临床诊断意义的特征区域或外部植入器械所在区域;
在手术中获取所述ROI的多张不同角度的二维投影平片;
根据所述多张不同角度的二维投影平片、所述医学三维影像以及深度学习网络模型,生成所述ROI与每一二维投影平片对应的区域掩模,其中,所述区域掩膜包括ROI在二维投影平片中的二维轮廓定位信息;
对所述ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射,生成术中所述ROI的三维定位信息;
根据所述ROI的三维定位信息进行预设区域的定位。
2.如权利要求1所述的定位方法,其特征在于,在对所述医学三维影像进行感兴趣区域ROI的提取之后,还包括:获得所述ROI的三维空间掩膜。
3.如权利要求2所述的定位方法,其特征在于,所述根据多张不同角度的二维投影平片、所述医学三维影像以及深度学习网络模型,生成所述ROI与每一二维投影平片对应的区域掩模包括:
根据所述多张不同角度的二维投影平片生成术中的三维体位模型,所述三维体位模型包括皮肤、软组织、骨骼的分布信息和外轮廓信息;
根据所述三维体位模型以及所述深度学习网络模型,对所述医学三维影像和三维空间掩膜进行形变,生成模拟的术中医学三维影像和模拟的术中三维空间掩膜;
对所述模拟的术中医学三维影像和所述模拟的术中三维空间掩膜对应ROI进行二维空间映射,生成所述ROI与每一二维投影平片对应的区域掩模以及模拟的二维影像。
4.如权利要求3所述的定位方法,其特征在于,所述根据所述多张不同角度的二维投影平片生成术中的三维体位模型包括:
对所述多张不同角度的二维投影平片进行预处理;
对预处理后的二维投影平片进行重采样,根据重采样后的二维投影平片生成术中的三维体位模型。
5.如权利要求3所述的定位方法,其特征在于,所述生成所述ROI与每一二维投影平片对应的区域掩模之后,对ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射之前,还包括:
根据所述区域掩膜以及与所述区域掩膜对应的二维投影平片,确定特定区域;
通过软组织增强和边缘信息提取算法,对所述特定区域进行图像结构特征提取;
根据特征提取后的数据,对所述模拟的二维影像进行相关性分析,确定精确的二维定位掩模。
6.一种术中预设区域的定位装置,其特征在于,包括:
提取模块,用于获取预设区域的医学三维影像,并对所述医学三维影像进行感兴趣区域ROI的提取,所述医学三维影像为CT或MR影像,所述ROI为人体内具备临床诊断意义的特征区域或外部植入器械所在区域;
获取模块,用于在手术中获取所述ROI的多张不同角度的二维投影平片;
深度学习模块,用于根据所述多张不同角度的二维投影平片、所述医学三维影像以及深度学习网络模型,生成所述ROI与每一二维投影平片对应的区域掩模,其中,所述区域掩膜包括ROI在二维投影平片中的二维轮廓定位信息;
映射模块,用于对所述ROI在多张不同角度的二维投影平片的区域掩模以及多张不同角度的二维投影平片进行三维空间映射,生成术中所述ROI的三维定位信息;
定位模块,用于根据所述ROI的三维定位信息进行预设区域的定位。
7.如权利要求6所述的定位装置,其特征在于,
所述提取模块还用于在对所述医学三维影像进行感兴趣区域ROI的提取之后,获得所述ROI的三维空间掩膜;
所述深度学习模块还用于根据所述多张不同角度的二维投影平片生成术中的三维体位模型,所述三维体位模型包括皮肤、软组织、骨骼的分布信息和外轮廓信息;根据所述三维体位模型以及所述深度学习网络模型,对所述医学三维影像和三维空间掩膜进行形变,生成模拟的术中医学三维影像和模拟的术中三维空间掩膜;对所述模拟的术中医学三维影像和所述模拟的术中三维空间掩膜对应ROI进行二维空间映射,生成所述ROI与每一二维投影平片对应的区域掩模以及模拟的二维影像。
8.如权利要求7所述的定位装置,其特征在于,还包括精确模块,用于根据所述区域掩膜以及与所述区域掩膜对应的二维投影平片,确定特定区域;通过软组织增强和边缘信息提取算法,对所述特定区域进行图像结构特征提取;根据特征提取后的数据,对所述模拟的二维影像进行相关性分析,确定精确的二维定位掩模。
9.一种存储介质,其特征在于,所述存储介质存储有计算机可执行指令,所述计算机可执行指令用于执行权利要求1-5任意一项所述的术中预设区域的定位方法。
10.一种电子设备,其特征在于,包括:
至少一个处理器;以及,
与所述至少一个处理器通信连接的存储器;其中,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行权利要求1-5任意一项所述的术中预设区域的定位方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010180471.0A CN111415404B (zh) | 2020-03-16 | 2020-03-16 | 术中预设区域的定位方法、装置、存储介质及电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010180471.0A CN111415404B (zh) | 2020-03-16 | 2020-03-16 | 术中预设区域的定位方法、装置、存储介质及电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111415404A true CN111415404A (zh) | 2020-07-14 |
CN111415404B CN111415404B (zh) | 2021-06-29 |
Family
ID=71494310
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010180471.0A Active CN111415404B (zh) | 2020-03-16 | 2020-03-16 | 术中预设区域的定位方法、装置、存储介质及电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111415404B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112201335A (zh) * | 2020-07-23 | 2021-01-08 | 中国人民解放军总医院 | 一种线阵超声内镜下识别腹腔内结构系统及其方法 |
CN113011418A (zh) * | 2021-02-09 | 2021-06-22 | 杭州海康慧影科技有限公司 | 确定图像中待处理区域的方法、装置、设备 |
CN114842004A (zh) * | 2022-07-04 | 2022-08-02 | 真健康(北京)医疗科技有限公司 | 基于神经网络模型的穿刺位置验证方法及设备 |
CN114886445A (zh) * | 2022-07-15 | 2022-08-12 | 康达洲际医疗器械有限公司 | 一种基于多叶光栅动态可调的双c臂三维成像方法与系统 |
CN114947907A (zh) * | 2022-08-01 | 2022-08-30 | 康达洲际医疗器械有限公司 | 一种基于动态多叶光栅的三维c臂成像方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080269600A1 (en) * | 2007-04-24 | 2008-10-30 | Medtronic, Inc. | Flexible Array For Use In Navigated Surgery |
CN102651145A (zh) * | 2012-04-06 | 2012-08-29 | 哈尔滨工业大学 | 股骨三维模型可视化方法 |
CN102999902A (zh) * | 2012-11-13 | 2013-03-27 | 上海交通大学医学院附属瑞金医院 | 基于ct配准结果的光学导航定位系统及其导航方法 |
CN108601630A (zh) * | 2015-12-22 | 2018-09-28 | 斯栢迈股份公司 | 在脊柱区域和邻近胸廓、骨盆或头部的区域的外科手术期间的术中图控导航装置 |
CN110706336A (zh) * | 2019-09-29 | 2020-01-17 | 上海昊骇信息科技有限公司 | 一种基于医疗影像数据的三维重建方法及系统 |
CN110840561A (zh) * | 2019-11-11 | 2020-02-28 | 郑健青 | 一种基于人工智能与图论算法的手术导航辅助系统 |
-
2020
- 2020-03-16 CN CN202010180471.0A patent/CN111415404B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080269600A1 (en) * | 2007-04-24 | 2008-10-30 | Medtronic, Inc. | Flexible Array For Use In Navigated Surgery |
CN102651145A (zh) * | 2012-04-06 | 2012-08-29 | 哈尔滨工业大学 | 股骨三维模型可视化方法 |
CN102999902A (zh) * | 2012-11-13 | 2013-03-27 | 上海交通大学医学院附属瑞金医院 | 基于ct配准结果的光学导航定位系统及其导航方法 |
CN108601630A (zh) * | 2015-12-22 | 2018-09-28 | 斯栢迈股份公司 | 在脊柱区域和邻近胸廓、骨盆或头部的区域的外科手术期间的术中图控导航装置 |
CN110706336A (zh) * | 2019-09-29 | 2020-01-17 | 上海昊骇信息科技有限公司 | 一种基于医疗影像数据的三维重建方法及系统 |
CN110840561A (zh) * | 2019-11-11 | 2020-02-28 | 郑健青 | 一种基于人工智能与图论算法的手术导航辅助系统 |
Non-Patent Citations (2)
Title |
---|
ZHE LUO ET AL.: "Intra-Operative 2-D Ultrasound and Dynamic 3-D Aortic Model Registration for Magnetic Navigation of Transcatheter Aortic Valve Implantation", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 * |
赵岩: "计算机辅助手术导航中的跟踪定位方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112201335A (zh) * | 2020-07-23 | 2021-01-08 | 中国人民解放军总医院 | 一种线阵超声内镜下识别腹腔内结构系统及其方法 |
CN112201335B (zh) * | 2020-07-23 | 2023-05-26 | 中国人民解放军总医院 | 一种线阵超声内镜下识别腹腔内结构系统及其方法 |
CN113011418A (zh) * | 2021-02-09 | 2021-06-22 | 杭州海康慧影科技有限公司 | 确定图像中待处理区域的方法、装置、设备 |
CN113011418B (zh) * | 2021-02-09 | 2024-02-23 | 杭州海康慧影科技有限公司 | 确定图像中待处理区域的方法、装置、设备 |
CN114842004A (zh) * | 2022-07-04 | 2022-08-02 | 真健康(北京)医疗科技有限公司 | 基于神经网络模型的穿刺位置验证方法及设备 |
CN114842004B (zh) * | 2022-07-04 | 2022-10-21 | 真健康(北京)医疗科技有限公司 | 基于神经网络模型的穿刺位置验证方法及设备 |
CN114886445A (zh) * | 2022-07-15 | 2022-08-12 | 康达洲际医疗器械有限公司 | 一种基于多叶光栅动态可调的双c臂三维成像方法与系统 |
CN114886445B (zh) * | 2022-07-15 | 2022-12-13 | 康达洲际医疗器械有限公司 | 一种基于多叶光栅动态可调的双c臂三维成像方法与系统 |
CN114947907A (zh) * | 2022-08-01 | 2022-08-30 | 康达洲际医疗器械有限公司 | 一种基于动态多叶光栅的三维c臂成像方法及系统 |
WO2024027124A1 (zh) * | 2022-08-01 | 2024-02-08 | 康达洲际医疗器械有限公司 | 一种基于动态多叶光栅的三维c臂成像方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111415404B (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111415404B (zh) | 术中预设区域的定位方法、装置、存储介质及电子设备 | |
Markelj et al. | A review of 3D/2D registration methods for image-guided interventions | |
US8897514B2 (en) | Imaging method for motion analysis | |
US8942455B2 (en) | 2D/3D image registration method | |
US8233692B2 (en) | Method of suppressing obscuring features in an image | |
JP2022517581A (ja) | 冠動脈の動的ロードマップを提供するための方法およびシステム | |
JP2022507622A (ja) | 拡張現実ディスプレイでの光学コードの使用 | |
US20070167784A1 (en) | Real-time Elastic Registration to Determine Temporal Evolution of Internal Tissues for Image-Guided Interventions | |
EP2056255B1 (en) | Method for reconstruction of a three-dimensional model of an osteo-articular structure | |
JP2010246883A (ja) | 患者位置決めシステム | |
JP2009536857A (ja) | 画像誘導放射線医療のための画像の変形可能な登録 | |
US10614598B2 (en) | Systems and methods for generating 2D projection from previously generated 3D dataset | |
JP2021520236A (ja) | 解剖学的構造に基づくフレームレス二次元/三次元画像レジストレーション | |
TWI836491B (zh) | 註冊二維影像資料組與感興趣部位的三維影像資料組的方法及導航系統 | |
Zollei | 2D-3D rigid-body registration of X-ray Fluoroscopy and CT images | |
Puggelli et al. | Accuracy assessment of CT-based 3D bone surface reconstruction | |
US20230316550A1 (en) | Image processing device, method, and program | |
De Silva et al. | Robust 2-D–3-D registration optimization for motion compensation during 3-D TRUS-guided biopsy using learned prostate motion data | |
Al Abboodi et al. | Supervised Transfer Learning for Multi Organs 3D Segmentation With Registration Tools for Metal Artifact Reduction in CT Images | |
CN114283179A (zh) | 基于超声图像的骨折远近端空间位姿实时获取与配准系统 | |
Al-Shayea et al. | An efficient approach to 3d image reconstruction | |
TWI842001B (zh) | 註冊二維影像資料組與感興部位的三維影像資料組的方法及導航系統 | |
Rizqie et al. | 3D coordinate reconstruction from 2D X-ray images for guided lung biopsy | |
WO2024002221A1 (zh) | 介入手术影像辅助方法、系统、装置及存储介质 | |
TWI773491B (zh) | 肺功能鑑別方法 |
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 |