半径方向リプル・アーティファクトを低減するＭＲ撮像の方法及びシステム
【課題】磁気共鳴撮像において、ＰＲＯＰＥＬＬＥＲ法又は類似の撮像プロトコルの応用を、頭部アキシャル走査以外のスライス配向及び他の解剖学的部位まで拡張することを可能にする。
【解決手段】スライス配向又は関心のある解剖学的構造とは独立にＰＲＯＰＥＬＬＥＲ撮像を現実的に実行することを可能にするＭＲ撮像の方法及びシステムを開示する。本発明は、ｋ空間の中央領域（８８、９１）の周りに回転されるＭＲデータのブレード（８２、９０）の取得、及びＭＲデータのブレード（８２、９０）から、スライス平面の視野の外部の組織からの位相エイリアシングに起因する半径方向リプル・アーティファクトを実質的に含まない任意のスライス配向の画像を再構成すること（９４）に関する。
【解決手段】スライス配向又は関心のある解剖学的構造とは独立にＰＲＯＰＥＬＬＥＲ撮像を現実的に実行することを可能にするＭＲ撮像の方法及びシステムを開示する。本発明は、ｋ空間の中央領域（８８、９１）の周りに回転されるＭＲデータのブレード（８２、９０）の取得、及びＭＲデータのブレード（８２、９０）から、スライス平面の視野の外部の組織からの位相エイリアシングに起因する半径方向リプル・アーティファクトを実質的に含まない任意のスライス配向の画像を再構成すること（９４）に関する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は一般的には、診断撮像に関し、さらに具体的には、磁気共鳴（ＭＲ）撮像において位相エイリアシング（半径方向リプル）アーティファクトを低減する方法及び装置に関する。本発明はまた、ｋ空間の部分的に充填されたブレードを充填する方法に関する。
【背景技術】
【０００２】
人体組織のような物質が一様な磁場（分極磁場Ｂ０）を受けると、組織内のスピンの個々の磁気モーメントはこの分極磁場と整列しようとするが、スピン固有のラーモア周波数で乱雑な秩序で磁場の周りで歳差運動する。物質すなわち組織が、ｘｙ平面に位置しておりラーモア周波数に近い周波数の磁場（励起磁場Ｂ１）を受けると、整列後の正味のモーメントすなわち「縦磁化」ＭＺが、ｘｙ平面に向かって回転し又は「傾斜」して正味の横磁気モーメントＭｔを発生する。励起信号Ｂ１が停止された後に、励起したスピンによって信号が放出され、この信号を受信し処理して画像を形成することができる。
【０００３】
これらの信号を用いて画像を形成するときには、磁場勾配（Ｇｘ、Ｇｙ及びＧｚ）が用いられる。典型的には、被撮像領域は、用いられている特定の局在化方法に従って上述の各勾配を変化させる一連の測定サイクルによって走査される。結果として得られるＮＭＲ受信信号集合をディジタル化し処理して、多くの周知の再構成手法のいずれかを用いて画像を再構成する。
【０００４】
高速スピン・エコー（ＦＳＥ）撮像は、アーティファクトを最小限にしてＭＲＩデータを収集する効率的な方法として広く用いられているＭＲ撮像手法である。しかしながら、ＦＳＥ画像でも、画像取得は通常数分間を要するため随意又は不随意の患者の運動から生ずるゴースト・アーティファクトを生ずる。
【０００５】
ＦＳＥ画像の運動アーティファクトを低減する多くの撮像手法が開発されている。かかるＦＳＥ手法の一つはPeriodically Rotated Overlapping Parallel Lines with Enhanced Reconstruction（周期的に回転する重畳平行線による改良型再構成、ＰＲＯＰＥＬＬＥＲ）撮像と呼ばれ、エコー鎖においてｋ空間の中心を通る矩形の短冊すなわち「ブレード」を測定するようにデータを収集することにより、ＭＲ信号を符号化する。この短冊は、続いてのエコー鎖では原点を中心としてｋ空間内で漸進的に回転され、これにより、所望の分解でにｋ空間の所要の領域の十分な測定を可能にしている。
【発明の開示】
【発明が解決しようとする課題】
【０００６】
周期的に回転する重畳平行線による改良型再構成は、頭部の近アキシャル走査において被検体の不慮の並進運動及び回転運動に関連したアーティファクトを低減させる場合には極めて実効的であることが判明している。ｋ空間の中心は取得中に多数回サンプリングされるため、ＰＲＯＰＥＬＬＥＲは運動に対する感受性が比較的低い。加えて、回転及び移動に対する明示的な補正を用いて運動アーティファクトをさらに低減している。これらのＰＲＯＰＥＬＬＥＲ撮像の利点にも拘わらず、この手法の応用は、主にサジタル走査、コロナル走査及び斜方（アキシャル方向から著しく変化した配向）走査では半径方向リプル・アーティファクトが観察される結果として、頭部アキシャル走査に限られている。これらのアーティファクトは、肝臓及び脊椎のＰＲＯＰＥＬＬＥＲ撮像でも観察される。ＰＲＯＰＥＬＬＥＲの応用を他のスライス配向及び他の解剖学的部位まで拡張するためには、これらの半径方向リプル・アーティファクトを低減しなければならない。
【０００７】
従って、半径方向リプル・アーティファクトを低減してＰＲＯＰＥＬＬＥＲ又は類似の撮像プロトコルを具現化したＭＲ撮像のシステム及び方法を提供することができると望ましい。
【課題を解決するための手段】
【０００８】
本発明は、ＭＲ画像が半径方向リプル・アーティファクトを実質的に含まないようにして、改良型ＰＲＯＰＥＬＬＥＲ取得を用いて取得されたＭＲデータからＭＲ画像を再構成するシステム及び方法での上述の欠点を克服するものである。
【０００９】
スライス配向又は関心のある解剖学的構造とは独立にＰＲＯＰＥＬＬＥＲ撮像を現実的に実行することを可能にするＭＲ撮像の方法及びシステムを開示する。本発明は、ｋ空間の中央領域の周りに回転されるＭＲデータのブレードの取得、及びＭＲデータのブレードから、従来は所与のＦＯＶの外部に解剖学的構造を有する被検体のコロナル又はサジタルＰＲＯＰＥＬＬＥＲ画像に見受けられていた位相エイリアシング・アーティファクトを実質的に含まない任意のスライス配向の画像を再構成することに関する。
【００１０】
従って、本発明の一観点では、ＭＲ撮像装置が、磁石のボア（中孔）の周りに配置されており分極磁場を印加する複数の勾配コイルを含んでいる。ＲＦ送受信器システム及びＲＦスイッチが、パルス・モジュールによって制御されて、ＲＦコイル・アセンブリへ及びＲＦコイル・アセンブリからＲＦ信号を送受してＭＲ画像を取得する。ＭＲＩ装置はまた、ＦＯＶから任意のスライス配向に沿ってＭＲデータのブレードを取得し、隣り合ったブレードがスライス平面において互いから回転されるようにしてＭＲデータのブレードでｋ空間を充填するようにプログラムされているコンピュータを含んでいる。さらにもう一つの観点では、コンピュータはまた、ｋ空間データの各々のブレードを画像空間へ変換し、幾つかの補正を施して、画像空間をｋ空間へ変換して戻し、このとき半径方向リプル・アーティファクトを実質的に含まない画像が再構成されるようにプログラムされている。
【００１１】
もう一つの観点では、本発明は、ＭＲデータ取得の方法を含んでおり、この方法は、ＦＯＶからのＭＲデータ取得のためのサンプリングを行なうようにエコー鎖を励起するステップと、ｋ空間の中央領域を通って延在する複数のｋ空間データ線を有するｋ空間ブレードを部分的に充填するようにエコー鎖をサンプリングするステップとを含んでいる。この方法はさらに、部分的に充填されたｋ空間ブレードの残りを、エコー鎖からサンプリングされたデータから決定されるデータで充填するステップを含んでいる。励起するステップ、サンプリングするステップ、及び充填するステップは、各々のｋ空間ブレードがｋ空間の中央領域を通って延在し、且つ隣り合ったｋ空間ブレードがｋ空間の中央領域の周りに互いから回転されるようにして、複数のｋ空間ブレードがｋ空間において充填されるまで繰り返される。
【００１２】
もう一つの観点では、本発明は、コンピュータ読み取り可能な記憶媒体に記憶されたコンピュータ・プログラムとして具現化され、このコンピュータ・プログラムは、コンピュータによって実行されると、コンピュータに、エンコードされたＦＯＶを超えて延在する部分を有する対象からｋ空間の複数のブレードを取得して、これらのｋ空間のブレードを撮像する、ことを行なわせる命令を有する。コンピュータはさらに、ＦＯＶの位相エンコード方向での中心部分を記述するブレード画像を単離して、ＦＯＶの位相エンコード方向での中心部分の外部に位置する位相エンコード位置からデータを除去するようにブレード画像を切り取る、ことを行なう。命令はまた、コンピュータに、切り取られたブレード画像をｋ空間の切り取られたブレードへ変換して、単離されたブレード画像に対応するｋ空間のブレードをｋ空間の切り取られたブレードで置き換える、ことを行なわせる。
【００１３】
本発明のその他様々な特徴及び利点は、以下の詳細な説明及び図面から明らかとなろう。
【００１４】
図面は、本発明を実施するのに現状で想到される好適な一実施形態を示す。
【発明を実施するための最良の形態】
【００１５】
図１に、本発明を組み入れた好ましいＭＲイメージング・システム１０の主要な構成要素を示す。このシステムの動作は操作コンソール１２から制御され、操作コンソール１２は、キーボード又は他の入力装置１３、制御パネル１４及び表示画面１６を含んでいる。コンソール１２はリンク１８を介して別個のコンピュータ・システム２０と交信し、システム２０は、操作者が表示画面１６での画像の形成及び表示を制御することを可能にする。コンピュータ・システム２０は、バックプレイン２０ａを介して互いに交信する幾つかのモジュールを含んでいる。これらのモジュールには、画像プロセッサ・モジュール２２、ＣＰＵモジュール２４、及び画像データ配列を記憶するフレーム・バッファとして当技術分野で公知のメモリ・モジュール２６が含まれる。コンピュータ・システム２０は、画像データ及びプログラムの記憶のためのディスク記憶装置２８及びテープ・ドライブ３０に結合されており、高速シリアル・リンク３４を介して別個のシステム制御３２と交信する。入力装置１３は、マウス、ジョイスティック、キーボード、トラックボール、指触起動式画面、光ワンド、音声制御、又は任意の類似の若しくは等価の入力装置を含むことができ、これらの入力装置を用いて対話型での幾何学的構成指定を行なうことができる。
【００１６】
システム制御３２は、バックプレイン３２ａによって共に接続された一組のモジュールを含んでいる。これらのモジュールには、ＣＰＵモジュール３６及びパルス発生器モジュール３８が含まれており、パルス発生器モジュール３８はシリアル・リンク４０を介して操作コンソール１２に接続している。リンク４０を介して、システム制御３２は、実行すべき走査シーケンスを指示する操作者からの指令を受け取る。パルス発生器モジュール３８はシステムの各構成要素を動作させて所望の走査シーケンスを実行すると共に、生成されるＲＦパルスのタイミング、強度及び形状、並びにデータ取得ウィンドウのタイミング及び長さを指示するデータを生成する。パルス発生器モジュール３８は、走査中に発生される勾配パルスのタイミング及び形状を指示する一組の勾配増幅器４２に接続している。パルス発生器モジュール３８はまた、患者に取り付けられた電極からの心電図（ＥＣＧ）信号のように患者に接続された多くの様々なセンサからの信号を受け取る生理学的取得制御器４４から患者データを受け取ることができる。最後に、パルス発生器モジュール３８は、患者及び磁石システムの状態に関連する様々なセンサから信号を受け取る走査室インタフェイス回路４６に接続している。また、走査室インタフェイス回路４６を介して、患者配置システム４８は、走査に望ましい位置まで患者を移動させる命令を受け取る。
【００１７】
パルス発生器モジュール３８によって発生された勾配波形は、Ｇｘ増幅器、Ｇｙ増幅器及びＧｚ増幅器を有する勾配増幅器システム４２に印加される。各々の勾配増幅器は、参照番号５０で全体的に示す勾配コイル・アセンブリ内の対応する物理的な勾配コイルを励起して、取得された信号を空間的に符号化（エンコード）するのに用いられる磁場勾配を生成する。勾配コイル・アセンブリ５０は、分極用磁石５４及び全身型ＲＦコイル５６を含む磁石アセンブリ５２の一部を成す。システム制御３２の送受信器モジュール５８がパルスを発生すると、パルスはＲＦ増幅器６０によって増幅されて、送受信（Ｔ／Ｒ）スイッチ６２によってＲＦコイル５６に結合される。患者の体内の励起した核によって放出される結果信号は同じＲＦコイル５６によって感知されて、送受信スイッチ６２を介して前置増幅器６４に結合される。増幅されたＭＲ信号は、送受信器５８の受信器部で復調、濾波及びディジタル化を施される。送受信スイッチ６２はパルス発生器モジュール３８からの信号によって制御されて、送信モード時にはＲＦ増幅器６０をコイル５６に電気的に接続し、受信モード時には前置増幅器６４をコイル５６に接続する。送受信スイッチ６２はまた、送信モード又は受信モードのいずれでも分離型ＲＦコイル（例えば表面コイル）を用いることを可能にする。
【００１８】
ＲＦコイル５６によって捕えられたＭＲ信号は、送受信器モジュール５８によってディジタル化されて、システム制御３２のメモリ・モジュール６６へ転送される。生（未処理）のｋ空間データの配列がメモリ・モジュール６６に取得されると走査が完了する。この生のｋ空間データを構成し直して、再構成したい各々の画像毎に別個のｋ空間データ配列とすることができ、これらの配列の各々をアレイ・プロセッサ６８に入力すると、アレイ・プロセッサ６８はデータをフーリエ変換するように動作して画像データの配列を生成する。この画像データはシリアル・リンク３４を介してコンピュータ・システム２０へ伝送され、ここで、ディスク記憶装置２８のようなメモリに記憶される。操作コンソール１２から受け取った命令に応じて、この画像データをテープ・ドライブ３０のような長期記憶装置に保管してもよいし、或いは画像データを画像プロセッサ２２によってさらに処理して操作コンソール１２へ伝送し、表示器１６に表示してもよい。
【００１９】
ＰＲＯＰＥＬＬＥＲに基づく取得に関して本発明を説明するが、本発明は他の撮像手法にも応用可能であり得るものと思量される。この点で、走査時間を短縮し、ＳＮＲ及び空間分解能を高めて、位相エイリアシングを低減するのに有利な本発明は、ＰＲＯＰＥＬＬＥＲ以外の多くの撮像手法にも応用可能であり得る。また、本発明が、分割ブレード式ＰＲＯＰＥＬＬＥＲ及びＴＵＲＢＯＰＲＯＰのようなＰＲＯＰＥＬＬＥＲの変形にも応用可能であり得ることが認められよう。
【００２０】
図２について説明する。同図には、ＰＲＯＰＥＬＬＥＲパルス・シーケンス７０の一部が示されている。パルス・シーケンス７０は、図示の例では、関心領域から１２本のスピン・エコー７２を取得するように設計されている。これらのスピン・エコーは全て単一の軸、例えばＧｘに対して収集される。各々のスピン・エコー７２は、ＲＦ再収束パルス７４に続いて周波数エンコード・パルス７６の間に取得され、これら一連のパルスは、定常状態の条件で印加されている。このスピン・エコー・データを用いてｋ空間を充填し、これを図３に模式的に示している。パルス・シーケンスはｋ空間が充填されるまで繰り返されるため、ｋ空間の中央領域では過剰サンプリングが生ずる。
【００２１】
図３は、図２に関して説明したＰＲＯＰＥＬＬＥＲパルス・シーケンスの反復印加の結果として取得されたＭＲデータで充填されたｋ空間７８を模式的に示す。ＰＲＯＰＥＬＬＥＲでは、取得される各々のエコーはｋ空間７６の単一の線８０に対応する。このようなものとして、１２本のスピン・エコー・データ取得では、ｋ空間の各々のブレード８２が１２本のデータの線を含む。ＰＲＯＰＥＬＬＥＲに基づく撮像は、各々のエコー鎖でのｋ空間データのブレードの回転を具現化している。この点で、ｋ空間のブレードは、ｋ空間が充填されるまで各々のエコー鎖取得によってｋ空間の中心８８の周りで漸進的に回転される。
【００２２】
図３について説明したように、ｋ空間７８はＰＲＯＰＥＬＬＥＲ取得を用いて充填される。この取得を通じて、ｋ空間の中央領域を通って延在する所与のブレード全ての位相エンコード位置は個々にサンプリングされる。すなわち、ｋ空間の各々の縦列が位相エンコード位置に対応し、各々の横列が周波数エンコード位置に対応する。従って、ｋ空間の各々の点を位相エンコード位置及び周波数エンコード位置によって識別することができる。位相エンコード位置（又はｋ空間内の縦列）の数は、１本のエコーのサンプリング中に印加される位相エンコード勾配の数の関数である。位相エンコード位置の数が増すにつれて、必要な位相エンコード・ステップの数も増す。このようなものとして、位相エンコード位置の数が増加すると、走査時間も必然的に増加する。位相エンコード位置の数が減少すると走査時間は短縮するが、信号及び分解能の損失を伴う。従って、一実施形態では、本発明は、ｋ空間のブレードのホモダイン処理を含んでいる。この処理を図４に示す。説明の目的でｋ空間の単一のブレードについて説明するが、当業者は、ｋ空間の各々の取得されたブレードに対して本書で説明するホモダイン処理手法を応用し得ることを認められよう。
【００２３】
ｋ空間ブレード９０は、ｋ空間の中央領域９１又は原点を通って延在しており、ブレード９０の連続的な回転時に中央領域９１は過剰サンプリングされる。ｋ空間ブレード９０は好ましくは、所定の数のｋ空間データ線９２によって画定され、各々の線９２がサンプリング済み（実線）又は未サンプリング（破線）のいずれかとして特徴付けられる。さらに、図４には、ｋ空間ブレード９０の全てに満たない位相エンコード位置（ビュー）が、サンプリング済みのＭＲデータによって充填されていることが示されている。全てに満たない位相エンコード位置がサンプリングされるので、位相エンコード・ステップの数を減少させることができる。代替的には、各々のｋ空間データ線の部分的充填の後に、位相エンコード・ステップの総数が、ホモダイン処理を行なっていない取得のステップ数と等しくなるようにして、位相エンコード位置（従って位相エンコード・ステップ）の数を増加させることもできる。後に改めて説明するが、未サンプリングの位相エンコード位置は、サンプリング済みの位相エンコード位置に基づくデータで充填される。
【００２４】
多くの部分的フーリエ再構成手法を用いて未サンプリングの位相エンコード位置のデータを決定することができるが、好ましい一手法はホモダイン処理である。ホモダイン処理は、エルミート共役対称性のための一対のフィルタを用いて、サンプリング済みの位相エンコード位置のデータに基づいて未サンプリングの位相エンコード位置のデータを決定する。未サンプリングの位相エンコード位置のデータを提供することに加え、ホモダイン処理はまた、サンプリングされたデータに含まれて共鳴周波数、流れ及び運動の変動にしばしば帰属可能な位相の誤差又はばらつきを補正する。ホモダイン処理又は他の部分的フーリエ再構成手法では、他の場合であれば要求されるような位相エンコード・ステップを必要とせずに、未サンプリングの位相エンコード位置（図４では「Ｕ」として示す）を「実効的に」サンプリングすることができる。加えて、１本のエコーをサンプリングするのに必要な時間（ＴＥ）が実効的に短縮され、これによっても走査時間を短縮して被検体スループットを増大させることができる。
【００２５】
好ましくは、サンプリング済みの位相エンコード位置は、ｋ空間の横列１列の二分の一よりも僅かに多い部分を含む。この点で、ｋ空間データ線の中心の周りの位相エンコード位置について空間周波数データが取得される。このようにして、位相エンコード・ステップの数を二分の一よりも僅かに少ない分だけ減少し得るものと思量される。一方、印加される位相エンコード勾配の数を減少させずに定義し直して、殆ど全てのサンプリング済みの位相エンコード位置がｋ空間データ線の二分の一よりも僅かに多い部分に配置されるようにする。結果として、走査時間は短縮せず、信号対雑音が増大して空間分解能が高まる。事実上、位相エンコード・ステップの数を倍増せずに位相エンコード位置の数を倍増することができる。実際に、本発明によって、位相エンコード位置（位相ＦＯＶ）の数の倍増を上回ることが可能になるものと思量される。
【００２６】
図５には、ＰＲＯＰＥＬＬＥＲ取得によってＭＲデータを取得して、半径方向リプル・アーティファクトを実質的に含まない画像を再構成する手法の各ステップが掲げられている。工程９４（コンピュータで実行可能な言語で具現化されており、１台又は複数のコンピュータによって実行され得る）はステップ９６で開始して、ＰＲＯＰＥＬＬＥＲ取得の設定を行なう。ステップ９６では、ＰＲＯＰＥＬＬＥＲ取得のパラメータを、所定数の利用者入力から読み込むか又は自動的に決定する。これらのパラメータとしては、繰り返し時間（ＴＲ）、視野（ＦＯＶ）の大きさ、及びスライス配向、例えばアキシャル、頭部、コロナル又は斜方等の諸設定がある。上述のように、本発明は特に、コロナル、サジタル又は斜方ＰＲＯＰＥＬＬＥＲ走査で取得されたＭＲデータから再構成される画像での半径方向リプル（位相エイリアシング）アーティファクトを低減することに関する。
【００２７】
ＰＲＯＰＥＬＬＥＲ走査のパラメータを入力し、設定し、また他の場合には決定した後に、ステップ９８においてＭＲデータを取得する。エンコードされて励起されたエコーをサンプリングして用いて、互いから回転されておりｋ空間の中央領域を通って延在するｋ空間のブレードを部分的に充填する。次いで、ブレード・ホモダイン処理ステップ９００を実行して、図４に関して説明したように、サンプリング済みの位相エンコード位置から、未サンプリングの位相エンコード位置のデータを決定する。ブレード・ホモダイン処理１００の後に、ステップ１０２において画像空間位相補正を行なって、画像空間データの位相誤差を補正する。位相誤差補正１０２の後に、ステップ１０４において逆高速フーリエ変換（ＦＦＴ）を用いて、位相を補正した画像空間を下降（ダウン）サンプリングする。ステップ１０４において下降サンプリングが実行されて、ＦＯＶの位相エンコード方向での中心部分のみを記述するｋ空間のブレードを決定する。このステップは、中心部分以外からのデータが画像再構成時に除外されるように、ＦＯＶの位相エンコード方向での中心部分に対応するｋ空間のブレードを切り取ることにより達成される。このように、ブレード画像空間を切り取った後に、得られたブレード画像に逆ＦＦＴを施してｋ空間ブレードを復元する。次いで、このｋ空間ブレードをＰＲＯＰＥＬＬＥＲ再構成工程の残りの部分に通すと、工程が実行されて全てのｋ空間ブレードからの画像を再構成する。
【００２８】
明確に述べると、下降サンプリング１０４の後に、所定数のＰＲＯＰＥＬＬＥＲ再構成ステップを実行して運動を補正する。これらの運動補正ステップは、回転補正１０６、並進補正１０８、相関加重１１０、及び最後に、画像再構成１１２を含んでいる。画像再構成１１２の後に、ステップ１１４において工程は終了する。以上の各ステップの結果として、得られた再構成画像は半径方向リプル（位相エイリアシング）アーティファクトを実質的に含まない。
【００２９】
本発明の一つの利点は、ブレード・ホモダイン・ステップ及び下降サンプリング・ステップを従来のＰＲＯＰＥＬＬＥＲ再構成手法にシームレスに挿入し、しかもＳＮＲ及び空間分解能を高めて、各々の半径方向ＰＲＯＰＥＬＬＥＲブレードに沿った位相エイリアシングを除去し得ることである。このように、典型的には位相ＦＯＶを超えて延在する解剖学的構造から画像に現われる半径方向リプル・アーティファクトが低減する。すなわち、本発明は、ＰＲＯＰＥＬＬＥＲブレード同士の間の位相折り返しの組み合わせの結果として現われる得る位相エイリアシング・アーティファクトを低減する。
【００３０】
図６及び図７では、２枚のブレード画像が、図４に関して上で述べたブレード・ホモダイン処理ステップを示している。図６は、ｋ空間の半部で８枚のビュー及びｋ空間の残り半部で２枚のビューを示す１０枚のビューでのブレード画像である。ホモダイン処理の後に、１０枚のビューを、ｋ空間の半部で８枚のビュー及びｋ空間の残り半部で８枚のビューを示す１６枚のビューでのブレードへ変換する。ブレード・ホモダイン処理の結果として、６枚の追加のデータのビューが画像再構成に用いられ、これにより、データ取得時に６回分の追加で時間を費やす位相エンコード・ステップを行なわないでもＳＮＲ及び空間分解能が高まる。
【００３１】
図８には、従来のＰＲＯＰＥＬＬＥＲ取得を用いてＦＯＶを超えて延在する部分を有する球形ファントムから再構成される画像を示す。図示のように、位相エンコード方向に沿ってエイリアシングが生じている。一方、図９では、同じ球形ファントムを本発明に従って撮像したときに、図８で優勢であった位相エイリアシング・アーティファクト又は半径方向リプル・アーティファクトが除去されている。この点で、対象がＦＯＶを位相エンコード方向に超えて延在するにも拘わらず、位相エンコード方向には位相折り返しが生じない。尚、図８の画像は、ＰＲＯＰＥＬＬＥＲに基づく取得を用いて、エコー鎖長さ（ＥＴＬ）を２４、ｆｕｌｌ ＮＥＸ、及び位相ＦＯＶを２０ｃｍとして撮像されたことを特記しておく。図９の画像は、ＰＲＯＰＥＬＬＥＲを用いて、ＥＴＬ＝２４、ｈａｌｆ ＮＥＸ、及び位相ＦＯＶを４０ｃｍとして取得した。このように、図５に関して説明したホモダイン工程を用いて、ＦＯＶを位相エンコード方向に実効的に倍増させた。
【００３２】
図１０及び図１１はさらに、本発明によって達成される利点を示す。明確に述べると、図１０に示すように、ＦＯＶを位相エンコード方向に超えて延在する対象から取得される全てのブレードの従来のＰＲＯＰＥＬＬＥＲ再構成は、画像にリプルを生ずる。これらのリプルは、位相折り返しエイリアシングの結果である。一方、図１１は図１０で撮像した同じ対象の画像であるが、本書で説明した改良型ＰＲＯＰＥＬＬＥＲ取得を用いており、位相エイリアシング・リプルが存在しない。図１１の画像は、図１０の画像よりも遥かに一様であり、従って診断的価値が高い。
【００３３】
図１２及び図１３はさらに、本発明によって達成される利点を示す。明確に述べると、図１２に示すように、サジタル平面に沿って得られた脊椎画像の従来のＴ２−ＰＲＯＰＥＬＬＥＲ再構成は、図面では脊柱領域の縞として現われている半径方向リプル・アーティファクトを生じ得る。一方、本書で開示している改良型ＰＲＯＰＥＬＬＥＲ取得は、本発明によれば、図１３に示すようにこれらの半径方向リプル・アーティファクトを低減する。明確に述べると、図１２で見られる「縞」が図１３には存在しない。
【００３４】
従って、ＭＲＩ装置が、磁石のボアの周りに配置されており分極磁場を印加する複数の勾配コイルを含んでいる。ＲＦ送受信器システム及びＲＦスイッチが、パルス・モジュールによって制御されて、ＲＦコイル・アセンブリへ及びＲＦコイル・アセンブリからＲＦ信号を送受してＭＲ画像を取得する。ＭＲＩ装置はまた、ＦＯＶからコロナル・スライス配向及びサジタル・スライス配向のいずれかに沿ってＭＲデータのブレードを取得し、隣り合ったブレードが互いから回転されるようにＭＲデータのブレードでｋ空間を充填するようにプログラムされているコンピュータを含んでいる。コンピュータはまた、ｋ空間データを画像空間へ変換して、画像空間から半径方向リプル・アーティファクトを実質的に含まない画像を再構成するようにプログラムされている。
【００３５】
また、ＭＲデータ取得の方法が開示されており、この方法は、ＦＯＶからのＭＲデータ取得のためのサンプリングを行なうようにエコー鎖を励起するステップと、ｋ空間の中央領域を通って延在する複数のｋ空間データ線を有するｋ空間ブレードを部分的に充填するようにエコー鎖をサンプリングするステップとを含んでいる。この方法はさらに、部分的に充填されたｋ空間ブレードの残りを、エコー鎖からサンプリングされたデータから決定されるデータで充填するステップを含んでいる。励起するステップ、サンプリングするステップ、及び充填するステップは、複数のｋ空間ブレードがｋ空間において充填されるまで繰り返されて、各々のｋ空間ブレードがｋ空間の中央領域を通って延在し、且つ隣り合ったｋ空間ブレードがｋ空間の中央領域の周りで互いから回転されるようにｋ空間を充填する。
【００３６】
本発明はまた、コンピュータ読み取り可能な記憶媒体に記憶されたコンピュータ・プログラムとして具現化され、このコンピュータ・プログラムは、コンピュータによって実行されると、コンピュータに、ＦＯＶを超えて延在する部分を有する対象からｋ空間の複数のブレードを取得して、これらのｋ空間のブレードを撮像する、ことを行なわせる命令を有する。コンピュータはさらに、ＦＯＶの位相エンコード方向での中心部分を記述するブレード画像を単離して、ＦＯＶの位相エンコード方向での中心部分の外部に位置する位相エンコード位置からデータを除去するようにブレード画像を切り取る、ことを行なう。命令はまた、コンピュータに、切り取られたブレード画像をｋ空間の切り取られたブレードへ変換して、単離されたブレード画像に対応するｋ空間のブレードをｋ空間の切り取られたブレードで置き換える、ことを行なわせる。
【００３７】
本発明は好適実施形態について説明されており、明示的に記載した以外の均等構成、代替構成及び改変が可能であり特許請求の範囲内にあることが認められよう。また、図面の符号に対応する特許請求の範囲中の符号は、単に本願発明の理解をより容易にするために用いられているものであり、本願発明の範囲を狭める意図で用いられたものではない。そして、本願の特許請求の範囲に記載した事項は、明細書に組み込まれ、明細書の記載事項の一部となる。
【図面の簡単な説明】
【００３８】
【図１】本発明と共に用いられるＭＲイメージング・システムの模式的ブロック図である。
【図２】本発明で応用可能なＰＲＯＰＥＬＬＥＲパルス・シーケンスの例の一部を示す図である。
【図３】図２のＰＲＯＰＥＬＬＥＲパルス・シーケンスの例を用いて取得されたＭＲデータで充填されるｋ空間の模式図である。
【図４】本発明のもう一つの観点に従ってホモダイン処理を用いて充填されるｋ空間の単一のブレードの模式図である。
【図５】本発明のもう一つの観点によるＰＲＯＰＥＬＬＥＲ再構成の各ステップを示す流れ図である。
【図６】ｋ空間の部分的に充填されたブレードを示す画像である。
【図７】元々部分的に充填されており、続いてホモダイン処理を用いて充填されたｋ空間のブレードを示す画像である。
【図８】位相ＦＯＶを超えて延在する部分を有する球形ファントムから従来のＰＲＯＰＥＬＬＥＲ取得を用いて再構成される単一のブレード画像である。
【図９】本発明の一観点に従って位相エイリアシングを低減した図８で撮像された球形ファントムから再構成される単一のブレード画像である。
【図１０】様々な配向で取得されたデータのブレードから従来のＰＲＯＰＥＬＬＥＲ再構成手法を用いて再構成されたファントムの画像である。
【図１１】図１０で撮像されたファントムの画像であるが、本発明の一観点による改良型ＰＲＯＰＥＬＬＥＲ再構成手法によって様々な配向で取得されたデータのブレードからの画像である。
【図１２】従来のＰＲＯＰＥＬＬＥＲ取得を用いてサジタル平面に沿って得られた脊椎画像である。
【図１３】本発明の一観点による改良型ＰＲＯＰＥＬＬＥＲ取得を用いてやはりサジタル平面に沿って得られた脊椎画像であって、図１２の画像に見られる半径方向リプル（位相エイリアシング）アーティファクトの存在しない画像である。
【符号の説明】
【００３９】
１０ ＭＲイメージング・システム
１２ 操作コンソール
１３ 入力装置
１４ 制御パネル
１６ 表示画面
１８ リンク
２０ コンピュータ・システム
２０ａ、３２ａ バックプレイン
２２ 画像プロセッサ・モジュール
２４、３６ ＣＰＵモジュール
２６、６６ メモリ・モジュール
２８ ディスク記憶装置
３０ テープ・ドライブ
３２ システム制御
３４ 高速シリアル・リンク
３８ パルス発生器モジュール
４０ シリアル・リンク
４２ 勾配増幅器
４４ 生理学的取得制御器
４６ 走査室インタフェイス回路
４８ 患者配置システム
５０ 勾配コイル・アセンブリ
５２ 磁石アセンブリ
５４ 分極用磁石
５６ 全身型ＲＦコイル
５８ 送受信器モジュール
６０ ＲＦ増幅器
６２ 送受信（Ｔ／Ｒ）スイッチ
６４ 前置増幅器
６８ アレイ・プロセッサ
７０ ＰＲＯＰＥＬＬＥＲパルス・シーケンス
７２ スピン・エコー
７４ ＲＦ再収束パルス
７６ 周波数エンコード・パルス
７８ ｋ空間
８０ データ線
８８ ｋ空間の中心
９０ ｋ空間ブレード
９１ 中央領域
９２ ｋ空間データ線
【技術分野】
【０００１】
本発明は一般的には、診断撮像に関し、さらに具体的には、磁気共鳴（ＭＲ）撮像において位相エイリアシング（半径方向リプル）アーティファクトを低減する方法及び装置に関する。本発明はまた、ｋ空間の部分的に充填されたブレードを充填する方法に関する。
【背景技術】
【０００２】
人体組織のような物質が一様な磁場（分極磁場Ｂ０）を受けると、組織内のスピンの個々の磁気モーメントはこの分極磁場と整列しようとするが、スピン固有のラーモア周波数で乱雑な秩序で磁場の周りで歳差運動する。物質すなわち組織が、ｘｙ平面に位置しておりラーモア周波数に近い周波数の磁場（励起磁場Ｂ１）を受けると、整列後の正味のモーメントすなわち「縦磁化」ＭＺが、ｘｙ平面に向かって回転し又は「傾斜」して正味の横磁気モーメントＭｔを発生する。励起信号Ｂ１が停止された後に、励起したスピンによって信号が放出され、この信号を受信し処理して画像を形成することができる。
【０００３】
これらの信号を用いて画像を形成するときには、磁場勾配（Ｇｘ、Ｇｙ及びＧｚ）が用いられる。典型的には、被撮像領域は、用いられている特定の局在化方法に従って上述の各勾配を変化させる一連の測定サイクルによって走査される。結果として得られるＮＭＲ受信信号集合をディジタル化し処理して、多くの周知の再構成手法のいずれかを用いて画像を再構成する。
【０００４】
高速スピン・エコー（ＦＳＥ）撮像は、アーティファクトを最小限にしてＭＲＩデータを収集する効率的な方法として広く用いられているＭＲ撮像手法である。しかしながら、ＦＳＥ画像でも、画像取得は通常数分間を要するため随意又は不随意の患者の運動から生ずるゴースト・アーティファクトを生ずる。
【０００５】
ＦＳＥ画像の運動アーティファクトを低減する多くの撮像手法が開発されている。かかるＦＳＥ手法の一つはPeriodically Rotated Overlapping Parallel Lines with Enhanced Reconstruction（周期的に回転する重畳平行線による改良型再構成、ＰＲＯＰＥＬＬＥＲ）撮像と呼ばれ、エコー鎖においてｋ空間の中心を通る矩形の短冊すなわち「ブレード」を測定するようにデータを収集することにより、ＭＲ信号を符号化する。この短冊は、続いてのエコー鎖では原点を中心としてｋ空間内で漸進的に回転され、これにより、所望の分解でにｋ空間の所要の領域の十分な測定を可能にしている。
【発明の開示】
【発明が解決しようとする課題】
【０００６】
周期的に回転する重畳平行線による改良型再構成は、頭部の近アキシャル走査において被検体の不慮の並進運動及び回転運動に関連したアーティファクトを低減させる場合には極めて実効的であることが判明している。ｋ空間の中心は取得中に多数回サンプリングされるため、ＰＲＯＰＥＬＬＥＲは運動に対する感受性が比較的低い。加えて、回転及び移動に対する明示的な補正を用いて運動アーティファクトをさらに低減している。これらのＰＲＯＰＥＬＬＥＲ撮像の利点にも拘わらず、この手法の応用は、主にサジタル走査、コロナル走査及び斜方（アキシャル方向から著しく変化した配向）走査では半径方向リプル・アーティファクトが観察される結果として、頭部アキシャル走査に限られている。これらのアーティファクトは、肝臓及び脊椎のＰＲＯＰＥＬＬＥＲ撮像でも観察される。ＰＲＯＰＥＬＬＥＲの応用を他のスライス配向及び他の解剖学的部位まで拡張するためには、これらの半径方向リプル・アーティファクトを低減しなければならない。
【０００７】
従って、半径方向リプル・アーティファクトを低減してＰＲＯＰＥＬＬＥＲ又は類似の撮像プロトコルを具現化したＭＲ撮像のシステム及び方法を提供することができると望ましい。
【課題を解決するための手段】
【０００８】
本発明は、ＭＲ画像が半径方向リプル・アーティファクトを実質的に含まないようにして、改良型ＰＲＯＰＥＬＬＥＲ取得を用いて取得されたＭＲデータからＭＲ画像を再構成するシステム及び方法での上述の欠点を克服するものである。
【０００９】
スライス配向又は関心のある解剖学的構造とは独立にＰＲＯＰＥＬＬＥＲ撮像を現実的に実行することを可能にするＭＲ撮像の方法及びシステムを開示する。本発明は、ｋ空間の中央領域の周りに回転されるＭＲデータのブレードの取得、及びＭＲデータのブレードから、従来は所与のＦＯＶの外部に解剖学的構造を有する被検体のコロナル又はサジタルＰＲＯＰＥＬＬＥＲ画像に見受けられていた位相エイリアシング・アーティファクトを実質的に含まない任意のスライス配向の画像を再構成することに関する。
【００１０】
従って、本発明の一観点では、ＭＲ撮像装置が、磁石のボア（中孔）の周りに配置されており分極磁場を印加する複数の勾配コイルを含んでいる。ＲＦ送受信器システム及びＲＦスイッチが、パルス・モジュールによって制御されて、ＲＦコイル・アセンブリへ及びＲＦコイル・アセンブリからＲＦ信号を送受してＭＲ画像を取得する。ＭＲＩ装置はまた、ＦＯＶから任意のスライス配向に沿ってＭＲデータのブレードを取得し、隣り合ったブレードがスライス平面において互いから回転されるようにしてＭＲデータのブレードでｋ空間を充填するようにプログラムされているコンピュータを含んでいる。さらにもう一つの観点では、コンピュータはまた、ｋ空間データの各々のブレードを画像空間へ変換し、幾つかの補正を施して、画像空間をｋ空間へ変換して戻し、このとき半径方向リプル・アーティファクトを実質的に含まない画像が再構成されるようにプログラムされている。
【００１１】
もう一つの観点では、本発明は、ＭＲデータ取得の方法を含んでおり、この方法は、ＦＯＶからのＭＲデータ取得のためのサンプリングを行なうようにエコー鎖を励起するステップと、ｋ空間の中央領域を通って延在する複数のｋ空間データ線を有するｋ空間ブレードを部分的に充填するようにエコー鎖をサンプリングするステップとを含んでいる。この方法はさらに、部分的に充填されたｋ空間ブレードの残りを、エコー鎖からサンプリングされたデータから決定されるデータで充填するステップを含んでいる。励起するステップ、サンプリングするステップ、及び充填するステップは、各々のｋ空間ブレードがｋ空間の中央領域を通って延在し、且つ隣り合ったｋ空間ブレードがｋ空間の中央領域の周りに互いから回転されるようにして、複数のｋ空間ブレードがｋ空間において充填されるまで繰り返される。
【００１２】
もう一つの観点では、本発明は、コンピュータ読み取り可能な記憶媒体に記憶されたコンピュータ・プログラムとして具現化され、このコンピュータ・プログラムは、コンピュータによって実行されると、コンピュータに、エンコードされたＦＯＶを超えて延在する部分を有する対象からｋ空間の複数のブレードを取得して、これらのｋ空間のブレードを撮像する、ことを行なわせる命令を有する。コンピュータはさらに、ＦＯＶの位相エンコード方向での中心部分を記述するブレード画像を単離して、ＦＯＶの位相エンコード方向での中心部分の外部に位置する位相エンコード位置からデータを除去するようにブレード画像を切り取る、ことを行なう。命令はまた、コンピュータに、切り取られたブレード画像をｋ空間の切り取られたブレードへ変換して、単離されたブレード画像に対応するｋ空間のブレードをｋ空間の切り取られたブレードで置き換える、ことを行なわせる。
【００１３】
本発明のその他様々な特徴及び利点は、以下の詳細な説明及び図面から明らかとなろう。
【００１４】
図面は、本発明を実施するのに現状で想到される好適な一実施形態を示す。
【発明を実施するための最良の形態】
【００１５】
図１に、本発明を組み入れた好ましいＭＲイメージング・システム１０の主要な構成要素を示す。このシステムの動作は操作コンソール１２から制御され、操作コンソール１２は、キーボード又は他の入力装置１３、制御パネル１４及び表示画面１６を含んでいる。コンソール１２はリンク１８を介して別個のコンピュータ・システム２０と交信し、システム２０は、操作者が表示画面１６での画像の形成及び表示を制御することを可能にする。コンピュータ・システム２０は、バックプレイン２０ａを介して互いに交信する幾つかのモジュールを含んでいる。これらのモジュールには、画像プロセッサ・モジュール２２、ＣＰＵモジュール２４、及び画像データ配列を記憶するフレーム・バッファとして当技術分野で公知のメモリ・モジュール２６が含まれる。コンピュータ・システム２０は、画像データ及びプログラムの記憶のためのディスク記憶装置２８及びテープ・ドライブ３０に結合されており、高速シリアル・リンク３４を介して別個のシステム制御３２と交信する。入力装置１３は、マウス、ジョイスティック、キーボード、トラックボール、指触起動式画面、光ワンド、音声制御、又は任意の類似の若しくは等価の入力装置を含むことができ、これらの入力装置を用いて対話型での幾何学的構成指定を行なうことができる。
【００１６】
システム制御３２は、バックプレイン３２ａによって共に接続された一組のモジュールを含んでいる。これらのモジュールには、ＣＰＵモジュール３６及びパルス発生器モジュール３８が含まれており、パルス発生器モジュール３８はシリアル・リンク４０を介して操作コンソール１２に接続している。リンク４０を介して、システム制御３２は、実行すべき走査シーケンスを指示する操作者からの指令を受け取る。パルス発生器モジュール３８はシステムの各構成要素を動作させて所望の走査シーケンスを実行すると共に、生成されるＲＦパルスのタイミング、強度及び形状、並びにデータ取得ウィンドウのタイミング及び長さを指示するデータを生成する。パルス発生器モジュール３８は、走査中に発生される勾配パルスのタイミング及び形状を指示する一組の勾配増幅器４２に接続している。パルス発生器モジュール３８はまた、患者に取り付けられた電極からの心電図（ＥＣＧ）信号のように患者に接続された多くの様々なセンサからの信号を受け取る生理学的取得制御器４４から患者データを受け取ることができる。最後に、パルス発生器モジュール３８は、患者及び磁石システムの状態に関連する様々なセンサから信号を受け取る走査室インタフェイス回路４６に接続している。また、走査室インタフェイス回路４６を介して、患者配置システム４８は、走査に望ましい位置まで患者を移動させる命令を受け取る。
【００１７】
パルス発生器モジュール３８によって発生された勾配波形は、Ｇｘ増幅器、Ｇｙ増幅器及びＧｚ増幅器を有する勾配増幅器システム４２に印加される。各々の勾配増幅器は、参照番号５０で全体的に示す勾配コイル・アセンブリ内の対応する物理的な勾配コイルを励起して、取得された信号を空間的に符号化（エンコード）するのに用いられる磁場勾配を生成する。勾配コイル・アセンブリ５０は、分極用磁石５４及び全身型ＲＦコイル５６を含む磁石アセンブリ５２の一部を成す。システム制御３２の送受信器モジュール５８がパルスを発生すると、パルスはＲＦ増幅器６０によって増幅されて、送受信（Ｔ／Ｒ）スイッチ６２によってＲＦコイル５６に結合される。患者の体内の励起した核によって放出される結果信号は同じＲＦコイル５６によって感知されて、送受信スイッチ６２を介して前置増幅器６４に結合される。増幅されたＭＲ信号は、送受信器５８の受信器部で復調、濾波及びディジタル化を施される。送受信スイッチ６２はパルス発生器モジュール３８からの信号によって制御されて、送信モード時にはＲＦ増幅器６０をコイル５６に電気的に接続し、受信モード時には前置増幅器６４をコイル５６に接続する。送受信スイッチ６２はまた、送信モード又は受信モードのいずれでも分離型ＲＦコイル（例えば表面コイル）を用いることを可能にする。
【００１８】
ＲＦコイル５６によって捕えられたＭＲ信号は、送受信器モジュール５８によってディジタル化されて、システム制御３２のメモリ・モジュール６６へ転送される。生（未処理）のｋ空間データの配列がメモリ・モジュール６６に取得されると走査が完了する。この生のｋ空間データを構成し直して、再構成したい各々の画像毎に別個のｋ空間データ配列とすることができ、これらの配列の各々をアレイ・プロセッサ６８に入力すると、アレイ・プロセッサ６８はデータをフーリエ変換するように動作して画像データの配列を生成する。この画像データはシリアル・リンク３４を介してコンピュータ・システム２０へ伝送され、ここで、ディスク記憶装置２８のようなメモリに記憶される。操作コンソール１２から受け取った命令に応じて、この画像データをテープ・ドライブ３０のような長期記憶装置に保管してもよいし、或いは画像データを画像プロセッサ２２によってさらに処理して操作コンソール１２へ伝送し、表示器１６に表示してもよい。
【００１９】
ＰＲＯＰＥＬＬＥＲに基づく取得に関して本発明を説明するが、本発明は他の撮像手法にも応用可能であり得るものと思量される。この点で、走査時間を短縮し、ＳＮＲ及び空間分解能を高めて、位相エイリアシングを低減するのに有利な本発明は、ＰＲＯＰＥＬＬＥＲ以外の多くの撮像手法にも応用可能であり得る。また、本発明が、分割ブレード式ＰＲＯＰＥＬＬＥＲ及びＴＵＲＢＯＰＲＯＰのようなＰＲＯＰＥＬＬＥＲの変形にも応用可能であり得ることが認められよう。
【００２０】
図２について説明する。同図には、ＰＲＯＰＥＬＬＥＲパルス・シーケンス７０の一部が示されている。パルス・シーケンス７０は、図示の例では、関心領域から１２本のスピン・エコー７２を取得するように設計されている。これらのスピン・エコーは全て単一の軸、例えばＧｘに対して収集される。各々のスピン・エコー７２は、ＲＦ再収束パルス７４に続いて周波数エンコード・パルス７６の間に取得され、これら一連のパルスは、定常状態の条件で印加されている。このスピン・エコー・データを用いてｋ空間を充填し、これを図３に模式的に示している。パルス・シーケンスはｋ空間が充填されるまで繰り返されるため、ｋ空間の中央領域では過剰サンプリングが生ずる。
【００２１】
図３は、図２に関して説明したＰＲＯＰＥＬＬＥＲパルス・シーケンスの反復印加の結果として取得されたＭＲデータで充填されたｋ空間７８を模式的に示す。ＰＲＯＰＥＬＬＥＲでは、取得される各々のエコーはｋ空間７６の単一の線８０に対応する。このようなものとして、１２本のスピン・エコー・データ取得では、ｋ空間の各々のブレード８２が１２本のデータの線を含む。ＰＲＯＰＥＬＬＥＲに基づく撮像は、各々のエコー鎖でのｋ空間データのブレードの回転を具現化している。この点で、ｋ空間のブレードは、ｋ空間が充填されるまで各々のエコー鎖取得によってｋ空間の中心８８の周りで漸進的に回転される。
【００２２】
図３について説明したように、ｋ空間７８はＰＲＯＰＥＬＬＥＲ取得を用いて充填される。この取得を通じて、ｋ空間の中央領域を通って延在する所与のブレード全ての位相エンコード位置は個々にサンプリングされる。すなわち、ｋ空間の各々の縦列が位相エンコード位置に対応し、各々の横列が周波数エンコード位置に対応する。従って、ｋ空間の各々の点を位相エンコード位置及び周波数エンコード位置によって識別することができる。位相エンコード位置（又はｋ空間内の縦列）の数は、１本のエコーのサンプリング中に印加される位相エンコード勾配の数の関数である。位相エンコード位置の数が増すにつれて、必要な位相エンコード・ステップの数も増す。このようなものとして、位相エンコード位置の数が増加すると、走査時間も必然的に増加する。位相エンコード位置の数が減少すると走査時間は短縮するが、信号及び分解能の損失を伴う。従って、一実施形態では、本発明は、ｋ空間のブレードのホモダイン処理を含んでいる。この処理を図４に示す。説明の目的でｋ空間の単一のブレードについて説明するが、当業者は、ｋ空間の各々の取得されたブレードに対して本書で説明するホモダイン処理手法を応用し得ることを認められよう。
【００２３】
ｋ空間ブレード９０は、ｋ空間の中央領域９１又は原点を通って延在しており、ブレード９０の連続的な回転時に中央領域９１は過剰サンプリングされる。ｋ空間ブレード９０は好ましくは、所定の数のｋ空間データ線９２によって画定され、各々の線９２がサンプリング済み（実線）又は未サンプリング（破線）のいずれかとして特徴付けられる。さらに、図４には、ｋ空間ブレード９０の全てに満たない位相エンコード位置（ビュー）が、サンプリング済みのＭＲデータによって充填されていることが示されている。全てに満たない位相エンコード位置がサンプリングされるので、位相エンコード・ステップの数を減少させることができる。代替的には、各々のｋ空間データ線の部分的充填の後に、位相エンコード・ステップの総数が、ホモダイン処理を行なっていない取得のステップ数と等しくなるようにして、位相エンコード位置（従って位相エンコード・ステップ）の数を増加させることもできる。後に改めて説明するが、未サンプリングの位相エンコード位置は、サンプリング済みの位相エンコード位置に基づくデータで充填される。
【００２４】
多くの部分的フーリエ再構成手法を用いて未サンプリングの位相エンコード位置のデータを決定することができるが、好ましい一手法はホモダイン処理である。ホモダイン処理は、エルミート共役対称性のための一対のフィルタを用いて、サンプリング済みの位相エンコード位置のデータに基づいて未サンプリングの位相エンコード位置のデータを決定する。未サンプリングの位相エンコード位置のデータを提供することに加え、ホモダイン処理はまた、サンプリングされたデータに含まれて共鳴周波数、流れ及び運動の変動にしばしば帰属可能な位相の誤差又はばらつきを補正する。ホモダイン処理又は他の部分的フーリエ再構成手法では、他の場合であれば要求されるような位相エンコード・ステップを必要とせずに、未サンプリングの位相エンコード位置（図４では「Ｕ」として示す）を「実効的に」サンプリングすることができる。加えて、１本のエコーをサンプリングするのに必要な時間（ＴＥ）が実効的に短縮され、これによっても走査時間を短縮して被検体スループットを増大させることができる。
【００２５】
好ましくは、サンプリング済みの位相エンコード位置は、ｋ空間の横列１列の二分の一よりも僅かに多い部分を含む。この点で、ｋ空間データ線の中心の周りの位相エンコード位置について空間周波数データが取得される。このようにして、位相エンコード・ステップの数を二分の一よりも僅かに少ない分だけ減少し得るものと思量される。一方、印加される位相エンコード勾配の数を減少させずに定義し直して、殆ど全てのサンプリング済みの位相エンコード位置がｋ空間データ線の二分の一よりも僅かに多い部分に配置されるようにする。結果として、走査時間は短縮せず、信号対雑音が増大して空間分解能が高まる。事実上、位相エンコード・ステップの数を倍増せずに位相エンコード位置の数を倍増することができる。実際に、本発明によって、位相エンコード位置（位相ＦＯＶ）の数の倍増を上回ることが可能になるものと思量される。
【００２６】
図５には、ＰＲＯＰＥＬＬＥＲ取得によってＭＲデータを取得して、半径方向リプル・アーティファクトを実質的に含まない画像を再構成する手法の各ステップが掲げられている。工程９４（コンピュータで実行可能な言語で具現化されており、１台又は複数のコンピュータによって実行され得る）はステップ９６で開始して、ＰＲＯＰＥＬＬＥＲ取得の設定を行なう。ステップ９６では、ＰＲＯＰＥＬＬＥＲ取得のパラメータを、所定数の利用者入力から読み込むか又は自動的に決定する。これらのパラメータとしては、繰り返し時間（ＴＲ）、視野（ＦＯＶ）の大きさ、及びスライス配向、例えばアキシャル、頭部、コロナル又は斜方等の諸設定がある。上述のように、本発明は特に、コロナル、サジタル又は斜方ＰＲＯＰＥＬＬＥＲ走査で取得されたＭＲデータから再構成される画像での半径方向リプル（位相エイリアシング）アーティファクトを低減することに関する。
【００２７】
ＰＲＯＰＥＬＬＥＲ走査のパラメータを入力し、設定し、また他の場合には決定した後に、ステップ９８においてＭＲデータを取得する。エンコードされて励起されたエコーをサンプリングして用いて、互いから回転されておりｋ空間の中央領域を通って延在するｋ空間のブレードを部分的に充填する。次いで、ブレード・ホモダイン処理ステップ９００を実行して、図４に関して説明したように、サンプリング済みの位相エンコード位置から、未サンプリングの位相エンコード位置のデータを決定する。ブレード・ホモダイン処理１００の後に、ステップ１０２において画像空間位相補正を行なって、画像空間データの位相誤差を補正する。位相誤差補正１０２の後に、ステップ１０４において逆高速フーリエ変換（ＦＦＴ）を用いて、位相を補正した画像空間を下降（ダウン）サンプリングする。ステップ１０４において下降サンプリングが実行されて、ＦＯＶの位相エンコード方向での中心部分のみを記述するｋ空間のブレードを決定する。このステップは、中心部分以外からのデータが画像再構成時に除外されるように、ＦＯＶの位相エンコード方向での中心部分に対応するｋ空間のブレードを切り取ることにより達成される。このように、ブレード画像空間を切り取った後に、得られたブレード画像に逆ＦＦＴを施してｋ空間ブレードを復元する。次いで、このｋ空間ブレードをＰＲＯＰＥＬＬＥＲ再構成工程の残りの部分に通すと、工程が実行されて全てのｋ空間ブレードからの画像を再構成する。
【００２８】
明確に述べると、下降サンプリング１０４の後に、所定数のＰＲＯＰＥＬＬＥＲ再構成ステップを実行して運動を補正する。これらの運動補正ステップは、回転補正１０６、並進補正１０８、相関加重１１０、及び最後に、画像再構成１１２を含んでいる。画像再構成１１２の後に、ステップ１１４において工程は終了する。以上の各ステップの結果として、得られた再構成画像は半径方向リプル（位相エイリアシング）アーティファクトを実質的に含まない。
【００２９】
本発明の一つの利点は、ブレード・ホモダイン・ステップ及び下降サンプリング・ステップを従来のＰＲＯＰＥＬＬＥＲ再構成手法にシームレスに挿入し、しかもＳＮＲ及び空間分解能を高めて、各々の半径方向ＰＲＯＰＥＬＬＥＲブレードに沿った位相エイリアシングを除去し得ることである。このように、典型的には位相ＦＯＶを超えて延在する解剖学的構造から画像に現われる半径方向リプル・アーティファクトが低減する。すなわち、本発明は、ＰＲＯＰＥＬＬＥＲブレード同士の間の位相折り返しの組み合わせの結果として現われる得る位相エイリアシング・アーティファクトを低減する。
【００３０】
図６及び図７では、２枚のブレード画像が、図４に関して上で述べたブレード・ホモダイン処理ステップを示している。図６は、ｋ空間の半部で８枚のビュー及びｋ空間の残り半部で２枚のビューを示す１０枚のビューでのブレード画像である。ホモダイン処理の後に、１０枚のビューを、ｋ空間の半部で８枚のビュー及びｋ空間の残り半部で８枚のビューを示す１６枚のビューでのブレードへ変換する。ブレード・ホモダイン処理の結果として、６枚の追加のデータのビューが画像再構成に用いられ、これにより、データ取得時に６回分の追加で時間を費やす位相エンコード・ステップを行なわないでもＳＮＲ及び空間分解能が高まる。
【００３１】
図８には、従来のＰＲＯＰＥＬＬＥＲ取得を用いてＦＯＶを超えて延在する部分を有する球形ファントムから再構成される画像を示す。図示のように、位相エンコード方向に沿ってエイリアシングが生じている。一方、図９では、同じ球形ファントムを本発明に従って撮像したときに、図８で優勢であった位相エイリアシング・アーティファクト又は半径方向リプル・アーティファクトが除去されている。この点で、対象がＦＯＶを位相エンコード方向に超えて延在するにも拘わらず、位相エンコード方向には位相折り返しが生じない。尚、図８の画像は、ＰＲＯＰＥＬＬＥＲに基づく取得を用いて、エコー鎖長さ（ＥＴＬ）を２４、ｆｕｌｌ ＮＥＸ、及び位相ＦＯＶを２０ｃｍとして撮像されたことを特記しておく。図９の画像は、ＰＲＯＰＥＬＬＥＲを用いて、ＥＴＬ＝２４、ｈａｌｆ ＮＥＸ、及び位相ＦＯＶを４０ｃｍとして取得した。このように、図５に関して説明したホモダイン工程を用いて、ＦＯＶを位相エンコード方向に実効的に倍増させた。
【００３２】
図１０及び図１１はさらに、本発明によって達成される利点を示す。明確に述べると、図１０に示すように、ＦＯＶを位相エンコード方向に超えて延在する対象から取得される全てのブレードの従来のＰＲＯＰＥＬＬＥＲ再構成は、画像にリプルを生ずる。これらのリプルは、位相折り返しエイリアシングの結果である。一方、図１１は図１０で撮像した同じ対象の画像であるが、本書で説明した改良型ＰＲＯＰＥＬＬＥＲ取得を用いており、位相エイリアシング・リプルが存在しない。図１１の画像は、図１０の画像よりも遥かに一様であり、従って診断的価値が高い。
【００３３】
図１２及び図１３はさらに、本発明によって達成される利点を示す。明確に述べると、図１２に示すように、サジタル平面に沿って得られた脊椎画像の従来のＴ２−ＰＲＯＰＥＬＬＥＲ再構成は、図面では脊柱領域の縞として現われている半径方向リプル・アーティファクトを生じ得る。一方、本書で開示している改良型ＰＲＯＰＥＬＬＥＲ取得は、本発明によれば、図１３に示すようにこれらの半径方向リプル・アーティファクトを低減する。明確に述べると、図１２で見られる「縞」が図１３には存在しない。
【００３４】
従って、ＭＲＩ装置が、磁石のボアの周りに配置されており分極磁場を印加する複数の勾配コイルを含んでいる。ＲＦ送受信器システム及びＲＦスイッチが、パルス・モジュールによって制御されて、ＲＦコイル・アセンブリへ及びＲＦコイル・アセンブリからＲＦ信号を送受してＭＲ画像を取得する。ＭＲＩ装置はまた、ＦＯＶからコロナル・スライス配向及びサジタル・スライス配向のいずれかに沿ってＭＲデータのブレードを取得し、隣り合ったブレードが互いから回転されるようにＭＲデータのブレードでｋ空間を充填するようにプログラムされているコンピュータを含んでいる。コンピュータはまた、ｋ空間データを画像空間へ変換して、画像空間から半径方向リプル・アーティファクトを実質的に含まない画像を再構成するようにプログラムされている。
【００３５】
また、ＭＲデータ取得の方法が開示されており、この方法は、ＦＯＶからのＭＲデータ取得のためのサンプリングを行なうようにエコー鎖を励起するステップと、ｋ空間の中央領域を通って延在する複数のｋ空間データ線を有するｋ空間ブレードを部分的に充填するようにエコー鎖をサンプリングするステップとを含んでいる。この方法はさらに、部分的に充填されたｋ空間ブレードの残りを、エコー鎖からサンプリングされたデータから決定されるデータで充填するステップを含んでいる。励起するステップ、サンプリングするステップ、及び充填するステップは、複数のｋ空間ブレードがｋ空間において充填されるまで繰り返されて、各々のｋ空間ブレードがｋ空間の中央領域を通って延在し、且つ隣り合ったｋ空間ブレードがｋ空間の中央領域の周りで互いから回転されるようにｋ空間を充填する。
【００３６】
本発明はまた、コンピュータ読み取り可能な記憶媒体に記憶されたコンピュータ・プログラムとして具現化され、このコンピュータ・プログラムは、コンピュータによって実行されると、コンピュータに、ＦＯＶを超えて延在する部分を有する対象からｋ空間の複数のブレードを取得して、これらのｋ空間のブレードを撮像する、ことを行なわせる命令を有する。コンピュータはさらに、ＦＯＶの位相エンコード方向での中心部分を記述するブレード画像を単離して、ＦＯＶの位相エンコード方向での中心部分の外部に位置する位相エンコード位置からデータを除去するようにブレード画像を切り取る、ことを行なう。命令はまた、コンピュータに、切り取られたブレード画像をｋ空間の切り取られたブレードへ変換して、単離されたブレード画像に対応するｋ空間のブレードをｋ空間の切り取られたブレードで置き換える、ことを行なわせる。
【００３７】
本発明は好適実施形態について説明されており、明示的に記載した以外の均等構成、代替構成及び改変が可能であり特許請求の範囲内にあることが認められよう。また、図面の符号に対応する特許請求の範囲中の符号は、単に本願発明の理解をより容易にするために用いられているものであり、本願発明の範囲を狭める意図で用いられたものではない。そして、本願の特許請求の範囲に記載した事項は、明細書に組み込まれ、明細書の記載事項の一部となる。
【図面の簡単な説明】
【００３８】
【図１】本発明と共に用いられるＭＲイメージング・システムの模式的ブロック図である。
【図２】本発明で応用可能なＰＲＯＰＥＬＬＥＲパルス・シーケンスの例の一部を示す図である。
【図３】図２のＰＲＯＰＥＬＬＥＲパルス・シーケンスの例を用いて取得されたＭＲデータで充填されるｋ空間の模式図である。
【図４】本発明のもう一つの観点に従ってホモダイン処理を用いて充填されるｋ空間の単一のブレードの模式図である。
【図５】本発明のもう一つの観点によるＰＲＯＰＥＬＬＥＲ再構成の各ステップを示す流れ図である。
【図６】ｋ空間の部分的に充填されたブレードを示す画像である。
【図７】元々部分的に充填されており、続いてホモダイン処理を用いて充填されたｋ空間のブレードを示す画像である。
【図８】位相ＦＯＶを超えて延在する部分を有する球形ファントムから従来のＰＲＯＰＥＬＬＥＲ取得を用いて再構成される単一のブレード画像である。
【図９】本発明の一観点に従って位相エイリアシングを低減した図８で撮像された球形ファントムから再構成される単一のブレード画像である。
【図１０】様々な配向で取得されたデータのブレードから従来のＰＲＯＰＥＬＬＥＲ再構成手法を用いて再構成されたファントムの画像である。
【図１１】図１０で撮像されたファントムの画像であるが、本発明の一観点による改良型ＰＲＯＰＥＬＬＥＲ再構成手法によって様々な配向で取得されたデータのブレードからの画像である。
【図１２】従来のＰＲＯＰＥＬＬＥＲ取得を用いてサジタル平面に沿って得られた脊椎画像である。
【図１３】本発明の一観点による改良型ＰＲＯＰＥＬＬＥＲ取得を用いてやはりサジタル平面に沿って得られた脊椎画像であって、図１２の画像に見られる半径方向リプル（位相エイリアシング）アーティファクトの存在しない画像である。
【符号の説明】
【００３９】
１０ ＭＲイメージング・システム
１２ 操作コンソール
１３ 入力装置
１４ 制御パネル
１６ 表示画面
１８ リンク
２０ コンピュータ・システム
２０ａ、３２ａ バックプレイン
２２ 画像プロセッサ・モジュール
２４、３６ ＣＰＵモジュール
２６、６６ メモリ・モジュール
２８ ディスク記憶装置
３０ テープ・ドライブ
３２ システム制御
３４ 高速シリアル・リンク
３８ パルス発生器モジュール
４０ シリアル・リンク
４２ 勾配増幅器
４４ 生理学的取得制御器
４６ 走査室インタフェイス回路
４８ 患者配置システム
５０ 勾配コイル・アセンブリ
５２ 磁石アセンブリ
５４ 分極用磁石
５６ 全身型ＲＦコイル
５８ 送受信器モジュール
６０ ＲＦ増幅器
６２ 送受信（Ｔ／Ｒ）スイッチ
６４ 前置増幅器
６８ アレイ・プロセッサ
７０ ＰＲＯＰＥＬＬＥＲパルス・シーケンス
７２ スピン・エコー
７４ ＲＦ再収束パルス
７６ 周波数エンコード・パルス
７８ ｋ空間
８０ データ線
８８ ｋ空間の中心
９０ ｋ空間ブレード
９１ 中央領域
９２ ｋ空間データ線
【特許請求の範囲】
【請求項１】
磁石（５４）のボアの周りに配置されており分極磁場を印加する複数の勾配コイル（５０）と、無線周波数（ＲＦ）コイル・アセンブリ（５６）へ無線周波数信号を送信して磁気共鳴画像を取得するようにパルス・モジュール（３８）により制御される無線周波数送受信器システム（５８）及び無線周波数スイッチ（６２）と、
任意のスライス配向に沿って視野（ＦＯＶ）から磁気共鳴データのブレード（８２、９０）を取得し、
隣り合ったブレード（８２、９０）が互いから回転するようにして磁気共鳴データの前記ブレード（８２、９０）でｋ空間（７８）を充填し、
前記ｋ空間（７８）を画像空間へ変換して、
前記画像空間から半径方向リプル・アーティファクトを実質的に含まない画像を再構成する
ようにプログラムされているコンピュータ（２０）と、
を備えた磁気共鳴（ＭＲ）撮像装置。
【請求項２】
前記コンピュータ（２０）はさらに、各々のブレード（８２、９０）がｋ空間（７８）の中央領域（８８、９１）を通って延在するように磁気共鳴データの前記ブレード（８２、９０）を配向させるようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項３】
前記コンピュータ（２０）はさらに、位相エイリアシング・アーティファクトを実質的に含まない前記画像を再構成するようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項４】
前記コンピュータ（２０）はさらに、磁気共鳴データのブレード（８２、９０）を部分的に取得するようにエコーを部分的にサンプリングして、該エコーから部分的に取得されたデータから磁気共鳴データの前記ブレード（８２、９０）の残りを決定するようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項５】
前記コンピュータ（２０）はさらに、
磁気共鳴データの各々のブレード（８２、９０）についてブレード画像を形成し（１０２）、
前記視野の位相エンコード方向での中央領域に対応する前記ブレード画像を単離し（１０４）、
前記視野の前記位相エンコード方向での前記中央領域の外部から位相エンコード・データを除去するように前記ブレード画像を切り取り（１０４）、
前記切り取られたブレード画像を磁気共鳴データの切り取られたブレード（８２、９０）へ変換し（１４）、
前記単離したブレード画像を形成した元となった磁気共鳴データの前記ブレード（８２、９０）を磁気共鳴データの前記切り取られたブレード（８２、９０）で置き換えて、
磁気共鳴データの前記切り取られたブレード（８２、９０）を画像空間へ変換する（１０６〜１１２）
ようにプログラムされている、請求項４に記載の磁気共鳴撮像装置。
【請求項６】
前記コンピュータ（２０）はさらに、前記切り取られたブレード画像を磁気共鳴データの前記切り取られたブレード（８２、９０）へ変換する逆高速フーリエ変換（ＦＦＴ）を施す（１０４）ようにプログラムされている、請求項５に記載の磁気共鳴撮像装置。
【請求項７】
前記コンピュータ（２０）はさらに、画像再構成の前に運動アーティファクトを補正する（１０２）ようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項８】
前記コンピュータ（２０）はさらに、前記視野の外部に少なくとも部分的に延在する任意の解剖学的関心対象から磁気共鳴データを取得する（９８）ようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項９】
前記コンピュータ（２０）はさらに、ＰＲＯＰＥＬＬＥＲ取得により磁気共鳴データの前記ブレード（８２、９０）を取得するようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項１】
磁石（５４）のボアの周りに配置されており分極磁場を印加する複数の勾配コイル（５０）と、無線周波数（ＲＦ）コイル・アセンブリ（５６）へ無線周波数信号を送信して磁気共鳴画像を取得するようにパルス・モジュール（３８）により制御される無線周波数送受信器システム（５８）及び無線周波数スイッチ（６２）と、
任意のスライス配向に沿って視野（ＦＯＶ）から磁気共鳴データのブレード（８２、９０）を取得し、
隣り合ったブレード（８２、９０）が互いから回転するようにして磁気共鳴データの前記ブレード（８２、９０）でｋ空間（７８）を充填し、
前記ｋ空間（７８）を画像空間へ変換して、
前記画像空間から半径方向リプル・アーティファクトを実質的に含まない画像を再構成する
ようにプログラムされているコンピュータ（２０）と、
を備えた磁気共鳴（ＭＲ）撮像装置。
【請求項２】
前記コンピュータ（２０）はさらに、各々のブレード（８２、９０）がｋ空間（７８）の中央領域（８８、９１）を通って延在するように磁気共鳴データの前記ブレード（８２、９０）を配向させるようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項３】
前記コンピュータ（２０）はさらに、位相エイリアシング・アーティファクトを実質的に含まない前記画像を再構成するようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項４】
前記コンピュータ（２０）はさらに、磁気共鳴データのブレード（８２、９０）を部分的に取得するようにエコーを部分的にサンプリングして、該エコーから部分的に取得されたデータから磁気共鳴データの前記ブレード（８２、９０）の残りを決定するようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項５】
前記コンピュータ（２０）はさらに、
磁気共鳴データの各々のブレード（８２、９０）についてブレード画像を形成し（１０２）、
前記視野の位相エンコード方向での中央領域に対応する前記ブレード画像を単離し（１０４）、
前記視野の前記位相エンコード方向での前記中央領域の外部から位相エンコード・データを除去するように前記ブレード画像を切り取り（１０４）、
前記切り取られたブレード画像を磁気共鳴データの切り取られたブレード（８２、９０）へ変換し（１４）、
前記単離したブレード画像を形成した元となった磁気共鳴データの前記ブレード（８２、９０）を磁気共鳴データの前記切り取られたブレード（８２、９０）で置き換えて、
磁気共鳴データの前記切り取られたブレード（８２、９０）を画像空間へ変換する（１０６〜１１２）
ようにプログラムされている、請求項４に記載の磁気共鳴撮像装置。
【請求項６】
前記コンピュータ（２０）はさらに、前記切り取られたブレード画像を磁気共鳴データの前記切り取られたブレード（８２、９０）へ変換する逆高速フーリエ変換（ＦＦＴ）を施す（１０４）ようにプログラムされている、請求項５に記載の磁気共鳴撮像装置。
【請求項７】
前記コンピュータ（２０）はさらに、画像再構成の前に運動アーティファクトを補正する（１０２）ようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項８】
前記コンピュータ（２０）はさらに、前記視野の外部に少なくとも部分的に延在する任意の解剖学的関心対象から磁気共鳴データを取得する（９８）ようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【請求項９】
前記コンピュータ（２０）はさらに、ＰＲＯＰＥＬＬＥＲ取得により磁気共鳴データの前記ブレード（８２、９０）を取得するようにプログラムされている、請求項１に記載の磁気共鳴撮像装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【公開番号】特開２００６−２２３８６４（Ｐ２００６−２２３８６４Ａ）
【公開日】平成１８年８月３１日（２００６．８．３１）
【国際特許分類】
【出願番号】特願２００６−３６０１６（Ｐ２００６−３６０１６）
【出願日】平成１８年２月１４日（２００６．２．１４）
【出願人】（３９００４１５４２）ゼネラル・エレクトリック・カンパニイ (6,332)
【氏名又は名称原語表記】ＧＥＮＥＲＡＬ ＥＬＥＣＴＲＩＣ ＣＯＭＰＡＮＹ
【Ｆターム（参考）】
【公開日】平成１８年８月３１日（２００６．８．３１）
【国際特許分類】
【出願日】平成１８年２月１４日（２００６．２．１４）
【出願人】（３９００４１５４２）ゼネラル・エレクトリック・カンパニイ (6,332)
【氏名又は名称原語表記】ＧＥＮＥＲＡＬ ＥＬＥＣＴＲＩＣ ＣＯＭＰＡＮＹ
【Ｆターム（参考）】
[ Back to top ]