JP7451366B2 - 画像生成装置、画像生成方法、および画像生成プログラム - Google Patents
画像生成装置、画像生成方法、および画像生成プログラム Download PDFInfo
- Publication number
- JP7451366B2 JP7451366B2 JP2020157877A JP2020157877A JP7451366B2 JP 7451366 B2 JP7451366 B2 JP 7451366B2 JP 2020157877 A JP2020157877 A JP 2020157877A JP 2020157877 A JP2020157877 A JP 2020157877A JP 7451366 B2 JP7451366 B2 JP 7451366B2
- Authority
- JP
- Japan
- Prior art keywords
- sensitivity distribution
- magnetic resonance
- coil
- image
- sensitivity
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 35
- 230000035945 sensitivity Effects 0.000 claims description 311
- 238000009826 distribution Methods 0.000 claims description 225
- 238000003384 imaging method Methods 0.000 claims description 54
- 239000002131 composite material Substances 0.000 claims description 28
- 230000005484 gravity Effects 0.000 claims description 23
- 238000013519 translation Methods 0.000 claims description 5
- 230000006870 function Effects 0.000 description 214
- 238000012545 processing Methods 0.000 description 79
- 239000013598 vector Substances 0.000 description 36
- 238000002595 magnetic resonance imaging Methods 0.000 description 25
- 230000008569 process Effects 0.000 description 18
- 230000015572 biosynthetic process Effects 0.000 description 13
- 238000003786 synthesis reaction Methods 0.000 description 13
- 238000004891 communication Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 11
- 230000003068 static effect Effects 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 5
- 239000004065 semiconductor Substances 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000002603 single-photon emission computed tomography Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000001151 other effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/36—Electrical details, e.g. matching or coupling of the coil to the receiver
- G01R33/3607—RF waveform generators, e.g. frequency generators, amplitude-, frequency- or phase modulators or shifters, pulse programmers, digital to analog converters for the RF signal, means for filtering or attenuating of the RF signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
- G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
Description
本明細書及び図面に開示の実施形態は、画像生成装置、画像生成方法、および画像生成プログラムに関する。
従来、ｋ空間における非カーテシアンでの撮像により収集された磁気共鳴データに基づく磁気共鳴画像の生成において、当該磁気共鳴データの収集点を当該ｋ空間の格子点に配置に関するグリッディングが、実行される。グリッディングにおける感度は、通常、再構成領域の中央部分の信号対雑音比(ｓｉｇｎａｌ－ｔｏ－ｎｏｉｓｅ ｒａｔｉｏ：以下、Ｓ／Ｎと呼ぶ）が最大となる。このため、グリッディングを用いて再構成された画像
周囲の領域は、Ｓ／Ｎが低下し、シェーディングやストリークなどのアーチファクトが発生することがある。
周囲の領域は、Ｓ／Ｎが低下し、シェーディングやストリークなどのアーチファクトが発生することがある。
また、磁気共鳴画像の生成において、被検体に対する本スキャン前に収集された受信コイルの感度マップを用いて、当該受信コイルごとの感度補正に関する感度エンコーディング（以下、ＳＥＮＳＥ（ｓｅｎｓｉｔｉｖｉｔｙ ｅｎｃｏｄｉｎｇ）と呼ぶ）を行うことがある。
ｋ空間における非カーテシアンでの撮像により収集された磁気共鳴データに基づく磁気共鳴画像の生成において、グリッディングとＳＥＮＳＥとを組み合わせた場合、アーチファクトの発生やＳ／Ｎの低下が、再構成された磁気共鳴画像に現れるとがある。
本明細書及び図面に開示の実施形態が解決しようとする課題の一つは、アーチファクトの発生を抑制し、Ｓ／Ｎを向上させた画像を生成することにある。ただし、本明細書及び図面に開示の実施形態により解決しようとする課題は上記課題に限られない。後述する実施形態に示す各構成による各効果に対応する課題を他の課題として位置づけることもできる。
実施形態に係る画像生成装置は、取得部と、位置合わせ部と、生成部とを備える、取得部は、被検体の撮像に用いられる受信コイルの感度の分布を示すコイル感度分布と、ｋ空間における非カーテシアンでの前記撮像により収集された磁気共鳴データとを取得する。位置合わせ部は、前記コイル感度分布に基づいて、前記ｋ空間における前記磁気共鳴データの配置に関するグリッディングの感度の分布を示すグリッディング感度分布と前記コイル感度分布との位置合わせを実行する。生成部は、前記位置合わせの結果と前記磁気共鳴データと前記コイル感度分布と前記グリッディング感度分布とに基づいて、磁気共鳴画像を生成する。
以下、図面を参照しながら、画像生成装置、画像生成方法、および画像生成プログラムの実施形態について詳細に説明する。図１は、画像生成装置１の一例を示すブロック図である。画像生成装置１は、例えば、本画像生成装置１における各種機能が搭載された各種モダリティや、院内などにおけるサーバに搭載される。なお、画像生成装置１における各種機能は、医用画像管理システム（以下、ＰＡＣＳ（Ｐｉｃｔｕｒｅ Ａｒｃｈｉｖｉｎｇ ａｎｄ Ｃｏｍｍｕｎｉｃａｔｉｏｎ Ｓｙｓｔｅｍｓ）のサーバや、病院情報システム（以下、ＨＩＳ（Ｈｏｓｐｉｔａｌ Ｉｎｆｏｒｍａｔｉｏｎ Ｓｙｓｔｅｍ）と呼ぶ）のサーバなどに搭載されてもよい。
また、本画像生成装置１における各種機能が搭載されたとは、例えば、磁気共鳴イメージング（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ：以下、ＭＲＩと呼ぶ）装置、ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ：陽電子放出コンピュータ断層撮影）－ＭＲＩ装置、ＳＰＥＣＴ（ｓｉｎｇｌｅ ｐｈｏｔｏｎ
ｅｍｉｓｓｉｏｎ ｃｏｍｐｕｔｅｄ ｔｏｍｏｇｒａｐｈｙ：単一光子放出コンピュータ断層撮影）－ＭＲＩ装置などである。以下、説明を具体的にするために、画像生成装置１は、ＭＲＩ装置に搭載されているものとする。このとき、ＭＲＩ装置は、処理回路１５における各種機能を有することとなる。
ｅｍｉｓｓｉｏｎ ｃｏｍｐｕｔｅｄ ｔｏｍｏｇｒａｐｈｙ：単一光子放出コンピュータ断層撮影）－ＭＲＩ装置などである。以下、説明を具体的にするために、画像生成装置１は、ＭＲＩ装置に搭載されているものとする。このとき、ＭＲＩ装置は、処理回路１５における各種機能を有することとなる。
（第１の実施形態）
図２は、本実施形態に係るＭＲＩ装置１００の一例を示す図である。図２に示すように、ＭＲＩ装置１００における画像生成装置１は、入出力インターフェース１７をさらに有する。なお、画像生成装置１は、図１に示すように、入出力インターフェース１７を非搭載であってもよい。図２に示すように、ＭＲＩ装置１００は、静磁場磁石１０１と、傾斜磁場コイル１０３と、傾斜磁場電源１０５と、寝台１０７と、寝台制御回路（寝台制御部）１０９と、送信回路１１３と、送信コイル１１５と、受信コイル１１７と、受信回路１１９と、撮像制御回路（収集部）１２１と、システム制御回路（システム制御部）１２３と、記憶装置１２５と、画像生成装置１とを備える。
図２は、本実施形態に係るＭＲＩ装置１００の一例を示す図である。図２に示すように、ＭＲＩ装置１００における画像生成装置１は、入出力インターフェース１７をさらに有する。なお、画像生成装置１は、図１に示すように、入出力インターフェース１７を非搭載であってもよい。図２に示すように、ＭＲＩ装置１００は、静磁場磁石１０１と、傾斜磁場コイル１０３と、傾斜磁場電源１０５と、寝台１０７と、寝台制御回路（寝台制御部）１０９と、送信回路１１３と、送信コイル１１５と、受信コイル１１７と、受信回路１１９と、撮像制御回路（収集部）１２１と、システム制御回路（システム制御部）１２３と、記憶装置１２５と、画像生成装置１とを備える。
静磁場磁石１０１は、中空の略円筒状に形成された磁石である。静磁場磁石１０１は、
内部の空間に略一様な静磁場を発生する。静磁場磁石１０１としては、例えば、超伝導磁
石等が使用される。
内部の空間に略一様な静磁場を発生する。静磁場磁石１０１としては、例えば、超伝導磁
石等が使用される。
傾斜磁場コイル１０３は、中空の略円筒形状に形成されたコイルであり、円筒形の冷却
容器の内面側に配置される。傾斜磁場コイル１０３は、傾斜磁場電源１０５から個別に電
流供給を受けて、互いに直交するＸ、Ｙ、及びＺの各軸に沿って磁場強度が変化する傾斜
磁場を発生する。傾斜磁場コイル１０３によって発生されるＸ、Ｙ、Ｚ各軸の傾斜磁場は
、例えば、スライス選択用傾斜磁場、位相エンコード用傾斜磁場および周波数エンコード
用傾斜磁場（リードアウト傾斜磁場ともいう）を形成する。スライス選択用傾斜磁場は、
任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場は、空間的位置に
応じて磁気共鳴信号（以下、ＭＲ（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ）信号と呼ぶ
）の位相を変化させるために利用される。周波数エンコード用傾斜磁場は、空間的位置に
応じてＭＲ信号の周波数を変化させるために利用される。
容器の内面側に配置される。傾斜磁場コイル１０３は、傾斜磁場電源１０５から個別に電
流供給を受けて、互いに直交するＸ、Ｙ、及びＺの各軸に沿って磁場強度が変化する傾斜
磁場を発生する。傾斜磁場コイル１０３によって発生されるＸ、Ｙ、Ｚ各軸の傾斜磁場は
、例えば、スライス選択用傾斜磁場、位相エンコード用傾斜磁場および周波数エンコード
用傾斜磁場（リードアウト傾斜磁場ともいう）を形成する。スライス選択用傾斜磁場は、
任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場は、空間的位置に
応じて磁気共鳴信号（以下、ＭＲ（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ）信号と呼ぶ
）の位相を変化させるために利用される。周波数エンコード用傾斜磁場は、空間的位置に
応じてＭＲ信号の周波数を変化させるために利用される。
傾斜磁場電源１０５は、撮像制御回路１２１の制御により、傾斜磁場コイル１０３に電
流を供給する電源装置である。
流を供給する電源装置である。
寝台１０７は、被検体Ｐが載置される天板１０７１を備えた装置である。寝台１０７は
、寝台制御回路１０９による制御のもと、被検体Ｐが載置された天板１０７１を、ボア１
１１内へ挿入する。
、寝台制御回路１０９による制御のもと、被検体Ｐが載置された天板１０７１を、ボア１
１１内へ挿入する。
寝台制御回路１０９は、寝台１０７を制御する回路である。寝台制御回路１０９は、入
出力インターフェース１７を介した操作者の指示により寝台１０７を駆動することで、天
板１０７１を長手方向および上下方向、場合によっては左右方向へ移動させる。
出力インターフェース１７を介した操作者の指示により寝台１０７を駆動することで、天
板１０７１を長手方向および上下方向、場合によっては左右方向へ移動させる。
送信回路１１３は、撮像制御回路１２１の制御により、ラーモア周波数で変調された高
周波パルスを送信コイル１１５に供給する。例えば、送信回路１１３は、発振部や位相選
択部、周波数変換部、振幅変調部、ＲＦアンプなどを有する。発振部は、静磁場中におけ
る対象原子核に固有の共鳴周波数のＲＦパルスを発生する。位相選択部は、発振部によっ
て発生したＲＦパルスの位相を選択する。周波数変換部は、位相選択部から出力されたＲ
Ｆパルスの周波数を変換する。振幅変調部は、周波数変換部から出力されたＲＦパルスの
振幅を例えばｓｉｎｃ関数に従って変調する。ＲＦアンプは、振幅変調部から出力された
ＲＦパルスを増幅して送信コイル１１５に供給する。
周波パルスを送信コイル１１５に供給する。例えば、送信回路１１３は、発振部や位相選
択部、周波数変換部、振幅変調部、ＲＦアンプなどを有する。発振部は、静磁場中におけ
る対象原子核に固有の共鳴周波数のＲＦパルスを発生する。位相選択部は、発振部によっ
て発生したＲＦパルスの位相を選択する。周波数変換部は、位相選択部から出力されたＲ
Ｆパルスの周波数を変換する。振幅変調部は、周波数変換部から出力されたＲＦパルスの
振幅を例えばｓｉｎｃ関数に従って変調する。ＲＦアンプは、振幅変調部から出力された
ＲＦパルスを増幅して送信コイル１１５に供給する。
送信コイル１１５は、傾斜磁場コイル１０３の内側に配置されたＲＦ（Ｒａｄｉｏ Ｆ
ｒｅｑｕｅｎｃｙ）コイルである。送信コイル１１５は、送信回路１１３からの出力に応
じて、高周波磁場に相当するＲＦパルスを発生する。
ｒｅｑｕｅｎｃｙ）コイルである。送信コイル１１５は、送信回路１１３からの出力に応
じて、高周波磁場に相当するＲＦパルスを発生する。
受信コイル１１７は、傾斜磁場コイル１０３の内側に配置されたＲＦコイルである。受
信コイル１１７は、高周波磁場によって被検体Ｐから放射されるＭＲ信号を受信する。受
信コイル１１７は、受信されたＭＲ信号を受信回路１１９へ出力する。受信コイル１１７
は、例えば、１以上、典型的には複数のコイルエレメントを有するコイルアレイである。
以下、説明を具体的にするために、受信コイル１１７は、複数のコイルエレメントを有す
るコイルアレイとして説明する。
信コイル１１７は、高周波磁場によって被検体Ｐから放射されるＭＲ信号を受信する。受
信コイル１１７は、受信されたＭＲ信号を受信回路１１９へ出力する。受信コイル１１７
は、例えば、１以上、典型的には複数のコイルエレメントを有するコイルアレイである。
以下、説明を具体的にするために、受信コイル１１７は、複数のコイルエレメントを有す
るコイルアレイとして説明する。
なお、受信コイル１１７は、一つのコイルエレメントにより構成されてもよい。また、
図２において送信コイル１１５と受信コイル１１７とは別個のＲＦコイルとして記載され
ているが、送信コイル１１５と受信コイル１１７とは、一体化された送受信コイルとして
実施されてもよい。送受信コイルは、被検体Ｐの撮像部位に対応し、例えば、頭部コイル
のような局所的な送受信ＲＦコイルである。
図２において送信コイル１１５と受信コイル１１７とは別個のＲＦコイルとして記載され
ているが、送信コイル１１５と受信コイル１１７とは、一体化された送受信コイルとして
実施されてもよい。送受信コイルは、被検体Ｐの撮像部位に対応し、例えば、頭部コイル
のような局所的な送受信ＲＦコイルである。
受信回路１１９は、撮像制御回路１２１の制御により、受信コイル１１７から出力され
たＭＲ信号に基づいて、デジタルのＭＲ信号（以下、ＭＲデータと呼ぶ）を生成する。具
体的には、受信回路１１９は、受信コイル１１７から出力されたＭＲ信号に対して各種信
号処理を施した後、各種信号処理が施されたデータに対してアナログ／デジタル（Ａ／Ｄ
（Ａｎａｌｏｇ ｔｏ Ｄｉｇｉｔａｌ））変換して、ＭＲデータを生成する。受信回路
１１９は、生成されたＭＲデータを、撮像制御回路１２１に出力する。例えば、ＭＲデー
タは、コイルエレメントごとに生成され、コイルエレメントを識別するタグとともに、撮
像制御回路１２１に出力される。
たＭＲ信号に基づいて、デジタルのＭＲ信号（以下、ＭＲデータと呼ぶ）を生成する。具
体的には、受信回路１１９は、受信コイル１１７から出力されたＭＲ信号に対して各種信
号処理を施した後、各種信号処理が施されたデータに対してアナログ／デジタル（Ａ／Ｄ
（Ａｎａｌｏｇ ｔｏ Ｄｉｇｉｔａｌ））変換して、ＭＲデータを生成する。受信回路
１１９は、生成されたＭＲデータを、撮像制御回路１２１に出力する。例えば、ＭＲデー
タは、コイルエレメントごとに生成され、コイルエレメントを識別するタグとともに、撮
像制御回路１２１に出力される。
撮像制御回路１２１は、処理回路１５から出力された撮像プロトコルに従って、傾斜磁
場電源１０５、送信回路１１３及び受信回路１１９等を制御し、被検体Ｐに対する撮像を
行う。撮像プロトコルは、検査の種類に応じたパルスシーケンスを有する。撮像プロトコ
ルには、傾斜磁場電源１０５により傾斜磁場コイル１０３に供給される電流の大きさ、傾
斜磁場電源１０５により電流が傾斜磁場コイル１０３に供給されるタイミング、送信回路
１１３により送信コイル１１５に供給される高周波パルスの大きさや時間幅、送信回路１
１３により送信コイル１１５に高周波パルスが供給されるタイミング、受信コイル１１７
によりＭＲ信号が受信されるタイミング等が定義されている。撮像制御回路１２１は、傾
斜磁場電源１０５、送信回路１１３及び受信回路１１９等を駆動して被検体Ｐを撮像した
結果、受信回路１１９からＭＲデータを受信すると、受信したＭＲデータを画像生成装置
１等へ転送する。
場電源１０５、送信回路１１３及び受信回路１１９等を制御し、被検体Ｐに対する撮像を
行う。撮像プロトコルは、検査の種類に応じたパルスシーケンスを有する。撮像プロトコ
ルには、傾斜磁場電源１０５により傾斜磁場コイル１０３に供給される電流の大きさ、傾
斜磁場電源１０５により電流が傾斜磁場コイル１０３に供給されるタイミング、送信回路
１１３により送信コイル１１５に供給される高周波パルスの大きさや時間幅、送信回路１
１３により送信コイル１１５に高周波パルスが供給されるタイミング、受信コイル１１７
によりＭＲ信号が受信されるタイミング等が定義されている。撮像制御回路１２１は、傾
斜磁場電源１０５、送信回路１１３及び受信回路１１９等を駆動して被検体Ｐを撮像した
結果、受信回路１１９からＭＲデータを受信すると、受信したＭＲデータを画像生成装置
１等へ転送する。
撮像制御回路１２１は、例えば、ｋ空間における非カーテシアン（Ｎｏｎ－Ｃａｒｔｅ
ｓｉａｎ）での撮像に関するパルスシーケンスを実行することにより、ＭＲデータ（以下
、非カーテシアンデータと呼ぶ）を収集する。非カーテシアンでの撮像とは、例えば、ラ
ジアル収集などである。以下、説明を具体的にするために、非カーテシアンでの撮像は、
ラジアル収集であるものとして説明する。また、撮像制御回路１２１は、被検体Ｐの撮像
に用いられる受信コイル１１７の感度を示す分布（以下、コイル感度分布とよぶ）の生成
に関するＭＲデータ（以下、コイル感度データと呼ぶ）を収集する。コイル感度分布は、
コイル感度マップとも称され、複素数のデータで表現される。コイル感度データの収集は
、例えば、当該非カーテシアンデータの収集に先立って、プリスキャンやロケータスキャ
ンなどにおいて撮像制御回路１２１により実行される。このとき、コイル感度分布は、例
えば、処理回路１５における生成機能１５７により、コイル感度データに対するフーリエ
変換により生成される。撮像制御回路１２１は、例えばプロセッサにより実現される。
ｓｉａｎ）での撮像に関するパルスシーケンスを実行することにより、ＭＲデータ（以下
、非カーテシアンデータと呼ぶ）を収集する。非カーテシアンでの撮像とは、例えば、ラ
ジアル収集などである。以下、説明を具体的にするために、非カーテシアンでの撮像は、
ラジアル収集であるものとして説明する。また、撮像制御回路１２１は、被検体Ｐの撮像
に用いられる受信コイル１１７の感度を示す分布（以下、コイル感度分布とよぶ）の生成
に関するＭＲデータ（以下、コイル感度データと呼ぶ）を収集する。コイル感度分布は、
コイル感度マップとも称され、複素数のデータで表現される。コイル感度データの収集は
、例えば、当該非カーテシアンデータの収集に先立って、プリスキャンやロケータスキャ
ンなどにおいて撮像制御回路１２１により実行される。このとき、コイル感度分布は、例
えば、処理回路１５における生成機能１５７により、コイル感度データに対するフーリエ
変換により生成される。撮像制御回路１２１は、例えばプロセッサにより実現される。
上記説明では、「プロセッサ」が各機能に対応するプログラムをメモリ１３から読み出
して実行する例を説明したが、実施形態はこれに限定されない。「プロセッサ」という文
言は、例えば、ＣＰＵ、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）
、特定用途向け集積回路（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａ
ｔｅｄ Ｃｉｒｃｕｉｔ：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プロ
グラマブル論理デバイス（Ｓｉｍｐｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅ
ｖｉｃｅ：ＳＰＬＤ）、複合プログラマブル論理デバイス（Ｃｏｍｐｌｅｘ Ｐｒｏｇｒ
ａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＣＰＬＤ）、及びフィールドプログラマブ
ルゲートアレイ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ：ＦＰ
ＧＡ））等の回路を意味する。
して実行する例を説明したが、実施形態はこれに限定されない。「プロセッサ」という文
言は、例えば、ＣＰＵ、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）
、特定用途向け集積回路（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａ
ｔｅｄ Ｃｉｒｃｕｉｔ：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プロ
グラマブル論理デバイス（Ｓｉｍｐｌｅ Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅ
ｖｉｃｅ：ＳＰＬＤ）、複合プログラマブル論理デバイス（Ｃｏｍｐｌｅｘ Ｐｒｏｇｒ
ａｍｍａｂｌｅ Ｌｏｇｉｃ Ｄｅｖｉｃｅ：ＣＰＬＤ）、及びフィールドプログラマブ
ルゲートアレイ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ：ＦＰ
ＧＡ））等の回路を意味する。
プロセッサが例えばＣＰＵである場合、プロセッサはメモリ１３に保存されたプログラ
ムを読み出して実行することで機能を実現する。一方、プロセッサがＡＳＩＣである場合
、メモリ１３にプログラムを保存する代わりに、当該機能がプロセッサの回路内に論理回
路として直接組み込まれる。なお、本実施形態の各プロセッサは、プロセッサごとに単一
の回路として構成される場合に限らず、複数の独立した回路を組み合わせて１つのプロセ
ッサとして構成し、その機能を実現するようにしてもよい。また、単一の記憶回路が各処
理機能に対応するプログラムを記憶するものとして説明したが、複数の記憶回路を分散し
て配置して、処理回路は個別の記憶回路から対応するプログラムを読み出す構成としても
構わない。
ムを読み出して実行することで機能を実現する。一方、プロセッサがＡＳＩＣである場合
、メモリ１３にプログラムを保存する代わりに、当該機能がプロセッサの回路内に論理回
路として直接組み込まれる。なお、本実施形態の各プロセッサは、プロセッサごとに単一
の回路として構成される場合に限らず、複数の独立した回路を組み合わせて１つのプロセ
ッサとして構成し、その機能を実現するようにしてもよい。また、単一の記憶回路が各処
理機能に対応するプログラムを記憶するものとして説明したが、複数の記憶回路を分散し
て配置して、処理回路は個別の記憶回路から対応するプログラムを読み出す構成としても
構わない。
システム制御回路１２３は、ハードウェア資源として図示していないプロセッサ、ＲＯ
Ｍ（Ｒｅａｄ－Ｏｎｌｙ Ｍｅｍｏｒｙ）やＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅ
ｍｏｒｙ）等のメモリ等を有し、システム制御機能によりＭＲＩ装置１００を制御する。
具体的には、システム制御回路１２３は、記憶装置１２５に記憶されたシステム制御プロ
グラムを読み出してメモリ上に展開し、展開されたシステム制御プログラムに従って本Ｍ
ＲＩ装置１００の各回路を制御する。例えば、システム制御回路１２３は、入出力インタ
ーフェース１７を介して操作者から入力された撮像条件に基づいて、撮像プロトコルを記
憶装置１２５から読み出す。システム制御回路１２３は、撮像プロトコルを撮像制御回路
１２１に送信し、被検体Ｐに対する撮像を制御する。システム制御回路１２３は、例えば
プロセッサにより実現される。なお、システム制御回路１２３は、処理回路１５に組み込
まれてもよい。このとき、システム制御機能は処理回路１５により実行され、処理回路１
５は、システム制御回路１２３の代替として機能する。
Ｍ（Ｒｅａｄ－Ｏｎｌｙ Ｍｅｍｏｒｙ）やＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅ
ｍｏｒｙ）等のメモリ等を有し、システム制御機能によりＭＲＩ装置１００を制御する。
具体的には、システム制御回路１２３は、記憶装置１２５に記憶されたシステム制御プロ
グラムを読み出してメモリ上に展開し、展開されたシステム制御プログラムに従って本Ｍ
ＲＩ装置１００の各回路を制御する。例えば、システム制御回路１２３は、入出力インタ
ーフェース１７を介して操作者から入力された撮像条件に基づいて、撮像プロトコルを記
憶装置１２５から読み出す。システム制御回路１２３は、撮像プロトコルを撮像制御回路
１２１に送信し、被検体Ｐに対する撮像を制御する。システム制御回路１２３は、例えば
プロセッサにより実現される。なお、システム制御回路１２３は、処理回路１５に組み込
まれてもよい。このとき、システム制御機能は処理回路１５により実行され、処理回路１
５は、システム制御回路１２３の代替として機能する。
記憶装置１２５は、システム制御回路１２３において実行される各種プログラム、各種
撮像プロトコル、撮像プロトコルを規定する複数の撮像パラメータを含む撮像条件等を記
憶する。記憶装置１２５は、例えば、ＲＡＭ、フラッシュメモリ等の半導体メモリ素子、
ＨＤＤ（Ｈａｒｄ ｄｉｓｋ Ｄｒｉｖｅ）、ＳＳＤ（Ｓｏｌｉｄ Ｓｔａｔｅ Ｄｒｉ
ｖｅ）、光ディスク等である。また、記憶装置１２５は、ＣＤ（Ｃｏｍｐａｃｔ Ｄｉｓ
ｃ）－ＲＯＭドライブやＤＶＤ（Ｄｉｇｉｔａｌ Ｖｅｒｓａｔｉｌｅ Ｄｉｓｃ）ドラ
イブ、フラッシュメモリ等の可搬型記憶媒体との間で種々の情報を読み書きする駆動装置
等であってもよい。なお、記憶装置１２５に記憶されるデータは、メモリ１３に記憶され
てもよい。このとき、メモリ１３は、記憶装置１２５の代替として機能する。
撮像プロトコル、撮像プロトコルを規定する複数の撮像パラメータを含む撮像条件等を記
憶する。記憶装置１２５は、例えば、ＲＡＭ、フラッシュメモリ等の半導体メモリ素子、
ＨＤＤ（Ｈａｒｄ ｄｉｓｋ Ｄｒｉｖｅ）、ＳＳＤ（Ｓｏｌｉｄ Ｓｔａｔｅ Ｄｒｉ
ｖｅ）、光ディスク等である。また、記憶装置１２５は、ＣＤ（Ｃｏｍｐａｃｔ Ｄｉｓ
ｃ）－ＲＯＭドライブやＤＶＤ（Ｄｉｇｉｔａｌ Ｖｅｒｓａｔｉｌｅ Ｄｉｓｃ）ドラ
イブ、フラッシュメモリ等の可搬型記憶媒体との間で種々の情報を読み書きする駆動装置
等であってもよい。なお、記憶装置１２５に記憶されるデータは、メモリ１３に記憶され
てもよい。このとき、メモリ１３は、記憶装置１２５の代替として機能する。
画像生成装置１は、通信インターフェース１１と、メモリ１３と、処理回路１５とを有
する。図１および図２に示すように、画像生成装置１において、通信インターフェース１
１と、メモリ１３と、処理回路１５とはバスにより電気的に接続されている。図１および
図２に示すように、画像生成装置１は、通信インターフェース１１を介して、ネットワー
クに接続されている。ネットワークには、例えば、各種モダリティや、ＨＩＳ、放射線情
報システム（ＲＩＳ：Ｒａｄｉｏｌｏｇｙ Ｉｎｆｏｒｍａｔｉｏｎ Ｓｙｓｔｅｍ）等
の医療機関内の情報処理システムと互いに通信可能に接続される。なお、図１に示す画像
生成装置１には、ユーザの各種時を入力するための入力インターフェースや、生成機能１
５７により生成された医用画像を表示するディスプレイ（出力インターフェース）を、図
２に示すように入出力インターフェース１７を有していてもよい。
する。図１および図２に示すように、画像生成装置１において、通信インターフェース１
１と、メモリ１３と、処理回路１５とはバスにより電気的に接続されている。図１および
図２に示すように、画像生成装置１は、通信インターフェース１１を介して、ネットワー
クに接続されている。ネットワークには、例えば、各種モダリティや、ＨＩＳ、放射線情
報システム（ＲＩＳ：Ｒａｄｉｏｌｏｇｙ Ｉｎｆｏｒｍａｔｉｏｎ Ｓｙｓｔｅｍ）等
の医療機関内の情報処理システムと互いに通信可能に接続される。なお、図１に示す画像
生成装置１には、ユーザの各種時を入力するための入力インターフェースや、生成機能１
５７により生成された医用画像を表示するディスプレイ（出力インターフェース）を、図
２に示すように入出力インターフェース１７を有していてもよい。
通信インターフェース１１は、例えば、被検体Ｐに対する検査において当該被検体Ｐを
撮像する各種モダリティや、ＨＩＳ、ＰＡＣＳなどとの間でデータ通信を行う。通信イン
ターフェース１１と各種モダリティおよび病院情報システムとの通信の規格は、如何なる
規格であっても良いが、例えば、ＨＬ７（Ｈｅａｒｔｈ Ｌｅｖｅｌ ７）、ＤＩＣＯＭ
、又はその両方等が挙げられる。
撮像する各種モダリティや、ＨＩＳ、ＰＡＣＳなどとの間でデータ通信を行う。通信イン
ターフェース１１と各種モダリティおよび病院情報システムとの通信の規格は、如何なる
規格であっても良いが、例えば、ＨＬ７（Ｈｅａｒｔｈ Ｌｅｖｅｌ ７）、ＤＩＣＯＭ
、又はその両方等が挙げられる。
メモリ１３は、種々の情報を記憶する記憶回路により実現される。例えば、メモリ１３
は、ＨＤＤやＳＳＤ、集積回路記憶装置等の記憶装置である。メモリ１３は、記憶部に相
当する。なお、メモリ１３は、ＨＤＤやＳＳＤ等以外にも、ＲＡＭ（Ｒａｎｄｏｍ Ａｃ
ｃｅｓｓ Ｍｅｍｏｒｙ）、フラッシュメモリ等の半導体メモリ素子、ＣＤ（Ｃｏｍｐａ
ｃｔ Ｄｉｓｃ）およびＤＶＤ（Ｄｉｇｉｔａｌ Ｖｅｒｓａｔｉｌｅ Ｄｉｓｃ）など
の光学ディスク、可搬性記憶媒体や、ＲＡＭ等の半導体メモリ素子等との間で種々の情報
を読み書きする駆動装置であってもよい。
は、ＨＤＤやＳＳＤ、集積回路記憶装置等の記憶装置である。メモリ１３は、記憶部に相
当する。なお、メモリ１３は、ＨＤＤやＳＳＤ等以外にも、ＲＡＭ（Ｒａｎｄｏｍ Ａｃ
ｃｅｓｓ Ｍｅｍｏｒｙ）、フラッシュメモリ等の半導体メモリ素子、ＣＤ（Ｃｏｍｐａ
ｃｔ Ｄｉｓｃ）およびＤＶＤ（Ｄｉｇｉｔａｌ Ｖｅｒｓａｔｉｌｅ Ｄｉｓｃ）など
の光学ディスク、可搬性記憶媒体や、ＲＡＭ等の半導体メモリ素子等との間で種々の情報
を読み書きする駆動装置であってもよい。
メモリ１３は、処理回路１５により実現される取得機能１５１、調整機能１５３、位置
合わせ機能１５５、および生成機能１５７を、コンピュータによって実行可能なプログラ
ムの形態で記憶する。メモリ１３は、取得機能１５１により通信インターフェース１１を
介して受信された各種データを記憶する。具体的には、メモリ１３は、例えば、取得機能
１５１により取得された複数のコイル感度分布と非カーテシアンデータとグリッディング
感度分布とを記憶する。コイル感度分布は、複数のコイルエレメント各々に対応する。す
なわち、メモリ１３は、複数のコイルエレメントに対応する複数のコイル感度分布を記憶
する。また、メモリ１３は、生成機能１５７により生成された磁気共鳴画像（以下、ＭＲ
画像と呼ぶ）を記憶する。
合わせ機能１５５、および生成機能１５７を、コンピュータによって実行可能なプログラ
ムの形態で記憶する。メモリ１３は、取得機能１５１により通信インターフェース１１を
介して受信された各種データを記憶する。具体的には、メモリ１３は、例えば、取得機能
１５１により取得された複数のコイル感度分布と非カーテシアンデータとグリッディング
感度分布とを記憶する。コイル感度分布は、複数のコイルエレメント各々に対応する。す
なわち、メモリ１３は、複数のコイルエレメントに対応する複数のコイル感度分布を記憶
する。また、メモリ１３は、生成機能１５７により生成された磁気共鳴画像（以下、ＭＲ
画像と呼ぶ）を記憶する。
グリッディング感度分布は、ｋ空間におけるＭＲデータの配置に関するグリッディング
の感度分布である。グリッディングは、非カーテシアンによる撮像で収集された非カーテ
シアンデータを補間して、ｋ空間の格子点に配置するために用いられる各種処理に相当す
る。グリッディングの条件は、当該各種処理において用いられる各種データ等に対する調
整、変更、修正などを実現するための要件および制約等に対応する。
の感度分布である。グリッディングは、非カーテシアンによる撮像で収集された非カーテ
シアンデータを補間して、ｋ空間の格子点に配置するために用いられる各種処理に相当す
る。グリッディングの条件は、当該各種処理において用いられる各種データ等に対する調
整、変更、修正などを実現するための要件および制約等に対応する。
グリッディング感度分布は、複素数のデータで表現される。グリッディング感度分布は
、非カーテシアンデータを畳み込むために用いられるカーネル関数のフーリエ変換に相当
する。すなわち、カーネル関数は、グリッディングの感度の元となる。カーネル関数とし
て、例えば、複素数で表現されたガウス関数を用いる。なお、カーネル関数は、ガウス関
数に限定されず他の関数であってもよい。また、カーネル関数は、メモリ１３に記憶され
てもよい。当該カーネル関数を用いて非カーテシアンデータを畳み込む畳み込み演算にお
いて、カーネル関数は、ｋ空間における非カーテシアンデータの位置と非カーテシアンデ
ータの値とをぼかす（Ｂｌｕｒ）ことに寄与する。また、当該畳み込み演算において、カ
ーネル関数の半値幅は、当該ぼかすことの程度に対応する。また、カーネル関数の位相勾
配が、グリッディング感度分布における位置に対応する。
、非カーテシアンデータを畳み込むために用いられるカーネル関数のフーリエ変換に相当
する。すなわち、カーネル関数は、グリッディングの感度の元となる。カーネル関数とし
て、例えば、複素数で表現されたガウス関数を用いる。なお、カーネル関数は、ガウス関
数に限定されず他の関数であってもよい。また、カーネル関数は、メモリ１３に記憶され
てもよい。当該カーネル関数を用いて非カーテシアンデータを畳み込む畳み込み演算にお
いて、カーネル関数は、ｋ空間における非カーテシアンデータの位置と非カーテシアンデ
ータの値とをぼかす（Ｂｌｕｒ）ことに寄与する。また、当該畳み込み演算において、カ
ーネル関数の半値幅は、当該ぼかすことの程度に対応する。また、カーネル関数の位相勾
配が、グリッディング感度分布における位置に対応する。
処理回路１５は、画像生成装置１の全体の制御を行う。処理回路１５は、上述のプロセ
ッサなどにより実現される。処理回路１５は、取得機能１５１、調整機能１５３、位置合
わせ機能１５５、および生成機能１５７などを備える。取得機能１５１、調整機能１５３
、位置合わせ機能１５５、および生成機能１５７をそれぞれ実現する処理回路１５は、取
得部、調整部、位置合わせ部、生成部に相当する。取得機能１５１、調整機能１５３、位
置合わせ機能１５５、および生成機能１５７などの各機能は、コンピュータによって実行
可能なプログラムの形態でメモリ１３に記憶されている。処理回路１５は、プロセッサで
ある。例えば、処理回路１５は、プログラムをメモリ１３から読み出し、実行することで
各プログラムに対応する機能を実現する。換言すると、各プログラムを読み出した状態の
処理回路１５は、取得機能１５１、調整機能１５３、位置合わせ機能１５５、および生成
機能１５７などの各機能を有することとなる。
ッサなどにより実現される。処理回路１５は、取得機能１５１、調整機能１５３、位置合
わせ機能１５５、および生成機能１５７などを備える。取得機能１５１、調整機能１５３
、位置合わせ機能１５５、および生成機能１５７をそれぞれ実現する処理回路１５は、取
得部、調整部、位置合わせ部、生成部に相当する。取得機能１５１、調整機能１５３、位
置合わせ機能１５５、および生成機能１５７などの各機能は、コンピュータによって実行
可能なプログラムの形態でメモリ１３に記憶されている。処理回路１５は、プロセッサで
ある。例えば、処理回路１５は、プログラムをメモリ１３から読み出し、実行することで
各プログラムに対応する機能を実現する。換言すると、各プログラムを読み出した状態の
処理回路１５は、取得機能１５１、調整機能１５３、位置合わせ機能１５５、および生成
機能１５７などの各機能を有することとなる。
処理回路１５は、取得機能１５１により、被検体Ｐの撮像に用いられる受信コイル１１
７の感度の分布を示すコイル感度分布と、ｋ空間における非カーテシアンでの被検体Ｐに
対する撮像により収集されたＭＲデータとを取得する。例えば、図１に示すように、処理
回路１５が単独の画像生成装置１に搭載される場合、取得機能１５１は、ネットワークお
よび通信インターフェース１１を介して、ＰＡＣＳやモダリティなどから、複数のコイル
エレメントに対応する複数のコイル感度分布と、ラジアル収集によるＭＲデータとを、ネ
ットワークを介して取得する。また、例えば、図２に示すように、処理回路１５がＭＲＩ
装置１００に搭載される場合、取得機能１５１は、生成機能１５７により生成された複数
のコイル感度分布と、受信回路１１９により生成されたＭＲデータとを、取得する。また
、受信コイル１１７が複数のコイルエレメントを有する場合、取得機能１５１は、複数の
コイルエレメントごとに、コイル感度分布と磁気共鳴データとを取得する。取得機能１５
１は、取得されたコイル感度分布とＭＲデータとを、メモリ１３に記憶させる。
７の感度の分布を示すコイル感度分布と、ｋ空間における非カーテシアンでの被検体Ｐに
対する撮像により収集されたＭＲデータとを取得する。例えば、図１に示すように、処理
回路１５が単独の画像生成装置１に搭載される場合、取得機能１５１は、ネットワークお
よび通信インターフェース１１を介して、ＰＡＣＳやモダリティなどから、複数のコイル
エレメントに対応する複数のコイル感度分布と、ラジアル収集によるＭＲデータとを、ネ
ットワークを介して取得する。また、例えば、図２に示すように、処理回路１５がＭＲＩ
装置１００に搭載される場合、取得機能１５１は、生成機能１５７により生成された複数
のコイル感度分布と、受信回路１１９により生成されたＭＲデータとを、取得する。また
、受信コイル１１７が複数のコイルエレメントを有する場合、取得機能１５１は、複数の
コイルエレメントごとに、コイル感度分布と磁気共鳴データとを取得する。取得機能１５
１は、取得されたコイル感度分布とＭＲデータとを、メモリ１３に記憶させる。
処理回路１５は、調整機能１５３により、コイル感度分布の深さに応じてグリッディン
グ感度分布が最適となるように、グリッディング感度分布を調整する。換言すれば、調整
機能１５３は、コイル感度分布に基づいて、グリッディング感度分布の元となるカーネル
関数の包絡形状または半値幅を調整する。具体的には、調整機能１５３は、受信コイル１
１７における複数のコイルエレメントごとに、当該コイルエレメントに対応するコイル感
度分布に基づいて、当該コイルエレメントに対応するカーネル関数を調整する。
グ感度分布が最適となるように、グリッディング感度分布を調整する。換言すれば、調整
機能１５３は、コイル感度分布に基づいて、グリッディング感度分布の元となるカーネル
関数の包絡形状または半値幅を調整する。具体的には、調整機能１５３は、受信コイル１
１７における複数のコイルエレメントごとに、当該コイルエレメントに対応するコイル感
度分布に基づいて、当該コイルエレメントに対応するカーネル関数を調整する。
例えば、調整機能１５３は、コイル感度分布における感度の広がりを計算し、計算され
た感度の広がりを所定の基準値と比較して、カーネル関数の半値幅を調整する。より詳細
には、調整機能１５３は、カーネル関数の調整として、複数のコイルエレメントごとすな
わち複数の収集チャネルごとに、コイル感度分布における感度の広がりを示す指標を用い
てカーネル関数の半値幅を調整する。なお、一つの収集チャネルは、一つのコイルエレメ
ントに対応することに限定されず、複数のコイルエレメントに対応付けられてもよい。複
数の所定の基準値および調整機能１５３における具体的な処理内容は、画像生成処理にお
いて説明する。なお、本調整機能１５３は、処理回路１５において非搭載であってもよい
。このとき、図１および図２に示す調整機能１５３は不要となり、グリッディングに用い
られるカーネル関数は、複数のコイルエレメントに亘って同一の形状となる。
た感度の広がりを所定の基準値と比較して、カーネル関数の半値幅を調整する。より詳細
には、調整機能１５３は、カーネル関数の調整として、複数のコイルエレメントごとすな
わち複数の収集チャネルごとに、コイル感度分布における感度の広がりを示す指標を用い
てカーネル関数の半値幅を調整する。なお、一つの収集チャネルは、一つのコイルエレメ
ントに対応することに限定されず、複数のコイルエレメントに対応付けられてもよい。複
数の所定の基準値および調整機能１５３における具体的な処理内容は、画像生成処理にお
いて説明する。なお、本調整機能１５３は、処理回路１５において非搭載であってもよい
。このとき、図１および図２に示す調整機能１５３は不要となり、グリッディングに用い
られるカーネル関数は、複数のコイルエレメントに亘って同一の形状となる。
処理回路１５は、位置合わせ機能１５５により、コイル感度分布に基づいて、グリッデ
ィング感度分布とコイル感度分布との位置合わせを実行する。また、受信コイル１１７が
複数のコイルエレメントを有する場合、位置合わせ機能１５５は、コイル感度分布に基づ
いて、グリッディング感度分布とコイル感度分布との位置合わせを実行する。例えば、位
置合わせ機能１５５は、グリッディング感度分布における感度の重心（以下、グリッディ
ング感度重心と呼ぶ）とコイル感度分布における感度の重心（以下、コイル感度重心と呼
ぶ）との位置合わせを実行する。
ィング感度分布とコイル感度分布との位置合わせを実行する。また、受信コイル１１７が
複数のコイルエレメントを有する場合、位置合わせ機能１５５は、コイル感度分布に基づ
いて、グリッディング感度分布とコイル感度分布との位置合わせを実行する。例えば、位
置合わせ機能１５５は、グリッディング感度分布における感度の重心（以下、グリッディ
ング感度重心と呼ぶ）とコイル感度分布における感度の重心（以下、コイル感度重心と呼
ぶ）との位置合わせを実行する。
当該位置合わせにより、位置合わせ機能１５５は、グリッディング感度重心とコイル感
度重心との差を算出する。当該差は、例えば、コイル感度重心からグリッディング感度重
心へ向かうベクトル（距離および方向）に相当する。すなわち、当該ベクトルの算出は、
グリッディング感度分布とコイル感度分布との位置合わせの結果に対応する。位置合わせ
機能１５５は、受信コイル１１７おける複数のコイルエレメント各々に関して、当該位置
合わせを実行する。位置合わせ機能１５５における具体的な処理内容は、非カーテシアン
データを用いてＭＲ画像を生成する処理（以下、画像生成処理と呼ぶ）において説明する
。
度重心との差を算出する。当該差は、例えば、コイル感度重心からグリッディング感度重
心へ向かうベクトル（距離および方向）に相当する。すなわち、当該ベクトルの算出は、
グリッディング感度分布とコイル感度分布との位置合わせの結果に対応する。位置合わせ
機能１５５は、受信コイル１１７おける複数のコイルエレメント各々に関して、当該位置
合わせを実行する。位置合わせ機能１５５における具体的な処理内容は、非カーテシアン
データを用いてＭＲ画像を生成する処理（以下、画像生成処理と呼ぶ）において説明する
。
処理回路１５は、生成機能１５７により、位置合わせ機能１５５による位置合わせの結
果と、ＭＲデータである非カーテシアンデータと、コイル感度分布と、グリッディング感
度分布とに基づいて、ＭＲ画像を生成する。調整機能１５３によりカーネル関数が調整さ
れている場合、生成機能１５７において用いられるグリッディング感度分布は、調整機能
１５３により調整されたカーネル関数のフーリエ変換に対応する。このとき、生成機能１
５７は、位置合わせの結果と磁気共鳴データとコイル感度分布と調整されたカーネル関数
に対応するグリッディング感度分布とに基づいて、ＭＲ画像を生成する。なお、受信コイ
ル１１７が複数のコイルエレメントを有する場合、生成機能１５７は、複数のコイルエレ
メントごとの複数の位置合わせの結果と複数の非カーテシアンデータと複数のコイル感度
分布と複数のグリッディング感度分布とに基づいて、ＭＲ画像を生成する。
果と、ＭＲデータである非カーテシアンデータと、コイル感度分布と、グリッディング感
度分布とに基づいて、ＭＲ画像を生成する。調整機能１５３によりカーネル関数が調整さ
れている場合、生成機能１５７において用いられるグリッディング感度分布は、調整機能
１５３により調整されたカーネル関数のフーリエ変換に対応する。このとき、生成機能１
５７は、位置合わせの結果と磁気共鳴データとコイル感度分布と調整されたカーネル関数
に対応するグリッディング感度分布とに基づいて、ＭＲ画像を生成する。なお、受信コイ
ル１１７が複数のコイルエレメントを有する場合、生成機能１５７は、複数のコイルエレ
メントごとの複数の位置合わせの結果と複数の非カーテシアンデータと複数のコイル感度
分布と複数のグリッディング感度分布とに基づいて、ＭＲ画像を生成する。
例えば、処理回路１５は、生成機能１５７により、グリッディング感度重心とコイル感
度重心との差（べクトル）を位置合わせの結果として用いて、非カーテシアンデータに対
して位相変調を実行する。すなわち、生成機能１５７は、当該位相変調として、非カーテ
シアンデータにおける位相に当該ベクトルに応じた位相の傾斜（勾配）を付与する位相変
調を実行する。非カーテシアンデータに対する位相変調は、グリッディングの条件の変更
に包含される。次いで、生成機能１５７は、グリッディングの感度の元となるカーネル関
数を用いて、位相変調後の非カーテシアンデータ（以下、位相変調後データと呼ぶ）をｋ
空間に配置する。具体的には、生成機能１５７は、カーネル関数を用いて位相変調後デー
タを畳み込む畳み込み演算を実行し、畳み込み演算の結果をｋ空間に配置する。
度重心との差（べクトル）を位置合わせの結果として用いて、非カーテシアンデータに対
して位相変調を実行する。すなわち、生成機能１５７は、当該位相変調として、非カーテ
シアンデータにおける位相に当該ベクトルに応じた位相の傾斜（勾配）を付与する位相変
調を実行する。非カーテシアンデータに対する位相変調は、グリッディングの条件の変更
に包含される。次いで、生成機能１５７は、グリッディングの感度の元となるカーネル関
数を用いて、位相変調後の非カーテシアンデータ（以下、位相変調後データと呼ぶ）をｋ
空間に配置する。具体的には、生成機能１５７は、カーネル関数を用いて位相変調後デー
タを畳み込む畳み込み演算を実行し、畳み込み演算の結果をｋ空間に配置する。
処理回路１５は、生成機能１５７により、ｋ空間に配置されたＭＲデータを用いて、Ｍ
Ｒ画像を生成する。具体的には、生成機能１５７は、ｋ空間に配置された畳み込み演算の
結果に対してフーリエ変換を実行する。生成機能１５７は、当該フーリエ変換の結果に対
して、当該差に基づいて上記位相変調とは逆の位相に対応するように平行移動を実行する
ことで、変換画像を生成する。なお、生成機能１５７は、当該差に基づいて、ｋ空間に配
置された非カーテシアンデータを、上記位相変調とは逆の位相に変調し、当該逆の位相に
変調された非カーテシアンデータに対してフーリエ変換を行うことにより、変換画像を生
成してもよい。
Ｒ画像を生成する。具体的には、生成機能１５７は、ｋ空間に配置された畳み込み演算の
結果に対してフーリエ変換を実行する。生成機能１５７は、当該フーリエ変換の結果に対
して、当該差に基づいて上記位相変調とは逆の位相に対応するように平行移動を実行する
ことで、変換画像を生成する。なお、生成機能１５７は、当該差に基づいて、ｋ空間に配
置された非カーテシアンデータを、上記位相変調とは逆の位相に変調し、当該逆の位相に
変調された非カーテシアンデータに対してフーリエ変換を行うことにより、変換画像を生
成してもよい。
処理回路１５は、生成機能１５７により、変換画像を用いてＭＲ画像を生成する。具体
的には、処理回路１５は、生成機能１５７により、位置合わせの結果を用いてグリッディ
ング感度分布とコイル感度分布とを合成することにより、合成感度分布を生成する。合成
感度分布は、複素数のデータで表現される。具体的には、生成機能１５７は、上記ベクト
ルとは逆のベクトル（以下、逆ベクトルと呼ぶ）を用いて、グリッディング感度分布を平
行移動する。次いで、生成機能１５７は、平行移動されたグリッディング感度分布とコイ
ル感度分布とを合成し、合成感度分布を生成する。続いて、生成機能１５７は、変換画像
と合成感度分布とを用いた感度エンコーディング（以下、ＳＥＮＳＥ（ＳＥＮＳｉｔｉｖ
ｉｔｙ Ｅｎｃｏｄｉｎｇ）と呼ぶ）により、ＭＲ画像を生成する。なお、生成機能１５
７は、変換画像の平方和の平方根（以下、ＳＯＳ（Ｓｕｍ Ｏｆ Ｓｑｕａｒｅ）と呼ぶ
）により、ＭＲ画像を生成してもよい。
的には、処理回路１５は、生成機能１５７により、位置合わせの結果を用いてグリッディ
ング感度分布とコイル感度分布とを合成することにより、合成感度分布を生成する。合成
感度分布は、複素数のデータで表現される。具体的には、生成機能１５７は、上記ベクト
ルとは逆のベクトル（以下、逆ベクトルと呼ぶ）を用いて、グリッディング感度分布を平
行移動する。次いで、生成機能１５７は、平行移動されたグリッディング感度分布とコイ
ル感度分布とを合成し、合成感度分布を生成する。続いて、生成機能１５７は、変換画像
と合成感度分布とを用いた感度エンコーディング（以下、ＳＥＮＳＥ（ＳＥＮＳｉｔｉｖ
ｉｔｙ Ｅｎｃｏｄｉｎｇ）と呼ぶ）により、ＭＲ画像を生成する。なお、生成機能１５
７は、変換画像の平方和の平方根（以下、ＳＯＳ（Ｓｕｍ Ｏｆ Ｓｑｕａｒｅ）と呼ぶ
）により、ＭＲ画像を生成してもよい。
なお、処理回路１５は、生成機能１５７により、変換画像をグリッディング感度分布で
除算することにより、変換画像におけるシェーディングを低減（補正）したディシェーデ
ィング画像を生成してもよい。このとき、生成機能１５７は、ディシェーディング画像と
コイル感度分布とを用いたＳＥＮＳＥ、またはディシェーディング画像のＳＯＳにより、
ＭＲ画像を生成する。生成機能１５７における具体的な処理内容は、画像生成処理におい
て説明する。
除算することにより、変換画像におけるシェーディングを低減（補正）したディシェーデ
ィング画像を生成してもよい。このとき、生成機能１５７は、ディシェーディング画像と
コイル感度分布とを用いたＳＥＮＳＥ、またはディシェーディング画像のＳＯＳにより、
ＭＲ画像を生成する。生成機能１５７における具体的な処理内容は、画像生成処理におい
て説明する。
以上のように構成された本実施形態のＭＲＩ装置１００により実行される画像生成処理
について、図３乃至図５を用いて説明する。図３および図４は、画像生成処理の概要を示
す図である。図５は、画像生成処理の手順の一例を示すフローチャートである。なお、画
像生成処理は、パラレルイメージング（ｐａｒａｌｌｅｌ ｉｍａｇｉｎｇ）に関する再
構成においても適宜利用可能である。
について、図３乃至図５を用いて説明する。図３および図４は、画像生成処理の概要を示
す図である。図５は、画像生成処理の手順の一例を示すフローチャートである。なお、画
像生成処理は、パラレルイメージング（ｐａｒａｌｌｅｌ ｉｍａｇｉｎｇ）に関する再
構成においても適宜利用可能である。
（画像生成処理）
図５における画像生成処理の前段階として、撮像制御回路１２１は、被検体Ｐに対して
コイル感度データの収集を実行する。処理回路１５は、生成機能１５７により、コイル感
度データに対してフーリエ変換を実行することで、コイル感度分布Ｓｃを生成する。また
、撮像制御回路１２１は、被検体Ｐに対して非カーテシアンでの撮像を実行し、非カーテ
シアンデータＮＣＤを収集する。コイル感度分布Ｓｃと非カーテシアンデータＮＣＤとは
、記憶装置１２５や、ＰＡＣＳなどに記憶される。
図５における画像生成処理の前段階として、撮像制御回路１２１は、被検体Ｐに対して
コイル感度データの収集を実行する。処理回路１５は、生成機能１５７により、コイル感
度データに対してフーリエ変換を実行することで、コイル感度分布Ｓｃを生成する。また
、撮像制御回路１２１は、被検体Ｐに対して非カーテシアンでの撮像を実行し、非カーテ
シアンデータＮＣＤを収集する。コイル感度分布Ｓｃと非カーテシアンデータＮＣＤとは
、記憶装置１２５や、ＰＡＣＳなどに記憶される。
（ステップＳ５０１）
処理回路１５は、取得機能１５１により、非カーテシアンデータＮＣＤとコイル感度分
布Ｓｃとを取得する。例えば、取得機能１５１は、非カーテシアンデータＮＣＤとコイル
感度分布Ｓｃとを、メモリ１３から取得する。なお、画像生成装置１が単独で用いられる
場合、取得機能１５１は、ネットワークおよび通信インターフェース１１を介して、外部
の記憶装置やＰＡＣＳなどから、非カーテシアンデータＮＣＤとコイル感度分布Ｓｃとを
取得する。取得機能１５１は、取得された非カーテシアンデータＮＣＤとコイル感度分布
Ｓｃとを、メモリ１３に記憶させる。
処理回路１５は、取得機能１５１により、非カーテシアンデータＮＣＤとコイル感度分
布Ｓｃとを取得する。例えば、取得機能１５１は、非カーテシアンデータＮＣＤとコイル
感度分布Ｓｃとを、メモリ１３から取得する。なお、画像生成装置１が単独で用いられる
場合、取得機能１５１は、ネットワークおよび通信インターフェース１１を介して、外部
の記憶装置やＰＡＣＳなどから、非カーテシアンデータＮＣＤとコイル感度分布Ｓｃとを
取得する。取得機能１５１は、取得された非カーテシアンデータＮＣＤとコイル感度分布
Ｓｃとを、メモリ１３に記憶させる。
（ステップＳ５０２）
処理回路１５は、調整機能１５３により、コイル感度分布Ｓｃに基づいて、グリッディ
ング感度分布の元となるカーネル関数を調整する（カーネル調整：ＫＡ）。具体的には、
調整機能１５３は、コイル感度分布Ｓｃにおける感度の広がりを示す指標を用いて、カー
ネル関数の半値幅を調整する。当該指標は、例えば、コイル感度分布における感度の広が
りの程度に相当する。具体的には、調整機能１５３は、例えば、以下の式により、当該指
標を、コイルエレメントｃに対応するコイル感度分布Ｓｃごとに計算する。
処理回路１５は、調整機能１５３により、コイル感度分布Ｓｃに基づいて、グリッディ
ング感度分布の元となるカーネル関数を調整する（カーネル調整：ＫＡ）。具体的には、
調整機能１５３は、コイル感度分布Ｓｃにおける感度の広がりを示す指標を用いて、カー
ネル関数の半値幅を調整する。当該指標は、例えば、コイル感度分布における感度の広が
りの程度に相当する。具体的には、調整機能１５３は、例えば、以下の式により、当該指
標を、コイルエレメントｃに対応するコイル感度分布Ｓｃごとに計算する。
上式（１）におけるｄｅｖ（ｘ）ｃおよびｄｅｖ（ｙ）ｃは、コイル感度分布Ｓｃにお
ける画像上の位置を表すベクトル
におけるｘ成分の当該指標およびｙ成分の当該指標にそれぞれ対応する。式（１）にお
ける右辺は、ｘ方向およびｙ方向における感度を重みとし、以下のコイル感度重心
を基準とした感度の標準偏差を、コイル感度分布Ｓｃにおける感度の積分で規格化した値
を示している。
ける画像上の位置を表すベクトル
ける右辺は、ｘ方向およびｙ方向における感度を重みとし、以下のコイル感度重心
を示している。
コイルエレメントｃにおけるコイル感度重心は、例えば、以下の式（２）
により算出される。式（２）の左辺は、コイル感度重心の位置ベクトルを示している。右
辺の分母は、コイル感度分布Ｓｃの位置ｒにおけるコイルの感度Ｓｃ（ｒ）をコイル感度
分布Ｓｃの全域に亘って積分することを示している。右辺の分子は、コイル感度分布Ｓｃ
の位置ｒにおけるコイルの感度Ｓｃ（ｒ）を重みとして位置ｒに掛けて、コイル感度分布
Ｓｃの全域に亘って積分することを示している。
辺の分母は、コイル感度分布Ｓｃの位置ｒにおけるコイルの感度Ｓｃ（ｒ）をコイル感度
分布Ｓｃの全域に亘って積分することを示している。右辺の分子は、コイル感度分布Ｓｃ
の位置ｒにおけるコイルの感度Ｓｃ（ｒ）を重みとして位置ｒに掛けて、コイル感度分布
Ｓｃの全域に亘って積分することを示している。
所定の基準値ｄｅｖｒｅｆにおける調整前のカーネル関数（以下、基準のカーネル関数
と呼ぶ）をＫｒｅｆ（ｋｘ、ｋｙ）とすると、処理回路１５は、生成機能１５７により、
算出された当該指標（ｄｅｖ（ｘ）ｃ、ｄｅｖ（ｙ）ｃ）と所定の基準値ｄｅｖｒｅｆと
基準のカーネル関数Ｋｒｅｆ（ｋｘ、ｋｙ）とに基づいて、コイルエレメントｃに関する
調整後のカーネル関数を、以下の式を用いて決定する。
と呼ぶ）をＫｒｅｆ（ｋｘ、ｋｙ）とすると、処理回路１５は、生成機能１５７により、
算出された当該指標（ｄｅｖ（ｘ）ｃ、ｄｅｖ（ｙ）ｃ）と所定の基準値ｄｅｖｒｅｆと
基準のカーネル関数Ｋｒｅｆ（ｋｘ、ｋｙ）とに基づいて、コイルエレメントｃに関する
調整後のカーネル関数を、以下の式を用いて決定する。
る位相の変更を示している。また、上式における右辺において、ネイピア数ｅの指数関数
の乗数は、調整前のカーネル関数の包絡形状の変更、すなわちカーネル関数の半値幅の調
整に対応する。所定の基準値ｄｅｖｒｅｆは、例えば、コイル感度分布Ｓｃの画像サイズ
の半分の長さ（画像の中心から画像の端部までの距離）である。コイル感度分布Ｓｃら算
出された重み付き標準偏差（ｄｅｖ（ｘ）ｃ、ｄｅｖ（ｙ）ｃ）は、画像サイズの半分を
上回ることはないため、式（３）に示す調整後のカーネル関数ＫＦの形状は、予めメモリ
１３に記憶された基準のカーネル関数Ｋｒｅｆを引き伸ばした形状となる。
処理回路１５は、調整機能１５３により、調整後のカーネル関数ＫＦに対してフーリエ
変換を実行することにより、コイルエレメントｃに応じた調整後のグリッディング感度分
布Ｓｇを生成する。これにより、グリッディング感度分布Ｓｇに関するカーネル関数の調
整は完了する。調整機能１５３は、コイルエレメントｃに対応するコイル感度分布Ｓｃに
応じて調整されたカーネル関数ＫＦと調整後のグリッディング感度分布Ｓｇとを、当該コ
イルエレメントｃと対応付けて、メモリ１３に記憶する。本ステップにおけるカーネル関
数の調整は、グリッディング条件の変更に包含される。以降の処理におけるグリッディン
グ感度分布およびカーネル関数は、本ステップにより調整されたものであるものとする。
変換を実行することにより、コイルエレメントｃに応じた調整後のグリッディング感度分
布Ｓｇを生成する。これにより、グリッディング感度分布Ｓｇに関するカーネル関数の調
整は完了する。調整機能１５３は、コイルエレメントｃに対応するコイル感度分布Ｓｃに
応じて調整されたカーネル関数ＫＦと調整後のグリッディング感度分布Ｓｇとを、当該コ
イルエレメントｃと対応付けて、メモリ１３に記憶する。本ステップにおけるカーネル関
数の調整は、グリッディング条件の変更に包含される。以降の処理におけるグリッディン
グ感度分布およびカーネル関数は、本ステップにより調整されたものであるものとする。
なお、処理回路１５において、当該調整機能１５３が搭載されない場合、本ステップは
省略される。このとき、メモリ１３には、全てのコイルエレメントに共通するカーネル関
数とグリッディング感度分布とが記憶されることとなる。
省略される。このとき、メモリ１３には、全てのコイルエレメントに共通するカーネル関
数とグリッディング感度分布とが記憶されることとなる。
（ステップＳ５０３）
処理回路１５は、位置合わせ機能１５５により、グリッディング感度分布Ｓｇとコイル
感度分布Ｓｃとの位置合わせを実行する。位置合わせ機能１５５は、グリッディング感度
分布Ｓｇとコイル感度分布Ｓｃとの位置合わせにより、例えば、コイル感度重心からグリ
ッディング感度重心へ向かうベクトル（以下、シフトベクトルと呼ぶ）を算出する。
処理回路１５は、位置合わせ機能１５５により、グリッディング感度分布Ｓｇとコイル
感度分布Ｓｃとの位置合わせを実行する。位置合わせ機能１５５は、グリッディング感度
分布Ｓｇとコイル感度分布Ｓｃとの位置合わせにより、例えば、コイル感度重心からグリ
ッディング感度重心へ向かうベクトル（以下、シフトベクトルと呼ぶ）を算出する。
以下、説明を簡便にするために、グリッディング感度重心が画像の中心であって、かつ
画像の原点が画像の中心であるものとする。このとき、シフトベクトルは、コイル感度重
心の位置ベクトルの逆ベクトルとなる。すなわち、シフトベクトルは、式（２）に示すコ
イルの感度重心の位置ベクトルにマイナスを付与したベクトルとなる
画像の原点が画像の中心であるものとする。このとき、シフトベクトルは、コイル感度重
心の位置ベクトルの逆ベクトルとなる。すなわち、シフトベクトルは、式（２）に示すコ
イルの感度重心の位置ベクトルにマイナスを付与したベクトルとなる
（ステップＳ５０４）
処理回路１５は、生成機能１５７により、位置合わせの結果によるシフトベクトルに基
づいて、非カーテシアンデータＮＣＤに対して位相変調（シフト演算：ＳＣＧ）を実行す
る。画像空間においてシフトベクトルによる画像のシフトは、フーリエ変換における公式
により、ｋ空間では、複素数変調すなわちシフトベクトルを勾配とする位相の勾配の付与
に対応する。
処理回路１５は、生成機能１５７により、位置合わせの結果によるシフトベクトルに基
づいて、非カーテシアンデータＮＣＤに対して位相変調（シフト演算：ＳＣＧ）を実行す
る。画像空間においてシフトベクトルによる画像のシフトは、フーリエ変換における公式
により、ｋ空間では、複素数変調すなわちシフトベクトルを勾配とする位相の勾配の付与
に対応する。
このため、生成機能１５７は、非カーテシアンデータＮＣＤに対して、位置合わせの結
果によるシフトベクトルに基づく以下の位相の勾配
を乗算する。これにより、生成機能１５７は、位相変調後データＰＭＤを生成する。再構
成対象のＭＲデータである非カーテシアンデータＮＣＤに対して位相変調を行う本ステッ
プにおける処理は、コイル感度分布Ｓｃに対応する受信コイル１１７の感度情報に基づい
て、グリッディング条件を変更することに対応する。
果によるシフトベクトルに基づく以下の位相の勾配
成対象のＭＲデータである非カーテシアンデータＮＣＤに対して位相変調を行う本ステッ
プにおける処理は、コイル感度分布Ｓｃに対応する受信コイル１１７の感度情報に基づい
て、グリッディング条件を変更することに対応する。
（ステップＳ５０５）
処理回路１５は、生成機能１５７により、カーネル関数ＫＦを用いて位相変調後の非カ
ーテシアンデータすなわち位相変調後データＰＭＤを畳み込む畳み込み演算ＣＣを実行し
て、ｋ空間のカーテシアン座標へ位相変調後データＰＭＤを再配置する。図３におけるＧ
Ｄは、ｋ空間のカーテシアン座標へ再配置された位相変調後データＰＭＤを示している。
図３に示すように、ｋ空間に配置された畳み込み演算ＣＣの結果を示す再配置済みデータ
ＧＤは、絶対値および位相、すなわち複素数に対応する２種のｋ空間データである。再配
置済みデータＧＤは、カーネル関数ＫＦによる畳み込み演算ＣＣ、すなわちカーネル関数
ＫＦの半値幅等に応じたぼかしにより、動径方向に垂直な偏角方向および動径方向に対し
て半値幅に応じた所定の幅を有するものとなる。
処理回路１５は、生成機能１５７により、カーネル関数ＫＦを用いて位相変調後の非カ
ーテシアンデータすなわち位相変調後データＰＭＤを畳み込む畳み込み演算ＣＣを実行し
て、ｋ空間のカーテシアン座標へ位相変調後データＰＭＤを再配置する。図３におけるＧ
Ｄは、ｋ空間のカーテシアン座標へ再配置された位相変調後データＰＭＤを示している。
図３に示すように、ｋ空間に配置された畳み込み演算ＣＣの結果を示す再配置済みデータ
ＧＤは、絶対値および位相、すなわち複素数に対応する２種のｋ空間データである。再配
置済みデータＧＤは、カーネル関数ＫＦによる畳み込み演算ＣＣ、すなわちカーネル関数
ＫＦの半値幅等に応じたぼかしにより、動径方向に垂直な偏角方向および動径方向に対し
て半値幅に応じた所定の幅を有するものとなる。
（ステップ５０６）
処理回路１５は、生成機能１５７により、ｋ空間に配置された再配置済みデータＧＤに
対してフーリエ変換ＦＴを実行する。これにより、図３に示す再構成画像ＲｅｃｏｎＩが
生成される。図３に示すように、再構成画像ＲｅｃｏｎＩにおいて、断面像の画像領域の
周辺部が暗くなるシェーディングという現象が発生する。生成機能１５７は、ｋ空間に配
置された磁気共鳴データすなわち再配置済みデータＧＤのフーリエ変換である再構成画像
ＲｅｃｏｎＩを、位置合わせ機能１５５により算出された差すなわちシフトベクトルに基
づいて、ステップＳ５０４において実施された位相変調とは逆の位相に対応するように平
行移動ＴＲを実行する。これにより、生成機能１５７は、変換画像Ｉｆｔを生成する。
処理回路１５は、生成機能１５７により、ｋ空間に配置された再配置済みデータＧＤに
対してフーリエ変換ＦＴを実行する。これにより、図３に示す再構成画像ＲｅｃｏｎＩが
生成される。図３に示すように、再構成画像ＲｅｃｏｎＩにおいて、断面像の画像領域の
周辺部が暗くなるシェーディングという現象が発生する。生成機能１５７は、ｋ空間に配
置された磁気共鳴データすなわち再配置済みデータＧＤのフーリエ変換である再構成画像
ＲｅｃｏｎＩを、位置合わせ機能１５５により算出された差すなわちシフトベクトルに基
づいて、ステップＳ５０４において実施された位相変調とは逆の位相に対応するように平
行移動ＴＲを実行する。これにより、生成機能１５７は、変換画像Ｉｆｔを生成する。
具体的には、生成機能１５７は、シフトベクトルの逆ベクトル（以下、逆シフトベクト
ルと呼ぶ）すなわちコイル感度重心の位置ベクトルである式（２）を用いて、再構成画像
ＲｅｃｏｎＩに対して平行移動ＴＲを実行することにより、変換画像Ｉｆｔを生成する。
ルと呼ぶ）すなわちコイル感度重心の位置ベクトルである式（２）を用いて、再構成画像
ＲｅｃｏｎＩに対して平行移動ＴＲを実行することにより、変換画像Ｉｆｔを生成する。
なお、変換画像Ｉｆｔの生成は、上記手順に限定されない。例えば、生成機能１５７は
、ｋ空間に配置された再配置済みデータＧＤを、ステップＳ５０４において実行された位
相変調とは逆の位相に変調する。当該逆の位相は、例えば、ネイピア数ｅにおける指数に
おいて、シフトベクトルの逆シフトベクトルを位相勾配として用いて表現される。具体的
には、当該逆の位相は、以下の式で表される。
、ｋ空間に配置された再配置済みデータＧＤを、ステップＳ５０４において実行された位
相変調とは逆の位相に変調する。当該逆の位相は、例えば、ネイピア数ｅにおける指数に
おいて、シフトベクトルの逆シフトベクトルを位相勾配として用いて表現される。具体的
には、当該逆の位相は、以下の式で表される。
生成機能１５７は、上式で示された位相を、ｋ空間に配置された再配置済みデータＧＤ
に乗算する。当該乗算は、フーリエ変換における公式により、画像空間では、逆シフトベ
クトルでの平行移動に相当する。生成機能１５７は、当該逆の位相に変調された再配置済
みデータＧＤに対してフーリエ変換ＦＴを行うことにより、変換画像Ｉｆｔを生成する。
に乗算する。当該乗算は、フーリエ変換における公式により、画像空間では、逆シフトベ
クトルでの平行移動に相当する。生成機能１５７は、当該逆の位相に変調された再配置済
みデータＧＤに対してフーリエ変換ＦＴを行うことにより、変換画像Ｉｆｔを生成する。
（ステップＳ５０７）
処理回路１５は、生成機能１５７により、位置合わせの結果を用いてグリッディング感
度分布Ｓｇとコイル感度分布Ｓｃとを合成することにより、合成感度分布Ｓｓを生成する
。具体的には、グリッディング感度重心とコイル感度重心とを位置合わせして、グリッデ
ィング感度分布Ｓｇとコイル感度分布Ｓｃと合成する。当該合成は、コイルエレメントご
とに実行される。当該合成は、例えば、以下の式に示すように、グリッディング感度分布
Ｓｇとコイル感度分布Ｓｃとの、画素ごとの積である。
処理回路１５は、生成機能１５７により、位置合わせの結果を用いてグリッディング感
度分布Ｓｇとコイル感度分布Ｓｃとを合成することにより、合成感度分布Ｓｓを生成する
。具体的には、グリッディング感度重心とコイル感度重心とを位置合わせして、グリッデ
ィング感度分布Ｓｇとコイル感度分布Ｓｃと合成する。当該合成は、コイルエレメントご
とに実行される。当該合成は、例えば、以下の式に示すように、グリッディング感度分布
Ｓｇとコイル感度分布Ｓｃとの、画素ごとの積である。
Ｓｓを示している。また、Ｓｇｃ（ｒ）は、コイルエレメントｃにおけるｒの位置のグリ
ッディングの感度を示している。Ｓｃｃ（ｒ）は、ｒの位置のコイルエレメントｃの感度
を示している。
なお、本ステップの処理は、ステップＳ５０４以降であれば、いずれの段階においても
実行可能である。以上のステップＳ５０１乃至５０７における各処理は、複数のコイルエ
レメント各々に対して実行される。
実行可能である。以上のステップＳ５０１乃至５０７における各処理は、複数のコイルエ
レメント各々に対して実行される。
（ステップＳ５０８）
処理回路１５は、生成機能１５７により、合成感度分布Ｓｓと変換画像Ｉｆｔとを用い
たＳＥＮＳＥによる合成（以下、ＳＥＮＳＥ合成と呼ぶ）により、ＭＲ画像Ｉｍｅｒｇｅ
を生成する。合成感度分布Ｓｓの生成は、複数のコイルエレメント各々に対して実行され
ている。このため、ＳＥＮＳＥ合成は、複数のコイルエレメントに対応する複数の合成感
度分布Ｓｓｃおよび複数の変換画像Ｉｆｔを用いる。なお、図４に示すように、生成機能
１５７は、合成感度分布Ｓｓを示す領域（以下、合成感度領域と呼ぶ）より、変換画像Ｉ
ｆｔを示す領域（以下、画像領域と呼ぶ）が広ければ、生成機能１５７は、ＳＥＮＳＥ合
成の前に、入出力インターフェース１７を介してユーザにより指定された画像領域に、変
換画像Ｉｆｔの領域と合成感度領域とを合わせるように、変換画像Ｉｆｔと合成感度分布
Ｓｓとに対して、クロッピング（Ｃｒｏｐｐｉｎｇ）を実行する。
処理回路１５は、生成機能１５７により、合成感度分布Ｓｓと変換画像Ｉｆｔとを用い
たＳＥＮＳＥによる合成（以下、ＳＥＮＳＥ合成と呼ぶ）により、ＭＲ画像Ｉｍｅｒｇｅ
を生成する。合成感度分布Ｓｓの生成は、複数のコイルエレメント各々に対して実行され
ている。このため、ＳＥＮＳＥ合成は、複数のコイルエレメントに対応する複数の合成感
度分布Ｓｓｃおよび複数の変換画像Ｉｆｔを用いる。なお、図４に示すように、生成機能
１５７は、合成感度分布Ｓｓを示す領域（以下、合成感度領域と呼ぶ）より、変換画像Ｉ
ｆｔを示す領域（以下、画像領域と呼ぶ）が広ければ、生成機能１５７は、ＳＥＮＳＥ合
成の前に、入出力インターフェース１７を介してユーザにより指定された画像領域に、変
換画像Ｉｆｔの領域と合成感度領域とを合わせるように、変換画像Ｉｆｔと合成感度分布
Ｓｓとに対して、クロッピング（Ｃｒｏｐｐｉｎｇ）を実行する。
ＳＥＮＳＥ合成は、例えば、以下の式（４）
により、実現される。式（４）における左辺Ｉｍｅｒｇｅ（ｒ）は、画像上の位置ｒに
おいて、画像上の位置ｒにおいて、生成機能１５７により生成されるＭＲ画像Ｉｍｅｒｇ
ｅの画素値を示している。式（４）の右辺におけるＳｓｃ（ｒ）は、コイルエレメントｃ
における位置ｒでの合成感度分布Ｓｓの合成感度を示し、記号＊は、複素共役を示してい
る。
おいて、画像上の位置ｒにおいて、生成機能１５７により生成されるＭＲ画像Ｉｍｅｒｇ
ｅの画素値を示している。式（４）の右辺におけるＳｓｃ（ｒ）は、コイルエレメントｃ
における位置ｒでの合成感度分布Ｓｓの合成感度を示し、記号＊は、複素共役を示してい
る。
また、式（４）の右辺におけるＩｆｔｃ（ｒ）は、コイルエレメントｃにおける位置ｒ
での変換画像Ｉｆｔの画素値を示している。式（４）の右辺の分母は、位置ｒにおいて、
コイルエレメントごとの合成感度分布Ｓｓｃの絶対値の２乗を全てのコイルエレメントに
亘って加算した値を示している。また、式（４）の右辺の分子は、位置ｒにおいて、変換
画像Ｉｆｔとコイルエレメントごとの合成感度分布Ｓｓｃをとの積を全てのコイルエレメ
ントに亘って加算した値、換言すれば、位置ｒにおける変換画像Ｉｆｔに対する合成感度
分布Ｓｓｃの、全コイルに亘る重み付け加算を示している。すなわち、式（４）は、当該
重み付け加算を合成感度で規格化した値をＭＲ画像Ｉｍｅｒｇｅとして算出ことを意味し
ている。
での変換画像Ｉｆｔの画素値を示している。式（４）の右辺の分母は、位置ｒにおいて、
コイルエレメントごとの合成感度分布Ｓｓｃの絶対値の２乗を全てのコイルエレメントに
亘って加算した値を示している。また、式（４）の右辺の分子は、位置ｒにおいて、変換
画像Ｉｆｔとコイルエレメントごとの合成感度分布Ｓｓｃをとの積を全てのコイルエレメ
ントに亘って加算した値、換言すれば、位置ｒにおける変換画像Ｉｆｔに対する合成感度
分布Ｓｓｃの、全コイルに亘る重み付け加算を示している。すなわち、式（４）は、当該
重み付け加算を合成感度で規格化した値をＭＲ画像Ｉｍｅｒｇｅとして算出ことを意味し
ている。
本実施形態の変形例として、本ステップＳ５０８において、処理回路１５は、生成機能
１５７により、合成感度分布Ｓｓを用いずに、コイルエレメントｃにおける位置ｒでの変
換画像Ｉｆｔの画素値Ｉｆｔｃ（ｒ）を用いたＳＯＳにより、ＭＲ画像Ｉｍｅｒｇｅを生
成してもよい。このとき、ステップＳ５０７は不要となる。ＳＯＳによるＭＲ画像Ｉｍｅ
ｒｇｅの生成は、例えば以下の式により実現される。
１５７により、合成感度分布Ｓｓを用いずに、コイルエレメントｃにおける位置ｒでの変
換画像Ｉｆｔの画素値Ｉｆｔｃ（ｒ）を用いたＳＯＳにより、ＭＲ画像Ｉｍｅｒｇｅを生
成してもよい。このとき、ステップＳ５０７は不要となる。ＳＯＳによるＭＲ画像Ｉｍｅ
ｒｇｅの生成は、例えば以下の式により実現される。
Ｓｇとコイル感度分布Ｓｃとは用いることなく、ＭＲ画像Ｉｍｅｒｇｅが生成される。
なお、本変形例のさらなる変形例として、ＳＯＳの計算において用いられる画素値Ｉｆ
ｔｃ（ｒ）の代わりに、画素値Ｉｆｔｃ（ｒ）を、コイルエレメントｃにおける位置ｒで
のグリッディング感度分布Ｓｇの画素値Ｓｇｃ（ｒ）で除算したディシェーディング画像
が用いられてもよい。生成機能１７３は、コイルエレメントｃにおけるディシェーディン
グ画像Ｉｄｃ（ｒ）を、例えば以下の式により算出する。
ｔｃ（ｒ）の代わりに、画素値Ｉｆｔｃ（ｒ）を、コイルエレメントｃにおける位置ｒで
のグリッディング感度分布Ｓｇの画素値Ｓｇｃ（ｒ）で除算したディシェーディング画像
が用いられてもよい。生成機能１７３は、コイルエレメントｃにおけるディシェーディン
グ画像Ｉｄｃ（ｒ）を、例えば以下の式により算出する。
を用いたＳＯＳにより、ＭＲ画像Ｉｍｅｒｇｅを生成する。
なお、調整機能１５３によりグリッディング感度分布Ｓｇに関するカーネル関数が調整
されない場合、グリッディング感度分布Ｓｇは、コイルエレメントによらない定数（Ｓｇ
ｃ（ｒ）＝Ｓｇ（ｒ））となる。このとき、合成感度分布Ｓｓを用いたＳＥＮＳＥ合成に
より生成されるＭＲ画像Ｉｍｅｒｇｅと、ディシェーディング画像Ｉｄｃ（ｒ）とコイル
感度分布Ｓｃｃとを用いたＳＥＮＳＥ合成により生成されたＭＲ画像Ｉｄｃｍｅｒｇｅと
は、以下の式に示すように等価となる。
されない場合、グリッディング感度分布Ｓｇは、コイルエレメントによらない定数（Ｓｇ
ｃ（ｒ）＝Ｓｇ（ｒ））となる。このとき、合成感度分布Ｓｓを用いたＳＥＮＳＥ合成に
より生成されるＭＲ画像Ｉｍｅｒｇｅと、ディシェーディング画像Ｉｄｃ（ｒ）とコイル
感度分布Ｓｃｃとを用いたＳＥＮＳＥ合成により生成されたＭＲ画像Ｉｄｃｍｅｒｇｅと
は、以下の式に示すように等価となる。
以上に述べた第１の実施形態に係るＭＲＩ装置１００および画像生成装置１によれば、
コイル感度分布Ｓｃと非カーテシアンデータＮＣＤとを取得し、コイル感度分布Ｓｃに基
づいてグリッディング感度分布Ｓｇとコイル感度分布Ｓｃとの位置合わせを実行し、位置
合わせの結果と非カーテシアンデータＮＣＤとコイル感度分布Ｓｃとグリッディング感度
分布Ｓｇとに基づいてＭＲ画像Ｉｍｅｒｇｅを生成する。具体的には、本ＭＲＩ装置１０
０および本画像生成装置１によれば、グリッディング感度重心とコイル感度重心との差に
相当するシフトベクトルを用いて非カーテシアンデータＮＣＤに対して位相変調を実行し
、グリッディングの感度の元となるカーネル関数ＫＦを用いて位相変調後の磁気共鳴デー
タをｋ空間に配置し、ｋ空間に配置された磁気共鳴データ（再配置済みデータＧＤ）を用
いてＭＲ画像Ｉｍｅｒｇｅを生成する。
コイル感度分布Ｓｃと非カーテシアンデータＮＣＤとを取得し、コイル感度分布Ｓｃに基
づいてグリッディング感度分布Ｓｇとコイル感度分布Ｓｃとの位置合わせを実行し、位置
合わせの結果と非カーテシアンデータＮＣＤとコイル感度分布Ｓｃとグリッディング感度
分布Ｓｇとに基づいてＭＲ画像Ｉｍｅｒｇｅを生成する。具体的には、本ＭＲＩ装置１０
０および本画像生成装置１によれば、グリッディング感度重心とコイル感度重心との差に
相当するシフトベクトルを用いて非カーテシアンデータＮＣＤに対して位相変調を実行し
、グリッディングの感度の元となるカーネル関数ＫＦを用いて位相変調後の磁気共鳴デー
タをｋ空間に配置し、ｋ空間に配置された磁気共鳴データ（再配置済みデータＧＤ）を用
いてＭＲ画像Ｉｍｅｒｇｅを生成する。
また、本ＭＲＩ装置１００および本画像生成装置１によれば、ｋ空間に配置された磁気
共鳴データ（再配置済みデータＧＤ）のフーリエ変換である再構成画像ＲｅｃｏｎＩに対
して、コイル感度重心とグリッディング感度重心との差に基づいて当該位相変調とは逆の
位相に対応するように平行移動ＴＲを実行することで、変換画像Ｉｆｔを生成し、当該変
換画像Ｉｆｔを用いてＭＲ画像Ｉｍｅｒｇｅを生成する。また、本ＭＲＩ装置１００およ
び本画像生成装置１によれば、当該差に基づいて、ｋ空間に配置された磁気共鳴データ（
再配置済みデータＧＤ）を位相変調とは逆の位相に変調し、当該逆の位相に変調された磁
気共鳴データに対してフーリエ変換ＦＴを行うことにより変換画像Ｉｆｔを生成し、当該
変換画像Ｉｆｔを用いてＭＲ画像Ｉｍｅｒｇｅを生成する。
共鳴データ（再配置済みデータＧＤ）のフーリエ変換である再構成画像ＲｅｃｏｎＩに対
して、コイル感度重心とグリッディング感度重心との差に基づいて当該位相変調とは逆の
位相に対応するように平行移動ＴＲを実行することで、変換画像Ｉｆｔを生成し、当該変
換画像Ｉｆｔを用いてＭＲ画像Ｉｍｅｒｇｅを生成する。また、本ＭＲＩ装置１００およ
び本画像生成装置１によれば、当該差に基づいて、ｋ空間に配置された磁気共鳴データ（
再配置済みデータＧＤ）を位相変調とは逆の位相に変調し、当該逆の位相に変調された磁
気共鳴データに対してフーリエ変換ＦＴを行うことにより変換画像Ｉｆｔを生成し、当該
変換画像Ｉｆｔを用いてＭＲ画像Ｉｍｅｒｇｅを生成する。
また、本ＭＲＩ装置１００および本画像生成装置１によれば、位置合わせの結果を用い
てグリッディング感度分布Ｓｇとコイル感度分布Ｓｃとを合成することにより、合成感度
分布Ｓｓを生成し、変換画像Ｉｆｔと合成感度分布Ｓｓとを用いた感度エンコーディング
（ＳＥＮＳＥ）により、ＭＲ画像Ｉｍｅｒｇｅを生成する。また、本ＭＲＩ装置１００お
よび本画像生成装置１によれば、変換画像Ｉｆｔの平方和の平方根（ＳＯＳ）により、Ｍ
Ｒ画像Ｉｍｅｒｇｅを生成する。
てグリッディング感度分布Ｓｇとコイル感度分布Ｓｃとを合成することにより、合成感度
分布Ｓｓを生成し、変換画像Ｉｆｔと合成感度分布Ｓｓとを用いた感度エンコーディング
（ＳＥＮＳＥ）により、ＭＲ画像Ｉｍｅｒｇｅを生成する。また、本ＭＲＩ装置１００お
よび本画像生成装置１によれば、変換画像Ｉｆｔの平方和の平方根（ＳＯＳ）により、Ｍ
Ｒ画像Ｉｍｅｒｇｅを生成する。
図６は、本ＭＲＩ装置１００および本画像生成装置１による画像生成処理の結果と比較
例とにおいて、生成されたＭＲ画像と、当該ＭＲ画像と正解画像との差の一例を示す図で
ある。図６に示すように、本実施形態にかかる画像生成処理の結果は、比較例と比べて、
ＭＲ画像における断面像の画像領域の周囲におけるストリークアーチファクトが改善され
ている。加えて、図６に示すように、本実施形態にかかる画像生成処理の結果は、比較例
と比べて、当該画像領域におけるノイズの低減と、ＭＲ画像の４隅におけるアーチファク
ト除去がされている。
例とにおいて、生成されたＭＲ画像と、当該ＭＲ画像と正解画像との差の一例を示す図で
ある。図６に示すように、本実施形態にかかる画像生成処理の結果は、比較例と比べて、
ＭＲ画像における断面像の画像領域の周囲におけるストリークアーチファクトが改善され
ている。加えて、図６に示すように、本実施形態にかかる画像生成処理の結果は、比較例
と比べて、当該画像領域におけるノイズの低減と、ＭＲ画像の４隅におけるアーチファク
ト除去がされている。
以上のことから、本実施形態にかかるＭＲＩ装置１００および画像生成装置１によれば
、コイル感度分布Ｓｃとグリッディング感度分布Ｓｇの位置合わせにより、受信コイル１
１７の感度情報をグリッディングに反映させること、すなわちコイル感度分布Ｓｃにグリ
ッディング感度分布Ｓｇを合わせること、換言すればコイル感度分布Ｓｃとグリッディン
グ感度分布Ｓｇとにおける高感度の部分を合わせることができる。このため、本ＭＲＩ装
置１００および本画像生成装置１によれば、収集速度が速い非カーテシアンでの撮像に伴
うグリッディング再構成において、ストリークと呼ばれる放射状のアーチファクトなどを
低減すること、およびＳ／Ｎを向上することができる。すなわち、本実施形態にかかるＭ
ＲＩ装置１００および画像生成装置１によれば、グリッディングの感度と感度マップにお
ける受信コイル１１７の感度とのずれによるストリークなどのアーチファクトの発生範囲
を抑制し、かつ診断能の高いＭＲ画像Ｉｍｅｒｇｅを得ることができる。
、コイル感度分布Ｓｃとグリッディング感度分布Ｓｇの位置合わせにより、受信コイル１
１７の感度情報をグリッディングに反映させること、すなわちコイル感度分布Ｓｃにグリ
ッディング感度分布Ｓｇを合わせること、換言すればコイル感度分布Ｓｃとグリッディン
グ感度分布Ｓｇとにおける高感度の部分を合わせることができる。このため、本ＭＲＩ装
置１００および本画像生成装置１によれば、収集速度が速い非カーテシアンでの撮像に伴
うグリッディング再構成において、ストリークと呼ばれる放射状のアーチファクトなどを
低減すること、およびＳ／Ｎを向上することができる。すなわち、本実施形態にかかるＭ
ＲＩ装置１００および画像生成装置１によれば、グリッディングの感度と感度マップにお
ける受信コイル１１７の感度とのずれによるストリークなどのアーチファクトの発生範囲
を抑制し、かつ診断能の高いＭＲ画像Ｉｍｅｒｇｅを得ることができる。
（第２の実施形態）
第１の実施形態との相違は、非カーテシアンデータに対して位相変調を行うのではなく
、調整後のカーネル関数ＫＦに対して位相変調を実行し、位相変調後のカーネル関数を用
いて非カーテシアンデータＮＣＤをｋ空間に配置し、ｋ空間に配置されたデータを用いて
ＭＲ画像Ｉｍｅｒｇｅを生成することにある。なお、調整機能１５３によりグリッディン
グ感度分布に関するカーネル関数が調整されない場合、位相変調の対象は、例えばメモリ
１３に記憶されたカーネル関数またはグリッディング感度分布のフーリエ変換となる。以
下、説明を具体的にするために、本実施形態における位相変調の対象は、調整機能１５３
によりコイルエレメントごとに調整されたカーネル関数であるものとする。
第１の実施形態との相違は、非カーテシアンデータに対して位相変調を行うのではなく
、調整後のカーネル関数ＫＦに対して位相変調を実行し、位相変調後のカーネル関数を用
いて非カーテシアンデータＮＣＤをｋ空間に配置し、ｋ空間に配置されたデータを用いて
ＭＲ画像Ｉｍｅｒｇｅを生成することにある。なお、調整機能１５３によりグリッディン
グ感度分布に関するカーネル関数が調整されない場合、位相変調の対象は、例えばメモリ
１３に記憶されたカーネル関数またはグリッディング感度分布のフーリエ変換となる。以
下、説明を具体的にするために、本実施形態における位相変調の対象は、調整機能１５３
によりコイルエレメントごとに調整されたカーネル関数であるものとする。
以下、図７および図８を用いて、本実施形態に係る画像生成処理の手順について説明す
る。図７は、本実施形態における画像生成処理の一部の概要を示す図である。図７に示す
画像生成処理の後段の処理は、図４と同様なため、不図示としている。図８は、本実施形
態に係る画像生成処理の手順の一例を示すフローチャートである。
る。図７は、本実施形態における画像生成処理の一部の概要を示す図である。図７に示す
画像生成処理の後段の処理は、図４と同様なため、不図示としている。図８は、本実施形
態に係る画像生成処理の手順の一例を示すフローチャートである。
（画像生成処理）
（ステップＳ８０１）
処理回路１５は、取得機能１５１により、非カーテシアンデータＮＣＤとコイル感度分
布Ｓｃとを取得する。本ステップの処理内容は、ステップＳ５０１と同様なため、説明は
省略する。
（ステップＳ８０１）
処理回路１５は、取得機能１５１により、非カーテシアンデータＮＣＤとコイル感度分
布Ｓｃとを取得する。本ステップの処理内容は、ステップＳ５０１と同様なため、説明は
省略する。
（ステップＳ８０２）
処理回路１５は、調整機能１５３により、コイル感度分布Ｓｃに基づいて、グリッディ
ング感度分布の元となるカーネル関数を調整する。本ステップの処理内容は、ステップＳ
５０２と同様なため、説明は省略する。カーネル調整ＫＡに調整されたカーネル関数ＫＦ
は、図７において点線で示されている。
処理回路１５は、調整機能１５３により、コイル感度分布Ｓｃに基づいて、グリッディ
ング感度分布の元となるカーネル関数を調整する。本ステップの処理内容は、ステップＳ
５０２と同様なため、説明は省略する。カーネル調整ＫＡに調整されたカーネル関数ＫＦ
は、図７において点線で示されている。
（ステップＳ８０３）
処理回路１５は、位置合わせ機能１５５により、調整前のグリッディング感度分布とコ
イル感度分布Ｓｃとの位置合わせを実行する。位置合わせ機能１５５は、調整前のグリッ
ディング感度分布とコイル感度分布Ｓｃとの位置合わせにより、例えば、コイル感度重心
からグリッディング感度重心へ向かうシフトベクトルを算出する。本ステップの処理内容
は、ステップＳ５０３と同様なため、説明は省略する。
処理回路１５は、位置合わせ機能１５５により、調整前のグリッディング感度分布とコ
イル感度分布Ｓｃとの位置合わせを実行する。位置合わせ機能１５５は、調整前のグリッ
ディング感度分布とコイル感度分布Ｓｃとの位置合わせにより、例えば、コイル感度重心
からグリッディング感度重心へ向かうシフトベクトルを算出する。本ステップの処理内容
は、ステップＳ５０３と同様なため、説明は省略する。
（ステップＳ８０４）
処理回路１５は、生成機能１５７により、位置合わせの結果によるシフトベクトルの逆
ベクトル（逆シフトベクトル）に基づいて、調整後のカーネル関数ＫＦに対して位相変調
を実行する。調整後のカーネル関数ＫＦに対する位相変調は、当該カーネル関数ＫＦにお
けるピーク位置の平行移動に相当する。換言すれば、当該位相変調は、当該カーネル関数
に対応するグリッディング感度分布において、グリッディング感度重心をコイル感度重心
へシフト（平行移動）することに相当する。すなわち、画像であるグリッディング感度分
布のシフトは、フーリエ変換における公式により、調整後のカーネル関数が配置されるｋ
空間では、複素数変調すなわち逆シフトベクトルを勾配とする位相の勾配の付与に対応す
る。
処理回路１５は、生成機能１５７により、位置合わせの結果によるシフトベクトルの逆
ベクトル（逆シフトベクトル）に基づいて、調整後のカーネル関数ＫＦに対して位相変調
を実行する。調整後のカーネル関数ＫＦに対する位相変調は、当該カーネル関数ＫＦにお
けるピーク位置の平行移動に相当する。換言すれば、当該位相変調は、当該カーネル関数
に対応するグリッディング感度分布において、グリッディング感度重心をコイル感度重心
へシフト（平行移動）することに相当する。すなわち、画像であるグリッディング感度分
布のシフトは、フーリエ変換における公式により、調整後のカーネル関数が配置されるｋ
空間では、複素数変調すなわち逆シフトベクトルを勾配とする位相の勾配の付与に対応す
る。
すなわち、生成機能１５７は、調整後のカーネル関数ＫＦに対して、位置合わせの結果
による逆シフトベクトルに基づく以下の位相の勾配
を乗算する。これにより、生成機能１５７は、位相変調後のカーネル関数ＰＭＫＦを生成
する。本ステップにおけるカーネル関数に対する位相変調は、例えば、複数のコイルエレ
メントに対応する調整後の複数のカーネル関数ごとに実施される。生成機能１５７は、位
相変調後のカーネル関数ＰＭＫＦ対してフーリエ変換ＦＴを実行し、位相変調後のカーネ
ル関数に対応する調整後のグリッディング感度分布Ｓｇを生成する。本ステップにおける
カーネル関数ＫＦに対する位相変調は、グリッディングの条件の変更に包含される。
による逆シフトベクトルに基づく以下の位相の勾配
する。本ステップにおけるカーネル関数に対する位相変調は、例えば、複数のコイルエレ
メントに対応する調整後の複数のカーネル関数ごとに実施される。生成機能１５７は、位
相変調後のカーネル関数ＰＭＫＦ対してフーリエ変換ＦＴを実行し、位相変調後のカーネ
ル関数に対応する調整後のグリッディング感度分布Ｓｇを生成する。本ステップにおける
カーネル関数ＫＦに対する位相変調は、グリッディングの条件の変更に包含される。
（ステップＳ８０５）
処理回路１５は、生成機能１５７により、位相変調後のカーネル関数ＰＭＫＦを用いて
非カーテシアンデータＮＣＤを畳み込む畳み込み演算ＡＣＣを実行して、ｋ空間のカーテ
シアン座標へ非カーテシアンデータＮＣＤを再配置する。図７におけるＡＧＤは、ｋ空間
のカーテシアン座標へ再配置された磁気共鳴データを示している。図７に示すように、ｋ
空間に配置された畳み込み演算ＡＣＣの結果を示す再配置済みデータＡＧＤは、絶対値お
よび位相、すなわち複素数に対応する２つのｋ空間データである。再配置済みデータＡＧ
Ｄは、位相変調後のカーネル関数ＰＭＫＦによる畳み込み演算ＡＣＣ、すなわちカーネル
関数ＰＭＫＦの半値幅等に応じたぼかしにより、動径方向に垂直な偏角方向および動径方
向に対して半値幅に応じた所定の幅を有するものとなる。
処理回路１５は、生成機能１５７により、位相変調後のカーネル関数ＰＭＫＦを用いて
非カーテシアンデータＮＣＤを畳み込む畳み込み演算ＡＣＣを実行して、ｋ空間のカーテ
シアン座標へ非カーテシアンデータＮＣＤを再配置する。図７におけるＡＧＤは、ｋ空間
のカーテシアン座標へ再配置された磁気共鳴データを示している。図７に示すように、ｋ
空間に配置された畳み込み演算ＡＣＣの結果を示す再配置済みデータＡＧＤは、絶対値お
よび位相、すなわち複素数に対応する２つのｋ空間データである。再配置済みデータＡＧ
Ｄは、位相変調後のカーネル関数ＰＭＫＦによる畳み込み演算ＡＣＣ、すなわちカーネル
関数ＰＭＫＦの半値幅等に応じたぼかしにより、動径方向に垂直な偏角方向および動径方
向に対して半値幅に応じた所定の幅を有するものとなる。
（ステップＳ８０６）
処理回路１５は、生成機能１５７により、ｋ空間に配置された再配置済みデータＡＧＤ
に対してフーリエ変換ＦＴを実行する。当該フーリエ変換ＦＴにより、生成機能１５７は
、図７に示す変換画像Ｉｆｔを生成する。図７に示す変換画像Ｉｆｔと図３に示す変換画
像Ｉｆｔとは略同一な画像となる。
処理回路１５は、生成機能１５７により、ｋ空間に配置された再配置済みデータＡＧＤ
に対してフーリエ変換ＦＴを実行する。当該フーリエ変換ＦＴにより、生成機能１５７は
、図７に示す変換画像Ｉｆｔを生成する。図７に示す変換画像Ｉｆｔと図３に示す変換画
像Ｉｆｔとは略同一な画像となる。
（ステップＳ８０７）
処理回路１５は、生成機能１５７により、位相変調後のカーネル関数ＰＭＫＦのフーリ
エ変換であるグリッディング感度分布Ｓｇとコイル感度分布Ｓｃとを合成することにより
、合成感度分布Ｓｓを生成する。具体的には、生成機能１５７は、グリッディング感度分
布Ｓｇとコイル感度分布Ｓｃとを、画素ごとに乗算することにより、合成感度分布Ｓｓを
生成する。合成感度分布Ｓｓの生成に関する演算は、ステップＳ５０７と略同一なため、
説明は省略する。
処理回路１５は、生成機能１５７により、位相変調後のカーネル関数ＰＭＫＦのフーリ
エ変換であるグリッディング感度分布Ｓｇとコイル感度分布Ｓｃとを合成することにより
、合成感度分布Ｓｓを生成する。具体的には、生成機能１５７は、グリッディング感度分
布Ｓｇとコイル感度分布Ｓｃとを、画素ごとに乗算することにより、合成感度分布Ｓｓを
生成する。合成感度分布Ｓｓの生成に関する演算は、ステップＳ５０７と略同一なため、
説明は省略する。
（ステップＳ８０８）
処理回路１５は、生成機能１５７により、合成感度分布Ｓｓと変換画像Ｉｆｔとを用い
たＳＥＮＳＥ合成により、ＭＲ画像Ｉｍｅｒｇｅを生成する。なお、生成機能１５７は、
ＳＥＮＳＥ合成の代わりに、ＳＯＳなどによりＭＲ画像Ｉｍｅｒｇｅを生成してもよい。
本ステップの処理内容は、ステップＳ５０８と同様なため、説明は省略する。
処理回路１５は、生成機能１５７により、合成感度分布Ｓｓと変換画像Ｉｆｔとを用い
たＳＥＮＳＥ合成により、ＭＲ画像Ｉｍｅｒｇｅを生成する。なお、生成機能１５７は、
ＳＥＮＳＥ合成の代わりに、ＳＯＳなどによりＭＲ画像Ｉｍｅｒｇｅを生成してもよい。
本ステップの処理内容は、ステップＳ５０８と同様なため、説明は省略する。
以上に述べた第２の実施形態に係るＭＲＩ装置１００および画像生成装置１によれば、
グリッディング感度重心とコイル感度重心との差を位置合わせの結果として用いて、グリ
ッディング感度分布の元となるカーネル関数ＫＦに対して位相変調を実行し、位相変調後
のカーネル関数ＰＭＫＦを用いて磁気共鳴データ（非カーテシアンデータ）ＮＣＤをｋ空
間に配置し、ｋ空間に配置された磁気共鳴データ（再配置済みデータ）ＡＧＤを用いて、
ＭＲ画像Ｉｍｅｒｇｅを生成する。これにより、本実施形態におけるＭＲＩ装置１００お
よび画像生成装置１によれば、第１の実施形態における画像生成処理よりも簡便にＭＲ画
像Ｉｍｅｒｇｅを生成することができる。なお、本実施形態における他の効果は、第１の
実施形態と同様なため、説明は省略する。
グリッディング感度重心とコイル感度重心との差を位置合わせの結果として用いて、グリ
ッディング感度分布の元となるカーネル関数ＫＦに対して位相変調を実行し、位相変調後
のカーネル関数ＰＭＫＦを用いて磁気共鳴データ（非カーテシアンデータ）ＮＣＤをｋ空
間に配置し、ｋ空間に配置された磁気共鳴データ（再配置済みデータ）ＡＧＤを用いて、
ＭＲ画像Ｉｍｅｒｇｅを生成する。これにより、本実施形態におけるＭＲＩ装置１００お
よび画像生成装置１によれば、第１の実施形態における画像生成処理よりも簡便にＭＲ画
像Ｉｍｅｒｇｅを生成することができる。なお、本実施形態における他の効果は、第１の
実施形態と同様なため、説明は省略する。
第１の実施形態および第２の実施形態における技術的思想を画像生成方法で実現する場
合、当該画像生成方法は、被検体Ｐの撮像に用いられる受信コイル１１７の感度の分布を
示すコイル感度分布Ｓｃと、ｋ空間における非カーテシアンでの撮像により収集された磁
気共鳴データ（非カーテシアンデータ）ＮＣＤとを取得し、コイル感度分布Ｓｃに基づい
て、ｋ空間における磁気共鳴データの配置に関するグリッディングの感度の分布を示すグ
リッディング感度分布Ｓｇとコイル感度分布Ｓｃとの位置合わせを実行し、位置合わせの
結果と磁気共鳴データＮＣＤとコイル感度分布Ｓｃとグリッディング感度分布Ｓｇとに基
づいて、ＭＲ画像Ｉｍｅｒｇｅを生成する。本画像生成方法に関する画生成処理の手順お
よび効果は、第１の実施形態の記載および第２の実施形態の記載と同様なため、説明は省
略する。
合、当該画像生成方法は、被検体Ｐの撮像に用いられる受信コイル１１７の感度の分布を
示すコイル感度分布Ｓｃと、ｋ空間における非カーテシアンでの撮像により収集された磁
気共鳴データ（非カーテシアンデータ）ＮＣＤとを取得し、コイル感度分布Ｓｃに基づい
て、ｋ空間における磁気共鳴データの配置に関するグリッディングの感度の分布を示すグ
リッディング感度分布Ｓｇとコイル感度分布Ｓｃとの位置合わせを実行し、位置合わせの
結果と磁気共鳴データＮＣＤとコイル感度分布Ｓｃとグリッディング感度分布Ｓｇとに基
づいて、ＭＲ画像Ｉｍｅｒｇｅを生成する。本画像生成方法に関する画生成処理の手順お
よび効果は、第１の実施形態の記載および第２の実施形態の記載と同様なため、説明は省
略する。
第１の実施形態および第２の実施形態における技術的思想を画像生成プログラムで実現
する場合、当該画像生成プログラムは、コンピュータに、被検体Ｐの撮像に用いられる受
信コイル１１７の感度の分布を示すコイル感度分布Ｓｃと、ｋ空間における非カーテシア
ンでの撮像により収集された磁気共鳴データ（非カーテシアンデータ）ＮＣＤとを取得し
、コイル感度分布Ｓｃに基づいて、ｋ空間における前記磁気共鳴データの配置に関するグ
リッディングの感度の分布を示すグリッディング感度分布Ｓｇとコイル感度分布Ｓｃとの
位置合わせを実行し、位置合わせの結果と磁気共鳴データＮＣＤとコイル感度分布Ｓｃと
グリッディング感度分布Ｓｇとに基づいて、ＭＲ画像Ｉｍｅｒｇｅを生成すること、を実
現させる。
する場合、当該画像生成プログラムは、コンピュータに、被検体Ｐの撮像に用いられる受
信コイル１１７の感度の分布を示すコイル感度分布Ｓｃと、ｋ空間における非カーテシア
ンでの撮像により収集された磁気共鳴データ（非カーテシアンデータ）ＮＣＤとを取得し
、コイル感度分布Ｓｃに基づいて、ｋ空間における前記磁気共鳴データの配置に関するグ
リッディングの感度の分布を示すグリッディング感度分布Ｓｇとコイル感度分布Ｓｃとの
位置合わせを実行し、位置合わせの結果と磁気共鳴データＮＣＤとコイル感度分布Ｓｃと
グリッディング感度分布Ｓｇとに基づいて、ＭＲ画像Ｉｍｅｒｇｅを生成すること、を実
現させる。
例えば、ＭＲＩ装置１００などのモダリティやＰＡＣＳサーバなどにおけるコンピュー
タに画像生成プログラムをインストールし、これらをメモリ上で展開することによっても
、再構成処理を実現することができる。このとき、コンピュータに当該手法を実行させる
ことのできるプログラムは、磁気ディスク（ハードディスクなど）、光ディスク（ＣＤ－
ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記憶媒体に格納して頒布することも可能であ
る。画像生成プログラムにより画像生成処理の手順および効果は、第１の実施形態および
第２の実施形態と同様なため、説明は省略する。
タに画像生成プログラムをインストールし、これらをメモリ上で展開することによっても
、再構成処理を実現することができる。このとき、コンピュータに当該手法を実行させる
ことのできるプログラムは、磁気ディスク（ハードディスクなど）、光ディスク（ＣＤ－
ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記憶媒体に格納して頒布することも可能であ
る。画像生成プログラムにより画像生成処理の手順および効果は、第１の実施形態および
第２の実施形態と同様なため、説明は省略する。
以上説明した少なくとも１つの実施形態等によれば、アーチファクトの発生を抑制し、
Ｓ／Ｎを向上させた画像を生成することができる。
Ｓ／Ｎを向上させた画像を生成することができる。
いくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり
、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態
で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え
、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、
発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範
囲に含まれるものである。
、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態
で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え
、変更、実施形態同士の組み合わせを行うことができる。これら実施形態やその変形は、
発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範
囲に含まれるものである。
１ 画像生成装置
１１ 通信インターフェース
１３ メモリ
１５ 処理回路
１７ 入出力インターフェース
１００ 磁気共鳴イメージング装置
１０１ 静磁場磁石
１０３ 傾斜磁場コイル
１０５ 傾斜磁場電源
１０７ 寝台
１０９ 寝台制御回路
１１１ ボア
１１３ 送信回路
１１５ 送信コイル
１１７ 受信コイル
１１９ 受信回路
１２１ 撮像制御回路
１２３ システム制御回路
１２５ 記憶装置
１５１ 取得機能
１５３ 調整機能
１５５ 位置合わせ機能
１５７ 生成機能
１１ 通信インターフェース
１３ メモリ
１５ 処理回路
１７ 入出力インターフェース
１００ 磁気共鳴イメージング装置
１０１ 静磁場磁石
１０３ 傾斜磁場コイル
１０５ 傾斜磁場電源
１０７ 寝台
１０９ 寝台制御回路
１１１ ボア
１１３ 送信回路
１１５ 送信コイル
１１７ 受信コイル
１１９ 受信回路
１２１ 撮像制御回路
１２３ システム制御回路
１２５ 記憶装置
１５１ 取得機能
１５３ 調整機能
１５５ 位置合わせ機能
１５７ 生成機能
Claims (12)
- 被検体の撮像に用いられる受信コイルの感度の分布を示すコイル感度分布と、ｋ空間に
おける非カーテシアンでの前記撮像により収集された磁気共鳴データとを取得する取得部
と、
前記コイル感度分布に基づいて、前記ｋ空間における前記磁気共鳴データの配置に関す
るグリッディングの感度の分布を示すグリッディング感度分布と前記コイル感度分布との
位置合わせを実行する位置合わせ部と、
前記位置合わせの結果と前記磁気共鳴データと前記コイル感度分布と前記グリッディン
グ感度分布とに基づいて、磁気共鳴画像を生成する生成部と、
を備えた画像生成装置。 - 前記位置合わせ部は、前記位置合わせの結果として、前記グリッディング感度分布にお
ける感度の重心と前記コイル感度分布における感度の重心との差を算出し、
前記生成部は、
前記差を用いて、前記磁気共鳴データに対して位相変調を実行し、
前記グリッディングの感度の元となるカーネル関数を用いて、前記位相変調後の磁気
共鳴データを前記ｋ空間に配置し、
前記ｋ空間に配置された磁気共鳴データを用いて、前記磁気共鳴画像を生成する、
請求項１に記載の画像生成装置。 - 前記位置合わせ部は、前記位置合わせの結果として、前記グリッディング感度分布にお
ける感度の重心と前記コイル感度分布における感度の重心との差を算出し、
前記生成部は、
前記差を用いて、前記グリッディング感度分布の元となるカーネル関数に対して位相
変調を実行し、
前記位相変調後のカーネル関数を用いて、前記磁気共鳴データを前記ｋ空間に配置し
、
前記ｋ空間に配置された磁気共鳴データを用いて、前記磁気共鳴画像を生成する、
請求項１に記載の画像生成装置。 - 前記生成部は、
前記ｋ空間に配置された磁気共鳴データのフーリエ変換に対して、前記グリッディン
グ感度分布における感度の重心と前記コイル感度分布における感度の重心との差を用いて
前記磁気共鳴データに対して位相変調を実行し、前記差に基づいて前記位相変調とは逆の
位相に対応するように平行移動を実行することで、変換画像を生成し、
前記変換画像を用いて、前記磁気共鳴画像を生成する、
請求項１乃至３のうちいずれか一項に記載の画像生成装置。 - 前記生成部は、
前記グリッディング感度分布における感度の重心と前記コイル感度分布における感度
の重心との差を用いて前記磁気共鳴データに対して位相変調を実行し、前記差に基づいて
、前記ｋ空間に配置された磁気共鳴データを、前記位相変調とは逆の位相に変調し、
前記逆の位相に変調された磁気共鳴データに対してフーリエ変換を行うことにより、
変換画像を生成し、
前記変換画像を用いて、前記磁気共鳴画像を生成する、
請求項１乃至３のうちいずれか一項に記載の画像生成装置。 - 前記受信コイルは、複数のコイルエレメントを有し、
前記取得部は、前記複数のコイルエレメントごとに、前記コイル感度分布と前記磁気共
鳴データとを取得し、
前記位置合わせ部は、前記複数のコイルエレメントごとに、前記コイル感度分布に基づ
いて、前記コイル感度分布と前記グリッディング感度分布との位置合わせを実行し、
前記生成部は、前記複数の位置合わせの結果と前記複数の磁気共鳴データと前記複数の
コイル感度分布と前記グリッディング感度分布とに基づいて、前記磁気共鳴画像を生成す
る、
請求項１乃至５のうちいずれか一項に記載の画像生成装置。 - 前記生成部は、
前記ｋ空間に配置された磁気共鳴データのフーリエ変換に対して、前記グリッディン
グ感度分布における感度の重心と前記コイル感度分布における感度の重心との差を用いて
前記磁気共鳴データに対して位相変調を実行し、前記差に基づいて前記位相変調とは逆の
位相に対応するように平行移動を実行することで、変換画像を生成し、または、前記差に
基づいて、前記ｋ空間に配置された磁気共鳴データを、前記位相変調とは逆の位相に変調
して、前記逆の位相に変調された磁気共鳴データに対してフーリエ変換を行うことにより
、変換画像を生成し、
前記位置合わせの結果を用いて前記グリッディング感度分布と前記コイル感度分布と
を合成することにより、合成感度分布を生成し、
前記変換画像と前記合成感度分布とを用いた感度エンコーディングにより、前記磁気
共鳴画像を生成する、
請求項１乃至６のうちいずれか一項に記載の画像生成装置。 - 前記生成部は、
前記ｋ空間に配置された磁気共鳴データのフーリエ変換に対して、前記グリッディン
グ感度分布における感度の重心と前記コイル感度分布における感度の重心との差を用いて
前記磁気共鳴データに対して位相変調を実行し、前記差に基づいて前記位相変調とは逆の
位相に対応するように平行移動を実行することで、変換画像を生成し、または、前記差に
基づいて、前記ｋ空間に配置された磁気共鳴データを、前記位相変調とは逆の位相に変調
して、前記逆の位相に変調された磁気共鳴データに対してフーリエ変換を行うことにより
、変換画像を生成し、
前記変換画像の平方和の平方根により、前記磁気共鳴画像を生成する、
請求項１乃至６のうちいずれか一項に記載の画像生成装置。 - 前記生成部は、
前記ｋ空間に配置された磁気共鳴データのフーリエ変換に対して、前記グリッディン
グ感度分布における感度の重心と前記コイル感度分布における感度の重心との差を用いて
前記磁気共鳴データに対して位相変調を実行し、前記差に基づいて前記位相変調とは逆の
位相に対応するように平行移動を実行することで、変換画像を生成し、または、前記差に
基づいて、前記ｋ空間に配置された磁気共鳴データを、前記位相変調とは逆の位相に変調
して、前記逆の位相に変調された磁気共鳴データに対してフーリエ変換を行うことにより
、変換画像を生成し、
前記変換画像を前記グリッディング感度分布で除算することにより、前記変換画像に
おけるシェーディングを低減したディシェーディング画像を生成し、
前記ディシェーディング画像と前記コイル感度分布とを用いた感度エンコーディング
、または前記ディシェーディング画像の平方和の平方根により、前記磁気共鳴画像を生成
する、
請求項１乃至６のうちいずれか一項に記載の画像生成装置。 - 前記コイル感度分布に基づいて、前記グリッディング感度分布の元となるカーネル関数
の包絡形状または半値幅を調整する調整部をさらに具備し、
前記生成部は、
前記位置合わせの結果と前記磁気共鳴データと前記コイル感度分布と前記調整された
カーネル関数に対応するグリッディング感度分布とに基づいて、前記磁気共鳴画像を生成
する、
請求項１乃至９のうちいずれか一項に記載の画像生成装置。 - 被検体の撮像に用いられる受信コイルの感度の分布を示すコイル感度分布と、ｋ空間に
おける非カーテシアンでの前記撮像により収集された磁気共鳴データとを取得し、
前記コイル感度分布に基づいて、前記ｋ空間における前記磁気共鳴データの配置に関す
るグリッディングの感度の分布を示すグリッディング感度分布と前記コイル感度分布との
位置合わせを実行し、
前記位置合わせの結果と前記磁気共鳴データと前記コイル感度分布とグリッディング感
度分布とに基づいて、磁気共鳴画像を生成すること、
を備えた画像生成方法。 - コンピュータに、
被検体の撮像に用いられる受信コイルの感度の分布を示すコイル感度分布と、ｋ空間に
おける非カーテシアンでの前記撮像により収集された磁気共鳴データとを取得し、
前記コイル感度分布に基づいて、前記ｋ空間における前記磁気共鳴データの配置に関す
るグリッディングの感度の分布を示すグリッディング感度分布と前記コイル感度分布との
位置合わせを実行し、
前記位置合わせの結果と前記磁気共鳴データと前記コイル感度分布とグリッディング感
度分布とに基づいて、磁気共鳴画像を生成すること、
を実現させる画像生成プログラム。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020157877A JP7451366B2 (ja) | 2020-09-18 | 2020-09-18 | 画像生成装置、画像生成方法、および画像生成プログラム |
| US17/391,758 US11543480B2 (en) | 2020-09-18 | 2021-08-02 | Image generating apparatus, image generating method, and non-volatile computer-readable storage medium storing therein image generating program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020157877A JP7451366B2 (ja) | 2020-09-18 | 2020-09-18 | 画像生成装置、画像生成方法、および画像生成プログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2022051417A JP2022051417A (ja) | 2022-03-31 |
| JP7451366B2 true JP7451366B2 (ja) | 2024-03-18 |
Family
ID=80741497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020157877A Active JP7451366B2 (ja) | 2020-09-18 | 2020-09-18 | 画像生成装置、画像生成方法、および画像生成プログラム |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US11543480B2 (ja) |
| JP (1) | JP7451366B2 (ja) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002056767A1 (fr) | 2001-01-19 | 2002-07-25 | Kabushiki Kaisha Toshiba | Imagerie parallele par resonance magnetique utilisant une carte de sensibilite de bobine haute precision |
| US20060186941A1 (en) | 2003-07-03 | 2006-08-24 | Eth Zurich, Eth Transfer | Minimum norm sence reconstruction for optimal spatial response |
| WO2008111416A1 (ja) | 2007-03-09 | 2008-09-18 | Hitachi Medical Corporation | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
| JP2014503290A (ja) | 2010-12-22 | 2014-02-13 | コーニンクレッカ フィリップス エヌ ヴェ | 校正スキャン、コイル感度マップ及びナビゲータを使用する、剛体動き補償のためのパラレルｍｒｉ方法 |
| US20160104279A1 (en) | 2014-10-10 | 2016-04-14 | Cedars-Sinai Medical Center | Method and system for "push-button" comprehensive cardiac mr examination using continuous self-gated 3d radial imaging |
| JP2016123853A (ja) | 2014-12-29 | 2016-07-11 | 株式会社東芝 | 磁気共鳴イメージング装置及び画像処理装置 |
| JP2019162369A (ja) | 2018-03-20 | 2019-09-26 | 株式会社日立製作所 | 磁気共鳴撮影装置、感度分布算出方法、及び、感度分布算出プログラム |
| JP2019208990A (ja) | 2018-06-07 | 2019-12-12 | キヤノンメディカルシステムズ株式会社 | 医用画像診断装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080154115A1 (en) | 2004-02-10 | 2008-06-26 | Koninklijke Philips Electronics N.V. | Magnetic Resonance Imaging Method |
| US10634752B2 (en) * | 2015-09-14 | 2020-04-28 | Cedars-Sinai Medical Center | Motion corrected simultaneously acquired multiple contrast coronary MRI systems and methods |
| US10802096B2 (en) * | 2017-12-26 | 2020-10-13 | Uih America, Inc. | Methods and systems for magnetic resonance imaging |
-
2020
- 2020-09-18 JP JP2020157877A patent/JP7451366B2/ja active Active
-
2021
- 2021-08-02 US US17/391,758 patent/US11543480B2/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002056767A1 (fr) | 2001-01-19 | 2002-07-25 | Kabushiki Kaisha Toshiba | Imagerie parallele par resonance magnetique utilisant une carte de sensibilite de bobine haute precision |
| US20060186941A1 (en) | 2003-07-03 | 2006-08-24 | Eth Zurich, Eth Transfer | Minimum norm sence reconstruction for optimal spatial response |
| WO2008111416A1 (ja) | 2007-03-09 | 2008-09-18 | Hitachi Medical Corporation | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
| JP2014503290A (ja) | 2010-12-22 | 2014-02-13 | コーニンクレッカ フィリップス エヌ ヴェ | 校正スキャン、コイル感度マップ及びナビゲータを使用する、剛体動き補償のためのパラレルｍｒｉ方法 |
| US20160104279A1 (en) | 2014-10-10 | 2016-04-14 | Cedars-Sinai Medical Center | Method and system for "push-button" comprehensive cardiac mr examination using continuous self-gated 3d radial imaging |
| JP2016123853A (ja) | 2014-12-29 | 2016-07-11 | 株式会社東芝 | 磁気共鳴イメージング装置及び画像処理装置 |
| JP2019162369A (ja) | 2018-03-20 | 2019-09-26 | 株式会社日立製作所 | 磁気共鳴撮影装置、感度分布算出方法、及び、感度分布算出プログラム |
| JP2019208990A (ja) | 2018-06-07 | 2019-12-12 | キヤノンメディカルシステムズ株式会社 | 医用画像診断装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2022051417A (ja) | 2022-03-31 |
| US11543480B2 (en) | 2023-01-03 |
| US20220091212A1 (en) | 2022-03-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10215822B2 (en) | Image processing apparatus and magnetic-resonance imaging apparatus | |
| JP6073627B2 (ja) | 磁気共鳴イメージング装置及び画像処理装置 | |
| US20140126796A1 (en) | Mri reconstruction with motion-dependent regularization | |
| JP6753798B2 (ja) | 医用撮像装置、画像処理方法及びプログラム | |
| US11327134B2 (en) | MRI apparatus, image processing apparatus, and image processing method | |
| JP6548531B2 (ja) | 磁気共鳴イメージング装置及び画像処理装置 | |
| US11016154B2 (en) | Magnetic resonance imaging apparatus and image reconstruction program | |
| JP7451366B2 (ja) | 画像生成装置、画像生成方法、および画像生成プログラム | |
| US11896362B2 (en) | Medical image processing apparatus, method of medical image processing, and nonvolatile computer readable storage medium storing therein medical image processing program | |
| JP6996865B2 (ja) | 磁気共鳴イメージング装置 | |
| JP6618786B2 (ja) | 磁気共鳴イメージング装置及び画像処理装置 | |
| US11918335B2 (en) | Magnetic resonance imaging method, apparatus, and computer storage medium | |
| US20230221393A1 (en) | Image processing apparatus, image processing method, and non-volatile computer-readable storage medium storing therein image processing program | |
| US10459051B2 (en) | Field of view increase in magnetic resonance imaging using nonlinear gradients and generalized iterative reconstruction | |
| JP6809870B2 (ja) | 磁気共鳴イメージング装置、および、信号処理方法 | |
| US11693076B2 (en) | Information processing apparatus, information processing method, and non-transitory computer-readable storage medium storing information processing program | |
| US11896361B2 (en) | Image generating apparatus, image generation method, and nontransitory computer-readable storage medium storing therein image generation program | |
| JP7177715B2 (ja) | 画像処理装置及び画像処理方法 | |
| JP7223619B2 (ja) | 磁気共鳴イメージング装置、画像処理装置及び画像処理方法 | |
| KR102062766B1 (ko) | 물 지방 분리 영상을 획득하는 방법 및 그 자기 공명 영상 장치 | |
| US20240138700A1 (en) | Medical image processing apparatus, method of medical image processing, and nonvolatile computer readable storage medium storing therein medical image processing program | |
| US20240078636A1 (en) | Image processing apparatus, image processing method, and magnetic resonance imaging apparatus | |
| JP2022092151A (ja) | 画像生成装置、画像生成方法、および画像生成プログラム | |
| JP2023180392A (ja) | 磁気共鳴イメージング装置、位相補正方法、および撮像制御方法 | |
| WO2024006944A1 (en) | Methods and systems for super-resolution with progressive sub-voxel up-sampling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230802 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20231221 |
|
| 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: 20240206 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240306 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7451366 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |