CN104287764A - 一种ct灌注成像方法和设备 - Google Patents
一种ct灌注成像方法和设备 Download PDFInfo
- Publication number
- CN104287764A CN104287764A CN201410462793.9A CN201410462793A CN104287764A CN 104287764 A CN104287764 A CN 104287764A CN 201410462793 A CN201410462793 A CN 201410462793A CN 104287764 A CN104287764 A CN 104287764A
- Authority
- CN
- China
- Prior art keywords
- perfusion
- parameters
- target
- algorithm
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000010412 perfusion Effects 0.000 title claims abstract description 657
- 238000003384 imaging method Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 claims abstract description 45
- 150000001875 compounds Chemical class 0.000 claims description 98
- 238000004382 potting Methods 0.000 claims description 98
- 238000010801 machine learning Methods 0.000 claims description 45
- 230000003902 lesion Effects 0.000 claims description 19
- 230000004044 response Effects 0.000 claims description 13
- 238000005481 NMR spectroscopy Methods 0.000 claims description 12
- 230000001575 pathological effect Effects 0.000 claims description 5
- 239000002872 contrast media Substances 0.000 description 25
- 210000001519 tissue Anatomy 0.000 description 18
- 210000004556 brain Anatomy 0.000 description 12
- 238000004458 analytical method Methods 0.000 description 11
- 210000004369 blood Anatomy 0.000 description 11
- 239000008280 blood Substances 0.000 description 11
- 230000017531 blood circulation Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000002347 injection Methods 0.000 description 6
- 239000007924 injection Substances 0.000 description 6
- 239000010410 layer Substances 0.000 description 6
- 230000008901 benefit Effects 0.000 description 5
- 210000004204 blood vessel Anatomy 0.000 description 5
- 238000000205 computational method Methods 0.000 description 5
- 230000036285 pathological change Effects 0.000 description 5
- 231100000915 pathological change Toxicity 0.000 description 5
- 210000003462 vein Anatomy 0.000 description 5
- 210000004185 liver Anatomy 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 230000036770 blood supply Effects 0.000 description 3
- 230000002490 cerebral effect Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000035699 permeability Effects 0.000 description 3
- 210000004227 basal ganglia Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000002224 dissection Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 208000028867 ischemia Diseases 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005316 response function Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 208000003174 Brain Neoplasms Diseases 0.000 description 1
- 241000459479 Capsula Species 0.000 description 1
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 102000010834 Extracellular Matrix Proteins Human genes 0.000 description 1
- 108010037362 Extracellular Matrix Proteins Proteins 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 238000002835 absorbance Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000008499 blood brain barrier function Effects 0.000 description 1
- 210000001218 blood-brain barrier Anatomy 0.000 description 1
- 230000037396 body weight Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 206010008118 cerebral infarction Diseases 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 210000002744 extracellular matrix Anatomy 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000001802 infusion Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 230000004089 microcirculation Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 208000018389 neoplasm of cerebral hemisphere Diseases 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 210000000496 pancreas Anatomy 0.000 description 1
- 210000002381 plasma Anatomy 0.000 description 1
- 238000010992 reflux Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 210000000952 spleen Anatomy 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 210000001103 thalamus Anatomy 0.000 description 1
Classifications
-
- 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
- A61B6/507—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 for determination of haemodynamic parameters, e.g. perfusion CT
-
- 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/032—Transmission computed tomography [CT]
-
- 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/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
-
- 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/486—Diagnostic techniques involving generating temporal series of image data
-
- 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
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- 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
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5247—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0263—Measuring blood flow using NMR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56366—Perfusion imaging
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Human Computer Interaction (AREA)
- Pulmonology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Theoretical Computer Science (AREA)
- Physiology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Hematology (AREA)
- Cardiology (AREA)
Abstract
本发明公开了一种CT灌注成像方法和设备。该方法包括:获取CT设备对当前灌注物扫描得到的灌注数据,所述灌注数据表示时间密度曲线;利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数;按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数;以所述目标灌注参数形成所述当前灌注物的目标灌注图像。通过本发明的技术方案,不仅使得灌注图像体现的灌注状态更加准确,而且也简化了操作。
Description
技术领域
本发明涉及CT设备技术领域,特别是涉及一种CT灌注成像方法和设备。
背景技术
在医疗设备成像中,灌注成像可以提供活组织器官的微循环血流状态(即灌注状态),通过组织器官的灌注状态实现从细胞水平揭示如脑梗塞、肿瘤等病理改变或生理改变,从而有助于外科手术前、后治疗计划的制定。目前,由于成本低、病变区域定位较准确以及普遍适用于各类患者,CT灌注成像成为最常用的一种灌注成像方式。具体地说,CT灌注成像是从静脉团注对比剂之后,对选定层面进行同层动态扫描,以获得该层面每一像素点的时间密度曲线(Time Density Curve,简称TDC),根据TDC曲线利用灌注算法计算灌注参数,再以灌注参数形成灌注图像。对于某一灌注参数的灌注图像来说,其每个像素点的色阶表示该点对应在组织器官上的位置处该灌注参数的大小。例如,对于脑部的灌注形成脑部血容量的灌注图像,在该灌注图像中,每一像素点的色阶表示该点对应的脑部位置上血容量的大小。
目前,现有技术提供了多种灌注算法可以用于对TDC曲线计算灌注参数,而在CT灌注成像时,现有技术中灌注图像是以其中一种灌注算法计算出来的灌注参数来形成的。本发明的发明人经过研究发现,这些灌注算法分别是基于不同的假设条件和计算方法,但这些假设都不完全符合实际,具体地说,有些假设本身不符合实际或难以实现、有些假设仅在特定的时刻或环境下才能发生、有些假设仅在特定组织器官内才合理等等,因此,在现有技术中,无论单独采用哪一种灌注算法计算灌注参数来形成灌注图像,其灌注图像中体现的灌注状态与真实的灌注状态之间也都存在偏差。例如,有些灌注算法计算出的灌注参数在一些特定部位(如缺血的部位)偏差尤其大,有些灌注算法计算出的某个或某些灌注参数被高估或低估,有些灌注算法计算出的灌注参数仅在某个或某些组织器官准确等等。进一步而言,由于现有技术中各种灌注算法的灌注图像都不能体现准确的灌注状态,为了能够得到更准确的灌注状态,医生还需要对各种灌注算法的灌注图像进行对比分析,从而导致查看灌注图像时需要频繁地在不同灌注算法的灌注图像之间切换,造成繁琐的操作。
发明内容
本发明所要解决的技术问题是,提供一种CT灌注成像方法和设备,以解决按照现有技术中采用单一灌注算法计算出的灌注参数形成灌注图像而导致的灌注图像体现的灌注状态不准确以及频繁切换灌注图像的操作繁琐的技术问题。
为解决上述技术问题,本发明实施例提供了一种CT灌注成像方法,该方法包括:
获取CT设备对当前灌注物扫描得到的灌注数据,所述灌注数据表示时间密度曲线;
利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数;
按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数;
以所述目标灌注参数形成所述当前灌注物的目标灌注图像。
可选的,所述各种灌注算法的权重为预先采用机器学习方法对历史灌注物在各种所述灌注算法下计算出的灌注参数与所述历史灌注物的真实灌注参数之间对应关系进行机器学习而得到的。
可选的,所述历史灌注物的真实灌注参数为采用核磁共振对所述历史灌注物进行灌注成像而得到的灌注参数。
可选的,所述方法还包括:
响应于设置系数的输入操作,获取为各种所述灌注算法输入的设置系数;
按照各种所述灌注算法的设置系数,对各个所述初始灌注参数进行重新计算,重新得到目标灌注参数;
以重新得到的目标灌注参数重新形成所述当前灌注物的目标灌注图像。
可选的,所述方法还包括:
根据所述目标灌注图像中各像素点的目标灌注参数,在所述目标灌注图像中标识病变区域,并提供所述病变区域的病变信息。
可选的,所述方法还包括:
响应于对所述目标灌注图像中选定位置上灌注参数的显示请求操作,查找所述选定位置上的目标灌注参数以及各个初始灌注参数并显示。
可选的,所述方法还包括:
分别以各个所述初始灌注参数形成所述当前灌注物在各种所述灌注算法下的灌注图像作为各个初始灌注图像,并显示各个所述初始灌注图像与所述目标灌注图像。
此外,本发明实施例还提供了一种CT灌注成像设备,包括:
第一获取模块,用于获取CT设备对当前灌注物扫描得到的灌注数据,所述灌注数据表示时间密度曲线;
第一计算模块,用于利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数;
第二计算模块,用于按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数;
第一形成模块,用于以所述目标灌注参数形成所述当前灌注物的目标灌注图像。
可选的,所述各种灌注算法的权重为预先采用机器学习方法对历史灌注物在各种所述灌注算法下计算出的灌注参数与所述历史灌注物的真实灌注参数之间对应关系进行机器学习而得到的。
可选的,所述设备还包括:
第二获取模块,用于响应于对设置系数的输入操作,获取为各种所述灌注算法输入的设置系统;
第三计算模块,用于按照各种所述灌注算法的设置系数,对各个所述初始灌注参数进行重新计算,重新得到目标灌注参数;
第二形成模块,用于以重新得到的目标灌注参数重新形成所述当前灌注物的目标灌注图像。
可选的,所述设备还包括:
第二查找模块,用于响应于对所述目标灌注图像中选定位置上灌注参数的显示请求操作,查找所述选定位置上的目标灌注参数以及各个初始灌注参数;
第一显示模块,用于显示所述选定位置上的目标灌注参数以及各个所述初始灌注参数。
与现有技术相比,本发明具有以下优点:
根据本发明实施例的技术方案,可以预先采集一些已知真实灌注参数的历史灌注物,采用机器学习方法对历史灌注物在各种所述灌注算法下计算出的灌注参数与所述历史灌注物的真实灌注参数之间对应关系进行机器学习,从而得到机器学习方法中各种灌注算法的权重,当需要当前灌注物形成灌注图像时,可以先获取CT设备对当前灌注物扫描得到的表示时间密度曲线的灌注数据,然后利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数,再按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数,从而以所述目标灌注参数形成所述当前灌注物的目标灌注图像。因此,由于灌注图像是以目标灌注参数来形成的,并且目标灌注参数又是综合了各种灌注算法计算出的初始灌注参数而计算得到的,使得目标灌注图像的形成综合了各种灌注算法,因此,各种灌注算法不同的假设条件可以互相弥补彼此间不符合实际的情况,使得在目标灌注图像中目标灌注参数能够更加接近于真实的灌注参数,从而不仅使得灌注图像体现的灌注状态更加准确,也避免了医生为了对比各个灌注状态不准确的灌注图像而频繁地切换各种灌注算法的灌注图像,简化了操作。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明中CT灌注成像方法一实施例的流程图;
图2为本发明实施例中修改目标灌注图像中各灌注算法权重一实施方式的流程图;
图3为本发明中CT灌注成像设备一实施例的结构图;
图4为本发明中CT灌注成像设备另一实施例的结构图;
图5为本发明中CT灌注成像设备又一实施例的结构图;
图6为本发明中CT灌注成像设备再一实施例的结构图;
图7为本发明中CT灌注成像设备又再一实施例的结构图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的发明人经过长期的研究发现,现有技术中虽然提供了多种灌注算法用于计算灌注参数,但是,由于灌注图像是以单独一种灌注算法计算出来的灌注参数来形成的,而所有的灌注算法采用的假设条件都不完全符合实际,这就使得现有技术用于形成灌注图像的灌注参数都与真实的灌注参数之间存在偏差,从而导致灌注图像体现的灌注状态都与真实的状态之间存在偏差,因此,医生在使用现有的灌注图像时,需要为其利用多种灌注算法计算出的灌注参数形成多个灌注图像,使得医生通过对比分析多个灌注图像来确定更准确的灌注状态,这就造成了繁琐的操作。发明人进一步研究后发现,实际上,各种灌注算法的假设条件与实际不符合的情况是不同的,因此,各种灌注算法的假设条件能够互相弥补彼此间不符合实际的情况,也正是因此,医生能够通过综合考虑各种灌注算法下得到的灌注图像确定更接近于真实的灌注状态。由此可见,为了获得更接近于真实灌注状态的灌注图像而实现操作的简化,可以综合考虑多种灌注算法计算出的灌注参数来形成一个灌注图像。
基于上述分析,本发明的主要思想在于:为了实现综合多种灌注算法计算出的灌注参数来形成一个更接近于真实灌注状态的灌注图像,可以采用机器学习的方式来建立各种灌注算法的灌注参数与真实灌注参数之间的对应关系,以便利用该机器学习算法实现将各种灌注参数计算出的初始灌注参数综合成更接近于真实的目标灌注参数,从而以目标灌注参数形成的目标灌注图像可以体现出更接近于真实的灌注状态,并进一步避免了频繁更换灌注算法以切换各种灌注图像的操作,避免了操作误差,简化了操作。具体地说,可以预先采集一些已知真实灌注参数的历史灌注物,采用机器学习方法对历史灌注物在各种灌注算法下计算出的灌注参数与历史灌注物的真实灌注参数之间建立对应关系,从而得到机器学习方法中各种灌注算法的权重,当需要当前灌注物形成灌注图像时,对于当前灌注物扫描得到的灌注数据,可以利用至少两种灌注算法,分别对灌注数据进行计算,得到各种灌注算法下的灌注参数作为各个初始灌注参数,再按照各种灌注算法的权重,对各个初始灌注参数进行计算,得到目标灌注参数,从而以目标灌注参数形成所述当前灌注物的目标灌注图像。
在本发明实施例的一个示例性应用场景中,例如可以应用在CT设备上。具体地,在这一应用场景示例中,采用机器学习方法对历史灌注物在各种所述灌注算法下计算出的灌注参数与历史灌注物的真实灌注参数之间对应关系进行机器学习而为各种灌注算法得到的权重可以预置在CT设备中,在需要对当前灌注物进行灌注成像时,CT设备可以对当前灌注物扫描来得到灌注数据,利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数,然后,CT设备可以按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数,再后,CT设备以所述目标灌注参数形成所述当前灌注物的目标灌注图像。
在本发明实施例的又一个示例性应用场景中,例如可以应用在与CT设备相连的PC等计算设备上。具体地,在这一应用场景示例中,采用机器学习方法对历史灌注物在各种所述灌注算法下计算出的灌注参数与历史灌注物的真实灌注参数之间对应关系进行机器学习而为各种灌注算法得到的权重可以预置在该计算设备中,在需要对当前灌注物进行灌注成像时,CT设备对当前灌注物扫描得到的灌注数据之后,该计算设备可以从CT设备中接收该灌注数据,并利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数,然后,该计算设备可以按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数,再后,该计算设备以所述目标灌注参数形成所述当前灌注物的目标灌注图像。
在介绍了本发明的主要思想以后,下面结合附图,详细说明本发明的各种非限制性实施方式。
参见图1,示出了本发明中CT灌注成像方法一实施例的流程图。在本实施例中,例如具体可以包括如下步骤:
S101、获取CT设备对当前灌注物扫描得到的灌注数据,所述灌注数据表示时间密度曲线。
具体实现时,在一些实施方式中,采用CT设备平扫当前待灌注物时,确定可以做灌注扫描后,结合病人体重、健康程度,病变位置、类型等信息,可以由医生确定造影剂剂量及造影剂注射速率,设定扫描开始时间及扫描时间、每层图像的间隔时间等,然后根据以上参数,CT设备可以自动匹配最佳灌注计算方法,或者由医生手动选择最佳灌注计算方法。在另一些实施方式中,也可以将病人相关信息输入设备,由CT设备给出建议的造影剂剂量、造影剂注释速率、确定各项扫描参数及最佳灌注计算方法等。可以理解的是,最佳灌注计算方法的选择会影响到后续目标灌注参数计算过程中的权重。
在当前灌注物(如病患的脑部、肝部等)静脉团注射对比剂之后,CT设备在选定层面对当前灌注物进行动态扫描,获得灌输数据。其中,灌注数据具体表示的是选定层面内每1个像素点的时间密度曲线,也即,通过该灌注数据,可以绘制出选定层面内每1个像素点的时间密度曲线。此外,选定层面可以是一层或多层,通常情况下采用的是单层。可以理解的是,对于不同的组织器官进行灌注时,可以依据组织器官分布来选择选定层面,例如,对于脑灌注来说,可以选择基底节层面作为选定层面,这样就使得选定层面中较为全面地包括了丘脑、基底节、内囊以及大脑前、中、后动脉。对于选定多层的情况,可以在某一层选择组织器官的主要供血血管,在某一层选择血液流出该组织的主要血管,用于病变所在层的灌注参数的计算。其中,供血血管和血液流出组织的血管以及病变可以在同一层,也可以分别在不同层。对于血管及病变区域的选择,可以是由医生手动选择,或者也可以是由系统自动确定,或者还可以是在系统自动确定不够准确的情况下,由医生手动调整。在同时计算多层的灌注参数的情况下,可以形成空间上三维的灌注结果,向医生提供病变体积等信息。
S102、利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数。
在本实施例中,为形成灌注图像,所使用的灌注参数例如可以是血容量、血流量、平均通过时间、达峰时间、表面通透性等参数中的一个或多个。可以理解的是,所使用的灌注参数可以针对所要灌注的组织器官来进行选择。例如,对于脑灌注来说,所使用的灌注参数可以是脑部的血容量、血流量、平均通过时间或达峰时间等参数。又如,对于肝脏灌注来说,所使用的灌注参数可以是肝脏的血容量、血流量、平均通过时间或表面通透性等参数。再如,对于其他组织器官的灌注,所使用的灌注参数可以是该组织器官的血容量、血流量、平均通过时间或表面通透性等参数。
在本实施例中,用于计算初始灌注参数的灌注算法,例如可以是现有灌注算法中的任意两种或两种以上,具体地,可以包括以下算法的任意两种或两种以上:瞬间法、最大斜率法、去卷积法、patlak分析法等。可以理解的是,用于计算灌注参数的灌注算法可以是针对所要灌注的组织器官来选择。例如,对于脑灌注来说,用于计算灌注参数的灌注算法可以是最大斜率法、去卷积法等。又如,对于肝脏灌注来说,用于计算灌注参数的灌注算法可以是“一室两进”模式的去卷积法、将脾脏强化达峰时间作为肝动脉期与门静脉期的分界点的最大斜率法等。再如,对于其他组织器官的灌注,用于计算灌注参数的灌注算法可以是最大斜率法、去卷积法、patlak分析法等。
需要说明的是,不同的灌注算法是基于不同的假设条件,因此,利用不同灌注算法计算出来的初始灌注参数,也与真实的灌注参数存在不同的偏差。
例如,对于瞬间法,其理论基础是对比剂稀释理论,即在没有对比剂外渗和消除对比剂再循环的情况下,可以根据时间密度曲线计算出脑血容量。假定对比剂单一通过,没有对比剂再循环影响的推测值。瞬间法的优点是不需要假设无静脉引流,因此可以将血容量与平均通过时间之比当作血流量,计算简单。但是,在一般情况下其组织通过时间的方差与平均通过时间之比在大脑中为常数的这种假设是不符合实际的,尤其是在发生缺血的部位。
又如,对于斜率法,其是由Peters提出,属于单房室的一种灌注算法,并且其能够独立于具体的模型结构。假设血管内壁完整,对比剂经由输入动脉进入组织器官,然后快速渗透到毛细血管和细胞外间质中,并且不与组织发生代谢作用,经过很短的一段时间后开始从一条静脉流出。假设所有注入的对比剂均留在脑内血管内而无静脉流出,则组织血流量等于组织TDC的最大初始斜率与动脉TDC的峰值高度之比值。斜率法的优点是概念简单,便于理解,但其主要弊端是脑组织时间密度曲线的最大初始斜率处对比剂无静脉流出的假设很难实现。如果在达到最大初始斜率峰值之前对比剂有明显的静脉流出,血流量将被低估。另外,由于最大初始斜率取决于对比剂团注的速度和病人的心脏输出功能。因此,要缩短最大初始斜率由起始到达峰值的时间,就只有提高对比剂的注射速度和对比剂剂量。为保证计算的准确性而要求较高注射速率(10-20ml/s),必须使用静脉插管,但这种方式不适于心功能不良以及静脉比较脆弱的病人,所以斜率法在临床应用中受到限制。
再如,对于多时图解分析法,其采用的是将感兴趣组织区域和输入动脉血浆的对比剂时间密度曲线进行变换处理,在满足一定条件的情况下,新合成的曲线有一部分近似于一条直线,对数据曲线能够进行线性拟合。拟合后生成的直线斜率在不同的图解分析方法中,分别表征对比剂的吸收率或者分布容积。多时图解分析方法的优点在于,其能够独立于任何特定的模型结构。一般情况下,分析所采用的拟合方法在肿瘤组织区域能够根据观测值拟合得更好,观测值的离散度更小,即观测值更稳定。但是,多时图解分析法的缺点在于,其采用的假设条件是存在一个对比剂不可回流的反应或步骤,使得对比剂或者标记物不可逆向流出组织,但在实际情况下,这样的情况只能在初始时经历剧烈的对比剂浓度变化,且组织对比剂浓度曲线下降得足够缓慢的组织房室模型中才能够发生。此外,多时图解分析法还忽略了血管容积。
又再如,对于去卷积算法,其是由Cenic在1999年提出的,其是将所研究的系统看作线性时不变的,在理想状态下,根据瞬时注射单位质量的对比剂所产生的组织增强曲线,得到组织的脉冲响应函数。该脉冲响应函数可以认为反映了物质在组织中随着时间的分布情况。去卷积法的优点在于,其要求对比剂的注射速度不必很高,一般3.5~4.0ml/s即可,因此其在临床上得到了广泛的应用。但是,去卷积法的缺点在于,其采用的假设条件要求对比物是不扩散的,虽然在脑部对比物不扩散是合理的假设,但是对其他器官就不正确了,并且这也不符合脑部肿瘤时血脑屏障损坏的情况。此外,去卷积法对噪声特别敏感,因此,也必须用有效的方法抑制噪声才能保证去卷积法的计算值准确。
此外,对于延迟敏感的灌注算法,其假设条件是在动脉灌注开始前组织中不存在灌注,因此,相对于延迟不敏感的灌注算法,延迟敏感的灌注算法会扩大对有动脉阻塞组织的估计,比如,会扩大对脑部中风区域的估计。
S103、按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数。
其中,所述各种灌注算法的权重可以为预先采用机器学习方法对历史灌注物在各种所述灌注算法下计算出的灌注参数与所述历史灌注物的真实灌注参数之间对应关系进行机器学习而得到的。具体实现时,可以预先采用一些已知真实灌注参数的历史灌注物及利用各种灌注算法为历史灌注物计算出的历史灌注参数,然后对历史灌注物在各种灌注算法下计算出的各个灌注参数与历史灌注物的真实灌注参数之间的对应关系进行机器学习,从而建立出一个表示灌注物在各种灌注算法下计算出的灌注参数与其真实灌注参数之间对应关系的机器学习模型,而已建立的机器学习模型即是由各种灌注算法的权重所组成,而当需要为未知真实灌注参数的当前灌注物形成灌注图像时,可以采用从已建立的机器学习模型中获得的各种灌注算法的权重来对当前灌注物在各种灌注算法下的初始灌注参数进行计算,得到目标灌注参数,则该目标灌注参数相对于各初始灌注参数就更接近于当前灌注物的真实灌注参数。
具体地说,通过机器学习方法建立的机器学习模型是可以用于表示多个输入变量的变量值与多个输出变量的变量值之间对应关系的计算工具,具体到本实施例,机器学习模型表示的可以是灌注物在各种灌注算法下计算出的灌注参数与该灌注物的真实灌注参数之间的对应关系,也即,在该机器学习模型中,输入变量是灌注物在各种灌注算法下计算出的各个灌注参数,输出变量为灌注物的真实灌注参数。由于历史灌注物在各种灌注算法下的各个灌注参数以及历史灌注物的真实灌注参数均是已知的,通过对历史灌注物在各种灌注算法下的灌注参数与历史灌注物的真实灌注参数之间对应关系的训练学习,可以建立表示任意灌注物在各种灌注算法下的灌注参数与该灌注物的真实灌注参数间对应关系的机器学习模型,从而以当前灌注物在各种灌注算法下的各个初始灌注参数作为输入变量,该机器学习模型就可以计算出作为输出变量的目标灌注参数,使得该目标灌注参数更接近于当前灌注物的真实灌注参数。因此,该机器学习模型,可以是表示一个针对多个变量与一个输出变量的机器学习模型。
此外,在本实施例中,机器学习方法可以采用任意一种能够实现机器学习功能的算法,例如可以是神经元网络、支持向量机、AdaBoost算法、贝叶斯网络等,本实施例对此不做限定。
可以理解的是,为了通过机器学习得到各种灌注算法的权重,需要利用历史灌注物在各种灌注算法下的灌注参数及其真实灌注参数。对于历史灌注物在各种灌注算法下的灌注参数,例如可以利用各种灌注算法对历史灌注物进行CT灌注成像来获得。而对于历史灌注物的真实灌注参数,实际上通过对历史灌注物进行解剖分析才能够获得最准确的真实灌注参数。但是,通常机器学习需要大量的训练数据,而能够通过解剖分析的手段获知真实灌注参数的历史灌注物数量却十分有限。因此,为了采集到已知真实灌注参数的历史灌注物的数量足够建立机器学习模型,由于核磁共振对灌注成像非常接近于真实的灌注状态,在本实施例的一些实施方式中,可以采用核磁共振对历史灌注物进行灌注成像而得到的灌注参照作为其真实灌注参数,也即,历史灌注物的真实灌注参数例如可以为采用核磁共振对所述历史灌注物进行灌注成像而得到的灌注参数。
需要说明的是,针对不同的组织器官,可以建立同一机器学习模型来使用,或者,也可以分别建立不同的机器学习模型来使用。具体地,在本实施例的一些实施方式中,可以根据组织器官的供血特点来确定哪些组织器官可以采用同一机器学习模型。例如,对于肾脏和胰脏来说,两者供血特点比较接近,均为单系统动脉供血器官,因此两者可以采用同一机器学习模型来进行灌注成像,该机器学习模型中涉及的灌注算法可以采用单输入模型。此外,类似地,针对不同的疾病,可以建立同一机器学习模型来使用,或者,也可以建立不同的机器学习模型来使用。
S104、以所述目标灌注参数形成所述当前灌注物的目标灌注图像。
具体实现时,可以以选定层面每1个像素点的目标灌注参数来确定目标灌注图像中相应像素点的灰度或色阶,以此生成目标灌注图像。在目标灌注图像中,各位置的灰度或色阶即体现出该位置上更接近于真实灌注参数的目标灌注参数,从而使得该目标灌注图像可以体现出更接近于真实的灌注状态。此时,可以响应于用户对图像的显示请求操作,显示形成的目标灌注图像,从而使得医生能够从目标灌注图像上查看到更准确的灌注状态。
可以理解的是,以机器学习模型综合各种灌注算法的初始灌注参数得到的目标灌注参数来形成目标灌注图像,虽然可以使得目标灌注图像相对于以初始灌注参数形成的灌注图像来说更接近于真实的灌注状态,但是,由于各种实际情况之间存在的差异,有时目标灌注图像所体现的灌注状态可能还是不够理想。为了使得目标灌注图像能够更加理想地体现真实的灌注状态,在本实施例的一些实施方式中,还可以允许用户通过为各种灌注算法自行输入设置系数来实现对各初始灌注参数进行计算,以实现对目标灌注图像的调整,从而得到体现更准确灌注状态的目标灌注图像。具体地,如图2所示,本实施例中还可以包括:
S201、响应于设置系数的输入操作,获取为各种所述灌注算法输入的设置系数。
具体实现时,在向用户显示目标灌注图像的同时,可以为用户提供一种为调整目标灌注图像而输入设置系数的操作方式,当用户执行了对权重的修改操作时,可以获取用户对各种灌注算法输入的设置系数,以用于重新生成目标灌注图像。
S202、按照各种所述灌注算法修改过的设置系数,对各个所述初始灌注参数进行重新计算,重新得到目标灌注参数。
具体实现时,例如可以以用户对各灌注算法输入的设置系数作为各自初始灌注参数在目标灌注参数中占据的比例系数,将各个初始灌注参数与其设置系数的乘积之和作为目标灌注参数。
S203、以重新得到的目标灌注参数重新形成所述当前灌注物的目标灌注图像。
在重新形成目标灌注图像之后,可以将重新形成的目标灌注图像向用户显示。可以理解的是,如果重新形成的目标灌注图像依然不够理想,则用户可以再次对各灌注参数进行修改,从而再重新形成目标灌注图像,直至目标灌注图像足够理想。
此外,在一些实施方式中,对于为各灌注算法输入过的一组设置系数,可以将其保存下来,以便后续用户再次调整目标灌注参数时可以将曾经输入过的设置系数组提供给用户选择,从而简化用户的操作。此外,对于以设置系数对各初始灌注参数计算出的目标灌注参数,可以将其各初始灌注参数与目标灌注参数作为一组训练数据,以用于对机器学习模型的再次训练。
接着返回图1。
在本实施例的另一些实施方式中,为了帮助用户更容易地定位病变位置,还可以在目标灌注图像中依据各像素点的目标灌注参数来确定哪些位置可能已发生病变并呈现这些确定出的病变位置。具体地,本实施例例如还可以包括:根据所述目标灌注图像中各像素点的目标灌注参数,在所述目标灌注图像中将所述标识像素点组成的区域标识为病变区域,并提供所述病变区域的病变信息。更具体地,可以考虑多种不同的目标灌注参数,在目标灌注图像中查找这些目标灌注参数都位于其各自预设范围内的像素点作为标识像素点,再以这些标识像素点在目标灌注图像中组成的区域标识为病变区域。其中,为了标识不同疾病的病变区域,可以选用不同种类的目标灌注参数,和/或,同一种类的目标灌注参数也可以选用不同的预设范围。此外,所述病变信息,例如可以包括依据其病变区域所分析得到的病变情况、病变程度和/或病变区域面积等。
可以理解的是,用户有时需要分析比对分别以各灌注算法下初始灌注参数形成的初始灌注图像以及综合各灌注算法形成的目标灌注图像,这样就会导致需要频繁地为当前灌注物反复进行灌注成像,造成操作的繁琐以及系统的繁忙。为了在满足用户分析比对各灌注图像的同时避免频繁地为当前灌注物反复进行灌注成像,在本实施例的又一些实施方式中,还可以为当前灌注物并行显示各种灌注算法分别形成的初始灌注图像以及综合各灌注算法形成的目标灌注图像。具体地,本实施例例如还可以包括:分别以各个所述初始灌注参数形成所述当前灌注物在各种所述灌注算法下的灌注图像作为各个初始灌注图像,并显示各个所述初始灌注图像与所述目标灌注图像。其中,各个初始灌注图像与目标灌注图像的并行显示,例如可以是将各个初始灌注图像与目标灌注图像同时显示在显示界面上;又如可以是同时形成了各个初始灌注图像与目标灌注图像之后,为用户提供一种选择显示哪一或哪些灌注图像的操作方式,并响应用户对图像的选择操作来显示用户从各初始灌注图像及目标灌注图像中选定的灌注图像。
此外,用户有时需要分析比对同一位置上各种灌注算法下得到的初始灌注参数以及综合各灌注算法得到的目标灌注参数,如果灌注参数仅在其灌注图像上提供显示,这样就会导致需要频繁地切换各灌注图像,造成操作的繁琐以及系统的繁忙。为了在满足用户分析比对各灌注参数的同时避免频繁地切换灌注图像,在本实施例的再一实施方式中,还可以在一个灌注图像上并行显示所有灌注图像相同位置上的灌注参数。具体地说,以目标灌注图像上的并行显示为例,本实施例例如还可以包括:响应于对所述目标灌注图像中选定位置上灌注参数的显示请求操作,查找所述选定位置上的目标灌注参数以及各个初始灌注参数并显示。其中,在并行显示各灌注参数时,可以将选定位置上的目标灌注参数及各初始灌注参数置入一个表格中进行显示,从而可以便于用户比对各灌注参数。可以理解的是,虽然此处是以目标灌注图像上并行显示各灌注参数为例进行介绍的,但结合前述并行显示各初始灌注图像及目标灌注图像的实施方式,各灌注参数也可以是在各个初始灌注图像中的任意图像上并行显示。对于并行显示各灌注参数时用户是在哪一灌注图像上触发并显示的,本实施例不做限定。
通过本实施例的技术方案,由于灌注图像是以目标灌注参数来形成的,并且目标灌注参数又是综合了各种灌注算法计算出的初始灌注参数而计算得到的,使得目标灌注图像的形成综合了各种灌注算法,因此,各种灌注算法不同的假设条件可以互相弥补彼此间不符合实际的情况,使得在目标灌注图像中目标灌注参数能够更加接近于真实的灌注参数,从而不仅使得灌注图像体现的灌注状态更加准确,也避免了医生为了对比各个灌注状态不准确的灌注图像而频繁地切换灌注图像,简化了操作。
在介绍了本发明示例性实施方式的方法之后,接下来,对本发明示例性实施方式的、用于CT灌注成像的设备进行介绍。
参见图3,示出了本发明中CT灌注成像设备一实施例的结构图。在本实施例中,所述设备例如具体可以包括:
第一获取模块301,用于获取CT设备对当前灌注物扫描得到的灌注数据,所述灌注数据表示时间密度曲线;
第一计算模块302,用于利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数;
第二计算模块303,用于按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数;
第一形成模块304,用于以所述目标灌注参数形成所述当前灌注物的目标灌注图像。
在本实施例的一些实施方式中,所述各种灌注算法的权重可以为预先采用机器学习方法对历史灌注物在各种所述灌注算法下计算出的灌注参数与所述历史灌注物的真实灌注参数之间对应关系进行机器学习而得到的。
在本实施例的另一些实施方式中,所述历史灌注物的真实灌注数据例如可以为采用核磁共振对所述历史灌注物进行灌注成像而得到的灌注参数。
参见图4,示出了本发明中CT灌注成像设备另一实施例的结构图。在本实施例中,除了图3所示的所有结构外,所述设备例如还可以包括:
第二获取模块401,用于响应于对设置系数的输入操作,获取为各种所述灌注算法输入的设置系数;
第三计算模块402,用于按照各种所述灌注算法的设置系数,对各个所述初始灌注参数进行重新计算,重新得到目标灌注参数;
第二形成模块403,用于以重新得到的目标灌注参数重新形成所述当前灌注物的目标灌注图像。
参见图5,示出了本发明中CT灌注成像设备又一实施例的结构图。在本实施例中,除了图3所示的所有结构外,所述设备例如还可以包括:
标识模块501,用于根据所述目标灌注图像中各像素点的目标灌注参数,在所述目标灌注图像中标识病变区域;
提供模块502,用于提供所述病变区域的病变信息。
参见图6,示出了本发明中CT灌注成像设备再一实施例的结构图。在本实施例中,除了图3所示的所有结构外,所述设备例如还可以包括:
第二查找模块601,用于响应于对所述目标灌注图像中选定位置上灌注参数的显示请求操作,查找所述选定位置上的目标灌注参数以及各个初始灌注参数;
第一显示模块602,用于显示所述选定位置上的目标灌注参数以及各个所述初始灌注参数。
参见图7,示出了本发明中CT灌注成像设备又再一实施例的结构图。在本实施例中,除了图3所示的所有结构外,所述设备例如还可以包括:
第三形成模块701,用于分别以各个所述初始灌注参数形成所述当前灌注物在各种所述灌注算法下的灌注图像作为各个初始灌注图像;
第二显示模块702,用于显示各个所述初始灌注图像与所述目标灌注图像。
通过本发明设备实施例的技术方案,由于灌注图像是以目标灌注参数来形成的,并且目标灌注参数又是综合了各种灌注算法计算出的初始灌注参数而计算得到的,使得目标灌注图像的形成综合了各种灌注算法,因此,各种灌注算法不同的假设条件可以互相弥补彼此间不符合实际的情况,使得在目标灌注图像中目标灌注参数能够更加接近于真实的灌注参数,从而不仅使得灌注图像体现的灌注状态更加准确,也避免了医生为了对比各个灌注状态不准确的灌注图像而频繁地切换灌注图像,简化了操作。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
对于设备实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的设备实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上所述仅是本发明的具体实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (11)
1.一种CT灌注成像方法,其特征在于,包括:
获取CT设备对当前灌注物扫描得到的灌注数据,所述灌注数据表示时间密度曲线;
利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数;
按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数;
以所述目标灌注参数形成所述当前灌注物的目标灌注图像。
2.根据权利要求1所述的方法,其特征在于,所述各种灌注算法的权重为预先采用机器学习方法对历史灌注物在各种所述灌注算法下计算出的灌注参数与所述历史灌注物的真实灌注参数之间对应关系进行机器学习而得到的。
3.根据权利要求1所述的方法,其特征在于,所述历史灌注物的真实灌注参数为采用核磁共振对所述历史灌注物进行灌注成像而得到的灌注参数。
4.根据权利要求1所述的方法,其特征在于,还包括:
响应于设置系数的输入操作,获取为各种所述灌注算法输入的设置系数;
按照各种所述灌注算法的设置系数,对各个所述初始灌注参数进行重新计算,重新得到目标灌注参数;
以重新得到的目标灌注参数重新形成所述当前灌注物的目标灌注图像。
5.根据权利要求1所述的方法,其特征在于,还包括:
根据所述目标灌注图像中各像素点的目标灌注参数,在所述目标灌注图像中标识病变区域,并提供所述病变区域的病变信息。
6.根据权利要求1所述的方法,其特征在于,还包括:
响应于对所述目标灌注图像中选定位置上灌注参数的显示请求操作,查找所述选定位置上的目标灌注参数以及各个初始灌注参数并显示。
7.根据权利要求1所述的方法,其特征在于,还包括:
分别以各个所述初始灌注参数形成所述当前灌注物在各种所述灌注算法下的灌注图像作为各个初始灌注图像,并显示各个所述初始灌注图像与所述目标灌注图像。
8.一种CT灌注成像设备,其特征在于,包括:
第一获取模块,用于获取CT设备对当前灌注物扫描得到的灌注数据,所述灌注数据表示时间密度曲线;
第一计算模块,用于利用至少两种灌注算法,分别对所述灌注数据进行计算,得到各种所述灌注算法下的灌注参数作为各个初始灌注参数;
第二计算模块,用于按照各种灌注算法的权重,对各个所述初始灌注参数进行计算,得到目标灌注参数;
第一形成模块,用于以所述目标灌注参数形成所述当前灌注物的目标灌注图像。
9.根据权利要求8所述的设备,其特征在于,所述各种灌注算法的权重为预先采用机器学习方法对历史灌注物在各种所述灌注算法下计算出的灌注参数与所述历史灌注物的真实灌注参数之间对应关系进行机器学习而得到的。
10.根据权利要求8所述的设备,其特征在于,还包括:
第二获取模块,用于响应于对设置系数的输入操作,获取为各种所述灌注算法输入的设置系统;
第三计算模块,用于按照各种所述灌注算法的设置系数,对各个所述初始灌注参数进行重新计算,重新得到目标灌注参数;
第二形成模块,用于以重新得到的目标灌注参数重新形成所述当前灌注物的目标灌注图像。
11.根据权利要求7所述的设备,其特征在于,还包括:
第二查找模块,用于响应于对所述目标灌注图像中选定位置上灌注参数的显示请求操作,查找所述选定位置上的目标灌注参数以及各个初始灌注参数;
第一显示模块,用于显示所述选定位置上的目标灌注参数以及各个所述初始灌注参数。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410462793.9A CN104287764B (zh) | 2014-09-11 | 2014-09-11 | 一种ct灌注成像方法和设备 |
US14/852,386 US10307126B2 (en) | 2014-09-11 | 2015-09-11 | Computed tomography perfusion imaging |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410462793.9A CN104287764B (zh) | 2014-09-11 | 2014-09-11 | 一种ct灌注成像方法和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104287764A true CN104287764A (zh) | 2015-01-21 |
CN104287764B CN104287764B (zh) | 2017-05-31 |
Family
ID=52307905
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410462793.9A Expired - Fee Related CN104287764B (zh) | 2014-09-11 | 2014-09-11 | 一种ct灌注成像方法和设备 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10307126B2 (zh) |
CN (1) | CN104287764B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017148332A1 (zh) * | 2016-02-29 | 2017-09-08 | 上海联影医疗科技有限公司 | 一种灌注分析方法与设备 |
CN113034438A (zh) * | 2021-03-01 | 2021-06-25 | 北京安德医智科技有限公司 | 一种ct灌注参数预测方法及装置、电子设备和存储介质 |
CN113393433A (zh) * | 2021-06-10 | 2021-09-14 | 北京安德医智科技有限公司 | 一种通用的医学影像灌注参数预测方法及装置 |
WO2021184600A1 (zh) * | 2020-03-18 | 2021-09-23 | 上海商汤智能科技有限公司 | 一种图像分割方法及装置、设备及计算机可读存储介质 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678379B (zh) * | 2016-01-12 | 2020-08-07 | 腾讯科技(深圳)有限公司 | 一种cnn的处理方法和装置 |
JP6899666B2 (ja) * | 2017-02-28 | 2021-07-07 | 国立大学法人信州大学 | 造影ctスキャン装置、造影ctスキャン装置の作動方法、及び造影ctスキャン像を得るための方法 |
US10628659B2 (en) | 2017-11-27 | 2020-04-21 | International Business Machines Corporation | Intelligent tumor tracking system |
US20230218169A1 (en) * | 2018-03-30 | 2023-07-13 | A-Moy Limited | Brain imaging system and brain imaging method |
US11747421B2 (en) * | 2018-05-18 | 2023-09-05 | Case Western Reserve University | System and method for quantifying perfusion using a dictionary matching |
CN114121235A (zh) * | 2021-11-22 | 2022-03-01 | 上海商汤智能科技有限公司 | 图像检测方法及相关装置、设备和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1805712A (zh) * | 2003-06-12 | 2006-07-19 | 伯拉考开发股份有限公司 | 通过超声造影成像中的补给曲线拟合进行的血液流动估计 |
CN103186895A (zh) * | 2013-04-15 | 2013-07-03 | 武汉大学 | 一种基于神经网络模型的ct灌注医疗图像智能融合方法 |
US20140114618A1 (en) * | 2012-10-19 | 2014-04-24 | Timothy A. Fonte | Systems and methods for numerically evaluating vasculature |
CN103976750A (zh) * | 2012-11-14 | 2014-08-13 | 美国西门子医疗解决公司 | 用于观察脉管系统以及灌注组织的系统 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4930146A (en) | 1989-07-10 | 1990-05-29 | General Electric Company | X-ray tube current control with constant loop gain |
WO2004026007A2 (en) | 2002-09-10 | 2004-03-25 | Newton Scientific, Inc. | X-ray feedback stabilization of an x-ray tube |
JP2005305024A (ja) | 2004-04-26 | 2005-11-04 | Toshiba Corp | X線コンピュータ断層撮影装置 |
US7522744B2 (en) | 2004-08-31 | 2009-04-21 | University Of Iowa Research Foundation | System and method for adaptive bolus chasing computed tomography (CT) angiography |
JP2008000499A (ja) | 2006-06-26 | 2008-01-10 | Toshiba Corp | X線ct装置、x線ct装置の制御方法、及びx線ct装置の制御プログラム |
JP2009039330A (ja) | 2007-08-09 | 2009-02-26 | Toshiba Corp | X線コンピュータ断層撮影装置および撮影方法 |
US9414798B2 (en) | 2009-03-25 | 2016-08-16 | Siemens Medical Solutions Usa, Inc. | System and method for automatic trigger-ROI detection and monitoring during bolus tracking |
JP5604965B2 (ja) | 2010-05-10 | 2014-10-15 | 株式会社島津製作所 | 放射線透視・撮影装置 |
CN102612248A (zh) | 2011-01-25 | 2012-07-25 | 南京普爱射线影像设备有限公司 | 一种用于x光机灯丝电流和管电流双闭环控制装置 |
JP5978427B2 (ja) | 2011-07-27 | 2016-08-24 | 株式会社日立製作所 | X線ct装置及び造影撮影方法 |
US8964940B2 (en) | 2012-11-21 | 2015-02-24 | Thermo Scientific Portable Analytical Instruments Inc. | Dynamically adjustable filament control through firmware for miniature x-ray source |
-
2014
- 2014-09-11 CN CN201410462793.9A patent/CN104287764B/zh not_active Expired - Fee Related
-
2015
- 2015-09-11 US US14/852,386 patent/US10307126B2/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1805712A (zh) * | 2003-06-12 | 2006-07-19 | 伯拉考开发股份有限公司 | 通过超声造影成像中的补给曲线拟合进行的血液流动估计 |
US20140114618A1 (en) * | 2012-10-19 | 2014-04-24 | Timothy A. Fonte | Systems and methods for numerically evaluating vasculature |
CN103976750A (zh) * | 2012-11-14 | 2014-08-13 | 美国西门子医疗解决公司 | 用于观察脉管系统以及灌注组织的系统 |
CN103186895A (zh) * | 2013-04-15 | 2013-07-03 | 武汉大学 | 一种基于神经网络模型的ct灌注医疗图像智能融合方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017148332A1 (zh) * | 2016-02-29 | 2017-09-08 | 上海联影医疗科技有限公司 | 一种灌注分析方法与设备 |
WO2021184600A1 (zh) * | 2020-03-18 | 2021-09-23 | 上海商汤智能科技有限公司 | 一种图像分割方法及装置、设备及计算机可读存储介质 |
CN113034438A (zh) * | 2021-03-01 | 2021-06-25 | 北京安德医智科技有限公司 | 一种ct灌注参数预测方法及装置、电子设备和存储介质 |
CN113393433A (zh) * | 2021-06-10 | 2021-09-14 | 北京安德医智科技有限公司 | 一种通用的医学影像灌注参数预测方法及装置 |
CN113393433B (zh) * | 2021-06-10 | 2022-03-01 | 北京安德医智科技有限公司 | 一种通用的医学影像灌注参数预测方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104287764B (zh) | 2017-05-31 |
US20160073992A1 (en) | 2016-03-17 |
US10307126B2 (en) | 2019-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104287764A (zh) | 一种ct灌注成像方法和设备 | |
CN110364238B (zh) | 基于机器学习的对比剂管理 | |
Chen et al. | Remodeling of cardiac fiber structure after infarction in rats quantified with diffusion tensor MRI | |
JP5716238B2 (ja) | 急性脳卒中の組織状態をマッピングする方法及びシステム | |
US8571288B2 (en) | Image display apparatus and magnetic resonance imaging apparatus | |
US20170127972A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
EP3281135B1 (en) | System and method for vascular tree generation using patient-specific structural and functional data, and joint prior information | |
CN108122616B (zh) | 个体特异性的心血管模型的生成方法及其应用 | |
JP5976625B2 (ja) | 画像表示装置及び磁気共鳴イメージング装置 | |
CN110391016B (zh) | 一种动态对比增强磁共振图像的分析方法 | |
JP6448918B2 (ja) | 画像処理装置、方法、及び医用画像診断装置 | |
WO2001085230A2 (en) | A method and apparatus for targeting material delivery to tissue | |
WO2006067719A2 (en) | A method, a system and a computer program for integration of medical diagnostic information and a geometric model of a movable body | |
WO2018133098A1 (zh) | 血管壁应力应变状态获取方法及系统 | |
CN101658428A (zh) | 灌注成像技术中对位图进行处理的方法和系统 | |
CN109730677A (zh) | 体素内不相干运动成像的信号处理方法、装置及存储介质 | |
US10438364B2 (en) | System and method for registration of perfusion images | |
US8913812B2 (en) | Medical image processing apparatus | |
CN109498046A (zh) | 基于核素图像与ct冠脉造影融合的心肌梗死定量评估方法 | |
EP1478277B1 (en) | M-tracking for space-time imaging | |
Young et al. | Reperfused myocardial infarction in mice: 3D mapping of late gadolinium enhancement and strain | |
Coon et al. | Quantification of left ventricular size and function using contrast-enhanced real-time 3D imaging with power modulation: comparison with cardiac MRI | |
Belmont et al. | An open-source ultrasound software for diagnosis of fistula maturation | |
Woodall et al. | Model discovery approach enables noninvasive measurement of intra-tumoral fluid transport in dynamic MRI | |
Murta et al. | Segmental quantitative analysis of myocardial contrast echocardiography images using a bullseye representation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Patentee after: Shenyang Neusoft Medical Systems Co.,Ltd. Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16 Patentee before: SHENYANG NEUSOFT MEDICAL SYSTEMS Co.,Ltd. |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170531 |