(57)【要約】【課題】頭部画像から精度良く病変部を検出する。また、特定の血管部位に注目して観察することを可能とする。【解決手段】医用画像処理装置では、頭部を撮影したＭＲＡ画像から脳血管領域が抽出された画像を用いてＧＣフィルタバンクにより１次検出が行われる。ＧＣフィルタバンクでは分析バンクにおいて多重解像度解析が行われ、解像度レベル毎に部分画像からベクトル集中度が算出される。次に、ベクトル集中度が閾値を超えるボクセルを１、閾値以下のボクセルを０とする重み画像が作成される。次いで、再構成バンクにおいて各部分画像に重み画像が乗算された部分画像を用いて元の画像の再構成が行われる。この再構成画像に再現された候補領域について２次検出、３次検出が行われ、その検出結果が表示部に表示出力される。【選択図】図２
2008104886.tif 000002
【特許請求の範囲】
【請求項１】
３次元画像である頭部ＭＲＡ画像をフィルタ処理して多重解像度解析し、解像度レベル毎に分解した部分画像を用いて、解像度レベル毎にベクトル集中度を算出する分析手段と、
前記部分画像から元の頭部ＭＲＡ画像を再構成するにあたり、前記算出されたベクトル集中度が所定値となる領域を重み付けした前記部分画像の重み画像を作成し、当該重み画像を乗算した前記部分画像から元の画像を再現する再構成手段と、
前記再構成された頭部ＭＲＡ画像を用いて、当該頭部ＭＲＡ画像に再現された病変部の候補領域のうち、正常な血管である偽陽性候補領域を削除する削除手段と、
を備えることを特徴とする医用画像処理装置。
【請求項２】
前記分析手段は、前記頭部ＭＲＡ画像に対し、下記式１１、１２により表されるフィルタＨ（ｚ）、Ｇ（ｚ）によるフィルタ処理を施すことにより多重解像度解析を行うことを特徴とする請求項１に記載の医用画像処理装置。
【数１】
2008104886.tif 000013
なお、上記式におけるｚは、ｚ変換を示す。
【請求項３】
前記分析手段は、前記頭部ＭＲＡ画像の３次元を構成するｘ、ｙ、ｚ方向についてそれぞれフィルタ処理を施し、ｘ、ｙ、ｚ方向のそれぞれに係る部分画像を得ることを特徴とする請求項２に記載の医用画像処理装置。
【請求項４】
前記再構成手段は、前記重み画像を乗算した前記部分画像に対し、下記式１１、１３、１４で表されるフィルタＬ（ｚ）、Ｋ（ｚ）、Ｈ（ｚ）Ｌ（ｚ）によるフィルタ処理を施すことにより元の画像を再現することを特徴とする請求項２又は３に記載の医用画像処理装置。
【数２】
2008104886.tif 000014
なお、上記式におけるｚは、ｚ変換を示す。
【請求項５】
前記削除手段は、前記再構成された頭部ＭＲＡ画像に再現された病変部の候補領域について特徴量を算出し、当該特徴量に基づいて前記候補領域から偽陽性候補領域を削除することを特徴とする請求項１〜４の何れか一項に記載の医用画像処理装置。
【請求項６】
前記特徴量は、前記再構成された頭部ＭＲＡ画像に再現された病変部の候補領域の大きさ、球形度、当該候補領域におけるベクトル集中度の平均若しくはベクトル集中度の最大値のうち、少なくとも１以上であることを特徴とする請求項５に記載の医用画像処理装置。
【請求項７】
前記病変部は、血管において生じる瘤であり、
前記頭部ＭＲＡ画像から血管像を抽出する抽出手段を備え、
前記再構成手段は、前記血管像が抽出された画像を用いて多重解像度解析及び元の画像の再構成を行うことを特徴とする請求項１〜６の何れか一項に記載の医用画像処理装置。
【請求項８】
前記頭部ＭＲＡ画像の血管像に含まれる一又は複数の血管部位を判別する画像制御手段と、
前記頭部ＭＲＡ画像をＭＩＰ表示する表示手段と、
表示操作を行うための操作手段と、
前記頭部ＭＲＡ画像において表示対象の血管部位が前記操作手段を介して選択操作されると、前記選択された血管部位に対応する血管像のみを前記頭部ＭＲＡ画像から抽出し、当該抽出された血管像のみを前記表示手段上に表示させる表示制御手段と、
を備えることを特徴とする請求項１〜７の何れか一項に記載の医用画像処理装置。
【請求項９】
３次元画像である頭部ＭＲＡ画像をフィルタ処理して多重解像度解析し、解像度レベル毎に分解した部分画像を用いて、解像度レベル毎にベクトル集中度を算出する分析工程と、
前記部分画像から元の頭部ＭＲＡ画像を再構成するにあたり、前記算出されたベクトル集中度が所定値となる領域を重み付けした前記部分画像の重み画像を作成し、当該重み画像を乗算した前記部分画像から元の画像を再現する再構成工程と、
前記再構成された頭部ＭＲＡ画像を用いて、当該頭部ＭＲＡ画像に再現された病変部の候補領域のうち、正常な血管である偽陽性候補領域を削除する削除工程と、
を含むことを特徴とする画像処理方法。
【請求項１０】
前記分析工程は、前記頭部ＭＲＡ画像に対し、下記１１、１２により表されるフィルタＨ（ｚ）、Ｇ（ｚ）によるフィルタ処理を施すことにより多重解像度解析を行うことを特徴とする請求項９に記載の画像処理方法。
【数３】
2008104886.tif 000015
なお、上記式におけるｚは、ｚ変換を示す。
【請求項１１】
前記分析工程は、前記頭部ＭＲＡ画像の３次元を構成するｘ、ｙ、ｚ方向についてそれぞれフィルタ処理を施し、ｘ、ｙ、ｚ方向のそれぞれに係る部分画像を得ることを特徴とする請求項１０に記載の画像処理方法。
【請求項１２】
前記再構成工程は、前記重み画像を乗算した前記部分画像に対し、下記式１１、１３、１４で表されるフィルタＬ（ｚ）、Ｋ（ｚ）、Ｈ（ｚ）Ｌ（ｚ）によるフィルタ処理を施すことにより元の画像を再現することを特徴とする請求項１０又は１１に記載の画像処理方法。
【数４】
2008104886.tif 000016
なお、上記式におけるｚは、ｚ変換を示す。
【請求項１３】
前記削除工程は、前記再構成された頭部ＭＲＡ画像に再現された病変部の候補領域について特徴量を算出し、当該特徴量に基づいて前記候補領域から偽陽性候補領域を削除することを特徴とする請求項９〜１２の何れか一項に記載の画像処理方法。
【請求項１４】
前記特徴量は、前記再構成された頭部ＭＲＡ画像に再現された病変部の候補領域の大きさ、球形度、当該候補領域におけるベクトル集中度の平均若しくはベクトル集中度の最大値のうち、少なくとも１以上であることを特徴とする請求項１３に記載の画像処理方法。
【請求項１５】
前記病変部は、血管において生じる瘤であり、
前記頭部ＭＲＡ画像から血管像を抽出する抽出工程を含み、
前記再構成工程は、前記血管像が抽出された画像を用いて多重解像度解析及び元の画像の再構成を行うことを特徴とする請求項９〜１４の何れか一項に記載の画像処理方法。
【請求項１６】
前記頭部ＭＲＡ画像の血管像に含まれる一又は複数の血管部位を判別する画像制御工程と、
前記頭部ＭＲＡ画像を表示手段にＭＩＰ表示する表示工程と、
前記頭部ＭＲＡ画像において表示対象の血管部位が操作手段を介して選択操作されると、前記選択された血管部位に対応する血管像のみを前記頭部ＭＲＡ画像から抽出し、当該抽出された血管像のみを前記表示手段上に表示させる表示制御工程と、
を含むことを特徴とする請求項９〜１５の何れか一項に記載の画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、患者の頭部を撮影して得られた頭部画像の画像解析、画像処理を行う医用画像処理装置及び画像処理方法に関する。
【背景技術】
【０００２】
近年、磁気共鳴撮影装置（以下、ＭＲＩ；Magnetic Resonance Imagingという）の普及と高性能化に伴い、脳ドック検査の件数が急速に増加してきている。脳ドック検査の目的の一つは、血管において生じる未破裂脳動脈瘤を早期に発見し、適切な処置や治療を行うことによって、脳動脈瘤の破裂によるくも膜下出血等、重篤な疾患の発病を防止することにある。
【０００３】
医師による未破裂脳動脈瘤の検出は、ＭＲＩで血管内の血液の流れを画像化したＭＲＡ画像（Magnetic Resonance Angiography）を用いて行われる。通常、医師による読影時には３次元画像データを様々な角度からＭＩＰ処理（Maximum Intensity Projection；最大値輝度投影法）によって２次元化した画像が用いられるが、血管に生じる未破裂脳動脈瘤は小さいため、重なって表示される周囲の血管像との識別が必要となり、医師の疲労が激しい。また、この疲労による見落としの可能性も考えられる。
【０００４】
従来から、このような医師の診断支援を行うため、画像処理によって病変部分の画像領域を検出する装置が開発されている（例えば、特許文献１、非特許文献１、２参照）。この装置は、一般にＣＡＤ（Computer-Aided Diagnosis）と呼ばれている。
【０００５】
一方、読影時にＭＩＰによる画像を医用画像処理装置において生成し、表示するとき、視線方向によっては複数の血管が重複して表示される場合がある。この場合、小さな未破裂動脈瘤を検出するには、医師が注目したい血管像と周囲の血管像との識別が必要なために、多くの読影時間を必要とし、医師の疲労も激しい。疲労によって本来検出すべき未破裂動脈瘤を見落とす可能性がある。
【０００６】
動脈瘤が多発する血管部位は、中大脳動脈分岐部、前交通動脈、内頸動脈後交通動脈分岐部等であることが知られている。従って、これらの多発部位を詳しく観察するには、医師が関心のある血管部位の血管像のみに注目したターゲットＭＩＰ画像を生成すればよい。
【０００７】
また、複数の血管部位が交差する領域においては、視線方向に沿った深さ情報を用いて当該深さ情報が大きい方の血管像、つまり視線方向において後方側にある血管部位の交差部分の輝度を低下させて表示する技術も開示されている（例えば、特許文献２参照）。この方法によれば、血管像に遠近感をもたせることができ、手前側の血管部位を観察しやすくなる。
【特許文献１】特開２００２−１１２９８６号公報
【特許文献２】特開平５−２７７０９１号公報
【非特許文献１】林則夫他、モルフォルジー処理を利用した頭部ＭＲＩ画像における小脳および脳患部の自動抽出法、医用画像情報学会雑誌、vol21.no1.pp109-115,2004
【非特許文献２】横山龍二郎他、脳ＭＲ画像におけるラクナ梗塞領域の自動検出、日本放射線技術学会雑誌、58(3),399-405,2002
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかしながら、未破裂脳動脈瘤を血管組織と識別して精度良く検出するような手法はいまだ提案されていない。
【０００９】
また、従来の医用画像処理装置では、血管像中の各血管部位を自動判別することができないため、上記のターゲットＭＩＰ画像を自動で作成するためには、３次元の画像データから関心のある血管部位の血管像のみを手動で指定する操作が必要となる。このような操作は容易ではなく、作成に時間を要する。
【００１０】
また、上記特許文献２に記載の方法は観察者に近い側の血管部位を観察しやすくするものである。従って、医師が観察したい視線方向に応じて作成されたＭＩＰ画像において、医師の関心がある血管部位の手前側に別の血管部位が交差している場合はそれを除去することができず、所望の観察方向から所望の血管部位を詳細に観察できるとは限らない。
【００１１】
本発明の課題は、頭部画像から精度良く病変部を検出することである。また、特定の血管部位に注目して観察することを可能とすることである。
【課題を解決するための手段】
【００１２】
請求項１に記載の発明によれば、
３次元画像である頭部ＭＲＡ画像をフィルタ処理して多重解像度解析し、解像度レベル毎に分解した部分画像を用いて、解像度レベル毎にベクトル集中度を算出する分析手段と、
前記部分画像から元の頭部ＭＲＡ画像を再構成するにあたり、前記算出されたベクトル集中度が所定値となる領域を重み付けした前記部分画像の重み画像を作成し、当該重み画像を乗算した前記部分画像から元の画像を再現する再構成手段と、
前記再構成された頭部ＭＲＡ画像を用いて、当該頭部ＭＲＡ画像に再現された病変部の候補領域のうち、正常な血管である偽陽性候補領域を削除する削除手段と、
を備えることを特徴とする医用画像処理装置が提供される。
【００１３】
請求項２に記載の発明によれば、
前記分析手段は、前記頭部ＭＲＡ画像に対し、下記式１１、１２により表されるフィルタＨ（ｚ）、Ｇ（ｚ）によるフィルタ処理を施すことにより多重解像度解析を行うことを特徴とする請求項１に記載の医用画像処理装置が提供される。
【数１】
2008104886.tif 000003
なお、上記式におけるｚは、ｚ変換を示す。
【００１４】
請求項３に記載の発明によれば、
前記分析手段は、前記頭部ＭＲＡ画像の３次元を構成するｘ、ｙ、ｚ方向についてそれぞれフィルタ処理を施し、ｘ、ｙ、ｚ方向のそれぞれに係る部分画像を得ることを特徴とする請求項２に記載の医用画像処理装置が提供される。
【００１５】
請求項４に記載の発明によれば、
前記再構成手段は、前記重み画像を乗算した前記部分画像に対し、下記式１１、１３、１４で表されるフィルタＬ（ｚ）、Ｋ（ｚ）、Ｈ（ｚ）Ｌ（ｚ）によるフィルタ処理を施すことにより元の画像を再現することを特徴とする請求項２又は３に記載の医用画像処理装置が提供される。
【数２】
2008104886.tif 000004
なお、上記式におけるｚは、ｚ変換を示す。
【００１６】
請求項５に記載の発明によれば、
前記削除手段は、前記再構成された頭部ＭＲＡ画像に再現された病変部の候補領域について特徴量を算出し、当該特徴量に基づいて前記候補領域から偽陽性候補領域を削除することを特徴とする請求項１〜４の何れか一項に記載の医用画像処理装置が提供される。
【００１７】
請求項６に記載の発明によれば、
前記特徴量は、前記再構成された頭部ＭＲＡ画像に再現された病変部の候補領域の大きさ、球形度、当該候補領域におけるベクトル集中度の平均若しくはベクトル集中度の最大値のうち、少なくとも１以上であることを特徴とする請求項５に記載の医用画像処理装置が提供される。
【００１８】
請求項７に記載の発明によれば、
前記病変部は、血管において生じる瘤であり、
前記頭部ＭＲＡ画像から血管像を抽出する抽出手段を備え、
前記再構成手段は、前記血管像が抽出された画像を用いて多重解像度解析及び元の画像の再構成を行うことを特徴とする請求項１〜６の何れか一項に記載の医用画像処理装置が提供される。
【００１９】
請求項８に記載の発明によれば、
前記頭部ＭＲＡ画像の血管像に含まれる一又は複数の血管部位を判別する画像制御手段と、
前記頭部ＭＲＡ画像をＭＩＰ表示する表示手段と、
表示操作を行うための操作手段と、
前記頭部ＭＲＡ画像において表示対象の血管部位が前記操作手段を介して選択操作されると、前記選択された血管部位に対応する血管像のみを前記頭部ＭＲＡ画像から抽出し、当該抽出された血管像のみを前記表示手段上に表示させる表示制御手段と、
を備えることを特徴とする請求項１〜７の何れか一項に記載の医用画像処理装置が提供される。
【００２０】
請求項９に記載の発明によれば、
３次元画像である頭部ＭＲＡ画像をフィルタ処理して多重解像度解析し、解像度レベル毎に分解した部分画像を用いて、解像度レベル毎にベクトル集中度を算出する分析工程と、
前記部分画像から元の頭部ＭＲＡ画像を再構成するにあたり、前記算出されたベクトル集中度が所定値となる領域を重み付けした前記部分画像の重み画像を作成し、当該重み画像を乗算した前記部分画像から元の画像を再現する再構成工程と、
前記再構成された頭部ＭＲＡ画像を用いて、当該頭部ＭＲＡ画像に再現された病変部の候補領域のうち、正常な血管である偽陽性候補領域を削除する削除工程と、
を含むことを特徴とする画像処理方法が提供される。
【００２１】
請求項１０に記載の発明によれば、
前記分析工程は、前記頭部ＭＲＡ画像に対し、下記式１１、１２により表されるフィルタＨ（ｚ）、Ｇ（ｚ）によるフィルタ処理を施すことにより多重解像度解析を行うことを特徴とする請求項９に記載の画像処理方法が提供される。
【数３】
2008104886.tif 000005
なお、上記式におけるｚは、ｚ変換を示す。
【００２２】
請求項１１に記載の発明によれば、
前記分析工程は、前記頭部ＭＲＡ画像の３次元を構成するｘ、ｙ、ｚ方向についてそれぞれフィルタ処理を施し、ｘ、ｙ、ｚ方向のそれぞれに係る部分画像を得ることを特徴とする請求項１０に記載の画像処理方法が提供される。
【００２３】
請求項１２に記載の発明によれば、
前記再構成工程は、前記重み画像を乗算した前記部分画像に対し、下記式１１、１３、１４で表されるフィルタＬ（ｚ）、Ｋ（ｚ）、Ｈ（ｚ）Ｌ（ｚ）によるフィルタ処理を施すことにより元の画像を再現することを特徴とする請求項１０又は１１に記載の画像処理方法が提供される。
【数４】
2008104886.tif 000006
なお、上記式におけるｚは、ｚ変換を示す。
【００２４】
請求項１３に記載の発明によれば、
前記削除工程は、前記再構成された頭部ＭＲＡ画像に再現された病変部の候補領域について特徴量を算出し、当該特徴量に基づいて前記候補領域から偽陽性候補領域を削除することを特徴とする請求項９〜１２の何れか一項に記載の画像処理方法が提供される。
【００２５】
請求項１４に記載の発明によれば、
前記特徴量は、前記再構成された頭部ＭＲＡ画像に再現された病変部の候補領域の大きさ、球形度、当該候補領域におけるベクトル集中度の平均若しくはベクトル集中度の最大値のうち、少なくとも１以上であることを特徴とする請求項１３に記載の画像処理方法が提供される。
【００２６】
請求項１５に記載の発明によれば、
前記病変部は、血管において生じる瘤であり、
前記頭部ＭＲＡ画像から血管像を抽出する抽出工程を含み、
前記再構成工程は、前記血管像が抽出された画像を用いて多重解像度解析及び元の画像の再構成を行うことを特徴とする請求項９〜１４の何れか一項に記載の画像処理方法が提供される。
【００２７】
請求項１６に記載の発明によれば、
前記頭部ＭＲＡ画像の血管像に含まれる一又は複数の血管部位を判別する画像制御工程と、
前記頭部ＭＲＡ画像を表示手段にＭＩＰ表示する表示工程と、
前記頭部ＭＲＡ画像において表示対象の血管部位が操作手段を介して選択操作されると、前記選択された血管部位に対応する血管像のみを前記頭部ＭＲＡ画像から抽出し、当該抽出された血管像のみを前記表示手段上に表示させる表示制御工程と、
を含むことを特徴とする請求項９〜１５の何れか一項に記載の画像処理方法が提供される。
【発明の効果】
【００２８】
請求項１〜６、９〜１４に記載の発明によれば、頭部画像に対して多重解像度解析を行って得られた部分画像によって解像度レベル毎にベクトル集中度を算出する。解像度レベル毎にベクトル集中度を算出することにより、様々な大きさの動脈瘤の検出に対応することができ、より精度の高い検出処理が可能となる。また、瘤等の病変部はある一定の高い値のベクトル集中度を呈するが、本発明によればこのような病変部である可能性が高い領域のみ再現されるよう再構成を行うので、病変部の候補領域のみ画像化した再構成画像を得ることができる。再構成画像を用いて特徴量の算出を行うことにより、算出する特徴量から候補領域以外の画像要素を排除することができ、候補領域についての特徴量を正確に算出することが可能となる。従って、検出処理自体の精度をさらに向上させることが可能となる。
【００２９】
請求項７、１５に記載の発明によれば、実質的にベクトル集中度を算出する領域を、瘤の病変部が存在する血管領域に絞ることができ、演算時間の短縮化を図ることが可能となる。
【００３０】
請求項８、１６に記載の発明によれば、医師が観察したい特定の血管部位の血管像のみを対象画像から抽出して観察することが可能となる。血管像は複数の血管部位が重なって表示されることがあるため、その重複部分の観察がしづらくなる。よって、医師により選択された特定の血管部位のみを表示することにより、血管部位の重複を解消することができ、医師が読影しやすい環境を提供することができる。
【発明を実施するための最良の形態】
【００３１】
〈第１実施形態〉
まず、構成を説明する。
図１に、本実施形態における医用画像処理装置１０の構成を示す。
この医用画像処理装置１０は、検査撮影により得られた医用画像を画像解析することにより、当該医用画像から病変部の候補領域を検出するものである。
なお、この医用画像処理装置１０を、医用画像を生成する画像生成装置や医用画像を保存・管理するサーバ、医師の読影に付すため、サーバに保存された医用画像を取り寄せて表示手段上に表示する読影端末等、各種装置がネットワークを介して接続された医用画像システムに設けることとしてもよい。また、本実施形態では医用画像処理装置１０単体で本発明を実現する例を説明するが、医用画像処理装置１０における機能を上記医用画像システムの各構成装置に分散させて医用画像システム全体で本発明を実現することとしてもよい。
【００３２】
以下、医用画像処理装置１０の各部について説明する。
医用画像処理装置１０は、図１に示すように、制御部１１、操作部１２、表示部１３、通信部１４、記憶部１５、病変候補検出部１６を備えて構成されている。
【００３３】
制御部１１は、ＣＰＵ（Central Processing Unit）、ＲＡＭ（Random Access Memory
）等を備えて構成されており、記憶部１５に格納されている各種制御プログラムを読み出して各種演算を行うとともに、各部１２〜１６における処理動作を統括的に制御する。
【００３４】
操作部１２は、キーボードやマウス等を備え、オペレータによりこれらが操作されると、その操作に応じた操作信号を生成して制御部１１に出力する。なお、表示部１３におけるディスプレイと一体に構成したタッチパネルを備えることとしてもよい。
【００３５】
表示部１３は、ＬＣＤ（Liquid Crystal Display）等の表示手段を備え、制御部１１からの指示に応じてこの表示手段上に、各種操作画面や、医用画像、医用画像から検出された病変候補の検出結果、その検出情報等の各種表示情報を表示させる。
【００３６】
通信部１４は、通信用のインターフェイスを備え、ネットワーク上の外部装置と情報の送受信を行う。例えば、通信部１４は画像生成装置から生成された医用画像を受信する、医用画像処理装置１０における病変候補の検出情報を読影端末に送信する等の通信動作を行う。
【００３７】
記憶部１５は、制御部１１において用いられる制御プログラム、病変候補検出部１６において用いられる検出処理等の各種処理プログラムの他、各プログラムの実行に必要なパラメータやその処理結果等のデータが記憶されている。
また、記憶部１５は、病変候補の検出対象である医用画像やその検出結果の情報等を記憶している。
【００３８】
病変候補検出部１６は、記憶部１５に記憶されている処理プログラムとの協働により、処理対象画像に必要に応じて各種画像処理（階調変換処理、鮮鋭性調整処理、ダイナミックレンジ圧縮処理等）を施す。また、病変候補検出部１６は検出処理を実行し、その検出結果を出力する。検出処理の内容については、後述する。
【００３９】
次に、上記医用画像処理装置１０による病変候補の検出処理について説明する。
本実施形態では、ＭＲＩにより患者の頭部を撮影し、脳内血管の血流を画像化したＭＲＡ画像（３次元画像）から、未破裂の脳動脈瘤の病変候補を検出する例を説明する。脳動脈瘤は、動脈の壁内にできる膨隆（拡張）のことであり、動脈壁に血流による圧力がかかることにより生じるものである。脳動脈瘤の内部では血栓が生じやすく、この脳動脈瘤が破裂すると、くも膜下出血等の重篤な疾患を発症することとなる。
【００４０】
図２は、検出処理の流れを説明するフローチャートである。この検出処理は、前述したように、病変候補検出部１６が記憶部１５に記憶される検出処理プログラムを読み込むことにより実行される処理である。
図２に示すように、検出処理ではまずＭＲＡの３次元画像データの入力が行われる（ステップＳ１）。具体的には、病変候補検出部１６により記憶部１５に記憶されていた処理対象のＭＲＡ画像の読み込みが行われる。
【００４１】
ここで、ＭＲＩにより得られる画像について説明する。
ＭＲＩは、磁場内における核磁気共鳴（以下、ＮＭＲ；Nuclear Magnetic Resonance）を利用して画像を得る方法である。
ＮＭＲでは、被検体を静磁場中に置き、その後、被検体において検出対象とする原子核の共鳴周波数のＲＦパルス（電波）を照射する。医用上は、通常、人体に多く存在する水を構成する水素原子の共鳴周波数が用いられる。被検体にＲＦパルスが照射されると、励起現象が生じ、共鳴周波数に共鳴する原子の原子核スピンの位相がそろうとともに、原子核スピンがＲＦパルスのエネルギーを吸収する。この励起状態でＲＦパルスの照射を止めると、緩和現象を生じ、原子核スピンの位相が不均一化するとともに、原子核スピンがエネルギーを放出する。この位相の緩和の時定数がＴ２、エネルギーの緩和の時定数がＴ１である。
【００４２】
ＭＲＩでは、撮像方法を変えることにより、検査目的に応じた様々な種類の画像を得ることが可能である。例えば、繰り返し時間ＴＲ、エコー時間ＴＥについて、ＴＲ＝Ｔ１、ＴＥ＜＜Ｔ２と撮影条件を調整したものをＴ１強調画像、ＴＲ＞＞Ｔ１、ＴＥ＝Ｔ２と撮影条件を調整したものをＴ２強調画像というが、両者は画像のコントラストに差異があり、主にＴ１強調画像は解剖構造の検出に、Ｔ２強調画像は病変部の検出に用いられる。また、ＦＬＡＩＲ法により撮影された画像は水からの信号を減衰させたＴ２強調画像であるが、特にＦＬＡＩＲ画像と呼ばれる。
【００４３】
ＭＲＡは、ＭＲＩにおける血管撮像方法をいう。ＭＲＩでは、被検体の足から頭の方向（この方向を体軸という）に傾斜磁場をかけることにより特定のスライス（断層）にのみエネルギーを吸収させることが可能である。撮像時には、スライス内の血管ではＲＦパルスで血液が飽和された状態となっているが、血管は常に血液の流れがあるため経時により非飽和の血流が流入するとそのスライスにおける信号強度が増加する。ＭＲＡは、この高信号を画像化することにより血流のある血管を画像化する方法である。
【００４４】
図３Ａに、ＭＲＡ画像例を示す。
図３Ａに示すように、血流のある血管領域は高信号となるので、ＭＲＡ画像では血管領域は白く表れている。
【００４５】
このようなＭＲＡ画像の３次元データが入力されると、当該３次元画像データに対し、候補を検出する準備段階として前処理が施される（ステップＳ２）。前処理としては、画像データの正規化処理、階調変換処理が行われる。正規化は、ボクセルを構成する全ての辺が等サイズの３次元画像データとなるように、線形補間法により変換することにより行われる。そして、この等サイズのボクセルに変換された３次元画像データに対し、濃度階調変換処理が施され、各ボクセルの持つ信号値が０〜１０２４の濃度階調に線形変換される。このとき、高信号値であるほど濃度値１０２４に近く、低信号値であるほど濃度値０に近くなるように変換される。なお、濃度階調のレンジは、０〜１０２４に限らず適宜設定可能である。
【００４６】
前処理が終了すると、３次元ＭＲＡ画像から血管の画像領域が抽出される（ステップＳ３）。まず、閾値処理が施され、ＭＲＡ画像の２値化が行われる。一般的に、ＭＲＡ画像では図３Ａに示すように血管領域は白く、その他の領域は黒っぽく表れるため、２値化画像では血管領域はその他の領域と異なる値となる。よって、領域拡張法により血管領域の抽出を行う。まず、２値化画像を用いて始点となるボクセル（最も白く高濃度値のボクセル）を決定し、２値化処理前の３次元ＭＲＡ画像においてその始点と決定されたボクセルの近傍２６ボクセルを調べ、ある判定条件（例えば、濃度値が５００以上であること）を満たす近傍ボクセルを血管領域と判断する。そして、この血管領域と判断された近傍ボクセルについても上記と同様の処理を繰り返す。このように、領域を拡張しながら判定条件を満たすボクセルを順次抽出することにより、血管領域を抽出することができる。図３Ｂに、図３ＡのＭＲＡ画像から抽出した血管領域を白（濃度値１０２４）、その他の領域を黒（濃度値０）で２値化した図を示す。
【００４７】
次いで、抽出された血管領域の３次元ＭＲＡ画像に対し、図４に示すようなベクトル集中度フィルタを用いたフィルタ処理が行われ、フィルタ処理により出力された処理画像から脳動脈瘤の１次候補領域の検出が行われる（ステップＳ４）。ベクトル集中度フィルタは、各ボクセルにおけるベクトル集中度を算出し、当該算出されたベクトル集中度の値をそのボクセルのボクセル値として画像化し出力するものである。ベクトル集中度は、濃度変化の勾配ベクトルの向きに着目し、ある注目点に対して近傍領域の勾配ベクトルがどの程度集中しているかを評価するものである。
【００４８】
図５に、脳動脈瘤モデルと血管モデルを示す。
図５に示すように、脳動脈瘤の場合には線形状の血管上に球状の瘤が存在するため、勾配ベクトル（図中の矢印は、勾配ベクトルの向きを示す）は瘤の中心へ向かう傾向がある。一方、血管は線形状であるためそのような傾向は生じない。そのため、脳動脈瘤モデルの形状に近い領域は他の血管領域に比べてベクトル集中度の値が高い。従って、ベクトル集中度フィルタによりベクトル集中度が高い領域のみを出力することにより、脳動脈瘤の１次候補領域を検出することができる。
【００４９】
具体的には、血管領域の抽出画像上に図４に示すような注目ボクセルＰを走査させ、注目ボクセルＰを中心とする半径Ｒの球の範囲内に、抽出された血管領域が存在した場合にベクトル集中度を算出する。
ベクトル集中度は、下記式１により算出される。
【数５】
2008104886.tif 000007
ここで、角度θは注目ボクセルＰから周辺ボクセルＱｊまでの方向ベクトルと、周辺ボクセルＱｊにおけるベクトルの方向との間の角度を示し、Ｍは演算の対象となった周辺ボクセルＱｊの個数を示している（図４参照）。
【００５０】
血管領域を構成する各ボクセルについてベクトル集中度が求められると、当該ベクトル集中度をボクセル値とする、図６に示すようなフィルタ処理画像が出力される。ベクトル集中度は０−１の範囲で出力されるので、図６では、ベクトル集中度が大きいほど（１に近いほど）、白く高濃度に表れるように画像化している。
【００５１】
次に、フィルタ処理画像について例えば閾値０．５を用いて２値化する閾値処理が施され、１次候補領域の検出が行われる。すなわち、２値化により閾値０．５より大きいベクトル集中度を有するボクセルからなる領域が１次候補領域として抽出される。図７Ａは図６に示すフィルタ処理画像の一部分を示す図である。このフィルタ処理画像を２値化することにより、図７Ｂに示すような２値化画像が得られる。図７Ｂに示すように、２値化画像において白く表れている、つまりベクトル集中度が大きい画像領域が１次候補領域として出力される。
【００５２】
なお、２値化するための閾値は、予め脳動脈瘤の存在が判明している教師画像を用いて決定しておく。すなわち、教師画像を２値化することにより、その存在が既に判明している脳動脈瘤の画像領域のみを抽出することができるような閾値が求められる。
また、閾値はｐ−タイル法等により統計的に解析して得ることとしてもよい。ｐ−タイル法は、濃度ヒストグラムを求め、この濃度ヒストグラムにおいてある一定の面積比率ｐ％を占めるところの濃度値を閾値として求める方法である。本実施形態では、フィルタ処理画像の濃度ヒストグラムを求め、最高濃度値側から面積比率ｐ％を占めるところの濃度値を閾値として決定する。
【００５３】
次いで、２次検出を行うため、１次候補領域の特徴を示す特徴量が算出される（ステップＳ５）。脳動脈瘤はある一定の大きさを持ち、形状が球状であることから、本実施形態では特徴量として候補領域の大きさ、球形度、領域内の各ボクセルにおけるベクトル集中度の平均値を算出することとする。しかし、脳動脈瘤を特徴付けることができるのであればどのような特徴量を用いるかは特に限定せず、例えば各ボクセルのベクトル集中度の最大値を算出することとしてもよいし、各ボクセルの濃度値の標準偏差等を算出することとしてもよい。
【００５４】
大きさの特徴量としては、候補領域を構成する各ボクセルの体積が算出される。ここでは、処理時間の短縮のため実際の体積ではなく、ボクセルの個数を算出してこれを体積を示す指標値として以後の演算に用いることとする。また、球形度の特徴量は、図８に示すように１次候補領域の体積と同一体積の球を、１次候補領域の重心と球の重心が一致するように配置したときに、この球と一致する１次候補領域部分の体積と１次候補領域の全体積との比率から求める。
【００５５】
特徴量が算出されると、当該特徴量に基づいて２次検出が行われる（ステップＳ６）。２次検出では、脳動脈瘤と正常組織である血管との識別が行われる。識別はルールベース法を用いた識別器の例を示すが、これに限らず、例えば人工ニューラルネットワーク、サポートベクトルマシン、判別分析等、識別ができるのであればその手法は問わない。
【００５６】
ルールベース法の識別器では、図９Ａ、図９Ｂに示すように、まず大きさに対する球形度、大きさに対するベクトル集中度の平均の関係を示すプロファイルが作成される。識別器では予め２次候補として検出する範囲（図９Ａ、図９Ｂで実線により囲まれた範囲）が決定されており、識別対象の１次候補の各特徴量を変量データとして上記プロファイルを作成した際に、その１次候補の各特徴量の変量データが検出範囲内に存在すれば真陽性候補、存在しなければ偽陽性候補として識別される。すなわち、検出範囲内に特徴量の変量データが分布する１次候補のみが２次検出される。
【００５７】
なお、上記検出範囲は、予め脳動脈瘤、或いは正常な血管であると判明している教師画像を用いて決定される。まず、脳動脈瘤、正常な血管の教師画像からそれぞれ大きさ、球形度、ベクトル集中度の平均が求められ、これら特徴量から大きさに対する球形度、大きさに対するベクトル集中度の平均値の関係のプロファイルが作成される。図９Ａ、図９Ｂでは、○のマーカで示す変量データが真陽性、つまり脳動脈瘤と判明している教師データであり、●のマーカで示す変量データが偽陽性、つまり血管と判明している教師データである。このとき、脳動脈瘤と正常な血管では画像の特徴が異なるため、大きさ、球形度、ベクトル集中度の平均値の特徴量の分布に偏りが生じる。よって、プロファイルにおいて、全ての脳動脈瘤の教師データが含まれる範囲が最小限となるように、大きさ、球形度、ベクトル集中度の平均値について２つの閾値が決定される。この閾値によって囲まれる領域範囲が検出範囲である。
【００５８】
図９Ａ、図９Ｂの例では、閾値を示す４本の点線で囲まれた範囲（その範囲を囲む部分のみ実線で示す）が検出範囲であり、「大きさ−球形度」、「大きさ−ベクトル集中度の平均」のそれぞれについて検出範囲が決定されている。そして、何れのプロファイルにおいても検出範囲内に位置する１次候補（その変量データを△のマーカで示す）が脳動脈瘤候補として２次検出される。
【００５９】
このようにして２次検出が行われると、この２次検出候補についてさらに３次検出が行われる（ステップＳ７）。３次検出では、３つの特徴量を用いて判別分析が行われる。判別分析の手法としては、マハラノビスの距離、主成分分析、線形判別関数等、何れのものも適用可能であるが、２次検出時の手法とは異なるものを適用するものとする。
【００６０】
そして、３次検出された候補領域が最終的な脳動脈瘤候補であるとしてその検出結果が出力される（ステップＳ８）。
図１０に、検出結果として表示部１３に表示出力される検出結果例を示す。
図１０に示すように、表示部１３では、３次元ＭＲＡ画像から作成されたＭＩＰ画像において、３次検出された脳動脈瘤の候補領域を指し示すマーカ情報（図１０中の矢印のマーカ）が表示される。このような表示を行うことにより、脳動脈瘤候補領域と他の領域とを識別可能としている。ＭＩＰ画像は、３次元ＭＲＡ画像データにＭＩＰ処理を施して作成される２次元画像であり、画像中の構造物を３次元的に表示することが可能である。ＭＩＰ処理は最大輝度投影法と呼ばれ、ある方向から平行光線によって投影を行い、ボクセル中の最大の輝度（信号値）を投影面に反映させて３次元的な観察を可能とする２次元画像を作成する処理である。
【００６１】
なお、脳動脈瘤候補領域について算出されたベクトル集中度の平均値等、脳動脈瘤候補の検出に関する情報を出力し、医師の診断時における参考情報としてもよい。また、ベクトル集中度の程度に応じてマーカ情報の色を変えることとしてもよい。例えば、ベクトル集中度が０．８以上であれば赤、０．７〜０．８であれば黄、０．７〜０．５であれば青等のように矢印のマーカの色を変えることにより、ベクトル集中度、つまり瘤の球状の度合いが強いことを医師は視覚的に容易に把握することができる。
【００６２】
なお、図１０ではマーカ情報により脳動脈瘤候補の位置を識別可能とする例を示したが、脳動脈瘤候補領域が他の領域と識別可能に表示されるのであればこのような表示方法に限らない。例えば、ＭＩＰ処理ではなく、ボリュームレンダリング法を用いて３次元的に表示することとしてもよい。ボリュームレンダリング法は、部分領域毎にそれらのボクセルに色情報と不透明度を与えることにより３次元表示を行う手法であり、注目領域については不透明度を高く設定し、それ以外の領域については低く設定することによって注目領域を浮き立たせることができるものである。よって、表示時には領域毎に不透明度、その不透明度に応じた色情報の設定処理を行う。
【００６３】
また、ＭＩＰ画像等ではなく、ベクトル集中度フィルタにより得られたフィルタ処理画像（図６参照）を用いて検出結果を示すこととしてもよい。このとき、例えばベクトル集中度が低い値の領域は青色を、高い値の領域は赤色を付す等して色分けすることにより、医師が算出されたベクトル集中度を視覚的に把握できるようにしてもよい。さらに、ベクトル集中度によって血管領域の色を変えたフィルタ処理画像を前記ＭＩＰ画像の対応する位置に重ね合わせて表示することとしてもよい。このように、医師の脳動脈瘤の検出にあたっての参考情報としてベクトル集中度の情報を提供することが可能である。
【００６４】
次に、図１１を参照して血管部位判別処理について説明する。
この血管部位判別処理は、制御部１１と記憶部１５に記憶された血管部位判別処理の処理プログラムとの協働により実現されるソフトウェア処理である。血管部位判別処理では、頭部を撮影した３次元ＭＲＡ画像上に表れる血管像についてその血管像に含まれる一又は複数の血管部位を判別する。
【００６５】
ＭＲＡはＭＲＩの血管撮像方法の１種である。ＭＲＩでは、被写体の足から頭の方向（この方向を体軸という）に傾斜磁場をかけることによって特定のスライス（断層）にのみエネルギーを吸収させることが可能である。撮像時にはスライス内の血管ではＲＦパルスで血液が飽和された状態となっているが、血管は常に流れがあるため、経時により非飽和の血流が流入するとそのスライスにおける信号強度が増加する。ＭＲＡはこの高信号を画像化することにより血流のある血管を画像化する方法である。
【００６６】
最初に、血管部位の判別に必要な参照画像について説明する。
参照画像は、図１２Ａに示すように、３次元ＭＲＡ画像上の血管像について、一又は複数の血管部位の位置及びその名称が予め設定されたものである。ここで、血管部位とは解剖学上の血管の分類をいい、血管部位の位置とは当該血管部位に属するボクセルの位置をいう。
【００６７】
図１２Ａでは、血管像に含まれる８つの血管部位（前大脳動脈、右中大脳動脈、左中脳動脈、右内頸動脈、左内頸動脈、右後大脳動脈、左後大脳動脈、脳底動脈）についてその位置及び血管部位の名称を設定した例を示している。なお、図１２Ａでは８つの血管部位のうち、３つの血管部位（右中大脳動脈、前大脳動脈、脳底動脈）についてのみ名称を示しているが、８つ全ての血管部位について名称は設定される。
【００６８】
参照画像ｇ２は、参照画像用に選択された、図１２Ｂに示すような頭部ＭＲＡ画像ｇ１の３次元データから作成される。まず、この３次元データからある間隔毎にアキシャル画像（体軸に垂直な面でボクセルを切り出した２次元断層画像）を作成する。そして、一のアキシャル画像において医師の指摘に基づき、手動操作により各血管部位に属するボクセルを指定し、さらにその血管部位の名称を指定する。これを体軸方向に位置を変えてスライスした各アキシャル画像について繰り返すことにより、３次元データを構成する全ボクセルのうち、各血管部位に属するボクセルの位置及び血管部位の名称を設定することができる。
【００６９】
また、参照画像ｇ２には血管の屈曲点、終局点、血管部位同士の交差点等の特徴的な箇所においてランドマークのボクセルが設定される。ランドマークは対象画像と参照画像との位置合わせに使用されるものであるが、詳細な説明は後述する。ランドマークについても医師の指摘に基づく手動操作に応じて設定される。
【００７０】
以上のようにして作成された参照画像ｇ２は、記憶部１５に保存される。
なお、参照画像ｇ２の作成は医用画像処理装置１０の制御部１１で行ってもよいし、外部で作成されたものを記憶部１５に保存することとしてもよい。また、血管像に含まれる８つの血管部位を示すため、図１２Ａでは各血管部位を識別表示したが、実際の参照画像ｇ２は背景が黒（低信号値）、血管像が白（高信号値）と２値化された画像である。そして、各血管部位に属するボクセルの位置情報、名称の情報及びランドマークであるボクセルの位置情報は、参照画像に付帯されている、或いは参照画像と対応付けて別ファイルとして記憶部１５に保存されている。
【００７１】
次に、上記参照画像を用いた血管部位判別処理について具体的に説明する。
医用画像処理装置１０では、まず判別対象の３次元ＭＲＡ画像（以下、対象画像という）に対して制御部１１により正規化処理が施される（ステップＳ１１）。
対象画像として用いられるＭＲＡ画像は、血液の流れを画像化したものであるため、被写体や撮影条件によってボクセルが等サイズではない直方体となったり、ボクセル値の最大値、最小値にばらつきが生じたりする。そこで、対象画像に関する前提条件を統一するため、正規化処理を施す。
【００７２】
正規化処理では、まずボクセルを構成する全ての辺が等サイズとなるように線形補間法により対象画像が変換される。次に、対象画像の全てのボクセルのボクセル値についてヒストグラムが作成され、ヒストグラムの上位５％以上のボクセル値を１０２４、最小のボクセル値を０として、対象画像の全てのボクセル値が０〜１０２４の階調に線形変換される。このとき、ボクセル値が高信号値であるほど濃度値１０２４に近く、低信号値であるほど濃度値０に近くなるように変換される。なお、濃度階調のレンジは、０〜１０２４に限らず適宜設定可能である。
【００７３】
図１３Ａ、図１３Ｂに、正規化処理の一例を示す。
図１３Ａに示す対象画像ｇ３、図１３Ｂに示す対象画像ｇ４はそれぞれ異なる患者を被写体としたものである。そのため、対象画像ｇ３から得られたヒストグラムｈ１（図１３Ａ参照）と、対象画像ｇ４から得られたヒストグラムｈ３（図１３Ｂ参照）とでは、２つの極大点があるという特徴が共通するが、そのボクセル値の範囲にかなり相違があり、全体としてヒストグラム特性が異なるものとなっていることが分かる。このようなヒストグラム特性を有する対象画像ｇ３、ｇ４について上記の正規化処理を施した後に再度ヒストグラムを作成すると、図１３Ａに示すヒストグラムｈ２、図１３Ｂに示すヒストグラムｈ４がそれぞれ得られる。ヒストグラムｈ２、ｈ４から分かるように、正規化処理によって各対象画像ｇ３、ｇ４のヒストグラム特性がほぼ同じものとなっている。
【００７４】
正規化処理を終えると、制御部１１では当該正規化された対象画像から血管像が抽出される（ステップＳ１２）。
まず、対象画像について閾値処理が施され、２値化が行われる。一般的に、ＭＲＡ画像では、図１４Ａに示すように血管像は白く、その他の組織部分は黒っぽく表れるため、２値化画像では血管像はその他の領域とは異なる値となる。よって、領域拡張法により血管像と同程度の信号値を有する領域の抽出を行う。
【００７５】
領域拡張法では、２値化画像を用いて始点となるボクセル（最も白く高濃度値のボクセル）を決定し、２値化処理前の対象画像においてその始点と決定されたボクセルの近傍２６ボクセルを調べ、ある判定条件（例えば、濃度値が５００以上であること）を満たす近傍ボクセルを血管像と判断する。そして、この血管像と判断された近傍ボクセルについても上記と同様の処理を繰り返す。このように、領域を拡張しながら判定条件を満たすボクセルを順次抽出することにより、血管像の画像領域を抽出することができる。
【００７６】
図１４Ｂに、血管像を抽出した血管抽出画像ｇ６を示す。
血管抽出画像ｇ６は、図１４Ａに示す正規化後の対象画像ｇ５から血管像を抽出し、当該血管像の領域を白（濃度値１０２４）、その他の領域を黒（濃度値０）で２値化したものである。
【００７７】
次に、制御部１１では、血管抽出画像の血管像の位置と参照画像の血管像の位置とを略一致させるため、各画像の重心位置を元に位置合わせが行われる（ステップＳ１３）。重心位置は、血管像に属する全てのボクセルの重心となるボクセルの位置である。
【００７８】
図１５Ａ、図１５Ｂを参照して、具体的に説明する。
図１５Ａは、位置合わせ前の血管抽出画像と参照画像とを重ね合わせた図である。図１５Ａから単に血管抽出画像と参照画像を合わせただけではそれぞれの血管像の位置が一致していないことが分かる。
そこで、制御部１１では、図１５Ａに示すように血管抽出画像の重心Ｐ（ｘ１、ｙ１、ｚ１）、参照画像の重心Ｑ（ｘ２、ｙ２、ｙ３）の位置が求められる。次いで、この重心位置Ｐ、Ｑが一致するように、血管抽出画像又は参照画像が平行移動される。平行移動により各重心位置Ｐ、Ｑを一致させた結果が、図１５Ｂに示す図である。図１５Ｂから血管抽出画像の血管像と参照画像の血管像の位置が大まかに一致していることが分かる。
【００７９】
さらに、精度よく位置合わせを行うため、制御部１１では血管抽出画像について剛体変形が行われる（ステップＳ１４）。
まず、剛体変形の前処理として相互相関係数を用いた対応点の検索が行われる。これは、位置合わせを行う２つの画像についてそれぞれ複数の対応点を設定し、この２つの画像において設定された対応点がそれぞれ一致するように一方の画像を剛体変形するためである。ここでは、参照画像において予め定められているランドマークのボクセルと、局所的に画像特性が類似する血管抽出画像のボクセルが対応点として設定される。画像特性の類似性は、血管抽出画像と参照画像について相互相関係数が求められ、この相互相関係数に基づいて判断される。
【００８０】
具体的には、図１６Ａに示すように、予め参照画像ｇ７の血管像において設定されている１２点のランドマークに対応する対応点が血管抽出画像から検索される。対応点の検索時には、図１６Ｂに示すように、血管抽出画像ｇ８において、参照画像ｇ７の各ランドマークと対応する位置のボクセルを開始点とし、血管抽出画像ｇ８及び参照画像ｇ７において当該開始点及びランドマークのボクセルからＸ軸、Ｙ軸、Ｚ軸方向に−１０〜＋１０ボクセルの範囲（２１×２１×２１ボクセルの立方領域）内のボクセルが探索され、各ボクセルについて、下記式２により相互相関係数Ｃ（以下、相関値Ｃという）が算出される。
【００８１】
【数６】
2008104886.tif 000008
上記式２においてＡ（ｉ，ｊ，ｋ）は参照画像ｇ７のボクセル位置、Ｂ（ｉ，ｊ，ｋ）は血管抽出画像ｇ８のボクセル位置を示す。ＩＪＫは探索領域のサイズを示し、ＩＪＫ＝２１×２１×２１である。
【００８２】
また、α、βは、それぞれ参照画像ｇ７、血管抽出画像ｇ８における探索領域内のボクセル値の平均値であり、下記式３、４により示される。σＡ、σＢは、それぞれ参照画像ｇ７、血管抽出画像ｇ８における探索領域内のボクセル値の標準偏差であり、下記式５、６により示される。
【数７】
2008104886.tif 000009
【００８３】
相関値Ｃは−１．０〜１．０の値域を持ち、最大値１．０に近いほど、参照画像ｇ７と血管抽出画像ｇ８の画像特性が類似していることを示す。
そこで、最も大きな相関値Ｃをとるボクセルの位置が、参照画像ｇ７のランドマークに対応する血管抽出画像ｇ８の対応点として設定される。
【００８４】
対応点が設定されると、制御部１１ではこの対応点に基づき、血管抽出画像ｇ８に対して剛体変形を施すことにより、血管抽出画像ｇ８の血管像と参照画像ｇ７の血管像との位置合わせが行われる。剛体変形は、アフィン変換の一つであり、回転、平行移動により座標変換を行うものである。位置合わせは、最小二乗法を用いた剛体変形を複数回繰り返すＩＣＰ（Iterative Closest Point）アルゴリズムにより血管抽出画像ｇ８の対応点が参照画像ｇ７のランドマークに一致するように行われる。このアルゴリズムによれば、剛体変形を行う毎に、参照画像のランドマークと血管抽出画像の対応点における距離の最小二乗誤差を算出し、当該最小二乗誤差がある閾値を超える等の終了条件を満たすまで剛体変形が繰り返される。なお、剛体変形をアフィン変換（拡大縮小、回転、並行移動による座標変換）に代えることとしてもよい。
【００８５】
次いで、参照画像に基づき、剛体変形により位置合わせされた血管抽出画像の血管像に含まれる各血管部位が制御部１１において判別される（ステップＳ１５）。
まず、参照画像の各血管部位に属するボクセル全てを対象に（これを対象ボクセルという）、血管抽出画像におけるある注目ボクセルとのユークリッド距離の２乗が求められる。そして、その求めたユークリッド距離が最短となる対象ボクセルが属する血管部位が、注目ボクセルが属する血管部位であると判断される。またこのとき、対象ボクセルに設定されている血管部位の名称から、注目ボクセルの血管部位の名称が判断される。
【００８６】
血管抽出画像の血管像を構成する全てのボクセルに対し、注目ボクセルを設定して上記の処理を行い、全ボクセルについて対応する血管部位が判別されると、その判別された血管部位に属するボクセルの位置、血管部位の名称を示す血管部位情報が制御部１１により生成され、対象画像に付帯される（ステップＳ１６）。例えば、位置（ｘ３、ｙ３、ｚ３）のボクセルは前大脳動脈の血管部位であると判別された場合、そのボクセル位置「（ｘ３，ｙ３、ｚ３）」のボクセルは血管名称「前大脳動脈」の血管部位であることを示す血管部位情報が対象画像のヘッダ領域に書き込まれる等して付帯されることとなる。
【００８７】
図１７Ａ、図１７Ｂに、血管部位の判別結果を示す。
図１７Ａに示す血管抽出画像ｇ９と図１７Ｂに示す血管抽出画像ｇ１１は、それぞれ異なる被写体から得られた画像であり、図１７Ａに示す画像ｇ１０、図１７Ｂに示す画像ｇ１２はそれぞれ血管抽出画像ｇ９、ｇ１１から各血管部位が判別され、その血管部位毎に色を変えて識別表示した画像である。画像ｇ１０、ｇ１２から、剛体変形による位置合わせを行うことにより、画像ｇ１０、ｇ１２の血管像が異なる形態（血管の位置、大きさ、延在方向等）であるにも拘わらず、同じように血管部位を識別できていることが分かる。
【００８８】
以上が、対象画像における血管部位を判別し、その血管部位情報を付帯するまでの流れである。
そして、このような対象画像について操作部１２を介して表示指示操作がなされると、制御部１１により対象画像にＭＩＰ処理が施されてＭＩＰ画像が生成され、表示部１３上に表示される。以下、ＭＩＰ画像の表示をＭＩＰ表示という。
【００８９】
ＭＩＰ法はある方向から平行光線によって投影を行うとともに、この投影線上にあるボクセル中の最大の輝度（ボクセル値）を投影面に反映させて２次元画像を作成する方法である。この投影方向は医師が観察を所望する視線方向となる。ここでは、医師が自由にその観察方向を操作することができる構成とし、操作部１２を介して指示操作された観察方向に応じて制御部１１により対象画像からＭＩＰ画像が作成され、表示部１３に表示されるものとする。
【００９０】
対象画像がＭＩＰ表示された状態においてさらに血管部位の識別表示の指示操作がなされると、制御部１１では対象画像に付帯されている血管部位情報が参照され、この血管部位情報に基づき、ＭＩＰ画像中の血管像において各血管部位が識別可能となるように表示制御が行われる（ステップＳ１７）。
【００９１】
例えば、血管部位情報に基づいて各血管部位のボクセル位置及び名称が判別されると、制御部１１では、ＭＩＰ画像において前大脳動脈の血管部位に属するボクセルには青色、脳底動脈の血管部位に属するボクセルには緑色等、血管部位毎に判別された位置にあるボクセルに血管部位毎の色が設定される。そして、対象画像のＭＩＰ画像に当該設定色を反映させる。さらにその血管部位の名称を指し示すアノテーション画像が作成されてＭＩＰ画像の対応する血管部位に合成される。
【００９２】
図１８に、識別表示例を示す。
図１８においてＭＩＰ画像ｇ１３は対象画像が頭部上方向からＭＩＰ表示されたものである。このＭＩＰ画像ｇ１３が表示された状態で各血管部位の識別表示が指示操作されると、識別表示画像ｇ１４が表示される。識別表示画像ｇ１４は、血管抽出画像において８種類の血管部位にそれぞれ異なる色を付す等して、各血管部位を識別表示したものである。
識別表示画像ｇ１４において、医師によりある血管部位が選択操作されると、制御部１１の表示制御により当該選択された血管部位に関連づけて「脳底動脈」等の血管部位の名称を示すアノテーションｍが表示される。
【００９３】
上記頭部上方向のＭＩＰ画像ｇ１３に対して、側面方向からのＭＩＰ表示が指示されると、制御部１１により側面方向に応じたＭＩＰ画像ｇ１５が作成されて、表示される。ＭＩＰ画像ｇ１５では、手前側の血管像と後方側の血管像とが重なり合っており、観察しづらい表示となっている。このようなＭＩＰ画像ｇ１５においても、血管部位の識別表示が可能である。血管部位情報を参照することにより、何れのボクセルが何れの血管部位に該当するのかを判別することができるため、ＭＩＰ表示の観察方向の変動によらず、識別表示すべきボクセルを特定できるからである。
【００９４】
ＭＩＰ画像ｇ１５に対応する識別表示画像は画像ｇ１６である。この識別表示画像ｇ１６においては、何れかの血管部位のみを抽出して表示することが可能である。
制御部１１では、操作部１２を介して何れか一の血管部位が選択操作されると、当該選択された血管部位のボクセルの輝度のみを投影させたＭＩＰ画像、つまり対象画像のＭＩＰ画像ｇ１５から選択された血管部位のみが抽出された血管選択画像ｇ１７が表示される。血管選択画像ｇ１７では、選択された血管部位のみＭＩＰ表示され、その他の血管部位が非表示とされているため、医師は関心のある血管部位のみに注目して観察することができる。
【００９５】
また、これらの表示画像ｇ１３〜ｇ１７は、同一画面上に並べて表示することとしてもよいし、一画面一画像として切替表示することとしてもよい。前者の場合、ＭＩＰ表示画像ｇ１３と識別表示画像ｇ１４、血管選択画像ｇ１７等とを比較観察することが可能となるし、後者の場合は全画面表示で各画像ｇ１３〜ｇ１７を観察することができ、詳細を観察しやすくなる。
【００９６】
以下、検出処理の実施例を示す。
２０人の患者についての３次元ＭＲＡ画像データを得た。これらの画像データは、そのマトリクスサイズが２５６×２５６であり、空間分解能は０．６２５〜０．７８ｍｍ、スライス厚は０．５〜１．２ｍｍである。これらの画像データには７つの未破裂脳動脈瘤が含まれていることが分かっている。未破裂脳動脈瘤は、経験のある脳神経外科医により決定されたものである。
【００９７】
まず、これら３次元ＭＲＡ画像データを、線形補完法を用いて等ボクセルサイズの画像データに変換し、正規化を行った。この正規化により、全ての画像データは、マトリクスサイズが４００×４００×２００、空間分解能が０．５ｍｍの等ボクセルの画像データとなった。次に、階調変換処理を施し、濃度階調を０〜１０２４に線形変換した後、医師により決定された脳動脈瘤の領域（真陽性）について大きさ（体積）、球形度、ベクトル集中度の各特徴量が算出されるとともに、正常な血管領域（偽陽性）についても同様の特徴量が算出される。そして、これら特徴量を２次検出の識別器であるルールベース法の検出範囲の決定、判別分析の教師データとして使用した。同様に、３次検出時の識別器における教師画像として用いて識別器の学習を行った。
【００９８】
上記のように構築された識別器により、２０症例の３次元ＭＲＡ画像データについて図７の検出処理を行ったところ、１症例あたり１．８５個の偽陽性候補が含まれていたが、脳動脈瘤の検出正答率は１００％であり、精度良く検出できることが確認された。
【００９９】
以上のように、第１実施形態によれば、ベクトル集中度フィルタを用いて中心部分に勾配ベクトルが集中する特徴を有する脳動脈瘤候補を精度良く検出し、その検出情報を医師に提供することができる。よって、医師の読影作業時の疲労や見落としを防止することができ、診断精度の向上が期待される。
【０１００】
また、脳動脈瘤候補領域は他の領域と識別可能に表示するので、画像上において医師が検出結果を容易に把握することができ、読影作業の円滑化を図ることができる。
【０１０１】
また、ベクトル集中度フィルタは、ＭＲＡ画像そのものではなく血管領域を抽出した抽出画像に適用するため、フィルタ処理に要する処理時間の短縮化を図ることができる。
【０１０２】
なお、上記実施形態は、本発明を適用した好適な一例であり、これに限定されない。
例えば、上記説明では、３次元ＭＲＡ画像を用いて脳動脈瘤候補を検出していたが、２次元ＭＲＡ画像を用いて検出を行うこととしてもよい。この場合、脳動脈瘤候補の大きさはピクセル数、球形度は円形度として２次元の特徴量を算出することとなる。
【０１０３】
また、ＭＲＡ画像を用いて動脈瘤候補の検出を行う例を説明したが、造影剤を用いて血管領域を撮影した造影ＭＲＡ画像等、他の撮影方法によるＭＲＩ画像を用いてもよいし、ＣＴＡ（Computed Tomography Angiography）、ＤＳＡ（Digital Subtraction Angiography）等、他の撮影装置により血管領域が撮影された画像を用いることとしてもよい。また、検出対象も動脈瘤のみならず、球状の瘤の特徴を有する病変部であれば本発明を適用して検出することが可能である。
【０１０４】
また、上記実施形態によれば、画像上の血管像に含まれる８つの血管部位の位置及び名称を予め定めた参照画像と対象画像とを、その血管像について位置合わせすることにより、対象画像上の血管像に含まれる各血管部位の位置及び名称を判別する。そして、この位置及び名称の情報を血管部位情報として対象画像に付帯させるので、対象画像のＭＩＰ表示を行う際には当該血管部位情報に基づき、各血管部位を容易に判別することができ、各血管部位の識別表示が可能となる。従って、医師は対象画像における特定の血管部位に注目して対象画像を観察することが可能となり、読影効率の向上を図ることができる。
【０１０５】
ＭＩＰ表示を行う際、そのＭＩＰ表示を行う方向によっては複数の血管部位が重複し、注目したい血管部位に対する観察が困難となる場合がある。
しかしながら、本実施形態によれば各血管部位の識別表示が可能となるため、医師は各血管部位についての位置や名称を容易に特定することができ、読影作業の効率化を図ることができる。
【０１０６】
また、識別表示された各血管部位のうち、任意の血管部位の選択操作に応じて、当該選択された血管部位のみを抽出して表示したターゲットＭＩＰ画像を作成し、表示するので、医師は注目したい血管部位のみを観察することが可能となる。よって、複数の血管部位の重複を解消することができ、動脈瘤の多発部位等、ある血管部位に特定して観察を行うことができる。
【０１０７】
また、対象画像と参照画像の血管像の位置合わせを行う際には、まず重心位置に基づいて大まかな位置合わせを行った後、相互相関係数により対象画像に剛体変形を施して対象画像の血管像の特徴点を参照画像の血管像の対応する特徴点に一致するように細かな位置合わせを行う。
【０１０８】
異なる被写体（患者）であっても主要な血管部位の形態（血管の長さや延在方向、太さ等）は概ね共通するが、主要でない細い血管部位についてはその形態に個体差が生じるため、血管部位の形態は被写体によって様々なものとなる。
しかしながら、本実施形態のように最終的に剛体変形により血管像の主な屈曲点等の特徴点の位置を合わせこむことにより、対象画像中の各血管部位がどのような形態を呈していても、精度よく参照画像の血管部位と対応させることができる。従って、被写体によらず一律に血管部位を判別することができ、汎用性が高い。
【０１０９】
さらに、位置合わせは２段階（重心位置に基づくものと、剛体変形によるもの）で行うので、血管部位を判別する判別精度を向上させることができる。また、剛体変形の前に重心位置により位置合わせを行っておくことにより、剛体変形に係る処理時間の短縮化を図ることができ、処理効率が良い。
【０１１０】
なお、上述した実施形態は、本発明を適用した好適な一例であり、これに限定されない。
例えば、ＭＲＡ画像を用いた例を説明したが、造影剤を用いて血管を撮影した造影ＭＲＡ画像等、他の撮影方法によるＭＲＩ画像を用いてもよい。また、ＣＴＡ（Computed Tomography  Angiography）、ＤＳＡ（Digital Subtraction Angiography）等、他の撮影装置により血管が撮影された画像を用いることとしてもよい。
【０１１１】
〈第２実施形態〉
第２実施形態では、ＧＣフィルタバンクによる検出処理の例を説明する。
第２実施形態に係る医用画像処理装置は、第１実施形態に係る医用画像処理装置１０と同一構成であり、動作が異なるのみである。よって、第１実施形態に係る医用画像処理装置１０（図１参照）と同一構成部分には同一の符号を付し、以下、第２実施形態の医用画像処理装置１０の動作について説明する。
【０１１２】
図１９は、第２実施形態に係る検出処理を示すフローチャートである。
図１９に示す検出処理では、まずＭＲＡの３次元画像データの入力が行われ（ステップＳ１０１）、この３次元画像データに対して前処理が行われる（ステップＳ１０２）。前処理が終了すると、３次元画像データから血管の画像領域が抽出される（ステップＳ１０３）。なお、このステップＳ１０１〜Ｓ１０３は、第１実施形態において図２を参照して説明したステップＳ１〜３と同様の処理であるため、ここでは詳細な説明を省略する。
【０１１３】
次に、抽出された血管領域の３次元ＭＲＡ画像を用いてＧＣフィルタバンクにより、脳動脈瘤の１次候補領域の検出が行われる（ステップＳ１０４）。
以下、ＧＣフィルタバンクによる検出方法について説明する。
ＧＣフィルタバンクは、各種フィルタ処理を組み合わせたものであり、分析バンクと再構成バンクに分かれている。分析バンクは原画像（血管領域の３次元ＭＲＡ画像）の多重解像度解析を行って解像度レベルの異なる画像（以下、部分画像という）を作成し、この部分画像から重み画像を作成するものである。一方、再構成バンクは重み画像によって各部分画像を重み付けした後、この重み付けられた部分画像から原画像の再構成を行うものである。
【０１１４】
図２０は、分析バンクを示す図である。
図２０に示すように、分析バンクでは血管領域の３次元ＭＲＡ画像を原画像Ｓ０としてフィルタバンクＡ（ｚｊ）においてフィルタ処理が行われ、各解像度レベルｊの部分画像が順次作成される。ここでは、解像度レベル１〜３の部分画像を作成する例を説明する。
フィルタバンクＡ（ｚｊ）は、図２１に示すようにフィルタＨ（ｚｊ）、Ｇ（ｚｊ）によるフィルタ処理を経て、画像Ｓｊ−１を部分画像Ｓｊ、Ｗｚｊ、Ｗｙｊ、Ｗｘｊに分解するものである。ここで、Ｓｊはｘ、ｙ、ｚの各方向においてＳｊ−１に平滑化フィルタＨ（ｚｊ）を施したものである。平滑化フィルタＨ（ｚｊ）は下記式１１により表される。
【数８】
2008104886.tif 000010
なお、上記式１１において用いられているｚは、ｚ変換を示すものである（以下、フィルタを示す式１２〜１４において同じ）。
【０１１５】
部分画像Ｗｚｊ、Ｗｙｊ、Ｗｘｊは、画像Ｓｊ−１のｘ、ｙ、ｚ方向のそれぞれについてフィルタＧ（ｚｊ）が適用されて得られた１階差分画像である。フィルタＧ（ｚｊ）は下記式１２により表される。
【数９】
2008104886.tif 000011
部分画像Ｗｚｊ、Ｗｙｊ、Ｗｘｊは、解像度レベルが小さい（ｊ＝１）ほど、細かな変動に対する検出性が高い。一方、解像度レベルが大きい（ｊ＝３）ほど粗い変動に対する検出性が高くなる。従って、解像度レベルが小さい部分画像Ｗｚｊ、Ｗｙｊ、Ｗｘｊは大きさが小さい動脈瘤の検出に適し、解像度レベルが大きい部分画像Ｗｚｊ、Ｗｙｊ、Ｗｘｊは大きい動脈瘤の検出に適している。
【０１１６】
Ｓｊについてはさらに解像度レベルｊ＋１のフィルタバンクＡ（ｚｊ＋１）によるフィルタ処理が行われることとなる。このように、繰り返しフィルタ処理が行われることにより、解像度レベルｊ（ｊ＝１〜３）毎に部分画像Ｓｊ、Ｗｚｊ、Ｗｙｊ、Ｗｘｊが作成される。
【０１１７】
次に、図２０に示すように、各解像度レベルの部分画像Ｗｚｊ、Ｗｙｊ、Ｗｘｊに対してベクトル集中度フィルタＧＣｊによるフィルタ処理が施され、解像度レベル毎のベクトル集中度が算出される。ベクトル集中度の算出方法については前述したのでここでは説明を省略する。
【０１１８】
算出されたベクトル集中度はニューラルネットワークＮＮに入力される。
ニューラルネットワークＮＮでは、予め教師データを用いて脳動脈瘤の可能性が高いほど１の値を、逆に低ければ０の値で出力するように、０〜１の範囲で出力値を出力するように設計されている。
【０１１９】
ニューラルネットワークＮＮから出力値が得られると、判定部ＮＭにおいて当該出力値に基づいて重み画像Ｖが生成され、出力される。重み画像Ｖはある閾値（ここでは０．８とする）より大きい出力値を有するボクセルについてはボクセル値を１とし、閾値０．８以下の出力値を有するボクセルについてはボクセル値を０とすることにより、作成される。すなわち、閾値である０．８を超えるベクトル集中度が算出されたボクセルは、脳動脈瘤の領域を構成するボクセルである可能性が高い。一方、閾値０．８以下のベクトル集中度を持つボクセルは正常な血管領域を構成するボクセルである可能性が高い。よって、閾値を境界としてボクセル値を２値化することにより、重み画像を作成するものである。重み画像Ｖは再構成バンクに入力される。
【０１２０】
図２２は、再構成バンクを示す図である。
再構成バンクは、各部分画像Ｓｊ、Ｗｚｊ、Ｗｙｊ、Ｗｘｊに重み付け処理を施した後、フィルタバンクＳ（ｚｊ）を経て元の原画像Ｓ０を再構成するものである。
重み付け処理では、各部分画像Ｓｊ、Ｗｚｊ、Ｗｙｊ、Ｗｘｊに重み画像Ｖが乗算される。すなわち、重み画像Ｖにおいてボクセル毎に設定された１又は０の値は、重み付け処理において重み付け係数として用いられることとなる。重み画像Ｖが乗算された各部分画像はフィルタバンクＳ（ｚｊ）に入力される。
【０１２１】
フィルタバンクＳ（ｚｊ）は、図２３に示すようにフィルタＬ（ｚｊ）、Ｋ（ｚｊ）、Ｈ（ｚｊ）Ｌ（ｚｊ）によるフィルタ処理を経て、Ｓｊ、Ｗｚｊ、Ｗｙｊ、ＷｘｊからＳｊ−１を再構成するものである。
ここで、フィルタＬ（ｚｊ）、Ｋ（ｚｊ）は、下記式１３、１４により表される。
【数１０】
2008104886.tif 000012
【０１２２】
すなわち、Ｓｊに対してはｘ、ｙ、ｚ方向についてそれぞれフィルタＬ（ｚｊ）によるフィルタ処理が施される。また、Ｗｚｊについてはｘ方向についてはフィルタＫ（ｚｊ）、ｙ、ｚ方向についてはフィルタＨ（ｚｊ）Ｌ（ｚｊ）が適用される。Ｗｙｊに対してはｙ方向についてフィルタＫ（ｚｊ）、ｚ方向についてフィルタＨ（ｚｊ）Ｌ（ｚｊ）が適用され、Ｗｘｊに対してはｚ方向についてフィルタＫ（ｚｊ）が適用される。そして、それぞれのフィルタ処理後のＳｊ、Ｗｚｊ、Ｗｙｊ、Ｗｘｊを加算することにより、Ｓｊ−１が得られる。
【０１２３】
Ｓｊ−１についてはさらにフィルタバンクＳ（ｚｊ−１）によるフィルタ処理にかけられ、Ｓｊ−２が再構成されることとなる。このようなフィルタ処理を繰り返すことにより元の画像Ｓ０を再構成することができる。
【０１２４】
原画像Ｓ０の再構成は、フィルタバンクＡ（ｚｊ）、Ｓ（ｚｊ）における各フィルタＨ（ｚｊ）、Ｇ（ｚｊ）、Ｌ（ｚｊ）、Ｋ（ｚｊ）が、図２１、図２３に示すフィルタの組み合わせである場合のみ、完全に元の画像を再現できるものである。よって、再構成の結果、重み画像Ｖにおいてボクセル値が１に設定されたボクセルについては完全に元の原画像Ｓ０が再現されるのに対し、重み画像Ｖにおいてボクセル値が０に設定されたボクセルについては原画像Ｓ０が再現されず、そのボクセル値は０となってしまう。従って、ＧＣフィルタバンクにより出力された出力画像Ｓ０は、脳動脈瘤である可能性が高い画像領域のみが再現された画像となる。この出力画像Ｓ０において現れた画像が脳動脈瘤の１次候補である。
【０１２５】
出力画像Ｓ０が得られると、この出力画像Ｓ０を用いて特徴量の算出が行われる（ステップＳ１０５）。特徴量は、第１実施形態と同様に候補領域の大きさ、球形度、領域内の各ボクセルにおけるベクトル集中度の平均値を算出する。算出された特徴量を用いて２次検出が行われた後（ステップＳ１０６）、さらに３次検出が行われ（ステップＳ１０７）、最終的な検出結果として出力される（ステップＳ１０８）。ステップＳ１０５〜Ｓ１０８の処理は、図２を参照して説明したステップＳ５〜Ｓ８と同一であるので、詳細な説明は省略する。
【０１２６】
また、第２実施形態においても血管部位判別処理が行われるが、第１実施形態と同様の処理内容であるので、その処理内容及び効果についての説明は省略する。
【０１２７】
以上のように、第２実施形態によれば、ＧＣフィルタバンクにより頭部画像に対して多重解像度解析を行って得られた部分画像によって解像度レベル毎にベクトル集中度を算出する。解像度レベル毎にベクトル集中度を算出することにより、様々な大きさの動脈瘤の検出に対応することができ、より精度の高い検出処理が可能となる。
【０１２８】
また、ベクトル集中度が閾値以上の所定値となるボクセルのみ元の画像が再現されるようボクセル値を１に設定し、その他のボクセルについてはボクセル値を０に設定した重み画像を生成して各部分画像に乗算して重み付けを行い、この重み付けされた各部分画像から元の画像を再構成する。よって、ベクトル集中度が高い、すなわち動脈瘤である可能性が高い領域のみが再構成されることとなり、動脈瘤の候補領域のみ画像化した再構成画像を得ることができる。このような再構成画像を用いて特徴量の算出を行うことにより、算出する特徴量から候補領域以外の画像要素を排除することができ、候補領域についての特徴量を正確に算出することが可能となる。従って、検出処理自体の精度をさらに向上させることが可能となる。
【図面の簡単な説明】
【０１２９】
【図１】実施の形態に係る医用画像処理装置の内部構成を示す図である。
【図２】医用画像処理装置により実行される検出処理の流れを説明するフローチャートである。
【図３Ａ】ＭＲＡ画像例を示す図である。
【図３Ｂ】血管領域を抽出した抽出画像例を示す図である。
【図４】ベクトル集中度フィルタを示す図である。
【図５】脳動脈瘤モデル、血管領域モデルを示す図である。
【図６】ベクトル集中度フィルタの出力画像例を示す図である。
【図７Ａ】閾値処理前のフィルタ処理画像を示す図である。
【図７Ｂ】閾値処理後のフィルタ処理画像を示す図である。
【図８】球形度の算出方法を説明する図である。
【図９Ａ】ルールベース法による識別方法について説明する図である。
【図９Ｂ】ルールベース法による識別方法について説明する図である。
【図１０】脳動脈瘤候補の検出結果の出力例を示す図である。
【図１１】医用画像処理装置により実行される血管部位判別処理を説明するフローチャートである。
【図１２Ａ】参照画像例を示す図である。
【図１２Ｂ】参照画像の作成に用いた元画像を示す図である。
【図１３Ａ】対象画像及び正規化処理を施す前後における対象画像のヒストグラムを示す図である。
【図１３Ｂ】図１３Ａとは異なる被写体の対象画像及び正規化処理を施す前後における対象画像のヒストグラムを示す図である。
【図１４Ａ】対象画像を示す図である。
【図１４Ｂ】図１４Ａの対象画像から血管を抽出した血管抽出画像を示す図である。
【図１５Ａ】血管抽出画像と参照画像を示す図である。
【図１５Ｂ】図１５Ａに示す血管抽出画像と参照画像とを重ね合わせた図である。
【図１６Ａ】参照画像におけるランドマークを示す図である。
【図１６Ｂ】血管抽出画像における対応点を示す図である。
【図１７Ａ】血管抽出画像とその血管部位の判別結果を示す図である。
【図１７Ｂ】血管抽出画像とその血管部位の判別結果を示す図である。
【図１８】対象画像において判別された各血管部位の識別表示例を示す図である。
【図１９】第２実施形態に係る検出処理を示すフローチャートである。
【図２０】ＧＣフィルタバンクの分析バンクを示す図である。
【図２１】フィルタバンクＡ（ｚｊ）を示す図である。
【図２２】ＧＣフィルタバンクの再構成バンクを示す図である。
【図２３】フィルタバンクＳ（ｚｊ）を示す図である。
【符号の説明】
【０１３０】
１０  医用画像処理装置
１１  制御部
１２  操作部
１３  表示部
１４  通信部
１５  記憶部
１６  病変候補検出部
【図１】
2008104886.tif 000017
【図２】
2008104886.tif 000018
【図４】
2008104886.tif 000019
【図５】
2008104886.tif 000020
【図８】
2008104886.tif 000021
【図９Ａ】
2008104886.tif 000022
【図９Ｂ】
2008104886.tif 000023
【図１１】
2008104886.tif 000024
【図１２Ａ】
2008104886.tif 000025
【図１５Ａ】
2008104886.tif 000026
【図１５Ｂ】
2008104886.tif 000027
【図１６Ａ】
2008104886.tif 000028
【図１６Ｂ】
2008104886.tif 000029
【図１７Ａ】
2008104886.tif 000030
【図１７Ｂ】
2008104886.tif 000031
【図１９】
2008104886.tif 000032
【図２０】
2008104886.tif 000033
【図２１】
2008104886.tif 000034
【図２２】
2008104886.tif 000035
【図２３】
2008104886.tif 000036
【図３Ａ】
2008104886.tif 000037
【図３Ｂ】
2008104886.tif 000038
【図６】
2008104886.tif 000039
【図７Ａ】
2008104886.tif 000040
【図７Ｂ】
2008104886.tif 000041
【図１０】
2008104886.tif 000042
【図１２Ｂ】
2008104886.tif 000043
【図１３Ａ】
2008104886.tif 000044
【図１３Ｂ】
2008104886.tif 000045
【図１４Ａ】
2008104886.tif 000046
【図１４Ｂ】
2008104886.tif 000047
【図１８】
2008104886.tif 000048