CN108567422A - 一种血管成像方法和装置 - Google Patents
一种血管成像方法和装置 Download PDFInfo
- Publication number
- CN108567422A CN108567422A CN201810175899.9A CN201810175899A CN108567422A CN 108567422 A CN108567422 A CN 108567422A CN 201810175899 A CN201810175899 A CN 201810175899A CN 108567422 A CN108567422 A CN 108567422A
- Authority
- CN
- China
- Prior art keywords
- echo
- signal
- subgraph
- subspace
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 27
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 22
- 210000003462 vein Anatomy 0.000 claims abstract description 103
- 210000001367 artery Anatomy 0.000 claims abstract description 89
- 230000002792 vascular Effects 0.000 claims abstract description 48
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000012545 processing Methods 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 8
- 238000012935 Averaging Methods 0.000 claims description 6
- 239000008280 blood Substances 0.000 abstract description 25
- 210000004369 blood Anatomy 0.000 abstract description 25
- 210000001519 tissue Anatomy 0.000 description 19
- 238000002595 magnetic resonance imaging Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 10
- 230000005291 magnetic effect Effects 0.000 description 10
- 230000008520 organization Effects 0.000 description 7
- 238000004590 computer program Methods 0.000 description 6
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 229910052739 hydrogen Inorganic materials 0.000 description 5
- 239000001257 hydrogen Substances 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 238000002601 radiography Methods 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 208000029028 brain injury Diseases 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 1
- 208000003174 Brain Neoplasms Diseases 0.000 description 1
- 208000034710 Cerebral arteriovenous malformation Diseases 0.000 description 1
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 208000002263 Intracranial Arteriovenous Malformations Diseases 0.000 description 1
- 238000012307 MRI technique Methods 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 201000000034 arteriovenous malformations of the brain Diseases 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 210000001627 cerebral artery Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 210000003792 cranial nerve Anatomy 0.000 description 1
- 108010002255 deoxyhemoglobin Proteins 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000005298 paramagnetic effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 238000002627 tracheal intubation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- 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/5635—Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- 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/504—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 diagnosis of blood vessels, e.g. by angiography
-
- 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Medical Informatics (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Vascular Medicine (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Optics & Photonics (AREA)
- Dentistry (AREA)
- Hematology (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请提供一种血管成像方法和装置,该方法包括:在第一重复采集周期,获取第一时间区间内的经过流动补偿的第一回波信号,获取第二时间区间内的经过流动补偿的第二回波信号;在第二重复采集周期,获取第三时间区间内的经过流动补偿的第三回波信号,获取第四时间区间内的经过流动散相的第四回波信号;根据第二回波信号获取静脉血管图像,并根据第一回波信号、第三回波信号和第四回波信号获取动静脉血管图像;根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像。通过本申请的技术方案,可以提供动脉血管图像和静脉血管图像,即实现全血管组织的成像,提高用户使用感受。
Description
技术领域
本申请涉及医疗技术领域,尤其是涉及一种血管成像方法和装置。
背景技术
MRI(Magnetic Resonance Imaging,核磁共振成像)是一种利用射频波与原子核***在外磁场中的相互作用,产生可检测信号的成像方法,MRI的本质是能量级间跃迁的量子效应。MRI的工作原理是:将待检体(如患者)置于磁场环境中,使用无线电射频脉冲激发待检体内的氢原子核,引起氢原子核的共振,并吸收能量。在停止射频脉冲后,氢原子核按照特定频率发出射频信号,并将之前吸收的能量释放出来,外部的检测装置接收待检体释放出来的射频信号,将射频信号转换为图像信号,并利用图像信号生成图像。由于MRI没有电离辐射对待检体的损害、采集到的参数多、信息量大、可多方位成像、对软组织有高分辨力等特点,MRI被广泛应用于临床疾病的诊断,是有效的检查方法。
MRA(Magnetic Resonance Angiography,核磁共振血管成像)是MRI技术的典型应用,是一种无创伤性、不需要使用插管、不需要使用对比造影剂的血管成像方法。目前,可以采用TOF(Time Of Flight,时间飞越)或者PC(Phase Contrast,相位对比)技术实现MRA。但是,无论采用哪种技术,均只能提供部分血管的图像信息,而无法同时提供与背景组织之间对比度很高的动脉血管的图像信息和静脉血管的图像信息,因此,用户使用感受比较差。
发明内容
本申请提供一种血管成像方法,应用于医疗设备,所述方法包括:
在第一重复采集周期,获取第一时间区间内的经过流动补偿的第一回波信号,并获取第二时间区间内的经过流动补偿的第二回波信号;
在第二重复采集周期,获取第三时间区间内的经过流动补偿的第三回波信号,并获取第四时间区间内的经过流动散相的第四回波信号;
根据所述第二回波信号获取静脉血管图像,并根据所述第一回波信号、所述第三回波信号和所述第四回波信号获取动静脉血管图像;
根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像。
本申请提供一种血管成像装置,应用于医疗设备,所述装置包括:
第一获取模块,用于在第一重复采集周期获取第一时间区间内的经过流动补偿的第一回波信号,获取第二时间区间内的经过流动补偿的第二回波信号;
第二获取模块,用于在第二重复采集周期获取第三时间区间内的经过流动补偿的第三回波信号,获取第四时间区间内的经过流动散相的第四回波信号;
第三获取模块,用于根据所述第二回波信号获取静脉血管图像,并根据所述第一回波信号、所述第三回波信号和所述第四回波信号获取动静脉血管图像,并根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像。
基于上述技术方案,本申请实施例中,在第一重复采集周期,获取第一时间区间内的经过流动补偿的第一回波信号,并获取第二时间区间内的经过流动补偿的第二回波信号,在第二重复采集周期,获取第三时间区间内的经过流动补偿的第三回波信号,并获取第四时间区间内的经过流动散相的第四回波信号;这样,可以根据第二回波信号获取静脉血管图像,并可以根据第一回波信号、第三回波信号和第四回波信号获取动静脉血管图像,并根据动静脉血管图像和静脉血管图像得到动脉血管图像。上述方式可以提供动脉血管图像(动脉与背景组织的高对比度图)和静脉血管图像(静脉与背景组织的高对比度图),即实现高对比度全血管组织的成像,提高用户使用感受,提高应用价值。
附图说明
为了更加清楚地说明本申请实施例或者现有技术中的技术方案,下面将对本申请实施例或者现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,还可以根据本申请实施例的这些附图获得其他的附图。
图1是本申请一种实施方式中的MRI设备的一个剖面示意图;
图2是本申请一种实施方式中的空间编码方式的示意图;
图3是本申请一种实施方式中的血管成像方法的流程图;
图4A和图4B是本申请一种实施方式中的K空间的划分示意图;
图5是本申请一种实施方式中的医疗设备的硬件结构图;
图6是本申请一种实施方式中的血管成像装置的结构图。
具体实施方式
在本申请使用的术语仅仅是出于描述特定实施例的目的,而非限制本申请。本申请和权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其它含义。还应当理解,本文中使用的术语“和/或”是指包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解,尽管在本申请可能采用术语第一、第二、第三等来描述各种信息,但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如,在不脱离本申请范围的情况下,第一信息也可以被称为第二信息,类似地,第二信息也可以被称为第一信息。取决于语境,此外,所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
本申请实施例中提出一种血管成像方法,该方法可以应用于医疗设备,例如,MRI设备、CT(Computed Tomography,电子计算机断层扫描)设备、PET(Positron EmissionComputed Tomography,正电子发射型计算机断层显像)设备、DR(Digital Radiography,数字放射显影***)设备等,后续以MRI设备为例,可以理解的是,该方法也可以应用于其它医疗设备,对此不做限制。
参见图1所示,为MRI设备的一个剖面示意图,MRI设备可以包括:主磁场150、梯度线圈110、射频发射线圈120和射频接收线圈130,当然,MRI设备还可以包括其它器件,对此不做限制,以图1所示的结构为例进行说明。
其中,主磁场150提供了成像所需的磁场环境,在将待检体置于扫描床140后,待检体对应的扫描部位需要位于主磁场150提供的磁场环境内。梯度线圈110、射频发射线圈120和射频接收线圈130位于主磁场150提供的磁场环境内。
其中,射频发射线圈120可以在指定扫描位置(如待检体的头部、心脏等位置)发射脉冲信号(用于激发待检体内的氢原子核,并引起氢原子核的共振),待检体在接收到脉冲信号后,按照特定频率发出射频信号,该射频信号是上述脉冲信号的共振信号,且待检体产生的射频信号被射频接收线圈130接收到。
射频接收线圈130接收到射频信号后,将射频信号传输给频谱仪,频谱仪对射频信号进行分析,将射频信号转换为图像信号,将图像信号传输给计算机。计算机利用图像信号生成图像,将图像提供给医护人员,由医护人员进行诊断。
在射频发射线圈120发射脉冲信号时,梯度线圈110还可以向待检体提供空间编码信息,这样,待检体受激发的部位根据空间编码信息在脉冲信号激励作用下产生射频信号,即该射频信号可以携带空间编码信息。射频接收线圈130接收到的射频信号包括空间编码信息,并最终在重建图像中恢复出不同位置处的组织信息,即重建计算机可以利用包括空间编码信息的图像信号生成图像。
本申请实施例中,为了得到静脉血管图像(如脑静脉血管图像等)和动脉血管图像(如脑动脉血管图像等),还可以采用图2所示的空间编码方式。当然,图2所示的空间编码方式只是一个示例,后续以图2的空间编码方式为例。
参见图2所示,可以将所有的采集周期划分为第一重复采集周期和第二重复采集周期。例如,第1个采集周期进行第一重复采集周期,第2个采集周期进行第二重复采集周期,第3个采集周期在修改部分编码梯度后重复第一重复采集周期,第4个采集周期在修改部分编码梯度后重复第二重复采集周期,以此类推。进一步的,第一重复采集周期的时长与第二重复采集周期的时长保持相同。例如,第一重复采集周期的时长和第二重复采集周期的时长均为50毫秒。
参见图2所示,第一重复采集周期可以包括但不限于:第一时间区间、第二时间区间、时间区间A和时间区间B等。此外,第二重复采集周期可以包括但不限于:第三时间区间、第四时间区间、时间区间C和时间区间D等。
在时间区间A,射频发射线圈120可以发射脉冲信号,且梯度线圈110可以向待检体提供空间编码信息;在时间区间A,射频接收线圈130可以不进行工作;在第一时间区间,射频发射线圈120可以不进行工作,且待检体可以按照特定频率发出射频信号(即上述脉冲信号的共振信号),且射频接收线圈130可以接收到待检体产生的射频信号,这个射频信号包括空间编码信息。
在时间区间B,射频发射线圈120发射脉冲信号,梯度线圈110向待检体提供空间编码信息;在时间区间B,射频接收线圈130不进行工作;在第二时间区间,射频发射线圈120不进行工作,待检体按照特定频率发出射频信号,且射频接收线圈130接收待检体产生的射频信号,这个射频信号包括空间编码信息。
在时间区间C,射频发射线圈120发射脉冲信号,梯度线圈110向待检体提供空间编码信息;在时间区间C,射频接收线圈130不进行工作;在第三时间区间,射频发射线圈120不进行工作,待检体按照特定频率发出射频信号,且射频接收线圈130接收待检体产生的射频信号,这个射频信号包括空间编码信息。
在时间区间D,射频发射线圈120发射脉冲信号,梯度线圈110向待检体提供空间编码信息;在时间区间D,射频接收线圈130不进行工作;在第四时间区间,射频发射线圈120不进行工作,待检体按照特定频率发出射频信号,射频接收线圈130接收待检体产生的射频信号,该射频信号包括空间编码信息。
在上述实施例中,射频接收线圈130接收到的射频信号可以是回波信号。
参见图2所示,在时间区间A、时间区间B和时间区间C,梯度线圈110在向待检体提供空间编码信息时,该空间编码信息用于对待检体的流动组织(如血液等)进行流动补偿,从而使得待检体产生的回波信号是经过流动补偿的回波信号。在时间区间D,梯度线圈110在向待检体提供空间编码信息时,该空间编码信息用于对待检体的流动组织进行流动散相(也可以称为相位发散、主动相位发散等),从而使得待检体产生的回波信号是经过流动散相的回波信号。
其中,梯度线圈110可以通过选层梯度、编码梯度、读出梯度的不同形状,向待检体提供用于进行流动补偿或者流动散相的空间编码信息。例如,通过时间区间A内的选层梯度、编码梯度、读出梯度,提供用于进行流动补偿的空间编码信息;通过时间区间B内的选层梯度、编码梯度、读出梯度,提供用于进行流动补偿的空间编码信息;时间区间C内的各梯度与时间区间A相同,不再赘述;通过时间区间D内的选层梯度、编码梯度、读出梯度,提供用于进行流动散相的空间编码信息。当然,图2所示的各选层梯度、编码梯度、读出梯度,只是用于实现流动补偿功能或者流动散相功能的示例,对此不做限制。
综上所述,在时间区间A、时间区间B和时间区间C,通过向待检体提供图2所示的选层梯度、编码梯度、读出梯度,从而可以对待检体内的流动组织进行流动补偿,这样,使得待检体可以发出经过流动补偿的回波信号。与此类似的,在时间区间D,待检体可以发出经过流动散相的回波信号。
其中,以流动组织是血液为例,对待检体内的血液进行流动补偿是指:由于血液具有流动性,而血液的流动会导致成像像素内信号变弱,且考虑到血液可近似为匀速流动,因此,可以采用流动补偿保持匀速流动的血液信号,从而可以避免血液信号变弱,这样,经过流动补偿的回波信号可以包括所有血液信号。此外,对待检体内的血液进行流动散相是指:考虑到血液是匀速流动,因此,可以采用对应一定流速的编码梯度发散掉对应流速的血液信号,从而衰减血液信号,这样,经过流动散相的回波信号包括信号衰减的血液信号。
在上述实施例中,通过使用选层梯度、编码梯度、读出梯度,可以实现对待检体内的流动组织的流动补偿或者流动散相,当然,选层梯度、编码梯度、读出梯度还具有其它功能,对此不做限制。此外,图2所示的选层梯度、编码梯度、读出梯度,只是实现流动补偿或者流动散相的示例,对此不做限制。
在上述实施例中,第一时间区间的回波时间(即回波信号的采集时长)与第三时间区间的回波时间保持相同,第二时间区间的回波时间与第四时间区间的回波时间保持相同,第二时间区间的回波时间大于第一时间区间的回波时间,第四时间区间的回波时间大于第三时间区间的回波时间。
例如,第一时间区间的回波时间可以较短,如10毫秒等,第三时间区间的回波时间与第一时间区间的回波时间相同。第二时间区间的回波时间可以较长,如40毫秒等,第四时间区间的回波时间与第二时间区间的回波时间相同。
在上述应用场景下,参见图3所示,为本申请实施例中提出的血管成像方法的流程示意图,该方法可以应用于医疗设备(如MRI设备),该方法包括:
步骤301,在第一重复采集周期,获取第一时间区间内的经过流动补偿的第一回波信号,并获取第二时间区间内的经过流动补偿的第二回波信号。
其中,获取第一时间区间内的经过流动补偿的第一回波信号,包括:获取第一时间区间对应的K空间内的经过流动补偿的第一回波信号,第一回波信号可以由多个一维回波信号组成,不同一维回波信号对应不同的空间相位编码。
参见图2所示,在时间区间A,射频发射线圈120发射脉冲信号,梯度线圈110向待检体提供空间编码信息,该空间编码信息用于对待检体内的流动组织进行流动补偿。在第一时间区间,射频接收线圈130接收待检体产生的经过流动补偿的回波信号A,该回波信号A是一维回波信号。由于第一时间区间对应的K空间,可以包括大量一维回波信号,因此,通过不断调整空间相位编码,可以将回波信号A转换成其它一维回波信号,每次调整空间相位编码时就可以得到新的回波信号,这样将得到的所有一维回波信号组成第一回波信号。
参见图4A所示,假设第一时间区间对应的K空间包括101个一维回波信号,在得到回波信号A后,可以调整100次空间相位编码,对此调整过程不做限制,最终得到101个一维回波信号,这些一维回波信号组成第一回波信号。
在图4A中,每个直线表示一个一维回波信号,共有101个直线,当然,这种表示方式只是一个示例,实际采集时不限于采集101个回波信号。
其中,获取第二时间区间内的经过流动补偿的第二回波信号,包括:获取第二时间区间对应的K空间内的经过流动补偿的第二回波信号,第二回波信号可以由多个一维回波信号组成,不同一维回波信号对应不同的空间相位编码。
第二回波信号的获取方式与第一回波信号的获取方式类似,不再赘述。
步骤302,在第二重复采集周期,获取第三时间区间内的经过流动补偿的第三回波信号,并获取第四时间区间内的经过流动散相的第四回波信号。
其中,获取第三时间区间内的经过流动补偿的第三回波信号,可以包括:
方式一、获取第三时间区间对应的K空间内的经过流动补偿的第三回波信号,该第三回波信号由多个一维回波信号组成,不同一维回波信号对应不同的空间相位编码,获取方式与第一回波信号的获取方式类似,在此不再赘述。
方式二、将第三时间区间对应的K空间划分成第一子空间和第二子空间,第一子空间和第二子空间组成第三时间区间对应的K空间。获取第一子空间内的经过流动补偿的回波信号,但不采集第二子空间的回波信号,而是从第一回波信号中选择与第二子空间对应的回波信号,并将第一子空间的回波信号与第二子空间的回波信号结合组成第三回波信号。
其中,在将第三时间区间对应的K空间划分成第一子空间和第二子空间时,第一子空间的相位编码索引与K空间中间的距离小于第一阈值(可以根据经验进行配置,如20等),第二子空间的相位编码索引与K空间中心的距离大于等于第一阈值,且第一子空间和第二子空间能够完全等价于所述K空间。
例如,假设第三时间区间对应的K空间包括101个一维回波信号,第一阈值为20,则K空间中心、第一子空间和第二子空间的示意图参见图4B所示。
K空间实际上是频率空间,是图像空间在傅利叶转换下的共轭空间,应用在磁共振造影的成像分析。上述图4A和图4B的K空间、K空间中心只是为了方便描述的示例,实际应用中,K空间的维度大小和K空间中心是已知的。
在已知K空间的维度大小后,参见图4A所示,可以通过调整空间相位编码,就可以得到K空间的所有回波信号。在已知K空间的维度大小和K空间中心后,参见图4B所示,可以先确定第一子空间,第一子空间位于K空间中心附近,例如,K空间中心的相位编码索引是51,第一子空间的相位编码索引与K空间中心的距离小于第一阈值,因此,第一子空间的相位编码索引是31-50、以及52-71。然后,可以通过调整空间相位编码,得到K空间的第一子空间的所有回波信号,即相位编码索引是31-50、52-71的所有回波信号。
其中,针对获取第一子空间内的经过流动补偿的回波信号的方式,与方式一类似,不同之处在于:在方式一中,是获取K空间的所有回波信号,而方式二中,是获取该K空间的第一子空间的所有回波信号,即方式二中的回波信号的获取数量小于方式一中的回波信号的获取数量,对此获取过程不再赘述。
其中,从所述第一回波信号中选择与第二子空间对应的回波信号,可以包括:将第一时间区间对应的K空间划分成第一子空间和第二子空间,其划分方式与第三时间区间对应的K空间的划分方式相同。然后,从第一回波信号中选择与第二子空间对应的回波信号,也就是说,第一回波信号中的第二子空间对应的回波信号,可以作为第三回波信号中的第二子空间对应的回波信号。
例如,参见图4B所示,可以将第一回波信号中的第1个回波信号(即第1条直线)作为第三回波信号中的第1个回波信号,以此类推,将第一回波信号中的第30个回波信号作为第三回波信号中的第30个回波信号,将第一回波信号中的第72个回波信号作为第三回波信号中的第72个回波信号,以此类推。
经过上述处理,可以得到第一子空间的回波信号和第二子空间的回波信号,并将第一子空间的回波信号与第二子空间的回波信号组成第三回波信号。
显然,在采用方式二时,由于第一回波信号的获取时间与第三回波信号的获取时间比较接近,在相隔较近的时间内,待检体内的流动组织信号不会发生较大变化,而第一回波信号和第三回波信号对应的选层梯度、编码梯度、读出梯度相同,因此,第一回波信号和第三回波信号可以相同或近似。基于此,针对第三回波信号中的第二子空间的回波信号,可以直接复用第一回波信号中的第二子空间的回波信号,不用采用调整空间相位编码的方式来获取第三回波信号中的第二子空间的回波信号,这样,大幅减少第三回波信号的获取时间。
其中,获取第四时间区间内的经过流动散相的第四回波信号,可以包括:
方式A、获取第四时间区间对应的K空间内的经过流动散相的第四回波信号,该第四回波信号由多个一维回波信号组成,不同一维回波信号对应不同的空间相位编码,获取方式与第一回波信号的获取方式类似,在此不再赘述。
方式B、将第四时间区间对应的K空间划分成第三子空间和第四子空间,第三子空间和第四子空间组成第四时间区间对应的K空间。获取第三子空间内的经过流动散相的回波信号,并不采集第四子空间的回波信号,而是采用直接零值填充的策略获取第四子空间的回波信号,即第四子空间的回波信号为0,并将第三子空间的回波信号与第四子空间的零值信号结合组成第四回波信号。
其中,在将第四时间区间对应的K空间划分成第三子空间和第四子空间时,第三子空间的相位编码索引与K空间中间的距离小于第二阈值(可以根据经验进行配置,如20等),第四子空间的相位编码索引与K空间中心的距离大于等于第二阈值,且第三子空间和第四子空间能够完全等价于所述K空间。
其中,针对获取第三子空间内的经过流动散相的回波信号的方式,与方式A类似,不同之处在于:在方式A中,是获取K空间的所有回波信号,而方式B中,是获取K空间的第三子空间的所有回波信号,对此获取过程不再赘述。
其中,采用零值填充的策略获取第四子空间的回波信号,可以包括:将第四子空间的所有回波信号均填充为0,即第四子空间的所有回波信号均为0。
经过上述处理,可以得到第三子空间的回波信号和第四子空间的回波信号,并将第三子空间的回波信号与第四子空间的回波信号组成第四回波信号。
显然,在采用方式B时,针对第四回波信号中的第四子空间的回波信号,可以使用0值填充,不用采用调整空间相位编码的方式来获取第四回波信号中的第四子空间的回波信号,这样,可以大幅减少第四回波信号的获取时间。
由于已经获取第四回波信号中的第三子空间的回波信号(即K空间中心的信息),因此,主要的信息均在第四回波信号的实际采集的编码线中,将第四回波信号中的第四子空间的回波信号填充为0,不会造成主要信息的丢失,即第四回波信号是准确的。
步骤303,根据第二回波信号获取静脉血管图像,即静脉血管的图像。
其中,由于第二回波信号是经过流动补偿的第二回波信号,且第二回波信号对应的回波时间较长,因此,第二回波信号可以凸显磁敏感效应较强的静脉组织血管。在此基础上,采用SWI(Susceptibility Weighted Imaging,磁敏感加权成像)技术,就可以获取到静脉血管图像,对此获取过程不做限制。
其中,在上述静脉血管图像中,可以凸显出富含顺磁性的去氧血红蛋白的静脉血管,可以用于脑中铁质沉积、急性脑损伤等疾病的对比度增强显示。
步骤304,根据第一回波信号、第三回波信号和第四回波信号获取动静脉血管图像。其中,所述动静脉血管图像是包括动脉血管和静脉血管的图像,并且在动静脉血管图像中,可以通过重建处理将脑实质背景信号基本去除。
其中,根据第一回波信号、第三回波信号和第四回波信号获取动静脉血管图像,包括:根据第一回波信号获取第一子图像,根据第三回波信号获取第二子图像,并根据第一子图像和/或第二子图像获取第三子图像;根据第四回波信号获取第四子图像;根据第三子图像和第四子图像获取动静脉血管图像。
进一步的,根据第一子图像和/或第二子图像获取第三子图像,可以包括但不限于:将第一子图像确定为第三子图像;或者,将第二子图像确定为第三子图像;或者,对第一子图像与第二子图像进行平均处理,得到第三子图像。其中,若对第一子图像与第二子图像进行平均处理,即对两个子图像进行平均处理(累加处理),得到第三子图像,则可以提高图像的信噪比,增强血管信号。
此外,根据第三子图像和第四子图像获取动静脉血管图像,可以包括但不限于:对第三子图像与第四子图像进行相减处理,得到动静脉血管图像。
在一个例子中,由于第一回波信号是经过流动补偿的信号,且流动补偿具有回聚相位的功能,从而保持动脉血管和静脉血管的信息,因此,第一回波信号包括动脉血管、静脉血管和背景组织(动脉血管和静脉血管之外的其它器官)的信息,同理,第三回波信号也可以包括动脉血管、静脉血管和背景组织的信息。综上所述,第三子图像可以包括动脉血管、静脉血管和背景组织的信息。
由于第四回波信号是经过流动散相的信号,且流动散相具有发散相位的功能,从而可以得到对应流速血液信号降低的动脉血管和静脉血管的信息,因此,第四回波信号可以包括血液信号降低但背景组织基本不受影响的信息。综上所述,第四子图像可以包括背景组织不受影响信息,但动脉血管和静脉血管信号降低的图像信息。基于此,在对第三子图像与第四子图像进行相减处理时,就可以得到没有背景组织的动静脉血管图像,即动静脉血管图像包括动脉血管和静脉血管的信息,从而使得背景组织得到抑制,消除动静脉血管图像中的背景组织,凸显血管组织信号。
具体的,假设第三子图像为SRP,SRP=SA+SB,SA是动脉血管和静脉血管的信息,而SB是背景组织的信息;假设第四子图像为SDP,SDP=SA’+SB’,SA’是动脉血管和静脉血管的信息,而SB’是背景组织的信息。由于进行流动补偿或者流动散相时,对背景组织基本没有影响,因此,SB和SB’近似相同,即二者相差不大,后续以SB和SB’相同为例。由于进行流动补偿时,能够保持动脉血管和静脉血管的信息,而进行流动散相时,能够降低对应流速的血液信号,因此SA远远大于SA’,SA’可近似为0,SA是实际的动脉血管和静脉血管的信息。
综上所述,在对第三子图像与第四子图像进行相减处理时,得到的血管信号是SA和SA’的差异值,即没有背景组织的动静脉血管图像,增强了动静脉血管与背景组织的对比度。
步骤305,根据静脉血管图像和动静脉血管图像获取动脉血管图像。如可以对动静脉血管图像与静脉血管图像进行相减处理,得到动脉血管图像。
其中,由于上述静脉血管图像是包括静脉血管的图像,而上述动静脉血管图像是包括动脉血管和静脉血管的图像,因此,在对动静脉血管图像与静脉血管图像进行相减处理后,得到的就是包括动脉血管的动脉血管图像。
在上述实施例中,针对根据回波信号获取图像的过程,射频接收线圈130可以将回波信号传输给频谱仪,频谱仪可以将回波信号转换为图像信号,将图像信号传输给计算机,由计算机利用图像信号生成图像。例如,射频接收线圈130可以将第二回波信号传输给频谱仪,频谱仪将第二回波信号转换为图像信号,将图像信号传输给计算机,由计算机利用该图像信号获取静脉血管图像。
基于上述技术方案,本申请实施例中,在第一重复采集周期,获取第一时间区间内的经过流动补偿的第一回波信号,并获取第二时间区间内的经过流动补偿的第二回波信号,在第二重复采集周期,获取第三时间区间内的经过流动补偿的第三回波信号,并获取第四时间区间内的经过流动散相的第四回波信号;这样,可以根据第二回波信号获取静脉血管图像,并可以根据第一回波信号、第三回波信号和第四回波信号获取动静脉血管图像,并根据动静脉血管图像和静脉血管图像得到动脉血管图像。上述方式可以提供动脉血管图像(动脉与背景组织的高对比度图)和静脉血管图像(静脉与背景组织的高对比度图),实现高对比度全血管组织的成像,提高用户使用感受,提高应用价值,在脑神经科学研究及临床应用方面具有重要价值。例如,在脑动静脉畸形、中风、急性脑损伤和肿瘤等诊断中,上述血管成像技术对疾病的辅助诊断非常重要。
如图5所示,为本申请实施例中的医疗设备的硬件结构图,该医疗设备可以包括扫描***510和控制台***520。该扫描***510可以包括线圈***511和扫描床512,该线圈***511是医疗设备的重要组成部分,线圈***511包含梯度线圈、射频发射线圈和射频接收线圈,当然,还可以有其它器件,如处理器5111和存储器5112等,对此不做限制。该扫描床512是配合线圈***511完成扫描任务的工具,用于支撑被检体。该控制台***520可以包括计算机521和输入面板等,对此不做限制。存储器5112中的血管成像装置作为一个逻辑意义上的装置,可以由处理器5111从存储器5112中将该装置对应的计算机程序指令读取到内存中运行,以实现本申请上述实施例中的血管成像方法。
处理器5111通过读取存储器5112中对应指令执行以下操作:在第一重复采集周期,获取第一时间区间内的经过流动补偿的第一回波信号,并获取第二时间区间内的经过流动补偿的第二回波信号;在第二重复采集周期,获取第三时间区间内的经过流动补偿的第三回波信号,并获取第四时间区间内的经过流动散相的第四回波信号;根据所述第二回波信号获取静脉血管图像,并根据所述第一回波信号、所述第三回波信号和所述第四回波信号获取动静脉血管图像;根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像。
参见图6所示,为本申请实施例的血管成像装置的结构图,所述装置包括:
第一获取模块611,用于在第一重复采集周期获取第一时间区间内的经过流动补偿的第一回波信号、第二时间区间内的经过流动补偿的第二回波信号;
第二获取模块612,用于在第二重复采集周期获取第三时间区间内的经过流动补偿的第三回波信号、第四时间区间内的经过流动散相的第四回波信号;
第三获取模块613,用于根据所述第二回波信号获取静脉血管图像,并根据所述第一回波信号、所述第三回波信号和所述第四回波信号获取动静脉血管图像,并根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像。
在一个例子中,所述第二获取模块612在获取第三时间区间内的经过流动补偿的第三回波信号时具体用于:将第三时间区间对应的K空间划分成第一子空间和第二子空间;获取所述第一子空间内的经过流动补偿的回波信号;从所述第一回波信号中选择与所述第二子空间对应的回波信号;将第一子空间的回波信号与第二子空间的回波信号组成第三回波信号;所述第二获取模块612在获取第四时间区间内的经过流动散相的第四回波信号时具体用于:将第四时间区间对应的K空间划分成第三子空间和第四子空间;获取所述第三子空间内的经过流动散相的回波信号;采用零值填充的策略获取第四子空间的回波信号;将第三子空间的回波信号与第四子空间的回波信号组成第四回波信号。
在一个例子中,所述第三获取模块613在根据所述第一回波信号、所述第三回波信号和所述第四回波信号获取动静脉血管图像时具体用于:根据所述第一回波信号获取第一子图像,根据所述第三回波信号获取第二子图像,根据所述第一子图像和/或所述第二子图像获取第三子图像;根据所述第四回波信号获取第四子图像;根据所述第三子图像和所述第四子图像获取动静脉血管图像。
在一个例子中,所述第三获取模块613在根据所述第一子图像和/或所述第二子图像获取第三子图像时具体用于:将第一子图像确定为所述第三子图像;或者,将第二子图像确定为所述第三子图像;或者,对所述第一子图像与所述第二子图像进行平均处理,得到所述第三子图像;所述第三获取模块613在根据所述第三子图像和所述第四子图像获取动静脉血管图像时具体用于:对所述第三子图像与所述第四子图像进行相减处理,得到动静脉血管图像。
在一个例子中,所述第三获取模块613在根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像时具体用于:对所述动静脉血管图像与静脉血管图像进行相减处理,得到动脉血管图像。
上述实施例阐明的***、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机,计算机的具体形式可以是个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件收发设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任意几种设备的组合。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可以由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其它可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其它可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
而且,这些计算机程序指令也可以存储在能引导计算机或其它可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或者多个流程和/或方框图一个方框或者多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其它可编程数据处理设备上,使得在计算机或者其它可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其它可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (15)
1.一种血管成像方法,其特征在于,应用于医疗设备,所述方法包括:
在第一重复采集周期,获取第一时间区间内的经过流动补偿的第一回波信号,并获取第二时间区间内的经过流动补偿的第二回波信号;
在第二重复采集周期,获取第三时间区间内的经过流动补偿的第三回波信号,并获取第四时间区间内的经过流动散相的第四回波信号;
根据所述第二回波信号获取静脉血管图像,并根据所述第一回波信号、所述第三回波信号和所述第四回波信号获取动静脉血管图像;
根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像。
2.根据权利要求1所述的方法,其特征在于,
所述第一时间区间的回波时间与所述第三时间区间的回波时间相同;
所述第二时间区间的回波时间与所述第四时间区间的回波时间相同;
所述第二时间区间的回波时间大于所述第一时间区间的回波时间;
所述第四时间区间的回波时间大于所述第三时间区间的回波时间。
3.根据权利要求1所述的方法,其特征在于,
所述获取第三时间区间内的经过流动补偿的第三回波信号,包括:
将所述第三时间区间对应的K空间划分成第一子空间和第二子空间;
获取所述第一子空间内的经过流动补偿的回波信号;
从所述第一回波信号中选择与所述第二子空间对应的回波信号;
将第一子空间的回波信号与第二子空间的回波信号组成第三回波信号。
4.根据权利要求3所述的方法,其特征在于,
从所述第一回波信号中选择与所述第二子空间对应的回波信号,包括:
将所述第一时间区间对应的K空间划分成第一子空间和第二子空间;
从所述第一回波信号中选择与所述第二子空间对应的回波信号。
5.根据权利要求1所述的方法,其特征在于,
所述获取第四时间区间内的经过流动散相的第四回波信号,包括:
将所述第四时间区间对应的K空间划分成第三子空间和第四子空间;
获取所述第三子空间内的经过流动散相的回波信号;
采用零值填充的策略获取所述第四子空间的回波信号;
将第三子空间的回波信号与第四子空间的回波信号组成第四回波信号。
6.根据权利要求3-5任一项所述的方法,其特征在于,
所述第一子空间的相位编码索引与K空间中间的距离小于第一阈值;
所述第二子空间的相位编码索引与K空间中心的距离大于等于第一阈值;
所述第三子空间的相位编码索引与K空间中间的距离小于第二阈值;
所述第四子空间的相位编码索引与K空间中心的距离大于等于第二阈值。
7.根据权利要求1所述的方法,其特征在于,所述根据所述第一回波信号、所述第三回波信号和所述第四回波信号获取动静脉血管图像,包括:
根据所述第一回波信号获取第一子图像,根据所述第三回波信号获取第二子图像,并根据所述第一子图像和/或所述第二子图像获取第三子图像;
根据所述第四回波信号获取第四子图像;
根据所述第三子图像和所述第四子图像获取动静脉血管图像。
8.根据权利要求7所述的方法,其特征在于,
所述根据所述第一子图像和/或所述第二子图像获取第三子图像,包括:
将所述第一子图像确定为所述第三子图像;或者,
将所述第二子图像确定为所述第三子图像;或者,
对所述第一子图像与所述第二子图像进行平均处理,得到第三子图像。
9.根据权利要求7所述的方法,其特征在于,
所述根据所述第三子图像和所述第四子图像获取动静脉血管图像,包括:对所述第三子图像与所述第四子图像进行相减处理,得到动静脉血管图像。
10.根据权利要求1所述的方法,其特征在于,所述根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像,包括:
对动静脉血管图像与静脉血管图像进行相减处理,得到动脉血管图像。
11.一种血管成像装置,其特征在于,应用于医疗设备,所述装置包括:
第一获取模块,用于在第一重复采集周期获取第一时间区间内的经过流动补偿的第一回波信号,获取第二时间区间内的经过流动补偿的第二回波信号;
第二获取模块,用于在第二重复采集周期获取第三时间区间内的经过流动补偿的第三回波信号,获取第四时间区间内的经过流动散相的第四回波信号;
第三获取模块,用于根据所述第二回波信号获取静脉血管图像,并根据所述第一回波信号、所述第三回波信号和所述第四回波信号获取动静脉血管图像,并根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像。
12.根据权利要求11所述的装置,其特征在于,所述第二获取模块获取第三时间区间内的经过流动补偿的第三回波信号时具体用于:将第三时间区间对应的K空间划分成第一子空间和第二子空间;获取所述第一子空间内的经过流动补偿的回波信号;从所述第一回波信号中选择与所述第二子空间对应的回波信号;将第一子空间的回波信号与第二子空间的回波信号组成第三回波信号;
在获取第四时间区间内的经过流动散相的第四回波信号时具体用于:将第四时间区间对应的K空间划分成第三子空间和第四子空间;获取所述第三子空间内的经过流动散相的回波信号;采用零值填充的策略获取第四子空间的回波信号;将第三子空间的回波信号与第四子空间的回波信号组成第四回波信号。
13.根据权利要求11所述的装置,其特征在于,
所述第三获取模块在根据所述第一回波信号、所述第三回波信号和所述第四回波信号获取动静脉血管图像时具体用于:根据所述第一回波信号获取第一子图像,根据所述第三回波信号获取第二子图像,并根据所述第一子图像和/或所述第二子图像获取第三子图像;根据所述第四回波信号获取第四子图像;根据所述第三子图像和所述第四子图像获取动静脉血管图像。
14.根据权利要求13所述的装置,其特征在于,所述第三获取模块在根据所述第一子图像和/或所述第二子图像获取第三子图像时具体用于:将第一子图像确定为所述第三子图像;或者,将第二子图像确定为所述第三子图像;或者,对所述第一子图像与所述第二子图像进行平均处理,得到第三子图像;
在根据所述第三子图像和所述第四子图像获取动静脉血管图像时具体用于:对所述第三子图像与第四子图像进行相减处理,得到动静脉血管图像。
15.根据权利要求11所述的装置,其特征在于,所述第三获取模块在根据所述静脉血管图像和所述动静脉血管图像获取动脉血管图像时具体用于:对所述动静脉血管图像与静脉血管图像进行相减处理,得到动脉血管图像。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810175899.9A CN108567422B (zh) | 2018-03-02 | 2018-03-02 | 一种血管成像方法和装置 |
US16/291,749 US20190271752A1 (en) | 2018-03-02 | 2019-03-04 | Magnetic resonance angiography |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810175899.9A CN108567422B (zh) | 2018-03-02 | 2018-03-02 | 一种血管成像方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108567422A true CN108567422A (zh) | 2018-09-25 |
CN108567422B CN108567422B (zh) | 2021-03-05 |
Family
ID=63576740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810175899.9A Active CN108567422B (zh) | 2018-03-02 | 2018-03-02 | 一种血管成像方法和装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20190271752A1 (zh) |
CN (1) | CN108567422B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110148134A (zh) * | 2019-04-02 | 2019-08-20 | 佛山瑞加图医疗科技有限公司 | 血管成像方法及装置 |
CN113805130A (zh) * | 2020-06-17 | 2021-12-17 | 西门子(深圳)磁共振有限公司 | 快速磁化率敏感性成像方法、装置及磁共振成像*** |
CN113970717A (zh) * | 2020-07-24 | 2022-01-25 | 西门子(深圳)磁共振有限公司 | 磁敏感加权成像方法、装置及磁共振成像*** |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08266503A (ja) * | 1995-03-31 | 1996-10-15 | Shimadzu Corp | Mrイメージング装置 |
CN102238909A (zh) * | 2008-12-04 | 2011-11-09 | 株式会社日立医疗器械 | 磁共振成像装置以及同步摄像方法 |
US8060180B2 (en) * | 2005-01-14 | 2011-11-15 | New York University | Multi-echo magnetic resonance imaging method and system |
CN103110420A (zh) * | 2013-01-31 | 2013-05-22 | 深圳先进技术研究院 | 磁共振血管成像方法和*** |
CN103371820A (zh) * | 2012-04-20 | 2013-10-30 | 西门子公司 | 借助磁共振技术采集呼吸的检查对象的测量数据组的方法 |
CN103389481A (zh) * | 2012-05-11 | 2013-11-13 | 上海联影医疗科技有限公司 | 磁共振频率相位双编码采样方法及其图像重建方法 |
CN103961097A (zh) * | 2013-02-04 | 2014-08-06 | 上海联影医疗科技有限公司 | 一种磁共振扫描短te成像方法及磁共振扫描*** |
CN104545918A (zh) * | 2014-12-29 | 2015-04-29 | 中国科学院深圳先进技术研究院 | 非对比增强磁共振静脉造影成像方法 |
CN104635189A (zh) * | 2013-11-13 | 2015-05-20 | 上海联影医疗科技有限公司 | 磁共振多对比度图像重建方法及*** |
CN104688228A (zh) * | 2015-03-20 | 2015-06-10 | 中国科学院深圳先进技术研究院 | 一种磁共振成像方法及设备 |
CN105259525A (zh) * | 2015-10-28 | 2016-01-20 | 北京大学 | 基于邻域共享压缩感知的动态对比增强磁共振快速成像方法 |
CN106597333A (zh) * | 2016-12-30 | 2017-04-26 | 上海联影医疗科技有限公司 | 一种磁共振并行成像方法及磁共振成像*** |
CN106646302A (zh) * | 2015-10-29 | 2017-05-10 | 上海联影医疗科技有限公司 | 一种磁共振成像方法 |
CN106821381A (zh) * | 2017-01-06 | 2017-06-13 | 上海交通大学 | 一种血管成像方法 |
US20170242087A1 (en) * | 2016-02-19 | 2017-08-24 | Samsung Electronics Co., Ltd. | Magnetic resonance imaging apparatus and method of obtaining magnetic resonance image by using multiple excitation with delayed spin-echoes |
CN107126211A (zh) * | 2017-06-07 | 2017-09-05 | 上海联影医疗科技有限公司 | 磁共振成像方法和*** |
-
2018
- 2018-03-02 CN CN201810175899.9A patent/CN108567422B/zh active Active
-
2019
- 2019-03-04 US US16/291,749 patent/US20190271752A1/en not_active Abandoned
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08266503A (ja) * | 1995-03-31 | 1996-10-15 | Shimadzu Corp | Mrイメージング装置 |
US8060180B2 (en) * | 2005-01-14 | 2011-11-15 | New York University | Multi-echo magnetic resonance imaging method and system |
CN102238909A (zh) * | 2008-12-04 | 2011-11-09 | 株式会社日立医疗器械 | 磁共振成像装置以及同步摄像方法 |
CN103371820A (zh) * | 2012-04-20 | 2013-10-30 | 西门子公司 | 借助磁共振技术采集呼吸的检查对象的测量数据组的方法 |
CN103389481A (zh) * | 2012-05-11 | 2013-11-13 | 上海联影医疗科技有限公司 | 磁共振频率相位双编码采样方法及其图像重建方法 |
CN103110420A (zh) * | 2013-01-31 | 2013-05-22 | 深圳先进技术研究院 | 磁共振血管成像方法和*** |
CN103961097A (zh) * | 2013-02-04 | 2014-08-06 | 上海联影医疗科技有限公司 | 一种磁共振扫描短te成像方法及磁共振扫描*** |
CN104635189A (zh) * | 2013-11-13 | 2015-05-20 | 上海联影医疗科技有限公司 | 磁共振多对比度图像重建方法及*** |
CN104545918A (zh) * | 2014-12-29 | 2015-04-29 | 中国科学院深圳先进技术研究院 | 非对比增强磁共振静脉造影成像方法 |
CN104688228A (zh) * | 2015-03-20 | 2015-06-10 | 中国科学院深圳先进技术研究院 | 一种磁共振成像方法及设备 |
CN105259525A (zh) * | 2015-10-28 | 2016-01-20 | 北京大学 | 基于邻域共享压缩感知的动态对比增强磁共振快速成像方法 |
CN106646302A (zh) * | 2015-10-29 | 2017-05-10 | 上海联影医疗科技有限公司 | 一种磁共振成像方法 |
US20170242087A1 (en) * | 2016-02-19 | 2017-08-24 | Samsung Electronics Co., Ltd. | Magnetic resonance imaging apparatus and method of obtaining magnetic resonance image by using multiple excitation with delayed spin-echoes |
CN106597333A (zh) * | 2016-12-30 | 2017-04-26 | 上海联影医疗科技有限公司 | 一种磁共振并行成像方法及磁共振成像*** |
CN106821381A (zh) * | 2017-01-06 | 2017-06-13 | 上海交通大学 | 一种血管成像方法 |
CN107126211A (zh) * | 2017-06-07 | 2017-09-05 | 上海联影医疗科技有限公司 | 磁共振成像方法和*** |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110148134A (zh) * | 2019-04-02 | 2019-08-20 | 佛山瑞加图医疗科技有限公司 | 血管成像方法及装置 |
CN113805130A (zh) * | 2020-06-17 | 2021-12-17 | 西门子(深圳)磁共振有限公司 | 快速磁化率敏感性成像方法、装置及磁共振成像*** |
CN113805130B (zh) * | 2020-06-17 | 2024-01-30 | 西门子(深圳)磁共振有限公司 | 快速磁化率敏感性成像方法、装置及磁共振成像*** |
CN113970717A (zh) * | 2020-07-24 | 2022-01-25 | 西门子(深圳)磁共振有限公司 | 磁敏感加权成像方法、装置及磁共振成像*** |
US20220026515A1 (en) * | 2020-07-24 | 2022-01-27 | Siemens Healthcare Gmbh | 3d flow compensated interleaved epi with a centric reordering scheme for fast high-resolution susceptibility-weighted imaging |
US11726158B2 (en) * | 2020-07-24 | 2023-08-15 | Siemens Healthcare Gmbh | 3D flow compensated interleaved EPI with a centric reordering scheme for fast high-resolution susceptibility-weighted imaging |
CN113970717B (zh) * | 2020-07-24 | 2023-11-24 | 西门子(深圳)磁共振有限公司 | 磁敏感加权成像方法、装置及磁共振成像*** |
Also Published As
Publication number | Publication date |
---|---|
US20190271752A1 (en) | 2019-09-05 |
CN108567422B (zh) | 2021-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ehman et al. | Magnetic resonance imaging with respiratory gating: techniques and advantages | |
van Geuns et al. | MR coronary angiography with breath-hold targeted volumes: preliminary clinical results | |
FI86775B (fi) | Kaernmagnetisk resonans anvaendande projektionsavbildningssystem foer en blodaoder. | |
CN110031786B (zh) | 磁共振图像重建方法、磁共振成像方法、设备及介质 | |
US6505064B1 (en) | Diagnostic imaging systems and methods employing temporally resolved intensity tracing | |
Wang et al. | Vision 20/20: Simultaneous CT‐MRI—Next chapter of multimodality imaging | |
EP3545837A1 (en) | Medical image providing apparatus and medical image processing method of the same | |
US10921405B2 (en) | Image reconstructing method and reconstructing apparatus | |
US10123761B2 (en) | Device and method for spatiotemporal reconstruction of a moving vascular pulse wave in the brain and other organs | |
EP2711896A2 (en) | Method and apparatus for image enhancement in magnetic resonance imaging using motion corrupted data | |
KR101967245B1 (ko) | 자기 공명 이미징 시스템 및 자기 공명 이미징 방법 | |
CN108567422A (zh) | 一种血管成像方法和装置 | |
CN109917315A (zh) | 磁共振成像扫描方法、装置、计算机设备和存储介质 | |
CN107865659A (zh) | 磁共振成像装置及获取磁共振图像的方法 | |
JP6625353B2 (ja) | 磁気共鳴イメージング装置 | |
CN104814737B (zh) | 磁共振成像装置以及磁共振成像方法 | |
CN105030238B (zh) | 在多个桌台位置处执行磁共振检查 | |
KR101826702B1 (ko) | 자기 공명 영상 장치 및 그 방법 | |
US20170131377A1 (en) | Magnetic resonance imaging apparatus and method | |
JP2016530020A (ja) | 押しボタン式血管壁イメージング | |
KR101649275B1 (ko) | 자기 공명 영상 생성 방법 및 그에 따른 의료 영상 장치 | |
US20120197103A1 (en) | Method For Calculation and Generation of Spatially-Tailored Parallel Radio Frequency Saturation Fields | |
US20120010496A1 (en) | medical image diagnosis apparatus and a controlling method | |
US20150374247A1 (en) | Method of measuring blood flow velocity performed by medical imaging apparatus, and the medical imaging apparatus | |
US20180100906A1 (en) | Magnetic resonance imaging apparatus and method of operating the same |
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 |