CN114723845B - 以spect为标准的盐水造影肺灌注图像重建方法、系统及设备 - Google Patents
以spect为标准的盐水造影肺灌注图像重建方法、系统及设备 Download PDFInfo
- Publication number
- CN114723845B CN114723845B CN202210649287.5A CN202210649287A CN114723845B CN 114723845 B CN114723845 B CN 114723845B CN 202210649287 A CN202210649287 A CN 202210649287A CN 114723845 B CN114723845 B CN 114723845B
- Authority
- CN
- China
- Prior art keywords
- saline
- spect
- patient
- angiography
- contrast
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 title claims abstract description 234
- 239000011780 sodium chloride Substances 0.000 title claims abstract description 232
- 238000002603 single-photon emission computed tomography Methods 0.000 title claims abstract description 110
- 230000010412 perfusion Effects 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 title claims abstract description 65
- 210000004072 lung Anatomy 0.000 title claims abstract description 51
- 239000012895 dilution Substances 0.000 claims abstract description 110
- 238000010790 dilution Methods 0.000 claims abstract description 110
- 238000002583 angiography Methods 0.000 claims abstract description 82
- 238000013136 deep learning model Methods 0.000 claims abstract description 70
- 238000004422 calculation algorithm Methods 0.000 claims description 19
- 210000000038 chest Anatomy 0.000 claims description 14
- 230000005540 biological transmission Effects 0.000 claims description 11
- 238000003384 imaging method Methods 0.000 claims description 7
- 238000012549 training Methods 0.000 claims description 7
- 230000004927 fusion Effects 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 238000007477 logistic regression Methods 0.000 claims description 3
- 238000007637 random forest analysis Methods 0.000 claims description 3
- 238000004579 scanning voltage microscopy Methods 0.000 claims description 3
- 150000003839 salts Chemical class 0.000 claims 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 3
- 238000012512 characterization method Methods 0.000 claims 1
- 238000012937 correction Methods 0.000 description 7
- 230000002685 pulmonary effect Effects 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 210000000115 thoracic cavity Anatomy 0.000 description 7
- 238000012544 monitoring process Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000009423 ventilation Methods 0.000 description 5
- 206010006322 Breath holding Diseases 0.000 description 4
- 210000004204 blood vessel Anatomy 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004070 electrodeposition Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 description 3
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 108010076039 Polyproteins Proteins 0.000 description 1
- GKLVYJBZJHMRIY-OUBTZVSYSA-N Technetium-99 Chemical compound [99Tc] GKLVYJBZJHMRIY-OUBTZVSYSA-N 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003601 intercostal effect Effects 0.000 description 1
- 210000004731 jugular vein Anatomy 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 210000001321 subclavian vein Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- 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/10108—Single photon emission computed tomography [SPECT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
-
- 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
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Optics & Photonics (AREA)
- Radiology & Medical Imaging (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- High Energy & Nuclear Physics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Nuclear Medicine (AREA)
Abstract
本发明涉及以SPECT为标准的盐水造影肺灌注图像重建方法、系统及设备。所述方法包括：获取患者盐水造影电阻抗稀释曲线；将所述患者盐水造影电阻抗稀释曲线输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；所述预先训练的深度学习模型采用下列方法构建：获取盐水造影电阻抗稀释曲线及同时期的SPECT图；采用深度学习模型将所述盐水造影电阻抗稀释曲线生成基于盐水造影重建的SPECT图；通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。本发明实现了盐水造影肺灌注图像重建。
Description
技术领域
本发明涉及智能医疗技术领域，更具体地，涉及一种以SPECT为标准的盐水造影肺灌注图像重建方法、系统、设备和计算机可读存储介质。
背景技术
EIT作为一种床旁无创、连续、动态、无辐射的先进肺通气监测技术，主要通过局部电极施加微弱电流，感应通气过程中胸腔生物电阻抗变化，再利用相应的成像算法来监测肺通气功能状态。本课题组在国内最早地开展了盐水造影EIT肺灌注成像相关临床应用研究及研发，区域的阻抗稀释曲线最大下降斜率是目前用来建构肺灌注图像的主要特征参数，但是仍有很大的优化空间，有必要进一步提高成像的准确性。
发明内容
基于上述问题，本研究首次提出以SPECT图像作为标准，开发了一种以SPECT为标准的盐水造影肺灌注图像重建方法、系统、设备和计算机可读存储介质，用以提高成像的准确性。
本申请的目的在于提供一种以SPECT为标准的盐水造影肺灌注图像重建方法，所述方法包括：
获取患者盐水造影电阻抗稀释曲线；
将所述患者盐水造影电阻抗稀释曲线输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；
可选的，所述预先训练的深度学习模型采用下列方法构建：
获取盐水造影电阻抗稀释曲线及同时期的SPECT图；
采用深度学习模型将所述盐水造影电阻抗稀释曲线生成基于盐水造影重建的SPECT图；
通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
进一步，所述方法还包括提取所述患者盐水造影电阻抗稀释曲线多参数特征，将所述患者盐水造影电阻抗稀释曲线提取的多参数特征输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；
可选的，所述采用深度学习模型将所述盐水造影电阻抗稀释曲线生成基于盐水造影重建的SPECT图为提取所述盐水造影电阻抗稀释曲线多参数特征，采用深度学习模型将所述提取所述盐水造影电阻抗稀释曲线多参数特征生成基于盐水造影重建的SPECT图；
可选的，所述多参数特征包括下列特征中的一种或几种：曲线最大阻抗下降幅值、曲线下面积、最大斜率、平均传输时间、胸廓形态校正。
所述轮廓形态校正为通过测量EIT监测断层的胸廓前后径、最大横径、胸围大小来进行胸廓形态校正。目前EIT图像主要通过默认的椭圆形轮廓特征来模拟前向模型，通过估算电极位置和电流分布来计算相对阻抗的信息分布，但每个患者胸廓存在个体差异，在此提出通过测量EIT监测断层的胸廓前后径、最大横径、胸围大小等参数特征来进行胸廓形态校正。
优选的，所述多参数特征包括下列特征中的一种或几种：曲线最大阻抗下降幅值、曲线下面积、最大斜率、平均传输时间。
进一步，所述患者盐水造影电阻抗稀释曲线由各个像素点（区域）的阻抗稀释曲线组成，可以分为：患者盐水造影全局阻抗稀释曲线、患者盐水造影心脏区域阻抗稀释曲线、患者盐水造影肺区域阻抗稀释曲线。
进一步，所述深度学习模型为DeepLab或unet，更优选的，所述深度学习模型为DeepLab V3。
进一步，所述提取所述患者盐水造影电阻抗稀释曲线多参数特征为采用机器学习算法提取所述患者盐水造影电阻抗稀释曲线多参数特征，可选的，所述提取所述患者盐水造影电阻抗稀释曲线多参数特征为采用XGBoost、SVM、逻辑回归、随机森林、GBDT、LightGBM或ANN算法中的一种或几种提取所述患者盐水造影电阻抗稀释曲线多参数特征；
可选的，所述方法还包括提取所述患者盐水造影电阻抗稀释曲线多参数特征，将所述多参数特征进行特征融合得到特征融合后的特征，将所述特征融合后的特征输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图。
进一步，所述方法还包括将所述患者盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像，将所述盐水造影EIT灌注图像输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；
可选的，所述预先训练的深度学习模型采用下列方法构建：
获取盐水造影电阻抗稀释曲线及同时期的SPECT图；
将所述盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像；
采用深度学习模型将所述盐水造影EIT灌注图像生成基于盐水造影重建的SPECT图；
通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
进一步，所述患者盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像为利用患者盐水造影电阻抗稀释曲线结合图像重建算法获取盐水造影EIT灌注图像；
可选的，所述利用患者盐水造影电阻抗稀释曲线结合最大斜率法获取盐水造影EIT灌注图像；
优选的，获取患者盐水造影电阻抗稀释曲线，提取所述患者盐水造影电阻抗稀释曲线多参数特征，结合图像重建算法获取盐水造影EIT灌注图像，可选的，所述多参数特征包括下列特征中的一种或几种：曲线最大下降斜率、曲线下面积、平均传输时间、电阻下降幅度。
本申请的目的在于提供一种以SPECT为标准的盐水造影肺灌注图像重建设备，所述设备包括：存储器和处理器；
所述存储器用于存储程序指令；
所述处理器用于调用程序指令，当程序指令被执行时，用于执行实现上述的以SPECT为标准的盐水造影肺灌注图像重建方法。
本申请的目的在于提供一种以SPECT为标准的盐水造影肺灌注图像重建系统，所述系统包括：
获取单元，用于获取患者盐水造影电阻抗稀释曲线；
生成单元，用于将所述患者盐水造影电阻抗稀释曲线输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；
模型训练单元，用于预先训练深度学习模型，采用下列方法构建深度学习模型：获取盐水造影电阻抗稀释曲线及同时期的SPECT图；采用深度学习模型将所述盐水造影电阻抗稀释曲线生成基于盐水造影重建的SPECT图；通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
本申请的目的在于提供一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现上述的以SPECT为标准的盐水造影肺灌注图像重建方法。
本申请的优点：
1.目前盐水造影EIT法虽然具有较高的准确性但还有很大优化空间，本申请以SPECT为标准对盐水造影电阻抗稀释曲线重建，获得基于盐水造影重建的SPECT图；
2.考虑到从盐水造影电阻抗稀释曲线数据维度及映射到SPECT图的困难，本申请还提供一种更优的解决方案，即先结合传统图像重建算法，将患者盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像，再采用深度学习模型将所述盐水造影EIT灌注图像生成基于盐水造影重建的SPECT图；
3.考虑现有技术中传统图像重建一般采用最大斜率法获取盐水造影EIT灌注图像，但是发明人研究发现盐水增强造影产生的阻抗稀释曲线包含了许多肺灌注相关信息，如曲线最大下降斜率、曲线下面积、最大振幅以及平均传输时间等，故而，本申请还提供一种更优的解决方案，即在盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像时，提取所述患者盐水造影电阻抗稀释曲线多参数特征，提高生成盐水造影EIT灌注图像的准确性，进而再采用深度学习模型将所述盐水造影EIT灌注图像生成基于盐水造影重建的SPECT图。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获取其他的附图。
图1是本发明实施例提供的以SPECT为标准的盐水造影肺灌注图像重建方法示意流程图；
图2是本发明实施例提供的以SPECT为标准的盐水造影肺灌注图像重建设备示意流程图；
图3是本发明实施例提供的以SPECT为标准的盐水造影肺灌注图像重建系统示意图；
图4是本发明实施例提供的以SPECT为标准的盐水造影肺灌注图像重建方法过程示意图；
图5是本发明实施例提供的盐水造影肺电阻抗稀释时间曲线反映肺区域灌注的主要特征参数图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。
在本发明的说明书和权利要求书及上述附图中的描述的一些流程中，包含了按照特定顺序出现的多个操作，但是应该清楚了解，这些操作可以不按照其在本文中出现的顺序来执行或并行执行，操作的序号如S101、S102等，仅仅是用于区分开各个不同的操作，序号本身不代表任何的执行顺序。另外，这些流程可以包括更多或更少的操作，并且这些操作可以按顺序执行或并行执行。需要说明的是，本文中的“第一”、“第二”等描述，是用于区分不同的消息、设备、模块等，不代表先后顺序，也不限定“第一”和“第二”是不同的类型。
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有作出创造性劳动前提下所获取的所有其他实施例，都属于本发明保护的范围。
图1是本发明实施例提供的以SPECT为标准的盐水造影肺灌注图像重建方法示意流程图，具体地，所述方法包括如下步骤：
S101：获取患者盐水造影电阻抗稀释曲线；
在一个实施例中，所述患者盐水造影电阻抗稀释曲线包括患者盐水造影全局阻抗稀释曲线、患者盐水造影心脏区域阻抗稀释曲线、患者盐水造影肺区域阻抗稀释曲线；在一个具体实施例中，获取患者盐水造影肺灌注区域阻抗稀释曲线，所述肺灌注区域根据不同划分方法划分成不同子肺灌注区域；可选的，所述子肺灌注区域包括右腹侧、左腹侧、右背侧、左背侧；可选的，所述子肺灌注区域还包括腹外侧、腹内侧、背内侧、背外侧；可选的，所述肺灌注区域，可分为1024（32 X 32 阵矩）个区域（像素点）。
在一个实施例中，所述患者盐水造影电阻抗稀释曲线的获取方法为：首先，呼吸屏气测试，要求最小8秒以上(呼吸机机械通气时，先予适当镇静，调整呼吸机为完全控制通气模式，按呼气或吸气屏气按键10s；自住呼吸患者嘱屏气8秒) ；呼吸屏气测试通过后，方可行盐水造影EIT检查。然后，病人连接肺电阻抗监测仪器，准备10％NaCl 10ml，确认患者已建立中心静脉导管(颈内静脉或锁骨下静脉导管均可)。盐水注射：一般要求2位操作者一起完成，其中一位在确认EIT机器工作正常，启动患者呼吸屏气同时，发号注射盐水的指令；另外一位操作者得到确认指令后，即刻从中心静脉导管快速注射10％NaCL 10ml到患者体内；在整个操作期间EIT监测仪开启记录模式，在注射盐水前2分钟开始持续采集胸部电阻抗信号数据，全过程要求至少持续5 分钟，完整记录屏气期间注射盐水导致肺电阻下降的过程。屏气期间整体电阻曲线开始下降作为盐水进入体内的起点(T0) ，一个动周期后作为盐水进入肺血管的起点(T1) ，整体电阻的最低点作为盐水通过肺血管的终点(T2) ，T0-T1时间段的电阻曲线反映盐水进入右心，不反映肺血管灌注；应用T1-T2时间段的各个肺区域的电阻-时间变化曲线(最大斜率拟合)进行盐水造影EIT肺灌注图像构建。
S102：将所述患者盐水造影电阻抗稀释曲线输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图。
在一个实施例中，所述预先训练的深度学习模型采用下列方法构建：
获取盐水造影电阻抗稀释曲线及同时期的SPECT图；
采用深度学习模型将所述盐水造影电阻抗稀释曲线生成基于盐水造影重建的SPECT图；
通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
在一个实施例中，所述方法还包括提取所述患者盐水造影电阻抗稀释曲线多参数特征，将所述患者盐水造影电阻抗稀释曲线提取的多参数特征输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图。
在一个实施例中，所述同时期的SPECT图获取方法为：首先吸入密闭系统的Xe放射性气体进行肺通气显像，然后通过静脉注射锝99聚合蛋白进行肺血流显像，通过多探头r相机通过SPECT图像，选取与EIT相同平面（第4-5肋间）通气与灌注的SPECT断层图像，以期作为训练的目标。
在一个实施例中，所述采用深度学习模型将所述盐水造影电阻抗稀释曲线生成基于盐水造影重建的SPECT图为提取所述盐水造影电阻抗稀释曲线多参数特征，采用深度学习模型将所述提取所述盐水造影电阻抗稀释曲线多参数特征生成基于盐水造影重建的SPECT图。
在一个实施例中，所述多参数特征包括下列特征中的一种或几种：曲线最大阻抗下降幅值、曲线下面积、最大斜率、平均传输时间、胸廓形态校正。考虑目前EIT图像主要通过默认的椭圆形轮廓特征来模拟前向模型，通过估算电极位置和电流分布来计算相对阻抗的信息分布，但每个患者胸廓存在个体差异，在此提出通过测量EIT监测断层的胸廓前后径、最大横径、胸围大小来进行胸廓形态校正；优选的，所述多参数特征包括下列特征中的一种或几种：曲线最大阻抗下降幅值、曲线下面积、最大斜率、平均传输时间、胸廓形态校正，所述曲线最大阻抗下降幅值又称电阻下降幅度，所述最大斜率又称曲线最大下降斜率，具体见图5。
在一个实施例中，所述深度学习模型为DeepLab或unet，更优选的，所述深度学习模型为DeepLab V3。
在一个实施例中，所述提取所述患者盐水造影电阻抗稀释曲线多参数特征为采用机器学习算法提取所述患者盐水造影电阻抗稀释曲线多参数特征，可选的，所述提取所述患者盐水造影电阻抗稀释曲线多参数特征为采用XGBoost、SVM、逻辑回归、随机森林、GBDT、LightGBM或ANN算法中的一种或几种提取所述患者盐水造影电阻抗稀释曲线多参数特征。
在一个实施例中，所述方法还包括提取所述患者盐水造影电阻抗稀释曲线多参数特征，将所述多参数特征进行特征融合得到特征融合后的特征，将所述特征融合后的特征输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图。
在一个实施例中，获取患者盐水造影电阻抗稀释曲线，将所述患者盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像，将所述盐水造影EIT灌注图像输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图。
在一个实施例中，所述预先训练的深度学习模型采用下列方法构建：
获取盐水造影电阻抗稀释曲线及同时期的SPECT图；
将所述盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像；
采用深度学习模型将所述盐水造影EIT灌注图像生成基于盐水造影重建的SPECT图；
通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
在一个实施例中，所述患者盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像为利用患者盐水造影电阻抗稀释曲线结合图像重建算法获取盐水造影EIT灌注图像；可选的，所述利用患者盐水造影电阻抗稀释曲线结合最大斜率法获取盐水造影EIT灌注图像。
在一个实施例中，获取患者盐水造影电阻抗稀释曲线，提取所述患者盐水造影电阻抗稀释曲线多参数特征，结合图像重建算法获取盐水造影EIT灌注图像，可选的，所述多参数特征包括下列特征中的一种或几种：曲线最大下降斜率、曲线下面积、平均传输时间、电阻下降幅度。
在一个具体实施例中，获取患者盐水造影电阻抗稀释曲线，将所述患者盐水造影电阻抗稀释曲线结合传统图像重建算法生成盐水造影EIT灌注图像，将所述盐水造影EIT灌注图像输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；所述预先训练的深度学习模型采用下列方法构建： 获取盐水造影电阻抗稀释曲线及同时期的SPECT图；将所述盐水造影电阻抗稀释曲线结合传统图像重建算法生成盐水造影EIT灌注图像；采用深度学习模型将所述盐水造影EIT灌注图像生成基于盐水造影重建的SPECT图；通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
在一个实施例中，获取患者盐水造影电阻抗稀释曲线，提取所述患者盐水造影电阻抗稀释曲线多参数特征，结合图像重建算法获取盐水造影EIT灌注图像，将所述盐水造影EIT灌注图像输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；所述预先训练的深度学习模型采用下列方法构建： 获取盐水造影电阻抗稀释曲线及同时期的SPECT图；提取所述患者盐水造影电阻抗稀释曲线多参数特征，结合图像重建算法获取盐水造影EIT灌注图像；采用深度学习模型将所述盐水造影EIT灌注图像生成基于盐水造影重建的SPECT图；通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。基于多参数特征获得的盐水造影EIT灌注图像是一种更为准确反映肺灌注的图像。在一个具体实施例中，所述多参数特征包括下列特征中的一种或几种：曲线最大阻抗下降幅值、曲线下面积、最大斜率、平均传输时间、胸廓形态校正。考虑目前EIT图像主要通过默认的椭圆形轮廓特征来模拟前向模型，通过估算电极位置和电流分布来计算相对阻抗的信息分布，但每个患者胸廓存在个体差异，在此提出通过测量EIT监测断层的胸廓前后径、最大横径、胸围大小来优化模拟的前向模型；优选的，所述多参数特征包括下列特征中的一种或几种：曲线最大阻抗下降幅值、曲线下面积、最大斜率、平均传输时间、胸廓形态校正。
在一个实施例中，见图4，获取患者盐水造影电阻抗稀释曲线，采用XGBoost算法提取所述患者盐水造影电阻抗稀释曲线多参数特征，结合图像重建算法获取盐水造影EIT灌注图像，将所述盐水造影EIT灌注图像输入预先训练的深度学习模型中，融合所述提取所述患者盐水造影电阻抗稀释曲线多参数特征，生成患者的基于盐水造影重建的SPECT图，所述深度学习模型为DeepLab。
图3是本发明实施例提供的一种以SPECT为标准的肺灌注图像重建设备，所述设备包括：存储器和处理器；
所述存储器用于存储程序指令；
所述处理器用于调用程序指令，当程序指令被执行时，用于执行实现上述的以SPECT为标准的盐水造影肺灌注图像重建方法。
图3是本发明实施例提供的一种以SPECT为标准的肺灌注图像重建系统，所述系统包括：
获取单元301，用于获取患者盐水造影电阻抗稀释曲线；
生成单元302，用于将所述患者盐水造影电阻抗稀释曲线输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；
模型训练单元，用于预先训练深度学习模型，采用下列方法构建深度学习模型：获取盐水造影电阻抗稀释曲线及同时期的SPECT图；采用深度学习模型将所述盐水造影电阻抗稀释曲线生成基于盐水造影重建的SPECT图；通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
本发明的一个目的在于提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述的以SPECT为标准的肺灌注图像重建方法。
本验证实施例的验证结果表明，为适应症分配固有权重相对于默认设置来说可以适度改善本方法的性能。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统，装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统，装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：只读存储器（ROM，Read Only Memory）、随机存取存储器（RAM，RandomAccess Memory）、磁盘或光盘等。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上对本发明所提供的一种计算机设备进行了详细介绍，对于本领域的一般技术人员，依据本发明实施例的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (15)
1.一种以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，所述方法包括：
获取患者盐水造影电阻抗稀释曲线；
提取所述患者盐水造影电阻抗稀释曲线多参数特征，所述多参数特征包括下列特征：曲线下面积、最大斜率和平均传输时间；
将所述患者盐水造影电阻抗稀释曲线提取的多参数特征输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；
所述预先训练的深度学习模型采用下列方法构建：
获取盐水造影电阻抗稀释曲线及同时期的SPECT图；
采用深度学习模型将所述盐水造影电阻抗稀释曲线生成基于盐水造影重建的SPECT图；
通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
2.根据权利要求1所述的以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，所述多参数特征还包括下列特征中的一种或几种：曲线最大阻抗下降幅值、胸廓形态校正。
3.根据权利要求1所述的以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，所述患者盐水造影电阻抗稀释曲线包括患者盐水造影全局阻抗稀释曲线、患者盐水造影心脏区域阻抗稀释曲线、患者盐水造影肺区域阻抗稀释曲线。
4.根据权利要求1所述的以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，所述深度学习模型为DeepLab或unet。
5.根据权利要求1所述的以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，所述深度学习模型为DeepLab V3。
6.根据权利要求1或2所述的以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，所述提取所述患者盐水造影电阻抗稀释曲线多参数特征为采用XGBoost、SVM、逻辑回归、随机森林、GBDT、LightGBM或ANN算法中的一种或几种提取所述患者盐水造影电阻抗稀释曲线多参数特征。
7.根据权利要求1或2所述的以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，所述方法还包括提取所述患者盐水造影电阻抗稀释曲线多参数特征，将所述多参数特征进行特征融合得到特征融合后的特征，将所述特征融合后的特征输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图。
8.根据权利要求1-5任意一项所述的以SPECT为标准的肺灌注图像重建方法，其特征在于，所述方法还包括将所述患者盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像，将所述盐水造影EIT灌注图像输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图。
9.根据权利要求1所述的以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，所述预先训练的深度学习模型采用下列方法构建：
获取盐水造影电阻抗稀释曲线及同时期的SPECT图；
将所述盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像；
采用深度学习模型将所述盐水造影EIT灌注图像生成基于盐水造影重建的SPECT图；
通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
10.根据权利要求8所述的以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，所述患者盐水造影电阻抗稀释曲线生成盐水造影EIT灌注图像为利用患者盐水造影电阻抗稀释曲线结合图像重建算法获取盐水造影EIT灌注图像。
11.根据权利要求10所述的以SPECT为标准的盐水造影肺灌注图像重建方法，其特征在于，获取患者盐水造影电阻抗稀释曲线，提取所述患者盐水造影电阻抗稀释曲线多参数特征，结合图像重建算法获取盐水造影EIT灌注图像。
12.一种以SPECT为标准的盐水造影肺灌注图像重建设备，其特征在于，所述设备包括：存储器和处理器；
所述存储器用于存储程序指令；
所述处理器用于调用程序指令，当程序指令被执行时，用于执行实现权利要求1-11任意一项所述的以SPECT为标准的盐水造影肺灌注图像重建方法。
13.一种以SPECT为标准的盐水造影肺灌注图像重建系统，其特征在于，所述系统包括：
获取单元，用于获取患者盐水造影电阻抗稀释曲线，提取所述患者盐水造影电阻抗稀释曲线多参数特征，所述多参数特征包括下列特征：曲线下面积、最大斜率和平均传输时间；
生成单元，用于将所述患者盐水造影电阻抗稀释曲线提取的多参数特征输入预先训练的深度学习模型中，生成患者的基于盐水造影重建的SPECT图；
模型训练单元，用于预先训练深度学习模型，采用下列方法构建深度学习模型：获取盐水造影电阻抗稀释曲线及同时期的SPECT图；采用深度学习模型将所述盐水造影电阻抗稀释曲线生成基于盐水造影重建的SPECT图；通过所述生成基于盐水造影重建的SPECT图与同时期的SPECT图像计算损失，反向传播，更新参数，获得预先训练的深度学习模型。
14.根据权利要求13所述的以SPECT为标准的盐水造影肺灌注图像重建系统，其特征在于，所述多参数特征还包括下列特征中的一种或几种：曲线最大阻抗下降幅值、胸廓形态校正。
15.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现权利要求1-11任意一项所述的以SPECT为标准的盐水造影肺灌注图像重建方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210649287.5A CN114723845B (zh) | 2022-06-10 | 2022-06-10 | 以spect为标准的盐水造影肺灌注图像重建方法、系统及设备 |
| PCT/CN2023/099037 WO2023237034A1 (zh) | 2022-06-10 | 2023-06-08 | 盐水造影电阻抗肺灌注及心脏显像的图像重建及分析方法、系统及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210649287.5A CN114723845B (zh) | 2022-06-10 | 2022-06-10 | 以spect为标准的盐水造影肺灌注图像重建方法、系统及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114723845A CN114723845A (zh) | 2022-07-08 |
| CN114723845B true CN114723845B (zh) | 2022-09-02 |
Family
ID=82233092
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210649287.5A Active CN114723845B (zh) | 2022-06-10 | 2022-06-10 | 以spect为标准的盐水造影肺灌注图像重建方法、系统及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114723845B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023237034A1 (zh) * | 2022-06-10 | 2023-12-14 | 中国医学科学院北京协和医院 | 盐水造影电阻抗肺灌注及心脏显像的图像重建及分析方法、系统及设备 |
| CN115530791B (zh) * | 2022-12-02 | 2023-03-10 | 中国医学科学院北京协和医院 | 基于盐水造影的心内分流电阻抗图像分析方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104582566A (zh) * | 2012-08-20 | 2015-04-29 | 德尔格医疗有限责任公司 | 用于确定肺灌注的量的区域分布的装置 |
| CN105212929A (zh) * | 2014-06-25 | 2016-01-06 | 德尔格制造股份两合公司 | 用于处理断层摄影数据的装置和方法 |
| CN111275786A (zh) * | 2020-02-25 | 2020-06-12 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 结合ML-EM与Unet算法的SPECT断层影像重建方法、装置和设备 |
| CN111340906A (zh) * | 2020-02-25 | 2020-06-26 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种结合ART与UNet算法的SPECT断层影像重建方法、装置和设备 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7117026B2 (en) * | 2002-06-12 | 2006-10-03 | Koninklijke Philips Electronics N.V. | Physiological model based non-rigid image registration |
| EP2293248A1 (en) * | 2009-09-08 | 2011-03-09 | Koninklijke Philips Electronics N.V. | Motion monitoring system for monitoring motion within a region of interest |
-
2022
- 2022-06-10 CN CN202210649287.5A patent/CN114723845B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104582566A (zh) * | 2012-08-20 | 2015-04-29 | 德尔格医疗有限责任公司 | 用于确定肺灌注的量的区域分布的装置 |
| CN105212929A (zh) * | 2014-06-25 | 2016-01-06 | 德尔格制造股份两合公司 | 用于处理断层摄影数据的装置和方法 |
| CN111275786A (zh) * | 2020-02-25 | 2020-06-12 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 结合ML-EM与Unet算法的SPECT断层影像重建方法、装置和设备 |
| CN111340906A (zh) * | 2020-02-25 | 2020-06-26 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种结合ART与UNet算法的SPECT断层影像重建方法、装置和设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114723845A (zh) | 2022-07-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114723845B (zh) | 以spect为标准的盐水造影肺灌注图像重建方法、系统及设备 | |
| JP7475376B2 (ja) | Ｍｒｉのストレイン測定値による心機能の迅速な定量評価 | |
| CN105380598B (zh) | 用于针对动脉狭窄的自动治疗规划的方法和系统 | |
| KR101910233B1 (ko) | 맥관구조를 수치 평가하는 시스템들 및 방법들 | |
| KR101611805B1 (ko) | 환자별 혈류 모델링 방법 및 시스템 | |
| KR101774095B1 (ko) | 제어 방법 및 제어 시스템 | |
| KR102402628B1 (ko) | 인간의 심장 및 심방을 모델링하는 방법 및 시스템 | |
| Johnson et al. | Coronary artery flow measurement using navigator echo gated phase contrast magnetic resonance velocity mapping at 3.0 T | |
| US11039804B2 (en) | Apparatus and method for determining a fractional flow reserve | |
| US11710569B2 (en) | Coronary artery disease metric based on estimation of myocardial microvascular resistance from ECG signal | |
| JP2010221060A (ja) | 壁運動及び心臓の血流のｍｒｉイメージング並びに心臓診断用のシステムを用いる心臓診断 | |
| JP5318855B2 (ja) | 撮像技術 | |
| CN111449657A (zh) | 一种床旁的基于盐水造影肺通气-血流灌注电阻抗断层成像方法 | |
| van Deel et al. | In vivo quantitative assessment of myocardial structure, function, perfusion and viability using cardiac micro-computed tomography | |
| Maï et al. | Effects of breathing and cardiac motion on spatial resolution in the microscopic imaging of rodents | |
| WO2009127827A1 (en) | Imaging technique | |
| CN114723844B (zh) | 经盐水造影校正的搏动灌注图像重建方法、系统及设备 | |
| CN115035208B (zh) | 肺灌注与区域v/q无创成像方法、系统及设备 | |
| RU2016117767A (ru) | Система для управления током крови и способ получения изображений в живом организме | |
| JP6981807B2 (ja) | 医用情報処理装置、ｘ線ｃｔ装置、医用情報処理プログラム、医用情報処理方法及び医用情報処理システム | |
| CN114052764A (zh) | 获取血流储备分数的方法、装置、系统和计算机存储介质 | |
| Pelgrim et al. | Validation of myocardial perfusion quantification by dynamic CT in an ex-vivo porcine heart model | |
| WO2023237034A1 (zh) | 盐水造影电阻抗肺灌注及心脏显像的图像重建及分析方法、系统及设备 | |
| Campbell | Improving Electrical Impedance Tomography Imaging of the Brain | |
| Kircher et al. | Influence of background lung tissue conductivity on the cardiosynchronous EIT signal components: A sensitivity study |
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 | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20220708 Assignee: Shenzhen Security Medical Technology Co.,Ltd. Assignor: PEKING UNION MEDICAL COLLEGE Hospital Contract record no.: X2023980038972 Denomination of invention: Method, system, and equipment for reconstructing saline contrast pulmonary perfusion images using SPECT as the standard Granted publication date: 20220902 License type: Common License Record date: 20230802 |
|
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20220708 Assignee: Dianqi biomedical technology (Suzhou) Co.,Ltd. Assignor: PEKING UNION MEDICAL COLLEGE Hospital Contract record no.: X2023980040775 Denomination of invention: Method, system, and equipment for reconstructing saline contrast pulmonary perfusion images using SPECT as the standard Granted publication date: 20220902 License type: Common License Record date: 20230829 |