(57)【特許請求の範囲】
【請求項１】
静磁場中におかれた被検体に対してＲＦパルスおよび傾斜磁場パルスを印加し、関心領域内の所望代謝物からの磁気共鳴信号を収集する磁気共鳴診断装置において、
前記被検体中で不要な信号を抑圧させるべくＲＦパルスを印加するＲＦパルス印加手段と、
前記ＲＦパルスにて励起された磁気共鳴信号の位相分散をするべく傾斜磁場パルスを印加する傾斜磁場パルス印加手段と、
前記ＲＦパルスの強度時間積分値を、前記関心領域を含む磁気共鳴画像の画像値から算出する計算機とを備え、
前記磁気共鳴画像には前記関心領域を設定するために収集されたロケータ画像とＴ１強調画像とＴ２強調画像とが含まれることを特徴とする磁気共鳴診断装置。
【請求項２】
静磁場中におかれた被検体に対してＲＦパルスおよび傾斜磁場パルスを印加し、関心領域内の所望代謝物からの磁気共鳴信号を収集する磁気共鳴診断装置において、
前記被検体中で不要な信号を抑圧させるべくＲＦパルスを印加するＲＦパルス印加手段と、
前記ＲＦパルスにて励起された磁気共鳴信号の位相分散をするべく傾斜磁場パルスを印加する傾斜磁場パルス印加手段と、
前記ＲＦパルスの強度時間積分値を、前記関心領域を含む磁気共鳴画像の画像値から算出する計算機とを備え、
前記計算機は、前記関心領域の画素値と緩和時間が既知であるリファレンス領域の画素値とに基づいて、前記リファレンス領域の緩和時間に対する前記関心領域の緩和時間の変化及びプロトン密度の変化を計算することを特徴とする磁気共鳴診断装置。
【請求項３】
静磁場中におかれた被検体に対してＲＦパルスおよび傾斜磁場パルスを印加し、関心領域内の所望代謝物からの磁気共鳴信号を収集する磁気共鳴診断装置において、
前記被検体中で診断の対象領域外から発生する磁気共鳴共鳴信号を飽和させるべくＲＦパルスを印加するＲＦパルス印加手段と、
前記ＲＦパルスにて励起された磁気共鳴信号の位相分散をするべく傾斜磁場パルスを印加する傾斜磁場パルス印加手段と、
前記ＲＦパルスの強度時間積分値を、前記関心領域を含む磁気共鳴画像の画像値から算出する計算機とを備え、
前記磁気共鳴画像には前記関心領域を設定するために収集されたロケータ画像とＴ１強調画像とＴ２強調画像とが含まれることを特徴とする磁気共鳴診断装置。
【請求項４】
静磁場中におかれた被検体に対してＲＦパルスおよび傾斜磁場パルスを印加し、関心領域内の所望代謝物からの磁気共鳴信号を収集する磁気共鳴診断装置において、
前記被検体中で診断の対象領域外から発生する磁気共鳴共鳴信号を飽和させるべくＲＦパルスを印加するＲＦパルス印加手段と、
前記ＲＦパルスにて励起された磁気共鳴信号の位相分散をするべく傾斜磁場パルスを印加する傾斜磁場パルス印加手段と、
前記ＲＦパルスの強度時間積分値を、前記関心領域を含む磁気共鳴画像の画像値から算出する計算機とを備え、
前記計算機は、前記関心領域の画素値と緩和時間が既知であるリファレンス領域の画素値とに基づいて、前記リファレンス領域の緩和時間に対する前記関心領域の緩和時間の変化及びプロトン密度の変化を計算することを特徴とする磁気共鳴診断装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、被検体の代謝物からの磁気共鳴信号を検出する磁気共鳴診断装置に係り、特に、不要水信号を抑制するためのＲＦパルスのフリップ角および局所励起領域外の信号を飽和するためのＲＦパルスのフリップ角を正確に、且つ調整時間を省略して計算機内部で求めることが可能な磁気共鳴診断装置に関する。
【０００２】
【従来の技術】
近年、疾病の早期診断（特に腫瘍の進行度の判定）に有効な方法として、１Ｈ−ＭＲＳ（ＭＲスペクトロスコピー）が行われている。１Ｈ−ＭＲＳは、ある基準集からの共鳴周波数の隔たりをｐｐｍで表した磁気共鳴スペクトルから、１Ｈの核を含む生体内の様々な物質の分子構造、化学環境、濃度等の情報を得る方法であり、典型的な利用法としては、代謝物としてのＮＡＡ（Ｎ−アセチルアスパラギン酸）、Ｃｈｏ（コリン）、ＰＣｒ／Ｃｒ（クレアチンリン酸/クレアチン）、Ｇｌｘ（グルタミン酸およびグルタミン）、Ｌａｃ（乳酸）、ｍＩ（ミオイノシトール）等を観測することに使われている。
【０００３】
この１Ｈ−ＭＲＳにおいて、生体内に存在するＮＡＡ等のプロトン代謝物の濃度は、生体水の濃度に比べて４桁程度小さく、通常のＡ／Ｄ変換器のビット数では代謝物からの信号を十分検出することができないため、通常、水信号を選択的に抑圧することが必須とされる。
【０００４】
この水信号を抑圧する方法として、さまざまな方法が提案されているが、一般的に、ＣＨＥＳＳ（化学シフト選択）パルスと呼ばれるＲＦパルスを使用する方法が実施されている。これは、図４の水信号抑圧区間ｔｃｈｅｓｓに示されているように、水のプロトンの周波数帯域のみを選択的に励起する数十ｍｓｅｃのＲＦパルス（ＣＨＥＳＳパルス）Ｐ１，Ｐ２，Ｐ３を印加し、水のプロトンの磁化のみを横平面に倒し、その状態で、傾斜磁場パルスＧ１，Ｇ２，Ｇ３を印加し、横磁化をスポイルすることで水信号を抑圧する方法である。
【０００５】
これら３本のＣＨＥＳＳパルスＰ１，Ｐ２，Ｐ３のうち、最初の２本のＣＨＥＳＳパルスＰ１，Ｐ２の強度時間積分値は水のプロトンの磁化を横平面に倒すためにフリップ角が９０°になるように設定され、３本目のＣＨＥＳＳパルスＰ３のフリップ角αは、水プロトンからの信号が実効的にゼロ又は最小化するようにプリスキャンで厳密に決定される。このように水抑制を実施することで、水信号に代謝物信号が埋もれることなく、代謝物の信号観測は可能となる。
【０００６】
また、脳の頭表の近傍あるいは、前立腺などの腹部領域で脂肪が大量に含まれている領域の近傍では、水信号と同様に、観測領域内への脂肪信号の混入が大きな問題となる。脂肪信号低減のために、通常、プリサチュレーション法が水信号抑圧法と共に併用されている。
【０００７】
このパルス系列は、図４の区間ｔｐｒｅｓａｔに示している。ＭＲＳの関心領域（ＲＯＩ）以外の領域の磁化を選択的に倒すために、周波数帯域△Ｆを持つスライス選択ＲＦパルス（以下、プリサットパルスと称する）Ｐ４，Ｐ５を傾斜磁場Ｇ４，Ｇ６とともに印加し、それぞれの後に、倒れた磁化をスポイルするためのスポイル傾斜磁場Ｇ５，Ｇ７を印加する。この脂肪抑圧区間ｔｐｒｅｓａｔに、代謝物の磁化を励起する局所励起区間ｔｌｏｃと、信号を観測する区間ｔａｃｑとが続く。
【０００８】
区間ｔｌｏｃは、フリップ角がそれぞれ９０°，１８０°，１８０°に成るように時間強度積分値が調整されているＲＦパルスＬ１，Ｌ２，Ｌ３によりエコー信号を発生させるダブルＳＥ（スピンエコー）法である。印加タイミングとしては、ｔｌｏｃ２＝ｔｌｏｃ１＋ｔｌｏｃ３である。また、局所励起用のシーケンスとして３つの９０°パルスからなるＳＴＥ（Stimulated Echo）法も使用される。また、ＲＦパルスＬ１，Ｌ２，Ｌ３と同時に印加されるＧｓ１，Ｇｓ２，Ｇｓ３は磁気共鳴信号を観測する領域を決定するためのスライス選択傾斜磁場である。Ｇsp1 ，Ｇsp2 が不要信号をスポイルするための傾斜磁場である。
【０００９】
またＧenc1，Ｇenc2の位相エンコード傾斜磁場を印加し、データ収集ごとに強度を変化させることでＭＲＳＩが可能となる。Ｇenc1，Ｇenc2を印加しない場合はＧsl1 ，Ｇsl2 ，Ｇsl3 により決まる単一の領域からの磁気共鳴信号を観測するシングルボクセルのＭＲＳが可能となる。ｓ１、ｓ２、ｓ３はｘ、ｙ、ｚの空間軸のいずれかを取得部位に応じて選ぶことができ、スライス選択傾斜磁場を複数軸選択することでオブリーク収集も可能となる。またスポイル傾斜磁場も不要信号が十分スポイルされるように、ＭＲＩ装置ごとに印加強度、印加時間、印加軸を選んで使用する。
【００１０】
以上の方法により１Ｈ−ＭＲＳにおいて、脂肪信号が対象領域の近傍にある場合でも良好なスペクトルが取得可能である。水信号抑制区間ｔｃｈｅｓｓの３本目のＣＨＥＳＳパルスＰ３のフリップ角αを、水の緩和時間を既知として調整を行わないで固定する方法があるが、組織（特に病変部位）により緩和時間が異なるために、正確なフリップ角を求めるためには調整が必要である。
【００１１】
送信パワーを変えながら水信号を最小とする手動の方法があるが、調整時間が長くかかることと操作者により差が出るなどの問題点がある。水抑制の３本目のパルスＰ３のフリップ角αを振りながら図４のシーケンスにより１Ｈ−ＭＲＳのデータ収集を行い、横軸がフリップ角で縦軸が水信号強度のグラフからαを求める方法があるが、この場合も「フリップ角を振る回数」×「繰り返し時間ＴＲ」という比較的長い時間を調整に要してしまう。
【００１２】
また、自動の場合には、Ｔ１が長い組織では、短い繰り返し時間ＴＲ（通常１Ｈ−ＭＲＳの計測の場合には１，５００〜３，０００ｍｓｅｃである）だと、縦磁化の緩和が不十分のため正確なフリップ角αが求まらないという問題がある。Ｔ１が長い組織に対して正確なフリップ角を求めるためには調整のときのみＴＲを６，０００以上に延長する必要があるが調整時間が余分にかかってしまう。また通常、１Ｈ−ＭＲＳの検査はシングルボクセルの場合は健側と患側の２種類以上行っているが、その度に水抑制のフリップ角調整を行っている。
【００１３】
このように、従来の磁気共鳴画像診断装置において１Ｈ−ＭＲＳを行う場合、水抑制のＣＨＥＳＳパルスＰ３のフリップ角αを正確に求めるには調整時間がかかるという欠点があった。
【００１４】
【発明が解決しようとする課題】
本発明の目的は、水抑制のためのＣＨＥＳＳパルスのフリップ角の決定を効率的に行うことにある。
【００１５】
【課題を解決するための手段】
本発明は、静磁場中におかれた被検体に対してＲＦパルスおよび傾斜磁場パルスを印加し、関心領域内の所望代謝物からの磁気共鳴信号を収集する磁気共鳴診断装置において、前記被検体中で不要な信号を抑圧させるべくＲＦパルスを印加するＲＦパルス印加手段と、前記ＲＦパルスにて励起された磁気共鳴信号の位相分散をするべく傾斜磁場パルスを印加する傾斜磁場パルス印加手段と、前記ＲＦパルスの強度時間積分値を、前記関心領域を含む磁気共鳴画像の画像値から算出する計算機とを備え、前記磁気共鳴画像には前記関心領域を設定するために収集されたロケータ画像とＴ１強調画像とＴ２強調画像とが含まれることを特徴とする。
【００１６】
【発明の実施の形態】
以下、図面を参照して本発明による装置を好ましい実施形態により説明する。
【００１７】
図１は、本実施形態に係る磁気共鳴映像装置の構成を示すブロック図である。同図において、静磁場磁石２０１は、励磁用電源２０２に駆動され、主勾配コイル群２０３およびシールドコイル群２０４は勾配コイル用電源２０５にて駆動される。これらにより、被検体２０６には一様な静磁場とそれと同一方向で互いに直交する３方向に線形磁場勾配を持つ傾斜磁場が印加される。主勾配コイル群２０３とシールドコイル群２０４は直列接続され、共通の勾配コイル用電源２０５に接続されてもよいし、上下左右の各コイルエレメントごとに複数の勾配コイル用電源２０５に接続され分割駆動されてもよい。
【００１８】
渦磁場時間応答補償用入力信号は、渦補償回路２０７にて生成される。シムコイル群２０８は、シムコイル用電源２０９に駆動され、静磁場の均一性が調整される。送信部２１０は高周波信号を出力するものであり、この高周波信号はプローブ２１１に送られ、被検体２０６に高周波磁場が印加される。このとき、プローブ２１１は送受両用でも送受信別々に設けてもよい。また、プローブ２１１と主コイル群２０３との間には、高周波シールド２１２が配設されている。
【００１９】
プローブ２１１で受信された磁気共鳴信号は、受信部２１３で検波された後、データ収集部２１４に転送され、計算機２１５に送られる。そして、上述した勾配コイル用電源２０５、シムコイル用電源２０９、送信部２１０、受信部２１３、データ収集部２１４は全てシーケンスコントローラ２１６により制御されており、システムコントローラ２１６は計算機２１５に、また、計算機２１５は、入力装置２１７により制御されている。計算機２１５は、データ収集部２１４から送られたデータのフーリエ変換等が行われ、被検体内部の所望化合物のＭＲＳを生成すると共に、水信号抑圧区間の３本目のＣＨＥＳＳパルスＰ３のフリップ角αを決定するための処理も行う。この処理についての詳細は後述する。そして、得られたスペクトルは画像ディスプレイ２１８に表示される。
【００２０】
図２（ａ）、図２（ｂ）、図２（ｃ）は、本発明の第１実施形態で水抑制のためのＣＨＥＳＳパルスのフリップ角αを決定するために使う１Ｈ−ＭＲＳの撮影プラン画像である。本実施形態は、ＣＨＥＳＳパルスのフリップ角が、１Ｈ−ＭＲＳの撮影以前に撮影された異なる３種類のシーケンスで得られた３種類の画像から求めることを特徴としている。
【００２１】
通常、１Ｈ−ＭＲＳの検査を行う際には、事前に、１Ｈ−ＭＲＳでの関心領域（以下、ＲＯＩと称する）を特定するために、イメージングを行う。通常は、グラディエントエコーによるロケータ画像１０１（図１（ａ））を取り、スピンエコー法によるＴ１強調画像１０２（図１（ｂ））およびＴ２強調画像１０３（図１（ｃ））を取得する。画像１０２，１０３は、例えば脳室の部分をスライスするアキシャル画像とし、１０１をサジタル画像とする。これらの画像からＣＨＥＳＳパルスのフリップ角αを決める。
【００２２】
１０４、１０５、１０６は、１Ｈ−ＭＲＳを取得するＲＯＩであり、つまり図４のＧsl1 ，Ｇsl2，Ｇsl3で励起される領域を示している。１０７、１０８、１０９は、フリップ角を計算するためのリファレンス領域を示す。ＲＯＩを決定する際には、１０１〜１０３のうちのいずれか一つの画像に対して、例えば、画像１０３に対して、ＲＯＩ１０６を腫瘍等の病変部位に入力装置２１７によりユーザーが設定する。
【００２３】
リファレンス領域１０９は、他の２枚の画像１０１，１０２が共にこのＲＯＩ１０６を含むように、ユーザーが入力装置２１７により設定する。または計算機２１５でリファレンス領域１０９が他の２枚の画像１０１，１０２に含まれるように制御しても良い。このようにしてリファレンス領域１０９は、緩和時間が既知である脳実質部に設定される。
【００２４】
ここで、ＲＯＩ１０４、１０５、１０６内のピクセルの平均値をそれぞれａ１，ａ２，ａ３とし、リファレンス領域１０７、１０８、１０９内のピクセルの平均値をそれぞれｂ１，ｂ２，ｂ３とする。画像１０１〜１０３が、複数枚の断面で収集されている場合はそれぞれに隣接する画像（ただしこれらの画像がＲＯＩおよびリファレンス領域内に含まれていることが条件である）も使用してピクセルの平均値を求めると、以下で求めるフリップ角αの精度が上がる。
【００２５】
以下に、ピクセル値ａ１，ａ２，ａ３，ｂ１，ｂ２，ｂ３を使用して、フリップ角αを求める方法を説明する。エコー時間ＴＥ，繰り返し時間ＴＲのスピンエコー法の信号強度Ｉは、プロトン密度をＭ、組織の縦緩和時間をＴ１，横緩和時間をＴ２とすると次式（１）で表される。
【００２６】
【数１】
0004718698.tif 000002
【００２７】
基準となる組織の緩和時間Ｔ１，Ｔ２およびプロトン密度Ｍに対して、別の組織においてそれぞれｄＴ１，ｄＴ２，ｄＭだけ変化した場合の信号強度の変化ｄＩは、Ｅｌｓｔｅｒ  ＡＤ（J Comput Assist Tomogr 12:130,1988）に記述されているように、次式（２）で表される。
【００２８】
【数２】
0004718698.tif 000003
【００２９】
ここでＳＴ１、ＳＴ２、ＳＭは次式（３）〜（５）で表される。
【００３０】
【数３】
0004718698.tif 000004
【００３１】
またフリップ角θ、エコー時間ＴＥ，繰り返し時間ＴＲ，グラディエントエコーの信号強度Ｉは、プロトン密度をＭ、組織の縦緩和時間をＴ１，横緩和時間をＴ２とすると次式（６）で表される。
【００３２】
【数４】
0004718698.tif 000005
【００３３】
Ｔ１，Ｔ２，ＭがそれぞれｄＴ１，ｄＴ２，ｄＭだけ変化した場合の信号強度の変化ｄＩは、スピンエコー法と同様にして次式（７）で表される。
【００３４】
【数５】
0004718698.tif 000006
【００３５】
ここでＧＴ１，ＧＴ２，ＧＭは、次式（８）〜（１０）で表される。
【００３６】
【数６】
0004718698.tif 000007
【００３７】
上述したようにリファレンス領域１０７〜１０９を脳実質部に設定しており、これらの緩和時間Ｔ１，Ｔ２は、多くの文献で述べられているように、既知である。ＲＯＩ１０４〜１０６のリファレンス領域１０７〜１０９に対する緩和時間、プロトン密度の変化が求める解である。
【００３８】
求める解ｄＴ１／Ｔ１，ｄＴ２／Ｔ２，ｄＭ／Ｍをそれぞれｘｘ、ｙｙ、ｚｚと置く。基準となるＴ１，Ｔ２を決定したため、画像１０１に対するＧＴ１，ＧＴ２，ＧＭ、およびＲＯＩ１０２、１０３に対するＳＴ１，ＳＴ２，ＳＭが求まり、これらをＡ11，Ａ12，Ａ13，Ａ21，Ａ22，Ａ23，Ａ31，Ａ32，Ａ33とし、式（７）および（２）に代入すると、次式（１１）〜（１３）の１次元連立方程式が成立する。
【００３９】
【数７】
0004718698.tif 000008
【００４０】
上記連立方程式を解くことでＲＯＩの縦緩和時間Ｔ１'をＴ１+ｄＴ１として求めることができる。
【００４１】
従って、水信号抑制区間の３本目のＣＨＥＳＳパルスＰ３の中心印加時刻と局所励起パルスＬ１の中心印加時刻までの時間をｔとすると、３本目のＣＨＥＳＳパルスＰ３のフリップ角αは、次式（１４）で与えられる。
【００４２】
【数８】
0004718698.tif 000009
【００４３】
またパルス強度と時間の関係をリニアで近似した場合αは次式（１５）で表される。
【００４４】
【数９】
0004718698.tif 000010
【００４５】
このように３枚の画像からフリップ角の最適値又はそれに近い値を求めることができるので、従来のように、水抑制の３本目のパルスＰ３のフリップ角αを振りながら１Ｈ−ＭＲＳのデータ収集を繰り返すよりも、水抑制のためのＣＨＥＳＳパルスＰ３のフリップ角αの決定を効率的に行うことができる。
【００４６】
図３は本発明の第２の実施形態を示す１Ｈ−ＭＲＳの撮影プラン画像を示す図である。プリサットを印加した場合のＣＨＥＳＳパルスのフリップ角、およびプリサットパルスのフリップ角が、１Ｈ−ＭＲＳの撮影以前に撮影された異なる３種類のシーケンスで得られた画像から求められることを示す。図４は従来例のパルスシーケンス図であるが本発明の説明において使用する。
【００４７】
次に第２の実施形態に関して詳細に説明する。図３（ａ）の３０１がロケータ画像、図３（ｂ）の３０２がスピンエコー法によるＴ１強調画像、図３（ｃ）の３０３がスピンエコー法によるＴ２強調画像である。３０４、３０５、３０６が図４のシーケンスのＧsl1 ，Ｇsl2 ，Ｇsl3 で励起される関心領域であり、３０７、３０８、３０９がフリップ角を計算するためのリファレンス領域を示している。３１０〜３１５の点線部分が区間ｔｐｒｅｓａｔのスライス選択用のＲＦパルスＰ４，Ｐ５と傾斜磁場Ｇ４，Ｇ６とで選択されるプリサット領域である。３１０〜３１５の実線領域が、以下の計算でピクセルの平均値を取るための領域である。３０４〜３０９および３１０〜３１５の実線部分がグラフィカルなインターフェースでユーザーにより指定される。リファレンス領域３０７、３０８、３０９のピクセルの平均値をそれぞれｂ１，ｂ２，ｂ３とし、ＲＯＩ３１０、３１１、３１２のピクセルの平均値をａ１，ａ２，ａ３とする。
【００４８】
リファレンス領域３０７〜３０９の脳実質部に関する既知の緩和時間Ｔ１，Ｔ２を使用し、ＲＯＩ３１０〜３１２の緩和時間をＴ１+ｄＴ１，Ｔ２+ｄＴ２、プロトン密度をＭ＋ｄＭと置き、画像３０１〜３０３が第１の実施形態の説明で用いたＴＥ，ＴＲ，θで得られたとすると、プリサットフリップ角αは第１の実施形態と同様に、式（１４）ないし式（１５）で求められる。ただし式（１４）、（１５）のtはプリサツトパルスＰ４の中心印加時刻から局所励起パルスＬ１の中心印加時刻までの時間を示す。
【００４９】
２本目のプリサットフリップ角の強度も同様に求めることができる。また水抑制パルスのフリップ角も第１の実施形態と同様に求めることができる。第２の実施形態では、１Ｈ−ＭＲＳに関して説明を行ったが、腹部等でＣＨＥＳＳパルスを使用した脂肪抑制を行う場合にも、脂肪抑制パルスのフリップ角を求める際にも適用できる。この場合にはリファレンスとなる物質として腹部の実質部を指定して、緩和時間も腹部実質部の緩和時間を使用する必要がある。上記２つの実施形態においてＣＨＥＳＳないしプリサツトフリップ角をグラディエントエコー法によるサジタル画像、スピンエコー法によるＴ１強調アキシャル画像，Ｔ２強調アキシャル画像の３種類の画像から求めているが、ＴＲ，ＴＥ（グラディエントエコー法ではＴＲ，ＴＥ，θ）が異なっていればいかなるシーケンスで得られたいかなるスライスの画像を使用しても良い。
【００５０】
ただし画像の中にＲＯＩおよびリファレンス領域が含まれている必要がある。また４種類以上のシーケンスを使用して画像を取得した場合はＳＮＲの高い画像から３つを使用すれば決定精度が高くなるため望ましい。また上記２つの実施形態において、１Ｈ−ＭＲＳに関する実施形態ではスピンエコー法によるシングルボクセル法のみで説明を行ったが、マルチボクセル法に関しても適用することができる。またＣＨＥＳＳ以降のパルスシーケンスがいかなるシーケンスに対しても適用可能である。
【００５１】
また３種類の画像でなく２種類の画像あるいは１種類の画像からもＣＨＥＳＳないしプリサットフリツプ角を求めることができる。例えば図２の場合で１０１、１０２の画像から求める場合を説明する。１０１がＴＥの短いグラディエントエコー画像であれば式（９）を０と近似しても良く、式（１１）のＡ１２の項は０となる。また１０２がＴＥが十分短いスピンエコー法によるＴ１強調画像であれば、式（４）を０と近似しても良く、式（１２）のＡ２２の項は０となる。つまり式（１１）、（１２）はｘｘ，ｚｚをパラメータとした２次元の連立方程式となり、所望の緩和時間Ｔ１'を求めることができ、式（１４）ないし（１５）により所望のαを求めることができる。このようにＴＥの短い任意の２つのシーケンスによる画像からフリップ角αを求めることができる。またＴＥの短いシーケンスを使用して、プロトン密度は脳実質部と関心領域でほぼ一定であるとすると単一の画像からαを求めることができる。つまり式（１２）の第２項、第３項を０と近似すれば１種類の画像からＴ１'を求めることができ、式（１４）ないし（１５）により所望のαを求めることができる。またこれらの１つまたは２つの画像からαを求める方法において、リファレンス領域を指定して、リファレンス領域の緩和時間を既知としているが、ある領域の水信号が最小となるフリップ角を従来の方法により求めて、ＲＯＩのαを上記のいずれかの方法により求めても良い。例えば１０７〜１０９の領域に対して従来のフリップ角を数度刻みで振る方法によりα'を求めると、この領域の緩和時間Ｔ１は次式で表される。
Ｔ１＝ｔ／ｌｎ（１−cos（α'））
この値を用いて上記と同様に１つまたは２つの画像からαを求めることができる。また上記２つの実施形態において、ＲＯＩに対してリファレンス領域を指定しているが、リファレンス領域を指定せず、画像全体をリファレンス領域として、ピクセル値を求めても良い。この場合は輪郭抽出を行い画像の平均値を求める必要がある。
【００５２】
本発明は、上述した実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々変形して実施することが可能である。さらに、上記実施形態には種々の段階が含まれており、開示される複数の構成要件における適宜な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されてもよい。
【００５３】
【発明の効果】
本発明により、１Ｈ−ＭＲＳにおいて水抑制パルスのフリップ角を事前に得られた画像から計算機内部で正確に得ることが可能となり、またフリップ角を求める調整時間を省くことが可能となる。
【図面の簡単な説明】
【図１】本発明の第１実施形態に係わる磁気共鳴映像装置の構成を示すブロック図。
【図２】本発明の第１の実施形態のフリップ角計算法で使用する３種の１Ｈ−ＭＲＳプラン画像を示す図。
【図３】本発明の第２の実施形態フリップ角計算法で使用する３種の  １Ｈ−ＭＲＳプラン画像を示す図。
【図４】ＭＲＳの一般的なパルスシーケンスを示す図。
【符号の説明】
１０１…頭部サジタル画像、
１０２…頭部アキシャル画像（Ｔ１強調画像）、
１０３…頭部アキシャル画像（Ｔ２強調画像）、
１０４…１Ｈ−ＭＲＳのＲＯＩ、
１０５…１Ｈ−ＭＲＳのＲＯＩ、
１０６…１Ｈ−ＭＲＳのＲＯＩ、
１０７…リファレンス領域、
１０８…リファレンス領域、
１０９…リファレンス領域、
２０１…静磁場磁石、
２０２…励磁用電源、
２０３…傾斜磁場コイル、
２０４…シールドコイル、
２０５…傾斜磁場コイル用電源、
２０６…被検体、
２０７…渦補償回路、
２０８…シムコイル、
２０９…シムコイル用電源、
２１０…送信部、
２１１…ＲＦコイル、
２１２…高周波シールド、
２１３…受信部、
２１４…データ収集部、
２１５…計算機、
２１６…シーケンスコントローラ、
２１７…入力装置、
２１８…ディスプレイ、
３０１…頭部サジタル画像、
３０２…頭部アキシャル画像（Ｔ１強調画像）、
３０３…頭部アキシヤル画像（Ｔ２強調画像）、
３０４…１Ｈ−ＭＲＳのＲＯＩ、
３０５…１Ｈ−ＭＲＳのＲＯＩ、
３０６…１Ｈ−ＭＲＳのＲＯＩ、
３０７…リファレンス領域、
３０８…リファレンス領域、
３０９…リファレンス領域、
３１０…プリサット領域、
３１１…プリサット領域、
３１２…プリサット領域、
３１３…プリサット領域、
３１４…プリサット領域、
３１５…プリサット領域。
【図１】
0004718698.tif 000011
【図２】
0004718698.tif 000012
【図３】
0004718698.tif 000013
【図４】
0004718698.tif 000014