CN1265762C - 线圈灵敏图生成方法和平行成像方法以及磁共振成像装置 - Google Patents

线圈灵敏图生成方法和平行成像方法以及磁共振成像装置 Download PDF

Info

Publication number
CN1265762C
CN1265762C CNB031386695A CN03138669A CN1265762C CN 1265762 C CN1265762 C CN 1265762C CN B031386695 A CNB031386695 A CN B031386695A CN 03138669 A CN03138669 A CN 03138669A CN 1265762 C CN1265762 C CN 1265762C
Authority
CN
China
Prior art keywords
coil
sensitivity map
sensitivity
phased array
pixel
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.)
Expired - Fee Related
Application number
CNB031386695A
Other languages
English (en)
Other versions
CN1468583A (zh
Inventor
池崎吉和
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.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Publication of CN1468583A publication Critical patent/CN1468583A/zh
Application granted granted Critical
Publication of CN1265762C publication Critical patent/CN1265762C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • 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
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Signal Processing (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

为了提高利用一般称为SENSE的平行成像方法进行自校准所得到的图像的均匀性，利用通过对相控阵列线圈线圈n(n＝1-N，N≥2)进行整个FOV校准扫描得到的复合图像Cn的求和图像Ab＝∑|Cn|，计算灵敏度因数Sn(p)(步骤V2，V3)。此外，利用加权最小二乘法对绝对值求和图像Ab的像素值Ab(p)的平方进行曲线拟合，生成灵敏度图Sn(步骤V4)。

Description

线圈灵敏图生成方法和平行成像方法以及磁共振成像装置
发明领域
本发明涉及线圈灵敏度图生成方法、平行成像方法以及磁共振成像(MRI)装置，更具体讲，涉及能够利用一般称为SENSE(灵敏度编码)的平行成像方法提高在自校准中得到的图像的均匀性的线圈灵敏图生成方法、平行成像方法以及MRI装置。
背景技术
例如，在下列文献中，对涉及SENSE的常规技术进行了描述：
(1)Pruessmann，K.P.，Weiger，M.，Scheidegger，M.B.，和Boesiger P.，Magn.Reson.Med.，42，952(1999)；和
(2)P.B.Roemer，W.A.Edelstein，C.E.Hayes，S.P.Souza，和O.M.Mueller，Magn.Reson.Med.，16，192-225(1990)。
图5为表示常规的MRI装置中进行平行成像处理操作的流程图。
在步骤J1，利用相控阵列线圈线圈n(n＝1-N，N≥2)对整个FOV(视场)进行校准扫描，以得到与相控阵列线圈线圈n相应的复合图像Cn。
在步骤J2，生成每个符合图像Cn的平方图像Cn 2，并且生成平方图像Cn 2的求和平方根图像As＝sqrt{∑Cn 2}。“sqrt{}”为求平方根的函数。
在步骤J3，利用以下公式对相控阵列线圈线圈n中的像素p的灵敏度因数Sn(p)′进行计算：
Sn(p)′＝Cn(p)/As(p)，
式中，将复合图像Cn中的像素p的像素值表示为Cn(p)，将求和平方根图像As中的像素p的像素值表示为As(p)。
在步骤J4，根据灵敏度因数Sn(p)′生成灵敏度图Sn′。
在步骤J5，利用相控阵列线圈线圈1-线圈N进行较小FOV的扫描，以得到相控阵列线圈线圈1-线圈N的相应的复合图像A1-AN。
在步骤J6，根据灵敏度图Sn′和复合图像A1-AN生成整个FOV的复合图像。
在步骤J7，重复进行步骤J5和J6，以得到理想数量的图像。
由于图5所示的常规的平行成像处理只利用相控阵列线圈线圈n(“自校准”)生成灵敏度图Sn′，因此不需要利用体线圈获得单独的校准图像，由此具有使操作简化的优点。
然而，尽管利用复合图像Cn的平方图像Cn 2的求和平方根图像As＝sqrt{∑Cn 2}计算的灵敏度因数Sn(p)′提供了良好的S/N(信噪)比，但是，就图像均匀性来说，它们不是最佳的。即，通过常规的自校准得到的图像具有均匀性极差的缺点。
发明概述
因此，本发明的一个目的是提供一种能够提高自校准中得到的图像的均匀性的线圈灵敏图生成方法、平行成像方法以及MRI装置。
按照其第一方面，本发明提供了一种线圈灵敏图生成方法，其特征在于包括：当将通过对每个相控阵列线圈线圈n(n＝1-N，N≥2)进行整个FOV校准扫描而得到的复合图像Cn中的像素p的像素值表示为Cn(p)，并且将通过将根据所述复合图像Cn生成的绝对值图像|Cn|相加而得到的求和图像Ab＝∑|Cn|中的像素p的像素值表示为Ab(p)时，利用下式计算相控阵列线圈线圈n中的像素p的灵敏度因数Sn(p)
Sn(p)＝Cn(p)/Ab(p)；
并且根据所述灵敏度因数Sn(p)，生成相控阵列线圈线圈n的相应的灵敏度图Sn。
第一方面的线圈灵敏度图生成方法利用复合图像Cn的绝对值图像|Cn|的求和图像Ab＝∑|Cn|计算灵敏度因数Sn(p)，因此，与利用复合图像Cn的平方图像Cn 2的求和平方根图像As＝sqrt{∑Cn 2}相比，提高了图像的均匀性。
按照其第二方面，本发明提供了具有前述结构的线圈灵敏度图生成方法，其特征在于包括：通过对所述灵敏度因数Sn(p)进行曲线拟合生成灵敏度图Sn。
利用复合图像Cn的绝对值图像|Cn|的求和图像Ab＝∑|Cn|计算的灵敏度因数Sn(P)具有比利用复合图像Cn的平方图像Cn 2的求和平方根图像As＝sqrt{∑Cn 2}更差的S/N比。
但是，第二方面的线圈灵敏度图生成方法通过对灵敏度因数Sn(p)进行曲线拟合生成灵敏度图Sn，因此，能够抑制噪声的影响，从而提高S/N比。
按照其第三方面，本发明提供了具有前述结构的线圈灵敏度图生成方法，其特征在于包括：利用二维多项式进行二维曲线拟合。
第三方面的线圈灵敏度图生成方法利用二维多项式进行二维曲线拟合，因此，能够生成具有平滑表面的灵敏度图。
按照其第四方面，本发明提供了具有前述结构的线圈灵敏度图生成方法，其特征在于包括：利用二维二次多项式a·x2+b·y2+c·x+d·y+e·x·y+f进行二维曲线拟合。
第四方面的线圈灵敏度图生成方法利用二维二次多项式进行二维曲线拟合，因此，能够生成具有平滑表面的灵敏度图。
按照其第五方面，本发明提供了具有前述结构的线圈灵敏度图生成方法，其特征在于包括：利用由所述绝对值求和图像Ab的像素值Ab(p)的平方加权的最小二乘法进行曲线拟合。
绝对值求和图像Ab的像素值Ab(p)的平方的值在信号部分中较大而在噪声部分中较小。
因此，第五方面的线圈灵敏度图生成方法利用由所述绝对值求和图像Ab的像素值Ab(p)的平方加权的最小二乘法进行曲线拟合。由此，可以进一步抑制噪声的影响，以进一步提高S/N比。
按照其第六方面，本发明提供了一种平行成像方法，其特征在于包括：通过对相控阵列线圈线圈n(n＝1-N，N≥2)进行整个FOV校准扫描得到相应的复合图像Cn；根据每个所述复合图像Cn生成绝对值图像|Cn|；通过将所述绝对值图像|Cn|相加生成求和图像Ab＝∑|Cn|；利用下式计算在相控阵列线圈线圈n中的像素p的灵敏度因数Sn(p)：
Sn(p)＝Cn(p)/Ab(p)，
式中，将所述复合图像Cn中的像素p的像素值表示为Cn(p)，并且将所述求和图像Ab中的像素p的像素值表示为Ab(p)；根据所述灵敏度因数Sn(p)生成相控阵列线圈线圈n的相应的灵敏度图Sn；通过对相控阵列线圈线圈n进行具有较小FOV的扫描得到相应的复数图An；并且根据所述灵敏度图Sn以及所述复数图An，生成整个FOV的复合图像。
第六方面的平行成像方法利用复合图像Cn的绝对值图像|Cn|的求和图像Ab＝∑|Cn|计算灵敏度因数Sn(p)，因此，与利用复合图像Cn的平方图像Cn 2的求和平方根图像As＝sqrt{∑Cn 2}相比，提高了图像的均匀性。
按照其第七方面，本发明提供了具有前述结构的平行成像方法，其特征在于包括：通过对所述灵敏度因数Sn(p)进行曲线拟合生成灵敏度图Sn。
利用复合图像Cn的绝对值图像|Cn|的求和图像Ab＝∑|Cn|计算的灵敏度因数Sn(p)具有比利用复合图像Cn的平方图像Cn 2的求和平方根图像As＝sqrt{∑Cn 2}差的S/N比。
但是，第七方面的平行成像方法通过对灵敏度因数Sn(p)进行曲线拟合生成灵敏度图Sn，因此，能够抑制噪声的影响，从而提高S/N比。
按照其第八方面，本发明提供了具有前述结构的平行成像方法，其特征在于包括：利用二维多项式进行二维曲线拟合。
第八方面的平行成像方法利用二维多项式进行二维曲线拟合，因此，能够生成具有平滑表面的灵敏度图。
按照其第九方面，本发明提供了具有前述结构的平行成像方法，其特征在于包括：利用二维二次多项式a·x2+b·y2+c·x+d·y+e·x·y+f进行二维曲线拟合。
第九方面的平行成像方法利用二维二次多项式进行二维曲线拟合，因此，能够生成具有平滑表面的灵敏度图。
按照其第十方面，本发明提供了具有前述结构的平行成像方法，其特征在于包括：利用由所述像素值Ab(p)的平方加权的最小二乘法进行曲线拟合。
绝对值求和图像Ab的像素值Ab(p)的平方的值在信号部分中较大而在噪声部分中较小。
因此，第十方面的平行成像方法利用由绝对值求和图像Ab的像素值Ab(p)的平方加权的最小二乘法进行曲线拟合。由此，可以进一步抑制噪声的影响，以进一步提高S/N比。
按照其第十一方面，本发明提供了一种MRI装置，其特征在于包括：一个用于发射RF(射频)脉冲的发射线圈；一个用于施加梯度磁场的梯度线圈；用于接收MR信号的相控阵列线圈线圈n(n＝1-N，N≥2)；用于驱动所述发射线圈、梯度线圈和相控阵列线圈线圈n，以进行整个FOV校准扫描并且得到与相控阵列线圈线圈n相应的复合图像Cn的校准扫描装置；用于根据每个所述复合图像Cn生成绝对值图像|Cn|并且通过将所述绝对值图像|Cn|相加生成求和图像Ab＝∑|Cn|的求和图像生成装置；用于利用以下公式计算在相控阵列线圈线圈n中的像素p的灵敏度因数Sn(p)的灵敏度因数计算装置：
Sn(p)＝Cn(p)/Ab(p)，
式中，将所述复合图像Cn中的像素p的像素值表示为Cn(p)，并且将所述求和图像Ab中的像素p的像素值表示为Ab(p)；以及用于根据所述灵敏度因数Sn(p)生成相控阵列线圈线圈n的相应灵敏度图Sn的灵敏度图生成装置。
第十一方面的MRI装置能够适当地实施如第一方面所描述的线圈灵敏度图生成方法。
按照其第十二方面，本发明提供了具有前述结构的MRI装置，其特征在于：所述灵敏度图生成装置通过对所述灵敏度因数Sn(p)进行曲线拟合，生成灵敏度图Sn。
第十二方面的MRI装置能够适当地实施如第二方面所描述的线圈灵敏度图生成方法。
按照其第十三方面，本发明提供了具有前述结构的MRI装置，其特征在于：所述灵敏度图生成装置利用二维多项式进行二维曲线拟合
第十三方面的MRI装置能够适当地实施如第三方面所描述的线圈灵敏度图生成方法。
按照其第十四方面，本发明提供了具有前述结构的MRI装置，其特征在于：所述灵敏度图生成装置利用二维二次多项式a·x2+b·y2+c·x+d·y+e·x·y+f进行二维曲线拟合。
第十四方面的MRI装置能够适当地实施如第四方面所描述的线圈灵敏度图生成方法。
按照其第十五方面，本发明提供了具有前述结构的MRI装置，其特征在于：所述灵敏度图生成装置利用由所述像素值Ab(p)的平方加权的最小二乘法进行曲线拟合。
第十五方面的MRI装置能够适当地实施如第五方面所描述的线圈灵敏度图生成方法。
按照其第十六方面，本发明提供了一种MRI装置，其特征在于包括：一个用于发射RF脉冲的发射线圈；一个用于施加梯度磁场的梯度线圈；用于接收MR信号的相控阵列线圈线圈n(n＝1-N，N≥2)；用于驱动所述发射线圈、梯度线圈和相控阵列线圈线圈n，以对整个FOV进行校准扫描并且得到与相控阵列线圈线圈n相应的复合图像Cn的校准扫描装置；用于根据每个所述复合图像Cn生成绝对值图像|Cn|并且通过将所述绝对值图像|Cn|相加生成求和图像Ab＝∑|Cn|的绝对值求和图像生成装置；用于利用以下公式计算在相控阵列线圈线圈n中的像素p的灵敏度因数Sn(p)的灵敏度因数计算装置：
Sn(p)＝Cn(p)/Ab(p)，
式中，将所述复合图像Cn中的像素p的像素值表示为Cn(p)，并且将所述求和图像Ab中的像素p的像素值表示为Ab(p)；用于根据所述灵敏度因数Sn(p)生成相控阵列线圈线圈n的相应灵敏度图Sn的灵敏度图生成装置；用于驱动所述发射线圈、梯度线圈和相控阵列线圈线圈n，以进行具有较小FOV的扫描并且得到相应的复合图像An的平行扫描装置；以及用于根据所述灵敏度图Sn和所述复合图像An生成整个FOV的复合图像的整个FOV图像生成装置。
第十六方面的MRI装置能够适当地实施如第六方面所描述的平行成像方法。
按照其第十七方面，本发明提供了具有前述结构的MRI装置，其特征在于：所述灵敏度图生成装置通过对灵敏度因数Sn(p)进行曲线拟合生成灵敏度图Sn。
第十七方面的MRI装置能够适当地实施如第七方面所描述的平行成像方法。
按照其第十八方面，本发明提供了具有前述结构的MRI装置，其特征在于：所述灵敏度图生成装置利用二维多项式进行二维曲线拟合。
第十八方面的MRI装置能够适当地实施如第八方面所描述的平行成像方法。
按照其第十九方面，本发明提供了具有前述结构的MRI装置，其特征在于：所述灵敏度图生成装置利用二维二次多项式a·x2+b·y2+c·x+d·y+e·x·y+f进行二维曲线拟合。
第十九方面的MRI装置能够适当地实施如第九方面所描述的平行成像方法。
按照其第二十方面，本发明提供了具有前述结构的MRI装置，其特征在于：所述灵敏度图生成装置利用由所述像素值Ab(p)的平方加权的最小二乘法进行曲线拟合。
第二十方面的MRI装置能够适当地实施如第十方面所描述的平行成像方法。
按照本发明的线圈灵敏度图生成方法、平行成像方法和MRI装置能够提高利用一般称为SENSE的平行成像方法进行自校准而得到的图像的均匀性。
根据以下对如附图所示的本发明的优选实施例的描述，本发明的其它目的和优点将会更加清楚。
附图简述
图1为按照本发明一个实施例的MRI装置的框图。
图2为表示按照本发明一个实施例的平行成像处理操作的流程图。
图3为表示按照本发明的绝对值求和图像Ab与常规的平方根求和图像As之间的差异的概念性示意图。
图4为示出了按照本发明的灵敏度图S1和S2与常规的灵敏度图S1′和S2′之间的差异的概念性示意图。
图5为示出了常规的平行成像处理的示意性操作的流程图。
本发明的详细描述
以下将参照在附图中示出的实施例，对本发明进行更详细的描述。
在MRI装置100中，磁铁部件1具有一个空腔B，用于将躺在台T上的受检者H放入其中，并且磁铁部件1包括围绕空腔B的一个用于生成恒定静磁场的静磁场线圈M、一个用于生成沿着X、Y和Z轴线的梯度磁场的梯度磁场线圈G以及一个用于发射RF脉冲的体线圈BC。
在受检者H的躯体周围布置用于接收MR信号的相控阵列线圈线圈1和线圈2。
将静磁场线圈M连接到静磁场电源电路2；将梯度磁场线圈G连接到梯度磁场驱动电路3；将体线圈BC连接到RF功率放大器4；并且将相控阵列线圈线圈1和线圈2连接到前置放大器5。
可以用永久磁铁代替静磁场线圈M。
顺序存储器电路8响应来自计算机7的指令，根据储存的脉冲序列，对梯度磁场驱动电路3进行操作，由此从梯度磁场线圈G生成梯度磁场；并且它还对栅极调制电路9进行操作，将来自RF振荡电路10的载波输出信号调制为预定周期和包络线形状的脉冲信号，在RF功率放大器4中对这个脉冲信号进行功率放大，然后提供给人体线圈BC，以发射用于有选择地激励在受检者H内的图像区域的RF脉冲。
前置放大器5对在相控阵列线圈线圈1和线圈2接收到的来自受检者H的MR信号进行放大，并且将该信号输入到相位检测器12。相位检测器12利用来自RF振荡电路10的作为参考信号的载波输出信号对来自前置放大器5的MR信号进行相位检测并且将相位检测的信号提供给A-D转换器11。A-D转换器11将经过相位检测的模拟信号转换为数字数据并且将其输入到计算机7。
计算机7对输入的数字数据进行累加，并且生成相控阵列线圈线圈1和线圈2的线圈灵敏度图。它还进行图像构成处理，从而生成图像。此外，计算机7负责总体控制，如接收从操作员控制台13提供的信息。
显示器设备显示图像和消息。
图2为示出了在MRI装置100进行的平行成像处理操作的流程图。
在步骤V1，利用相控阵列线圈线圈n(n＝1，2)对整个FOV进行校准扫描，从而得到与相控阵列线圈线圈n相应的复合图像Cn。
在步骤V2，生成每个复合图像Cn的绝对值图像|Cn|并且生成绝对值图像|Cn|的求和图像Ab＝∑|Cn|。
在步骤V3，利用以下公式计算在相控阵列线圈线圈n中的像素p的灵敏度因数Sn(p)：
Sn(p)＝Cn(p)/Ab(p)，
式中，将所述复合图像Cn中的像素p的像素值表示为Cn(p)，并且将绝对值求和图像Ab中的像素p的像素值表示为Ab(p)。
在步骤V4，如果操作员命令进行曲线拟合，则将程序转到步骤V5；而如果操作员命令不进行曲线拟合，则将程序转到步骤V6。
在步骤V5，通过将灵敏度因数Sn(p)与二维二次多项式Z＝a·x2+b·y2+c·x+d·y·x·y+f进行二维曲线拟合，生成灵敏度图Sn。这是通过对绝对值求和图像Ab的像素值Ab(p)的平方使用加权最小二乘法完成的。确切地说，计算使评价函数En最小的Zn：
En＝∑Ab(p)2·(Sn(p)-Zn(p)2
并且将Zn定义为灵敏度图Sn。然后，程序转到步骤V7。
在步骤V6，根据灵敏度因数Sn(p)生成灵敏度图Sn。然后程序转到步骤V7。
在步骤V7，利用相控阵列线圈线圈1-线圈N进行具有较小FOV的扫描，从而得到相控阵列线圈线圈1-线圈N的相应复合图像A1-AN。
在步骤V8，根据灵敏度图Sn和复合图像A1-AN生成整个FOV的复合图像。
在步骤V9，重复进行步骤V7和V8，以得到理想数量的图像。
图3为示出了按照本发明的绝对值求和图像Ab与常规的平方根求和图像As之间的差异的概念性示意图。
按照本发明的绝对值求和图像Ab具有比常规的平方根求和图像As更高的均匀性。
图4为示出了按照本发明的灵敏度图S1和S2与常规的灵敏度图S1′和S2′之间的差异的概念性示意图。
按照本发明的灵敏度图S1和S2比常规的灵敏度图S1′和S2′变化得更平缓。
根据MRI装置100，可以得到以下效果：
(1)由于利用绝对值求和图像Ab＝∑|Cn|计算灵敏度因数Sn(p)，因此，与利用求和平方根图像As＝sqrt{∑Cn 2}相比，提高了图像的均匀性。
(2)由于利用二维二次多项式进行二维曲线拟合，因此，可以生成具有平滑表面的灵敏度图，并且可以抑制噪声的影响，从而提高了S/N比。此外，同时利用加权最小二乘法对绝对值求和图像Ab的像素值Ab(p)的平方进行曲线拟合，可以进一步抑制噪声的影响，从而进一步提高了S/N比。
在不脱离本发明的精神和范围的情况下，可以设计本发明的许多不同的实施例。应该理解，除了在所附权利要求中所规定的以外，本发明并不局限于本说明中所描述的特定实施例。

Claims (15)

1.一种线圈灵敏度图生成方法，该方法包括：当将通过对每个相控阵列线圈线圈n(n＝1-N，N≥2)进行整个视场校准扫描而得到的复合图像Cn中的像素p的像素值表示为Cn(p)，将根据所述复合图像Cn生成的绝对值图像|Cn|相加而得到的求和图像Ab＝∑|Cn|中的像素p的像素值表示为Ab(p)时，利用下式计算相控阵列线圈线圈n中的像素p的灵敏度因数Sn(p)
Sn(p)＝Cn(p)/Ab(p)；
并且根据所述灵敏度因数Sn(p)，生成相控阵列线圈线圈n的相应的灵敏度图Sn。
2.如权利要求1所述的线圈灵敏度图生成方法，该方法包括：通过对所述灵敏度因数Sn(p)进行曲线拟合生成灵敏度图Sn。
3.如权利要求2所述的线圈灵敏度图生成方法，该方法包括：利用二维多项式进行二维曲线拟合。
4.如权利要求3所述的线圈灵敏度图生成方法，该方法包括：利用二维二次多项式a·x2+b·y2+c·x+d·y+e·x·y+f进行二维曲线拟合。
5.如权利要求2所述的线圈灵敏度图生成方法，该方法包括：利用由所述绝对值求和图像Ab的像素值Ab(p)的平方加权的最小二乘法进行曲线拟合。
6.一种平行成像方法，该方法包括：通过对相控阵列线圈线圈n(n＝1-N，N≥2)进行整个视场校准扫描得到相应的复合图像Cn；根据每个所述复合图像Cn生成绝对值图像|Cn|；通过将所述绝对值图像|Cn|相加生成求和图像Ab＝∑|Cn|；利用下式计算在相控阵列线圈线圈n中的像素p的灵敏度因数Sn(p)：
Sn(p)＝Cn(p)/Ab(p)，式中，将所述复合图像Cn中的像素p的像素值表示为Cn(p)，并且将所述求和图像Ab中的像素p的像素值表示为Ab(p)；根据所述灵敏度因数Sn(p)生成相控阵列线圈线圈n的相应的灵敏度图Sn；通过对相控阵列线圈线圈n进行具有较小视场的扫描得到相应的复数图An；并且根据所述灵敏度图Sn以及所述复数图An，生成整个视场的复合图像。
7.如权利要求6所述的平行成像方法，该方法包括：通过对所述灵敏度因数Sn(p)进行曲线拟合生成灵敏度图Sn。
8.如权利要求7所述的平行成像方法，该方法包括：利用二维多项式进行二维曲线拟合。
9.如权利要求8所述的平行成像方法，该方法包括：利用二维二次多项式a·x2+b·y2·x+d·y+e·x·y+f进行二维曲线拟合。
10.如权利要求7所述的平行成像方法，该方法包括：利用由所述像素值Ab(p)的平方加权的最小二乘法进行曲线拟合。
11.一种磁共振成像装置，该装置包括：一个用于发射RF脉冲的发射线圈；一个用于施加梯度磁场的梯度线圈；用于接收磁共振信号的相控阵列线圈线圈n(n＝1-N，N≥2)；用于驱动所述发射线圈、梯度线圈和相控阵列线圈线圈n，以对整个视场进行校准扫描并且得到与相控阵列线圈线圈n相应的复合图像Cn的校准扫描装置；用于根据每个所述复合图像Cn生成绝对值图像|Cn|并且通过将所述绝对值图像|Cn|相加生成求和图像Ab＝∑|Cn|的绝对值求和图像生成装置；用于利用以下公式计算在相控阵列线圈线圈n中的像素p的灵敏度因数Sn(p)的灵敏度因数计算装置：
Sn(p)＝Cn(p)/Ab(p)，
式中，将所述复合图像Cn中的像素p的像素值表示为Cn(p)，并且将所述求和图像Ab中的像素p的像素值表示为Ab(p)；用于根据所述灵敏度因数Sn(p)生成相控阵列线圈线圈n的相应灵敏度图Sn的灵敏度图生成装置；用于驱动所述发射线圈、梯度线圈和相控阵列线圈线圈n，以进行具有较小视场的扫描并且得到相应的复合图像An的平行扫描装置；以及用于根据所述灵敏度图Sn和所述复合图像An生成整个视场的复合图像的整个视场图像生成装置。
12.如权利要求11所述的磁共振成像装置，其中，所述灵敏度图生成装置通过对灵敏度因数Sn(p)进行曲线拟合生成灵敏度图Sn。
13.如权利要求12所述的磁共振成像装置，其中，所述灵敏度图生成装置利用二维多项式进行二维曲线拟合。
14.如权利要求13所述的磁共振成像装置，其中，所述灵敏度图生成装置利用二维二次多项式a·x2+b·y2+c·x+d·y+e·x·y+f进行二维曲线拟合。
15.如权利要求12所述的磁共振成像装置，其中，所述灵敏度图生成装置利用由所述像素值Ab(p)的平方加权的最小二乘法进行曲线拟合。
CNB031386695A 2002-06-07 2003-06-06 线圈灵敏图生成方法和平行成像方法以及磁共振成像装置 Expired - Fee Related CN1265762C (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP166994/02 2002-06-07
JP166994/2002 2002-06-07
JP2002166994A JP4030805B2 (ja) 2002-06-07 2002-06-07 コイル感度マップ作成方法およびｍｒｉ装置

Publications (2)

Publication Number Publication Date
CN1468583A CN1468583A (zh) 2004-01-21
CN1265762C true CN1265762C (zh) 2006-07-26

Family

ID=29706743

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB031386695A Expired - Fee Related CN1265762C (zh) 2002-06-07 2003-06-06 线圈灵敏图生成方法和平行成像方法以及磁共振成像装置

Country Status (5)

Country Link
US (1) US7187791B2 (zh)
JP (1) JP4030805B2 (zh)
KR (1) KR100562003B1 (zh)
CN (1) CN1265762C (zh)
DE (1) DE10325886B4 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8878931B2 (en) 2009-03-04 2014-11-04 Honeywell International Inc. Systems and methods for managing video data
US8941464B2 (en) 2005-10-21 2015-01-27 Honeywell International Inc. Authorization system and a method of authorization
US9344684B2 (en) 2011-08-05 2016-05-17 Honeywell International Inc. Systems and methods configured to enable content sharing between client terminals of a digital video management system

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070182410A1 (en) * 2004-01-15 2007-08-09 Niemi Anne K Coil sensitivity estimation for parallel imaging
CN102385044B (zh) * 2004-02-26 2014-12-10 株式会社东芝 磁共振成象装置及磁共振成象装置的数据处理方法
JP4607928B2 (ja) * 2004-04-23 2011-01-05 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 画像生成方法およびｍｒｉ装置
JP4037385B2 (ja) * 2004-04-23 2008-01-23 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 画像生成方法およびｍｒｉ装置
JP4647959B2 (ja) * 2004-09-10 2011-03-09 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｍｒｉ装置、画像均一度評価方法および画像均一度評価装置
JP2006075380A (ja) * 2004-09-10 2006-03-23 Ge Medical Systems Global Technology Co Llc Ｍｒ画像生成方法およびｍｒｉ装置
JP5127116B2 (ja) * 2005-02-15 2013-01-23 株式会社東芝 磁気共鳴イメージング装置
JP2007090001A (ja) * 2005-09-30 2007-04-12 Ge Medical Systems Global Technology Co Llc Ｍｒスキャン方法およびｍｒｉ装置
US8320647B2 (en) * 2007-11-20 2012-11-27 Olea Medical Method and system for processing multiple series of biological images obtained from a patient
CN101498773B (zh) * 2008-01-29 2012-05-23 西门子（中国）有限公司 一种提高图像亮度均匀度的方法和装置
JP5443695B2 (ja) * 2008-03-05 2014-03-19 株式会社東芝 磁気共鳴イメージング装置
CN101545960B (zh) * 2008-03-26 2012-07-18 西门子（中国）有限公司 筒状表面线圈及用于其的射频信号处理方法和系统
EP2233941A1 (en) * 2009-03-25 2010-09-29 Koninklijke Philips Electronics N.V. MR imaging using parallel signal acquisition
JP5290868B2 (ja) * 2009-05-28 2013-09-18 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー コイル感度推定装置、磁気共鳴イメージング装置、コイル感度推定方法、およびプログラム
CN101872001B (zh) * 2010-06-29 2012-11-07 中国科学院生物物理研究所 并行发射接收射频接口电路和相控阵发射接收头线圈
CN102008306B (zh) * 2010-12-22 2012-11-14 中国科学院深圳先进技术研究院 磁共振并行成像方法
US20130289912A1 (en) * 2012-03-30 2013-10-31 Siemens Aktiengesellschaft Eigen-vector approach for coil sensitivity maps estimation
JP5897415B2 (ja) * 2012-06-29 2016-03-30 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
CN103091656B (zh) * 2012-12-25 2015-06-03 深圳先进技术研究院 基于正则化约束多项式拟合磁共振线圈灵敏度的计算方法
JP6430417B2 (ja) * 2013-03-13 2018-11-28 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． パラレルイメージング加速パラメータの自動最適化
EP3692381B1 (en) 2017-10-06 2023-12-06 University of Cincinnati Systems and methods for estimating complex radiofrequency fields in a magnetic resonance imaging
CN108287325B (zh) * 2018-01-03 2020-08-11 上海东软医疗科技有限公司 一种图像重建方法、装置及设备
US10866292B2 (en) * 2018-05-31 2020-12-15 General Electric Company Methods and systems for coil selection in magnetic resonance imaging
KR102584250B1 (ko) * 2022-06-22 2023-10-05 주식회사 에어스메디컬 데이터 처리 방법, 컴퓨터 프로그램 및 장치

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5659629A (en) 1995-04-07 1997-08-19 General Electric Company Visualization of a multi-dimensional data set acquired with a surface receive coil
JP3600351B2 (ja) * 1996-03-12 2004-12-15 株式会社東芝 磁気共鳴イメージング装置
DE19653476A1 (de) 1996-12-20 1998-07-02 Siemens Ag Verfahren zur Intensitätskorrektur eines Bildes
US6289232B1 (en) 1998-03-30 2001-09-11 Beth Israel Deaconess Medical Center, Inc. Coil array autocalibration MR imaging
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
JP2004504910A (ja) * 2000-07-31 2004-02-19 コーニンクレッカ　フィリップス　エレクトロニクス　エヌ　ヴィ 高速動的撮像を形成する磁気共鳴方法
US6483308B1 (en) 2000-08-31 2002-11-19 Ge Medical Systems Global Technology Company, Llc Method and apparatus for processing MRI data acquired with a plurality of coils using dixon techniques
CN1491095A (zh) * 2001-01-19 2004-04-21 株式会社东芝 使用高精度线圈感度图的并行mr成象
US6559642B2 (en) 2001-05-09 2003-05-06 Ge Medical Systems Global Technology Company, Llc Calibration method for use with sensitivity encoding MRI acquisition

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8941464B2 (en) 2005-10-21 2015-01-27 Honeywell International Inc. Authorization system and a method of authorization
US8878931B2 (en) 2009-03-04 2014-11-04 Honeywell International Inc. Systems and methods for managing video data
US9344684B2 (en) 2011-08-05 2016-05-17 Honeywell International Inc. Systems and methods configured to enable content sharing between client terminals of a digital video management system

Also Published As

Publication number Publication date
KR20030095298A (ko) 2003-12-18
CN1468583A (zh) 2004-01-21
KR100562003B1 (ko) 2006-03-22
DE10325886B4 (de) 2009-12-03
JP2004008533A (ja) 2004-01-15
US7187791B2 (en) 2007-03-06
JP4030805B2 (ja) 2008-01-09
US20030228043A1 (en) 2003-12-11
DE10325886A1 (de) 2004-02-12

Similar Documents

Publication Publication Date Title
CN1265762C (zh) 线圈灵敏图生成方法和平行成像方法以及磁共振成像装置
CN100339047C (zh) 磁共振成像设备和磁共振成像方法
CN1251645C (zh) 磁共振成像装置
CN1268285C (zh) 磁共振成像装置
CN1299121C (zh) 磁共振成像设备
CN1217622C (zh) 减少快速自旋回波磁共振图象中麦克斯韦项假象的方法
CN1231177C (zh) 磁共振成像方法及实施该方法的装置
CN1745705A (zh) 磁共振影像生成方法和磁共振成像装置
CN1203808C (zh) 核磁共振成像装置和方法
CN1491095A (zh) 使用高精度线圈感度图的并行mr成象
CN1891150A (zh) 具有减少的fse尖点伪影的mr成像方法及系统
CN1374069A (zh) 磁共振成像装置及其磁场变化测量方法和补偿方法
CN1851494A (zh) 改进发送侧加速的体积选择磁共振成像的方法和设备
CN101051074A (zh) 磁共振成像装置和磁共振成像方法
CN1530072A (zh) 核磁共振成像设备和核磁共振成像方法
JP2012522563A (ja) ｋ空間に依存するＲＦパルス選択によるパラレル送信におけるＳＡＲ低減
CN1828329A (zh) 磁共振检测器及方法
CN1251644C (zh) 多片层磁共振成象方法和装置
CN1576875A (zh) 磁场非单调变化时在自旋回波图像中避免外围干扰的方法
CN1264023C (zh) 数据采集方法、补偿磁场偏移的方法以及磁共振成像装置
CN1192743C (zh) Mr成像方法,相位误差测量方法和mri系统
CN1975395A (zh) 改进核磁共振频谱分析中弱敏感原子核类探测的方法和装置
JP6552804B2 (ja) 磁気共鳴イメージング装置
US10775469B2 (en) Magnetic resonance imaging apparatus and method
CN109975343A (zh) 磁共振成像方法、系统及存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20060726

Termination date: 20140606