CN111640106A - Multimode medical image conversion method based on artificial intelligence - Google Patents
Multimode medical image conversion method based on artificial intelligence Download PDFInfo
- 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
Links
Images
Classifications
-
- 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
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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/10—Segmentation; Edge detection
- G06T7/194—Segmentation; Edge detection involving foreground-background segmentation
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
-
- 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/10104—Positron emission tomography [PET]
-
- 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/10132—Ultrasound image
-
- 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/20048—Transform domain processing
- G06T2207/20052—Discrete cosine transform [DCT]
-
- 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/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
-
- 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/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- 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
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图像中,并在融合的多模态图像中有效的增强真腔、假腔和血栓的特征。
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.
Description
技术领域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:
所述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,
所述表达式2:said Expression 2:
逆DCT变换公式3:Inverse DCT transform formula 3:
所述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:
第三,计算分块图像的方差,每一个分块图像的平均值μ和方差σ2公式如下:Third, calculate the variance of the block image, the average μ and variance σ 2 of each block image are as follows:
公式5:Formula 5:
公式6:Formula 6:
第四,基于方差结果比较对融合矩阵更新,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:
所述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.
所述表达式如下:said The expression is as follows:
表达式2:Expression 2:
逆DCT变换定义如下:The inverse DCT transform is defined as follows:
公式3:Formula 3:
所述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:
第三,计算分块图像的方差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:
公式6:Formula 6:
第四,基于方差结果比较对融合矩阵更新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)
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)
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)
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 |
-
2020
- 2020-06-02 CN CN202010487136.5A patent/CN111640106A/en active Pending
Patent Citations (5)
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)
Title |
---|
徐梦佳: "基于深度卷积神经网络的多模态医学影像分析方法研究" * |
Cited By (5)
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 |