CN105308473A - 平行成像加速参数的自动优化 - Google Patents

平行成像加速参数的自动优化 Download PDF

Info

Publication number
CN105308473A
CN105308473A CN201480014017.0A CN201480014017A CN105308473A CN 105308473 A CN105308473 A CN 105308473A CN 201480014017 A CN201480014017 A CN 201480014017A CN 105308473 A CN105308473 A CN 105308473A
Authority
CN
China
Prior art keywords
reduction factor
factor
coil
snr
reduction
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
Application number
CN201480014017.0A
Other languages
English (en)
Other versions
CN105308473B (zh
Inventor
P·博尔纳特
M·P·J·于里森
M·I·多内瓦
A·J·W·杜杰达姆
K·内尔克
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN105308473A publication Critical patent/CN105308473A/zh
Application granted granted Critical
Publication of CN105308473B publication Critical patent/CN105308473B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34—Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/34092—RF coils specially adapted for NMR spectrometers
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • General Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

一种平行磁共振成像系统(1)包括：具有多个线圈元件的至少一个射频(RF)线圈(10，12)、智能选择单元(24)、平行成像参数单元(28)以及序列控制(16)。智能选择单元(24)根据利用所述至少一个RF线圈的对对象的预扫描或先前扫描，构建(60)信号图和基于不同的减小因子集合的多个噪声图。平行成像参数单元(28)选择包括最高信噪比(SNR)噪声图相对应的减小因子的集合。序列控制(16)基于所选择的减小因子来执行对所述对象的磁共振成像扫描。

Description

平行成像加速参数的自动优化
技术领域
下文总体上涉及医学成像。其具体结合磁共振成像获得应用，并且特别参考磁共振成像进行描述。然而，将要理解的是，其也在其他场景使用而不必限制于上述应用。
背景技术
在利用平行成像的3维(3-D)磁共振(MR)扫描或利用在切片方向的平行成像的多切片扫描中，使用具有多个线圈元件和/或线圈阵列的局部线圈和/或全身线圈。例如，全身线圈可以包括平行操作的8、16、32等个线圈元件。线圈元件提供用于重建图像的交叠的空间信息。取决于扫描的几何配置和在扫描中使用的线圈的设计，方向不同，平行成像的性能不同。
针对每个扫描，操作者可以选择针对一个或多个编码方向(诸如切片编码方向、相位编码方向或频率编码方向)的加速因子。加速因子通过减小共振激励和充分采样k空间所需要的采样周期来减小扫描的时间。加速的量通常称为减小因子，其中，R＝1是指没有减小，R＝2是指减小了因子2等。编码方向可以包括，例如，一个或多个相位编码方向、读取方向或切片方向。减小因子可以是分数或者是整数。例如，减小因子可以包括1、8/7、4/3、1.5、8/5、2、2.5、8/3、3、3.5、4等。差的或不正确的选择导致低的信噪比(SNR)，而低的SNR导致差的图像质量。对减小因子的差的或不正确的选择只有在扫描完成和结果图像被查看之后才被操作者意识到。另外，各个加速因子的量提供了针对每个扫描重复的复杂的设置，这使得系统不友好并且操作者错误的可能性很大。
其他方法包括选择线圈元件的“最优”集合来执行对给定的感兴趣区域的MR扫描。选择线圈元件的最优集合的过程包括调查扫描或预扫描，其可以提供噪声信息。备选地，对象的先前扫描可以提供噪声信息。
下文公开了新的且改进的对平行加速参数的自动优化，其解决了上面提到的问题，以及其他问题。
发明内容
根据一个方面，平行磁共振成像系统包括具有多个线圈元件的至少一个射频(RF)线圈、智能选择单元、平行成像参数单元以及序列控制。智能选择单元，根据利用所述至少一个RF线圈的对对象预扫描或先前扫描，构建信号图和基于沿着一个或多个减小方向的减小因子的不同集合的多个噪声图。平行成像参数单元选择与优化信噪比(SNR)和扫描速度的噪声图相对应的减小因子的一个集合。序列控制基于所选择的减小因子来执行对对象的磁共振成像扫描。
根据另一方面，平行磁共振成像的方法包括根据对对象的预扫描基于接收到的总体减小因子来构建信号图和多个噪声图，每个噪声图基于沿着一个或多个减小方向的减小因子的不同集合而被确定。选择与满足信噪比(SNR)标准和扫描速度标准的噪声图相对应的减小因子。基于所选择的减小因子来执行对对象的磁共振成像扫描。
根据另一方面，平行磁共振成像系统包括具有多个线圈元件的至少一个RF线圈、用户接口和一个或多个处理器。用户接口接收最小图像质量和/或生成图像数据的最大时间中的至少一个。一个或多个处理器被配置为根据利用所述至少一个RF线圈的对对象预扫描或先前扫描，构建信号图和针对沿着一个或多个减小方向的多个减小因子中的每个的多个噪声图。所述一个或多个处理器还被配置为选择减小因子，其生成具有至少接收到的最小图像质量和/或在接收到的最大时间内的图像数据，并且使得成像系统基于所选择的减小因子来执行对对象的磁共振成像扫描。
一个优点是减小了要被选择的参数的数量。
另一优点存在于基于诸如总体减小因子、最小信噪比(SNR)等可量化的参数来自动选择优化的编码方向和/或减小因子。
另一优点存在于使用从常规线圈调查扫描获得的信息，其包括线圈敏感度信息和线圈的患者负载，以用于优化编码方向和/或减小因子。
另一优点存在于简化用户接口。
在阅读和理解以下的详细说明之后，本领域普通技术人员将会意识到再其他的优点。
附图说明
本发明可以采取各种部件和各部件布置以及各种步骤和各步骤安排的形式。附图仅用于说明优选实施例的目的，而不应被理解为限制本发明。
图1示意性地示出了对平行加速参数的进行自动优化的系统的实施例。
图2是智能选择处理的一个实施例的流程图。
图3是对平行加速参数进行自动优化的一种方法的流程图。
具体实施方式
参考图1，示意性地示出了对平行加速参数进行自动优化的系统1的实施例。系统1包括诸如水平膛扫描器、开放系统、C型扫描器等的MR扫描器2。扫描器包括定义检查区域的开口或膛，对象4被放置于所述检查区域中以进行成像。MR扫描器2包括主磁体6、一个或多个梯度线圈8和多线圈元件、全身射频(RF)线圈10。主磁体6可以包括超导线圈或常导线圈，并且产生静态B0场，例如水平静态场或垂直静态场。梯度线圈应用磁场梯度来对对象组织中的共振进行方向性编码，并且产生梯度回波。方向性编码可以定义相位编码方向、准备方向、切片编码方向等。在三维(3D)扫描中，梯度线圈通常在正交于读取方向的两个正交方向进行相位编码，而不是定义和重建2D图像的堆叠。扫描器2可以包括一个或多个局部RF线圈12，诸如躯干线圈、头线圈、前列腺线圈、腿线圈等。全身RF线圈8和/或局部RF线圈12包括具有多个线圈元件或线圈阵列的至少一个线圈，其基于诸如敏感度编码(SENSE)的平行成像技术提供平行成像。全身RF线圈和/或局部RF线圈生成RF脉冲，所述RF脉冲引起和操控对象4组织中的磁共振。全身RF线圈和/或局部线圈根据共振衰变接收RF信号作为磁共振数据，所述磁共振数据被传送到RF接收器14。线圈元件或阵列可以单独地或可选择地发送或接收组以提供所选择的发送和/或接收减小因子。
系统1包括序列控制器16，其控制对诸如T1扫描、T2*扫描等的成像扫描、预扫描和/或线圈调查扫描的操作。线圈调查扫描或预扫描可以包括例如梯度回波或快速场回波(FFE)序列。线圈调查扫描或预扫描使用k空间的笛卡尔采样或六边形采样。在一个实施例中，采样可以包括对k空间的压缩平行采样或感测采样。序列控制器16控制RF发射器单元18和梯度控制器20，RF发射器单元18控制对全身RF线圈10和/或局部RF线圈12的操作，梯度控制器20控制对梯度线圈8的操作。RF发射器单元18被布置为独立地驱动线圈元件或阵列，例如利用多个发射器。控制单元和对应线圈之间的通信可以是无线的或有线的。
系统1包括用户接口22、智能选择单元24、信噪比确定单元26和平行成像参数单元28。用户接口被适当地实施在工作站36中，被配置为与医疗保健从业者或扫描器2的操作者进行交互，诸如接收总体减小因子、最小SNR、最大扫描时间等，显示加速因子以便操作者审查等。工作站36包括电子处理器或电子处理设备30、显示设备32以及至少一个输入设备34。显示设备32显示菜单、面板和用户控制，而输入设备34接收来自医疗保健医师的输入，包括总体减小因子或加速因子或最小SNR。输入设备也接收命令以开始预扫描或线圈调查扫描。工作站36可以是桌上计算机、膝上计算机、平板电脑、移动计算设备、智能手机等。显示设备30可以包括计算机监视器、触摸屏、阴极射线管(CRT)、存储管、平面直角显示器、发光二极管(LED)显示器、电致发光(ELD)显示器、等离子显示板(PDP)、液晶显示器(LCD)、有机发光二极管(OLED)、投影仪等。输入设备32可以是键盘、鼠标、麦克风等。
智能选择单元24从RF接收器14接收线圈调查数据并且将调查数据存储在存储器38中。存储器可以是诸如随机存取存储器(RAM)的计算机存储器，或者诸如硬盘、盘阵列、固态存储器等盘存储器。智能选择单元生成满足接收到的总体减小因子和/或最小SNR的多个减小因子。减小因子包括相位编码。例如，总体减小因子4可以包括减小因子(4x1)、(3.5x8/7)、(2.5x8/5)、(2x2)、(1.5x8/3)、(1x4)等。针对被编码方向的减小因子的乘积近似于总体减小因子的值。这些值是基于粗网格的离散的值，其可以在迭代中被进一步细化。智能选择单元构建信号图40和多个噪声图42。
SNR确定单元计算针对不同方向的不同平行减小因子的SNR。SNR通过信号图除以噪声图中的一个来确定，其提供针对每个体素的SNR值。SNR由体积的平均来确定，例如其中N是体素的数量，而snr(r)是第r个体素的信噪比(SNR)。预期计算平均SNR或从噪声图获得标准的其他方法。在一个实施例中，该平均可以仅包括感兴趣区域中的那些体素。
平行成像参数单元28选择满足具有最高可接受SNR的接收到的总体减小因子的减小因子和方向。减小因子是针对被编码的方向用于构建具有最高SNR的噪声图。平行成像参数单元可以在迭代中将所选择的减小因子提供给智能选择单元利用较小的步骤或所选择的减小因子周围的较小的增量离散值以进行进一步的细化。在另一实施例中，减小因子在离散步骤中在编码方向上展开，并且以较少步骤在最高SNR周围进一步细化。
用户接口可以将SNR值显示为对比度图像、针对每个减小因子集合的平均SNR值和/或最高SNR。用户接口提供所选择的减小因子到序列控制以基于所选择的减小因子和/或编码方向来执行对对象的磁共振成像扫描。扫描可以包括对k空间的笛卡尔采样、六边形采样或压缩感测采样。
各种单元或模块22、24、26和28由诸如工作站36的电子处理器或电子处理设备30的电子数据处理设备，或者由能通过网络与工作站36连接的基于网络的服务器计算机等适当地实现。此外，公开的用户接口、信号和噪声图构建以及参数选择技术使用执行所公开的技术的存储指令(例如，软件)的非暂时性存储介质来适当地实施，所述指令能够由电子数据处理设备读取并且能够电子数据处理设备执行。
参考图2，用流程图表示了智能选择处理的一个实施例。简要概括，自动地计算在每个方向的SENSE减小因子。预扫描确定每个线圈元件或线圈元件的组的噪声量。利用不同的SENSE减小因子针对不同方向计算针对每个体素的SNR。操作者可以指定总体减小因子，而系统告知操作者针对所选择的SENSE因子的SNR惩罚。备选地，操作者指定最大可接受SNR并且系统选择减小因子。更详细地查看该方法，步骤或单元44，执行预扫描或线圈调查扫描，其接收来自N个数据通道的数据和Q身体数据。每个通道对应于一个(或多个)线圈元件。噪声协方差矩阵Ψ根据接收到的数据而被构建。
步骤或单元46在接收到的数据上执行多平面重构，其将数据重构成具有临床扫描的几何配置的斜片。步骤或单元48计算线圈敏感度和调控数据以确定由感兴趣区域的患者所占据的成像体的部分中的敏感度。
信号图和噪声图由步骤或单元50根据线圈敏感度和调控数据来构建。信号图针对每个减小因子的集合是相同的，而噪声图基于减小因子而不同。利用不同的减小因子重复构建噪声图以识别最佳SNR的方向。在一个实施例中，使用矩阵Η＝(SΗΨ-1S+R-1)-1SΗΨ-1来计算噪声图，其中，S是线圈敏感度图，R是由患者或感兴趣区域所占据的图像体积中的最大预期MR信号的图，而Ψ是协方差矩阵。噪声(r)被计算为diag(ΗΨΗΗ)。
通过步骤或单元52将信号图除以噪声图，其创建SNR图。SNR图包括针对每个体素的SNR值。步骤或单元54计算SNR图的平均值或其他标准，诸如其中N是体素的数量，而snr(r)是第r个体素的信噪比(SNR)。与噪声图相对应的所有体素的平均SNR被用于确定每个噪声图的最高SNR，如参考图1所讨论的那样。基于该信息，包括按每个方向的减小因子的减小因子被确定，并且操作者可以选择在图像质量(SNR)和成像速度(减小因子)之间的最佳平衡。
图3用流程图表示了对平行加速参数进行自动优化的方法，例如沿着减小方向的加速因子，其可以由一个或多个处理器或离散的单元或模块来实施。在步骤56，总体减小因子或SNR最小值由用户接口接收。在步骤58，执行诸如快速场回波(FFE)扫描的线圈调查扫描。线圈调查可以利用对k空间的笛卡尔采样、六边形采样或压缩感测采样来执行。备选地，执行先前扫描。
在步骤59，生成满足接收到的总体减小因子的多个减小因子。总体减小因子可以基于最小SNR。例如，利用总体减小因子4，集合(4，1)、(3.5，8/7)、(3，4/3)等满足总体减小因子4。减小因子加速扫描时间或减小对k空间的采样。减小因子包括针对两个编码方向的减小因子。在一个实施例中，所生成的减小因子的集合包括粗网格或离散值，其可以基于粗值的初始最高SNR以较小的增量被进一步细化。
在步骤60构建信号图和多个噪声图。该步骤包括参考图2所描述的处理。在步骤60针对所生成的减小因子的集合和/或基于最小SNR，来构建信号图和噪声图。继续如参考图2所描述的智能选择处理，在步骤62中基于减小因子的集合针对每个噪声图的每个体素来确定SNR。通过信号图处理噪声图来计算SNR。还计算在感兴趣对象或整个成像体积上的平均SNR。
在决定步骤64，迭代被执行，其重复对噪声图的构建，以及关于初始最高SNR和相对应的减小因子利用减小因子的较小增量的SNR确定。
在步骤66基于对应的噪声图的最高SNR来选择减小因子。对小于总体减小因子的减小因子的搜索沿着减小方向。在步骤68，利用所选择的减小因子来执行磁共振成像扫描。可以例如利用对使用所选择的减小因子压缩的k空间的SENSE采样来执行该扫描。
根据本发明的某种示例性实施例，示例性MR系统可以通过对误差传播的快速评估、本质上执行相对低分辨率的G因子评估，针对在开始可以提供的给定成像几何配置、对象范围、接收线圈配置和/或其他信息来优化依赖于SENSE因子的SENSE性能和采样方法方案。SENSE扫描参数采用参数空间的形式。例如，这些参数可以经由检查参数呈现给用户。基于可提供的信息，示例性MR系统可以确定最优(或最好的)参数集合，或者其可以通过将特定SENSE参数改变为某些有用的参数来满足一些预定义的检查参数特异的约束(例如，SNR损失、G因子分布/统计、特别感兴趣区域(ROI)等)。还可能的是，通知用户关于他/她的选择的结果，并且提供解决问题的选项。
通过最佳敏感度编码过程和对展开问题(G因子(1))的良好的适应来获得最终图像质量。在减小的FOV中单个位置的SENSE问题可以被认为成例如Sx＝y。
在该实例中，落在彼此之上(由于均一的k空间欠采样的折叠)来自不同位置的N个未知体素信号贡献可以被概括为向量x。接收到的信号被敏感度编码，并且利用M个接收线圈被测量(M≥N)。其对应的信号被组合在向量y。矩阵S表示(NxM)的敏感度矩阵。所选择的减小因子R(在2DSENSE的情况下Rtotal＝R1+R2)和所选择的k空间特定子采样方案定义在测量过程期间被适当地进行敏感度编码r全FOV中的那些体素将被折叠在彼此之上。这意味着Sx＝y的实际形式取决于所选择的经验参数，其应当被概括在向量Δ中。
对Sx＝y中将该矩阵求逆得出向量x，其包含对应的N个体素特定异的信号。例如，〖〖〖(S〗^HS)〗^(-1)S〗^Hy＝x。
矩阵〖〖〖(S〗^HS)〗^(-1)S〗^H是S的伪逆，并且其范数描述了从测量到最终图像的误差传播。一个潜在的范数g因子可以给出针对好的数值条件的标准化测量。例如，g＝√([〖(S^HS)〗^(-1))〖(S^HS)〗^]≥1。
g因子取决于矩阵S，矩阵S取决于所选择的采样参数Δ。因此，g因子可以被用于判断所选择的SENSE测试参数Δ是否最优或者是否应该差别地选择他们以做的更好。为此，人们不得不在所有可能的设置Δ中将g最小化。然而，g＝√([〖(S^HS)〗^(-1))〖(S^HS)〗^]≥1描述了在约减的FOV中的仅一个体素，并且不得不找到在约减的FOV(在2D成像情况下：矩阵LxK，在3D更高维度)中几乎所有体素的折衷。这可能是数值上的挑战并且评估可能被隐藏在背景中，但是可以想象针对加速的若干度量。
如果不是在rFOV中的所有LxK个体素都被用于优化Δ参数，则该问题的尺度可以减小。为了这个目的，约减的FOV可以不统一地或者随机地“欠采样地”(低分辨率方法)约减，因此g因子元素的数量将被考虑。这也意味着线圈敏感度是空间中相当平滑的函数。另外，该评估还可以被集中在具有特别的临床感兴趣(ROI)的图像中的重要区域。
减小问题的尺度的第二种方法是变换在虚拟线圈系统中所有涉及的线圈(使用SVD或类似的方法)。因此，只有N个显著线圈在该评估中必须被考虑。
然而，使用所有线圈直接计算针对每个位置的S^HS，作为伪逆的基本计算的元素，也是有效率的，并且避免忽略线圈信息的部分。
或者可以采用更简单的范数来代替g＝√([〖(S^HS)〗^(-1))〖(S^HS)〗^]≥1以用于评估。
可以应用统计过程来分析在减小的FOV中的g因子分布以允许比较不同的参数集Δ。此外，对最优Δ的搜索可以以分层的方式执行，从Δ空间和rFOV空间中的粗分辨率开始，其可以在优化过程期间被调节/增加。
可以被典型地应用在SENSE展开中的其他信息(噪声相关矩阵、调控图像，或在最简单情况下将线圈敏感度限制于所期望的信号来自的区域的二值掩码)也可以被包括在优化中的g因子的计算中。
此外，关于编码问题的结构的信息可以被并入。这能够在例如k-tSENSE应用中特别有趣，其中训练数据集可以用于展开。如果该信息在测量之前已经可用，或者作为粗略估计(例如，来自库)或者来自测量的模型，则这可以用作到适当的k-t采样方案的选择中的另外的输入以及专用于该患者或成像问题的减小因子。
例如，使用给定阵列线圈执行具有固定加速因子(R总体＝4)3DFFE扫描。例如，这可能由于临床需要而保证扫描短于20秒。用户定义扫描几何配置和读出方向以及最终成像矩阵和R总体。被给予这些约束和SENSE参考扫描信息，扫描器确定应当使用哪些潜在的欠采样因子R1和R2以及哪个采样方案。在没有系统帮助的情况下做出该决定是复杂的。相比之下，根据本发明的示例性实施例，可能的是由系统提供该决定，增加用户友好和简单性的同时帮助改善图像质量。
将会意识到，结合本文所呈现的特定示例性实施例将某种结构和/或功能特征描述为被并入定义的元件和/或组件。然而，可以预期的是，为了相同或相似的益处，在合适时这些特征也可以同样地被并入其他元件和/或组件。也将意识到，示例性实施例的不同的方面可以合适地选择性地被采用于以完成适于所期望的应用的其他替代实施例，其他替代实施例通过其来实现并入其中的方面的各个优点。
还将意识到，本文所描述的特定元件或组件可以经由硬件、软件、固件或其组合适当地实施其功能性。另外，将会意识到，本文所描述的合并在一起的特定元件可以在适当的环境下成为独立的元件或以其他方式分离。类似地，描述为由一个特定元件实施的多个特定功能可以由独立动作的多个不同的元件来实施以实现各个功能，或者某些单独的功能可以被分离并且由合作动作的多个不同的元件来实施。备选地，在合适时彼此不同的本文中以其他方式描述和/或示出的一些元件或部件可以在物理上或功能上组合。
简言之，参考优选的实施例阐述了本说明。显然，他人在阅读和理解本说明之后将进行修改和变动。目的是，本发明被理解为包括所有这些修改和变化，只要这些修改和变化落入所附权利要求或其等价方案的范围内。也就是说，将会意识到，上文所公开的或其他特征和功能或其替代中的各种，可以按期望组合到很多其他不同的系统或应用中，并且同样可以随后由本领域技术人员进行各种目前未预见或意料之外的其中的替代、修改、变型或改进，同样旨在由权利要求书所包含。

Claims (20)

1.一种平行磁共振成像系统(1)，包括：
具有多个线圈元件的至少一个射频(RF)线圈(10，12)；
智能选择单元(24)，其根据利用所述至少一个RF线圈的对对象的预扫描或先前扫描，构建(60)信号图和基于沿着一个或多个减小方向的不同的小减因子的集合的多个噪声图；
平行成像参数单元(28)，其选择与对信噪比(SNR)和扫描速度进行优化的噪声图相对应的减小因子的集合；
序列控件(16)，其基于所选择的减小因子来执行对所述对象的磁共振成像扫描。
2.根据权利要求1所述的系统，其中，所述减小因子的集合包括针对两个减小方向的两个减小因子。
3.根据权利要求1和2中的任一项所述的系统，还包括：
用户接口(22)，其接收最大成像时间或总体减小因子的其中一个，并且
其中，所述智能选择单元(24)基于所述最大成像时间或所述总体减小因子来生成多个减小因子的集合。
4.根据权利要求1-3中的任一项所述的系统，其中，所述减小因子的集合包括多个减小因子，每个减小因子与被编码的方向相对应。
5.根据权利要求求1和2中的任一项所述的系统，还包括：
用户接口(22)，其接收最小SNR，并且
其中，所述智能选择单元(24)生成满足接收到的最小SNR的总体减小因子。
6.根据权利要求5所述的系统，其中，所述平行成像参数单元(28)选择所述总体减小因子，所述总体减小因子使生成具有至少接收到的最小SNR的图像的成像时间最小化。
7.根据权利要求1-6中的任一项所述的系统，其中，对减小因子的所述选择包括沿着每个编码方向的搜索。
8.根据权利要求1-7中的任一项所述的系统，其中，所述智能选择单元(24)被配置为：
根据所述预扫描或所述先前扫描来构建(44)3维(3D)数据的集合；
执行(46)对所述3D数据的多平面重构以确定具有所述成像扫描几何配置的通道和Q身体数据；
基于所确定的通道数据和Q身体数据来计算(48)线圈敏感度和调控数据；以及
基于所确定的线圈敏感度和调控数据来生成(50)信号图和噪声图。
9.根据权利要求1所述的系统，其中，最高SNR包括全部体素的平均SNR。
10.一种平行磁共振成像的方法，包括：
基于从对对象的先前扫描接收到的总体减小因子来构建(60)信号图和多个噪声图，每个噪声图基于沿着一个或多个减小方向的不同的减小因子的集合被确定；
选择(66)与满足信噪比(SNR)标准和扫描速度标准的噪声图相对应的减小因子；
基于所选择的减小因子来执行对所述对象的磁共振成像扫描。
11.根据权利要求10所述的方法，其中，选择所述减小因子包括选择两个减小因子，每个减小因子针对一个减小方向。
12.根据权利要求10和11中的任一项所述的方法，还包括：
接收最大成像时间或总体减小因子中的一个；
基于所述最大成像时间或所述总体减小因子来生成多个减小方向的集合。
13.根据权利要求10-12中的任一项所述的方法，其中，所述减小因子的每个集合包括多个减小因子，每个减小因子与编码方向相对应。
14.根据权利要求10和11中的任一项所述的方法，还包括：
接收最小SNR标准，以及
生成满足接收到的最小SNR标准的总体减小因子。
15.根据权利要求14所述的方法，还包括：
选择所述减小因子和/或减小方向，所述减小因子和/或减小方向将成像时间最小化以生成至少满足接收到的最小SNR标准的图像。
16.根据权利要求10-15中的任一项所述的方法，还包括：
以离散的步骤沿着每个编码方向搜索所述减小因子；
以更小的步骤在最高SNR周围细化所述减小因子。
17.根据权利要求10-16中的任一项所述的方法，还包括：
根据所述先前扫描来构建(44)3维数据的集合；
执行(46)对所述3D数据的多平面重构以确定具有所述成像扫描几何配置的通道和Q身体数据；
基于所确定的通道数据和Q身体数据来计算(48)线圈敏感度和调控数据；并且
基于所确定的线圈敏感度和调控数据来生成(50)信号图和噪声图。
18.一种承载软件的非暂时性计算机可读存储介质，所述软件控制一个或多个电子数据处理设备(30)以执行根据权利要求10-17中的任一项所述的方法。
19.一种电子数据处理设备(30)，其被配置为执行根据权利要求10-17中的任一项所述的方法。
20.一种平行磁共振成像系统(1)，包括：
具有多个线圈元件的至少一个RF线圈(10，12)；
用户接口(22)，其接收最小图像质量和/或生成图像数据的最大时间中的至少一个；
一个或多个处理器(30)，被配置为：
根据利用所述至少一个RF线圈对对象的先前扫描，构建(60)信号图和针对沿着一个或多个减小方向的多个减小因子中的每个的多个噪声图；
选择(66)减小因子，所述减小因子生成具有至少接收到的最小图像质量和/或在接收到的最大时间内的图像数据；
令所述成像系统基于所选择的减小因子来执行(68)对所述对象的磁共振成像扫描。
CN201480014017.0A 2013-03-13 2014-03-12 平行成像加速参数的自动优化 Active CN105308473B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201361778582P 2013-03-13 2013-03-13
US61/778,582 2013-03-13
PCT/IB2014/059667 WO2014141086A1 (en) 2013-03-13 2014-03-12 Automatic optimization of parallel imaging acceleration parameters

Publications (2)

Publication Number Publication Date
CN105308473A true CN105308473A (zh) 2016-02-03
CN105308473B CN105308473B (zh) 2018-11-27

Family

ID=50440712

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480014017.0A Active CN105308473B (zh) 2013-03-13 2014-03-12 平行成像加速参数的自动优化

Country Status (5)

Country Link
US (1) US10107884B2 (zh)
EP (1) EP2972450B1 (zh)
JP (1) JP6430417B2 (zh)
CN (1) CN105308473B (zh)
WO (1) WO2014141086A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111830449A (zh) * 2019-04-19 2020-10-27 西门子(深圳)磁共振有限公司 呼吸信号提取方法及装置、磁共振成像系统、存储介质

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012023098A1 (en) * 2010-08-20 2012-02-23 Koninklijke Philips Electronics N.V. Virtual coil emulation in parallel transmission mri
US10429475B2 (en) * 2013-03-12 2019-10-01 The General Hospital Corporation Method for increasing signal-to-noise ratio in magnetic resonance imaging using per-voxel noise covariance regularization
US10955508B2 (en) * 2016-08-25 2021-03-23 Koninklijke Philips N.V. BO-corrected sensitivity encoding magnetic resonance imaging
EP3470868A1 (en) 2017-10-16 2019-04-17 Koninklijke Philips N.V. Quantitative measurement of relaxation times in magnetic resonance imaging
US11023785B2 (en) * 2018-07-23 2021-06-01 International Business Machines Corporation Sparse MRI data collection and classification using machine learning
EP3910358A1 (en) * 2020-05-15 2021-11-17 Koninklijke Philips N.V. Automated adjustment of the undersampling factor
WO2022032133A1 (en) * 2020-08-06 2022-02-10 University Of Pittsburgh - Of The Commonwealth System Of Higher Education Optimizing magnetic resonance imaging protocols for increased speed
CN112630709B (zh) * 2020-11-03 2024-08-27 杭州微影医疗科技有限公司 金属植入物磁共振成像方法、系统、终端及存储介质
US20220343564A1 (en) * 2021-04-27 2022-10-27 Siemens Healthcare Gmbh Method for automated regularization of hybrid k-space combination using a noise adjustment scan

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020153890A1 (en) * 2001-04-20 2002-10-24 Bruno Madore Combining unfold with parallel magnetic resonance imaging
CN1468583A (zh) * 2002-06-07 2004-01-21 Ge医疗系统环球技术有限公司 线圈灵敏图生成方法和平行成像方法以及磁共振成像装置
CN1572246A (zh) * 2003-05-08 2005-02-02 株式会社东芝 基于扩大展开法的平行成像
US20050083051A1 (en) * 2003-08-25 2005-04-21 Swen Campagna Method for determining an acceleration factor of a parallel image acquisition in magnetic resonance imaging
CN1930485A (zh) * 2004-03-12 2007-03-14 皇家飞利浦电子股份有限公司 用于优化mri扫描参数的预扫描
US20080068016A1 (en) * 2006-09-15 2008-03-20 Ajeetkumar Gaddipati System and method of accelerated mr propeller imaging
US7449888B1 (en) * 2005-07-27 2008-11-11 General Electric Company Method and apparatus for multi-dimensional parallel MR imaging
US20120092009A1 (en) * 2010-10-19 2012-04-19 The Board Of Trustees Of The Leland Stanford Junior University Autocalibrating parallel imaging reconstruction method from arbitrary k-space sampling with reduced noise

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
BR9810032A (pt) * 1998-04-17 2000-09-19 Koninkl Philips Electronics Nv Processo e aparelho de ressonância magnética para obter imagens por intermédio de ressonância magnética
US6841998B1 (en) 2001-04-06 2005-01-11 Mark Griswold Magnetic resonance imaging method and apparatus employing partial parallel acquisition, wherein each coil produces a complete k-space datasheet
JP2005525185A (ja) * 2002-05-13 2005-08-25 コーニンクレッカ　フィリップス　エレクトロニクス　エヌ　ヴィ 磁場強度の関数としてのパラレルイメージングにおけるリダクションファクタの固有の限界
CN1977181A (zh) 2004-06-28 2007-06-06 皇家飞利浦电子股份有限公司 平行磁共振成像
US7429862B2 (en) 2006-04-18 2008-09-30 Siemens Aktiengesellschaft PPA MRT method and apparatus with scaled contribution of k-space data to the reconstructed image
JP4464373B2 (ja) 2006-07-12 2010-05-19 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｍｒｉ装置
US7394252B1 (en) 2007-05-03 2008-07-01 The General Hospital Corporation Regularized GRAPPA reconstruction
US8354844B2 (en) 2009-12-23 2013-01-15 Universitaetsklinikum Freiburg Method for data acquisition acceleration in magnetic resonance imaging (MRI) with N-dimensional spatial encoding using two or more receiver coil arrays and non-linear phase distributions
US10551461B2 (en) * 2010-09-01 2020-02-04 Commissariat A L'energie Atomique Et Aux Energies Alternatives Method for performing parallel magnetic resonance imaging
JP5926285B2 (ja) * 2010-12-22 2016-05-25 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 任意ｋ空間トラジェクトリの高速並行再構成

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020153890A1 (en) * 2001-04-20 2002-10-24 Bruno Madore Combining unfold with parallel magnetic resonance imaging
CN1468583A (zh) * 2002-06-07 2004-01-21 Ge医疗系统环球技术有限公司 线圈灵敏图生成方法和平行成像方法以及磁共振成像装置
CN1572246A (zh) * 2003-05-08 2005-02-02 株式会社东芝 基于扩大展开法的平行成像
US20050083051A1 (en) * 2003-08-25 2005-04-21 Swen Campagna Method for determining an acceleration factor of a parallel image acquisition in magnetic resonance imaging
CN1930485A (zh) * 2004-03-12 2007-03-14 皇家飞利浦电子股份有限公司 用于优化mri扫描参数的预扫描
US7449888B1 (en) * 2005-07-27 2008-11-11 General Electric Company Method and apparatus for multi-dimensional parallel MR imaging
US20080068016A1 (en) * 2006-09-15 2008-03-20 Ajeetkumar Gaddipati System and method of accelerated mr propeller imaging
US20120092009A1 (en) * 2010-10-19 2012-04-19 The Board Of Trustees Of The Leland Stanford Junior University Autocalibrating parallel imaging reconstruction method from arbitrary k-space sampling with reduced noise

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
PAUL T. WEAVERS ET AL: "Acceleration Apportionment: A Method of Improved 2D SENSE Acceleration Applied to 3D Contrast-Enhanced MR Angiography", 《PAUL T. WEAVERS ET AL》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111830449A (zh) * 2019-04-19 2020-10-27 西门子(深圳)磁共振有限公司 呼吸信号提取方法及装置、磁共振成像系统、存储介质
CN111830449B (zh) * 2019-04-19 2022-10-18 西门子(深圳)磁共振有限公司 呼吸信号提取方法及装置、磁共振成像系统、存储介质

Also Published As

Publication number Publication date
CN105308473B (zh) 2018-11-27
EP2972450A1 (en) 2016-01-20
WO2014141086A1 (en) 2014-09-18
US10107884B2 (en) 2018-10-23
JP2016509940A (ja) 2016-04-04
US20160018498A1 (en) 2016-01-21
JP6430417B2 (ja) 2018-11-28
EP2972450B1 (en) 2021-02-17

Similar Documents

Publication Publication Date Title
CN105308473A (zh) 平行成像加速参数的自动优化
CN104781685B (zh) 用于使用非相干采样和冗余haar子波的动态磁共振成像的图像重建
CN111513716B (zh) 使用扩展灵敏度模型和深度神经网络进行磁共振图像重建的方法和系统
CN102159965B (zh) 用于mri的b1映射和b1l匀场
CN101971045B (zh) 用于并行磁共振成像的线圈选择
US11762049B2 (en) Medical data processing apparatus, medical data processing method, and magnetic resonance imaging apparatus
CN104765011B (zh) 磁共振原始数据的重建方法和装置及磁共振系统
CN104101849A (zh) 磁共振系统控制序列的确定
CN103154761B (zh) 在并行发射mri中的虚拟线圈仿真
US20110241670A1 (en) Accelerated dynamic magnetic resonance imaging system and method
CN103380384B (zh) 使用交错的水参考扫描的具有自动相位和b0校正的磁共振波谱成像方法和装置
CN108957375A (zh) 磁共振成像方法、装置、设备及存储介质
CN103323801A (zh) 磁共振系统控制序列的确定
CN103477238A (zh) 采用来自先验采集的约束的压缩感测mr图像重建
CN102007423A (zh) 用于患者安全性和经改进的扫描性能的实时局部和全身sar估计
CN103027681B (zh) 用于重构并行获取的mri图像的系统
US11815581B2 (en) Medical information processing apparatus, magnetic resonance imaging apparatus, and medical information processing method
CN105051564A (zh) 用于压缩感测MRI中的经改进的k-空间采样的方法
US9971007B2 (en) Method and apparatus for accelerated magnetic resonance imaging
US20100264923A1 (en) System for Improved MR Image Reconstruction
CN106793971A (zh) 磁共振成像装置
RU2764643C2 (ru) Управляемая потоком данных коррекция фазозависимых артефактов в системе магнитно-резонансной томографии
US7956611B2 (en) Magnetic resonance imaging apparatus and method
CN110542872A (zh) 磁共振成像方法及设备、相位信息的获取方法及装置
US9459335B2 (en) System and method for parallel magnetic resonance imaging with optimally selected in-plane acceleration

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