JP2018093989A - Magnetic resonance imaging device and reception coil sensitivity map calculation method - Google Patents

Magnetic resonance imaging device and reception coil sensitivity map calculation method Download PDF

Info

Publication number
JP2018093989A
JP2018093989A JP2016239822A JP2016239822A JP2018093989A JP 2018093989 A JP2018093989 A JP 2018093989A JP 2016239822 A JP2016239822 A JP 2016239822A JP 2016239822 A JP2016239822 A JP 2016239822A JP 2018093989 A JP2018093989 A JP 2018093989A
Authority
JP
Japan
Prior art keywords
sensitivity
reception
map
coil
sensitivity map
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
JP2016239822A
Other languages
Japanese (ja)
Other versions
JP6783642B2 (en
Inventor
眞次 黒川
Shinji Kurokawa
眞次 黒川
康弘 鎌田
Yasuhiro Kamata
康弘 鎌田
直哉 坂口
Naoya Sakaguchi
直哉 坂口
致知 濱田
Yoshitomo Hamada
致知 濱田
秀太 羽原
Hideta Habara
秀太 羽原
毅倫 村瀬
Takemichi Murase
毅倫 村瀬
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.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2016239822A priority Critical patent/JP6783642B2/en
Publication of JP2018093989A publication Critical patent/JP2018093989A/en
Application granted granted Critical
Publication of JP6783642B2 publication Critical patent/JP6783642B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

PROBLEM TO BE SOLVED: To accurately calculate a reception coil sensitivity distribution (sensitivity map) in the presence of a subject, and thereby improve quality of a subject image created by using sensitivity information.SOLUTION: Regarding a reception coil having a plurality of reception channels, relationships among sensitivity maps of a plurality of reception channels acquired outside of a test object are obtained first as an index (first index). Relationships among images (reference images) of a plurality of reception channels acquired by setting the test object are obtained as a second index. A coordinate conversion formula is calculated, which causes the map of the second index to coincide with the map of the first index. The sensitivity map of the reception coil acquired outside of the test object is corrected by a coordinate conversion formula, and a reception sensitivity map which is substantially equivalent to the reception sensitivity map when the subject is set is created.SELECTED DRAWING: Figure 4

Description

本発明は、磁気共鳴イメージング装置（以下、ＭＲＩ装置という）に関し、特にＭＲＩ装置の受信コイルの正確な感度分布を算出する技術に関する。   The present invention relates to a magnetic resonance imaging apparatus (hereinafter referred to as an MRI apparatus), and more particularly to a technique for calculating an accurate sensitivity distribution of a receiving coil of the MRI apparatus.

ＭＲＩ装置では、高周波磁場である核磁気共鳴信号を受信するために受信用ＲＦコイル（以下、受信コイル）が用いられる。受信コイルには、広い感度分布を持ち全身撮像用として用いられるボディーコイルと、狭い領域で高い感度分布を持ち局所撮像に好適な表面コイル（局所コイルともいう）とがある。一般的なＭＲＩ装置には、標準でガントリーにボディーコイルが内蔵されており、撮像条件を決めるためのプリスキャンなどでは、ボディーコイルを用いた撮像で得られた情報も利用される。   In the MRI apparatus, a receiving RF coil (hereinafter referred to as a receiving coil) is used to receive a nuclear magnetic resonance signal that is a high-frequency magnetic field. The receiving coil includes a body coil that has a wide sensitivity distribution and is used for whole body imaging, and a surface coil (also referred to as a local coil) that has a high sensitivity distribution in a narrow region and is suitable for local imaging. In general MRI apparatuses, a body coil is built in a gantry as standard, and information obtained by imaging using the body coil is also used in pre-scanning for determining imaging conditions.

ＭＲＩ装置で得られる画像は、受信コイルの感度分布の影響を受ける。特に表面コイルは感度が高いものの感度分布の均一性に劣る。このため、撮像後の画像は、受信コイルの感度分布の影響を排除するための受信コイルの感度情報を用いた補正（感度補正）が行われる。   An image obtained by the MRI apparatus is affected by the sensitivity distribution of the receiving coil. In particular, the surface coil has high sensitivity but is inferior in uniformity of sensitivity distribution. For this reason, the captured image is subjected to correction (sensitivity correction) using sensitivity information of the receiving coil to eliminate the influence of the sensitivity distribution of the receiving coil.

感度補正の精度を高めるためには正確な感度情報が必要である。受信コイルの感度情報は、例えば、均質なファントムの画像から得られる（即ち、当該画像自体が感度情報である）が、受信コイルの感度は、被検体（例えば人体）が存在するときと、ファントムを用いたときとでは異なる。被検体が存在する条件で受信コイルの感度を求めるため、被検体を撮像して得たデータを用い、低周波領域のデータのみを用いた再構成やフィルター処理等の画像処理によって受信コイルの感度分布を求める技術（自己参照型補正）がある（例えば特許文献１〜３）。特許文献３には、自己参照型補正の改良技術として、Ｂ１マップを用いて補正したプロトン強調画像をフィルタリングして受信感度分布を得る技術が開示されている。   In order to increase the accuracy of sensitivity correction, accurate sensitivity information is required. The sensitivity information of the receiving coil is obtained from, for example, a homogeneous phantom image (that is, the image itself is sensitivity information), but the sensitivity of the receiving coil is determined when the subject (for example, a human body) exists and when the phantom is present. It is different from when using. In order to determine the sensitivity of the receiving coil in the presence of the subject, the sensitivity of the receiving coil is obtained by image processing such as reconstruction or filtering using only data in the low frequency region, using data obtained by imaging the subject. There is a technique (self-reference type correction) for obtaining a distribution (for example, Patent Documents 1 to 3). Patent Document 3 discloses a technique for obtaining a reception sensitivity distribution by filtering a proton-weighted image corrected using a B1 map as an improved technique of self-reference type correction.

また表面コイルは感度が不均一なのでリファレンススキャンを行い、比較的均一なボディコルとの比で感度を補正することも提案されている（特許文献４）。
表面コイルの感度分布とボディーコイルの感度分布の比は被検体のリファレンススキャンで求めることができるので、表面コイルで撮像したときに比較的均一なボディーコイルの感度むら程度に補正することはできる。しかし、ボディーコイルの感度分布が均一でないときにはボディーコイルの感度分布補正も必要になる。その際、自己参照のような手法でボディーコイルの感度分布を求めると、参照画像のコントラストによっては失敗することがある。 In addition, since the surface coil has non-uniform sensitivity, it is also proposed to perform a reference scan and correct the sensitivity in a ratio with a relatively uniform body col (Patent Document 4).
Since the ratio between the sensitivity distribution of the surface coil and the sensitivity distribution of the body coil can be obtained by the reference scan of the subject, it can be corrected to a relatively uniform body coil sensitivity unevenness when imaged by the surface coil. However, when the sensitivity distribution of the body coil is not uniform, it is also necessary to correct the sensitivity distribution of the body coil. At this time, if the sensitivity distribution of the body coil is obtained by a method such as self-reference, it may fail depending on the contrast of the reference image.

さらに、ボディーコイルを用いた補正後に残るボディーコイルの感度不均一を補正するためにさまざまな技術が提案されている（特許文献５〜８）。例えば、特許文献５には、複数のコイルから得た画像を用いて画像の初期値を推定し、これにｈｏｍｏｍｏｒｐｈｉｃ ｆｉｌｔｅｒと多チャンネルブラインドデコンボリューションを作用させて感度分布を得る手法が提案されている。   Furthermore, various techniques have been proposed to correct non-uniform sensitivity of the body coil remaining after correction using the body coil (Patent Documents 5 to 8). For example, Patent Document 5 proposes a method of estimating an initial value of an image using images obtained from a plurality of coils, and obtaining a sensitivity distribution by applying a homomorphic filter and multichannel blind deconvolution to the image. .

また特許文献６には、表面コイルと同じ領域において感度分布がほぼ均一と見做せるボディーコイルで得た画像を用いて、表面コイルの感度分布情報を補正する技術が採用されている。   Patent Document 6 employs a technique for correcting the sensitivity distribution information of the surface coil using an image obtained with a body coil that can be considered to have a substantially uniform sensitivity distribution in the same region as the surface coil.

さらに特許文献６に記載された技術では、ボディーコイルの不均一性を画像処理で補正している。また、参照とするボディーコイルの感度分布を補正した後、補正後のボディーコイルの感度分布を用いて表面コイルの感度分布を求める技術も提案されている。例えば、特許文献７には、送信感度と受信感度は等しいという前提のもと、送信コイルのＢ１分布をシミングで均一にして、受信コイルの受信感度を補正する技術が開示されている。
特許文献８では、複数ＴＥの撮像でＴ２を求め、Ｔ２から組織とそのプロトン密度を推定してＴ２とプロトン密度から撮像時のＴＥでの参照画像を作成し、感度マップを求める技術が開示されているが、この技術では、Ｔ２の計測が必要になる上に、様々な誤差（Ｔ２とプロトン密度の対応、計算されるＴ２値）に影響される恐れがある。 Furthermore, in the technique described in Patent Document 6, the non-uniformity of the body coil is corrected by image processing. In addition, a technique has been proposed in which the sensitivity distribution of the body coil to be referenced is corrected, and then the sensitivity distribution of the surface coil is obtained using the corrected sensitivity distribution of the body coil. For example, Patent Document 7 discloses a technique for correcting the reception sensitivity of the reception coil by making the B1 distribution of the transmission coil uniform by shimming on the assumption that the transmission sensitivity and the reception sensitivity are equal.
Patent Document 8 discloses a technique for obtaining a sensitivity map by obtaining T2 by imaging a plurality of TEs, estimating a tissue and its proton density from T2 and creating a reference image at TE at the time of imaging from T2 and proton density. However, this technique requires measurement of T2, and may be affected by various errors (correspondence between T2 and proton density, calculated T2 value).

特表２００７−５０３２３９号公報（米国特許出願公開第２００６／０２６１８０９号公報）Japanese translation of PCT publication No. 2007-503239 (US Patent Application Publication No. 2006/0261809) 米国特許第６６８０６１０号公報US Pat. No. 6,680,610 米国特許第９３１６７０７号公報U.S. Pat. No. 9,316,707 ＥＰ０６９５９４Ａ１（１９９６／０２／０７）公報EP069594A1 (1996/02/07) publication 米国特許出願公開第２０１３／０２８１８２２号US Patent Application Publication No. 2013/0281822 特開２００９−２０７７５６号公報（米国特許出願公開第２００９／０２２４７５６号公報）JP 2009-207756 A (US Patent Application Publication No. 2009/0224756) 特許第５７０６８９９号（米国特許第９０３６８８４号）Patent No. 5706899 (U.S. Pat. No. 9,036,884) 特許第４１３０４０５号Japanese Patent No. 4130405

本発明は、ボディーコイル等の感度分布補正について参照画像のコントラストに影響されず、被検体存在下での受信コイルの感度分布（感度マップ）を精度よく算出することができ、それにより感度情報を利用して作成される被検体画像の画質を向上することを課題とする。本発明の感度分布算出方法はボディーコイルだけでなく表面コイルにも適用することができる。   The present invention can accurately calculate the sensitivity distribution (sensitivity map) of the receiving coil in the presence of the subject without being affected by the contrast of the reference image with respect to the sensitivity distribution correction of the body coil or the like. It is an object of the present invention to improve the image quality of a subject image created by use. The sensitivity distribution calculation method of the present invention can be applied not only to a body coil but also to a surface coil.

上記課題を解決するため、本発明は、複数の受信チャンネルを有する受信コイルについて、検査対象以外で取得した複数の受信チャンネルの感度マップ間の関係を指標（第一の指標）として求めておき、検査対象に対して取得した複数の受信チャンネルの画像（参照画像）間の関係を第二の指標として求めて、第二の指標のマップと第一の指標のマップを似たような分布に変形する座標変換式を算出する。検査対象以外で取得した受信コイルの感度マップを、前記座標変換式で補正し、検査対象に対する受信感度マップと実質的に等価な受信感度マップを作成する。撮像に用いる受信コイル、即ち受信感度マップを作成する受信コイルは、座標変換式を算出する際に用いた複数の受信チャンネルを有する受信コイルの一部或いは全部でもよいし、それとは異なるコイルでもよい。   In order to solve the above-described problem, the present invention obtains, as an index (first index), a relationship between sensitivity maps of a plurality of reception channels acquired other than the inspection target for a reception coil having a plurality of reception channels, The relationship between the images (reference images) of multiple reception channels acquired for the inspection target is obtained as a second index, and the second index map and the first index map are transformed into similar distributions. A coordinate conversion formula is calculated. The sensitivity map of the receiving coil acquired outside the inspection target is corrected by the coordinate conversion formula, and a reception sensitivity map substantially equivalent to the reception sensitivity map for the inspection target is created. The receiving coil used for imaging, that is, the receiving coil for creating the receiving sensitivity map, may be a part or all of the receiving coil having a plurality of receiving channels used when calculating the coordinate conversion formula, or may be a different coil. .

また広域用受信コイルの受信感度マップを作成した場合には、さらに被検体の撮像に用いる異なる受信コイルの感度マップを、作成した広域用受信コイルの受信感度マップで補正してもよい。   When the reception sensitivity map of the wide-area reception coil is created, the sensitivity map of different reception coils used for imaging the subject may be corrected with the created reception sensitivity map of the wide-area reception coil.

すなわち、本発明のＭＲＩ装置は、複数の受信チャンネルを有する受信コイルを備え、検査対象の核磁気共鳴信号を受信し、前記検査対象の画像を取得する撮像部と、前記受信コイルの感度情報を用いて前記検査対象の画像作成に必要な計算を行う計算部と、を備え、前記検査対象以外で予め取得した前記受信コイルの感度マップ及び前記複数の受信チャンネルそれぞれの感度マップまたは前記複数の受信チャンネルそれぞれの感度マップから導出したマップを記憶する記憶部を更に備え、前記計算部は、前記複数の受信チャンネルの感度マップ間の関係を表す第一の指標マップと、前記撮像部が前記検査対象について前記受信コイルの受信チャンネルごとに取得した複数の参照画像間の関係を表す第二の指標マップとの差を最小化するように前記第一の指標マップを座標変換する座標変換式を算出する座標変換式算出部と、前記座標変換式を用いて、前記検査対象存在時の受信感度マップを作成する感度マップ作成部と、を備えることを特徴とする。
なお前記複数の受信チャンネルそれぞれの感度マップから導出したマップは、前記第一の指標マップのほか、感度マップから第一の指標マップを算出する途中のデータも含んでもよい。 That is, the MRI apparatus of the present invention includes a receiving coil having a plurality of receiving channels, receives an nuclear magnetic resonance signal to be examined, acquires an image of the examination target, and sensitivity information of the receiving coil. And a calculation unit that performs calculations necessary for creating an image of the inspection target, and a sensitivity map of the reception coil and a sensitivity map of each of the plurality of reception channels or the plurality of receptions acquired in advance other than the inspection target A storage unit that stores a map derived from the sensitivity map of each channel; and the calculation unit includes a first index map that represents a relationship between the sensitivity maps of the plurality of reception channels, and the imaging unit includes the inspection target. So as to minimize a difference from a second index map representing a relationship between a plurality of reference images acquired for each reception channel of the reception coil. A coordinate conversion formula calculation unit that calculates a coordinate conversion formula for converting the coordinates of the first index map, and a sensitivity map generation unit that generates a reception sensitivity map when the inspection target exists using the coordinate conversion formula; It is characterized by providing.
The map derived from the sensitivity map of each of the plurality of reception channels may include data in the middle of calculating the first index map from the sensitivity map in addition to the first index map.

本発明によれば、複数の受信チャンネルを有する受信コイルのあらかじめ用意した感度マップ及び参照画像を用いて算出した感度マップの座標変換式を用いることにより、参照画像のコントラストに影響されず、精度よく受信感度マップを作成することができる。   According to the present invention, by using a sensitivity map coordinate conversion formula calculated using a sensitivity map prepared in advance and a reference image of a reception coil having a plurality of reception channels, the reference image is not affected by the contrast of the reference image and is accurately obtained. A reception sensitivity map can be created.

本発明が適用されるＭＲＩ装置の一例の全体構成を示す図。The figure which shows the whole structure of an example of the MRI apparatus with which this invention is applied. 実施形態のＭＲＩ装置の処理の概要を示すフロー図。The flowchart which shows the outline | summary of the process of the MRI apparatus of embodiment. 第一実施形態のＭＲＩ装置の信号処理部の機能ブロック図。The functional block diagram of the signal processing part of the MRI apparatus of 1st embodiment. 第一実施形態の信号処理部の処理の一例を示すフロー図。The flowchart which shows an example of the process of the signal processing part of 1st embodiment. 第一実施形態における座標変換処理の概要を模式的に示す図。The figure which shows typically the outline | summary of the coordinate transformation process in 1st embodiment. 第一実施形態の信号処理部の処理の他の例を示すフロー図。The flowchart which shows the other example of the process of the signal processing part of 1st embodiment. 第一実施形態により処理の効果を説明する図。The figure explaining the effect of a process by 1st embodiment. 第二実施形態の信号処理部の機能ブロック図。The functional block diagram of the signal processing part of 2nd embodiment. 第二実施形態の信号処理部の処理の流れを示すフロー図。The flowchart which shows the flow of a process of the signal processing part of 2nd embodiment. （ａ）シミュレーションで得た初期感度マップ、（ｂ）ファントム撮像で得た感度マップ、（ｃ）補正後の初期感度マップを示す図。(A) The initial sensitivity map obtained by simulation, (b) The sensitivity map obtained by phantom imaging, (c) The figure which shows the initial sensitivity map after correction | amendment. 第三実施形態の信号処理部の機能ブロック図。The functional block diagram of the signal processing part of 3rd embodiment. 第三実施形態の信号処理部の処理の流れを示すフロー図。The flowchart which shows the flow of a process of the signal processing part of 3rd embodiment. 第四実施形態における受信コイルの組み合わせ例を示す図。The figure which shows the example of a combination of the receiving coil in 4th embodiment. 第四実施形態の信号処理部の処理の流れを示すフロー図。The flowchart which shows the flow of a process of the signal processing part of 4th embodiment.

最初に、本発明が適用されるＭＲＩ装置の全体概要を、図面を参照して説明する。図１は、本発明に係るＭＲＩ装置の一実施形態の全体構成を示すブロック図である。このＭＲＩ装置は、ＮＭＲ現象を利用して被検体の断層画像を得るもので、静磁場発生部２と、傾斜磁場発生部３と、送信部５と、受信部６と、信号処理部７と、シーケンサ４と、中央処理装置（ＣＰＵ）８とを備えている。   First, an overall outline of an MRI apparatus to which the present invention is applied will be described with reference to the drawings. FIG. 1 is a block diagram showing the overall configuration of an embodiment of an MRI apparatus according to the present invention. This MRI apparatus uses a NMR phenomenon to obtain a tomographic image of a subject, and includes a static magnetic field generation unit 2, a gradient magnetic field generation unit 3, a transmission unit 5, a reception unit 6, and a signal processing unit 7. A sequencer 4 and a central processing unit (CPU) 8 are provided.

静磁場発生部２は、被検体１の周りの空間に均一な静磁場を発生させるもので、永久磁石方式、常電導方式あるいは超電導方式の静磁場発生装置を備える。発生する磁場の方向により、垂直磁場方式、水平磁場方式などがあり、本発明はいずれの方式にも適用できる。   The static magnetic field generation unit 2 generates a uniform static magnetic field in the space around the subject 1 and includes a permanent magnet type, normal conduction type or superconductivity type static magnetic field generation apparatus. Depending on the direction of the generated magnetic field, there are a vertical magnetic field method, a horizontal magnetic field method, and the like, and the present invention can be applied to any method.

傾斜磁場発生部３は、ＭＲＩ装置の座標系（静止座標系）であるＸ、Ｙ、Ｚの３軸方向に傾斜磁場を印加する傾斜磁場コイル９と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源１０とから成り、シーケンサ４からの命令に従ってそれぞれのコイルの傾斜磁場電源１０を駆動することにより、Ｘ、Ｙ、Ｚの３軸方向に傾斜磁場Ｇｘ、Ｇｙ、Ｇｚを印加する。これら３軸方向の傾斜磁場の組み合わせにより、任意の方向の傾斜磁場を発生することができる。例えば、撮影時には、スライス面（撮影断面）に直交する方向にスライス方向傾斜磁場パルス（Ｇｓ）を印加して被検体１に対するスライス面を設定し、そのスライス面に直交して且つ互いに直交する残りの２つの方向に位相エンコード方向傾斜磁場パルス（Ｇｐ）と周波数エンコード方向傾斜磁場パルス（Ｇｆ）を印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。   The gradient magnetic field generator 3 includes a gradient magnetic field coil 9 that applies gradient magnetic fields in the three-axis directions of X, Y, and Z, which are coordinate systems (stationary coordinate systems) of the MRI apparatus, and gradient magnetic fields that drive the respective gradient magnetic field coils. A gradient power supply Gx, Gy, Gz is applied in the three axial directions of X, Y, and Z by driving the gradient magnetic field power supply 10 of each coil according to a command from the sequencer 4. By combining these three-axis gradient magnetic fields, a gradient magnetic field in any direction can be generated. For example, at the time of imaging, a slice direction gradient magnetic field pulse (Gs) is applied in a direction orthogonal to the slice plane (imaging cross section) to set a slice plane for the subject 1, and the remaining orthogonal to the slice plane and orthogonal to each other The phase encoding direction gradient magnetic field pulse (Gp) and the frequency encoding direction gradient magnetic field pulse (Gf) are applied in the two directions, and position information in each direction is encoded in the echo signal.

シーケンサ４は、高周波磁場パルス（以下、「ＲＦパルス」という）と傾斜磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段で、ＣＰＵ８の制御で動作し、被検体１の断層画像のデータ収集に必要な種々の命令を送信部５、傾斜磁場発生部３、および受信部６に送る。   The sequencer 4 is a control unit that repeatedly applies a high-frequency magnetic field pulse (hereinafter referred to as “RF pulse”) and a gradient magnetic field pulse in a predetermined pulse sequence. The sequencer 4 operates under the control of the CPU 8 and collects tomographic image data of the subject 1. Various commands necessary for the transmission are sent to the transmission unit 5, the gradient magnetic field generation unit 3, and the reception unit 6.

送信部５は、被検体１の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体１にＲＦパルスを照射するもので、高周波発振器１１と変調器１２と高周波増幅器１３と送信側の高周波コイル（送信コイル）１４ａとを備える。高周波発振器１１から出力されたＲＦパルスをシーケンサ４からの指令によるタイミングで変調器１２により振幅変調し、この振幅変調されたＲＦパルスを高周波増幅器１３で増幅した後に被検体１に近接して配置された送信コイル１４ａに供給することにより、ＲＦパルスが被検体１に照射される。   The transmitter 5 irradiates the subject 1 with an RF pulse in order to cause nuclear magnetic resonance to occur in the nuclear spins of the atoms constituting the living tissue of the subject 1, and includes a high frequency oscillator 11, a modulator 12, and a high frequency amplifier. 13 and a high frequency coil (transmission coil) 14a on the transmission side. The RF pulse output from the high-frequency oscillator 11 is amplitude-modulated by the modulator 12 at a timing according to a command from the sequencer 4, and after the amplitude-modulated RF pulse is amplified by the high-frequency amplifier 13, the RF pulse is arranged close to the subject 1. By supplying to the transmission coil 14a, the subject 1 is irradiated with the RF pulse.

受信部６は、被検体１の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号（ＮＭＲ信号）を検出するもので、受信側の高周波コイル（受信コイル） １４ｂと信号増幅器１５と直交位相検波器１６とＡ／Ｄ変換器１７とを備える。送信コイル１４ａから照射された電磁波によって誘起された被検体１の応答のＮＭＲ信号が被検体１に近接して配置された受信コイル１４ｂで検出され、信号増幅器１５で増幅された後、シーケンサ４からの指令によるタイミングで直交位相検波器１６により直交する二系統の信号に分割され、それぞれがＡ／Ｄ変換器１７でディジタル量に変換されて、信号処理部７に送られる。   The receiving unit 6 detects an echo signal (NMR signal) emitted by nuclear magnetic resonance of nuclear spins constituting the biological tissue of the subject 1, and includes a receiving-side high-frequency coil (receiving coil) 14 b and a signal amplifier 15. And a quadrature detector 16 and an A / D converter 17. The NMR signal of the response of the subject 1 induced by the electromagnetic wave irradiated from the transmission coil 14a is detected by the reception coil 14b arranged close to the subject 1, amplified by the signal amplifier 15, and then from the sequencer 4. Are divided into two orthogonal signals by the quadrature detector 16, converted into digital quantities by the A / D converter 17, and sent to the signal processor 7.

信号処理部７は、各種データ処理と処理結果の表示及び保存等を行うもので、その機能の一部はＣＰＵ８により実現される。信号処理部７（ＣＰＵ８）が行う処理には、被検体１の画像再構成や被検体の特性を示す数値の演算、受信コイル感度等の装置の特性を用いた信号や処理結果の補正などが含まれる。   The signal processing unit 7 performs various data processing, processing result display and storage, and the like, and a part of its functions is realized by the CPU 8. The processing performed by the signal processing unit 7 (CPU 8) includes image reconstruction of the subject 1, calculation of numerical values indicating the characteristics of the subject, correction of signals and processing results using device characteristics such as reception coil sensitivity, and the like. included.

信号処理部７は、光ディスク、磁気ディスク等の外部記憶装置１８と、ＣＲＴ等からなるディスプレイ１９と、トラックボール又はマウス、キーボード等から成る操作部２０とを備えている。受信部６からのデータがＣＰＵ８に入力されると、ＣＰＵ８が信号処理、画像再構成等の処理を実行し、その結果である被検体１の断層画像をディスプレイ１９に表示すると共に、外部記憶装置１８の磁気ディスク等に記録する。ＣＰＵ８は、上述した処理を行う計算部として機能するとともにシーケンサ４や装置全体の制御を行う制御部として機能する。   The signal processing unit 7 includes an external storage device 18 such as an optical disk or a magnetic disk, a display 19 including a CRT, and an operation unit 20 including a trackball, a mouse, a keyboard, or the like. When data from the receiving unit 6 is input to the CPU 8, the CPU 8 executes processing such as signal processing and image reconstruction, displays a tomographic image of the subject 1 as a result on the display 19, and an external storage device. Recorded on 18 magnetic disks and the like. The CPU 8 functions as a calculation unit that performs the above-described processing, and also functions as a control unit that controls the sequencer 4 and the entire apparatus.

操作部２０は、ＭＲＩ装置の各種制御情報や信号処理部７で行う処理の制御情報を入力するもので、ディスプレイ１９に近接して配置され、操作者がディスプレイ１９を見ながら操作部２０を通してインタラクティブにＭＲＩ装置の各種処理を制御する。   The operation unit 20 is used to input various control information of the MRI apparatus and control information of processing performed by the signal processing unit 7. The operation unit 20 is disposed in the vicinity of the display 19, and the operator interactively through the operation unit 20 while looking at the display 19. In addition, various processes of the MRI apparatus are controlled.

なお、図１において、送信側の高周波コイル１４ａと傾斜磁場コイル９は、被検体１が挿入される静磁場発生部２の静磁場空間内に、垂直磁場方式であれば被検体１に対向して、水平磁場方式であれば被検体１を取り囲むようにして設置されている。また、受信側の高周波コイル１４ｂは、被検体１に対向して、或いは取り囲むように設置されている。   In FIG. 1, the high-frequency coil 14a and the gradient magnetic field coil 9 on the transmission side are opposed to the subject 1 in the static magnetic field space of the static magnetic field generating unit 2 in which the subject 1 is inserted if the vertical magnetic field method is used. If the horizontal magnetic field method is used, the subject 1 is installed so as to surround it. The high-frequency coil 14b on the receiving side is installed so as to face or surround the subject 1.

送信コイル１４ａと受信コイル１４ｂとは、別個の高周波コイルでもよいが、一つの高周波コイルが送信コイル１４ａと受信コイル１４ｂを兼ねる場合もある。また送受信を兼ねる全身用コイルと、それとは別の局所コイル等の高周波コイルを組み合わせて用いる場合も有る。本発明が適用されるＭＲＩ装置では、少なくとも一つの高周波コイルとして、複数の受信チャンネルを持つ受信コイルを備える。「複数の受信チャンネルを持つ」ということは、複数の小型受信コイルを組み合わせた受信コイルであって、各小型受信コイルがそれぞれＮＭＲ信号を受信する受信コイル、例えばマルチアレイコイルであること、受信コイルとしては一体であるが受信ポートを複数持ち、各受信ポートがそれぞれＮＭＲ信号受信する受信コイル、例えばバードケージ型コイル等であること、のいずれも含む趣旨である。以下、これらを総称して「多チャンネル受信コイル」或いは「多チャンネルコイル」という。   The transmission coil 14a and the reception coil 14b may be separate high-frequency coils, but one high-frequency coil may also serve as the transmission coil 14a and the reception coil 14b. In some cases, a whole-body coil that also performs transmission and reception and a high-frequency coil such as another local coil are used in combination. The MRI apparatus to which the present invention is applied includes a reception coil having a plurality of reception channels as at least one high-frequency coil. “Having a plurality of receiving channels” means a receiving coil in which a plurality of small receiving coils are combined, each small receiving coil receiving a NMR signal, for example, a multi-array coil, a receiving coil Although it is integrated, it has the meaning of including any of receiving coils having a plurality of receiving ports, each receiving port receiving a NMR signal, such as a birdcage type coil. Hereinafter, these are collectively referred to as “multi-channel receiving coil” or “multi-channel coil”.

次に、上述した構成のＭＲＩ装置の撮像及び信号処理の概要を、図２のフローを参照して説明する。なお以下の説明では、静磁場発生部２、傾斜磁場発生部３、シーケンサ４、送信部５、及び受信部６を総括して撮像部１００ともいう。   Next, an overview of imaging and signal processing of the MRI apparatus configured as described above will be described with reference to the flow of FIG. In the following description, the static magnetic field generation unit 2, the gradient magnetic field generation unit 3, the sequencer 4, the transmission unit 5, and the reception unit 6 are collectively referred to as an imaging unit 100.

まずＭＲＩ装置の静磁場空間に被検体１を配置し、撮像部１００を駆動し撮像を行う。撮像は位置決め等のためのリファレンススキャン（Ｓ１０１）と、本撮像（Ｓ１０３）とを含み、少なくともリファレンススキャンは、多チャンネルコイル、例えば全身用コイルを用いる。信号処理部７は、リファレンススキャンで取得したＮＭＲ信号を用いて参照画像を作成し（Ｓ１０３）、本撮像により得られたＮＭＲ信号を用いて、目的とする画像を作成する（Ｓ１０４）。目的とする画像には、プロトン密度画像、拡散強調画像等の強調画像、それら画像から派生する計算値画像などが含まれる。画像作成ステップＳ１０４には、予め取得或いは算出した、受信コイルの感度マップを用いた感度むらの補正（以下、単に感度補正という）（Ｓ１０４２）が含まれる。   First, the subject 1 is placed in the static magnetic field space of the MRI apparatus, and the imaging unit 100 is driven to perform imaging. Imaging includes a reference scan (S101) for positioning and the like, and main imaging (S103), and at least the reference scan uses a multi-channel coil, for example, a whole body coil. The signal processing unit 7 creates a reference image using the NMR signal acquired by the reference scan (S103), and creates a target image using the NMR signal obtained by the main imaging (S104). The target images include enhanced images such as proton density images and diffusion enhanced images, and calculated value images derived from these images. The image creation step S104 includes correction of sensitivity unevenness (hereinafter simply referred to as sensitivity correction) (S1042) using the sensitivity map of the receiving coil that is acquired or calculated in advance.

以下、詳述する実施形態のＭＲＩ装置では、受信コイルの感度マップの算出において、予め取得した多チャンネルコイルの感度マップとリファレンススキャンのデータが用いられる。具体的には、各チャンネルの感度マップとリファレンススキャンで得られた各チャンネルの参照画像を用いて座標変換式を導出し、この座標変換式を適用して、被検体存在時の受信コイルの受信感度マップを作成する（Ｓ１０４１）。信号処理部７は、この受信感度マップを用いて感度補正（Ｓ１０４２）や画像再構成に必要な演算を行い、画像を作成する。   Hereinafter, in the MRI apparatus of the embodiment described in detail, the sensitivity map of the multi-channel coil and the reference scan data acquired in advance are used in the calculation of the sensitivity map of the receiving coil. Specifically, a coordinate transformation formula is derived using the sensitivity map of each channel and the reference image of each channel obtained by the reference scan, and this coordinate transformation formula is applied to receive the receiving coil when the subject exists. A sensitivity map is created (S1041). The signal processing unit 7 performs an operation necessary for sensitivity correction (S1042) and image reconstruction using the reception sensitivity map to create an image.

上述した信号処理部７の計算部として機能は、ＣＰＵ８に備えられたメモリ（記憶部）やＣＰＵ８に組み込まれたプログラム等のソフトウェアを実行することで実現することができる。またその一部の機能をＡＳＩＣやＦＰＧＡなどのハードウェアで実現してもよい。
以上の概要を踏まえ、以下、信号処理部７が行う処理の各実施形態を説明する。 The function as the calculation unit of the signal processing unit 7 described above can be realized by executing software such as a memory (storage unit) provided in the CPU 8 or a program incorporated in the CPU 8. A part of the functions may be realized by hardware such as ASIC or FPGA.
Based on the above outline, each embodiment of the processing performed by the signal processing unit 7 will be described below.

＜第一実施形態＞
本実施形態は、多チャンネルコイルからで得たデータを用い座標変換式を作成し、被検体存在時の感度マップを作成する。本実施形態は、さらに、多チャンネルコイルについて作成された感度マップを用いて、本撮像に用いられた表面コイルの感度を補正する。 <First embodiment>
In the present embodiment, a coordinate conversion equation is created using data obtained from a multi-channel coil, and a sensitivity map when the subject is present is created. In the present embodiment, the sensitivity of the surface coil used for the main imaging is corrected using the sensitivity map created for the multi-channel coil.

本実施形態における信号処理部７は、図３の機能ブロック図に示すように、受信コイルの感度情報に関する種々の計算を行う感度計算部３０と、ＮＭＲ信号から画像を作成する画像作成部４０とを含む。感度計算部３０は、さらに、あらかじめ取得した多チャンネルコイルの各チャンネルに対する感度マップに相当する情報を保存する感度マップ保存部（記憶部）３１と、各チャンネルの感度マップの分布の違いを被検体に対するリファレンススキャンで得られる参照画像に内在する感度マップの分布の違いと揃えるような座標変換を決定する感度マップ座標変換決定部３２と、あらかじめ取得したコイルの感度マップをその座標変換で変換することで被検体存在時の感度マップ（予め取得した感度マップと区別して受信感度マップと呼ぶ）を作成する感度マップ作成部３３と、を備える。画像作成部４０には、感度マップ作成部３３が作成した受信感度マップを用いて感度補正を行う感度補正部４１が含まれる。感度マップ座標変換決定部３２は、指標算出部３２１と座標変換式算出部３２２を含む。   As shown in the functional block diagram of FIG. 3, the signal processing unit 7 in the present embodiment includes a sensitivity calculation unit 30 that performs various calculations related to sensitivity information of the receiving coil, an image generation unit 40 that generates an image from the NMR signal, including. The sensitivity calculation unit 30 further compares the sensitivity map storage unit (storage unit) 31 that stores information corresponding to the sensitivity map for each channel of the multi-channel coil acquired in advance, and the difference in sensitivity map distribution between the channels. A sensitivity map coordinate transformation determination unit 32 for determining coordinate transformation to align with the difference in distribution of sensitivity maps inherent in the reference image obtained by reference scanning with respect to the reference image, and transforming the previously acquired coil sensitivity map by the coordinate transformation And a sensitivity map creation unit 33 for creating a sensitivity map when the subject is present (referred to as a reception sensitivity map in distinction from a previously acquired sensitivity map). The image creation unit 40 includes a sensitivity correction unit 41 that performs sensitivity correction using the reception sensitivity map created by the sensitivity map creation unit 33. The sensitivity map coordinate conversion determination unit 32 includes an index calculation unit 321 and a coordinate conversion formula calculation unit 322.

感度マップ保存部３１に保存されている多チャンネルコイルの感度マップは、図４に示すように多チャンネルコイル全体の感度マップＳｗと、多チャンネルコイルを構成する各チャンネルの感度マップＳｎ（ｎはｎ番目のチャンネルであることを示し、チャンネル数を最大値とする整数）を含む。これら感度分布は、本撮像の対象である被検体とは異なるファントムなどで取得または計算機シミュレーションなどで作成されたものであり、被検体がセットされた時の感度マップとは異なる。   The sensitivity map of the multi-channel coil stored in the sensitivity map storage unit 31 includes a sensitivity map Sw for the entire multi-channel coil and a sensitivity map Sn (n is n for each channel constituting the multi-channel coil) as shown in FIG. Indicating the first channel and an integer with the maximum number of channels). These sensitivity distributions are obtained by a phantom or the like different from the subject that is the target of the main imaging or created by computer simulation or the like, and are different from the sensitivity map when the subject is set.

次に信号処理部７の処理の詳細を説明する。図４は信号処理部７の処理のフローを示す図であり、ここでは、前提として図２のステップＳ１０１において、感度マップが感度マップ保存部３０に保存されている多チャンネルコイルを用いたリファレンススキャンが行われ、リファレンススキャン（Ｓ１０１）により各チャンネルで得られたＮＭＲ信号から参照画像が作成されているものとする。   Next, details of the processing of the signal processing unit 7 will be described. FIG. 4 is a diagram showing a processing flow of the signal processing unit 7. Here, as a premise, a reference scan using a multi-channel coil in which the sensitivity map is stored in the sensitivity map storage unit 30 in step S101 of FIG. It is assumed that a reference image is created from the NMR signal obtained in each channel by the reference scan (S101).

信号処理部７（感度マップ座標変換決定部３２）は、この参照画像を用いて、各チャンネルの感度マップＳｎの分布の違いを被検体に対するリファレンススキャンと揃えるような座標変換を決定する。具体的には、まず、指標算出部３２１が、感度マップ保存部３１に保存されている各チャンネルの感度マップの対応画素間の関係（第一の指標マップ）を算出する（Ｓ４０１）。対応画素間の関係は、各チャンネルの感度分布の関係の指標となるものであり、例えば、各チャンネルの画素値の比、画素値をベクトル要素とするベクトル、当該ベクトルの偏角などを用いることができる。感度マップの全画素について指標を算出することで、指標を画素値とするマップ（指標マップ）が得られる。なお図４では第一の指標マップを算出するステップＳ４０１を、参照画像作成後に行っているものとして記載しているが、第一の指標マップは各チャンネルの感度マップがあれば算出できるので、撮像に先立って指標算出部３２１が算出し、予め感度マップ保存部３０に保存しておいてもよい。   The signal processing unit 7 (sensitivity map coordinate conversion determination unit 32) uses this reference image to determine coordinate conversion that aligns the difference in sensitivity map Sn distribution of each channel with the reference scan for the subject. Specifically, first, the index calculation unit 321 calculates the relationship (first index map) between corresponding pixels of the sensitivity map of each channel stored in the sensitivity map storage unit 31 (S401). The relationship between corresponding pixels serves as an index of the sensitivity distribution relationship of each channel. For example, a ratio of pixel values of each channel, a vector having pixel values as vector elements, a declination of the vector, and the like are used. Can do. By calculating the index for all the pixels of the sensitivity map, a map (index map) having the index as a pixel value is obtained. In FIG. 4, step S401 for calculating the first index map is described as being performed after the reference image is created. However, since the first index map can be calculated if there is a sensitivity map for each channel, imaging is performed. Prior to the calculation, the index calculation unit 321 may calculate and store the sensitivity map storage unit 30 in advance.

指標算出部３２１は、Ｓ１０２で作成された各チャンネルの参照画像を用いて、各チャンネルの参照画像の対応する画素間の関係（第二の指標マップ）を算出する（Ｓ４０２）。ここで用いる指標は、第一の指標マップの指標と同じものとする。即ち、第一の指標がベクトルであればベクトル、ベクトルの偏角であればベクトルの偏角とする。続いて座標変換式算出部３２２が、Ｓ４０１で算出した第一の指標マップを、Ｓ４０２で算出した第二の指標マップに一致させる座標変換式を求める（Ｓ４０３）。   The index calculation unit 321 calculates the relationship (second index map) between the corresponding pixels of the reference image of each channel using the reference image of each channel created in S102 (S402). The index used here is the same as the index of the first index map. That is, if the first index is a vector, the vector is used. If the first index is a vector, the vector is used. Subsequently, the coordinate conversion formula calculation unit 322 obtains a coordinate conversion formula that matches the first index map calculated in S401 with the second index map calculated in S402 (S403).

座標変換式算出方法について、２チャンネルコイルの場合を例にして説明する。図５に座標変換式算出処理を模式的に示す。感度マップ保存部３０に保存されているチャンネル１（ｃｈ１）の感度マップをＳ１（ｒ）、チャンネル２（ｃｈ２）の感度マップをＳ２（ｒ）とする。「ｒ」は、画素の座標を示す。これら各チャンネルの画素値を要素とするベクトルＳ１２（ここでは２次元ベクトル）の偏角φ１２（ｒ）を算出する。偏角φは画素毎に求められる。これを感度マップの対応画素間の関係を表す偏角分布（第一の指標マップ）とする。 The coordinate conversion formula calculation method will be described taking the case of a two-channel coil as an example. FIG. 5 schematically shows the coordinate conversion formula calculation process. The sensitivity map of channel 1 (ch1) stored in the sensitivity map storage unit 30 is S 1 (r), and the sensitivity map of channel 2 (ch2) is S 2 (r). “R” indicates the coordinates of the pixel. A deviation angle φ 12 (r) of a vector S 12 (two-dimensional vector here) having the pixel value of each channel as an element is calculated. The declination φ is obtained for each pixel. This is a declination distribution (first index map) representing the relationship between corresponding pixels in the sensitivity map.

またチャンネル１（ｃｈ１）の参照画像をＲ１（ｒ）、チャンネル２（ｃｈ２）の参照画像をＲ２（ｒ）とし、各チャンネルの画素値を要素とするベクトルＲ１２の偏角θ１２（ｒ）を算出する。偏角θは画素毎に求められる。これを感度マップの画素間の関係を表す偏角分布（第二の指標マップ）とする。 The Channel 1 (ch1) of the reference image R 1 (r), the reference picture of the channel 2 (ch2) and R 2 (r), the deflection angle theta 12 vector R 12 to the pixel value of each channel element ( r) is calculated. The declination θ is obtained for each pixel. This is the declination distribution (second index map) representing the relationship between the pixels of the sensitivity map.

次に次式に従い、偏角φと偏角θとの差を最小とする座標変換式を求める。具体的には、例えばｘｙ平面での平行移動と拡大縮小の場合、平行移動と拡大縮小の４つのパラメータ（ｘ方向及びｙ方向の移動量、ｘ方向及びｙ方向の拡大率）を１画素に相当する量、変化させながら探索を行い、式（１）を最小とする位置を決定する。
［数１］
Σ｜（φ１２（ｒ）−θ１２（ｒ’））２ ｜ （１） Next, a coordinate conversion equation that minimizes the difference between the deflection angle φ and the deflection angle θ is obtained according to the following equation. Specifically, for example, in the case of parallel movement and enlargement / reduction on the xy plane, the four parameters of parallel movement and enlargement / reduction (the amount of movement in the x and y directions and the enlargement ratio in the x and y directions) are set to one pixel. A search is performed while changing the corresponding amount, and the position that minimizes Equation (1) is determined.
[Equation 1]
Σ | (φ 12 (r) −θ 12 (r ′)) 2 | (1)

座標変換式は、例えば式（２）或いは式（２’）で表すことができる。

Figure 2018093989
式（２’）はｒをｘ−ｙ座標で表した時のアフィン変換を示し、拡大・縮小、平行移動、及び回転を含む。 The coordinate conversion expression can be expressed by, for example, Expression (2) or Expression (2 ′).
Figure 2018093989
 Expression (2 ′) represents an affine transformation when r is expressed by xy coordinates, and includes enlargement / reduction, translation, and rotation.

以上の説明は、２チャンネルの例であるが、３チャンネル以上でもよく、ｎ（ｎは２以上の整数）チャンネルの場合は、各チャンネルの感度を要素とするｎ次元ベクトルの成す角が０になるように座標変換をすれば良い。その際、チャンネル感度の全体的なずれを吸収するために、各チャンネルの感度にかける係数も変化させて成す角を０に近づけても良い。また上述した座標変換は一例であり、その他、座標変換には公知のレジストレーション技術、剛体レジストレーションや非剛体レジストレーションなどを採用してもよい。例えばチャンネル間の画素値の比やベクトルの偏角を求めてからレジストレーションを行ってもよいし、ベクトルのまま複数チャンネルを持つデータとしてレジストレーションを行ってもよい。さらにベクトルの大きさを規格化してからレジストレーションを行うなど、対象とするデータ形式も任意に変更できる。   The above explanation is an example of 2 channels, but 3 channels or more may be used. In the case of n (n is an integer of 2 or more) channels, the angle formed by an n-dimensional vector having the sensitivity of each channel as an element is 0. Coordinate conversion may be performed so that At this time, in order to absorb the overall shift in channel sensitivity, the angle formed by changing the coefficient applied to the sensitivity of each channel may be close to zero. The coordinate transformation described above is an example, and other known registration techniques, rigid registration, non-rigid registration, and the like may be employed for coordinate transformation. For example, registration may be performed after obtaining the ratio of pixel values between channels or vector declination, or registration as data having a plurality of channels as vectors. Furthermore, the target data format can be arbitrarily changed, such as performing registration after normalizing the vector size.

ステップＳ４０３により、座標変換式Ｆ（ｒ→ｒ’）が定まったならば、この座標変換式を用いて、感度マップ保存部３０に保存されている多チャンネルコイルの感度マップＳｗを座標変換し、被検体の受信感度を含む感度マップＳ’ｗ（受信感度マップ）を得る（Ｓ４０４）。   If the coordinate conversion formula F (r → r ′) is determined in step S403, the multi-channel coil sensitivity map Sw stored in the sensitivity map storage unit 30 is coordinate-converted using this coordinate conversion formula, A sensitivity map S′w (reception sensitivity map) including the reception sensitivity of the subject is obtained (S404).

次に、表面コイルについて予め取得した座標変換式を導出した多チャンネルコイルの感度マップに対する比としての感度マップＣに、Ｓ４０４で得た受信感度マップＳ’ｗをかけることで表面コイルの感度マップを求め、本撮像で取得した画像を補正する（Ｓ４０５：感度補正）。   Next, the sensitivity map of the surface coil is obtained by multiplying the sensitivity map C as a ratio with respect to the sensitivity map of the multi-channel coil from which the coordinate transformation formula obtained in advance for the surface coil is derived by applying the reception sensitivity map S′w obtained in S404. Then, the image acquired by the main imaging is corrected (S405: sensitivity correction).

本撮像で得た画像Ｉの感度補正（Ｓ４０５）は、この補正後感度マップＳ’ｗと、本撮像に用いた受信コイルの座標変換式を導出した多チャンネルコイルの感度マップに対する比としての感度マップＣとを用いて次式により行う。
［数３］
Ｉ’（ｘ、ｙ）＝Ｉ（ｘ、ｙ）／（Ｃ×Ｓ’ｗ） （３） Sensitivity correction (S405) of the image I obtained by the main imaging is the sensitivity as a ratio between the sensitivity map S′w after correction and the sensitivity map of the multi-channel coil derived from the coordinate conversion equation of the receiving coil used for the main imaging. Using the map C, the following equation is used.
[Equation 3]
I ′ (x, y) = I (x, y) / (C × S′w) (3)

なお式（３）における「Ｃ×Ｓ’ｗ」は、式（３）の計算の前に算出しておいてもよいし、一度に式（３）補正後の画像Ｉ’ を算出してもよい。また式（３）が発散しないために、右辺の分母をゼロにしないための処理、例えば感度がゼロになる領域を計算から排除する処理や一定値を挿入する処理などの公知技術、を行ってもよい。   Note that “C × S′w” in the equation (3) may be calculated before the calculation of the equation (3), or the corrected image I ′ may be calculated at once in the equation (3). Good. Further, since the expression (3) does not diverge, a known technique such as processing for preventing the right side denominator from being zero, for example, processing for excluding a region where sensitivity is zero or processing for inserting a constant value is performed. Also good.

上述の例では本撮像には表面コイルを用い、座標変換式を導出した多チャンネルコイルの感度マップに対する比としての表面コイルの感度マップＣをリファレンススキャンで求めたが、本撮像に用いた受信コイルが、座標変換式を導出した多チャンネルコイルと同一の受信コイルの場合には、図６のように、Ｓ４０４で得た受信感度マップＳ’ｗを用いて、従来の感度補正と同様に本撮像で得た画像の感度補正を行う（Ｓ４０６）。   In the above-described example, the surface coil is used for the main imaging, and the sensitivity map C of the surface coil as a ratio to the sensitivity map of the multi-channel coil from which the coordinate conversion formula is derived is obtained by the reference scan. However, in the case of the same receiving coil as the multi-channel coil from which the coordinate conversion formula is derived, as shown in FIG. 6, using the receiving sensitivity map S′w obtained in S404, the main imaging is performed as in the conventional sensitivity correction. The sensitivity correction of the image obtained in step S406 is performed (S406).

以上の例では、感度マップ各点での偏角の差の２乗和を最小にしたが、式（３）のかわりに式（４）のように比較する領域で平均値がそろうようにオフセットを加えてから感度マップ各点での偏角の差の２乗和を最小にする方法などを採用してもよい。
［数４］
Σ｜((φ１２(ｒ)−Avg(φ１２(ｒ)))−(θ１２(ｒ’)−Avg(θ１２(ｒ’)))）２｜
（４）
ここで、Avgは比較する領域での平均を表す。 In the above example, the sum of squares of the difference in the declination at each point of the sensitivity map is minimized, but the offset is set so that the average value is equal in the area to be compared as in equation (4) instead of equation (3). Or the like, a method of minimizing the sum of squares of the difference in declination at each point of the sensitivity map may be employed.
[Equation 4]
Σ | ((φ 12 (r) −Avg (φ 12 (r))) − (θ 12 (r ′) − Avg (θ 12 (r ′)))) 2 |
(4)
Here, Avg represents the average in the area to be compared.

また上述の例では、各チャンネルの感度マップ間或いは参照画像間の関係性を示す指標として偏角を用いたが、二つの指標（第一の指標と第二の指標）が同じであれば、偏角に限らず、例えば、各チャンネルの感度マップ或いは参照画像の画素値の比、差分など別の指標を用いることも可能である。   In the above example, the declination is used as an index indicating the relationship between the sensitivity maps of each channel or between the reference images. If the two indexes (the first index and the second index) are the same, For example, another index such as a sensitivity map of each channel or a ratio or difference of pixel values of a reference image can be used.

次に、第一実施形態による補正の効果について図７を用いて説明する。図７はリファレンススキャンの撮像条件をコントラストが生じるような条件にした頭部の撮像結果で、下側は撮像した画像であり、上側は画像上のプロファイルライン上のプロファイルを示す図である。また図７において、左（ａ）が補正なしの画像、中央（ｂ）は従来法（ボディーコイルの参照画像にローパスフィルターをかけたものをボディーコイルの感度マップとする方法）で感度補正した画像、右側（ｃ）は第一実施形態の手法（あらかじめファントムで取得した感度マップを変形してボディーコイルの感度マップとする方法）で感度補正した画像である。   Next, the effect of the correction according to the first embodiment will be described with reference to FIG. FIG. 7 shows the imaging result of the head under the condition that the imaging condition of the reference scan is such that a contrast is generated. The lower side is a captured image, and the upper side is a diagram showing a profile on a profile line on the image. In FIG. 7, the left (a) is an uncorrected image, and the center (b) is an image whose sensitivity is corrected by a conventional method (a method in which a body coil reference image is subjected to a low-pass filter is used as a body coil sensitivity map). The right side (c) is an image whose sensitivity has been corrected by the method of the first embodiment (a method in which a sensitivity map acquired in advance with a phantom is transformed into a body coil sensitivity map).

これらの結果からわかるように、補正なしの結果（ａ）は、ラインプロファイルを見ると本来左右対称のところ、輝度値に差があることがわかる。従来法の結果（ｂ）では、リファレンススキャンで生じた影響が補正されず、ラインプロファイルを見ると左右の差が残っていることがわかる。一方、本実施形態の結果（ｃ）では、従来法（ｂ）と同じリファレンススキャンを使用しているが、左右の差がなくなっていることがわかる。   As can be seen from these results, the result (a) without correction shows that there is a difference in luminance value when the line profile is essentially symmetrical. In the result (b) of the conventional method, the effect caused by the reference scan is not corrected, and it can be seen that the difference between the left and right remains when the line profile is viewed. On the other hand, in the result (c) of this embodiment, the same reference scan as in the conventional method (b) is used, but it can be seen that there is no difference between left and right.

本実施形態によれば、多チャンネルコイルの各チャンネルの感度分布間の関係性を示す第一の指標と、同コイルの各チャンネルで得た参照画像間の関係性を示す第二の指標とを一致させる座標変換式を求めておき、当該座標変換式を用いて被検体の感度を含む受信コイルの受信感度マップを作成する。これにより、感度補正の精度を向上することができる。   According to the present embodiment, the first index indicating the relationship between the sensitivity distributions of each channel of the multi-channel coil and the second index indicating the relationship between the reference images obtained in each channel of the coil. A coordinate conversion formula to be matched is obtained, and a reception sensitivity map of the receiving coil including the sensitivity of the subject is created using the coordinate conversion formula. Thereby, the accuracy of sensitivity correction can be improved.

また本実施形態によれば、多チャンネルコイルの感度マップを作成する際に、感度マップがリファレンススキャンのコントラストに影響されにくくなるという点がある。なぜなら、例で用いた指標のチャンネル間の比ではコントラストが割り算により消えるというように、コントラストを打ち消すような指標をとれるからである。従ってコントラストを気にせず、リファレンススキャンの撮像パラメータを設定できるようになるので、リファレンススキャンのＴＥを調整して撮像音を低減するということも可能になる。   In addition, according to the present embodiment, when creating a sensitivity map of a multi-channel coil, there is a point that the sensitivity map is not easily influenced by the contrast of the reference scan. This is because an index that cancels the contrast can be taken such that the contrast disappears by division in the ratio between the channels of the index used in the example. Accordingly, the imaging parameter of the reference scan can be set without worrying about the contrast, and it is also possible to reduce the imaging sound by adjusting the TE of the reference scan.

なお以上の実施形態では、画像全体を座標変換式により補正された受信コイルの感度マップ又はそれにより補正された別の受信コイルの感度マップを用いて感度補正を行う例を説明したが、他の感度補正と組み合わせても良い。たとえば、２Ｄ撮像でも３Ｄ撮像でも、体軸に垂直な方向の感度マップとして本実施形態の被検体感度を含む受信感度マップを用いて感度補正を行い、体軸方向の感度補正は従来の感度補正で行っても良い。   In the above embodiment, an example has been described in which sensitivity correction is performed using the sensitivity map of the receiving coil corrected by the coordinate conversion formula for the entire image or the sensitivity map of another receiving coil corrected thereby. It may be combined with sensitivity correction. For example, in both 2D imaging and 3D imaging, sensitivity correction is performed using the reception sensitivity map including the subject sensitivity of the present embodiment as a sensitivity map in the direction perpendicular to the body axis, and the sensitivity correction in the body axis direction is the conventional sensitivity correction. You can go there.

また以上の実施形態では、受信感度補正について説明したが、従来の受信感度補正技術と同様に、送信感度補正やフィルター処理などさらに他の補正や画像処理を行っても良いことは言うまでもない。   In the above embodiment, the reception sensitivity correction has been described, but it goes without saying that other corrections and image processing such as transmission sensitivity correction and filter processing may be performed as in the conventional reception sensitivity correction technique.

＜第二実施形態＞
本実施形態においても、感度マップ保存部３１に保存された各受信チャンネルの感度マップから算出した第一指標と、各チャンネルの参照画像を用いて算出した第二指標とを一致させる座標変換式を用いて、基準コイルの感度分布を補正して、感度補正に用いることは第一実施形態と同じである。本実施形態は、その前提である多チャンネル受信コイルの感度分布（感度マップ保存部に格納された感度マップ）の取得手法に特徴がある。 <Second embodiment>
Also in the present embodiment, a coordinate conversion equation for matching the first index calculated from the sensitivity map of each reception channel stored in the sensitivity map storage unit 31 with the second index calculated using the reference image of each channel is used. It is the same as in the first embodiment that the sensitivity distribution of the reference coil is corrected and used for sensitivity correction. The present embodiment is characterized by a technique for obtaining the sensitivity distribution (sensitivity map stored in the sensitivity map storage unit) of the multi-channel receiving coil, which is the premise thereof.

感度分布の取得手法には、大きく分けて、ファントムを用いる方法と、シミュレーションによる方法とがあり、いずれの方法を採用してもよい。本実施形態では、シミュレーションで感度マップを作成し、それをファントム撮像で得た感度マップで補正する手法を説明する。   The sensitivity distribution acquisition method is roughly classified into a method using a phantom and a method using simulation, and either method may be adopted. In the present embodiment, a method of creating a sensitivity map by simulation and correcting it with a sensitivity map obtained by phantom imaging will be described.

本実施形態における信号処理部７の、主として感度計算部の機能ブロック図を図８に示す。図示するように感度計算部３０は、感度マップ保存部３１、感度マップ座標変換決定部３２及び感度マップ作成部３３に加えて、受信コイルの初期感度マップ（被検体が存在しないときの感度マップ）を算出するシミュレーション部３５と、シミュレーション部３５が算出した初期感度マップをファントム撮像で得た画像に合わせる初期感度マップ調整部３６とを有する。   FIG. 8 shows a functional block diagram mainly of the sensitivity calculation unit of the signal processing unit 7 in the present embodiment. As shown in the figure, the sensitivity calculation unit 30 includes an initial sensitivity map of the receiving coil (sensitivity map when no subject exists) in addition to the sensitivity map storage unit 31, the sensitivity map coordinate conversion determination unit 32, and the sensitivity map creation unit 33. And an initial sensitivity map adjustment unit 36 that matches the initial sensitivity map calculated by the simulation unit 35 with an image obtained by phantom imaging.

以下、図９を参照して、本実施形態の信号処理部７の処理を説明する。
まずシミュレーション部３５は、感度マップをＦＤＴＤ法など公知のシミュレーション技術により算出する（Ｓ８０１）。シミュレーションでは、ファントムのサイズ、物性値、対象とする受信コイルのサイズ、キャパシタやインダクタの値、素子間のサイズ等を設定し、磁気共鳴周波数の信号を流したときの感度分布（感度マップ）を算出する。本実施形態では、図１０に示すように、大きな被検体にも対応できるように、シミュレーションではファントムで得られる感度マップＤ（図１０（ｂ））よりも大きなファントムを設定することもできる。これにより、図１０（ａ）に示すような初期感度マップＳを得る。 Hereinafter, the processing of the signal processing unit 7 of the present embodiment will be described with reference to FIG.
First, the simulation unit 35 calculates a sensitivity map by a known simulation technique such as the FDTD method (S801). In the simulation, the sensitivity distribution (sensitivity map) when the magnetic resonance frequency signal is sent is set by setting the size of the phantom, the physical property value, the size of the target receiving coil, the value of the capacitor or inductor, the size between elements, etc. calculate. In the present embodiment, as shown in FIG. 10, a phantom larger than the sensitivity map D (FIG. 10B) obtained by the phantom can be set in the simulation so that a large subject can be handled. Thereby, an initial sensitivity map S as shown in FIG.

次に初期感度マップ調整部３６が、シミュレーション結果である初期感度マップＳを、実際にファントムを撮像して得た感度マップＤに近づける修正を行う（Ｓ８０２）。ここで、ファントム撮像はアディアバティックパルスを用いて送信不均一の影響を取り除くなどして受信感度だけを取り出せるようにする。Ｂ１マップを求めて送信不均一の影響を取り除いても良い。また、表面コイルを置いて表面コイルの影響を含めた受信感度にしても良い。   Next, the initial sensitivity map adjustment unit 36 corrects the initial sensitivity map S, which is a simulation result, to be close to the sensitivity map D obtained by actually imaging the phantom (S802). Here, in phantom imaging, only the reception sensitivity can be extracted by removing the influence of transmission non-uniformity using an adiabatic pulse. The influence of non-uniform transmission may be removed by obtaining the B1 map. Alternatively, a reception coil including the influence of the surface coil may be provided by placing a surface coil.

初期感度マップＳをファントム撮像の感度マップＤに近づける方法として、具体的には、各チャンネルの初期感度マップに対し式（５）による座標変換（回転、平行移動、拡大縮小）と、チャンネルｉ毎に異なる定数Ｋｉの掛け算とを繰り返し、Ｓ’ｉ（ｒ）とＤ（ｒ）との差を最小化する座標変換式及び定数を算出する。ここで初期感度マップとファントム画像との差分は、それぞれ被検体（ファントム）が存在する領域について比較する。
［数５］
Ｓ’ｉ（ｒ）＝Ki×Ｓｉ(f（ｒ）） （５）
式中、ｆは、座標変換式ｆ（ｒ→ｒ’）を表す。 As a method for bringing the initial sensitivity map S closer to the sensitivity map D for phantom imaging, specifically, coordinate conversion (rotation, parallel movement, enlargement / reduction) according to the equation (5) is performed on the initial sensitivity map of each channel, and for each channel i. And a multiplication of different constants Ki is repeated to calculate a coordinate conversion formula and constants that minimize the difference between S′i (r) and D (r). Here, the difference between the initial sensitivity map and the phantom image is compared for each region where the subject (phantom) exists.
[Equation 5]
S′i (r) = Ki × Si (f (r)) (5)
In the formula, f represents a coordinate conversion formula f (r → r ′).

最小化の手法は、シミュレーションで求めた初期感度マップとファントム画像との差の２乗の平均を用いる方法など公知の手法を用いることができるが、ここでは両者を厳密に一致させる必要はなく、おおまかな修正でよい。   As a minimization method, a known method such as a method using an average of the square of the difference between the initial sensitivity map obtained by simulation and the phantom image can be used, but here it is not necessary to strictly match both, Rough correction is sufficient.

最後に大まかに修正したシミュレーション感度マップＳ’（ｒ）でファントム撮像画像Ｄを割って係数（画素毎の係数）βを求め、この分布を持つ係数βを修正後のシミュレーション感度Ｓ’（ｒ）に掛ける（Ｓ８０３）。この際、係数βの分布を滑らかにするためにメディアンフィルターやローパスフィルターなどで平滑化する。またファントム撮像画像Ｄにおいて被検体（ファントム）が存在しない領域の画素については、被検体が存在する領域で算出した係数βを外挿する（Ｓ８０４）。
外挿は公知の技術を用いることができるが、係数βは、シミュレーションが完全であれば恒常的に１となる性質のものであるため、単純な処理で安定に外挿できる。例えば、外挿すべき点に最も距離が近い位置の点の係数で置き換えるという処理でも良い。 Finally, the coefficient (pixel-specific coefficient) β is obtained by dividing the phantom captured image D by the roughly corrected simulation sensitivity map S ′ (r), and the coefficient β having this distribution is corrected to the simulation sensitivity S ′ (r). (S803). At this time, in order to smooth the distribution of the coefficient β, it is smoothed by a median filter, a low-pass filter, or the like. For the pixels in the region where the subject (phantom) does not exist in the phantom captured image D, the coefficient β calculated in the region where the subject exists is extrapolated (S804).
A known technique can be used for extrapolation. However, since the coefficient β has a property of constantly being 1 if the simulation is complete, it can be extrapolated stably by simple processing. For example, a process of replacing with a coefficient of a point closest to the point to be extrapolated may be performed.

以上の処理により、シミュレーションで求めた初期感度マップの修正が完了する。修正された感度マップ（図１０（ｃ））は、第一実施形態で説明した受信感度マップ（被検体に対する感度マップ）の作成に用いられる。   With the above processing, the correction of the initial sensitivity map obtained by the simulation is completed. The corrected sensitivity map (FIG. 10C) is used to create the reception sensitivity map (sensitivity map for the subject) described in the first embodiment.

本実施形態のＭＲＩ装置は、記憶部（感度マップ保存部）に記憶する複数の受信チャンネルの感度マップが、シミュレーションで得た感度画像をファントム画像で補正した感度マップである。本実施形態によれば、感度マップ作成に用いる被検体以外での感度マップをファントムサイズや物性値などを自由に設定して作成できる。   In the MRI apparatus of this embodiment, the sensitivity maps of a plurality of reception channels stored in the storage unit (sensitivity map storage unit) are sensitivity maps obtained by correcting sensitivity images obtained by simulation with phantom images. According to this embodiment, it is possible to create a sensitivity map other than the subject used for creating the sensitivity map by freely setting the phantom size, the physical property value, and the like.

なお以上の説明では、初期感度マップをファントム撮像結果におおまかに近づくように修正する処理Ｓ８０２で、チャンネルごとに同じ座標変換（回転と平行移動と拡大縮小）と、チャンネルごとに異なる定数の掛け算を行ったが、この方法に限らない。たとえば、チャンネルごとに異なる定数の掛け算のみを行っても良いし、定数の掛け算を省略してもよい。   In the above description, the same coordinate transformation (rotation, parallel movement, and enlargement / reduction) is performed for each channel, and multiplication of constants that are different for each channel is performed in the process S802 for correcting the initial sensitivity map so as to approximate the phantom imaging result. Although it went, it is not restricted to this method. For example, only different constant multiplication may be performed for each channel, or constant multiplication may be omitted.

また以上の説明では、シミュレーションで大きなファントムを設定して、ファントム撮像結果に合うように補正したが、ファントム撮像と同じファントムのシミュレーション結果から、補正（回転、拡大縮小、平行移動、増幅、係数）を求めたうえで、その補正を大きなファントムのシミュレーション結果に適用して感度マップ（感度マップ保存部に格納すべき感度マップ）を作成してもよい。   In the above description, a large phantom was set in the simulation and corrected to match the phantom imaging result. However, correction (rotation, enlargement / reduction, translation, amplification, coefficient) was made based on the same phantom simulation result as phantom imaging. Then, the sensitivity map (sensitivity map to be stored in the sensitivity map storage unit) may be created by applying the correction to the simulation result of a large phantom.

また、以上の説明では一つの感度マップ（感度マップ保存部に格納すべき感度マップ）を作成する場合を説明したが、被検体の条件（位置やサイズなど）に応じて使い分けられるように、条件の異なるあらかじめ感度マップを複数用意しても良い。   In the above description, the case where one sensitivity map (sensitivity map to be stored in the sensitivity map storage unit) is created has been described. However, it is necessary to use different conditions depending on the condition (position, size, etc.) of the subject. A plurality of sensitivity maps may be prepared in advance.

＜第三実施形態＞
本実施形態では、撮像用の受信コイルの感度マップとして位相情報を含む感度マップを作成することが特徴である。位相情報を含む感度マップを得ることで、感度むらの補正（いわゆる感度補正）のみならず、受信コイル感度マップをパラレルイメージングに利用することも可能となる。 <Third embodiment>
The present embodiment is characterized in that a sensitivity map including phase information is created as a sensitivity map of a receiving coil for imaging. By obtaining a sensitivity map including phase information, it is possible not only to correct sensitivity unevenness (so-called sensitivity correction) but also to use the reception coil sensitivity map for parallel imaging.

本実施形態では、感度マップ保存部３１は、複数チャンネルの受信コイルの各チャンネルの感度マップを複素数で保存している。例えば、均質なファントムを撮像することで受信コイルが得た信号は、信号強度と位相を含む複素データであるが、感度むらのみを補正する場合には、感度の絶対値である信号強度のみがあればよく、信号強度が感度マップである。しかし、本実施形態は、位相情報も利用することで、例えば画像再構成演算に、位相を含む受信感度情報を用いるパラレルイメージング等への適用を可能にする。   In the present embodiment, the sensitivity map storage unit 31 stores the sensitivity map of each channel of the reception coils of a plurality of channels as a complex number. For example, the signal obtained by the receiving coil by imaging a homogeneous phantom is complex data including the signal strength and phase, but when correcting only the sensitivity unevenness, only the signal strength that is the absolute value of the sensitivity is obtained. The signal intensity is a sensitivity map. However, this embodiment makes it possible to apply to, for example, parallel imaging using reception sensitivity information including a phase for image reconstruction calculation by using phase information as well.

本実施形態の感度マップ作成部が受信感度マップを作成する対象となる受信コイルは、別の撮像用表面コイルの感度マップを補正するための多チャンネルコイルでもよいが、好適には、それ自体が撮像用の多チャンネルコイルとなる。以下の説明では、撮像用の多チャンネルコイルである場合を説明する。感度マップ保存部３１以外の信号処理部７の構成は、図１１に示すように、第一実施形態において図３に例示した機能ブロック図の構成に、複素数を扱うために絶対値／位相補正部（位相補正部）３７が追加される。以下、図１２を参照して、第一実施形態と異なる点を中心に本実施形態の信号処理部７の動作を説明する。   The receiving coil for which the sensitivity map creating unit of the present embodiment creates the receiving sensitivity map may be a multi-channel coil for correcting the sensitivity map of another imaging surface coil, but it is preferably itself. It becomes a multi-channel coil for imaging. In the following description, a case of a multi-channel coil for imaging will be described. As shown in FIG. 11, the configuration of the signal processing unit 7 other than the sensitivity map storage unit 31 is the same as the configuration of the functional block diagram illustrated in FIG. 3 in the first embodiment as an absolute value / phase correction unit for handling complex numbers. A (phase correction unit) 37 is added. Hereinafter, with reference to FIG. 12, the operation of the signal processing unit 7 of the present embodiment will be described focusing on differences from the first embodiment.

まず座標変換決定処理Ｓ５０１では、リファレンススキャンによって、基準コイルの各チャンネルで得られた参照画像の絶対値（各画素の絶対値）をｎ次元ベクトルとし、感度マップ保存部３１に保存されている各チャンネルの感度分布の絶対値（各画素の絶対値）をｎ次元ベクトルとする。これら参照画像のベクトルと、感度マップのベクトルとがなす角が全体的に小さくなるように、感度マップを座標変換する座標変換式を算出する。この手法は、第一実施形態のＳ４０３と同様であり、拡大・縮小、平行移動、及び回転を含む。   First, in the coordinate conversion determination process S501, the absolute value of the reference image (absolute value of each pixel) obtained in each channel of the reference coil by the reference scan is set as an n-dimensional vector, and is stored in the sensitivity map storage unit 31. An absolute value of channel sensitivity distribution (absolute value of each pixel) is an n-dimensional vector. A coordinate conversion formula for coordinate conversion of the sensitivity map is calculated so that the angle formed by the reference image vector and the sensitivity map vector is reduced as a whole. This method is the same as S403 of the first embodiment, and includes enlargement / reduction, translation, and rotation.

次いでＳ５０１で得られた座標変換式を用いて、感度マップ保存部３１に保存されている各チャンネルの感度マップ（複素数）を座標変換する（Ｓ５０２）。次いで、座標変換後の感度マップの絶対値と位相値を補正する（Ｓ５０３）。このため、まず画素毎に、感度分布の絶対値と参照画像の絶対値から、式（６）を満たす係数ｋｉを算出する。   Next, using the coordinate conversion formula obtained in S501, the sensitivity map (complex number) of each channel stored in the sensitivity map storage unit 31 is subjected to coordinate conversion (S502). Next, the absolute value and phase value of the sensitivity map after coordinate conversion are corrected (S503). For this reason, first, for each pixel, a coefficient ki that satisfies Expression (6) is calculated from the absolute value of the sensitivity distribution and the absolute value of the reference image.

［数６］
｜Ｒｉ｜＝ｋｉ×｜Ｓｉ｜ （６）
式中、Ｓｉ（ｉは、ｉ番目のチャンネルであることを示し、チャンネル数を最大値とする整数）は座標変換されたチャンネルｉの感度マップ、Ｒｉは、チャンネルｉの参照画像である。この係数ｋｉを用いて、次式（７）から、各チャンネルの感度マップ（複素数）Ｒ’ｉを算出する。 [Equation 6]
| Ri | = ki × | Si | (6)
In the equation, Si (i indicates an i-th channel and an integer having the maximum number of channels) is a sensitivity map of channel i after coordinate conversion, and Ri is a reference image of channel i. Using this coefficient ki, the sensitivity map (complex number) R′i of each channel is calculated from the following equation (7).

［数７］
Ｒ’ｉ＝ｋｉ×Ｓｉ×exp(j×di)×（｜Ｓ｜／｜Ｒ｜） （７）
式中、Ｓは各チャンネルの感度マップの画素値を要素とするベクトル（チャンネル数がｎの場合、ｎ次元ベクトル）、Ｒは参照画像の画素値を要素とするベクトルである。jは虚数単位で、ｄiは、感度マップＲiの偏角（ａｒｇ（Ｒi））と参照画像Ｓiの偏角ａｒｇ（Ｓi）との差（位相差）で、式（８）で定義される。 [Equation 7]
R′i = ki × Si × exp (j × di) × (| S | / | R |) (7)
In the equation, S is a vector having a pixel value of the sensitivity map of each channel as an element (when the number of channels is n, an n-dimensional vector), and R is a vector having a pixel value of the reference image as an element. j is an imaginary unit, and di is a difference (phase difference) between the deviation angle (arg (Ri)) of the sensitivity map Ri and the deviation angle arg (Si) of the reference image Si, and is defined by Expression (8).

［数８］
ｄi＝ａｒｇ（Ｒi）−ａｒｇ（Ｓi） （８）
なお「ａｒｇ（Ｒｉ）」及び「ａｒｇ（Ｓｉ）は、いずれもアンラップ処理後の偏角を表す。
式（７）により得られる各チャンネルの感度マップＲ’ｉは、チャンネル間の比は参照画像におけるチャンネル間の比と同じであり、且つ感度マップＲｉを合成して得られる基準コイル全体の感度マップはフラットなものとなる。 [Equation 8]
di = arg (Ri) -arg (Si) (8)
Note that “arg (Ri)” and “arg (Si)” both represent declination angles after unwrapping.
The sensitivity map R′i of each channel obtained by Expression (7) has the same ratio between channels as the ratio between channels in the reference image, and the sensitivity map of the entire reference coil obtained by synthesizing the sensitivity map Ri. Will be flat.

各チャンネルの参照画像が、感度マップ保存部３１に保存されているチャンネル毎の感度マップよりも低解像度の場合或いは金属の存在などの理由で正確な参照画像が得られない領域がある場合、それら欠落した領域或いは足りない画素については、ｋｉとｄｉは、それぞれ補間によって求め、その後、式（７）よりＲ’ｉが求められる。これにより低分解能の参照画像から高分解能の感度マップを得たり、正確な参照画像が得られない領域の感度マップを得たりすることができる。   If the reference image of each channel has a lower resolution than the sensitivity map for each channel stored in the sensitivity map storage unit 31 or if there is an area where an accurate reference image cannot be obtained due to the presence of metal, etc. For missing regions or missing pixels, ki and di are obtained by interpolation, and then R′i is obtained from equation (7). As a result, a high-resolution sensitivity map can be obtained from a low-resolution reference image, or a sensitivity map in a region where an accurate reference image cannot be obtained can be obtained.

こうして絶対値及び位相が決定した多チャンネル受信コイルの感度マップは、パラレルイメージングの本撮像が行われたのであれば、当該多チャンネル受信コイルの各チャンネルで受信したＮＭＲ信号を用いた公知のパラレルイメージングに基く演算、例えばｋ空間データを用いた折り返し展開演算や実空間データを用いた折り返し除去演算、に用いることができる（Ｓ５０４）。   The sensitivity map of the multi-channel receiving coil whose absolute value and phase are determined in this way is known parallel imaging using NMR signals received by each channel of the multi-channel receiving coil if the main imaging of parallel imaging is performed. For example, a folding expansion calculation using k-space data and a folding removal calculation using real space data (S504).

また本撮像が、上述の処理で感度マップを算出した多チャンネル受信コイルとは別の表面コイルを用いた撮像であれば、第一実施形態と同様にこの表面コイルの感度分布を補正した上で、表面コイルで得た画像に対し感度補正を行ってもよい。   If the imaging is an imaging using a surface coil different from the multi-channel receiving coil for which the sensitivity map has been calculated by the above processing, the sensitivity distribution of the surface coil is corrected as in the first embodiment. Sensitivity correction may be performed on the image obtained by the surface coil.

本実施形態のＭＲＩ装置は、記憶部（感度マップ保存部）に記憶される複数の受信チャンネルの感度マップ、及び、複数の参照画像は、画素値が複素数であり、感度マップ作成部は、複数の受信チャンネルの感度マップ及び前記複数の参照画像の位相を揃える位相補正部をさらに備える。また計算部（感度計算部）は、感度マップ作成部により、位相が補正された、各受信チャンネルの受信感度マップを用いて、パラレルイメージングに基く演算を行い、画像を作成することができる。   In the MRI apparatus of the present embodiment, the sensitivity map of the plurality of reception channels stored in the storage unit (sensitivity map storage unit) and the plurality of reference images have complex pixel values, and the sensitivity map creation unit includes a plurality of sensitivity map creation units. And a phase correction unit for aligning the phases of the plurality of reference images. The calculation unit (sensitivity calculation unit) can perform an operation based on parallel imaging using the reception sensitivity map of each reception channel whose phase is corrected by the sensitivity map generation unit, and can generate an image.

本実施形態によれば、各チャンネルの感度マップを複素数として扱い、チャンネル間の関係が、被検体存在時の画像（参照画像）のチャンネル間の関係と同様になるように座標変換し、絶対値及び位相を持つ感度情報を得るので、精度のよい感度補正が可能であり、また感度情報を利用したパラレルイメージングも可能となる。   According to this embodiment, the sensitivity map of each channel is treated as a complex number, coordinate conversion is performed so that the relationship between channels is the same as the relationship between channels of an image (reference image) when the subject is present, and absolute values Since the sensitivity information having the phase and the phase is obtained, the sensitivity can be corrected with high accuracy, and parallel imaging using the sensitivity information is also possible.

さらに本実施形態によれば、感度マップを作成する際に、低分解能或いは情報が不足した参照画像であっても、高分解能の感度マップを得ることができる。従って、リファレンススキャンの撮像時間を短縮することも可能となる。   Furthermore, according to the present embodiment, when creating a sensitivity map, a high-resolution sensitivity map can be obtained even for a reference image with low resolution or lack of information. Therefore, it is possible to shorten the imaging time of the reference scan.

＜第四実施形態＞
第一〜第三実施形態では、複数チャンネルを有する受信コイルについてチャンネル毎に予め取得した感度マップとチャンネル毎の参照画像とから導出した座標変換式を用いて、当該受信コイルの感度マップを作成（補正）した実施形態であるが、本実施形態は、座標変換式を導出した複数チャンネルを有する受信コイル（第一の受信コイル）とは別の受信コイル（第二の受信コイル）の感度マップを、その座標変換を用いて作成する。 <Fourth embodiment>
In the first to third embodiments, a sensitivity map of the receiving coil is created using a coordinate conversion formula derived from a sensitivity map acquired in advance for each channel for a receiving coil having a plurality of channels and a reference image for each channel ( In this embodiment, a sensitivity map of a receiving coil (second receiving coil) different from the receiving coil (first receiving coil) having a plurality of channels from which a coordinate conversion equation is derived is used. It is created using the coordinate transformation.

第一の受信コイルと第二の受信コイルとの組み合わせは、第一の受信コイルが複数チャンネルを有する受信コイルであって且つ両受信コイルの感度分布が重なっていれば、特に限定されず、例えば、図１３に示すような、１組のちょう型コイル（２チャンネルコイル）６１とソレノイドコイル６２との組み合わせや、第一の受信コイルにさらにＣＲＣコイル６３を追加する場合などがある。また本実施形態では、感度マップ保存部３１に、第一の受信コイルの感度マップ、第一の受信コイルの各チャンネルの感度マップ、及び第二の受信コイルの感度マップが予め保存されているものとする。感度マップ保存部３１に予め保存されている感度マップの取得方法は特に限定されない。例えば、ファントム撮像、シミュレーション、或いは第二実施形態のようにシミュレーションで算出したものをファントム撮像結果で補正する、などの取得方法を採用できる。   The combination of the first receiving coil and the second receiving coil is not particularly limited as long as the first receiving coil is a receiving coil having a plurality of channels and the sensitivity distributions of both receiving coils overlap. As shown in FIG. 13, there is a combination of a pair of butterfly coils (two-channel coils) 61 and a solenoid coil 62, or a case where a CRC coil 63 is further added to the first receiving coil. In the present embodiment, the sensitivity map storage unit 31 stores in advance a sensitivity map of the first receiving coil, a sensitivity map of each channel of the first receiving coil, and a sensitivity map of the second receiving coil. And The method for obtaining the sensitivity map stored in advance in the sensitivity map storage unit 31 is not particularly limited. For example, it is possible to employ an acquisition method such as phantom imaging, simulation, or correction of a phantom imaging result obtained by simulation as in the second embodiment.

以下、図１３に示すような２チャンネル受信コイル６１とソレノイドコイル６２を組み合わせた受信コイル６０を例に、本実施形態の信号処理部の処理を説明する。図１４に処理の流れを示す。   Hereinafter, the processing of the signal processing unit of the present embodiment will be described by taking as an example a receiving coil 60 that combines a two-channel receiving coil 61 and a solenoid coil 62 as shown in FIG. FIG. 14 shows the flow of processing.

まず２チャンネル受信コイル６１を用いてリファレンススキャンを行い、各チャンネルＣｈ１、Ｃｈ２の参照画像を取得する（Ｓ７０１）。各チャンネルの参照画像の対応画素の画素値を要素とする２次元ベクトルＲを設定し、その偏角θを求める。一方、２チャンネル受信コイル６１について予め求めた各チャンネルＣｈ１、Ｃｈ２の感度分布について対応画素の画素値を要素とする２次元ベクトルＳを設定し、その偏角φを求める。全ての画素或いは予め設定した領域内の画素について、偏角θ、φをそれぞれ算出し、２つの偏角分布を得る（Ｓ７０２）。第一実施形態と同様に、感度マップから求めた偏角分布が参照画像から求めた偏角分布と同じような分布となるような座標変換式を決定する（Ｓ７０３）。   First, a reference scan is performed using the two-channel receiving coil 61, and reference images of the channels Ch1 and Ch2 are acquired (S701). A two-dimensional vector R having the pixel value of the corresponding pixel of the reference image of each channel as an element is set, and the deviation angle θ is obtained. On the other hand, a two-dimensional vector S having the pixel value of the corresponding pixel as an element is set for the sensitivity distribution of each channel Ch1 and Ch2 obtained in advance for the two-channel receiving coil 61, and the deviation angle φ is obtained. The declination angles θ and φ are calculated for all the pixels or the pixels in the preset area, and two declination distributions are obtained (S702). Similar to the first embodiment, a coordinate conversion formula is determined such that the declination distribution obtained from the sensitivity map becomes the same distribution as the declination distribution obtained from the reference image (S703).

この座標変換式を用いて、ソレノイドコイル６２について予め求めた感度分布を座標変換し、ソレノイドコイル６２の感度マップとする（Ｓ７０４）。これにより被検体に対する感度マップが得られるので、ソレノイドコイル６２を用いて撮像した画像の感度補正の精度を向上することができる。   Using this coordinate conversion formula, the sensitivity distribution obtained in advance for the solenoid coil 62 is coordinate-converted to obtain a sensitivity map of the solenoid coil 62 (S704). As a result, a sensitivity map for the subject is obtained, so that the accuracy of sensitivity correction of an image captured using the solenoid coil 62 can be improved.

なお以上の例ではソレノイドコイル６２の感度マップを作成する例であるが、座標変換式を用いて、２チャンネル受信コイル６１の各チャンネルＣｈ１、Ｃｈ２の感度マップを作成してもよいことは言うまでもない。   In the above example, the sensitivity map of the solenoid coil 62 is generated. However, it goes without saying that the sensitivity map of each channel Ch1 and Ch2 of the two-channel receiving coil 61 may be generated using a coordinate conversion formula. .

さらに以上の例では第一の受信コイル（２チャンネル受信コイル６１）の感度マップ（感度マップ保存部に格納すべき感度マップ）と参照画像をもとに作成した座標変換式を用いて第二の受信コイル（ソレノイドコイル）の感度マップ（被検体に対する感度マップ）を作成したが、第一の受信コイルと第二の受信コイルが一部重複していてもよい。例えば、上述した２チャンネル受信コイル６１とソレノイドコイル６２の組み合わせの例で、受信コイル６１のチャンネルとソレノイドコイルの感度マップ（感度マップ保存部に格納すべき感度マップ）の画素値を要素とするベクトルと参照画像の画素値を要素とするベクトルについて、それらの分布（或いは偏角の分布）が等しくなるような座標変換式を算出し、この座標変換式を用いてソレノイドコイルの感度マップ（感度マップ保存部に格納すべき感度マップ）を変換して、被検体に対する感度マップを作成してもよい。   Furthermore, in the above example, the second conversion is made using the coordinate conversion formula created based on the sensitivity map (sensitivity map to be stored in the sensitivity map storage unit) of the first reception coil (2-channel reception coil 61) and the reference image. Although the sensitivity map of the receiving coil (solenoid coil) (sensitivity map for the subject) is created, the first receiving coil and the second receiving coil may partially overlap. For example, in the example of the combination of the two-channel receiving coil 61 and the solenoid coil 62 described above, a vector having the pixel values of the sensitivity map (sensitivity map to be stored in the sensitivity map storage unit) of the channel of the receiving coil 61 and the solenoid coil as elements. And a vector whose elements are the pixel values of the reference image, a coordinate conversion formula is calculated so that their distributions (or distributions of declination) are equal. Using this coordinate conversion formula, a sensitivity map (sensitivity map) of the solenoid coil is calculated. The sensitivity map to be stored in the storage unit may be converted to create a sensitivity map for the subject.

本実施形態によれば、複数のコイルから構成される受信コイルにおいて、ソレノイドコイルのように感度補正の基準となるエレメントがあるときに、基準の感度むらを補正することができ効率よい感度補正を行うことができる。   According to the present embodiment, when there is an element serving as a reference for sensitivity correction, such as a solenoid coil, in a receiving coil composed of a plurality of coils, it is possible to correct the sensitivity variation of the reference and perform efficient sensitivity correction. It can be carried out.

以上、本発明の実施形態を説明したが、本発明は個々の実施形態に限定されるものではなく、例えば技術的に矛盾しない限り、各実施形態の構成や手順或いは変形例を他の実施形態に組み合わせたり、本発明において必須ではない構成を省くことも可能である。   As mentioned above, although embodiment of this invention was described, this invention is not limited to each embodiment, For example, as long as there is no technical contradiction, a structure, a procedure, or a modification of each embodiment is another embodiment. It is also possible to omit the configuration that is not essential in the present invention.

２：静磁場発生部、３：傾斜磁場発生部、４：シーケンサ、５：送信部、６：受信部、７：信号処理部、８：中央処理装置（ＣＰＵ）、１４ａ：送信コイル、１４ｂ：受信コイル、１８：外部記憶装置、１９：ディスプレイ、２０：操作部、３０：感度計算部、３１：感度マップ保存部（記憶部）、３２：感度マップ座標変換決定部、３２１：指標算出部、３２２：座標変換式算出部、３３：感度マップ作成部、３５：シミュレーション部、３６：初期感度マップ調整部、３７：絶対値／位相補正部、４０：画像作成部、４１：感度補正部、６０：受信コイル、６１：２チャンネル受信コイル、６２：ソレノイドコイル、６３：ＣＲＣコイル。

2: static magnetic field generation unit, 3: gradient magnetic field generation unit, 4: sequencer, 5: transmission unit, 6: reception unit, 7: signal processing unit, 8: central processing unit (CPU), 14a: transmission coil, 14b: Receiving coil, 18: external storage device, 19: display, 20: operation unit, 30: sensitivity calculation unit, 31: sensitivity map storage unit (storage unit), 32: sensitivity map coordinate conversion determination unit, 321: index calculation unit, 322: Coordinate conversion formula calculation unit, 33: Sensitivity map creation unit, 35: Simulation unit, 36: Initial sensitivity map adjustment unit, 37: Absolute value / phase correction unit, 40: Image creation unit, 41: Sensitivity correction unit, 60 : Reception coil, 61: 2-channel reception coil, 62: solenoid coil, 63: CRC coil.

Claims (13)

複数の受信チャンネルを有する受信コイルを備え、検査対象の核磁気共鳴信号を受信し、前記検査対象の画像を取得する撮像部と、
前記受信コイルの感度情報を用いて前記検査対象の画像作成に必要な計算を行う計算部と、
を備えた磁気共鳴イメージング装置であって、
前記検査対象以外で予め取得した前記受信コイルの感度マップ、及び、前記複数の受信チャンネルの感度マップまたは前記複数の受信チャンネルの感度マップから導出したマップを記憶する記憶部を更に備え、
前記計算部は、
前記複数の受信チャンネルの感度マップ間の関係を表す第一の指標マップと、前記撮像部が前記検査対象について前記受信コイルの受信チャンネルごとに取得した複数の参照画像間の関係を表す第二の指標マップとの差を最小化するように前記第一の指標マップを座標変換する座標変換式を算出する座標変換式算出部と、
前記座標変換式を用いて、前記検査対象存在時の受信感度マップを作成する感度マップ作成部と、を備えることを特徴とする磁気共鳴イメージング装置。 An imaging unit comprising a receiving coil having a plurality of receiving channels, receiving a nuclear magnetic resonance signal to be examined, and acquiring the image to be examined;
A calculation unit for performing calculation necessary for creating an image of the inspection object using sensitivity information of the receiving coil;
A magnetic resonance imaging apparatus comprising:
A storage unit that stores a sensitivity map of the reception coil acquired in advance other than the inspection target, and a sensitivity map of the plurality of reception channels or a map derived from the sensitivity map of the plurality of reception channels;
The calculator is
A first index map representing a relationship between sensitivity maps of the plurality of reception channels, and a second representing a relationship between a plurality of reference images acquired for each reception channel of the reception coil by the imaging unit for the inspection target. A coordinate conversion formula calculation unit that calculates a coordinate conversion formula for converting the coordinates of the first index map so as to minimize a difference from the index map;
A magnetic resonance imaging apparatus, comprising: a sensitivity map creating unit that creates a reception sensitivity map when the inspection target exists using the coordinate conversion formula. 請求項１に記載の磁気共鳴イメージング装置であって、
前記第一の指標マップは、各受信チャンネルの感度マップの画素値の比、前記画素値を要素とするベクトル、及び、前記ベクトルの偏角のいずれか一つを、画素値とするマップであり、
前記第二の指標マップは、各受信チャンネルの参照画像の画素値の比、前記画素値を要素とするベクトル、及び、前記ベクトルの偏角のいずれか一つを、画素値とするマップであることを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1,
The first index map is a map in which any one of a pixel value ratio of a sensitivity map of each reception channel, a vector having the pixel value as an element, and a declination of the vector is a pixel value. ,
The second index map is a map in which one of a pixel value ratio of a reference image of each reception channel, a vector having the pixel value as an element, and a declination of the vector is a pixel value. A magnetic resonance imaging apparatus. 請求項１に記載の磁気共鳴イメージング装置であって、
前記計算部は、前記検査対象の画像を、前記感度マップ作成部が作成した受信感度マップを用いて補正する感度補正部を備えることを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1,
The magnetic resonance imaging apparatus, wherein the calculation unit includes a sensitivity correction unit that corrects the examination target image using the reception sensitivity map created by the sensitivity map creation unit. 請求項１に記載の磁気共鳴イメージング装置であって、
前記感度マップ作成部は、前記複数の受信チャンネルを有する受信コイルについて複数チャンネルを組み合わせた全体としての前記受信感度マップを作成することを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1,
The magnetic resonance imaging apparatus, wherein the sensitivity map creation unit creates the reception sensitivity map as a whole by combining a plurality of channels with respect to the reception coil having the plurality of reception channels. 請求項１に記載の磁気共鳴イメージング装置であって、
前記複数の受信チャンネルを有する受信コイルは、広域用受信コイルであり、
前記計算部は、前記感度マップ作成部が作成した前記広域用受信コイルの前記受信感度マップと、前記広域用受信コイルとは異なる撮像用受信コイルについて予め求められた感度マップとを用いて、前記撮像用受信コイルの感度マップを補正する感度マップ補正部をさらに備えることを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1,
The receiving coil having the plurality of receiving channels is a wide area receiving coil,
The calculation unit uses the reception sensitivity map of the wide-area reception coil created by the sensitivity map creation unit, and a sensitivity map obtained in advance for an imaging reception coil different from the wide-area reception coil, A magnetic resonance imaging apparatus, further comprising a sensitivity map correction unit that corrects a sensitivity map of the imaging receiving coil. 請求項１に記載の磁気共鳴イメージング装置であって、
前記記憶部に記憶する第一の指標マップまたは前記複数の受信チャンネルの感度マップは、前記複数の受信チャンネルを有する受信コイルを用いてファントムを撮像することにより取得したファントム画像及びシミュレーションで得た感度画像の少なくとも一方を用いて作成された第一の指標マップまたは感度マップであることを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1,
The first index map stored in the storage unit or the sensitivity map of the plurality of reception channels is a phantom image acquired by imaging a phantom using a reception coil having the plurality of reception channels and a sensitivity obtained by simulation. A magnetic resonance imaging apparatus, wherein the magnetic resonance imaging apparatus is a first index map or sensitivity map created using at least one of images. 請求項６に記載の磁気共鳴イメージング装置であって、
前記記憶部に記憶する第一の指標マップまたは前記複数の受信チャンネルの感度マップは、前記シミュレーションで得た感度画像を前記ファントム画像で補正して作成されることを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 6,
The magnetic resonance imaging apparatus, wherein the first index map or the sensitivity map of the plurality of reception channels stored in the storage unit is created by correcting the sensitivity image obtained by the simulation with the phantom image. 請求項１に記載の磁気共鳴イメージング装置であって、
前記記憶部に記憶される感度マップは前記複数の受信チャンネルの感度マップであり、前記複数の受信チャンネルの感度マップ、及び、前記複数の参照画像は、画素値が複素数であり、
前記感度マップ作成部は、前記複数の受信チャンネルの感度マップ及び前記複数の参照画像の位相を揃える位相補正部をさらに備えることを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1,
The sensitivity map stored in the storage unit is a sensitivity map of the plurality of reception channels, the sensitivity map of the plurality of reception channels, and the plurality of reference images, pixel values are complex numbers,
The magnetic resonance imaging apparatus, wherein the sensitivity map creation unit further includes a phase correction unit that aligns the phases of the sensitivity maps of the plurality of reception channels and the plurality of reference images. 請求項８に記載の磁気共鳴イメージング装置であって、
前記計算部は、前記感度マップ作成部により、位相が補正された、各受信チャンネルの受信感度マップを用いて、パラレルイメージングに基く演算を行い、画像を作成することを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 8,
A magnetic resonance imaging apparatus characterized in that the calculation unit performs an operation based on parallel imaging using the reception sensitivity map of each reception channel whose phase is corrected by the sensitivity map generation unit, and generates an image. . 請求項１に記載の磁気共鳴イメージング装置であって、
前記感度マップ作成部は、前記複数の受信チャンネルを有する受信コイルの一部の受信チャンネルについて、又は、前記複数の受信チャンネルを有する受信コイルと感度領域が重なる別のコイルについて、前記受信感度マップを作成することを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1,
The sensitivity map creating unit is configured to obtain the reception sensitivity map for a part of reception channels of the reception coil having the plurality of reception channels, or for another coil having a sensitivity region overlapping with the reception coil having the plurality of reception channels. A magnetic resonance imaging apparatus comprising: 請求項１０に記載の磁気共鳴イメージング装置であって、
前記複数の受信チャンネルを有する受信コイルは、１つの受信チャンネルがソレノイドコイルで、他の受信チャンネルが前記ソレノイドコイルと重なる領域に感度を持つ複数チャンネルコイルであり、
前記座標変換式算出部は、少なくとも前記複数チャンネルコイルについて、前記第一の指標マップと前記第二の指標マップとを用いて座標変換式を作成し、
前記感度マップ作成部は、前記座標変換式を用いて、前記ソレノイドコイルの受信感度マップを作成することを特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 10,
The reception coil having the plurality of reception channels is a multiple channel coil having sensitivity in a region where one reception channel is a solenoid coil and another reception channel overlaps the solenoid coil,
The coordinate conversion formula calculation unit creates a coordinate conversion formula using at least the first index map and the second index map for at least the plurality of channel coils,
The magnetic resonance imaging apparatus, wherein the sensitivity map creating unit creates a reception sensitivity map of the solenoid coil using the coordinate conversion formula. 検査対象以外で取得した受信コイルの感度マップを用いて、検査対象に対する受信コイルの受信感度マップを算出する方法であって、
複数の受信チャンネルを有する受信コイルについて、検査対象以外で取得した前記複数の受信チャンネルの感度マップの画素間の関係を第一の指標として算出し、第一の指標マップを得、
検査対象に対して取得した前記複数の受信チャンネルの参照画像の画素間の関係を第二の指標として算出し、第二の指標マップを得、
前記第一の指標マップと第二の指標マップとが一致するように導出した座標変換式を検査対象以外で取得した受信コイルの感度マップに適用し、前記受信感度マップを算出することを特徴とする、受信コイルの感度マップ算出方法。 A method for calculating a reception sensitivity map of a reception coil with respect to an inspection target using a sensitivity map of the reception coil acquired outside the inspection target,
For a receiving coil having a plurality of receiving channels, the relationship between the pixels of the sensitivity map of the plurality of receiving channels acquired outside the inspection target is calculated as a first index, obtaining a first index map,
Calculating a relationship between pixels of the reference images of the plurality of reception channels acquired for the inspection target as a second index, obtaining a second index map;
Applying a coordinate transformation formula derived so that the first index map and the second index map match to a sensitivity map of a receiving coil acquired outside the inspection target, and calculating the reception sensitivity map, A receiving coil sensitivity map calculation method. 請求項１２に記載の受信コイルの受信感度マップを算出する方法であって、
前記座標変換式を、前記複数の受信チャンネルとは異なる受信コイルに適用し、前記複数の受信チャンネルとは異なる受信コイルについて、前記受信感度マップを作成することを特徴とする、受信コイルの感度マップ算出方法。

A method for calculating a reception sensitivity map of a reception coil according to claim 12, comprising:
The coordinate map is applied to a reception coil different from the plurality of reception channels, and the reception sensitivity map is created for a reception coil different from the plurality of reception channels. Calculation method.

JP2016239822A 2016-12-09 2016-12-09 Sensitivity map calculation method for magnetic resonance imaging device and receiving coil Active JP6783642B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016239822A JP6783642B2 (en) 2016-12-09 2016-12-09 Sensitivity map calculation method for magnetic resonance imaging device and receiving coil

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016239822A JP6783642B2 (en) 2016-12-09 2016-12-09 Sensitivity map calculation method for magnetic resonance imaging device and receiving coil

Publications (2)

Publication Number Publication Date
JP2018093989A true JP2018093989A (en) 2018-06-21
JP6783642B2 JP6783642B2 (en) 2020-11-11

Family

ID=62632494

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016239822A Active JP6783642B2 (en) 2016-12-09 2016-12-09 Sensitivity map calculation method for magnetic resonance imaging device and receiving coil

Country Status (1)

Country Link
JP (1) JP6783642B2 (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008006117A (en) * 2006-06-30 2008-01-17 Hitachi Medical Corp Magnetic resonance imaging apparatus
US20090224756A1 (en) * 2008-03-05 2009-09-10 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and magnetic resonance imaging method
US20120163692A1 (en) * 2009-09-17 2012-06-28 Koninklijke Philips Electronics N.V. Image intensity correction for magnetic resonance imaging
US20150241537A1 (en) * 2014-02-27 2015-08-27 Toshiba Medical Systems Corporation Magnetic Resonance Imaging with Consistent Geometries
US20160334492A1 (en) * 2015-05-13 2016-11-17 Medizinische Universität Wien Method for Magnetic Resonance Imaging

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008006117A (en) * 2006-06-30 2008-01-17 Hitachi Medical Corp Magnetic resonance imaging apparatus
US20090224756A1 (en) * 2008-03-05 2009-09-10 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP2009207756A (en) * 2008-03-05 2009-09-17 Toshiba Corp Magnetic resonance imaging apparatus
US20120163692A1 (en) * 2009-09-17 2012-06-28 Koninklijke Philips Electronics N.V. Image intensity correction for magnetic resonance imaging
JP2013505045A (en) * 2009-09-17 2013-02-14 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Image intensity correction for magnetic resonance imaging
US20150241537A1 (en) * 2014-02-27 2015-08-27 Toshiba Medical Systems Corporation Magnetic Resonance Imaging with Consistent Geometries
JP2015160141A (en) * 2014-02-27 2015-09-07 株式会社東芝 Magnetic resonance imaging system
US20160334492A1 (en) * 2015-05-13 2016-11-17 Medizinische Universität Wien Method for Magnetic Resonance Imaging

Also Published As

Publication number Publication date
JP6783642B2 (en) 2020-11-11

Similar Documents

Publication Publication Date Title
NL2003804C2 (en) System and method for automated scan planning using symmetry detection and image registration.
JP6464088B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP2016504144A (en) Metal resistant MR imaging
US20090259120A1 (en) Magnetic resonance imaging apparatus and image generating method
US20230260088A1 (en) Medical information processing apparatus and medical information processing method
JP5683984B2 (en) Magnetic resonance imaging apparatus and nonlinear distortion correction method
US20160154079A1 (en) Method and magnetic resonance apparatus for acquiring a sensitivity map for at least one local coil in a magnetic resonance scanner
JP5154751B2 (en) Medical image processing device
US9470772B2 (en) Magnetic resonance image distortion correction with a reference volume
JP2006255046A (en) Magnetic resonance imaging method and image processing apparatus
JP5181235B2 (en) Image processing apparatus, medical image photographing apparatus, program, and image processing method
US20100130849A1 (en) System and method for patient motion compensation during magnetic resonance scans
US11650281B2 (en) Excitation region setting method and magnetic resonance imaging apparatus
JP7419145B2 (en) Medical information processing device and medical information processing method
JP6783642B2 (en) Sensitivity map calculation method for magnetic resonance imaging device and receiving coil
JP4901420B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP6113012B2 (en) Magnetic resonance imaging apparatus and method for calculating B1 map for correction
US11698432B2 (en) Magnetic resonance imaging system, and main magnetic field correction method therefor and storage medium
JP5484001B2 (en) Magnetic resonance imaging apparatus and image correction method
JP4901627B2 (en) Magnetic resonance imaging device
US20240074671A1 (en) Image processing apparatus, image processing method, and magnetic resonance imaging apparatus
US11041921B2 (en) MRI apparatus, image processing apparatus, and MRI method
JP2018186892A (en) Magnetic resonance imaging apparatus
JP5638324B2 (en) Magnetic resonance imaging apparatus and image correction method
JP2018186942A (en) Magnetic resonance imaging apparatus and pulse designing method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190701

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20200520

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200616

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200728

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20201006

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20201022

R150 Certificate of patent or registration of utility model

Ref document number: 6783642

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250