(57)【要約】【課題】  ＣＰＲ画像による画像診断の効率向上を実現する画像処理装置の提供。【解決手段】  記憶部１４は、管状構造物に関するボリュームデータを記憶する。芯線抽出部１６２は、ボリュームデータに含まれる管状構造物の芯線を抽出する。ＣＰＲ画像発生部１６４は、抽出された芯線を含む複数の曲面に関する複数のＣＰＲ画像のデータをボリュームデータに基づいて発生する。最長画像特定部１８４は、発生された複数のＣＰＲ画像の中から、最長の画像長を有するＣＰＲ画像を特定する。最小傾き画像特定部１８８は、複数のＣＰＲ画像にそれぞれ含まれる複数の芯線について計算された複数の最大傾きのうち、最小傾きの芯線を含むＣＰＲ画像を特定する。表示部２０は、特定されたＣＰＲ画像を表示する。【選択図】      図１
2010051730.tif 000002
【特許請求の範囲】
【請求項１】
管状構造物に関するボリュームデータを記憶する記憶部と、
前記ボリュームデータに含まれる管状構造物の芯線を抽出する抽出部と、
前記抽出された芯線を含む複数の曲面に関する複数のＣＰＲ画像のデータを前記ボリュームデータに基づいて発生する画像発生部と、
前記発生された複数のＣＰＲ画像の中から、所定の画像長を有する第１ＣＰＲ画像又は所定の傾きを有する前記芯線を含む第２ＣＰＲ画像を特定する画像特定部と、
前記特定された第１ＣＰＲ画像又は第２ＣＰＲ画像を表示する表示部と、
を具備する画像処理装置。
【請求項２】
前記画像特定部は、
前記複数のＣＰＲ画像の複数の画像長をそれぞれ計算する画像長計算部と、
前記計算された複数の画像長のうちの最長の画像長を有する前記第１ＣＰＲ画像を前記複数のＣＰＲ画像から特定する最長画像特定部と、
を備える請求項１記載の画像処理装置。
【請求項３】
前記ＣＰＲ画像発生部は、前記芯線上の点列を一定の向きで通過する複数の切断線により前記曲面を構成し、
前記画像長計算部は、前記切断線の直交方向に関する前記ＣＰＲ画像の長さを前記画像長として計算する、
請求項２記載の画像処理装置。
【請求項４】
前記画像特定部は、
前記複数のＣＰＲ画像上の複数の芯線の最大傾きをそれぞれ計算する傾き計算部と、
前記計算された複数の最大傾きのうちの最小傾きを有する芯線を含む前記第２ＣＰＲ画像を前記複数のＣＰＲ画像から特定する最小傾き画像特定部と、
を備える請求項１記載の画像処理装置。
【請求項５】
前記ＣＰＲ画像発生部は、前記芯線上の点列を一定の向きで通過する複数の切断線により前記曲面を構成し、
前記傾き計算部は、前記芯線の複数個所のそれぞれについて、前記切断線の直交方向に対する傾きをそれぞれ計算し、前記計算された各傾きのうちの最大値を、前記最大傾きとして特定する、
請求項４記載の画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＣＰＲ（Ｃｕｒｖｅｄ  ｍｕｌｔｉ―Ｐｌａｎａｒ  Ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ）画像のデータを発生する画像処理装置に関する。
【背景技術】
【０００２】
Ｘ線コンピュータ断層撮影装置や磁気共鳴イメージング装置等で得られたボリュームデータに含まれる管状構造物を観察する方法としてＣＰＲ法が知られている（例えば非特許文献１参照）。ＣＰＲ法では、Ｘ軸、Ｙ軸、あるいはＺ軸に平行な切断線を初期的に設定する。そして、設定された切断線で管状構造物の芯線に沿ってボリュームデータを切断することにより形成される曲断面に関するＣＰＲ画像を初期的に表示している。しかし、切断線に平行な方向に走行する部分をＣＰＲ画像で観察する場合、その位置関係が把握しづらく、また、平行な部分の画像が乱れてしまう。この際ユーザは、切断線の向きを観察に適した向きに手動で設定し直さなければならない。このため、画像診断時間の遅れが生じている。
【非特許文献１】“ＣＰＲ―ｃｕｒｖｅｄ  ｐｌａｎａｒ  ｒｅｆｏｒｍａｔｉｏｎ”  ｖｉｓｕａｌｉｚａｔｉｏｎ  ２００２  ＶＩＳ２００２  ＩＥＥＥ  Ｖｏｌ  Ｉｓｓｕｅ，  １―１  Ｎｏｖ．  ２００２．  Ｐ３７―４４
【発明の開示】
【発明が解決しようとする課題】
【０００３】
本発明の目的は、ＣＰＲ画像による画像診断の効率向上を実現する画像処理装置を提供することにある。
【課題を解決するための手段】
【０００４】
本発明のある局面に係る画像処理装置は、管状構造物に関するボリュームデータを記憶する記憶部と、前記ボリュームデータに含まれる管状構造物の芯線を抽出する抽出部と、前記抽出された芯線を含む複数の曲面に関する複数のＣＰＲ画像のデータを前記ボリュームデータに基づいて発生する画像発生部と、前記発生された複数のＣＰＲ画像の中から、所定の画像長を有する第１ＣＰＲ画像又は所定の傾きを有する前記芯線を含む第２ＣＰＲ画像を特定する画像特定部と、前記特定された第１ＣＰＲ画像又は第２ＣＰＲ画像を表示する表示部と、を具備する。
【発明の効果】
【０００５】
本発明によれば、ＣＰＲ画像による画像診断の効率向上が実現する。
【発明を実施するための最良の形態】
【０００６】
以下、図面を参照しながら本発明の実施形態に係わる画像処理装置を説明する。
【０００７】
図１は、本発明の実施形態に係わる画像処理装置の構成を示す図である。図１に示すように画像処理装置１は、制御部１０を中枢として、入力部１２、記憶部１４、ＣＰＲ処理部１６、画像特定部１８、及び表示部２０を備える。
【０００８】
入力部１２は、ユーザからの各種指令や情報入力を受け付ける。入力部１２としては、マウスやトラックボールなどのポインティングデバイス、スイッチボタン等の選択デバイス、あるいはキーボード等の入力デバイスが適宜利用可能である。
【０００９】
記憶部１４は、Ｘ線コンピュータ断層撮影装置（Ｘ線ＣＴ装置）や磁気共鳴イメージング装置（ＭＲＩ装置）、超音波診断装置、核医学診断装置等によって発生された管状構造物に関するボリュームデータを記憶する。管状構造物は、血管、腸、食道等のあらゆる管状の構造物を指す。以下、説明を具体的に行なうため管状構造物は、造影剤により造影された血管であるとする。ボリュームデータは、スライス方向に沿って配列される複数のスライス画像のデータにより構成される。
【００１０】
ＣＰＲ処理部１６は、記憶部１４に記憶されているボリュームデータを読み込んでＣＰＲ処理をする。図１に示すように、ＣＰＲ処理部１６は、芯線抽出部１６２とＣＰＲ画像発生部１６４とを有する。芯線抽出部１６２は、既存の技術を用いて、ボリュームデータ内の血管の芯線を抽出する。ＣＰＲ画像発生部１６４は、芯線抽出部１６２により抽出された芯線を含む複数の曲面に関する複数のＣＰＲ画像のデータをボリュームデータから発生する。ＣＰＲ処理部１６によるＣＰＲ処理の詳細については、後述する。
【００１１】
画像特定部１８は、ＣＰＲ処理部１６により発生された複数のＣＰＲ画像の中から最も画像の乱れの少ないＣＰＲ画像を特定する。図１に示すように、画像特定部１８は、画像長計算部１８２、最長画像特定部１８４、傾き計算部１８６、及び最小傾き画像特定部１８８を有する。画像長計算部１８２は、ＣＰＲ処理部１６により発生された複数のＣＰＲ画像の複数の画像長をそれぞれ計算する。最長画像特定部１８４は、複数のＣＰＲ画像の中から最長の画像長を有するＣＰＲ画像を特定する。以下、最長の画像長を有するＣＰＲ画像を、最長ＣＰＲ画像と呼ぶことにする。傾き計算部１８６は、複数のＣＰＲ画像に含まれる複数の芯線の最大傾きをそれぞれ計算する。最小傾き画像計算部１８８は、複数のＣＰＲ画像の中から、複数の最大傾きのうちの最小傾きを有する芯線を含むＣＰＲ画像を特定する。以下、複数の最大傾きのうちの最小傾きを有する芯線を含むＣＰＲ画像を、最小傾きＣＰＲ画像と呼ぶことにする。画像特定部１８による画像特定処理の詳細については、後述する。
【００１２】
表示部２０は、画像特定部１８により特定された最長ＣＰＲ画像や最小傾きＣＰＲ画像を表示する。表示部２０としては、例えばＣＲＴディスプレイや、液晶ディスプレイ、有機ＥＬディスプレイ、プラズマディスプレイ等が利用可能である。
【００１３】
制御部１０は、画像処理装置１の全体を制御する。制御部１０は、入力部１２からの最長ＣＰＲ画像の表示処理の開始要求信号を受けると、画像処理装置１の各部を制御して、最長ＣＰＲ画像の表示処理を実行する。また、制御部１０は、入力部１２からの最小傾きＣＰＲ画像の表示処理の開始要求信号を受けると、画像処理装置１の各部を制御して、最小傾きＣＰＲ画像の表示処理を実行する。
【００１４】
次にＣＰＲ画像発生処理について説明する。図２は、ＣＰＲ画像発生処理の原理を説明するための図である。なお図２に示すＸＹＺ直交座標系のＺ軸は、スライス方向（被検体の体軸）に規定されているものとする。図２に示すように、ＣＰＲ画像発生処理は、ボリュームデータ３０内を蛇行する血管３１の全体を、血管３１の芯線３２を利用して重なり無く一平面状に可視化することを特徴としている。
【００１５】
図２に示すように、ＣＰＲ画像４０の断面３３は、曲面である。この曲断面３３は、芯線３２を構成する複数の点をそれぞれ貫く複数の切断線３４により構成される。このように、曲断面３３は、ボリュームデータ３０内に、芯線３２に沿うように設定される。一枚のＣＰＲ画像４０に関する全切断線３４の向きは、同一である。切断線３４の向きは、Ｘ軸、Ｙ軸、及びＺ軸のそれぞれに対する角度で規定される。典型的には、切断線３４は、Ｚ軸に直交する平面（ＸＹ平面）内に設定される。図２における切断線３４は、ＸＹ平面であって、Ｘ軸に平行に設定されている。これら複数の切断線３４から構成される曲断面３３が平面上に展開されることでＣＰＲ画像４０が得られる。
【００１６】
ここで、図２に示すように、曲断面３３上にＡＢ座標系を規定する。Ａ軸は、切断線３４に対して平行に規定される。Ｂ軸は、曲断面３３上に沿ってＡ軸に対して直交に規定される。すなわち、ＡＢ座標系をＸＹＺ座標系から観察すると、Ｂ軸は、芯線３２と同じようにＹ軸に関して曲がって見える。このＡ軸とＢ軸とを用いると、ＣＰＲ画像４０の座標系は、Ｂ軸に平行なｘ軸と、Ａ軸をＺ軸に平行に直線化したｚ軸とにより規定される。
【００１７】
次に図３と図４とを参照しながら、ＣＰＲ画像の乱れについて説明する。図３に示すように、血管径Ｒを備える血管３１の芯線３２上の複数の点についてそれぞれ複数の切断線３４が設定されているとする。各切断線３４は、Ｘ軸に平行に設定されているとする。また、切断線３４は、芯線３２上に一定間隔ｐ毎に設定される。血管３１は、切断線３４に対して平行になるように曲がっている。ここで、切断線３４に直交する血管部分３１１における芯線３４上の点Ｃ１と、切断線３４に平行するように曲がっている血管部分３１２における芯線３４上の点Ｃ２とを考える。点Ｃ２における芯線３２の走行方向ｄ２は、点Ｃ１における芯線３２の走行方向ｄ１に比して、切断線３４に対する角度が鋭角である。従って、点Ｃ１に設定された切断線３４における血管３１の幅Ｄ１は、血管径Ｒと同等であるが、点Ｃ２に設定された切断線３４における血管３１の幅Ｄ２は、血管径Ｒよりも大きい。
【００１８】
図４は、図３において示されている複数の切断線３４により構成される断曲面に関するＣＰＲ画像４０を示す。図４に示すように、図３の切断線３４に直交する血管部分３１１に由来するＣＰＲ画像４０上の血管部分３１３は、画像に乱れは生じない。しかし、切断線３４に平行するように曲がっている血管部分３１２に由来するＣＰＲ画像４０上の血管部分３１４は、膨張して見えてしまう。すなわち、芯線３２の走行方向が切断線３４に対して鋭角になればなるほど、画像に乱れが生じてしまう。
【００１９】
本実施形態に係わる画像処理装置１は、互いに切断線の向きの異なる複数のＣＰＲ画像のデータを発生し、発生された複数のＣＰＲ画像の中から、切断線に対して非直交方向に走行する血管部分に由来する画像の乱れが最も少ないＣＰＲ画像を表示する。
【００２０】
以下、制御部１０の制御により行なわれる、最長ＣＰＲ画像の表示処理を説明する。図５は、最長ＣＰＲ画像の表示処理の流れを示す図である。
【００２１】
ユーザにより入力部１２を介して最長ＣＰＲ画像の表示処理の開始指示がなされると、制御部１０は、ＣＰＲ処理部１６の芯線抽出部１６２に芯線抽出処理を行なわせる（ステップＳＡ１）。
【００２２】
芯線抽出処理において芯線抽出部１６２は、既存の技術を用いて、ボリュームデータから血管の芯線を抽出する。この際、芯線は、点列として抽出される。芯線を抽出する技術としては、例えば、ベッセルトラッキング法（例えばＯｎｎｏ  Ｗｉｎｋ，  Ｗｉｒｏ  Ｊ．  Ｎｉｅｓｓｅｎ，  “Ｆａｓｔ  Ｄｅｌｉｎａｔｉｏｎ  ａｎｄ  Ｖｉｓｕａｌｉｚａｔｉｏｎ  ｏｆ  Ｖｅｓｓｅｌ  ｉｎ  ３―Ｄ  Ａｎｇｉｏｇｒａｐｈｉｃ  Ｉｍａｇｅｓ”，  ＩＥＥＥ  Ｔｒａｎｓ．  Ｍｅｄ．  Ｉｍａｇ．  Ｖｏｌ．１９，  Ｎｏ．４，  ２０００）や、管腔臓器の内部領域を細線化する方法（例えばＧ．  Ｄ．  Ｒｕｂｉｎ，  Ｄ．  Ｓ．  Ｐａｉｋ，  Ｐ．  Ｃ．  Ｊｏｈｎｓｔｏｎ，  Ｓ．  Ｎａｐｅｌ，  “Ｍｅａｓｕｒｍｅｎｔ  ｏｆ  ｔｈｅ  Ａｏｒｔａ  ａｎｄ  Ｉｔｓ  Ｂｒａｎｃｈｅｓ  ｗｉｔｈ  Ｈｅｌｉｃａｌ  ＣＴ，”  Ｒａｄｉｏｌｏｇｙ，  Ｖｏｌ．２０６，  Ｎｏ．３，  ｐｐ．８２３―９，  Ｍａｒ．，  １９９８）等が適用可能である。
【００２３】
芯線が抽出されると制御部１０は、ＣＰＲ処理部１６のＣＰＲ画像発生部１６４に初期の切断線の向きの設定処理を行なわせる（ステップＳＡ２）。初期の切断線の向きの設定処理においてＣＰＲ画像発生部１６４は、ユーザにより入力部１２を介して指定された向きを、又は予め設定された向きを初期の切断線の向きθ１に設定する。
【００２４】
初期の切断線の向きθ１が設定されると制御部１０は、ＣＰＲ画像発生部１６４にＣＰＲ画像の発生処理を行なわせる（ステップＳＡ３）。ＣＰＲ画像の発生処理においてＣＰＲ画像発生部１６４は、設定された向きθ１の複数の切断線が芯線を構成する複数の点をそれぞれ通過するように、ボリュームデータ内に曲断面を設定する。そしてＣＰＲ画像発生部１６４は、設定した曲断面を平面上に展開してＣＰＲ画像のデータを発生する。なお、切断線は、芯線を構成する全点について設定されるとしても良いし、数点おきに設定されるとしても良い。ＣＰＲ画像のｚ軸方向に関する解像度が不足している場合、不足分は、切断線上の各画素の画素値に基づいて補間可能である。
【００２５】
ＣＰＲ画像のデータが発生されると制御部１０は、画像特定部１８の画像長計算部１８２に画像長の計算処理を行なわせる（ステップＳＡ４）。画像長の計算処理において画像長計算部１８２は、発生されたＣＰＲ画像の画像長を計算する。ＣＰＲ画像のデータは、その画像長に関連付けられて記憶部１４に記憶される。
【００２６】
このステップＳＡ３とステップＳＡ４との処理は、切断線の向きを変えながら複数回繰り返される。切断線は、図６に示すように、所定の回転軸ＲＡ回りに向きθ１から＋１８０°まで一定の角度間隔（例えば５°）で順番に設定される。この角度間隔は、予め設定されていても良いし、入力部１２を介して任意に設定可能としてもよい。
【００２７】
ステップＳＡ４において画像長が計算されると制御部１０は、向き（θ１＋１８０°）の切断線に関するＣＰＲ画像について画像長が計算されたか否かを判定する（ステップＳＡ５）。向き（θ１＋１８０°）の切断線に関するＣＰＲ画像について画像長が計算されていない（ステップＳＡ５：ＮＯ）と判定すると制御部１０は、ＣＰＲ画像発生部１６４に、切断線を所定の回転軸回りに一定角度回転して次の切断線を設定させる（ステップＳＡ６）。そしてステップＳＡ６において設定された切断線で、ＣＰＲ画像発生部１６４によるステップＳＡ３と、画像長計算部１８２によるステップＳＡ４とが繰り返される。
【００２８】
一方ステップＳＡ５において、向き（θ１＋１８０°）の切断線に関するＣＰＲ画像について画像長が計算された（ステップＳＡ５：ＹＥＳ）と判定すると制御部１０は、画像特定部１８の最長画像特定部１８４に画像特定処理を行なわせる（ステップＳＡ６）。
【００２９】
画像特定処理において最長画像特定部１８４は、発生された複数のＣＰＲ画像の中から、最長ＣＰＲ画像を特定する。図７は、切断線の向きが互いに異なる２つのＣＰＲ画像４０Ａ，４０Ｂとその画像長ＬＡ，ＬＢとを示す図である。図７に示すように、切断線の向きθＡに関するＣＰＲ画像４０Ａに描出されている血管３１Ａは、切断線の向きθＢに関するＣＰＲ画像４０Ｂに描出されている血管３１Ｂに比して、画像上で大きく曲がっている。従って、ＣＰＲ画像４０Ａの画像長３１Ａは、ＣＰＲ画像４０Ｂの画像長３１Ｂに比して小さい。このように、ＣＰＲ画像の画像長は、切断線の向きに応じて変化する。画像長が長いＣＰＲ画像は、画像長が短いＣＰＲ画像に比して、切断線に対して非直交方向に走行する血管部分が少ないといえる。すなわち、画像長が短ければ短いほど、切断線に対して非直交方向に走行する血管部分に由来する画像の乱れが少ない。つまり、最小の画像長を有する最長ＣＰＲ画像は、発生された複数のＣＰＲ画像のうち、画像の乱れが最も少ない。
【００３０】
最長ＣＰＲ画像が特定されると制御部１０は、表示部２０に表示処理を行なわせる（ステップＳＡ７）。この表示処理において表示部２０は、特定された最長ＣＰＲ画像を記憶部１４から読み出して表示する。
【００３１】
次に、制御部１０の制御により行なわれる、最小傾きＣＰＲ画像の表示処理を説明する。図８は、最小傾きＣＰＲ画像の表示処理の流れを示す図である。
【００３２】
ユーザにより入力部１２を介して最小傾きＣＰＲ画像の表示処理の開始指示がなされると、制御部１０は、ＣＰＲ処理部１６の芯線抽出部１６２に芯線抽出処理を行なわせる（ステップＳＢ１）。
【００３３】
芯線が抽出されると制御部１０は、ＣＰＲ処理部１６のＣＰＲ画像発生部１６４に初期の切断線の向きの設定処理を行なわせる（ステップＳＢ２）。初期の切断線の向きの設定処理においてＣＰＲ画像発生部１６４は、ユーザにより入力部１２を介して指定された向きを、又は予め設定された向きを初期の切断線の向きθ１に設定する。
【００３４】
初期の切断線の向きθ１が設定されると制御部１０は、ＣＰＲ画像発生部１６４にＣＰＲ画像の発生処理を行なわせる（ステップＳＢ３）。ＣＰＲ画像の発生処理においてＣＰＲ画像発生部１６４は、設定された向きθ１の複数の切断線が構成する曲断面に関するＣＰＲ画像のデータを発生する。
【００３５】
ＣＰＲ画像のデータが発生されると制御部１０は、画像特定部１８の傾き計算部１８６に傾き計算処理を行なわせる（ステップＳＢ４）。傾き計算処理において傾き計算部１８６は、ＣＰＲ画像上の芯線を微小間隔で区切り、各微小区間における芯線の傾きを計算する。より詳細には、微小区間の両端点を結ぶ線（接線）の傾きが、当該微小区間における芯線の傾きとして計算される。なお、傾きの計算は、上記方法に限定されない。例えば、傾き計算部は、芯線上の複数点（複数の芯点）に基づいて複数点を曲線に近似し、この近似曲線を微小間隔で区切って、各区間における近似曲線の接線の傾きを、芯線の接線の傾きとして計算してもよい。この曲線近似は、スプライン補間やラグランジュ補間等の既存のどの方法を用いてもよい。
【００３６】
各微小区間における芯線の傾きが計算されると制御部１０は、傾き計算部１８６に最大傾きの特定処理を行なわせる（ステップＳＢ５）。最大傾きの特定処理において傾き特定部１８６は、計算された複数の傾きの中から最大の傾きを、当該ＣＰＲ画像の芯線の最大傾きとして特定する。ＣＰＲ画像のデータは、その最大傾きに関連付けられて記憶部１４に記憶される。
【００３７】
ステップＳＢ５において最大傾きが特定されると制御部１０は、向き（θ１＋１８０°）の切断線に関するＣＰＲ画像について最大傾きが特定されたか否かを判定する（ステップＳＢ６）。向き（θ１＋１８０°）の切断線に関するＣＰＲ画像について最大傾きが特定された特定されていない（ステップＳＢ６：ＮＯ）と判定すると制御部１０は、ＣＰＲ画像発生部１６４に、切断線を所定の回転軸回りに一定角度回転して次の切断線を設定させる（ステップＳＢ７）。そしてステップＳＢ７において設定された切断線で、ＣＰＲ画像発生部１６４によるステップＳＢ３と、傾き計算部１８６によるステップＳＢ４と、傾き計算部１８６によるステップＳＢ６とが順番に繰り返される。
【００３８】
一方ステップＳＢ６において、向き（θ１＋１８０°）の切断線に関するＣＰＲ画像について最大傾きが特定された（ステップＳＢ６：ＹＥＳ）と判定すると制御部１０は、画像特定部１８の最小傾き画像特定部１８８に画像特定処理を行なわせる（ステップＳＢ８）。画像特定処理において最小傾き画像特定部１８８は、複数の最大傾きのうちの最小傾きの芯線を含む最小傾きＣＰＲ画像を、複数のＣＰＲ画像の中から特定する。
【００３９】
図９は、切断線の向きが互いに異なる２つのＣＰＲ画像の芯線３２Ａ，３２Ｂと傾きとを示す図である。図９に示すように、傾きは、切断線の直交方向の点間距離と平行方向の点間距離との比で計算される。より詳細には、傾きは、切断線の直交方向の点間距離に対する平行方向の点間距離で計算される。切断線の向きθＡに関するＣＰＲ画像上の芯線３２Ａの最大傾きｍａｘ［ＴＡ（ｍ）］＝ＴＡ（５）は、点ＰＡ（５）と点ＰＡ（６）との間隔で計算される。また、切断線の向きθＢに関するＣＰＲ画像上の芯線３２Ｂの最大傾きｍａｘ［ＴＢ（ｍ）］＝ＴＢ（１）は、点ＰＢ（１）と点ＰＢ（２）との間隔で計算される。また、切断線の向きθＡに関するＣＰＲ画像上の芯線３２Ａは、切断線の向きθＢに関するＣＰＲ画像上の芯線３２Ｂに比して、ｘ軸方向に大きく曲がっている。従って、芯線３２Ａの最大傾きＴＡ（５）は、芯線３２Ｂの最大傾きＴＢ（５）に比して大きい。
【００４０】
最大傾きＴＡ（５）が計算される点ＰＡ（５）と点ＰＡ（６）との間の芯線３２Ａ部分を観察すれば明らかなように、傾きの鋭い芯線部分は、傾きの緩い芯線部分に比して、よりｘ軸すなわち切断線に対して鋭角に走行する。すなわち、含まれる芯線の最大傾きが小さければ小さいほど、切断線に対して非直交方向に走行する血管部分に由来する画像の乱れは少ない。つまり、最小の最大傾きの芯線を有するＣＰＲ画像は、発生された複数のＣＰＲ画像のうち、画像の乱れが最も少ない。
【００４１】
最小の芯線の最大傾きを有するＣＰＲ画像が特定されると制御部１０は、表示部２０に表示処理を行なわせる（ステップＳＢ８）。この表示処理において表示部２０は、特定されたＣＰＲ画像を記憶部１４から読み出して表示する。
【００４２】
このように、本実施形態に係わる画像処理装置１は、複数の曲断面に関する複数のＣＰＲ画像のデータを発生し、発生された複数のＣＰＲ画像の中から最長ＣＰＲ画像又は最小傾きＣＰＲ画像を自動的に特定して初期的に表示する。この最長ＣＰＲ画像又は最小傾きＣＰＲ画像は、複数のＣＰＲ画像のうち、切断線の非直交方向に走行する血管部分に由来する画像の乱れが最も少ない。かくして本実施形態によれば、ＣＰＲ画像による画像診断の効率向上が実現する。
【００４３】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【図面の簡単な説明】
【００４４】
【図１】本発明の実施形態に係わる画像処理装置の構成を示す図。
【図２】ＣＰＲ画像発生処理の原理を説明するための図。
【図３】ＣＰＲ画像の乱れを説明するための図。
【図４】ＣＰＲ画像の乱れを説明するための他の図。
【図５】図１の制御部の制御のもとに行なわれる、最長ＣＰＲ画像の表示処理の流れを示す図。
【図６】図５のステップＳＡ２及びステップＳＡ６における切断線の設定処理を説明するための図。
【図７】図５のステップＳＡ７における最長ＣＰＲ画像の特定処理を説明するための図。
【図８】図１の制御部の制御のもとに行なわれる、最小傾きＣＰＲ画像の表示処理の流れを示す図。
【図９】図８のステップＳＢ８における最小傾きＣＰＲ画像の特定処理を説明するための図。
【符号の説明】
【００４５】
１…画像処理装置、１０…制御部、１２…入力部、１４…記憶部、１６…ＣＰＲ処理部、１６２…芯線抽出部、１６４…ＣＰＲ画像発生部、１８…画像特定部、１８２…画像長計算部、１８４…最長画像特定部、１８６…傾き計算部、１８８…最小傾き画像特定部、２０…表示部
【図１】
2010051730.tif 000003
【図２】
2010051730.tif 000004
【図３】
2010051730.tif 000005
【図４】
2010051730.tif 000006
【図５】
2010051730.tif 000007
【図６】
2010051730.tif 000008
【図７】
2010051730.tif 000009
【図８】
2010051730.tif 000010
【図９】
2010051730.tif 000011