JP2015013206A - Ｍｒｉ装置及びデータ処理装置 - Google Patents

Abstract

【課題】動脈瘤の破裂リスクを正確かつ容易に判定可能にするＭＲＩ装置を提供する。
【解決手段】被検体の動脈瘤を含む診断対象部位に対し所定のＭＲＩ撮影を行なって３次元形態データと３次元流速データを収集し、輪郭抽出部１５及び剪断応力計測部１６は、３次元形態データに基づいた動脈瘤の輪郭抽出と３次元流速データに基づいた動脈瘤内壁における剪断応力の計測を行なう。次いで、パラメータ算出部１７は、動脈瘤の輪郭情報及び剪断応力の計測結果に基づいて各種の診断パラメータを算出し、破裂リスク判定部１８及び高リスク部位検出部１９は、前記診断パラメータと所定閾値との比較により動脈瘤の破裂リスク判定と高リスク部位の位置検出を行なう。そして、表示部２０は、前記３次元形態データを用いて生成された３次元画像データに高リスク部位の位置情報と破裂リスクの判定結果を付加して自己のモニタに表示する。
【選択図】図１

Description

本発明は、ＭＲＩ装置及びデータ処理装置に係り、特に、動脈瘤や静脈瘤に対する破裂
リスクの判定を可能にするＭＲＩ装置及びデータ処理装置に関する。

磁気共鳴イメージング法（ＭＲＩ）は、静磁場中に置かれた物体の原子核スピンを、そ
のラーモア周波数をもつ高周波パルス（ＲＦパルス）で励起し、この励起に伴って発生す
る磁気共鳴信号（ＭＲ信号）に基づいて画像データを生成するイメージング法である。

ＭＲＩ装置は、生体内から検出されるＭＲ信号に基づいて画像データを生成する画像診
断装置であり、解剖学的診断情報のみならず生化学情報や機能診断情報など多くの情報を
得ることができるため、今日の画像診断の分野では不可欠なものとなっている。

このようなＭＲＩ装置を使用して血管内の狭窄部位や血管壁に生じた動脈瘤等の診断を
行なう場合、血管内を流れている血液を画像化することによって血管の形状を観測する方
法が一般に行なわれており、例えば、Gd-DTPA等の造影剤を注入することにより血液から
のＭＲ信号を高いコントラスト比で収集する造影ＭＲＡ撮影やＴＯＦ（time of flight）
法のパルスシーケンスを適用した非造影ＭＲＡ(magnetic resonance angiography)撮影等
が用いられている（例えば、特許文献１参照。）。

又、傾斜磁場中を流れる物体が発生したＭＲ信号の位相変化量はその移動速度に略比例
するという事実に基づき、血管内を流れる血液が発生したＭＲ信号の位相変化量を計測す
ることによりその速度情報を画像化するＰＣ（phase contrast）法も開発されている（例
えば、特許文献２参照。）。

そして、上述の造影ＭＲＡ撮影あるいは非造影ＭＲＡ撮影により３ｍｍ以上のサイズを
有した動脈瘤の検出が可能となり、特に、５ｍｍ以上の動脈瘤は早期治療の対象とされて
きた。

特開平０８−２０６０９１号公報 特開２００２−１６５７７１号公報

上述のように動脈瘤に対する治療は、ＭＲＩ撮影によって得られた形態画像データにお
いて５ｍｍ以上の大きさを有するものに対して行なわれてきた。しかしながら５ｍｍ以下
の動脈瘤であっても極めて高い破裂リスクを有している場合があり、従来のような形態情
報に基づいたＭＲＩ診断ではその破裂リスクを正確に判定することができないという問題
点を有していた。

本発明は、上述の問題点に鑑みてなされたものであり、その目的は、動脈瘤等の破裂リ
スクを正確かつ容易に判定することが可能なＭＲＩ装置及びデータ処理装置を提供するこ
とにある。

上記課題を解決するために、請求項１に係る本発明のＭＲＩ装置は、所定のパルスシー
ケンスを適用したＭＲＡ撮影により被検体の動脈瘤を含む診断対象部位が発生するＭＲ信
号を収集するＭＲ信号収集手段と、前記ＭＲ信号を処理して前記診断対象部位の３次元形
態データを生成する形態データ生成手段と、前記３次元形態データに基づいて前記動脈瘤
の輪郭を抽出し、輪郭データを生成する輪郭抽出手段と、前記ＭＲ信号を処理して前記診
断対象部位の３次元流速データを生成する流速データ生成手段と、前記輪郭データと前記
３次元流速データに基づいて前記動脈瘤における剪断応力を計測する剪断応力計測手段と
、前記剪断応力の計測結果あるいはこの剪断応力の計測結果と前記輪郭データに基づいて
前記動脈瘤の破裂リスクを判定する破裂リスク判定手段と、前記破裂リスクの判定結果を
表示する表示手段とを備えたことを特徴としている。

又、請求項２に係る本発明のデータ処理装置は、ＭＲＩ装置によって収集された動脈瘤
を含む診断対象部位の３次元形態データ及び３次元流速データを保管するデータ保管手段
と、前記３次元形態データに基づいて前記動脈瘤の輪郭を抽出し、輪郭データを生成する
輪郭抽出手段と、前記輪郭データと前記３次元流速データに基づいて前記動脈瘤における
剪断応力を計測する剪断応力計測手段と、前記剪断応力の計測結果あるいはこの剪断応力
の計測結果と前記輪郭データに基づいて前記動脈瘤の破裂リスクを判定する破裂リスク判
定手段と、前記破裂リスクの判定結果を表示する表示手段とを備えたことを特徴としてい
る。

本発明によれば、動脈瘤等の破裂リスクを正確かつ容易に判定することができる。この
ため、破裂リスクの高い動脈瘤等を早期に治療することが可能となる。

本発明の第１の実施例におけるＭＲＩ装置の全体構成を示すブロック図。 同実施例のＭＲＩ装置が備えるＭＲ信号収集部の具体的な構成を示すブロック図。 同実施例における剪断応力の算出方法を示す図。 同実施例のＭＲＩ装置が備えるパラメータ算出部の機能ブロック図。 同実施例における動脈瘤の破裂リスク判定基準と破裂リスクのレベルを説明するための図。 同実施例における破裂リスク判定結果の表示手順を示すフローチャート。 本発明の第２の実施例におけるデータ解析装置の全体構成を示すブロック図。

以下、図面を参照して本発明の実施例を説明する。

以下に述べる第１の実施例のＭＲＩ装置は、先ず、被検体の動脈瘤を含む診断対象部位
（脳動脈）に対し所定のパルスシーケンスを適用したＭＲＩ撮影を行なって３次元形態デ
ータと３次元流速データを収集し、３次元形態データに基づいた動脈瘤の輪郭抽出と３次
元流速データに基づいた動脈瘤内壁における剪断応力の計測を行なう。次いで、動脈瘤の
輪郭情報及び剪断応力の計測結果に基づいて各種の診断パラメータを算出し、これらの診
断パラメータと所定の閾値との比較によって動脈瘤の破裂リスク判定と高リスク部位の位
置検出を行なう。そして、前記３次元形態データを処理して生成した３次元画像データに
高リスク部位の位置情報を重畳し、更に、破裂リスクの判定結果を付加して表示部に表示
する。

尚、以下の実施例では、ＴＯＦ法のパルスシーケンスを適用した非造影ＭＲＡ撮影によ
って動脈瘤を含む血管の３次元形態データを収集し、ＰＣ法のパルスシーケンスを適用し
た非造影ＭＲＡ撮影によって動脈瘤内の３次元流速データを収集する場合について述べる
が、他の方法によって３次元形態データや３次元流速データを収集してもよい。例えば、
血管の３次元形態データは、Gd-DTPA等の造影剤を注入した状態で行なわれる造影ＭＲＡ
撮影や磁化の定常状態を利用したＴrueＳＳＦＰ(true steady-state free precession)法
あるいは特開２０００−５１４４号公報等に記載されているＦＢＩ（fresh blood imagin
g）法のパルスシーケンスを適用した非造影ＭＲＡ撮影によって収集してもよい。

又、頭部の脳動脈を診断対象部位とし、この脳動脈に発生した動脈瘤に対して破裂リス
クの判定を行なう場合について述べるが、他の部位に発生した動脈瘤あるいは静脈瘤に対
する破裂リスクの判定であっても構わない。

（装置の構成）
本発明の第１の実施例におけるＭＲＩ装置の構成につき図１乃至図５を用いて説明する
。尚、図１は、本実施例におけるＭＲＩ装置の全体構成を示すブロック図であり、図２及
び図４は、このＭＲＩ装置が備えるＭＲ信号収集部及びパラメータ算出部の具体的な構成
を示すブロック図である。

図１に示すＭＲＩ装置１００は、被検体の動脈瘤を含む頭部の診断対象部位（脳動脈）
に対し静磁場及び各種の傾斜磁場を印加し、このとき前記診断対象部位から発生するＭＲ
信号を収集するＭＲ信号収集部１０と、ＴＯＦ法のパルスシーケンスを適用した傾斜磁場
を用いて得られる形態計測モードのＭＲ信号を処理して前記診断対象部位の３次元形態デ
ータを生成する形態データ生成部１１と、ＰＣ法のパルスシーケンスを適用した傾斜磁場
を用いて得られる血流速度計測モードのＭＲ信号を処理することにより３次元流速データ
を生成する流速データ生成部１２と、得られた３次元形態データ及び３次元流速データを
保存するデータ記憶部１３と、このデータ記憶部１３に一旦保存された３次元形態データ
をレンダリング処理して３次元画像データを生成する画像データ生成部１４を備えている
。

又、ＭＲＩ装置１００は、３次元形態データにおける動脈瘤の輪郭を抽出する輪郭抽出
部１５と、データ記憶部１３に保存された３次元流速データ及び３次元形態データに基づ
いて動脈瘤の内壁における剪断応力を計測する剪断応力計測部１６と、剪断応力の計測結
果及び動脈瘤の輪郭情報に基づいて各種の診断パラメータを算出するパラメータ算出部１
７と、これらの診断パラメータと予め設定された閾値との比較により動脈瘤に対する破裂
リスクを判定する破裂リスク判定部１８と、上述の診断パラメータに基づいて前記動脈瘤
における高リスク部位を検出する高リスク部位検出部１９を備え、更に、画像データ生成
部１４が生成した診断対象部位の３次元画像データに高リスク部位検出部１９が検出した
高リスク部位の位置情報と破裂リスク判定部１８による破裂リスクの判定結果を付加して
表示する表示部２０と、被検体情報の入力、ＭＲＩ撮影条件の設定、各種コマンド信号の
入力等を行なう入力部２１と、被検体の心拍時相を計測する心拍時相計測部２２と、ＭＲ
Ｉ装置１００が備えた上述の各ユニットを統括的に制御するシステム制御部２３を備えて
いる。

次に、上述のＭＲ信号収集部１０の具体的な構成につき図２のブロック図を用いて説明
する。

図２に示すＭＲ信号収集部１０は、被検体１５０の診断対象部位に対して静磁場を発生
する静磁場発生部１と、ＴＯＦ法のパルスシーケンスを適用した形態計測モードの傾斜磁
場及びＰＣ法のパルスシーケンスを適用した血流速度計測モードの傾斜磁場を発生する傾
斜磁場発生部２と、被検体１５０に対してＲＦパルスの照射とＭＲ信号の受信を行なう送
受信部３と、形態計測モード及び血流速度計測モードにおける上述のパルスシーケンスを
制御するシーケンス制御部４と、被検体１５０を載置する天板５とを備えている。

静磁場発生部１は、常伝導磁石あるいは超電導磁石等によって構成される主磁石１０１
と、この主磁石１０１を駆動するための静磁場電源１０２を備え、静磁場電源１０２は、
主磁石１０１に対して所定の電流を供給することにより図示しないガントリの撮影野に配
置された被検体１５０の診断対象部位に対して強力な静磁場を形成する。尚、上述の主磁
石１０１は、永久磁石によって構成されていてもよい。

一方、傾斜磁場発生部２は、互いに直交するＸ軸方向、Ｙ軸方向及びＺ軸方向に対して
傾斜磁場を形成する複数の傾斜磁場コイル２０１と、傾斜磁場コイル２０１の各々に対し
てパルス電流を供給する傾斜磁場電源２０２を備えている。

傾斜磁場電源２０２は、被検体１５０の診断対象部位が配置されたガントリの撮影野に
対して符号化を行なう。即ち、傾斜磁場電源２０２は、シーケンス制御部４から供給され
る形態計測モード及び血流速度計測モードのシーケンス制御信号に基づいてＸ軸方向，Ｙ
軸方向及びＺ軸方向の傾斜磁場コイル２０１に供給するパルス電流を制御することにより
各々の方向に対して傾斜磁場を形成する。そして、Ｘ軸方向，Ｙ軸方向及びＺ軸方向の傾
斜磁場は合成されて互いに直交するスライス選択傾斜磁場、位相エンコード傾斜磁場及び
読み出し（周波数エンコード）傾斜磁場が所望の方向に形成される。

又、ＰＣ法が適用された本実施例の血流速度計測モードでは、移動体からのＭＲ信号に
対して位相変化を発生させるためのフローエンコード傾斜磁場が上述のスライス選択傾斜
磁場、位相エンコード傾斜磁場及び読み出し傾斜磁場の方向に対して更に印加される。即
ち、スライス選択傾斜磁場、位相エンコード傾斜磁場、読み出し傾斜磁場及びフローエン
コード傾斜磁場は、主磁石１０１によって形成された静磁場に重畳されて被検体１５０の
診断対象部位に印加される。

次に、送受信部３は、被検体１５０に対してＲＦパルスを照射する送信コイル３０１及
び送信部３０２と、被検体１５０にて発生したＭＲ信号を受信する受信コイル３０３及び
受信部３０４を有している。但し、送信コイル３０１の機能と受信コイル３０３の機能を
１つのコイルで兼ね備えた送受信コイルを用いてもよい。

送信部３０２は、シーケンス制御部４から供給されるシーケンス制御信号に基づき、主
磁石１０１の静磁場強度によって決定される水素原子核スピンの磁気共鳴周波数（ラーモ
ア周波数）と同じ周波数の搬送波を有し所定の選択励起波形で変調されたパルス電流を生
成する。即ち、送信部３０２は、例えば、水組織の磁気共鳴周波数を中心周波数とし所定
の周波数帯域を有するＲＦパルスを発生させるためのパルス電流を生成する。

送信コイル３０１は、送信部３０２から供給されるパルス電流によって駆動され、被検
体１５０の診断対象部位に対してＲＦパルスを照射する。一方、受信コイル３０３は、前
記ＲＦパルスの照射によって被検体１５０の診断対象部位における組織から発生したＭＲ
信号を検出する。又、ＴＯＦ法が適用された本実施例の形態計測モードでは、撮影スライ
ス断面に隣接した領域に対する前飽和パルス（pre-saturation pulse）の照射がＭＲ信号
の収集に先立って行なわれる。尚、受信コイル３０３は、通常、ＭＲ信号を高感度で検出
するために小口径のコイルが複数個（Ｎ個）配列された所謂アレイコイルによって構成さ
れている。

受信部３０４は、図示しないＮチャンネルの増幅回路、中間周波変換回路、検波回路、
フィルタリング回路及びＡ／Ｄ変換器を備え、受信コイル３０３が検出した微小なＭＲ信
号に対し増幅、中間周波変換、位相検波、フィルタリング等の信号処理を行なった後Ａ／
Ｄ変換する。但し、前記増幅回路は、受信コイル３０３が検出したＭＲ信号を高いＳ／Ｎ
で増幅するために受信コイル３０３の近傍に設けられている。

尚、本実施例の形態計測モードでは、例えば、特許文献１に記載されているようなＴＯ
Ｆ法のパルスシーケンスに基づいた傾斜磁場の印加とＲＦパルスの照射によってＭＲ信号
の収集を行ない、血流速度計測モードでは、例えば、特許文献２に記載されているような
ＰＣ法のパルスシーケンスに基づいた傾斜磁場の印加とＲＦパルスの照射によってＭＲ信
号の収集を行なう。

そして、上述の主磁石１０１、傾斜磁場コイル２０１、送信コイル３０１及び受信コイ
ル３０３は、ＭＲＩ装置１００の図示しないガントリに設けられ、このガントリの中央部
には撮影野が形成される。即ち、ガントリの中心部には天板４と共に被検体１５０が挿入
される撮影野が設けられ、この撮影野の周囲には受信コイル３０３、送信コイル３０１、
傾斜磁場コイル２０１及び主磁石１０１がＺ軸（体軸）を共軸として同心円状に配置され
ている。

次に、シーケンス制御部４は、図示しないＣＰＵを備え、システム制御部２３から供給
されるシーケンス情報に基づいてＴＯＦ法が適用された形態計測モードのシーケンス制御
信号及びＰＣ法が適用された血流速度計測モードのシーケンス制御信号を生成する。そし
て、これらのシーケンス制御信号を傾斜磁場発生部２の傾斜磁場電源２０２及び送受信部
３の送信部３０２へ供給することにより、傾斜磁場コイル２０１及び送信コイル３０１に
供給されるパルス電流を制御する。

一方、天板５は、ガントリの近傍に設置された図示しない寝台の上面においてＺ軸方向
にスライド自在に取り付けられ、天板５に載置された被検体１５０を体軸方向（Ｚ軸方向
）に移動することにより被検体１５０の診断対象部位を撮影野の所望位置に配置する。こ
の場合、診断対象部位が撮影野の近傍に設けられた受信コイル３０３に対向するような天
板５の移動が図示しない天板移動機構部及び天板移動制御部によって行なわれる。

図１へ戻って、形態データ生成部１１は、図示しないＭＲ信号記憶部と高速演算部を備
えている。そして、ＴＯＦ法が適用された形態計測モードのパルスシーケンスに基づいて
診断対象部位から収集されたＭＲ信号は、受信部３０４によって中間周波変換、位相検波
、更には、Ａ／Ｄ変換された後シーケンス制御部４から供給される撮像位置情報を付帯情
報として前記ＭＲ信号記憶部に保存される。一方、前記高速演算部は、前記ＭＲ信号記憶
部に一旦保存されたＭＲ信号及び撮像位置情報を読み出し、フーリエ変換による再構成処
理を行なって３次元形態データを生成する。

同様にして、流速データ生成部１２は、図示しないＭＲ信号記憶部と高速演算部を備え
、ＰＣ法が適用された血流速度計測モードのパルスシーケンスに基づいて前記診断対象部
位から収集されたＭＲ信号は、受信部３０４において上述の信号処理が行なわれた後撮像
位置情報を付帯情報として前記ＭＲ信号記憶部に保存される。そして、前記高速演算部は
、前記ＭＲ信号記憶部に一旦保存されたＭＲ信号及び撮像位置情報を読み出し、フーリエ
変換による再構成処理を行なうことにより動脈瘤を含む診断対象部位の内部を流れる血流
の３次元流速データを生成する。

尚、ＰＣ法では、血流のような移動体から得られるＭＲ信号に対して位相変化を生じさ
せるために、互いに逆方向の勾配を有した１対のフローエンコード傾斜磁場をスライス選
択傾斜磁場方向、位相エンコード傾斜磁場方向及び読み出し傾斜磁場方向の各々に対して
印加し、順方向のフローエンコード傾斜磁場を印加して得られたＭＲ信号と逆方向のフロ
ーエンコード傾斜磁場を印加して得られたＭＲ信号との差分処理によって新たに得られる
差分ＭＲ信号を再構成処理することによりフローエンコード傾斜磁場方向に対する血流の
流速成分を得ることができる。

この場合、１対からなる上述のフローエンコード傾斜磁場をスライス選択傾斜磁場方向
、位相エンコード傾斜磁場方向及び読み出し傾斜磁場方向の各々に対して順次印加して得
られたＭＲ信号を差分処理する。そして、得られた差分ＭＲ信号を再構成処理することに
より３次元形態データの各ボクセルに対応したスライス選択傾斜磁場方向の流速成分Ｉｓ
、位相エンコード傾斜磁場方向の流速成分Ｉｐ及び読み出し傾斜磁場方向の流速成分Ｉｒ
を算出し、これらの流速成分に基づいて３次元流速データを生成する。即ち、前記３次元
流速データは、３次元形態データの各ボクセルに対応する複数の位置にて算出されたスラ
イス選択傾斜磁場方向の流速成分Ｉｓ、位相エンコード傾斜磁場方向の流速成分Ｉｐ及び
読み出し傾斜磁場方向の流速成分Ｉｒによって構成される。

データ記憶部１３は、形態データ記憶領域と流速データ記憶領域を有し、上述の形態デ
ータ生成部１１から供給される３次元形態データは前記形態データ記憶領域に、又、流速
データ生成部１２から供給される３次元流速データは前記流速データ記憶領域に夫々保存
される。この場合、心拍時相計測部２２において計測された当該被検体の心拍時相情報も
３次元形態データ及び３次元流速データの付帯情報として上述の記憶領域に保存される。

一方、画像データ生成部１４は、図示しない不透明度・色調設定部とレンダリング処理
部を備えている。前記不透明度・色調設定部は、データ記憶部１３の形態データ記憶領域
に保存された３次元形態データを読み出し、この３次元形態データのボクセル値に基づい
て不透明度や色調をボクセル単位で設定する。一方、前記レンダリング処理部は、前記不
透明度・色調設定部が設定した不透明度や色調の情報に基づいて上述の３次元形態データ
をレンダリング処理し、ボリュームレンダリング画像データ等の３次元画像データを生成
する。

輪郭抽出部１５は、図示しない演算処理部を備え、例えば、領域拡張（region growing
）法を適用して動脈瘤内壁の輪郭を抽出する。即ち、前記演算処理部は、データ記憶部１
３の形態データ記憶領域に保存された３次元形態データを読み出し、更に、表示部２０に
表示された前記３次元画像データの動脈瘤に対して入力部２１が設定する中心位置の情報
を受信する。そして、３次元形態データに対し動脈瘤の中心位置を基準とした領域拡張法
を適用させて動脈瘤の輪郭を抽出する。例えば、前記演算処理部は、動脈瘤の中心位置を
基準とする所定サイズの３次元関心領域を３次元形態データに対して設定し、この３次元
関心領域を順次拡大させながらその表面と交叉する３次元形態データのボクセル値と所定
の閾値γとを比較する。そして、閾値γより大きなボクセル値を有する３次元形態データ
のボクセルを動脈瘤の内壁として認識することによりその輪郭を抽出する。

尚、ここでは、領域拡張法を用いて動脈瘤の輪郭を抽出する場合について述べたが、輪
郭抽出法は上述の領域拡張法に限定されるものではなく、例えば、marching-cube法のよ
うな他の輪郭抽出法を用いて動脈瘤の輪郭を抽出してもよい。

次に、剪断応力計測部１６は、図示しない流速値算出部と剪断応力算出部を備えている
。前記流速データ算出部は、先ず、上述の３次元形態データのボクセルＢ（ｘ、ｙ、ｚ）
に対応した３次元流速データを構成するスライス選択傾斜磁場方向の流速成分Ｉｓ（ｘ、
ｙ、ｚ）、位相エンコード傾斜磁場方向の流速成分Ｉｐ（ｘ、ｙ、ｚ）及び読み出し傾斜
磁場方向の流速成分Ｉｒ（ｘ、ｙ、ｚ）をデータ記憶部１３の流速データ記憶領域から読
み出す。次いで、これらの流速成分Ｉｓ（ｘ、ｙ、ｚ）、Ｉｐ（ｘ、ｙ、ｚ）及びＩｒ（
ｘ、ｙ、ｚ）を下式（１）へ代入して３次元流速データの絶対値（以下では、３次元流速
値と呼ぶ。）Ｉｏ（ｘ、ｙ、ｚ）を算出する。

一方、前記剪断応力算出部は、動脈瘤の内壁に沿って流れる血液の流速値を上述の３次
元流速値Ｉｏ（ｘ、ｙ、ｚ）に基づいて計測し、更に、この内壁に対して垂直な方向にお
ける前記流速値の変化率に基づいて動脈瘤の内壁における剪断応力を算出する。図３は剪
断応力の算出方法を示したものであり、動脈流の内壁に与える剪断応力τｗは次式（２）
によって表すことができる。

ここで、αは血液の粘性係数、ｕは動脈瘤の内壁に沿った血液の流速値、ｈは前記内壁
からの距離を示している。

次に、図１に示したパラメータ算出部１７の具体例につき図４の機能ブロック図を用い
て説明する。このパラメータ算出部１７は、アスペクト比、平均剪断応力、最小剪断応力
及び最大剪断応力を診断パラメータとして算出あるいは抽出する機能を有し、図４に示す
ようにアスペクト比算出部１７１、平均剪断応力算出部１７２、最小剪断応力抽出部１７
３及び最大剪断応力抽出部１７４を備えている。

アスペクト比算出部１７１は、輪郭抽出部１５から供給される動脈瘤の輪郭データに基
づき、動脈瘤の最大内径と最小内径との比（アスペクト比）を算出する。一方、平均剪断
応力算出部１７２は、輪郭抽出部１５から供給される前記動脈瘤の輪郭データと剪断応力
計測部１６から供給される剪断応力の算出結果を受信し、動脈瘤の輪郭データに囲まれた
３次元領域（即ち、動脈瘤の内部）にて得られた複数からなる剪断応力の算出結果を加算
平均して平均剪断応力を算出する。又、最小剪断応力抽出部１７３は、前記複数からなる
剪断応力の算出結果の中から最小剪断応力を抽出し、最大剪断応力抽出部１７４は、前記
複数からなる剪断応力の算出結果の中から最大剪断応力を抽出する。

尚、上述の最小剪断応力は、通常、図３に示すように動脈瘤の血液流入口Ｐａに対向し
た頂点Ｐｂの近傍において発生し、この領域における最小剪断応力が小さい程他の領域（
動脈瘤の内壁）において大きな剪断応力が発生することが知られている。

再び図１へ戻って、破裂リスク判定部１８は、パラメータ算出部１７によって算出ある
いは抽出されたアスペクト比、平均剪断応力及び最小剪断応力と所定の閾値とを比較する
ことにより動脈瘤に対する破裂リスクを判定する。図５（ａ）は、破裂リスクの判定に用
いる判定基準１乃至判定基準３を示したものであり、これらの判定基準の何れかに該当す
る動脈瘤は高リスク候補として判定される。即ち、判定基準１に示すように動脈瘤のアス
ペクト比が閾値β１（β１＝１．６）以上である場合、判定基準２に示すように心臓収縮
末期（peak systole）の動脈瘤内部にて算出される平均剪断応力が閾値β２（β２＝２．
９６Ｎ／ｍ２：低リスク動脈瘤が有する平均剪断応力１．４８Ｎ／ｍ２の２倍）以上であ
る場合、あるいは、判定基準３に示すように動脈瘤内壁の頂点近傍にて抽出される最小剪
断応力が閾値β３（β３＝０．５Ｎ／ｍ２）以下である場合、この動脈瘤は高リスク候補
として判定される。

そして、図５（ｂ）に示すように、上述のアスペクト比、平均剪断応力及び最小剪断応
力と閾値β１乃至閾値β３との比較において判定基準１乃至判定基準３の何れにも該当し
ない場合、この動脈瘤の破裂リスクは低レベル（即ち、破裂の可能性は極めて小さいレベ
ル）にあると判定される。

一方、判定基準１乃至判定基準３の何れか１つに該当する場合の破裂リスクは中レベル
にあると判定され、２つ以上の判定基準に該当する場合の破裂リスクは高レベルにあると
判定される。

図１へ戻って、高リスク部位検出部１９は、動脈瘤の内壁にて得られた複数からなる剪
断応力の算出結果の中からパラメータ算出部１７の最大剪断応力抽出部１７４が抽出した
最大剪断応力の位置情報に基づいて動脈瘤における高リスク部位を検出する。

一方、表示部２０は、図示しない表示データ生成部、データ変換部及びモニタを備え、
前記表示データ生成部は、画像データ生成部１４にて生成された当該診断対象部位の３次
元画像データに高リスク部位検出部１９が検出した最大剪断応力を有する高リスク部位の
位置情報を重畳し、更に、破裂リスク判定部１８によって判定された破裂リスクのレベル
情報や入力部２１からシステム制御部２３を介して供給される当該被検体の被検体情報等
を付加して表示データを生成する。そして、前記データ変換部は、得られた表示データを
所定の表示フォーマットに変換し前記モニタに表示する。

尚、上述の３次元画像データに重畳される高リスク部位の位置情報は、破裂リスク判定
部１８による破裂リスクの低レベル、中レベル及び高レベルに対応させてカラー表示され
る。カラー表示された上述の位置情報を観測することにより、操作者は、破裂リスクが高
い部位とその程度を容易に把握することが可能となる。

次に、入力部２１は、被検体情報や個人健康情報の入力、計測モードの選択、３次元形
態データ収集条件及び３次元流速データ収集条件の設定、３次元画像データ生成条件の設
定、高リスク候補の判定に用いる閾値β１乃至閾値β３の設定、輪郭抽出に用いる閾値γ
の設定、３次元画像データの動脈瘤に対する中心位置の設定、表示データ生成条件の設定
、更には、各種コマンド信号の入力等を行なう。尚、上述の被検体情報として、被検体の
氏名、ＩＤ、性別、年齢等があり、個人健康情報として、被検体の血圧、喫煙暦及び病歴
や家族の病歴等がある。

心拍時相計測部２２は、被検体の心電波形に基づいて心拍時相を設定する機能を有し、
被検体体表面に装着され心電波形を計測する計測用電極と、この計測用電極によって計測
された心電波形を所定の振幅に増幅する増幅回路と、増幅された心電波形をデジタル信号
に変換するＡ／Ｄ変換器と、デジタル変換された心電波形に基づいて、例えば、心臓の拡
張末期を基準とする１心拍周期に対して心拍時相を設定する心拍時相設定部（何れも図示
せず）を備えている。

システム制御部２３は、図示しないＣＰＵと記憶回路を備え、ＭＲＩ装置１００を統括
して制御する機能を有している。前記記憶回路には、形態計測モードのＴＯＦ法及び血流
速度計測モードのＰＣ法に適用可能なパルスシーケンスデータが保管され、更に、入力部
２１にて入力／設定／選択された各種の情報が保存される。そして、前記ＣＰＵは、前記
記憶回路に保存された上述の情報に基づいてＭＲＩ装置１００の各ユニットを制御し、当
該被検体１５０の動脈瘤に対する破裂リスクの判定を行なう。

（破裂リスク判定結果の表示手順）
次に、本実施例における破裂リスク判定結果の表示手順につき図６のフローチャートを
用いて説明する。

被検体１５０のＭＲＩ撮影に先立ちＭＲＩ装置１００の操作者は、天板５に載置した被
検体１５０をＺ軸方向へ移動することによりその診断対象部位（脳動脈）をガントリの撮
影野に配置し、入力部２１において被検体情報や個人健康情報の入力、３次元形態データ
収集条件及び３次元流速データ収集条件の設定、３次元画像データ生成条件の設定、高リ
スク候補の判定に用いる閾値β１乃至閾値β３の設定、表示データ生成条件の設定等を行
なう（図６のステップＳ１）。

次いで、操作者は、入力部２１において形態計測モードを選択し撮影開始コマンドを入
力する。そして、このコマンド信号を受信したシステム制御部２３は、ＭＲ信号収集部１
０のシーケンス制御部４に対して形態計測モードに適用されるＴＯＦ法のシーケンス情報
を供給し、シーケンス制御部４は、このシーケンス情報に基づいて生成したシーケンス制
御信号を傾斜磁場発生部２の傾斜磁場電源２０２及び送受信部３の送信部３０２へ供給す
ることにより、当該診断対象部位に対するＭＲ信号の収集を、例えば、１心拍周期におい
て行なう。

一方、形態データ生成部１１は、自己のＭＲ信号記憶部に一旦保存された形態計測モー
ドのＭＲ信号を再構成処理して３次元形態データを生成し、得られた１心拍周期分の時系
列的な３次元形態データは心拍時相計測部２２から供給される当該被検体の心拍時相情報
と共にデータ記憶部１３の形態データ記憶領域に保存される。（図６のステップＳ２）。

診断対象部位に対する３次元形態データの生成と保存が終了したならば、操作者は、入
力部２１において血流速度計測モードを選択し撮影開始コマンドを再度入力する。そして
、このコマンド信号を受信したシステム制御部２３は、ＭＲ信号収集部１０のシーケンス
制御部４に対して血流速度計測モードに適用されるＰＣ法のシーケンス情報を供給し、シ
ーケンス制御部４は、このシーケンス情報に基づいて生成したシーケンス制御信号を傾斜
磁場発生部２の傾斜磁場電源２０２及び送受信部３の送信部３０２へ供給することにより
１心拍周期におけるＭＲ信号の収集を行なう。そして、流速データ生成部１２は、自己の
ＭＲ信号記憶部に一旦保存された血流速度計測モードのＭＲ信号を再構成処理して３次元
流速データを生成し、得られた１心拍周期分の時系列的な３次元流速データは当該被検体
の心拍時相情報と共にデータ記憶部１３の流速データ記憶領域に保存される。（図６のス
テップＳ３）。

診断対象部位に対する時系列的な３次元形態データ及び３次元流速データの収集が終了
したならば、画像データ生成部１４は、データ記憶部１３の形態データ記憶領域から読み
出した所定心拍時相の３次元形態データをレンダリング処理して３次元画像データを生成
し、表示部２０のモニタに表示する。

表示部２０に表示された３次元画像データを観測した操作者は、この３次元画像データ
に示された動脈瘤の中心位置を入力部２１に設けられた入力デバイスを用いて設定する。
次いで、輪郭抽出部１５は、データ記憶部１３の形態データ記憶領域に保存されている３
次元形態データを読み出し、更に、入力部２１からシステム制御部２３を介して供給され
る動脈瘤の中心位置情報を受信する。そして、３次元形態データに対し動脈瘤の中心位置
を基準とする領域拡張法を適用させて動脈瘤の輪郭を抽出する（図６のステップＳ４）
一方、剪断応力計測部１６は、所定心拍時相における３次元形態データのボクセルに対
応した３次元流速データをデータ記憶部１３の流速データ記憶領域から読み出し、この３
次元流速データを構成するスライス選択傾斜磁場方向、位相エンコード傾斜磁場方向及び
読み出し傾斜磁場方向の流速成分に基づいて３次元流速値を算出する。そして、この３次
元流速値に基づき動脈瘤の内壁に沿って流れる血液の流速値を計測し、更に、この内壁に
対して垂直な方向における前記流速値の変化率に基づいて動脈瘤の内壁における剪断応力
を算出する（図６のステップＳ５）。

次に、パラメータ算出部１７は、輪郭抽出部１５から供給される動脈瘤の輪郭データに
基づいて動脈瘤のアスペクト比を算出し、前記輪郭データに囲まれた領域にて得られた複
数からなる剪断応力の算出結果を加算平均して平均剪断応力を算出する。更に、前記複数
からなる剪断応力の算出結果の中から最小剪断応力及び最大剪断応力を抽出する（図６の
ステップＳ６）。

そして、破裂リスク判定部１８は、パラメータ算出部１７によって算出あるいは抽出さ
れたアスペクト比、平均剪断応力及び最小剪断応力と所定の閾値とを比較することにより
動脈瘤に対する破裂リスクを判定し（図６のステップＳ７）、高リスク部位検出部１９は
、パラメータ算出部１７が複数からなる剪断応力の算出結果の中から抽出した最大剪断応
力の位置情報に基づいて動脈瘤における高リスク部位を検出する（図６のステップＳ８）
。

一方、画像データ生成部１４は、データ記憶部１３の形態データ記憶領域から読み出し
た１心拍周期分の３次元形態データをレンダリング処理して時系列的な３次元画像データ
を生成する（図６のステップＳ９）。

そして、表示部２０は、画像データ生成部１４にて生成された時系列的な３次元画像デ
ータの各々に高リスク部位検出部１９が検出した最大剪断応力を有する高リスク部位の位
置情報を重畳し、更に、破裂リスク判定部１８によって判定された破裂リスクのレベル情
報や入力部２１からシステム制御部２３を介して供給される当該被検体の被検体情報等を
付加して表示データを生成し自己のモニタに表示する（図６のステップＳ１０）。

以上述べた本発明の実施例によれば、動脈瘤に対するＭＲＩ撮影によって収集した３次
元形態データ及び３次元流速データに基づいて破裂リスクの判定を行なっているため、動
脈瘤の破裂リスクを正確かつ容易に判定することができる。このため、破裂リスクの高い
動脈瘤を早期に治療することが可能となる。

特に、上述の実施例では、３次元流速データの流速分布情報を用いて動脈瘤内壁に発生
する剪断応力を計測し、この剪断応力に基づいて破裂リスクの判定を行なっているため、
従来のような形態情報からは得られなかった動脈瘤の機能情報を正確に得ることができる
。又、最大剪断応力の発生部位を検出することにより高い破裂リスクを有する動脈瘤の内
壁を特定することが可能となる。

更に、上述の３次元形態データ及び３次元流速データは、ＴＯＦ法及びＰＣ法のパルス
シーケンスを適用した非造影ＭＲＡ撮影によって収集しているため、非侵襲的な検査が可
能となり被検体に対する負担を軽減することができる。

次に、本発明の第２の実施例におけるデータ解析装置について説明する。このデータ解
析装置は、先ず、別途設置されたＭＲＩ装置によって収集され、ネットワーク等を介して
供給された当該被検体の動脈瘤を含む診断対象部位（脳動脈）の３次元形態データ及び３
次元流速データを一旦保存する。次いで、３次元形態データに基づいた動脈瘤の輪郭抽出
と３次元流速データに基づいた動脈瘤内壁における剪断応力の計測を行ない、動脈瘤の輪
郭情報及び剪断応力の計測結果に基づいて各種の診断パラメータを算出する。そして、こ
れらの診断パラメータと所定の閾値との比較によって動脈瘤の破裂リスク判定と高リスク
部位の位置検出を行ない、前記３次元形態データを処理して生成した３次元画像データに
高リスク部位の位置情報と破裂リスクの判定結果を付加して表示部に表示する。

（装置の構成）
本発明の第２の実施例におけるデータ解析装置の構成につき図７を用いて説明する。尚
、図７は、前記データ解析装置の全体構成を示すブロック図であり、図１に示した第１の
実施例におけるＭＲＩ装置１００のユニットと同一の構成及び機能を有するユニットは同
一の符号を付加し詳細な説明は省略する。

即ち、図７に示すデータ解析装置２００は、当該被検体の動脈瘤を含む診断対象部位（
脳動脈）から予め収集された３次元形態データ及び３次元流速データを保管するデータ保
管部２４と、このデータ保管部２４に保管された３次元形態データをレンダリング処理し
て３次元画像データを生成する画像データ生成部１４と、３次元形態データにおける動脈
瘤の輪郭を抽出する輪郭抽出部１５と、データ記憶部１３に保存された３次元流速データ
及び３次元形態データに基づいて動脈瘤の内壁における剪断応力を計測する剪断応力計測
部１６と、剪断応力の計測結果及び動脈瘤の輪郭情報に基づいて各種の診断パラメータを
算出するパラメータ算出部１７を備えている。

更に、データ解析装置２００は、パラメータ算出部１７において算出された診断パラメ
ータと予め設定された閾値との比較により動脈瘤に対する破裂リスクを判定する破裂リス
ク判定部１８と、上述の診断パラメータに基づいて前記動脈瘤における高リスク部位を検
出する高リスク部位検出部１９と、画像データ生成部１４が生成した診断対象部位の３次
元画像データに高リスク部位検出部１９が検出した高リスク部位の位置情報を重畳し、更
に、破裂リスク判定部１８による破裂リスクの判定結果を付加して表示する表示部２０と
、被検体情報の入力や各種コマンド信号の入力等を行なう入力部２１ａと、データ解析装
置２００が備える上述の各ユニットを統括的に制御するシステム制御部２３ａを備えてい
る。

データ解析装置２００は、別途設置されたＭＲＩ装置によって予め収集され、ネットワ
ーク２５あるいは図示しない大容量の記憶媒体等を介して供給された当該被検体の診断対
象部位に対する３次元形態データ及び３次元流速データを保管する。この場合、３次元形
態データ及び３次元流速データの収集時における当該被検体の心拍時相情報も付帯情報と
して保存される。

一方、入力部２１ａは、３次元画像データ生成条件の設定、高リスク候補の判定に用い
る閾値β１乃至閾値β３の設定、輪郭抽出に用いる閾値γの設定、３次元画像データの動
脈瘤に対する中心位置の設定、表示データ生成条件の設定、更には、各種コマンド信号の
入力等を行なう。

システム制御部２３ａは、図示しないＣＰＵと記憶回路を備え、データ解析装置２００
を統括して制御する機能を有している。前記記憶回路には、入力部２１ａにて入力あるい
は設定された各種の情報が保存される。そして、前記ＣＰＵは、前記記憶回路に保存され
た上述の入力情報や設定情報に基づいてデータ解析装置２００の各ユニットを統括的に制
御し、当該被検体１５０の動脈瘤に対する破裂リスクの判定を行なう。

尚、本実施例における破裂リスク判定結果の表示手順は、図６のフローチャートに示し
たステップＳ４乃至ステップ１０の手順と同様であるため説明は省略する。

以上述べた本発明の実施例によれば、動脈瘤に対するＭＲＩ撮影によって予め収集され
た３次元形態データ及び３次元流速データに基づいて破裂リスクの判定を行なっているた
め、動脈瘤の破裂リスクを正確かつ容易に判定することができる。このため、破裂リスク
の高い動脈瘤を早期に治療することが可能となる。

特に、３次元流速データの流速分布情報を用いて動脈瘤内壁に発生する剪断応力を計測
し、この剪断応力に基づいて破裂リスクの判定を行なっているため、従来のような形態情
報からは得られなかった動脈瘤の機能情報を正確に得ることができる。又、最大剪断応力
の発生部位を検出することにより高い破裂リスクを有する動脈瘤の内壁を特定することが
可能となる。

更に、別途設置されたＭＲＩ装置からネットワークあるいは記憶媒体等を介して供給さ
れる３次元形態データ及び３次元流速データを用いて動脈瘤の破裂リスクを判定すること
ができるため、操作者は、時間や場所の制約をあまり受けることなく当該被検体の動脈瘤
に対する診断を効率よく行なうことができる。

以上、本発明の実施例について述べてきたが、本発明は、上述の実施例に限定されるも
のではなく変形して実施することが可能である。例えば、上述の実施例では、ＴＯＦ法の
パルスシーケンスを適用した非造影ＭＲＡ撮影によって動脈瘤を含む血管の３次元形態デ
ータを収集し、更に、ＰＣ法のパルスシーケンスを適用した非造影ＭＲＡ撮影によって動
脈瘤内の３次元流速データを収集する場合について述べたが、他の方法によって３次元形
態データや３次元流速データを収集してもよい。例えば、血管の３次元形態データは、Gd
-DTPA等の造影剤を注入した状態で行なわれる造影ＭＲＡ撮影や磁化の定常状態を利用し
たＴrueＳＳＦＰ法あるいは特開２０００−５１４４号公報等に記載されているＦＢＩ法
のパルスシーケンスを適用した非造影ＭＲＡ撮影によって収集してもよい。

又、頭部の脳動脈を診断対象部位とし、この脳動脈に発生した動脈瘤に対して破裂リス
クの判定を行なう場合について述べたが、他の部位に発生した動脈瘤あるいは静脈瘤に対
する破裂リスクの判定であっても構わない。

更に、領域拡張法を用いて動脈瘤の輪郭抽出を行なう場合について述べたが、例えば、
maching-cube法のような他の輪郭抽出法を用いてもよい。又、独立したユニット（即ち、
形態データ生成部１１と流速データ生成部１２）によって３次元形態データと３次元流速
データを生成する場合について述べたが、同一のユニットによってこれらのデータを生成
してもよい。

一方、上述の実施例では、３次元形態データの各ボクセルに対応させて３次元流速デー
タを生成する場合について述べたが、これに限定されるものではなく、３次元流速データ
は任意の位置において生成してもよい。又、破裂リスクの判定に使用される閾値β１乃至
β３の値は、図５に示した数値に限定されない。

更に、上述の実施例における画像データ生成部１４は、３次元形態データをレンダリン
グ処理して３次元画像データを生成する場合について述べたが、ＭＩＰ（maximum intens
ity projection）画像データであっても構わない。この場合、表示部２０は、画像データ
生成部１４が生成した診断対象部位のＭＩＰ画像データに高リスク部位検出部１９が検出
した高リスク部位の位置情報を重畳し、更に、破裂リスク判定部１８による破裂リスクの
判定結果を付加して表示する。

又、上述の実施例では、診断パラメータとしてパラメータ算出部１７が算出あるいは抽
出した動脈瘤のアスペクト比や動脈瘤内壁の平均剪断応力及び最小剪断応力と所定の閾値
β１乃至閾値β３とを比較することによって動脈瘤が有する破裂リスクのレベルを判定す
る場合について述べたが、破裂リスクのレベル判定に用いる診断パラメータの数は上述の
３つに限定されない。

１０…ＭＲ信号収集部
１１…形態データ生成部
１２…流速データ生成部
１３…データ記憶部
１４…画像データ生成部
１５…輪郭抽出部
１６…剪断応力計測部
１７…パラメータ算出部
１７１…アスペクト比算出部
１７２…平均剪断応力算出部
１７３…最小剪断応力抽出部
１７４…最大剪断応力抽出部
１８…破裂リスク判定部
１９…高リスク部位検出部
２０…表示部
２１、２１ａ…入力部
２２…心拍時相計測部
２３、２３ａ…システム制御部
２４…データ保管部
１００…ＭＲＩ装置
２００…データ解析装置

上記課題を解決するために、請求項１に係る本発明のＭＲＩ装置は、所定のパルスシー
ケンスを適用したＭＲＡ撮影により被検体の動脈瘤を含む診断対象部位が発生するＭＲ信
号を収集するＭＲ信号収集手段と、前記ＭＲ信号を処理して前記診断対象部位の３次元形
態データを生成する形態データ生成手段と、前記３次元形態データに基づいて前記動脈瘤
の輪郭データを生成する輪郭抽出手段と、前記ＭＲ信号を処理して前記診断対象部位の３
次元流速データを生成する流速データ生成手段と、前記輪郭データと前記３次元流速デー
タに基づいて前記動脈瘤の頂点の近傍における剪断応力を計測する剪断応力計測手段と、
前記動脈瘤の頂点の近傍における剪断応力の計測結果に基づいて前記動脈瘤の破裂リスク
を判定する破裂リスク判定手段と、前記破裂リスクの判定結果を表示する表示手段とを備
えたことを特徴としている。

又、請求項８に係る本発明のデータ処理装置は、ＭＲＩ装置によって収集された動脈瘤
を含む診断対象部位の３次元形態データ及び３次元流速データを保管するデータ保管手段
と、前記３次元形態データに基づいて前記動脈瘤の輪郭データを生成する輪郭抽出手段と
、前記輪郭データと前記３次元流速データに基づいて前記動脈瘤の頂点の近傍における剪
断応力を計測する剪断応力計測手段と、前記動脈瘤の頂点の近傍における剪断応力の計測
結果に基づいて前記動脈瘤の破裂リスクを判定する破裂リスク判定手段と、前記破裂リス
クの判定結果を表示する表示手段とを備えたことを特徴としている。

Claims (2)

1. 所定のパルスシーケンスを適用したＭＲＡ撮影により被検体の動脈瘤を含む診断対象部
   位が発生するＭＲ信号を収集するＭＲ信号収集手段と、
   前記ＭＲ信号を処理して前記診断対象部位の３次元形態データを生成する形態データ生成
   手段と、
   前記３次元形態データに基づいて前記動脈瘤の輪郭を抽出し、輪郭データを生成する輪郭
   抽出手段と、
   前記ＭＲ信号を処理して前記診断対象部位の３次元流速データを生成する流速データ生成
   手段と、
   前記輪郭データと前記３次元流速データに基づいて前記動脈瘤における剪断応力を計測す
   る剪断応力計測手段と、
   前記剪断応力の計測結果あるいはこの剪断応力の計測結果と前記輪郭データに基づいて前
   記動脈瘤の破裂リスクを判定する破裂リスク判定手段と、
   前記破裂リスクの判定結果を表示する表示手段とを
   備えたことを特徴とするＭＲＩ装置。
2. ＭＲＩ装置によって収集された動脈瘤を含む診断対象部位の３次元形態データ及び３次
   元流速データを保管するデータ保管手段と、
   前記３次元形態データに基づいて前記動脈瘤の輪郭を抽出し、輪郭データを生成する輪郭
   抽出手段と、
   前記輪郭データと前記３次元流速データに基づいて前記動脈瘤における剪断応力を計測す
   る剪断応力計測手段と、
   前記剪断応力の計測結果あるいはこの剪断応力の計測結果と前記輪郭データに基づいて前
   記動脈瘤の破裂リスクを判定する破裂リスク判定手段と、
   前記破裂リスクの判定結果を表示する表示手段とを
   備えたことを特徴とするデータ処理装置。

Images