CN102575995A - 用于检测气泡的磁共振成像系统和方法 - Google Patents
用于检测气泡的磁共振成像系统和方法 Download PDFInfo
- Publication number
- CN102575995A CN102575995A CN2010800458861A CN201080045886A CN102575995A CN 102575995 A CN102575995 A CN 102575995A CN 2010800458861 A CN2010800458861 A CN 2010800458861A CN 201080045886 A CN201080045886 A CN 201080045886A CN 102575995 A CN102575995 A CN 102575995A
- Authority
- CN
- China
- Prior art keywords
- pixel
- magnetic resonance
- resonance imaging
- bubble
- mri
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
- A61N7/02—Localised ultrasound hyperthermia
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/42—Details of probe positioning or probe attachment to the patient
- A61B8/4272—Details of probe positioning or probe attachment to the patient involving the acoustic interface between the transducer and the tissue
- A61B8/4281—Details of probe positioning or probe attachment to the patient involving the acoustic interface between the transducer and the tissue characterised by sound-transmitting media or devices for coupling the transducer to the tissue
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
-
- 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/4808—Multimodal MR, e.g. MR combined with positron emission tomography [PET], MR combined with ultrasound or MR combined with computed tomography [CT]
- G01R33/4814—MR combined with ultrasound
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
- A61B2018/00636—Sensing and controlling the application of energy
- A61B2018/00666—Sensing and controlling the application of energy using a threshold value
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
- A61B90/37—Surgical systems with images on a monitor during operation
- A61B2090/374—NMR or MRI
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/082—Measurement of solid, liquid or gas content
-
- 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
-
- 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/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20224—Image subtraction
-
- 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
Abstract
一种用于检测成像体积(108)之内的气泡(124,148,304,306,404,406)的磁共振成像系统(100),所述磁共振成像系统包括:磁体(102),其适于生成磁场,磁场用于对位于所述成像体积之内的受试者(104)的核的磁自旋进行取向;射频系统(110,112),其适于采集磁共振数据(160,164),其中所述射频系统包括射频收发器(112)和射频线圈(110);磁场梯度线圈(114),其适于对所述成像体积之内的核的磁自旋进行空间编码;磁场梯度线圈电源(116),其适于向所述磁场梯度线圈供应电流;以及计算机系统(132),其适于从所述磁共振成像数据构造图像并控制所述磁共振成像系统的操作,其中所述计算机系统适于利用从所述磁共振成像数据构造的磁共振图像(162,166,300,400)检测所述成像体积之内的气泡。
Description
技术领域
本发明涉及图像处理,尤其是在磁共振图像中识别气泡。
背景技术
可以使用来自聚焦超声换能器的超声有选择地处置身体内部的区域。发射超声作为高能量机械振动。这些振动在受到阻尼时诱发组织发热,并且作为高强度场,它们还可能导致空化。组织发热和空化都可用于在临床环境中有意地破坏组织。不过,利用超声加热组织比空化更容易控制,期望目标外部的空化可能导致非预期热损伤,例如在患者皮肤处。可以使用超声处置来消融组织并有选择地杀死组织肿瘤区域。这种技术已经用于处置子宫肌瘤,减少了对子宫切除术的需求,子宫切除术常常需要长时间恢复和住院治疗。
为了有选择地处置组织,可以使用聚焦超声换能器在特定处置体积上聚焦超声。换能器能够发射超声,并且其通常安装在诸如脱气水的介质内。然后使用致动器调整超声换能器的位置,由此调整利用超声处理的受试者的区域。可以使用电聚焦对超声场进行引导或成形以消融更大的组织体积。
通常使用标准医学成像技术来规划超声程序,也可用于引导处置。磁共振成像(MRI)、计算机断层摄影(CT)和超声成像已被用于规划和引导超声处置。聚焦超声换能器通常具有可以致动其的有限范围,因此必须相对于超声系统适当定位患者。超声不能通过空气发射到身体,因此使用超声耦合介质,例如超声凝胶、超声凝胶衬垫和/或脱气水从超声系统向受试者表面发射超声。通常使用诸如聚酯薄膜的膜在超声换能器浴器(bath)和用于向受试者传导超声的介质之间形成分界。
发明内容
本发明在独立权利要求中提供了一种磁共振成像系统、一种计算机程序产品和一种计算机实施的方法。在从属权利要求中给出了实施例。
超声程序的困难在于,在利用超声耦合介质将受试者声学耦合到超声系统时,在受试者表面和超声耦合介质之间,例如在凝胶衬垫-皮肤界面处会俘获气泡。如果向存在显著大于波长的一个或多个气泡的区域中施加高强度聚焦超声,可能会将超声反射回发射换能器,诱发换能器损伤或减少换能器寿命。如果向存在一个或多个小气泡的区域中施加高强度聚焦超声,超声可能诱发一个或多个气泡中的空化,导致受试者与气泡接触的表面损伤或烧伤。因此,有益的是在向受试者施加高强度超声之前识别气泡。本发明的实施例可以通过利用图像处理识别磁共振图像中的气泡以提供这个技术问题的解决方案。在一些实施例中,可以使用气泡敏感脉冲序列采集磁共振成像数据。
本发明提供了用于检测成像体积内的气泡的磁共振成像系统。还可以将本发明用于检测多个气泡。该磁共振成像系统包括适于生成磁场的磁体,磁场用于对位于成像体积之内的受试者的核的磁自旋进行取向。该磁体可以是超导磁体、永磁体、电磁体或前面三种磁体的任何组合。该磁共振成像系统还包括适于采集磁共振数据的射频系统。这里使用的磁共振数据是磁共振成像系统从磁共振成像系统的成像体积之内的受试者采集的数据。
磁共振成像数据包括从受试者的核接收的射频(radio)信号表示或记录。磁共振成像数据可以用于构造受试者的成像区域的图像或可视化。主要从成像区域接收射频信号。不过,利用磁共振成像数据的傅里叶变换构造磁共振图像。结果,成像体积外部的核可能对磁共振成像数据有贡献,因此对构造的任何图像或可视化有贡献。
射频系统包括射频收发器和射频线圈。射频线圈可以是能够发射和接收射频信号的单个线圈,或者射频线圈可以由发射线圈和接收线圈构成。类似地,射频收发器可以是收发器,或可以是独立的发射和接收射频装置。该磁共振成像系统还包括适于对成像体积之内的磁自旋进行空间编码的磁场梯度线圈。该磁共振成像系统还包括适于向磁场梯度线圈供应电流的磁场梯度线圈电源。
该磁共振成像系统还包括计算机系统,适于从磁共振成像数据构造图像并控制磁共振成像系统的操作。计算机系统适于利用从磁共振成像数据构造的磁共振成像检测成像体积之内的气泡。检测成像体积之内的气泡是有益的,因为可以使用气泡的知识和位置进行处置规划或诊断。
在另一个实施例中,计算机系统适于通过访问磁共振图像检测成像体积之内的气泡。计算机系统还适于通过执行如下步骤来检测气泡:通过从趋势(trend)检测图像减去磁共振图像来创建减影图像。通过从趋势检测图像减去磁共振成像图像创建减影图像。
计算机系统还适于通过执行如下步骤来检测成像体积之内的气泡:通过向减影图像应用阈值来创建二值图像。计算机系统还适于通过执行如下步骤来检测成像体积之内的气泡:通过对二值图像进行中值滤波来创建输出图像。如这里使用的,中值滤波器是通过计算像素的近邻的中值从图像去除斑点噪声和黑白点相间噪声的滤波器。中值滤波器在去除噪声同时保留边缘方面是有效的。如这里使用的,图像的低通滤波器是空间滤波器。图像的低通滤波器或低通滤波执行去除特定尺寸以下的对象的步骤。
在另一个实施例中,所述磁共振成像系统适于利用气泡敏感脉冲序列采集磁共振成像数据。可以利用气泡敏感磁共振脉冲序列检测气泡。这里使用的脉冲序列是磁共振成像系统为了采集磁共振成像数据而执行的操作序列。例如,脉冲序列包含针对根据时间施加到磁共振梯度线圈的电压以及根据时间施加到射频线圈的射频信号的指令。气泡敏感序列的范例是梯度回波脉冲序列,其中:回波时间(TE)为15毫秒,翻转角(FA)为60度,体素尺寸为1.11mm乘1.11mm乘2.5mm,重建分辨率最大。使用很多脉冲序列对质子或氢原子密度成像。由于气泡中氢原子密度比诸如水的液体中低得多,因此很多不同的脉冲序列能够用于采集可用于检测气泡的磁共振成像数据。
在另一个实施例中,所述气泡敏感脉冲序列为梯度回波脉冲序列。梯度回波序列具有5到25毫秒的回波时间,优选介于10和20毫秒。梯度回波序列的分辨率介于0.8和1.5毫米之间,优选介于1和1.25毫米。在低于1毫米的分辨率下,发现补偿越来越大的信噪比是有益的。在低于1毫米的分辨率以下,还发现高估了气泡数量。在高于1.25毫米的分辨率下,发现高估了气泡数量。梯度回波序列的切片厚度介于1和4毫米之间,优选介于2和3毫米之间。
梯度回波序列的翻转角介于10和80度之间,优选介于50和70度之间。调整磁共振图像的重建分辨率可用于改善气泡检测的质量。可以使用二维或三维梯度脉冲序列。如果使用二维梯度脉冲序列检测气泡,在切片堆叠平行于超声窗口时,受试者和超声窗口之间的气泡可视性性能良好。例如,如果超声窗口安装在充当床的受试者支撑中,那么受试者将躺在床上。可以由磁共振成像系统使用在受试者的冠状面中采集的图像堆叠检测受试者和超声窗口之间的气泡。
在另一实施例中,磁共振成像系统还包括用于生成超声波的超声系统。超声系统包括超声窗口。超声窗口适于接收受试者。超声窗口可以是与受试者直接接触的窗口,或者超声窗口可以利用超声耦合介质与受试者形成接触。超声耦合介质是适于发射超声波并能够顺应表面形状的介质。超声耦合介质的范例是凝胶衬垫或超声耦合凝胶或水。
在另一实施例中,磁共振成像系统适于检测超声窗口和受试者之间的气泡。本实施例尤其有利,因为气泡影响到超声的传播。气泡可能导致超声散射,或者超声可能导致气泡空化。
在另一实施例中,超声系统是高强度聚焦超声系统。本实施例尤其有利,因为高强度超声可能导致气泡空化。检测气泡允许通过气泡空化对受试者或受试者一部分造成意外损伤。
在另一实施例中,该计算机系统适于计算超声波路径。可以利用射线跟踪技术计算超声波的路径。该计算机系统适于计算检测到的气泡沿路径造成的超声衰减。可以通过估计气泡半径的空化或振荡来计算超声波的衰减或者可以使用查找表格完成计算,查找表格使用了试验数据。磁共振成像系统还适于在超声波衰减高于预定衰减阈值时用信号通知操作员。本实施例是有利的,因为如果气泡将引起超声衰减或反射过大,它允许操作员重新定位受试者。可以通过多种方式用信号通知操作员。可以由可听信号,可以有计算机显示器或监视器上的信号,或者可以有信号指示符,例如控制面板上的灯,其可以指出衰减将过大。
在另一实施例中,超声系统适于在超过预定衰减阈值时调整超声波的路径。本实施例是有利的,因为使用超声波无需对受试者重新定位,以便从超声波的路径移除气泡。可以通过几种不同方式实现超声的路径。超声系统可以具有能够机械移动的超声换能器。在这种情况下,可以重新定位超声系统。在其他实施例中,超声系统可以具有超声换能器,该超声换能器具有多个超声换能器元件。个体元件可以具有电源,使得可以改变超声换能器元件生成的超声相位和振幅。这样允许将超声换能器元件选择为使得超声遵循未受检测的气泡衰减的路径。
在另一实施例中,该计算机系统包括显示器。显示器适于显示输出图像。可以显示输出图像自身,或者可以在磁共振成像图像上显示或叠加于其上。也可以显示输出图像,使其适于吸引操作员的注意。例如,可以使用亮色,例如红色、黄色或橙色指示气泡位置。可以使用示出解剖结构的磁共振成像图像上叠加的闪烁图像突出所检测气泡的位置。
在另一方面中,本发明提供了一种包括机器可执行指令的计算机程序产品,机器可执行指令由计算机系统执行,以执行访问从成像体积的磁共振成像数据构造的磁共振图像的步骤。磁共振图像可以存储在例如硬盘驱动器或其他计算机存储介质上,或者可以位于计算机系统的随机存取存储器之内。换言之,磁共振图像可以是先前采集的,或一旦从磁共振成像数据构造出磁共振图像,就可以开始检测气泡的过程。这里使用的计算机系统是适于执行机器可执行指令的机器。计算机系统的范例是单一计算机、计算机网络、计算机群、嵌入式系统、诸如蜂窝电话或个人数字助理的移动计算装置和微控制器。机器可执行指令还包括利用磁共振图像执行成像体积之内气泡检测的步骤。
在另一实施例中,计算机程序产品还包括用于执行通过检测磁共振图像中对比度趋势来创建趋势检测图像的步骤的指令。计算机程序产品还包括用于执行通过从趋势检测图像减去磁共振图像来创建减影图像的步骤的指令。
计算机程序产品还包括通过向减影图像施加阈值创建二值图像的步骤。计算机程序产品还包括通过对二值图像进行中值滤波来执行创建输出图像的步骤。在本发明的计算机程序产品另一实施方式中,通过执行如下步骤检测对比度趋势:对磁共振图像进行中值滤波以确保第一多个像素的每个的数值都在预定范围之内;并对所得的中值滤波磁共振图像进行低通滤波。
在另一个实施例中,磁共振图像包括第一多个像素。趋势检测图像包括第二多个像素。减影图像包括第三多个像素。二值图像包括第四多个像素。第一多个像素的每个包括数值。第二多个像素的每个包括数值。第三多个像素的每个包括数值。对于第一多个像素的每个,在第一多个像素、第二多个像素、第三多个像素和第四多个像素中都有相应像素。通过检测第一多个像素中的对比度趋势计算第二多个像素的每个的数值。通过从第二多个像素的相应像素减去第一多个像素的相应像素或从第一多个像素的相应像素减去第二多个像素的相应像素,计算第三多个像素的每个的数值。如果第三多个像素的相应像素高于预定阈值,为第四多个像素的每个像素分配第一值。如果第三多个像素的相应像素低于预定阈值,为第四多个像素的每个像素分配第二值。
在由计算机实施时,第一、第二、第三和第四多个像素可以在一些实施方式中是相同的。例如,在不再需要多个像素时,可以对其进行重写并重新使用以创建另一多个像素。另一个范例是,在创建二值图像时,可以用二值图像重写第三多个像素,以便创建第四多个像素。
在另一实施例中,通过执行对磁共振图像的中值滤波的步骤以确保多个像素的每个的数值在预定范围之内从而检测对比度趋势。通过执行对所得中值滤波图像进行低通滤波的步骤进一步检测对比度趋势。实质上,通过首先从像素去除极值,然后施加初始去除小细节的低通滤波器,检测对比度趋势。在这两个步骤之后,图像中剩余的是图像对比度的大趋势。
在另一方面中,本发明提供了一种计算机实施的方法,用于检测磁共振图像中的气泡。该方法包括访问磁共振图像的步骤。该方法还包括通过检测磁共振图像中的对比度趋势来创建趋势检测图像。该方法还包括通过从趋势检测图像减去磁共振图像创建减影图像。
该方法还包括通过向减影图像应用阈值来创建二值图像。该方法还包括通过对二值图像进行中值滤波创建输出图像的步骤。
附图说明
在下文中将仅通过举例,并参考附图描述本发明的优选实施例,在附图中:
图1图示了根据本发明实施例的磁共振成像系统;
图2示出的方框图图示了根据本发明实施例的方法;
图3a示出了具有气泡的凝胶体模的磁共振图像;
图3b示出了从图3a生成的识别气泡的输出图像;
图4a示出了皮肤和超声窗口界面的磁共振图像;以及
图4b示出了从图4a生成的识别气泡的输出图像。
附图标记列表
100 磁共振成像系统
102 磁体
104 受试者
106 受试者支撑
108 成像区
110 射频线圈
112 射频收发器
114 磁场梯度线圈
116 磁场梯度线圈电源
118 超声系统
120 聚焦超声换能器
122 凝胶衬垫
124 气泡
126 超声路径
128 超声焦点
130 超声窗口
132 计算机
134 硬件接口
136 微处理器
138 存储器
140 存储设备
142 用户接口
144 显示器
146 磁共振图像
148 气泡
150 计算机程序产品
152 磁共振成像系统控制模块
154 超声系统控制模块
156 磁共振图像构造模块
158 气泡识别模块
160 磁共振成像数据
162 磁共振图像
164 磁共振成像数据存档
166 磁共振图像存档
168 计算机程序产品
300 磁共振图像
302 输出图像
304 气泡
306 白色区域
400 磁共振图像
402 输出图像
404 气泡
406 白色区域
具体实施方式
这些附图中的编号类似的元件是等价元件或执行相同功能。如果功能等价,先前论述过的元件未必会在后面的图中加以论述。
图1示出了根据本发明实施例的磁共振成像系统100。磁共振成像系统100具有磁体102。在本范例中,示出了具有通过中心的膛的圆柱形磁体。也可以使用其他各种磁体。例如,也可以使用类似两个亥姆霍兹线圈的开放式磁体。在磁体102的膛之内是定位于受试者支撑106上的受试者104。磁共振成像系统100适于采集成像区108之内的磁共振成像数据。射频天线110采集成像区108之内的磁共振成像数据。射频线圈110连接到射频收发器112。磁体的膛内还有磁场梯度线圈114。磁场梯度线圈114由磁场梯度电源116供电。
受试者支撑106下方是超声系统118。在本范例中,示出了用于高强度聚焦超声系统中的类型的聚焦超声收发器120。超声系统118具有超声窗口130,其适于允许超声波从超声系统118通过进入受试者104体内。在本范例中,在受试者支撑106之内有间隙,适于接收凝胶衬垫122,凝胶衬垫122将超声窗口130耦合到受试者104。在本范例中,在凝胶衬垫122和受试者104之间示出了气泡124。不过,气泡124位于成像区108之内。超声通过受试者的体积由标记为126的线划界。超声被聚焦在受试者之内的体积128中。在这幅图中可以看出,气泡124在超声通过的边界126之内。不过,气泡在成像区108之内,可以在磁共振图像146中被识别出。
磁共振成像系统100还包括计算机系统132。计算机系统具有硬件接口134,适于将计算机系统132连接到磁共振成像系统100的各其他部件。
射频收发器112、超声系统118和磁场梯度电源116都被示为连接到硬件接口134。计算机系统132之内是微处理器136。微处理器136连接到硬件接口134、用户接口142、计算机存储设备140和计算机存储器138。微处理器136适于执行机器可执行指令。计算机存储器138适于存储数据和机器可执行指令。
计算机存储器存储适于控制和运行磁共振成像系统100的计算机程序产品150。计算机程序产品150包括磁共振成像系统控制模块152,包含用于控制和操作磁共振成像系统的指令。计算机程序产品150还包括超声控制模块154,包含用于控制超声系统118的操作的机器可执行指令。计算机程序产品150还包括磁共振图像构造模块156。这个模块包含用于从磁共振成像数据160构造磁共振图像和透视图的例程和算法。计算机程序产品150还包括气泡识别模块158,适于对磁共振图像162执行图像分析,以识别气泡124。存储器138还具有磁共振成像数据160和磁共振图像162。
计算机系统132还具有计算机存储设备140。计算机存储设备的形式可以是固态硬盘驱动器、常规磁盘硬盘驱动器或其他适当的计算机存储介质。计算机存储设备140能够储存数据,以及计算机程序产品168的副本。在计算机存储设备中是磁共振成像数据164的存档以及磁共振图像166的存档。
还连接到微处理器136的是用户接口142。用户接口包括允许操作员观看数据或信息,并且还控制计算机系统132的操作的元件。例如,用户接口可以包括鼠标、触摸板、键盘或其他用户输入装置。图1中还示出了显示器144,其还包括用户接口142的一部分。显示器144上显示了磁共振图像146。叠加在磁共振图像146上的是利用气泡识别模块158识别的气泡148。
在工作中,将把受试者104放在受试者支撑106上。磁共振成像系统100可用于提供解剖学数据,解剖学数据可用于引导超声系统118的操作。在开始超声处理程序之前,磁共振成像系统在受试者104和凝胶衬垫122之间的界面区域中采集磁共振成像数据。气泡识别模块158然后识别受试者104和凝胶衬垫122之间的气泡124。利用显示器144向操作员通知这种情况。操作员然后能够决定移动受试者104并移除气泡124或调整超声换能器120的操作。这可以涉及到移动超声换能器120或控制在程序中激活超声换能器的哪些元件。
图2示出了根据本发明实施例的方法实施例。在步骤200中,访问磁共振成像图像。可以访问存储的图像,或可以使用刚从采集的磁共振成像数据构造的磁共振图像。在步骤202中,创建在其中检测到大对比度趋势的图像。在步骤204中,从趋势检测图像减去磁共振图像。在步骤206中,然后向减影图像应用阈值,最后在步骤208中构造输出图像。通过向二值图像应用中值滤波器构造输出图像。
图3a示出了其中嵌入了气泡的凝胶体模的磁共振图像300。图3b示出了输出图像302,其中的气泡已经在图3a所示的图像中识别出来。在图3b中气泡被识别为图像中的白色区域。例如,图3a中识别了气泡300,图3b中识别了与气泡300对应的白色区域302。
图4a示出了人皮肤和凝胶衬垫之间界面的磁共振图像400。在图4b中示出的输出图像402中将图4a中的气泡识别为白色区域。图4a和4b展示出,可以将方法实施例用于识别超声窗口和受试者表面之间的气泡。例如,图4a中识别了气泡400,图4b中识别了与气泡400对应的白色区域402。
Claims (15)
1.一种磁共振成像系统,其包括用于生成超声波的超声系统(118);其中,所述超声系统包括超声窗口(130);且其中,所述超声窗口适于接收受试者,其中,所述磁共振成像系统(100)被配置成检测成像体积(108)之内的气泡(124,148,304,306,404,406),所述磁共振成像系统包括:
-磁体(102),其适于生成磁场,该磁场用于对位于所述成像体积之内的所述受试者(104)的核的磁自旋进行取向;
-射频系统(110,112),其适于采集磁共振数据(160,164),其中,所述射频系统包括射频收发器(112)和射频线圈(110);
-磁场梯度线圈(114),其适于对所述成像体积之内的核的磁自旋进行空间编码;
-磁场梯度线圈电源(116),其适于向所述磁场梯度线圈供应电流;以及
-计算机系统(132),其适于从所述磁共振成像数据构造图像并控制所述磁共振成像系统的操作,其中,所述计算机系统适于利用从所述磁共振成像数据构造的磁共振图像(162,166,300,400)检测所述成像体积之内的所述气泡,且其中
所述超声系统被配置成基于所检测的气泡调整所述超声波的路径。
2.根据权利要求1所述的磁共振成像系统,其中,所述计算机系统适于通过执行如下步骤检测所述成像体积之内的所述气泡:
-访问(200)所述磁共振图像;
-通过检测所述磁共振图像中的对比度趋势创建(202)趋势检测图像;
-通过从所述趋势检测图像减去所述磁共振图像创建(204)减影图像;
-通过向第三多个像素的每个应用阈值创建(206)二值图像;以及
-通过对所述二值图像进行中值滤波创建(208)输出图像(146,302,402)。
3.根据权利要求1或2所述的磁共振成像系统,其中,所述磁共振成像系统适于利用气泡敏感脉冲序列采集磁共振成像数据。
4.根据权利要求3所述的磁共振成像系统,其中,所述气泡敏感脉冲序列为梯度回波序列,其中,所述梯度回波序列具有5和25毫秒之间的回波时间,其中,所述梯度回波序列具有0.8和1.5毫米之间的分辨率,其中,所述梯度回波序列具有1和4毫米之间的切片厚度,且其中,所述梯度回波序列具有10和80度之间的翻转角。
5.根据前述权利要求中的任一项所述的磁共振成像系统,其中,包括超声换能器的所述超声系统被配置成通过机械移动所述超声换能器来调整超声束的路径。
6.根据前述权利要求中的任一项所述的磁共振成像系统,其中,包括具有多个换能器元件的超声换能器阵列的所述超声系统被配置成通过调节由所述超声换能器的换能器元件生成的超声的相位和振幅来调整超声束的路径。
7.根据权利要求5或6所述的磁共振成像系统,其中,所述磁共振成像系统适于检测所述超声窗口和所述受试者之间的所述气泡。
8.根据权利要求5、6或7所述的磁共振成像系统,其中,所述超声系统是高强度聚焦超声系统(118,120,130)。
9.根据权利要求5到8中的任一项所述的磁共振成像系统,其中,所述计算机系统还适于计算所述超声波的路径(126);其中,所述计算机系统还适于计算所述气泡沿所述路径造成的所述超声波的衰减;且其中,所述磁共振成像系统还适于在所述超声波的衰减高于预定衰减阈值时用信号通知操作员。
10.根据权利要求9所述的磁共振成像系统,其中,所述超声系统适于在超过所述预定衰减阈值时调整所述超声波的路径。
11.根据前述权利要求中的任一项所述的磁共振成像系统,其中,所述计算机系统包括显示器(144);且其中,所述显示器适于显示所述输出图像。
12.一种用于检测磁共振图像(162,166,300,400)之内的气泡(124,148,304,306,404,406)的计算机程序产品,所述计算机程序产品包括用于由计算机系统(132)执行以执行以下步骤的机器可执行指令:
-访问(200)从成像体积(108)的磁共振成像数据(160,164)构造的磁共振图像(162,166,166,400);以及
-利用所述磁共振图像检测所述成像体积之内的所述气泡。
13.根据权利要求11所述的计算机程序产品,其中,利用所述磁共振图像检测所述成像体积之内的所述气泡包括如下步骤:
-通过检测所述磁共振图像中的对比度趋势创建(202)趋势检测图像;
-通过从所述趋势检测图像减去所述磁共振图像创建(204)减影图像;
-通过向第三多个像素的每个应用阈值创建(206)二值图像;以及
-通过对所述二值图像进行中值滤波创建(208)输出图像(146,302,402)。
14.根据权利要求13所述的计算机程序产品,其中,所述磁共振图像包括第一多个像素;其中,所述趋势检测图像包括第二多个像素;其中,所述减影图像包括第三多个像素;其中,所述二值图像包括第四多个像素;其中,所述第一多个像素的每个都包括数值;其中,所述第二多个像素的每个都包括数值;其中,所述第三多个像素的每个都包括数值;其中,对于所述第一多个像素的每个,在所述第一多个像素、所述第二多个像素、所述第三多个像素和所述第四多个像素中都有相应像素;其中,通过检测所述第一多个像素中的对比度趋势计算所述第二多个像素的每个的数值;其中,通过从所述第二多个像素的相应像素减去所述第一多个像素的相应像素或通过从所述第一多个像素的相应像素减去所述第二多个像素的相应像素计算所述第三多个像素的每个的数值;其中,如果所述第三多个像素的相应像素高于预定阈值,为所述第四多个像素的每个像素分配第一值;且其中,如果所述第三多个像素的相应像素低于所述预定阈值,为所述第四多个像素的每个分配第二值。
15.一种用于检测磁共振图像(162,166,300,400)中的气泡(124,148,304,306,404,406)的计算机实施的方法,所述方法包括:
-访问(200)磁共振图像;
-通过检测所述磁共振图像中的对比度趋势创建(202)趋势检测图像;
-通过从所述趋势检测图像减去所述磁共振图像创建(204)减影图像;
-通过向第三多个像素的每个应用阈值创建(206)二值图像;以及
-通过对所述二值图像进行中值滤波创建(208)输出图像(146,302,402)。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP09172765A EP2312303A1 (en) | 2009-10-12 | 2009-10-12 | Magnetic resonance imaging system and method for detecting a gas bubble |
EP09172765.1 | 2009-10-12 | ||
PCT/IB2010/054492 WO2011045708A1 (en) | 2009-10-12 | 2010-10-05 | Magnetic resonance imaging system and method for detecting a gas bubble |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102575995A true CN102575995A (zh) | 2012-07-11 |
Family
ID=41651139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010800458861A Pending CN102575995A (zh) | 2009-10-12 | 2010-10-05 | 用于检测气泡的磁共振成像系统和方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US20120203097A1 (zh) |
EP (2) | EP2312303A1 (zh) |
JP (1) | JP2013507165A (zh) |
CN (1) | CN102575995A (zh) |
BR (1) | BR112012008223A2 (zh) |
RU (1) | RU2012119445A (zh) |
WO (1) | WO2011045708A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107106872A (zh) * | 2014-09-02 | 2017-08-29 | 波尔多大学 | 通过超声用于心脏区域的切除或监测的系统和有关方法 |
CN108426908A (zh) * | 2018-02-09 | 2018-08-21 | 上海大学 | 一种检测岩心中泡沫气液比的方法 |
CN109959886A (zh) * | 2017-12-22 | 2019-07-02 | 西门子保健有限责任公司 | 用于确定磁共振成像装置的成像质量信息的方法、装置 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20130009138A (ko) * | 2011-07-14 | 2013-01-23 | 삼성전자주식회사 | 집속 초음파 치료 장치 및 이의 초점 제어 방법 |
BR112015003632A2 (pt) | 2012-08-22 | 2017-09-26 | Miitors Aps | medidor ultrassônico de fluxo, método de montagem de um medidor ultrassônico de fluxo. |
US10589129B2 (en) * | 2016-09-14 | 2020-03-17 | Insightec, Ltd. | Therapeutic ultrasound with reduced interference from microbubbles |
CN108919152B (zh) * | 2018-04-20 | 2024-02-06 | 福建省计量科学研究院(福建省眼镜质量检验站) | 振动与冲击传感器的磁灵敏度三维检定系统及方法 |
US11403760B2 (en) | 2019-02-21 | 2022-08-02 | Medtronic Navigation, Inc. | Method and apparatus for magnetic resonance imaging thermometry |
US11276174B2 (en) * | 2019-02-21 | 2022-03-15 | Medtronic Navigation, Inc. | Method and apparatus for magnetic resonance imaging thermometry |
US11426229B2 (en) | 2019-02-21 | 2022-08-30 | Medtronic Navigation, Inc. | Method and apparatus for magnetic resonance imaging thermometry |
CN112987017B (zh) * | 2019-11-29 | 2024-04-19 | 深圳迈瑞生物医疗电子股份有限公司 | 超声成像系统及其掉电控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0627206A2 (en) * | 1993-03-12 | 1994-12-07 | Kabushiki Kaisha Toshiba | Method and apparatus for ultrasound medical treatment |
US6128522A (en) * | 1997-05-23 | 2000-10-03 | Transurgical, Inc. | MRI-guided therapeutic unit and methods |
US20030083597A1 (en) * | 2001-10-29 | 2003-05-01 | Insightec-Txsonics Ltd. | System and method for sensing and locating disturbances in an emergy path of a focused ultrasound system |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH03136642A (ja) * | 1989-10-20 | 1991-06-11 | Olympus Optical Co Ltd | 超音波治療装置 |
JP3325300B2 (ja) * | 1992-02-28 | 2002-09-17 | 株式会社東芝 | 超音波治療装置 |
JP3860227B2 (ja) * | 1993-03-10 | 2006-12-20 | 株式会社東芝 | Mriガイド下で用いる超音波治療装置 |
JPH06315541A (ja) * | 1993-03-12 | 1994-11-15 | Toshiba Corp | 画像診断装置を用いた治療装置 |
JPH06261907A (ja) * | 1993-03-15 | 1994-09-20 | Toshiba Corp | 超音波治療装置 |
DK0711179T3 (da) * | 1993-07-30 | 2005-02-14 | Imcor Pharmaceutical Co | Stabiliserede sammensætninger med mikrobobler til ultralyd |
JPH0824267A (ja) * | 1994-07-15 | 1996-01-30 | Toshiba Corp | 超音波治療装置 |
JP3343031B2 (ja) * | 1996-06-27 | 2002-11-11 | 株式会社マキ製作所 | 核磁気共鳴を用いた青果物自動選別方法及びその装置 |
US5752515A (en) * | 1996-08-21 | 1998-05-19 | Brigham & Women's Hospital | Methods and apparatus for image-guided ultrasound delivery of compounds through the blood-brain barrier |
JP4322322B2 (ja) * | 1997-03-31 | 2009-08-26 | 株式会社東芝 | 超音波治療装置 |
JP2000065768A (ja) * | 1998-08-24 | 2000-03-03 | Natl Food Res Inst | 炊飯米の品質評価方法 |
JP4406107B2 (ja) * | 1999-03-09 | 2010-01-27 | 株式会社東芝 | 超音波治療装置 |
WO2001017455A2 (en) * | 1999-09-10 | 2001-03-15 | Transurgical, Inc. | Occlusion of tubular anatomical structures by energy application |
JP2001104358A (ja) * | 1999-10-08 | 2001-04-17 | Toshiba Corp | 超音波治療装置 |
US6692442B2 (en) * | 2001-12-13 | 2004-02-17 | Koninklijke Philips Electronics N.V. | Device for producing an on-line image of a body part into which a contrasting agent has been introduced |
WO2006129099A1 (en) * | 2005-06-02 | 2006-12-07 | Cancercure Technology As | Ultrasound treatment system |
US20070083120A1 (en) * | 2005-09-22 | 2007-04-12 | Cain Charles A | Pulsed cavitational ultrasound therapy |
EP1780651A1 (en) * | 2005-10-25 | 2007-05-02 | Bracco Imaging, S.P.A. | Method and system for automatic processing and evaluation of images, particularly diagnostic images |
CN100574829C (zh) * | 2006-08-24 | 2009-12-30 | 重庆融海超声医学工程研究中心有限公司 | 一种影像设备引导的高强度聚焦超声治疗系统 |
JP5551936B2 (ja) * | 2006-09-29 | 2014-07-16 | コーニンクレッカ フィリップス エヌ ヴェ | ハンズフリー超音波診断装置 |
JP5350275B2 (ja) * | 2007-02-23 | 2013-11-27 | コーニンクレッカ フィリップス エヌ ヴェ | 超音波装置、治療システム及びワークフローを向上させる方法 |
US8251908B2 (en) * | 2007-10-01 | 2012-08-28 | Insightec Ltd. | Motion compensated image-guided focused ultrasound therapy system |
-
2009
- 2009-10-12 EP EP09172765A patent/EP2312303A1/en not_active Ceased
-
2010
- 2010-10-05 CN CN2010800458861A patent/CN102575995A/zh active Pending
- 2010-10-05 EP EP10768585A patent/EP2488857A1/en not_active Withdrawn
- 2010-10-05 BR BR112012008223A patent/BR112012008223A2/pt not_active IP Right Cessation
- 2010-10-05 JP JP2012532704A patent/JP2013507165A/ja not_active Ceased
- 2010-10-05 WO PCT/IB2010/054492 patent/WO2011045708A1/en active Application Filing
- 2010-10-05 US US13/501,554 patent/US20120203097A1/en not_active Abandoned
- 2010-10-05 RU RU2012119445/14A patent/RU2012119445A/ru not_active Application Discontinuation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0627206A2 (en) * | 1993-03-12 | 1994-12-07 | Kabushiki Kaisha Toshiba | Method and apparatus for ultrasound medical treatment |
US6128522A (en) * | 1997-05-23 | 2000-10-03 | Transurgical, Inc. | MRI-guided therapeutic unit and methods |
US20030083597A1 (en) * | 2001-10-29 | 2003-05-01 | Insightec-Txsonics Ltd. | System and method for sensing and locating disturbances in an emergy path of a focused ultrasound system |
Non-Patent Citations (3)
Title |
---|
DC CHENG ET AL.: "Bubble recognition from image sequences", 《INVERSE PROBLEMS AND EXPERIMENTAL》 * |
JONATHAN HINDLEY ET AL.: "MRI Guidance of Focused Ultrasound Therapy of Uterine Fibroids:Early Results", 《AMERICAN JOURNAL OF ROENTGENOLOGY》 * |
MICHAEL LAMBA ET AL.: "Magnetic resonance imaging of microbubbles in a superheated emulsion chamber for brachytherapy dosimetry", 《MEDICAL PHYSICS》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107106872A (zh) * | 2014-09-02 | 2017-08-29 | 波尔多大学 | 通过超声用于心脏区域的切除或监测的系统和有关方法 |
CN109959886A (zh) * | 2017-12-22 | 2019-07-02 | 西门子保健有限责任公司 | 用于确定磁共振成像装置的成像质量信息的方法、装置 |
CN109959886B (zh) * | 2017-12-22 | 2021-09-10 | 西门子保健有限责任公司 | 用于确定磁共振成像装置的成像质量信息的方法、装置 |
CN108426908A (zh) * | 2018-02-09 | 2018-08-21 | 上海大学 | 一种检测岩心中泡沫气液比的方法 |
Also Published As
Publication number | Publication date |
---|---|
EP2312303A1 (en) | 2011-04-20 |
RU2012119445A (ru) | 2013-11-20 |
JP2013507165A (ja) | 2013-03-04 |
BR112012008223A2 (pt) | 2019-09-24 |
EP2488857A1 (en) | 2012-08-22 |
US20120203097A1 (en) | 2012-08-09 |
WO2011045708A1 (en) | 2011-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102575995A (zh) | 用于检测气泡的磁共振成像系统和方法 | |
US8725232B2 (en) | Therapeutic apparatus | |
EP2681576B1 (en) | Accelerated mr thermometry mapping involving an image ratio constrained reconstruction | |
CN102695963B (zh) | 用于核磁共振中计算局部比能量吸收率(sar)的方法 | |
US9420989B2 (en) | Apparatus, method and computer program product for calculating the speed of ultrasound in at least two tissue types | |
CN104125802B (zh) | 使用冷却周期期间采集的磁共振数据对处置计划的修改 | |
Mei et al. | Combining two‐dimensional spatially selective RF excitation, parallel imaging, and UNFOLD for accelerated MR thermometry imaging | |
CN101888876A (zh) | 运动补偿的影像引导聚焦超声治疗系统 | |
US20210204915A1 (en) | Focused ultrasound system with optimized monitoring of cavitation | |
CN102405078A (zh) | 用于控制治疗设备的控制设备 | |
JP2012130701A (ja) | データを伝送するためのシステム及び方法 | |
CN103180876A (zh) | 用于识别医学图像中的处置区域的医学成像系统、计算机实现方法和计算机程序产品 | |
JP2011224355A (ja) | 画像処理装置及び画像処理装置の制御プログラム | |
Winkler et al. | Direct SAR mapping by thermoacoustic imaging: A feasibility study | |
CN105072983B (zh) | 用于诊断成像的热疗 | |
US20180164392A1 (en) | Magnetic resonance imaging apparatus and controlling method thereof | |
CN105030238A (zh) | 在多个桌台位置处执行磁共振检查 | |
US20190170838A1 (en) | Coil apparatus, magnetic resonance imaging apparatus, and method of controlling the coil apparatus | |
CN106102576A (zh) | 动态医疗图像中运动的检测 | |
Bitton et al. | MR‐acoustic radiation force imaging (MR‐ARFI) and susceptibility weighted imaging (SWI) to visualize calcifications in ex vivo swine brain | |
JP2007260001A (ja) | 磁気共鳴イメージング装置および撮像方法 | |
JP4125134B2 (ja) | 磁気共鳴アコーストグラフィ | |
KR102468547B1 (ko) | 혈관 조영 영상을 생성하기 위한 자기 공명 영상 생성 장치 및 방법 | |
JP5484272B2 (ja) | 磁気共鳴イメージング装置および受信コイル接続状態の確認方法 | |
Kaye | Novel MRI Tools for Focused Ultrasound Surgery |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C05 | Deemed withdrawal (patent law before 1993) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120711 |