(57)【要約】      （修正有）【課題】大動脈などの血管の検出精度を向上させた血管検出装置を提供する。【解決手段】ＭＲ装置１００は、肺と肝臓と大動脈とを含む撮影部位を横切るｍ枚のスライスの画像を作成する画像作成手段９１と、スライスの画像に基づいて、スライスが肝臓を横切るか否かを判断するための指標を求める手段と、指標に基づいて、ｍ枚のスライスの中から、肝臓を横切る第１のスライスを選択する選択手段９５と、第１のスライスから検出された大動脈の位置に基づいて、肺を横切る第２のスライスにおける血管を検出する検出手段９６とを有する。【選択図】図１
2015123306.tif 000002
【特許請求の範囲】
【請求項１】
第１の部位と第２の部位と血管とを含む撮影部位を横切るｍ枚のスライスの画像を作成する画像作成手段と、
前記スライスの画像に基づいて、前記スライスが前記第１の部位を横切るか否かを判断するための指標を求める手段と、
前記指標に基づいて、前記ｍ枚のスライスの中から、前記第１の部位を横切る第１のスライスを選択する選択手段と、
前記第１のスライスにおける前記血管を検出する検出手段であって、前記第１のスライスにおける前記血管の位置に基づいて、前記第２の部位を横切る第２のスライスにおける血管を検出する検出手段と、
を有する血管検出装置。
【請求項２】
前記選択手段は、
前記指標に基づいて、前記ｍ枚のスライスの中から、前記第１の部位を横切る２枚以上のスライスを求め、前記２枚以上のスライスの中から、前記第１のスライスを選択する、請求項１に記載の血管検出装置。
【請求項３】
前記選択手段は、
前記２枚以上のスライスのうち、前記第２の部位に最も近いスライスを、前記第１のスライスとして選択する、請求項２に記載の血管検出装置。
【請求項４】
前記検出手段は、
前記スライスごとに前記血管の検索領域を設定し、前記検索領域の中から前記血管を検出する、請求項１〜３のうちのいずれか一項に記載の血管検出装置。
【請求項５】
前記ｍ枚のスライスには、それぞれスライス番号ｎ＝１，２，３，・・・ｍが対応付けられており、
前記第１のスライスのスライス番号ｎはｎ＝ｊ（１＜ｊ＜ｍ）であり、
前記検出手段は、
前記第１のスライスにおける前記血管を検出する第１のステップと、
前記第１のスライスに対してスライス番号が大きくなる方向に並ぶスライス番号ｎ＝ｊ＋１〜ｍのスライスにおける前記血管を検出する第２のステップと、
前記第１のスライスに対してスライス番号が小さくなる方向に並ぶスライス番号ｎ＝ｊ−１〜１のスライスにおける前記血管を検出する第３のステップと、
を実行する、請求項４に記載の血管検出装置。
【請求項６】
前記検出手段は、
前記第２のステップにおいて、
スライス番号ｎ（ｊ≦ｎ≦ｍ−１）のスライスにおける前記血管の位置に基づいて、スライス番号ｎ＋１（ｊ≦ｎ≦ｍ−１）のスライスに第１の検索領域を設定し、
前記第３のステップにおいて、
スライス番号ｎ（ｊ≧ｎ≧２）のスライスにおける前記血管の位置に基づいて、スライス番号ｎ−１（ｊ≧ｎ≧２）のスライスに第２の検索領域を設定する、請求項５に記載の血管検出装置。
【請求項７】
前記検出手段は、
前記第２のステップにおいて、
スライス番号ｎ（ｊ≦ｎ≦ｍ−１）のスライスにおける前記血管の位置を表す第１の角度に基づいて、前記第１の検索領域を設定し、
前記第３のステップにおいて、
スライス番号ｎ（ｊ≧ｎ≧２）のスライスにおける前記血管の位置を表す第２の角度に基づいて、前記第２の検索領域を設定する、請求項６に記載の血管検出装置。
【請求項８】
前記検出手段は、
前記第２のステップと前記第３のステップとを並列に実行する、請求項５〜７のうちのいずれか一項に記載の血管検出装置。
【請求項９】
前記指標を求める手段は、
前記第１の部位および前記第２の部位が第１の論理値を有し、背景領域が第２の論理値を有するように、前記ｍ枚のスライスの各々の画像を２値化する第１の２値化手段と、
前記第１の部位が前記第１の論理値を有し、前記第２の部位および前記背景領域が前記第２の論理値を有するように、前記ｍ枚のスライスの各々の画像を２値化する第２の２値化手段と、
前記第１の２値化手段により得られた第１の２値画像の特徴量と、前記第２の２値化手段により得られた第２の２値画像の特徴量とに基づいて、前記指標を算出する算出手段と、
を有する、請求項１〜８のうちのいずれか一項に記載の血管検出装置。
【請求項１０】
前記第１の２値画像の特徴量は、前記第１の２値画像における前記第１の論理値の面積であり、
前記第２の２値画像の特徴量は、前記第２の２値画像における前記第１の論理値の面積である、請求項９に記載の血管検出装置。
【請求項１１】
前記算出手段は、
前記第１の２値画像における前記第１の論理値の面積と、前記第２の２値画像における前記第１の論理値の面積とを求め、これらの面積の比に基づいて、前記指標を算出する、請求項１０に記載の血管検出装置。
【請求項１２】
前記指標を求める手段は、
前記第１の部位および前記第２の部位が第１の論理値を有し、背景領域が第２の論理値を有するように、前記ｍ枚のスライスの各々の画像を２値化する第１の２値化手段と、
前記第１の部位が前記第１の論理値を有し、前記第２の部位および前記背景領域が前記第２の論理値を有するように、前記ｍ枚のスライスの各々の画像を２値化する第２の２値化手段と、
前記第１の２値画像から所定の部位を削除する第１の削除処理と、前記第２の２値画像から前記所定の部位を削除する第２の削除処理とを実行する削除手段と、
前記第１の削除処理後の第１の２値画像の特徴量と、前記第２の削除処理後の第２の２値画像の特徴量とに基づいて、前記指標を算出する算出手段と、
を有する、請求項１〜８のうちのいずれか一項に記載の血管検出装置。
【請求項１３】
前記第１の削除処理後の第１の２値画像の特徴量は、前記第１の削除処理後の第１の２値画像における前記第１の論理値の面積であり、
前記第２の削除処理後の第２の２値画像の特徴量は、前記第２の削除処理後の第２の２値画像における前記第１の論理値の面積である、請求項１２に記載の血管検出装置。
【請求項１４】
前記算出手段は、
前記第１の削除処理後の第１の２値画像における前記第１の論理値の面積と、前記第２の削除処理後の第２の２値画像における前記第１の論理値の面積とを求め、これらの面積の比に基づいて、前記指標を算出する、請求項１３に記載の血管検出装置。
【請求項１５】
前記第１の２値化手段は、領域拡張法を用いた２値化であり、
前記第２の２値化手段は、閾値処理を用いた２値化である、請求項９〜１４のうちのいずれか一項に記載の血管検出装置。
【請求項１６】
前記第１の部位は肝臓を含み、前記第２の部位は肺を含む、請求項１〜１５のうちのいずれか一項に記載の血管検出装置。
【請求項１７】
請求項１〜１６のうちのいずれか一項に記載の血管検出装置を有する医用装置。
【請求項１８】
第１の部位と第２の部位と血管とを含む撮影部位を横切るｍ枚のスライスの画像を作成する画像作成処理と、
前記スライスの画像に基づいて、前記スライスが前記第１の部位を横切るか否かを判断するための指標を求める処理と、
前記指標に基づいて、前記ｍ枚のスライスの中から、前記第１の部位を横切る第１のスライスを選択する選択処理と、
前記第１のスライスにおける前記血管を検出する検出処理であって、前記第１のスライスにおける前記血管の位置に基づいて、前記第２の部位を横切る第２のスライスにおける血管を検出する検出処理と、
を計算機に実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、血管を検出する血管検出装置、被検体の画像を取得し、取得した画像から血管を検出する磁気共鳴装置、および血管を検出するためのプログラムに関する。
【背景技術】
【０００２】
従来より、造影剤を用いて被検体を撮影する磁気共鳴装置が知られている（特許文献１）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００９−２６１９０４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
造影剤を用いて被検体を撮影する方法の一例として、オペレータが、画像の中から大動脈の位置を見つけて、大動脈の位置に、造影剤を検出するためのトラッカー領域を設定する方法がある。この方法では、トラッカー領域に所定量の造影剤が到達したときに、撮影が実行される。しかし、経験の浅いオペレータは、大動脈の位置を見つけるのに長時間を要したり、大動脈とは別の場所にトラッカー領域を設定しまうことがある。したがって、大動脈を短時間で自動的に検出する手法の開発が望まれている。大動脈を短時間で自動的に検出する手法の一例が、論文“Med Imag Tech Vol.31 No.2 March 2013”で発表されている。
【０００５】
上記の論文では、勾配磁場のＺdephaserの面積を大きくすることより、大動脈内の血液の信号をできるだけ低減させている。したがって、各スライスから低信号の領域を検出することにより、スライスごとに大動脈を検出することができる。しかし、肺を横切るスライスでは、大動脈だけでなく、肺に含まれる空気の領域も低信号となるので、肺に含まれる空気の領域が大動脈として誤検出されることがある。したがって、大動脈などの血管の検出精度を向上させるための技術が望まれている。
【課題を解決するための手段】
【０００６】
本発明の第１の観点は、第１の部位と第２の部位と血管とを含む撮影部位を横切るｍ枚のスライスの画像を作成する画像作成手段と、
前記スライスの画像に基づいて、前記スライスが前記第１の部位を横切るか否かを判断するための指標を求める手段と、
前記指標に基づいて、前記ｍ枚のスライスの中から、前記第１の部位を横切る第１のスライスを選択する選択手段と、
前記第１のスライスにおける前記血管を検出する検出手段であって、前記第１のスライスにおける前記血管の位置に基づいて、前記第２の部位を横切る第２のスライスにおける血管を検出する検出手段と、
を有する血管検出装置である。
【０００７】
本発明の第２の観点は、上記の血管検出装置を有する磁気共鳴装置である。
【０００８】
本発明の第３の観点は、第１の部位と第２の部位と血管とを含む撮影部位を横切るｍ枚のスライスの画像を作成する画像作成処理と、
前記スライスの画像に基づいて、前記スライスが前記第１の部位を横切るか否かを判断するための指標を求める処理と、
前記指標に基づいて、前記ｍ枚のスライスの中から、前記第１の部位を横切る第１のスライスを選択する選択処理と、
前記第１のスライスにおける前記血管を検出する検出処理であって、前記第１のスライスにおける前記血管の位置に基づいて、前記第２の部位を横切る第２のスライスにおける血管を検出する検出処理と、
を計算機に実行させるためのプログラムである。
【発明の効果】
【０００９】
第２のスライスにおいて第２の部位を血管として誤検出する可能性を低減することができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の第１の形態の磁気共鳴装置の概略図である。
【図２】本形態で実行されるスキャンを示す図である。
【図３】撮影部位を概略的に示す図である。
【図４】ローカライザスキャンＬＳおよび本スキャンＭＳを実行するときのフローの一例を示す図である。
【図５】ローカライザスキャンＬＳの説明図である。
【図６】識別器Ｃを作成するために使用された教師データを概略的に示す図である。
【図７】アキシャル画像Ｄ１に設定された検索領域Ｒ１を示す図である。
【図８】検索領域Ｒ１に設定されたウィンドウＷを示す図である。
【図９】識別器Ｃを用いて、大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを求めるときの説明図である。
【図１０】アキシャル画像Ｄ１の大動脈ＡがウィンドウＷの先端部に内接しているときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）を示す図である。
【図１１】アキシャル画像Ｄ２に設定された検索領域Ｒ２を示す図である。
【図１２】検索領域Ｒ２から大動脈Ａが検出されたときの様子を概略的に示す図である。
【図１３】アキシャル画像Ｄ３に設定された検索領域Ｒ３を示す図である。
【図１４】検索領域Ｒ３から検出された大動脈Ａの位置を示す図である。
【図１５】スライスＳＬ１〜ＳＬ１０から検出された大動脈Ａの位置を表す回転角θ＝θ１〜θ１０を概略的に示す図である。
【図１６】誤検出の例を示す図である。
【図１７】スライスＳＬ２に設定された検索領域Ｒ２を示す図である。
【図１８】本形態におけるステップＳＴ３のフローを示す図である。
【図１９】選択されたスライスＳＬ１のアキシャル画像Ｄ１を概略的に示す図である。
【図２０】シード領域Ｒｓｅｅｄの一例を示す図である。
【図２１】領域拡張法でスライスＳＬ１の画像Ｄ１を２値化することにより得られた２値画像ＤＡ１を概略的に示す図である。
【図２２】通常の閾値処理でスライスＳＬ１の画像Ｄ１を２値化することにより得られた２値画像ＤＢ１を概略的に示す図である。
【図２３】ステップＳＴ３４０で実行される処理を説明するための図である。
【図２４】スライスＳＬ３の画像に対してステップＳＴ３０を実行するときの説明図である。
【図２５】スライスＳＬ７の画像に対してステップＳＴ３を実行するときの説明図である。
【図２６】各スライスと面積比Ｒとの関係を表すプロファイルを示す図である。
【図２７】スライスＳＬ１〜ＳＬ１０の中から、肝臓側に位置するスライスを求めるときの説明図である。
【図２８】アキシャル画像Ｄ５に設定された検索領域Ｒ５を示す図である。
【図２９】検索領域Ｒ５から大動脈Ａが検出されたときの様子を概略的に示す図である。
【図３０】アキシャル画像Ｄ６に設定された検索領域Ｒ６を示す図である。
【図３１】検索領域Ｒ６から検出された大動脈Ａの位置を示す図である。
【図３２】アキシャル画像Ｄ７に設定された検索領域Ｒ７を示す図である。
【図３３】検索領域Ｒ７から検出された大動脈Ａの位置を示す図である。
【図３４】スライスＳＬ６〜ＳＬ１０から検出された大動脈Ａの位置を表す回転角θ＝θ６〜θ１０を概略的に示す。
【図３５】アキシャル画像Ｄ４に設定された検索領域Ｒ４を示す図である。
【図３６】検索領域Ｒ４から検出された大動脈Ａの位置を示す図である。
【図３７】アキシャル画像Ｄ３に設定された検索領域Ｒ３を示す図である。
【図３８】検索領域Ｒ３から検出された大動脈Ａの位置を示す図である。
【図３９】スライスＳＬ４〜ＳＬ１から検出された大動脈Ａの位置を表す回転角θ＝θ４〜θ１を概略的に示す図である。
【図４０】スライスＳＬ１〜ＳＬ１０の大動脈Ａの中心を概略的に示す図である。
【図４１】トラッカー領域Ｒｔを概略的に示す図である。
【図４２】第２の形態のＭＲ装置２００を概略的に示す図である。
【図４３】第２の形態においてステップＳＴ３を実行するときのフローの説明図である。
【図４４】腕部が削除された２値画像ＤＣ１を概略的に示す図である。
【図４５】通常の閾値処理でスライスＳＬ１のアキシャル画像Ｄ１を２値化することにより得られた２値画像ＤＢ１を概略的に示す。
【図４６】腕部が削除された２値画像ＤＤ１を概略的に示す図である。
【図４７】ステップＳＴ３４０で実行される処理を説明するための図である。
【図４８】各スライスと面積比Ｒとの関係を表すプロファイルを示す図である。
【発明を実施するための形態】
【００１１】
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。
【００１２】
（１）第１の形態
図１は、本発明の第１の形態の磁気共鳴装置の概略図である。
磁気共鳴装置（以下、「ＭＲ装置」と呼ぶ。ＭＲ：Magnetic Resonance）１００は、マグネット２、テーブル３、受信コイル４などを有している。
【００１３】
マグネット２は、被検体１２が収容されるボア２１を有している。また、マグネット２は、超伝導コイル、勾配コイル、およびＲＦコイルなどが内蔵されている。
【００１４】
テーブル３は、クレードル３ａを有している。クレードル３ａは、ボア２１内に移動できるように構成されている。クレードル３ａによって、被検体１２はボア２１に搬送される。
【００１５】
受信コイル４は、被検体１２に取り付けられている。受信コイル４は、被検体１２からの磁気共鳴信号を受信する。
【００１６】
ＭＲ装置１００は、更に、造影剤注入装置５、送信器６、勾配磁場電源７、受信器８、制御部９、操作部１０、および表示部１１などを有している。
【００１７】
造影剤注入装置５は、被検体１２に造影剤を注入する。
送信器６はＲＦコイルに電流を供給し、勾配磁場電源７は勾配コイルに電流を供給する。
受信器８は、受信コイル４から受け取った信号に対して、検波などの信号処理を実行する。
【００１８】
制御部９は、表示部１１に必要な情報を伝送したり、受信器８から受け取ったデータに基づいて画像を再構成するなど、ＭＲ装置１００の各種の動作を実現するように、ＭＲ装置１００の各部の動作を制御する。制御部９は、画像作成手段９１〜設定手段９７などを有している。
【００１９】
画像作成手段９１は、被検体１２のアキシャル画像（例えば、図５参照）を作成する。
第１の２値化手段９２は領域拡張法を用いて画像を２値化する。
第２の２値化手段９３は閾値処理を用いて画像を２値化する。
算出手段９４は後述する面積比Ｒ（例えば、図２６参照）を算出する。
選択手段９５は大動脈を検出するときの基準となるスライスを選択する。
検出手段９６は大動脈を検出する。
設定手段９７は、造影剤を検出するためのトラッカー領域を設定する。
【００２０】
制御部９は、画像作成手段９１〜設定手段９７を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。制御部９は血管検出装置に相当する。
【００２１】
操作部１０は、オペレータにより操作され、種々の情報を制御部９に入力する。表示部１１は種々の情報を表示する。
ＭＲ装置１００は、上記のように構成されている。
【００２２】
図２は本形態で実行されるスキャンを示す図、図３は撮影部位を概略的に示す図である。
本形態では、ローカライザスキャンＬＳと本スキャンＭＳなどが実行される。
【００２３】
ローカライザスキャンＬＳは、スライス位置やトラッカー領域Ｒｔを設定するときに使用される画像を取得するためのスキャンである。トラッカー領域Ｒｔは、大動脈Ａに流れる造影剤を検出するために設定される領域である。
ローカライザスキャンＬＳの後に本スキャンＭＳが実行される。
【００２４】
本スキャンＭＳでは、被検体に造影剤が注入され、トラッカー領域Ｒｔから造影剤を検出するためのシーケンスが繰り返し実行される。そして、トラッカー領域Ｒｔに所定量の造影剤が注入したときに、肝臓の画像を取得するためのイメージングシーケンスが実行される。以下に、ローカライザスキャンＬＳおよび本スキャンＭＳを実行するときのフローの一例について説明する。
【００２５】
図４は、ローカライザスキャンＬＳおよび本スキャンＭＳを実行するときのフローの一例を示す図である。
ステップＳＴ１では、ローカライザスキャンＬＳ（図２参照）を実行する。
【００２６】
図５は、ローカライザスキャンＬＳの説明図である。
ローカライザスキャンＬＳは、被検体の肺および肝臓を含む部位の画像を取得するためのスキャンである。図５には、ローカライザスキャンＬＳを実行するときのスライスが示されている。ローカライザスキャンＬＳでは、アキシャル、サジタル、およびコロナルの画像が取得されるが、図５では、説明の便宜上、アキシャルの画像を取得するためのｍ枚のスライス（ここでは、ｍ＝１０枚のスライスＳＬ１〜ＳＬ１０）が示されている。スライスＳＬ１〜ＳＬ１０にはスライス番号ｎ（ｎ＝１〜１０）が対応付けられている。スライス番号ｎは、スライス位置がＳ側からＩ側に向かうに従って「１」づつ大きくなるように設定されている。画像作成手段９１（図１参照）は、ローカライザスキャンＬＳにより収集されたデータに基づいて、スライスＳＬ１〜ＳＬ１０の画像Ｄ１〜Ｄ１０を作成する（以下では、この画像を「アキシャル画像」と呼ぶ）。スライスＳＬ１〜ＳＬ１０は大動脈Ａを横切っているので、アキシャル画像Ｄ１〜Ｄ１０には、大動脈Ａの断面が描出される。アキシャル画像Ｄ１〜Ｄ１０の大動脈Ａの断面の位置を比較すると、大動脈Ａの断面は、Ｓ側からＩ側に近づくに従って、背骨を中心にして反時計方向ＡＮに回転していることが分かる。
【００２７】
ローカライザスキャンＬＳで使用されるシーケンスは、流速の速い血液のエコー信号ができるだけ収束しないように設計されている。したがって、流速の速い血液の信号は十分に低減することができる。例えば、大動脈Ａを流れる血液の流速について考えると、大動脈Ａを流れる血液の流速は速いので、大動脈Ａの血液の信号は十分に低減することができる。図５では、大動脈Ａを黒で表している。
アキシャル画像Ｄ１〜Ｄ１０を作成した後、ステップＳＴ２に進む。
【００２８】
ステップＳＴ２では、アキシャル画像Ｄ１〜Ｄ１０の各々から、脳脊髄液ＣＳＦを検出する。脳脊髄液ＣＳＦの検出方法は、例えば、“Med Imag Tech Vol.31 No.2 March 2013”に記載されている方法を用いることができる。アキシャル画像Ｄ１〜Ｄ１０ごとに脳脊髄液ＣＳＦを検出した後、ステップＳＴ３に進む。
【００２９】
ステップＳＴ３では、検出手段９６（図１参照）が、アキシャル画像Ｄ１〜Ｄ１０の各々から大動脈Ａを検出する。本形態では、大動脈Ａを識別するための識別器を用いて大動脈Ａを検出する。以下に、識別器の作成方法について説明する。
【００３０】
識別器は、被検体を撮影する前に予め準備されている。本形態では、識別器は、機械学習により作成されている。具体的には、教師データを用意し、機械学習で教師データを学習させることにより、大動脈Ａを検出するのに適した識別器Ｃを作成している（図６参照）。
【００３１】
図６は、識別器Ｃを作成するために使用された教師データの一例を概略的に示す図である。
教師データｖ１〜ｖｐは、実際の複数の人間を撮影することにより得られた腹部のアキシャル画像や胸部のアキシャル画像などを用いて作成することができる。
【００３２】
教師データｖ１〜ｖｐは、大動脈の断面が長方形領域の右側に内接するように規定されている。教師データｖ１〜ｖｐの大動脈は黒色で示されている。
【００３３】
これらの教師データｖ１〜ｖｐを機械学習で学習させることにより、アキシャル画像の中から大動脈を検出するのに適した識別器Ｃを作成する。
【００３４】
本形態では、検出手段９６が、上記のようにして作成された識別器Ｃを用いて、スライスＳＬ１〜ＳＬ１０の各々から大動脈Ａを検出する。以下に、本形態における大動脈Ａの検出方法について説明する。尚、本形態における大動脈Ａの検出方法の説明にあたっては、本形態における大動脈Ａの検出方法の効果を明確にするために、先ず、検出手段９６が、本形態とは別の方法で大動脈Ａを検出するやり方について説明する。
【００３５】
図７〜図１５は、本形態とは別の方法で大動脈Ａを検出するときの手順を示す図である。
先ず、肺側に位置するスライスＳＬ１のアキシャル画像Ｄ１に、大動脈Ａの検索範囲を規定する検索領域を設定する（図７参照）。
【００３６】
図７は、アキシャル画像Ｄ１に設定された検索領域Ｒ１を示す図である。
検出手段９６は、脳脊髄液ＣＳＦから水平に延びるラインＬ０と、脳脊髄液ＣＳＦを中心としてラインＬ０に対してθ＝θ０だけ傾いたラインＬ１とを規定する。ラインＬ０とラインＬ１とに挟まれる領域が、大動脈Ａの検索範囲を規定する検索領域Ｒ１として設定される。検索領域Ｒ１は、大動脈Ａが含まれるように設定される。θ０は、例えば、θ０＝４５°である。
【００３７】
次に、検出手段９６は、検索領域Ｒ１に、大動脈Ａを検出するときに使用するウィンドウを設定する（図８参照）。
【００３８】
図８は、検索領域Ｒ１に設定されたウィンドウＷを示す図である。
ウィンドウＷは、ｘ×ｙのピクセルサイズを有している。検出手段９６は、脳脊髄液ＣＳＦを基準にして、ラインＬ０上にウィンドウＷを設定する。ウィンドウＷを設定した後、検出手段９６は、検索領域Ｒ１内で、脳脊髄液ＣＳＦを中心としてウィンドウＷを回転させることによりウィンドウＷの回転角θを変更し、更にウィンドウＷのサイズ（ピクセル数ｘ、ピクセル数ｙ）も変更する。そして、識別器Ｃ（図６参照）を用いて、大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを求める（図９参照）。
【００３９】
図９は、識別器Ｃを用いて、大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘ、およびピクセル数ｙを求めるときの説明図である。
【００４０】
検出手段９６は、アキシャル画像Ｄ１に設定されたウィンドウＷの３つのパラメータ（θ，ｘ，ｙ）のうちのいずれかのパラメータを変更するたびに、ウィンドウＷ内のピクセルデータを抽出する。識別器Ｃは、抽出されたデータに基づいて、アキシャル画像Ｄ１の大動脈ＡがウィンドウＷの先端部に内接している可能性が高いか低いかを判断するための値を出力する。検出手段９６は、識別器Ｃの出力値Ｏｕｔに基づいて、アキシャル画像Ｄ１の大動脈ＡがウィンドウＷの先端部に内接している可能性が高いか低いかを判断する。本形態では、ウィンドウＷの先端部に大動脈Ａが内接している可能性が高いほど出力値Ｏｕｔが大きくなるように、識別器Ｃが構成されている。したがって、識別器Ｃの出力値Ｏｕｔが最大のときに、アキシャル画像Ｄ１の大動脈ＡがウィンドウＷの先端部に内接している可能性が最も高いと判断することができる。ここでは、回転角θ＝θ１、ピクセル数ｘ＝ｘ１、およびピクセル数ｙ＝ｙ１のときに、識別器Ｃの出力値Ｏｕｔが最大になったとする。したがって、ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）のときに、アキシャル画像Ｄ１の大動脈ＡがウィンドウＷの先端部に内接していると判断される。図１０に、アキシャル画像Ｄ１の大動脈ＡがウィンドウＷの先端部に内接しているときのウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）を示す。
【００４１】
ウィンドウＷの回転角θ１およびピクセル数ｙ１は、スライスＳＬ１における大動脈Ａの位置を表しており、ウィンドウＷのピクセル数ｘ１はスライスＳＬ１における大動脈Ａの内径を表している。したがって、ウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）を求めることによって、スライスＳＬ１の中から大動脈Ａを検出することができる。
【００４２】
スライスＳＬ１においてウィンドウＷのパラメータ（θ，ｘ，ｙ）＝（θ１，ｘ１，ｙ１）を求めた後、隣のスライスＳＬ２におけるアキシャル画像Ｄ２に、大動脈Ａの検索範囲を規定する検索領域を設定する（図１１参照）。
【００４３】
図１１は、アキシャル画像Ｄ２に設定された検索領域Ｒ２を示す図である。
スライスＳＬ２はスライスＳＬ１の隣に位置しているので、スライスＳＬ２における大動脈Ａの位置は、スライスＳＬ１における大動脈Ａの位置に近いと考えられる。そこで、検出手段９６は、スライスＳＬ１のアキシャル画像Ｄ１から検出した大動脈Ａの位置を表す回転角θ１を基準にして、回転角θ１を含む角度範囲θｓ２≦θ≦θｅ２を規定し、この角度範囲θｓ２≦θ≦θｅ２を検索領域Ｒ２として設定する。
【００４４】
次に、検出手段９６は、ウィンドウＷを用いて、検索領域Ｒ２の中から大動脈Ａを検出する。図１２に、検索領域Ｒ２から大動脈Ａが検出されたときの様子を概略的に示す。検出手段９６は、先に説明したように、ウィンドウＷの回転角θおよびウィンドウＷのサイズ（ピクセル数ｘ、ピクセル数ｙ）を変更し、アキシャル画像Ｄ３の大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘおよびｙを求める。図１２から、大動脈Ａは、θ＝θ２のライン上に位置していることがわかる。
【００４５】
スライスＳＬ２における大動脈Ａを検出した後、検出手段９６は、スライスＳＬ２に隣接するスライスＳＬ３のアキシャル画像Ｄ３に、大動脈Ａの検索範囲を規定する検索領域を設定する（図１３参照）。
【００４６】
図１３は、アキシャル画像Ｄ３に設定された検索領域Ｒ３を示す図である。
スライスＳＬ３はスライスＳＬ２の隣に位置しているので、スライスＳＬ３における大動脈Ａの位置は、スライスＳＬ２における大動脈Ａの位置に近いと考えられる。そこで、検出手段９６は、スライスＳＬ２のアキシャル画像Ｄ２から検出した大動脈Ａの位置を表す回転角θ２を基準にして、回転角θ２を含む角度範囲θｓ３≦θ≦θｅ３を規定し、この角度範囲θｓ３≦θ≦θｅ３を、検索領域Ｒ３として設定する。
【００４７】
次に、検出手段９６は、ウィンドウＷを用いて、検索領域Ｒ３の中から大動脈Ａを検出する。図１４に、検索領域Ｒ３から検出された大動脈Ａの位置を示す。大動脈Ａは、θ＝θ３のライン上に位置している。
【００４８】
スライスＳＬ３の大動脈Ａを検出した後も、同様の手法で、隣のスライスのアキシャル画像に検索領域を設定し、ウィンドウＷを用いて検索領域の中から大動脈Ａを検出する。したがって、全てのスライスＳＬ１〜ＳＬ１０における大動脈Ａを検出することができる。図１５に、スライスＳＬ１〜ＳＬ１０から検出された大動脈Ａの位置を表す回転角θ＝θ１〜θ１０を概略的に示す。
【００４９】
上記のように、図７〜図１５の方法では、スライスＳＬ１からスライスＳＬ１０に向かう順（スライス番号ｎの小さい順）に大動脈Ａが検出される。しかし、スライスＳＬ１は、肺を横切る位置に設定されているので、スライスＳＬ１のアキシャル画像Ｄ１には、大動脈Ａの他に、肺の領域も低信号で描出される。したがって、スライスＳＬ１のアキシャル画像Ｄ１から大動脈Ａを検出する場合、大動脈Ａと同様に低信号で描出された肺の領域を、大動脈として誤検出することがある。図１６に、誤検出の例を示す。図１６では、回転角θ＝θ１´において、肺の一部の領域を大動脈として誤検出している。したがって、隣のスライスＳＬ２に検索領域Ｒ２を設定する場合、スライスＳＬ１において大動脈として誤検出された領域の回転角θ＝θ１´に基づいて、スライスＳＬ２の検索領域Ｒ２が設定される。図１７に、スライスＳＬ２に設定された検索領域Ｒ２を示す。図１７では、スライスＳＬ２における大動脈Ａが検索領域Ｒ２に含まれていないので、スライスＳＬ２の大動脈Ａを検出することができない。したがって、残りのスライスＳＬ３〜ＳＬ１０でも大動脈を検出することができないという問題がある。この問題に対処する方法として、肺から最も離れた位置のスライスＳＬ１０から、肺を横切るスライスＳＬ１に向かう順（スライス番号ｎの大きい順）に大動脈Ａを検出することが考えられる。しかし、スライスＳＬ１０のように、肝臓の下端に近いスライスには、大動脈Ａとは別の動脈の断面も現れやすいので、スライスＳＬ１０から大動脈Ａの検出を開始すると、別の動脈を大動脈Ａとして検出する可能性がある。したがって、やはり、大動脈Ａを正しく検出することができないという問題がある。そこで、本形態では、各スライスの大動脈を正しく検出することができるように、ステップＳＴ３の処理を工夫している。以下に、本形態におけるステップＳＴ３の処理について説明する。
【００５０】
図１８は、本形態におけるステップＳＴ３のフローを示す図である。
ステップＳＴ３０では、スライスが肝臓を横切るか否かを判断するための指標を算出する。ステップＳＴ３０は、ステップＳＴ３１０〜ＳＴ３４０を有している。以下、各ステップについて説明する。
【００５１】
ステップＳＴ３１０では、スライスＳＬ１〜ＳＬ１０のアキシャル画像Ｄ１〜Ｄ１０の中から、１枚のアキシャル画像を選択する。ここでは、肺を横切るスライスＳＬ１のアキシャル画像Ｄ１が選択されたとする。図１９に、選択されたスライスＳＬ１のアキシャル画像Ｄ１を概略的に示す。アキシャル画像Ｄ１には、被検体の胴部だけでなく、左腕部および右腕部も描出されている。スライスＳＬ１のアキシャル画像Ｄ１を選択した後、ステップＳＴ３２０に進む。
【００５２】
ステップＳＴ３２０では、第１の２値化手段９２が、領域拡張法を用いて、スライスＳＬ１の２値化を実行する。
【００５３】
第１の２値化手段９２は、先ず、アキシャル画像Ｄ１の中の背景領域（被検体の体外領域）にシード領域を設定する。図２０にシード領域Ｒｓｅｅｄの一例を示す。シード領域Ｒｓｅｅｄを背景領域に設定する方法としては、クレードル３ａの位置情報を用いる方法がある。クレードル３ａは被検体の体外に位置しているので、クレードル３ａの位置にシード領域Ｒｓｅｅｄを設定することにより、シード領域Ｒｓｅｅｄを確実に背景領域（体外領域）に設定することができる。尚、オペレータがシード領域Ｒｓｅｅｄを手動で設定してもよい。
【００５４】
シード領域Ｒｓｅｅｄを設定した後、第１の２値化手段９２は、アキシャル画像Ｄ１の信号値に基づいて、被検体の体外と体内とを区別するときに使用する信号値の閾値ｖを求める。
【００５５】
被検体の体外の空間には、クレードル３ａや、ボア内の空気などが存在しているが、クレードル３ａの信号値や空気の信号値は、被検体の体内組織の信号値よりも十分に小さい値になる。このように、体外と体内組織との間には信号値に大きな差があるので、スライスＳＬ１のアキシャル画像Ｄ１の信号値に基づいて、被検体の体外と体内とを区別するための信号値の閾値ｖを求めることが可能となる。閾値ｖを求める方法の一例としては、シード領域Ｒｓｅｅｄ内の信号値を用いる方法がある。シード領域Ｒｓｅｅｄは被検体の体外に設定されているので、シード領域Ｒｓｅｅｄの信号値は、被検体の体内の信号値よりも十分に小さい値（ノイズ）を有している。したがって、シード領域Ｒｓｅｅｄの信号値を用いることにより、被検体の体外と体内とを区別するための信号値の閾値ｖを求めることができる。
【００５６】
次に、第１の２値化手段９２は、ピクセルの信号値と閾値ｖとを比較しながら、閾値ｖより小さい信号値を有するピクセルの領域を拡張し、スライスＳＬ１のアキシャル画像Ｄ１を２値化する。図２１に、領域拡張法でスライスＳＬ１の画像Ｄ１を２値化することにより得られた２値画像ＤＡ１を概略的に示す。被検体の体外における大部分のピクセルの信号値は閾値ｖよりも小さいが、被検体の体表面における大部分のピクセルの信号値は閾値ｖよりも大きい。したがって、ピクセルの領域は、被検体の体外領域の全体に渡って拡張されるが、被検体の体内領域の内側には拡張されないので、背景領域（被検体の体外）と被検体の体内とを区別することができる。図２１では、拡張された領域内のピクセルが論理値０（黒）で表され、他の領域内のピクセルが論理値１（白）で表されている。
【００５７】
領域拡張法を用いて２値化を実行した場合、背景領域（被検体の体外）には論理値０が割り当てられ、被検体の体内には論理値１が割り当てられる。尚、被検体の体内には論理値１が割り当てられるので、肺にも論理値１が割り当てられる。
【００５８】
領域拡張法を用いてスライスＳＬ１の２値画像ＤＡ１を得た後、ステップＳＴ３３０に進む。
【００５９】
ステップＳＴ３３０では、第２の２値化手段９３が、領域拡張法を用いずに、通常の閾値処理を用いて、スライスＳＬ１のアキシャル画像Ｄ１の２値化を実行する。具体的には、スライスＳＬ１のアキシャル画像Ｄ１の中で、ピクセルの信号値が閾値ＴＨより小さい領域に論理値０を割り当て、一方、ピクセルの信号値が閾値ＴＨより大きい領域に論理値１を割り当てることにより、スライスＳＬ１のアキシャル画像Ｄ１を２値化する。図２２に、通常の閾値処理でスライスＳＬ１のアキシャル画像Ｄ１を２値化することにより得られた２値画像ＤＢ１を概略的に示す。閾値ＴＨは、例えば、ステップＳＴ３２０の領域拡張法を用いた２値化を実行するときに求めた閾値ｖを使用することができる。
【００６０】
先に説明したように、体内組織と空気との間には信号値に大きな差があり、体内組織の信号値は大きくなるが、空気の信号値は小さくなる。したがって、体内組織に位置するピクセルの大部分は閾値ＴＨよりも大きい信号値を有するが、空気で占められる空間に位置するピクセルの大部分は閾値ＴＨよりも小さい信号値を有する。このような信号値の差異が生じているので、通常の閾値処理により得られた２値画像ＤＢ１は、被検体の体内組織の大部分に論理値１が割り当てられるが、背景領域（被検体の体外）や肺の内側（空気が含まれている空間）の大部分に論理値０が割り当てられる。２値化処理を実行した後、ステップＳＴ３４０に進む。
【００６１】
図２３は、ステップＳＴ３４０で実行される処理を説明するための図である。
ステップＳＴ３４０では、算出手段９４（図１参照）が、２値画像ＤＡ１の中の論理値１が割り当てられている領域の面積Ａと、２値画像ＤＢ１の中の論理値１が割り当てられている領域の面積Ｂとを計算する。面積ＡおよびＢの値を求めたら、算出手段９４は、面積比Ｒ＝Ｂ／Ａを計算する。この面積比Ｒが、スライスＳＬ１が肝臓を横切るか否かを判断するための指標として用いられる。面積比Ｒを用いてスライスが肝臓を横切るか否かを判断する具体的なやり方については後述する。
【００６２】
図２３に示すように、領域拡張法を用いた２値化を実行した場合、肺に論理値１を割り当てることができる。一方、閾値処理を用いた２値化を実行した場合、被検体の体内組織は論理値１が割り当てられるが、肺の内側は論理値０が割り当てられる。したがって、２つの２値画像ＤＡ１およびＤＢ１を比較すると、２値画像ＤＢ１の論理値１が割り当てられている領域の面積Ｂは、２値画像ＤＡ１の論理値１が割り当てられている領域の面積Ａよりも小さい値となる。したがって、スライスＳＬ１における面積比Ｒは、１よりもかなり小さい値になる。
【００６３】
このようにして、ステップＳＴ３０が実行される。尚、図１９〜図２３では、スライスＳＬ１のアキシャル画像Ｄ１に対してステップＳＴ３０を実行した例について説明されているが、他のスライスＳＬ２〜ＳＬ１０のアキシャル画像に対しても、スライスＳＬ１のアキシャル画像Ｄ１と同様に、ステップＳＴ３０を実行する。図２４および図２５に、他のスライスの画像に対してステップＳＴ３０を実行する例を示す。以下、図２４および図２５について順に説明する。
【００６４】
図２４は、スライスＳＬ３の画像に対してステップＳＴ３０を実行するときの説明図である。
先ず、ステップＳＴ３１０において、スライスＳＬ３のアキシャル画像Ｄ３を選択する。スライスＳＬ３は、肺と肝臓との境界付近に位置しており、肺と肝臓との両方の臓器を横切っている。スライスＳＬ３のアキシャル画像Ｄ３を選択した後、ステップＳＴ３２０に進む。
【００６５】
ステップＳＴ３２０では、第１の２値化手段９２が、アキシャル画像Ｄ３の中の背景領域（被検体の体外領域）にシード領域を設定し、領域拡張法を用いて２値化を実行する。したがって、スライスＳＬ３の２値画像ＤＡ３が得られる。領域拡張法を用いて２値化を実行した場合、背景領域（被検体の体外）には論理値０が割り当てられ、被検体の体内には論理値１が割り当てられる。被検体の体内は論理値１で表されるので、肺にも肝臓にも論理値１が割り当てられる。
領域拡張法を用いてスライスＳＬ３の２値画像ＤＡ３を得た後、ステップＳＴ３３０に進む。
【００６６】
ステップＳＴ３３０では、第２の２値化手段９３が、領域拡張法を用いずに、通常の閾値処理を用いて、スライスＳＬ３のアキシャル画像Ｄ３の２値化を実行する。これにより、スライスＳＬ３の２値画像ＤＢ３が得られる。通常の閾値処理により得られたスライスＳＬ３の２値画像ＤＢ３は、背景領域（被検体の体外）や肺の内側（空気が含まれている空間）の大部分に論理値０が割り当てられるが、肝臓の大部分には論理値１が割り当てられる。したがって、スライスＳＬ３の２値画像ＤＢ３は、スライスＳＬ１の２値画像ＤＢ１（図２３参照）と比較すると、論理値１が割り当てられている領域の面積が広くなる。２値化処理を実行した後、ステップＳＴ３４０に進む。
【００６７】
ステップＳＴ３４０では、算出手段９４が、２値画像ＤＡ３の中の論理値１が割り当てられている領域の面積Ａと、２値画像ＤＢ３の中の論理値１が割り当てられている領域の面積Ｂとを計算する。面積ＡおよびＢの値を求めたら、算出手段９４は、面積比Ｒ＝Ｂ／Ａを計算する。
【００６８】
２値画像ＤＡ３およびＤＢ３を比較すると、２値画像ＤＢ３の論理値１が割り当てられている領域の面積Ｂは、２値画像ＤＡ３の論理値１が割り当てられている領域の面積Ａよりも小さい値となる。したがって、スライスＳＬ３における面積比Ｒは、１よりも小さい値になる。
【００６９】
このようにして、スライスＳＬ３のアキシャル画像に対してステップＳＴ３が実行される。次に、スライスＳＬ７のアキシャル画像に対してステップＳＴ３を実行する場合について説明する。
【００７０】
図２５は、スライスＳＬ７の画像に対してステップＳＴ３を実行するときの説明図である。
先ず、ステップＳＴ３１０において、スライスＳＬ７のアキシャル画像Ｄ７を選択する。スライスＳＬ７は、肺を横切っておらず、肝臓を横切っている。スライスＳＬ７のアキシャル画像Ｄ７を選択した後、ステップＳＴ３２０に進む。
【００７１】
ステップＳＴ３２０では、第１の２値化手段９２が、アキシャル画像Ｄ７の中の背景領域（被検体の体外領域）にシード領域を設定し、領域拡張法を用いて２値化を実行する。したがって、スライスＳＬ７の２値画像ＤＡ７が得られる。領域拡張法を用いて２値化を実行した場合、背景領域（被検体の体外）には論理値０が割り当てられ、被検体の体内には論理値１が割り当てられる。被検体の体内は論理値１で表されるので、肝臓には論理値１が割り当てられる。
領域拡張法を用いてスライスＳＬ７の２値画像ＤＡ７を得た後、ステップＳＴ３３０に進む。
【００７２】
ステップＳＴ３３０では、第２の２値化手段９３が、領域拡張法を用いずに、通常の閾値処理を用いて、スライスＳＬ７のアキシャル画像Ｄ７の２値化を実行する。これにより、スライスＳＬ７の２値画像ＤＢ７が得られる。スライスＳＬ７は肺を横切っておらず、肝臓を横切っているので、２値画像ＤＢ７は、被検体の体内の大部分の領域に、論理値１が割り当てられている。２値化処理を実行した後、ステップＳＴ３４０に進む。
【００７３】
ステップＳＴ３４０では、算出手段９４が、２値画像ＤＡ７の中の論理値１が割り当てられている領域の面積Ａと、２値画像ＤＢ７の中の論理値１が割り当てられている領域の面積Ｂとを計算する。面積ＡおよびＢの値を求めたら、算出手段９４は、面積比Ｒ＝Ｂ／Ａを計算する。
【００７４】
２値画像ＤＡ７およびＤＢ７を比較すると、２値画像ＤＢ７の論理値１が割り当てられている領域の面積Ｂは、２値画像ＤＡ７の論理値１が割り当てられている領域の面積Ａとほぼ同じである。したがって、スライスＳＬ７における面積比Ｒは、１に近い値になる。
【００７５】
このように、算出手段９４は、各スライスの画像に対して面積比Ｒを算出する。したがって、各スライスと面積比Ｒとの関係を表すプロファイルを得ることができる。図２６に、各スライスと面積比Ｒとの関係を表すプロファイルを示す。
【００７６】
図２６のプロファイルを参照すると、スライスが肺側に位置しているか、肝臓側に位置しているかに応じて、面積比Ｒが異なることがわかる。以下に、スライスＳＬ１〜ＳＬ１０のうち、代表して、３つのスライスＳＬ１、ＳＬ３、およびＳＬ７を取り上げ、３つのスライスＳＬ１、ＳＬ３、およびＳＬ７の面積比Ｒの違いについて説明する。尚、説明の便宜上、肺を横切るスライスＳＬ１の面積比Ｒと、肝臓を横切るスライスＳＬ７の面積比Ｒを先に説明し、最後に、肝臓と肺との境界を横切るスライスＳＬ３の面積比Ｒについて説明する。
【００７７】
（１）肺を横切るスライスＳＬ１の面積比Ｒについて
スライスＳＬ１では、２値画像ＤＡ１の肺の領域には論理値１が割り当てられるが、２値画像ＤＢ１の肺の領域には論理値０が割り当てられる。したがって、２値画像ＤＢ１の論理値１の面積Ｂは、２値画像ＤＡ１の論理値１の面積Ａよりもかなり小さいので、スライスＳＬ１における面積比Ｒは１よりもかなり小さい値になる。
【００７８】
（２）肝臓を横切るスライスＳＬ７の面積比Ｒについて
スライスＳＬ７では、２値画像ＤＡ７の肝臓の領域には論理値１が割り当てられ、２値画像ＤＢ７の肝臓の領域にも論理値１が割り当てられる。したがって、２値画像ＤＢ７の論理値１の面積Ｂは、２値画像ＤＡ７の論理値１の面積Ａにほぼ等しいので、スライスＳＬ７における面積比Ｒは１に近い値になる。
【００７９】
（３）肝臓と肺との境界を横切るスライスＳＬ３の面積比Ｒについて
スライスＳＬ３の２値画像ＤＡ３では、肝臓の領域に論理値１が割り当てられ、肺の領域にも論理値１が割り当てられる。一方、スライスＳＬ３の２値画像ＤＢ３では、肝臓の領域には論理値１が割り当てられるが、肺の領域には論理値０が割り当てられる。したがって、スライスＳＬ３における面積比Ｒは、スライスＳＬ１における面積比Ｒよりも大きいが、スライスＳＬ７における面積比Ｒよりも小さくなる。
【００８０】
したがって、面積比Ｒが１に近いほど、スライスが肝臓を横切る可能性は高く、面積比Ｒが小さくなるに従って、スライスが肝臓を横切る可能性は低くなることがわかる。
スライスＳＬ１〜ＳＬ１０の面積比Ｒを求めた後、ステップＳＴ３１に進む。
【００８１】
ステップＳＴ３１では、選択手段９５が、面積比Ｒに基づいて、スライスＳＬ１〜ＳＬ１０の中から、大動脈Ａを検出するときの基準となるスライスＳＬｊ（スライス番号ｎ＝ｊ）を選択する。以下に、スライスＳＬｊを選択する方法について説明する。
【００８２】
選択手段９５は、先ず、スライスＳＬ１〜ＳＬ１０の中から、肝臓側に位置するスライスを求める（図２７参照）。
【００８３】
図２７は、スライスＳＬ１〜ＳＬ１０の中から、肝臓側に位置するスライスを求めるときの説明図である。
【００８４】
選択手段９５は、面積比Ｒの最大値ｍａｘと面積比Ｒの最小値ｍｉｎとの差ΔＲを求める。そして、ΔＲに基づいて、面積比Ｒが１に近いか否かを判断するための閾値ＴＨを設定する。閾値ＴＨは、以下の式で表すことができる。
ＴＨ＝ｍｉｎ＋ｋ・ΔＲ    ・・・（１）
ここで、ｍｉｎ：面積比Ｒの最小値
ｋ：係数
【００８５】
ｋは、例えば、ｋ＝０．９である。本形態では、閾値ＴＨ以上の面積比Ｒを有するスライスを、肝臓側に位置するスライスとして求める。図２７では、スライスＳＬ５〜ＳＬ１０が、肝臓側に位置するスライスとして求められる。先に説明したように、スライスの面積比Ｒが１に近いほど、スライスが肝臓側に位置している可能性は高いと考えられるので、閾値ＴＨ以上の面積比Ｒを有するスライスＳＬ５〜ＳＬ１０は、肝臓を横切っていると考えることができる。
【００８６】
スライスＳＬ５〜ＳＬ１０を求めた後、選択手段９５は、スライスＳＬ５〜ＳＬ１０の中から、一つのスライスを、大動脈を検出するときの基準となるスライスＳＬｊとして選択する。ここでは、ｊ＝５、すなわち、スライスＳＬ５〜ＳＬ１０のうち、一番肺側に位置するスライスＳＬ５（スライス番号ｎ＝５）が、大動脈Ａを検出するときの基準となるスライスＳＬｊとして選択される。スライスＳＬ５は面積比Ｒが１に近い値であるので、スライスＳＬ５は肺をほとんど横切っていないと考えることができる。したがって、スライスＳＬ５における大動脈Ａの周囲には、肺が原因で生じる低信号領域は存在していないと考えられる。スライスＳＬ５を選択した後、ステップＳＴ３２に進む。
【００８７】
ステップＳＴ３２では、検出手段９６が、ステップＳＴ３１で選択されたスライスＳＬ５の大動脈Ａを検出する。
【００８８】
図２８および図２９は、スライスＳＬ５の大動脈Ａを検出するときの説明図である。
検出手段９６は、先ず、スライスＳＬ５のアキシャル画像Ｄ５に、検索範囲を規定する検索領域を設定する。図２８に、アキシャル画像Ｄ５に設定された検索領域Ｒ５を示す。
【００８９】
検出手段９６は、脳脊髄液ＣＳＦから水平に延びるラインＬ０と、脳脊髄液ＣＳＦを中心としてラインＬ０に対してθ＝θ０だけ傾いたラインＬ１とを規定する。ラインＬ０とラインＬ１とに挟まれる領域が、大動脈Ａの検索範囲を規定する検索領域Ｒ５として設定される。検索領域Ｒ５は、大動脈Ａが含まれるように設定される。θ０は、例えば、θ０＝９０°である。
【００９０】
次に、検出手段９６は、ウィンドウＷを用いて、検索領域Ｒ５の中から大動脈Ａを検出する。図２９に、検索領域Ｒ５から大動脈Ａが検出されたときの様子を概略的に示す。検出手段９６は、先に説明したように、ウィンドウＷの回転角θおよびウィンドウＷのサイズ（ピクセル数ｘ、ピクセル数ｙ）を変更し、アキシャル画像Ｄ５の大動脈ＡがウィンドウＷの先端部に内接するときのウィンドウＷの回転角θ、ピクセル数ｘおよびｙを求める。図２９から、大動脈Ａは、θ＝θ５のライン上に位置していることがわかる。
スライスＳＬ５の大動脈Ａを検出した後、ステップＳＴ３３に進む。
【００９１】
ステップＳＴ３３では、検出手段９６が、スライスＳＬ５に対してスライス番号ｎが大きくなる方向に並ぶスライスＳＬ６〜ＳＬ１０の大動脈Ａを検出する。
【００９２】
図３０〜図３４は、スライスＳＬ６〜ＳＬ１０の大動脈Ａを検出するときの説明図である。
ステップＳＴ３３では、検出手段９６は、先ず、スライスＳＬ６〜ＳＬ１０のうち、スライスＳＬ５に隣接するスライスＳＬ６の大動脈Ａを検出する（図３０および図３１参照）。
【００９３】
図３０および図３１は、スライスＳＬ６の大動脈Ａを検出するときの説明図である。
検出手段９６は、先ず、スライスＳＬ６のアキシャル画像Ｄ６に、検索範囲を規定する検索領域を設定する。図３０に、アキシャル画像Ｄ６に設定された検索領域Ｒ６を示す。
【００９４】
スライスＳＬ６はスライスＳＬ５の隣に位置しているので、スライスＳＬ６における大動脈Ａの位置は、スライスＳＬ５における大動脈Ａの位置に近いと考えられる。そこで、検出手段９６は、スライスＳＬ５のアキシャル画像Ｄ５から検出した大動脈Ａの位置を表す回転角θ５を基準にして、回転角θ５を含む角度範囲θｓ６≦θ≦θｅ６を規定し、この角度範囲θｓ６≦θ≦θｅ６を検索領域Ｒ６として設定する。
【００９５】
尚、先に説明したように、大動脈Ａの断面は、Ｓ側からＩ側に近づくに従って、背骨を中心にして反時計方向ＡＮに回転している。したがって、アキシャル画像Ｄ６における大動脈Ａは、θ＝θ５に対して反時計方向ＡＮに位置している可能性が高いので、検索領域Ｒ６のうちのθｓ６≦θ≦θ５の角度範囲は狭くし、θ５≦θ≦θｅ６の角度範囲を広く設定してもよい。このように、大動脈Ａの走行方向に基づいて検索領域Ｒ６を設定することにより、大動脈Ａが存在している可能性が低い領域を、検索領域Ｒ６から除外することができる。
【００９６】
検索領域Ｒ６を設定した後、検出手段９６は、ウィンドウＷを用いて、検索領域Ｒ６の中から大動脈Ａを検出する。図３１に、検索領域Ｒ６から検出された大動脈Ａの位置を示す。大動脈Ａは、θ＝θ６のライン上に位置している。
【００９７】
スライスＳＬ６における大動脈Ａを検出した後、検出手段９６は、スライスＳＬ６に対してＩ方向（スライス番号ｎが大きくなる方向）に隣接するスライスＳＬ７のアキシャル画像Ｄ７に、大動脈Ａの検索範囲を規定する検索領域を設定する（図３２参照）。
【００９８】
図３２は、アキシャル画像Ｄ７に設定された検索領域Ｒ７を示す図である。
【００９９】
スライスＳＬ７はスライスＳＬ６の隣に位置しているので、スライスＳＬ７における大動脈Ａの位置は、スライスＳＬ６における大動脈Ａの位置に近いと考えられる。そこで、検出手段９６は、スライスＳＬ６のアキシャル画像Ｄ６から検出した大動脈Ａの位置を表す回転角θ６を基準にして、回転角θ６を含む角度範囲θｓ７≦θ≦θｅ７を規定し、この角度範囲θｓ７≦θ≦θｅ７を検索領域Ｒ７として設定する。尚、大動脈Ａの走行方向を考慮して、検索領域Ｒ７のうちのθｓ７≦θ≦θ６の角度範囲は狭くし、θ６≦θ≦θｅ７の角度範囲を広く設定してもよい。
【０１００】
検索領域Ｒ７を設定した後、検出手段９６は、ウィンドウＷを用いて、検索領域Ｒ７の中から大動脈Ａを検出する。図３３に、検索領域Ｒ７から検出された大動脈Ａの位置を示す。大動脈Ａは、θ＝θ７のライン上に位置している。
【０１０１】
スライスＳＬ７の大動脈Ａを検出した後も、同様の手法で、スライスＳＬ８、ＳＬ９、およびＳＬ１０の大動脈Ａを検出する。したがって、スライスＳＬ５に対してスライス番号ｎが大きい方向に並ぶスライスＳＬ６〜ＳＬ１０の大動脈Ａを検出することができる。図３４に、スライスＳＬ６〜ＳＬ１０から検出された大動脈Ａの位置を表す回転角θ＝θ６〜θ１０を概略的に示す。
【０１０２】
スライスＳＬ８の大動脈Ａを検出する場合、検出手段９６は、スライスＳＬ７から検出された大動脈の位置を表す回転角θ＝θ７に基づいて、スライスＳＬ８のアキシャル画像Ｄ８に検索領域Ｒ８を設定し、検索領域Ｒ８の中から大動脈Ａを検出する。図３４では、スライスＳＬ８の大動脈Ａは、θ＝θ８のライン上に位置している。スライスＳＬ８の大動脈Ａを検出した後、θ＝θ８に基づいて、隣のスライスＳＬ９のアキシャル画像Ｄ９に検索領域Ｒ９を設定し、検索領域Ｒ９の中から大動脈Ａを検出する。図３４では、スライスＳＬ９の大動脈Ａは、θ＝θ９のライン上に位置している。スライスＳＬ９の大動脈Ａを検出した後、θ＝θ９に基づいて、隣のスライスＳＬ１０のアキシャル画像Ｄ１０に検索領域Ｒ１０を設定し、検索領域Ｒ１０の中から大動脈Ａを検出する。このようにして、スライスＳＬ６〜ＳＬ１０の大動脈Ａを検出することができる。
スライスＳＬ６〜ＳＬ１０の大動脈Ａを検出した後、ステップＳＴ３４に進む。
【０１０３】
ステップＳＴ３４では、検出手段９６が、スライスＳＬ５に対してスライス番号ｎが小さくなる方向に並ぶスライスＳＬ４〜ＳＬ１の大動脈Ａを検出する。
【０１０４】
図３５〜図３９は、スライスＳＬ４〜ＳＬ１の大動脈Ａを検出するときの説明図である。
ステップＳＴ３４では、検出手段９６は、先ず、ＳＬ４〜ＳＬ１のうち、スライスＳＬ５に隣接するスライスＳＬ４の大動脈Ａを検出する（図３５および図３６参照）。
【０１０５】
図３５および図３６は、スライスＳＬ４の大動脈Ａを検出するときの説明図である。
【０１０６】
検出手段９６は、先ず、スライスＳＬ４のアキシャル画像Ｄ４に、検索範囲を規定する検索領域を設定する。図３５に、アキシャル画像Ｄ４に設定された検索領域Ｒ４を示す。
【０１０７】
スライスＳＬ４はスライスＳＬ５の隣に位置しているので、スライスＳＬ４における大動脈Ａの位置は、スライスＳＬ５における大動脈Ａの位置に近いと考えられる。そこで、検出手段９６は、スライスＳＬ５のアキシャル画像Ｄ５から検出した大動脈Ａの位置を表す回転角θ５を基準にして、回転角θ５を含む角度範囲θｓ４≧θ≧θｅ４を規定し、この角度範囲θｓ４≧θ≧θｅ４を検索領域Ｒ４として設定する。
【０１０８】
尚、図３５から、大動脈Ａの断面は、Ｉ側からＳ側に近づくに従って、背骨を中心にして時計方向ＣＲに回転している。したがって、アキシャル画像Ｄ４における大動脈Ａは、θ＝θ５に対して時計方向ＣＲに位置している可能性が高いので、検索領域Ｒ４のうちのθｓ４≧θ≧θ５の角度範囲は狭くし、θ５≧θ≧θｅ４の角度範囲を広く設定してもよい。このように、大動脈Ａの走行方向に基づいて検索領域Ｒ４を設定することにより、大動脈Ａが存在している可能性が低い領域を、検索領域Ｒ４から除外することができる。
【０１０９】
検索領域Ｒ４を設定した後、検出手段９６は、ウィンドウＷを用いて、検索領域Ｒ４の中から大動脈Ａを検出する。図３６に、検索領域Ｒ４から検出された大動脈Ａの位置を示す。大動脈Ａは、θ＝θ４のライン上に位置している。
【０１１０】
スライスＳＬ４における大動脈Ａを検出した後、検出手段９６は、スライスＳＬ４に対してＳ方向（スライス番号ｎが小さくなる方向）に隣接するスライスＳＬ３のアキシャル画像Ｄ３に、大動脈Ａの検索範囲を規定する検索領域を設定する（図３７参照）。
【０１１１】
図３７は、アキシャル画像Ｄ３に設定された検索領域Ｒ３を示す図である。
スライスＳＬ３はスライスＳＬ４の隣に位置しているので、スライスＳＬ３における大動脈Ａの位置は、スライスＳＬ４における大動脈Ａの位置に近いと考えられる。そこで、検出手段９６は、スライスＳＬ４のアキシャル画像Ｄ４から検出した大動脈Ａの位置を表す回転角θ４を基準にして、回転角θ４を含む角度範囲θｓ３≧θ≧θｅ３を規定し、この角度範囲θｓ３≧θ≧θｅ３を検索領域Ｒ３として設定する。尚、大動脈Ａの走行方向を考慮して、検索領域Ｒ３のうちのθｓ３≧θ≧θ４の角度範囲は狭くし、θ４≧θ≧θｅ３の角度範囲を広く設定してもよい。
【０１１２】
検索領域Ｒ３を設定した後、検出手段９６は、ウィンドウＷを用いて、検索領域Ｒ７の中から大動脈Ａを検出する。図３８に、検索領域Ｒ３から検出された大動脈Ａの位置を示す。大動脈Ａは、θ＝θ３のライン上に位置している。
【０１１３】
スライスＳＬ３の大動脈Ａを検出した後も、同様の手法で、スライスＳＬ２およびＳＬ１の大動脈Ａを検出する。したがって、スライスＳＬ５に対してスライス番号ｎが小さい方向に並ぶスライスＳＬ４〜ＳＬ１の大動脈Ａを検出することができる。図３９に、スライスＳＬ４〜ＳＬ１から検出された大動脈Ａの位置を表す回転角θ＝θ４〜θ１を概略的に示す。
【０１１４】
スライスＳＬ２の大動脈Ａを検出する場合、検出手段９６は、スライスＳＬ３から検出された大動脈の位置を表す回転角θ＝θ３に基づいて、スライスＳＬ２のアキシャル画像Ｄ２に検索領域Ｒ２を設定し、検索領域Ｒ２の中から大動脈Ａを検出する。図３９では、スライスＳＬ２の大動脈Ａは、θ＝θ２のライン上に位置している。スライスＳＬ２の大動脈Ａを検出した後、θ＝θ２に基づいて、隣のスライスＳＬ１のアキシャル画像Ｄ１に検索領域Ｒ１を設定し、検索領域Ｒ１の中から大動脈Ａを検出する。このようにして、スライスＳＬ４〜ＳＬ１の大動脈Ａを検出することができる。
【０１１５】
したがって、ステップＳＴ３２〜ＳＴ３４を実行することにより、全てのスライスＳＬ１〜ＳＬ１０の大動脈Ａを検出することができる。全てのスライスＳＬ１〜ＳＬ１０の大動脈Ａを検出した後、ステップＳＴ４（図４参照）に進む。
【０１１６】
ステップＳＴ４では、検出手段９６が、スライスＳＬ１〜ＳＬ１０から検出された大動脈Ａの中心を求める。図４０に、スライスＳＬ１〜ＳＬ１０の大動脈Ａの中心を概略的に示す。図４０には、紙面の制約上、スライスＳＬ１、ＳＬ２、ＳＬ３、およびＳＬ１０の大動脈Ａの中心のみを示してある。大動脈Ａの断面形状は略円形と見なすことができるので、ウィンドウＷに内接する円を考え、その円の中心を検出することにより、大動脈Ａの中心を求めることができる。大動脈Ａの中心を求めた後、ステップＳＴ５に進む。
【０１１７】
ステップＳＴ５では、設定手段９７（図１参照）が、検出された血管の中心に基づいて、造影剤を検出するためのトラッカー領域を設定する。図４１に、設定されたトラッカー領域Ｒｔを概略的に示す。トラッカー領域Ｒｔは、大動脈Ａの内側に位置するように設定される。トラッカー領域Ｒｔを設定した後、ステップＳＴ６に進む。
【０１１８】
ステップＳＴ６では、本スキャンＭＳ（図２参照）が実行される。本スキャンＭＳでは、被検体に造影剤が注入され、トラッカー領域Ｒｔから造影剤を検出するためのシーケンスが繰り返し実行される。そして、トラッカー領域Ｒｔに所定量の造影剤が注入したときに、肝臓の画像を取得するためのスキャンが実行され、フローが終了する。
【０１１９】
本形態では、スライスＳＬ１〜ＳＬ１０と面積比Ｒとの関係を表すプロファイルを求め（図２７参照）、このプロファイルに基づいて、大動脈Ａを検出するときの基準となるスライスＳＬ５を選択する。そして、スライスＳＬ１〜ＳＬ１０のうち、選択されたスライスＳＬ５の大動脈を最初に検出する。スライスＳＬ５の面積比Ｒは１に近い値であるので、スライスＳＬ５における大動脈Ａの周囲には、肺が原因で生じる低信号領域は存在していないと考えられる。したがって、スライスＳＬ５の大動脈を検出する場合、肺が原因で生じる低信号領域を大動脈として誤検出せずに、大動脈を正しく検出することができる。本形態では、スライスＳＬ５の大動脈を検出した後、スライスＳＬ５の大動脈の位置を基準にして、スライスＳＬ６〜ＳＬ１０の検索領域とスライスＳＬ４〜ＳＬ１の検索領域を設定する。したがって、スライスＳＬ１のように、肺を横切るスライスであっても、肺による低信号領域を大動脈として誤検出しないようにすることができる。
【０１２０】
本形態では、スライスＳＬ６〜ＳＬ１０の大動脈を検出した後（ステップＳＴ３３）、スライスＳＬ４〜ＳＬ１の大動脈を検出している（ステップＳＴ３４）。しかし、スライスＳＬ６〜ＳＬ１０の大動脈を検出するためのステップＳＴ３３と、スライスＳＬ４〜ＳＬ１の大動脈を検出するためのステップＳＴ３４とを並列で実行してもよい。ステップＳＴ３３とＳＴ３４とを並列で実行することにより、スライスＳＬ１〜ＳＬ１０の大動脈の検出に係る時間を短縮することができる。
【０１２１】
本形態では、スライスＳＬ５〜ＳＬ１０のうちのスライスＳＬ５を、大動脈を検出するときの基準となるスライスＳＬｊとして選択している。しかし、スライスＳＬ５〜ＳＬ１０の中から、スライスＳＬ５とは別のスライス（例えば、スライスＳＬ４）を、大動脈を検出するときの基準とスライスＳＬｊとして選択してもよい。ただし、スライスＳＬ１０のように、肝臓の下端に近いスライスには、大動脈Ａとは別の動脈の断面も現れやすいので、スライスＳＬ１０から大動脈Ａの検出を開始すると、別の動脈を大動脈Ａとして検出する可能性がある。したがって、スライスＳＬ５〜ＳＬ１０の中から、大動脈を検出するときの基準となるスライスＳＬｊを選択する場合、肝臓の下端から離れた位置のスライス（肺に近いスライス）を選択することが望ましい。
【０１２２】
本形態では、脳脊髄液ＣＳＦを基準にしてウィンドウＷを設定しているが、脳脊髄液ＣＳＦとは別の部位を基準にしてウィンドウＷを設定してもよい。また、本形態では大動脈Ａを検出しているが、本発明は、大動脈Ａとは別の血管を検出する場合にも適用することができる。
【０１２３】
本形態では、アキシャル画像に基づいて、大動脈を検出しているが、アキシャル面とは別の面（例えば、アキシャル面に対して斜めに交差するオブリーク面）の画像に基づいて大動脈を決定してもよい。
【０１２４】
本形態では、ウィンドウＷは矩形状であるが、別の形状（例えば、楕円形状）であってもよい。
【０１２５】
（２）第２の形態
第１の形態では、図１８に示すフローに従ってステップＳＴ３が実行されているが、第２の形態では、図１８に示すフローとは別のフローに従ってステップＳＴ３を実行する例について説明する。
【０１２６】
図４２は、第２の形態のＭＲ装置２００を概略的に示す図である。
尚、第２の形態のＭＲ装置２００は、第１の形態のＭＲ装置１００と比較すると、２値画像から腕部を削除するための削除手段９８を有する点が異なっているが、その他の構成は第１の形態のＭＲ装置１００と同じである。したがって、第２の形態のＭＲ装置２００の説明にあたっては、削除手段９８について主に説明する。
制御部９は、画像作成手段９１〜削除手段９８を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
【０１２７】
図４３は、第２の形態においてステップＳＴ３を実行するときのフローの説明図である。
ステップＳＴ３１０およびＳＴ３２０は第１の形態と同じであるので説明は省略する。
【０１２８】
ステップＳＴ３２０により、２値画像ＤＡ１（図２１参照）が得られる。尚、２値画像ＤＡ１には、被検体の胴部だけでなく、左腕部および右腕部も描出されている。第２の形態では、左腕部および右腕部を削除する処理をする。この処理を実行するために、ステップＳＴ３２１に進む。
【０１２９】
ステップＳＴ３２１では、削除手段９８が、２値画像ＤＡ１から、腕部を削除する処理を実行する（図４４参照）。
【０１３０】
図４４は、腕部が削除された２値画像ＤＣ１を概略的に示す図である。
ステップＳＴ３２１では、２値画像ＤＡ１から、左腕部と右腕部とを削除するための画像処理を実行する。これにより、腕部が削除された２値画像ＤＣ１を得ることができる。腕部を削除する方法としては、例えば、“Med Imag Tech Vol.31 No.2 March 2013”に記載されているにように、縮小処理、領域拡張処理、および拡大処理を実行する方法を用いることができる。腕部が削除された２値画像ＤＣ１を得た後、ステップＳＴ３３０に進む。
【０１３１】
ステップＳＴ３３０では、第２の２値化手段９３が、領域拡張法を用いずに、通常の閾値処理を用いて、スライスＳＬ１のアキシャル画像Ｄ１の２値化を実行する。この２値化は、第１の形態のステップＳＴ３３０で実行される２値化と同じである。図４５に、通常の閾値処理でスライスＳＬ１のアキシャル画像Ｄ１を２値化することにより得られた２値画像ＤＢ１を概略的に示す。２値画像ＤＢ１を得た後、ステップＳＴ３３１に進む。
【０１３２】
ステップＳＴ３３１では、削除手段９８が、２値画像ＤＢ１に対して、左腕部と右腕部とを削除するための画像処理を実行する。ステップＳＴ３３１で実行される画像処理は、ステップＳＴ３２１と同様に、縮小処理、領域拡張処理、および拡大処理を実行する方法を用いることができる。２値画像ＤＢ１は肺に論理値０が割り当てられているので、縮小処理、領域拡張処理、拡大処理を実行すると、腕部が削除されるだけでなく、胴部の領域が狭くなる。図４６に、腕部が削除された２値画像ＤＤ１を概略的に示す。腕部を削除した後、ステップＳＴ３４０に進む（図４７参照）。
【０１３３】
図４７は、ステップＳＴ３４０で実行される処理を説明するための図である。
ステップＳＴ３４０では、算出手段９４が、２値画像ＤＣ１の中の論理値１が割り当てられている領域の面積Ａと、２値画像ＤＤ１の中の論理値１が割り当てられている領域の面積Ｂとを計算する。面積ＡおよびＢの値を求めたら、算出手段９４は、面積比Ｒ＝Ｂ／Ａを計算する。
【０１３４】
図４７に示すように、面積Ｂは面積Ａよりもかなり小さいので、面積比Ｒは、１よりもかなり小さい値になることがわかる。
【０１３５】
このようにして、ステップＳＴ３０が実行される。尚、図４４〜図４７では、スライスＳＬ１のアキシャル画像Ｄ１に対してステップＳＴ３０を実行した例について説明されている。しかし、他のスライスＳＬ２〜ＳＬ１０の各々の画像に対しても、スライスＳＬ１のアキシャル画像Ｄ１と同様に、ステップＳＴ３０を実行し、面積比Ｒを求める。したがって、各スライスと面積比Ｒとの関係を表すプロファイルを得ることができる。図４８に、各スライスと面積比Ｒとの関係を表すプロファイルを示す。各スライスのアキシャル画像に対してステップＳＴ３０を実行した後、ステップＳＴ３１に進む。
ステップＳＴ３１〜ＳＴ３４は第１の形態と同じであるので、説明は省略する。
【０１３６】
第２の形態では、腕部が削除された２値画像を用いて面積比Ｒを求めている。面積比Ｒは、肝臓側では１に近い値になり、肺側では０に近い値になる。したがって、肝臓側の面積比と肺側の面積比との差が、第１の形態の場合（図２６参照）よりも大きくなるので、肺と肝臓との境界の検出精度を更に向上させることが可能となる。
【０１３７】
尚、第１および第２の形態では、スライスが肝臓を横切るか否かを判断するための指標として、面積比Ｒ＝Ａ／Ｂを用いている。しかし、大動脈Ａを検出するときの基準となるスライスＳＬｊを求めることができるのであれば、面積比Ｒとは別の指標を用いてもよい。例えば、面積比Ｒの代わりに、面積Ｂを、スライスが肝臓を横切るか否かを判断するための指標として用いてもよい。
【符号の説明】
【０１３８】
２  マグネット
３  テーブル
３ａ  クレードル
４  受信コイル
５  造影剤注入装置
６  送信器
７  勾配磁場電源
８  受信器
９  制御部
１０  操作部
１１  表示部
１２  被検体
２１  ボア
９１  画像作成手段
９２  第１の２値化手段
９３  第２の２値化手段
９４  算出手段
９５  選択手段
９６  検出手段
９７  設定手段
９８  削除手段
【図１】
2015123306.tif 000003
【図２】
2015123306.tif 000004
【図３】
2015123306.tif 000005
【図４】
2015123306.tif 000006
【図５】
2015123306.tif 000007
【図６】
2015123306.tif 000008
【図７】
2015123306.tif 000009
【図８】
2015123306.tif 000010
【図９】
2015123306.tif 000011
【図１０】
2015123306.tif 000012
【図１１】
2015123306.tif 000013
【図１２】
2015123306.tif 000014
【図１３】
2015123306.tif 000015
【図１４】
2015123306.tif 000016
【図１５】
2015123306.tif 000017
【図１６】
2015123306.tif 000018
【図１７】
2015123306.tif 000019
【図１８】
2015123306.tif 000020
【図１９】
2015123306.tif 000021
【図２０】
2015123306.tif 000022
【図２１】
2015123306.tif 000023
【図２２】
2015123306.tif 000024
【図２３】
2015123306.tif 000025
【図２４】
2015123306.tif 000026
【図２５】
2015123306.tif 000027
【図２６】
2015123306.tif 000028
【図２７】
2015123306.tif 000029
【図２８】
2015123306.tif 000030
【図２９】
2015123306.tif 000031
【図３０】
2015123306.tif 000032
【図３１】
2015123306.tif 000033
【図３２】
2015123306.tif 000034
【図３３】
2015123306.tif 000035
【図３４】
2015123306.tif 000036
【図３５】
2015123306.tif 000037
【図３６】
2015123306.tif 000038
【図３７】
2015123306.tif 000039
【図３８】
2015123306.tif 000040
【図３９】
2015123306.tif 000041
【図４０】
2015123306.tif 000042
【図４１】
2015123306.tif 000043
【図４２】
2015123306.tif 000044
【図４３】
2015123306.tif 000045
【図４４】
2015123306.tif 000046
【図４５】
2015123306.tif 000047
【図４６】
2015123306.tif 000048
【図４７】
2015123306.tif 000049
【図４８】
2015123306.tif 000050