CN111640106A - Multimode medical image conversion method based on artificial intelligence - Google Patents

Multimode medical image conversion method based on artificial intelligence Download PDF

Info

Publication number
CN111640106A
CN111640106A CN202010487136.5A CN202010487136A CN111640106A CN 111640106 A CN111640106 A CN 111640106A CN 202010487136 A CN202010487136 A CN 202010487136A CN 111640106 A CN111640106 A CN 111640106A
Authority
CN
China
Prior art keywords
image data
image
oct
ivus
fusion
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
CN202010487136.5A
Other languages
Chinese (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.)
Fuwai Shenzhen Hospital, CAMS&PUMC
Original Assignee
Fuwai Shenzhen Hospital, CAMS&PUMC
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 Fuwai Shenzhen Hospital, CAMS&PUMC filed Critical Fuwai Shenzhen Hospital, CAMS&PUMC
Priority to CN202010487136.5A priority Critical patent/CN111640106A/en
Publication of CN111640106A publication Critical patent/CN111640106A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • 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
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/194Segmentation; Edge detection involving foreground-background segmentation
    • 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/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • 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/10072Tomographic images
    • G06T2207/10101Optical tomography; Optical coherence tomography [OCT]
    • 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/10072Tomographic images
    • G06T2207/10104Positron emission tomography [PET]
    • 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/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20052Discrete cosine transform [DCT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20172Image enhancement details
    • G06T2207/20192Edge enhancement; Edge preservation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging
    • 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/30101Blood vessel; Artery; Vein; Vascular

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明公开了一种基于人工智能的多模态医学影像的转换方法,主要涉及医学影像领域。包括获取主动脉夹层的OCT图像数据;获取主动脉夹层的IVUS图像数据;分别对OCT图像数据和IVUS图像数据的背景采用全局阈值分割的方法,提取主动脉的外轮廓,并利用形态学孔洞填充和面积阈值方法将噪声信息进行清理;通过离散余弦变换的频域信息融合方法,将OCT图像数据和IVUS图像数据进行图像信息融合,获得融合图像。本发明的有益效果在于:它结合多种模态的胸主动脉夹层图像,将超声图像中的结构特征融合到OCT图像中,并在融合的多模态图像中有效的增强真腔、假腔和血栓的特征。

Figure 202010487136

The invention discloses a multi-modal medical image conversion method based on artificial intelligence, which mainly relates to the field of medical images. Including obtaining OCT image data of aortic dissection; obtaining IVUS image data of aortic dissection; using global threshold segmentation method for the background of OCT image data and IVUS image data respectively, extracting the outer contour of the aorta, and using morphological hole filling and area threshold method to clean up the noise information; through the frequency domain information fusion method of discrete cosine transform, the image information of OCT image data and IVUS image data is fused to obtain the fusion image. The beneficial effect of the present invention is that it combines the thoracic aortic dissection images of multiple modalities, fuses the structural features in the ultrasound image into the OCT image, and effectively enhances the true cavity and the false cavity in the fused multimodal image. and thrombus characteristics.

Figure 202010487136

Description

一种基于人工智能的多模态医学影像的转换方法A method for converting multimodal medical images based on artificial intelligence

技术领域technical field

本发明涉及医学影像领域,具体是一种基于人工智能的多模态医学影像的转换方法。The invention relates to the field of medical imaging, in particular to a multimodal medical imaging conversion method based on artificial intelligence.

背景技术Background technique

近年来,随着多种医学成像技术(CT.MRI.Utrasound.PET、OCT、Microscop等)的迅速发展,并已经被广泛应用于疾病的早期检测、诊断和治疗中,致使临床中,影像学专家每天需要面临大量的患者医学影像数据,单纯的依靠人工阅读和判定这些影像数据,不仅耗时费力,同时带有主观性。因此,计算机辅助诊断技术成为解决该问题的一个有力措施。尤其是,最近几年兴起的一种机器学习方法“深度学习”的引入,逐渐加快了现代计算机医学影像计算和分析迈入智能化和精确化的步伐。In recent years, with the rapid development of a variety of medical imaging technologies (CT, MRI, Utrasound, PET, OCT, Microscop, etc.), they have been widely used in the early detection, diagnosis and treatment of diseases. Experts need to face a large amount of patient medical image data every day. Simply relying on manual reading and judgment of these image data is not only time-consuming and labor-intensive, but also subjective. Therefore, computer aided diagnosis technology has become a powerful measure to solve this problem. In particular, the introduction of "deep learning", a machine learning method that has emerged in recent years, has gradually accelerated the pace of modern computer medical image computing and analysis toward intelligence and precision.

胸主动脉夹层是一个很险恶的、凶险的心血管疾病,它的致残和致死率很高。有统计数据表明,一旦发现急性胸部夹层动脉瘤,尤其是Ⅰ型的夹层动脉瘤,如果不予以积极的处理,它的一个月的死亡率能达到70%-80%。胸主动脉夹层是指在主动脉壁存在或不存在自身病变的基础上,并在一系列可能外因(如高血压、外伤等)的作用下导致主动脉内膜撕裂,血液由内膜撕裂口进入主动脉壁中层,造成主动脉中层沿长轴分离,管腔呈现真假两腔的一种病理状态。夹层好发部位在主动脉根部至左锁骨下动脉开口以远的弧形区域。该病起病突然,患者会出现剧烈疼痛、休克和压迫症状。少数患者由于心脏压塞、大量出血、恶性高血压、严重的主动脉瓣反流和心肌、中枢神经和肾脏的持续性缺血等,可致迅速死亡。故及时的诊断和治疗,对患者的生命和生活质量具有至关重要的作用。Thoracic aortic dissection is a very sinister and dangerous cardiovascular disease with high morbidity and mortality. Statistics show that once an acute thoracic dissecting aneurysm is found, especially a type I dissecting aneurysm, its one-month mortality rate can reach 70%-80% if it is not actively treated. Thoracic aortic dissection refers to the tear of the aortic intima under the action of a series of possible external factors (such as hypertension, trauma, etc.) on the basis of the presence or absence of its own pathological changes in the aortic wall. The rupture enters the middle layer of the aortic wall, causing the middle layer of the aorta to separate along the long axis, and the lumen presents a pathological state of true and false lumens. The prevalent site of dissection is the arcuate region from the aortic root to the opening of the left subclavian artery. The onset of the disease is sudden, with severe pain, shock, and compression. A small number of patients can die rapidly due to cardiac tamponade, massive hemorrhage, malignant hypertension, severe aortic regurgitation, and persistent ischemia of the myocardium, central nervous system, and kidneys. Therefore, timely diagnosis and treatment play a vital role in the patient's life and quality of life.

现在对于胸主动脉夹层的检查一般要辅助胸部的医学影像。例如CT、X线、核磁共振等。其中,OCT已经成为血管腔内成像分辨率最高的新兴成像方式,逐渐在胸主动脉夹层疾病的诊断中凸显其其重要地位。但是,考虑到胸主动脉夹层中会存在假腔、真腔和血栓多种复杂的解剖学结构,同时由于动脉壁三层组织结构的特点,在OCT成像中会存在反向散射,造成信号"的衰减和渗透,同时存在散斑噪声的干扰,最终形成弱边缘现象,这给OCT图像的自动化分割带来了研究的难点。At present, the examination of thoracic aortic dissection is generally supplemented by medical imaging of the chest. Such as CT, X-ray, nuclear magnetic resonance and so on. Among them, OCT has become an emerging imaging modality with the highest resolution in vascular intraluminal imaging, and has gradually highlighted its important position in the diagnosis of thoracic aortic dissection. However, considering that there are various complex anatomical structures of false lumen, true lumen and thrombus in thoracic aortic dissection, and due to the characteristics of the three-layer tissue structure of the arterial wall, there will be backscattering in OCT imaging, resulting in signal" At the same time, there is the interference of speckle noise, which eventually forms a weak edge phenomenon, which brings difficulties to the automatic segmentation of OCT images.

发明内容SUMMARY OF THE INVENTION

本发明的目的在于提供一种基于人工智能的多模态医学影像的转换方法,它结合多种模态的胸主动脉夹层图像,将超声图像中的结构特征融合到OCT图像中,并在融合的多模态图像中有效的增强真腔、假腔和血栓的特征。The purpose of the present invention is to provide a multi-modal medical image conversion method based on artificial intelligence, which combines multiple modalities of thoracic aortic dissection images, fuses the structural features in the ultrasound image into the OCT image, and fuses the structural features in the OCT image. Effective enhancement of true lumen, false lumen and thrombus features in multimodal images of

本发明为实现上述目的,通过以下技术方案实现:The present invention is achieved by the following technical solutions in order to achieve the above object:

一种基于人工智能的多模态医学影像的转换方法,包括:An artificial intelligence-based multimodal medical image conversion method, including:

获取主动脉夹层的OCT图像数据;Obtain OCT image data of aortic dissection;

获取主动脉夹层的IVUS图像数据;Obtain IVUS image data of aortic dissection;

分别对OCT图像数据和IVUS图像数据的背景采用全局阈值分割的方法,提取主动脉的外轮廓,并利用形态学孔洞填充和面积阈值方法将噪声信息进行清理;The background of the OCT image data and the IVUS image data was segmented by global threshold, the outer contour of the aorta was extracted, and the noise information was cleaned up by morphological hole filling and area threshold methods;

通过离散余弦变换的频域信息融合方法,将OCT图像数据和IVUS图像数据进行图像信息融合,获得融合图像。Through the frequency domain information fusion method of discrete cosine transform, the OCT image data and the IVUS image data are fused to obtain the fusion image.

所述获取主动脉夹层的OCT图像数据,为采用激光波长λ为924纳米,纵向分辨率为7微米,穿透能力1.725毫米进行获取。The OCT image data of the aortic dissection is obtained by using a laser with a wavelength λ of 924 nm, a longitudinal resolution of 7 μm, and a penetration capacity of 1.725 mm.

所述获取主动脉夹层的OCT图像数据,包括在轴旋转角度为0°、+90°、-90°、180°状态下进行扫描,获得四角度OCT图像,通过将四角度OCT图像配准获得主动脉夹层的二维OCT横截面图像,对于每个胸主动脉血管,提取200层OCT横截面图像数据,The obtaining of the OCT image data of the aortic dissection includes scanning in a state where the axis rotation angles are 0°, +90°, -90°, and 180° to obtain a four-angle OCT image, which is obtained by registering the four-angle OCT image. Two-dimensional OCT cross-sectional images of aortic dissection, for each thoracic aortic vessel, 200 slices of OCT cross-sectional image data were extracted,

以及,as well as,

所述获取主动脉夹层的IVUS图像数据包括每次胸主动脉血管提取200层IVUS图像数据。The acquiring IVUS image data of aortic dissection includes extracting 200 slices of IVUS image data of the thoracic aorta each time.

所述将OCT图像数据和IVUS图像数据进行图像信息融合,包括以下步骤:Described carrying out image information fusion with OCT image data and IVUS image data, comprises the following steps:

第一,对OCT图像数据和IVUS图像数据进行重叠子块划分,块大小采用8*8;First, the OCT image data and the IVUS image data are divided into overlapping sub-blocks, and the block size is 8*8;

第二,对获得的子块图像分别进行二维离散余弦变换,实现从空间域变换到频域空间中,DCT变换使用公式1:Second, two-dimensional discrete cosine transform is performed on the obtained sub-block images to realize the transformation from the space domain to the frequency domain space. The DCT transform uses formula 1:

Figure BDA0002519568090000031
Figure BDA0002519568090000031

所述h,p分别为频域信息中的坐标,取值范围为h,p=0,1,2……,N-1,The h, p are the coordinates in the frequency domain information respectively, and the value range is h, p=0, 1, 2..., N-1,

N代表划分的图像块大小,N represents the divided image block size,

所述

Figure BDA0002519568090000032
表达式2:said
Figure BDA0002519568090000032
Expression 2:

Figure BDA0002519568090000033
Figure BDA0002519568090000033

逆DCT变换公式3:Inverse DCT transform formula 3:

Figure BDA0002519568090000034
Figure BDA0002519568090000034

所述i,j是空间域坐标变量,取值范围为i,j=0,1,2……,N-1,N代表划分的图像块大小,The i, j are coordinate variables in the space domain, and the value range is i, j=0, 1, 2..., N-1, N represents the size of the divided image block,

根据公式1可以得到包含图像能量信息的主要部分的DC系数,即D(0,0),其余的D(h,p)取值为AC系数;According to formula 1, the DC coefficient containing the main part of the image energy information can be obtained, that is, D(0, 0), and the rest D(h, p) are AC coefficients;

归一化变换系数可由公式4:The normalized transform coefficient can be obtained from Equation 4:

Figure BDA0002519568090000041
Figure BDA0002519568090000041

第三,计算分块图像的方差,每一个分块图像的平均值μ和方差σ2公式如下:Third, calculate the variance of the block image, the average μ and variance σ 2 of each block image are as follows:

公式5:Formula 5:

Figure BDA0002519568090000042
Figure BDA0002519568090000042

公式6:Formula 6:

Figure BDA0002519568090000043
Figure BDA0002519568090000043

第四,基于方差结果比较对融合矩阵更新,Fourth, the fusion matrix is updated based on the variance result comparison,

选取方差较大的子块DCT变换结果的逆变换来更新融合矩阵,直至完成所有子块的融合,即可获得最终的融合图像矩阵。Select the inverse transform of the DCT transform result of the sub-block with larger variance to update the fusion matrix, until the fusion of all sub-blocks is completed, the final fusion image matrix can be obtained.

对比现有技术,本发明的有益效果在于:Compared with the prior art, the beneficial effects of the present invention are:

通过将OCT图像与超声图像进行融合,获取高分辨率的胸主动脉壁面信息,多模态融合可以有效利用OCT和IVUS对组织病灶区域的不同成像优势,有效弥补了单一的OCT图像中弱边缘结构信息的问题。我们针对融合后的结果图像采用深度神经网络的方法进行深层次的特征表征,从而实现最终的真腔、假和血栓的有效检测。有效解决了胸主动脉夹层中假腔和血栓的边缘模糊问题,不仅可以有效的提取胸主动脉夹层中假腔,还可对假腔中血栓进行准确的监测,为疾病的进一步分析奠定了良好的影像基础。By fusing OCT images with ultrasound images, high-resolution thoracic aortic wall information can be obtained. Multimodal fusion can effectively take advantage of the different imaging advantages of OCT and IVUS in tissue lesion areas, and effectively compensate for weak edges in a single OCT image. Problems with structural information. We use the deep neural network method to perform deep feature representation for the fused result images, so as to achieve the final effective detection of true lumen, false and thrombus. It effectively solves the blurred edge of false lumen and thrombus in thoracic aortic dissection. It can not only effectively extract false lumen in thoracic aortic dissection, but also accurately monitor thrombus in false lumen, laying a solid foundation for further analysis of the disease. image base.

附图说明Description of drawings

附图1是本发明图像融合的基本流程示意图。FIG. 1 is a schematic diagram of the basic flow of the image fusion of the present invention.

具体实施方式Detailed ways

下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所限定的范围。The present invention will be further described below in conjunction with specific embodiments. It should be understood that these examples are only used to illustrate the present invention and not to limit the scope of the present invention. In addition, it should be understood that after reading the teaching content of the present invention, those skilled in the art can make various changes or modifications to the present invention, and these equivalent forms also fall within the scope defined by the present application.

下述实施例中所涉及的仪器、试剂、材料等,若无特别说明,均为现有技术中已有的常规仪器、试剂、材料等,可通过正规商业途径获得。下述实施例中所涉及的实验方法,检测方法等,若无特别说明,均为现有技术中已有的常规实验方法,检测方法等。The instruments, reagents, materials, etc. involved in the following examples, unless otherwise specified, are all conventional instruments, reagents, materials, etc. existing in the prior art, and can be obtained through regular commercial channels. The experimental methods, detection methods, etc. involved in the following examples, unless otherwise specified, are all conventional experimental methods, detection methods, etc. in the prior art.

实施例1:一种基于人工智能的多模态医学影像的转换方法Example 1: A method for converting multimodal medical images based on artificial intelligence

简单的利用OCT图像单一的灰度信息,只能对动脉的外壁和内壁实现分割,但却无法实现对动脉内膜至外膜中间的部分细节进行清晰地描述,即无法准确分割出假腔以及血栓的边界。因此,本示例通过多模态信息融合对问题进行解决。Simply using the single grayscale information of the OCT image can only segment the outer and inner walls of the artery, but it cannot clearly describe some details from the intima to the adventitia of the artery, that is, it cannot accurately segment the false lumen and the inner wall. Thrombus borders. Therefore, this example solves the problem through multimodal information fusion.

其具体包括以下步骤:It specifically includes the following steps:

1)获取主动脉夹层的OCT图像数据1) Obtain OCT image data of aortic dissection

采用光学相干断层扫描(OCT)对胸主动脉夹层的血管进行体外成像,In vitro imaging of vessels with thoracic aortic dissection using optical coherence tomography (OCT),

OCT成像技术是目前分辨率最高的血管腔内成像技术,可以有效地显示血管的内壁结构,识别血管内膜的撕裂、斑块组织成分、支架贴合不良等现象,进而可以为动脉壁上发生的病变(如:夹层、微小血栓)以及支架贴合情况等提供强有力的检测工具。OCT imaging technology is currently the highest resolution intraluminal imaging technology, which can effectively display the inner wall structure of blood vessels, identify the tearing of the vascular intima, plaque tissue components, poor stent fit and other phenomena, and then can be used for the arterial wall. The occurrence of lesions (eg, dissection, microthrombosis) and stent fit provide powerful detection tools.

本示例的OCT影像采集方法是:The OCT image acquisition method for this example is:

采用激光波长λ为924纳米,纵向分辨率为7微米,穿透能力1.725毫米进行扫描,The laser wavelength λ is 924 nanometers, the longitudinal resolution is 7 microns, and the penetration ability is 1.725 mm for scanning.

由于主动脉直径较大(约位于2.5-3.5厘米区间),而OCT成像技术穿透力弱(大约在1毫米-2毫米之间),当血管壁由于夹层逐渐增厚时,难以清楚分辨血管外膜结构,即不能直接通过直径进行完整血管壁四周的成像,而需要在多个离散的轴旋转角度方向上成像多次,具体是:以90°为旋转间隔,Due to the large diameter of the aorta (about 2.5-3.5 cm) and the weak penetration of OCT imaging (about 1 mm-2 mm), it is difficult to clearly distinguish the blood vessels when the vessel wall is gradually thickened due to dissection The adventitial structure, that is, cannot be directly imaged around the complete blood vessel wall through the diameter, but needs to be imaged multiple times in multiple discrete axis rotation angles, specifically: with 90° as the rotation interval,

在轴旋转角度为0°、+90°、-90°、180°状态下进行扫描,获得四角度OCT图像,通过将四角度OCT图像配准获得主动脉夹层的二维OCT横截面图像,对于每个胸主动脉血管,提取200层OCT横截面图像数据。Scan under the state of axis rotation at 0°, +90°, -90°, and 180° to obtain a four-angle OCT image, and obtain a two-dimensional OCT cross-sectional image of aortic dissection by registering the four-angle OCT image. For each thoracic aortic vessel, 200 slices of OCT cross-sectional image data were extracted.

2)获取主动脉夹层的IVUS(超声)图像数据2) Obtain IVUS (ultrasound) image data of aortic dissection

IVUS具有高穿透能力,但是却显示分辨率相对低下,通常介于80-120微米。故其无法分辨80微米以下的组织结构,也就是主动脉内壁的内部结构。但是,IVUS对主动脉的管壁和管腔结构能够清晰辨识。故能够依靠其对假腔、血栓的位置进行粗略估计。IVUS has high penetrating power, but shows relatively low resolution, usually between 80-120 microns. Therefore, it cannot distinguish tissue structures below 80 microns, that is, the internal structure of the inner wall of the aorta. However, IVUS can clearly identify the wall and luminal structure of the aorta. Therefore, it can be used to roughly estimate the location of false lumen and thrombus.

每次胸主动脉血管提取200层IVUS图像数据。200 slices of IVUS image data were extracted for each thoracic aortic vessel.

3)分别对OCT图像数据和IVUS图像数据的背景采用全局阈值分割的方法,提取主动脉的外轮廓,并利用形态学孔洞填充和面积阈值方法将噪声信息进行清理。3) The background of the OCT image data and the IVUS image data was segmented by global threshold, the outer contour of the aorta was extracted, and the noise information was cleaned up by the morphological hole filling and area threshold methods.

其步骤是:OCT图像数据/IVUS图像数据——基于全局阈值的二值化分割——形态学孔洞填充——去除背景噪声区域——血管区域提取。The steps are: OCT image data/IVUS image data - binarization segmentation based on global threshold value - morphological hole filling - background noise removal area - blood vessel area extraction.

4)将OCT图像数据和IVUS图像数据进行图像信息融合,获得融合图像4) Fusion of OCT image data and IVUS image data to obtain a fused image

通过对两种成像图像进行信息融合,丰富血管夹层壁内的假腔和血栓区域的特征信息,提高后续环节对特征学习的准确性。By information fusion of the two imaging images, the feature information of the false lumen and the thrombus region in the vessel dissection wall is enriched, and the accuracy of feature learning in subsequent links is improved.

通过离散余弦变换的频域信息融合方法,对OCT图像数据和IVUS图像数据进行融合,具体步骤为:Through the frequency domain information fusion method of discrete cosine transform, the OCT image data and the IVUS image data are fused. The specific steps are:

第一,对OCT图像数据和IVUS图像数据进行重叠子块划分,块大小采用8*8;First, the OCT image data and the IVUS image data are divided into overlapping sub-blocks, and the block size is 8*8;

第二,对获得的子块图像分别进行二维离散余弦变换,实现从空间域变换到频域空间中,DCT变换的公式如下:Second, two-dimensional discrete cosine transform is performed on the obtained sub-block images to realize the transformation from the space domain to the frequency domain space. The formula of DCT transform is as follows:

公式1:Formula 1:

Figure BDA0002519568090000071
Figure BDA0002519568090000071

所述h,p分别为频域信息中的坐标,取值范围为h,p=0,1,2……,N-1,N代表划分的图像块大小。The h and p are the coordinates in the frequency domain information respectively, and the value range is h, p=0, 1, 2..., N-1, and N represents the size of the divided image blocks.

所述

Figure BDA0002519568090000072
表达式如下:said
Figure BDA0002519568090000072
The expression is as follows:

表达式2:Expression 2:

Figure BDA0002519568090000073
Figure BDA0002519568090000073

逆DCT变换定义如下:The inverse DCT transform is defined as follows:

公式3:Formula 3:

Figure BDA0002519568090000074
Figure BDA0002519568090000074

所述i,j是空间域坐标变量,取值范围为i,j=0,1,2……,N-1,N代表划分的图像块大小,根据公式1可以得到包含图像能量信息的主要部分的DC系数,即D(0,0),其余的D(h,p)取值为AC系数。The i, j are coordinate variables in the space domain, and the value range is i, j=0, 1, 2..., N-1, N represents the size of the divided image block, according to formula 1, the main elements including the image energy information can be obtained. Part of the DC coefficients, namely D(0, 0), and the rest D(h, p) are AC coefficients.

归一化变换系数可由以下公式表达:The normalized transform coefficients can be expressed by the following formula:

表达式4:Expression 4:

Figure BDA0002519568090000081
Figure BDA0002519568090000081

第三,计算分块图像的方差Third, calculate the variance of the tiled image

每一个分块图像的平均值μ和方差σ2公式如下:The formulas for the mean μ and variance σ 2 of each segmented image are as follows:

公式5:Formula 5:

Figure BDA0002519568090000082
Figure BDA0002519568090000082

公式6:Formula 6:

Figure BDA0002519568090000083
Figure BDA0002519568090000083

第四,基于方差结果比较对融合矩阵更新Fourth, update the fusion matrix based on the variance result comparison

通过计算得到的方差作为衡量图像对比度信息的指标,方差越大越清晰,图像包含的细节越清晰,特征越丰富。选取方差较大的子块DCT变换结果的逆变换来更新融合矩阵1_fusion中对应位置的元素值,依此类推,直至完成所有子块的融合,即可获得最终的融合图像矩阵。The variance obtained by calculation is used as an index to measure the contrast information of the image. The larger the variance, the clearer the image, the clearer the details and the richer the features. Select the inverse transform of the DCT transform result of the sub-block with large variance to update the element value of the corresponding position in the fusion matrix 1_fusion, and so on, until the fusion of all sub-blocks is completed, and the final fusion image matrix can be obtained.

Claims (4)

1.一种基于人工智能的多模态医学影像的转换方法,其特征在于,包括:1. a method for converting multimodal medical images based on artificial intelligence, is characterized in that, comprises: 获取主动脉夹层的OCT图像数据;Obtain OCT image data of aortic dissection; 获取主动脉夹层的IVUS图像数据;Obtain IVUS image data of aortic dissection; 分别对OCT图像数据和IVUS图像数据的背景采用全局阈值分割的方法,提取主动脉的外轮廓,并利用形态学孔洞填充和面积阈值方法将噪声信息进行清理;The background of the OCT image data and the IVUS image data was segmented by global threshold, the outer contour of the aorta was extracted, and the noise information was cleaned up by morphological hole filling and area threshold methods; 通过离散余弦变换的频域信息融合方法,将OCT图像数据和IVUS图像数据进行图像信息融合,获得融合图像。Through the frequency domain information fusion method of discrete cosine transform, the OCT image data and the IVUS image data are fused to obtain the fusion image. 2.根据权利要求1所述一种基于人工智能的多模态医学影像的转换方法,其特征在于,所述获取主动脉夹层的OCT图像数据,为采用激光波长λ为924纳米,纵向分辨率为7微米,穿透能力1.725毫米进行获取。2. the conversion method of a kind of artificial intelligence-based multimodal medical image according to claim 1, is characterized in that, the described acquisition of the OCT image data of aortic dissection is 924 nanometers for using laser wavelength λ, and the longitudinal resolution is 924 nanometers. The acquisition was performed at 7 microns with a penetration capacity of 1.725 mm. 3.根据权利要求1所述一种基于人工智能的多模态医学影像的转换方法,其特征在于,所述获取主动脉夹层的OCT图像数据,包括在轴旋转角度为0°、+90°、-90°、180°状态下进行扫描,获得四角度OCT图像,通过将四角度OCT图像配准获得主动脉夹层的二维OCT横截面图像,对于每个胸主动脉血管,提取200层OCT横截面图像数据,3. a kind of artificial intelligence-based multi-modal medical image conversion method according to claim 1, is characterized in that, described acquisition of the OCT image data of aortic dissection, including in the axis rotation angle is 0 °, +90 ° Scan at -90°, 180° to obtain four-angle OCT images, and obtain two-dimensional OCT cross-sectional images of aortic dissection by registering the four-angle OCT images. For each thoracic aortic vessel, extract 200 layers of OCT cross-sectional image data, 以及,as well as, 所述获取主动脉夹层的IVUS图像数据包括每次胸主动脉血管提取200层IVUS图像数据。The acquiring IVUS image data of aortic dissection includes extracting 200 slices of IVUS image data of the thoracic aorta each time. 4.根据权利要求1所述一种基于人工智能的多模态医学影像的转换方法,其特征在于,所述将OCT图像数据和IVUS图像数据进行图像信息融合,包括以下步骤:4. the conversion method of a kind of artificial intelligence-based multimodal medical image according to claim 1, is characterized in that, described carrying out image information fusion with OCT image data and IVUS image data, comprises the following steps: 第一,对OCT图像数据和IVUS图像数据进行重叠子块划分,块大小采用8*8;First, the OCT image data and the IVUS image data are divided into overlapping sub-blocks, and the block size is 8*8; 第二,对获得的子块图像分别进行二维离散余弦变换,实现从空间域变换到频域空间中,DCT变换使用公式1:Second, two-dimensional discrete cosine transform is performed on the obtained sub-block images to realize the transformation from the space domain to the frequency domain space. The DCT transform uses formula 1:
Figure FDA0002519568080000021
Figure FDA0002519568080000021
所述h,p分别为频域信息中的坐标,取值范围为h,p=0,1,2……,N-1,The h, p are the coordinates in the frequency domain information respectively, and the value range is h, p=0, 1, 2..., N-1, N代表划分的图像块大小,N represents the divided image block size, 所述
Figure FDA0002519568080000022
表达式2:
said
Figure FDA0002519568080000022
Expression 2:
Figure FDA0002519568080000023
Figure FDA0002519568080000023
逆DCT变换公式3:Inverse DCT transform formula 3:
Figure FDA0002519568080000024
Figure FDA0002519568080000024
所述i,j是空间域坐标变量,取值范围为i,j=0,1,2……,N-1,N代表划分的图像块大小,The i, j are coordinate variables in the space domain, and the value range is i, j=0, 1, 2..., N-1, N represents the size of the divided image block, 根据公式1可以得到包含图像能量信息的主要部分的DC系数,即D(0,0),其余的D(h,p)取值为AC系数;According to formula 1, the DC coefficient containing the main part of the image energy information can be obtained, that is, D(0, 0), and the rest D(h, p) are AC coefficients; 归一化变换系数可由公式4:The normalized transform coefficient can be obtained from Equation 4:
Figure FDA0002519568080000025
Figure FDA0002519568080000025
第三,计算分块图像的方差,每一个分块图像的平均值μ和方差σ2公式如下:Third, calculate the variance of the block image, the average μ and variance σ 2 of each block image are as follows: 公式5:Formula 5:
Figure FDA0002519568080000026
Figure FDA0002519568080000026
公式6:Formula 6:
Figure FDA0002519568080000031
Figure FDA0002519568080000031
第四,基于方差结果比较对融合矩阵更新,Fourth, the fusion matrix is updated based on the variance result comparison, 选取方差较大的子块DCT变换结果的逆变换来更新融合矩阵,直至完成所有子块的融合,即可获得最终的融合图像矩阵。Select the inverse transform of the DCT transform result of the sub-block with large variance to update the fusion matrix, until the fusion of all sub-blocks is completed, the final fusion image matrix can be obtained.
CN202010487136.5A 2020-06-02 2020-06-02 Multimode medical image conversion method based on artificial intelligence Pending CN111640106A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010487136.5A CN111640106A (en) 2020-06-02 2020-06-02 Multimode medical image conversion method based on artificial intelligence

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010487136.5A CN111640106A (en) 2020-06-02 2020-06-02 Multimode medical image conversion method based on artificial intelligence

Publications (1)

Publication Number Publication Date
CN111640106A true CN111640106A (en) 2020-09-08

Family

ID=72330328

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010487136.5A Pending CN111640106A (en) 2020-06-02 2020-06-02 Multimode medical image conversion method based on artificial intelligence

Country Status (1)

Country Link
CN (1) CN111640106A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113096056A (en) * 2021-04-06 2021-07-09 全景恒升(北京)科学技术有限公司 Intravascular image fusion method based on region complementation
CN114240935A (en) * 2022-02-24 2022-03-25 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) Space-frequency domain feature fusion medical image feature identification method and device
WO2022209652A1 (en) * 2021-03-29 2022-10-06 テルモ株式会社 Computer program, information processing method, and information processing device
CN116452484A (en) * 2023-06-20 2023-07-18 深圳英美达医疗技术有限公司 Fusion method, device, computer equipment and storage medium of different medical images

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104376549A (en) * 2014-11-20 2015-02-25 华北电力大学(保定) Intravascular ultrasound image and intravascular-OCT image fusing method
US20170224286A1 (en) * 2014-09-26 2017-08-10 Terumo Kabushiki Kaisha Tomographic image forming apparatus and control method
CN108010035A (en) * 2017-11-07 2018-05-08 深圳市金城保密技术有限公司 Finger vena image segmentation method and its system, terminal based on the detection of direction paddy shape
CN109091167A (en) * 2018-06-29 2018-12-28 东南大学 The prediction technique that Coronary Atherosclerotic Plaque increases
CN110288582A (en) * 2019-06-26 2019-09-27 电子科技大学 An Accurate Segmentation Method of Cell Nuclei

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170224286A1 (en) * 2014-09-26 2017-08-10 Terumo Kabushiki Kaisha Tomographic image forming apparatus and control method
CN104376549A (en) * 2014-11-20 2015-02-25 华北电力大学(保定) Intravascular ultrasound image and intravascular-OCT image fusing method
CN108010035A (en) * 2017-11-07 2018-05-08 深圳市金城保密技术有限公司 Finger vena image segmentation method and its system, terminal based on the detection of direction paddy shape
CN109091167A (en) * 2018-06-29 2018-12-28 东南大学 The prediction technique that Coronary Atherosclerotic Plaque increases
CN110288582A (en) * 2019-06-26 2019-09-27 电子科技大学 An Accurate Segmentation Method of Cell Nuclei

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
徐梦佳: "基于深度卷积神经网络的多模态医学影像分析方法研究" *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022209652A1 (en) * 2021-03-29 2022-10-06 テルモ株式会社 Computer program, information processing method, and information processing device
CN113096056A (en) * 2021-04-06 2021-07-09 全景恒升(北京)科学技术有限公司 Intravascular image fusion method based on region complementation
CN114240935A (en) * 2022-02-24 2022-03-25 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) Space-frequency domain feature fusion medical image feature identification method and device
CN116452484A (en) * 2023-06-20 2023-07-18 深圳英美达医疗技术有限公司 Fusion method, device, computer equipment and storage medium of different medical images
CN116452484B (en) * 2023-06-20 2023-09-26 深圳英美达医疗技术有限公司 Fusion method, device, computer equipment and storage medium of different medical images

Similar Documents

Publication Publication Date Title
CN111640106A (en) Multimode medical image conversion method based on artificial intelligence
JP7236768B2 (en) Calculation method of coronary flow reserve ratio based on myocardial blood flow and CT image
US7215802B2 (en) System and method for vascular border detection
CN112652052B (en) A method and system for three-dimensional reconstruction of coronary artery based on vascular classification criteria
US8805043B1 (en) System and method for creating and using intelligent databases for assisting in intima-media thickness (IMT)
Soepriatna et al. Cardiac and respiratory-gated volumetric murine ultrasound
WO2022105623A1 (en) Intracranial vascular focus recognition method based on transfer learning
Goudarzi et al. Segmentation of arm ultrasound images in breast cancer-related lymphedema: a database and deep learning algorithm
Olender et al. A mechanical approach for smooth surface fitting to delineate vessel walls in optical coherence tomography images
WO2020210948A1 (en) Method and apparatus for acquiring blood flow volume and blood flow velocity of coronary artery
CN110353639B (en) Blood supply area quantification method and system based on blood vessel enhanced radiography
WO2023151280A1 (en) Three-dimensional fusion method and fusion system for dual-mode coronary artery blood vessel images
Gindes et al. Comparison of ex‐vivo high‐resolution episcopic microscopy with in‐vivo four‐dimensional high‐resolution transvaginal sonography of the first‐trimester fetal heart
CN106419843A (en) Method and system for establishing a hypertrophic obstructive cardiomyopathy HOCM heart model
Franquet et al. Identification of the in vivo elastic properties of common carotid arteries from MRI: A study on subjects with and without atherosclerosis
Zheng et al. A deep learning method for motion artifact correction in intravascular photoacoustic image sequence
CN113470060A (en) Coronary artery multi-angle curved surface reconstruction visualization method based on CT image
Yue et al. Speckle tracking in intracardiac echocardiography for the assessment of myocardial deformation
CN114974518A (en) Pulmonary nodule image recognition method and device based on multimodal data fusion
Jansen et al. Local and global distensibility assessment of abdominal aortic aneurysms in vivo from probe tracked 2D ultrasound images
Saba et al. Analysis of carotid artery plaque and wall boundaries on CT images by using a semi-automatic method based on level set model
Katouzian et al. IVUS-histology image registration
Xiao Key Technologies of New Type of Intravascular Ultrasound Image Processing
CN113160116B (en) Method, system and equipment for automatically segmenting inner membrane and outer membrane of left ventricle
Raj et al. An insight into elasticity analysis of common carotid artery using ultrasonography

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: 20200908