(57)【要約】【課題】脳出血部位の体積計算を適正に行う方法および装置を実現する。【解決手段】脳出血を発症している頭部をＸ線ＣＴ装置で撮影して得られる連続スライスの複数の画像について正常な脳実質部に関するＲＯＩと脳出血部に関するＲＯＩをそれぞれ設定し（３０５，３０７）、正常な脳実質部に関するＲＯＩ内の画素のＣＴ値の平均値を複数の画像についてそれぞれ求め（３１７）、１つの画像の平均値を基準として他の画像の平均値のシフト量を画像ごとに求め（３１９）、画像ごとのシフト量を用いて脳出血部に関するＲＯＩ内の画素のＣＴ値を画像ごとに調整し（３２１）、調整後のＣＴ値について複数の画像を通じての最大値を特定し（３２３）、脳出血部位の体積を計算する（３２５）。【選択図】図３
2007152105.tif 000002
【特許請求の範囲】
【請求項１】
脳出血を発症している頭部をＸ線ＣＴ装置で撮影して得られる連続スライスの複数の画像について正常な脳実質部に関するＲＯＩと脳出血部に関するＲＯＩをそれぞれ設定し、
正常な脳実質部に関するＲＯＩ内の画素のＣＴ値の平均値を複数の画像についてそれぞれ求め、
１つの画像の平均値を基準として他の画像の平均値のシフト量を画像ごとに求め、
画像ごとのシフト量を用いて脳出血部に関するＲＯＩ内の画素のＣＴ値を画像ごとに調整し、
調整後のＣＴ値について複数の画像を通じての最大値を特定し、
次式によって脳出血部位の体積を計算する、
【数１】
2007152105.tif 000015
ここで、
ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての画素番号
Ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての総画素数
ｐｎ  ：脳出血部に関するＲＯＩ内の画素の調整後のＣＴ値
ＤＦＯＶ：デファインドＦＯＶの大きさ
ＳＦＯＶ：スキャンＦＯＶの大きさ
ＳｌｉｃｅＴｈｉｃｋｎｅｓｓ：スライス厚
Ｍａｇｎｉｆｉｃａｔｉｏｎ：倍率
ことを特徴とする脳出血部位体積計算方法。
【請求項２】
脳出血部に関するＲＯＩに含まれる骨部を除外する、
ことを特徴とする請求項１に記載の脳出血部位体積計算方法。
【請求項３】
脳出血部に関するＲＯＩに含まれる異物部を除外する、
ことを特徴とする請求項１または請求項２に記載の脳出血部位体積計算方法。
【請求項４】
脳出血を発症している頭部をＸ線ＣＴ装置で撮影して得られる連続スライスの複数の画像について正常な脳実質部に関するＲＯＩと脳出血部に関するＲＯＩをそれぞれ設定する設定手段と、
正常な脳実質部に関するＲＯＩ内の画素のＣＴ値の平均値を複数の画像についてそれぞれ求める平均値算出手段と、
１つの画像の平均値を基準として他の画像の平均値のシフト量を画像ごとに求めるシフト量算出手段と、
画像ごとのシフト量を用いて脳出血部に関するＲＯＩ内の画素のＣＴ値を画像ごとに調整する調整手段と、
調整後のＣＴ値について複数の画像を通じての最大値を特定する特定手段と、
次式によって脳出血部位の体積を計算する計算手段と、
【数２】
2007152105.tif 000016
ここで、
ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての画素番号
Ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての総画素数
ｐｎ  ：脳出血部に関するＲＯＩ内の画素の調整後のＣＴ値
ＤＦＯＶ：デファインドＦＯＶの大きさ
ＳＦＯＶ：スキャンＦＯＶの大きさ
ＳｌｉｃｅＴｈｉｃｋｎｅｓｓ：スライス厚
Ｍａｇｎｉｆｉｃａｔｉｏｎ：倍率
を具備することを特徴とする脳出血部位体積計算装置。
【請求項５】
脳出血部に関するＲＯＩに含まれる骨部を除外する骨部除外手段、
を具備することを特徴とする請求項４に記載の脳出血部位体積計算装置。
【請求項６】
脳出血部に関するＲＯＩに含まれる異物部を除外する異物部除外手段、
を具備することを特徴とする請求項４または請求項５に記載の脳出血部位体積計算装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脳出血部位体積計算（Ｃａｌｃｕｌａｔｉｎｇ  ３−Ｄ  Ｖｏｌｕｍｅ  ｏｆ  Ｃｅｒｅｂｒａｌ  Ｈｅｍｏｒｒｈａｇｅ）方法および装置に関し、特に、脳出血を発症している頭部のＸ線ＣＴ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）画像に基づいて脳出血部位の体積を計算する方法および装置に関する。
【背景技術】
【０００２】
脳出血の診断および治療のために、頭部のＸ線ＣＴ画像に基づいて脳出血部位の体積計算が行われる。脳出血部位体積計算は、先ず画像上で脳出血部位を特定し、次に特定された部位の体積を計算することによって行われる。計算の方式としては、脳出血部位の特定に専門医が関与する半自動方式（例えば、特許文献１参照）と、脳出血部位特定から体積計算まで全て自動で行う全自動方式がある。
【特許文献１】特開２００５−１１８５１０号公報（第６−７頁、図３−８）
【発明の開示】
【発明が解決しようとする課題】
【０００３】
半自動方式は、脳出血部位の特定（セグメンテーション：  ｓｅｇｍｅｎｔａｔｉｏｎ）に専門医が関与するので労力と時間がかかる。また、セグメンテーションの精度が担当者の技量に左右されるので、得られる結果は担当者ごとに異なり、同一担当者による再現性も良くない。
【０００４】
これに対して、全自動方式は省力的ではあるが、セグメンテーションの精度に問題があり、十分な成果は得られていない。すなわち、上部脳領域（ｕｐｐｅｒ  ｃｅｒｅｂｒａｌ  ｓｅｃｔｉｏｎ）では過小セグメンテーション（ｕｎｄｅｒ  ｓｅｇｍｅｎｔａｔｉｏｎ）となる傾向があり、中下部脳領域（ｍｉｄｄｌｅ  ｏｒ  ｌｏｗｅｒ  ｃｅｒｅｂｒａｌ  ｓｅｃｔｉｏｎ）では過剰セグメンテーション（ｏｖｅｒ  ｓｅｇｍｅｎｔａｔｉｏｎ）となる傾向がある。さらには、脳出血部でないところをセグメンテーションしたり、脳出血部をセグメンテーションし損なったりすることもある。
【０００５】
また、いずれの方式も、パーシャルヴォリューム・エフェクト（Ｐａｒｔｉａｌ  Ｖｏｌｕｍｅ  Ｅｆｆｅｃｔ）および複数画像間のＣＴ値シフト（ｓｈｉｆｔ）を無視しているので、計算結果は不正確なものとならざるを得ない。
【０００６】
そこで、本発明の課題は、脳出血部位の体積計算を適正に行う方法および装置を実現することである。また、本発明の課題は、ユーザーの技量に影響されない安定な結果が得られる脳出血部位計算方法および装置を実現することである。
【０００７】
また、本発明の課題は、ユーザーにとって操作しやすい脳出血部位計算方法および装置を実現することである。また、本発明の課題は、上中下いずれの脳領域にも適応し、パーシャルヴォリューム・エフェクトおよび複数画像間のＣＴ値シフトの問題も克服した脳出血部位計算方法および装置を実現することである。
【課題を解決するための手段】
【０００８】
上記の課題を解決するための第１の観点での発明は、脳出血を発症している頭部をＸ線ＣＴ装置で撮影して得られる連続スライスの複数の画像について正常な脳実質部に関するＲＯＩと脳出血部に関するＲＯＩをそれぞれ設定し、正常な脳実質部に関するＲＯＩ内の画素のＣＴ値の平均値を複数の画像についてそれぞれ求め、１つの画像の平均値を基準として他の画像の平均値のシフト量を画像ごとに求め、画像ごとのシフト量を用いて脳出血部に関するＲＯＩ内の画素のＣＴ値を画像ごとに調整し、調整後のＣＴ値について複数の画像を通じての最大値を特定し、次式によって脳出血部位の体積を計算する、
【０００９】
【数３】
2007152105.tif 000003
【００１０】
ここで、
ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての画素番号
Ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての総画素数
ｐｎ  ：脳出血部に関するＲＯＩ内の画素の調整後のＣＴ値
ＤＦＯＶ：デファインドＦＯＶの大きさ
ＳＦＯＶ：スキャンＦＯＶの大きさ
ＳｌｉｃｅＴｈｉｃｋｎｅｓｓ：スライス厚
Ｍａｇｎｉｆｉｃａｔｉｏｎ：倍率
ことを特徴とする脳出血部位体積計算方法である。
【００１１】
上記の課題を解決するための第２の観点での発明は、脳出血を発症している頭部をＸ線ＣＴ装置で撮影して得られる連続スライスの複数の画像について正常な脳実質部に関するＲＯＩと脳出血部に関するＲＯＩをそれぞれ設定する設定手段と、正常な脳実質部に関するＲＯＩ内の画素のＣＴ値の平均値を複数の画像についてそれぞれ求める平均値算出手段と、１つの画像の平均値を基準として他の画像の平均値のシフト量を画像ごとに求めるシフト量算出手段と、画像ごとのシフト量を用いて脳出血部に関するＲＯＩ内の画素のＣＴ値を画像ごとに調整する調整手段と、調整後のＣＴ値について複数の画像を通じての最大値を特定する特定手段と、次式によって脳出血部位の体積を計算する計算手段と、
【００１２】
【数４】
2007152105.tif 000004
【００１３】
ここで、
ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての画素番号
Ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての総画素数
ｐｎ  ：脳出血部に関するＲＯＩ内の画素の調整後のＣＴ値
ＤＦＯＶ：デファインドＦＯＶの大きさ
ＳＦＯＶ：スキャンＦＯＶの大きさ
ＳｌｉｃｅＴｈｉｃｋｎｅｓｓ：スライス厚
Ｍａｇｎｉｆｉｃａｔｉｏｎ：倍率
を具備することを特徴とする脳出血部位体積計算装置である。
【００１４】
脳出血部に関するＲＯＩに含まれる骨部を除外することが、正しい結果を得る点で好ましい。脳出血部に関するＲＯＩに含まれる異物部を除外することが、正しい結果を得る点で好ましい。
【発明の効果】
【００１５】
上記各観点での本発明によれば、脳出血を発症している頭部をＸ線ＣＴ装置で撮影して得られる連続スライスの複数の画像について正常な脳実質部に関するＲＯＩと脳出血部に関するＲＯＩをそれぞれ設定し、正常な脳実質部に関するＲＯＩ内の画素のＣＴ値の平均値を複数の画像についてそれぞれ求め、１つの画像の平均値を基準として他の画像の平均値のシフト量を画像ごとに求め、画像ごとのシフト量を用いて脳出血部に関するＲＯＩ内の画素のＣＴ値を画像ごとに調整し、調整後のＣＴ値について複数の画像を通じての最大値を特定し、次式によって脳出血部位の体積を計算するので、
【００１６】
【数５】
2007152105.tif 000005
【００１７】
ここで、
ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての画素番号
Ｎ：脳出血部に関するＲＯＩ内の画素の複数の画像を通じての総画素数
ｐｎ  ：脳出血部に関するＲＯＩ内の画素の調整後のＣＴ値
ＤＦＯＶ：デファインドＦＯＶの大きさ
ＳＦＯＶ：スキャンＦＯＶの大きさ
ＳｌｉｃｅＴｈｉｃｋｎｅｓｓ：スライス厚
Ｍａｇｎｉｆｉｃａｔｉｏｎ：倍率
脳出血部位の体積計算を適正に行う方法および装置を実現することができる。
【００１８】
また、上記各観点での本発明によれば、ユーザーの技量に影響されない安定な結果が得られる脳出血部位計算方法および装置を実現することができる。また、上記各観点での本発明によれば、ユーザーにとって操作しやすい脳出血部位計算方法および装置を実現することができる。
【００１９】
また、上記各観点での本発明によれば、上中下いずれの脳領域にも適応し、パーシャルヴォリューム・エフェクトおよび複数画像間のＣＴ値シフトの問題も克服した脳出血部位計算方法および装置を実現することができる。
【発明を実施するための最良の形態】
【００２０】
以下、図面を参照して発明を実施するための最良の形態を説明する。なお、本発明は、発明を実施するための最良の形態に限定されるものではない。図１に画像処理装置の構成をブロック（ｂｌｏｃｋ）図によって示す。
【００２１】
本装置は本発明を実施するための最良の形態の一例である。本装置の構成によって、脳出血部位体積計算装置に関する発明を実施するための最良の形態の一例が示される。本装置の動作によって、脳出血部位体積計算方法に関する発明を実施するための最良の形態の一例が示される。
【００２２】
図１に示すように、本装置は、データ（ｄａｔａ）処理部１０、表示部２０、操作部３０、記憶部４０および入出力部５０を有する。データ処理部１０は、表示部２０および操作部３０を通じてのユーザー（ｕｓｅｒ）によるインタラクティブ（ｉｎｔｅｒａｃｔｉｖｅ）な操作に基づいて、記憶部４０に記憶されたデータについて所定のデータ処理を行う。データ処理部１０は、また、入出力部５０を通じて外部装置に対するデータ入出力を行う。
【００２３】
脳出血部位体積計算に使用するＸ線ＣＴ画像は、入出力部５０を通じて外部装置から入力される。Ｘ線ＣＴ画像は、連続スライス（ｓｌｉｃｅ）の複数の画像である。以下、Ｘ線ＣＴ画像を単に画像ともいう。
【００２４】
外部装置は例えばＸ線ＣＴ装置や医用画像サーバー（ｓｅｒｖｅｒ）等である。なお、本装置は、Ｘ線ＣＴ装置や医用画像サーバーの一部であってもよい。その場合は、入出力部５０は必ずしも必要でない。
【００２５】
本装置によって行われる脳出血部位体積計算について説明する。図２に、脳出血部位体積計算を行うためのプロセス（ｐｒｏｃｅｓｓ）を示す。図２に示すように、脳出血部位体積計算は３つのプロセスＳ１，Ｓ２，Ｓ３によって行われる。
【００２６】
プロセスＳ１は選択プロセスである。選択プロセスＳ１では、平均値検出部の選択、脳出血部の選択、および、除外部の選択が行われる。なお、除外部の選択は必要に応じて行われる。選択プロセスＳ１の詳細については後にあらためて説明する。
【００２７】
プロセスＳ２は解析／調整プロセスである。解析／調整プロセスＳ２では、ＣＴ値の平均値解析とＣＴ値の調整が行われる。解析／調整プロセスＳ２の詳細については後にあらためて説明する。
【００２８】
プロセスＳ３は計算プロセスである。計算プロセスＳ３では体積計算が行われる。体積計算により、脳出血部位の体積が算出される。計算プロセスＳ３の詳細については後にあらためて説明する。
【００２９】
図３に、本装置の動作のフローチャート（ｆｌｏｗ  ｃｈａｒｔ）を示す。ステップ３０１（ｓｔｅｐ）で、画像呼出を行う。画像呼出は、ユーザーによる操作部３０の操作に基づいて行われる。画像呼出により、複数の画像の最初のものが表示部２０に表示される。
【００３０】
表示画像について、ステップ３０３で、脳出血の有無を判定する。脳出血の有無の判定はユーザーの目視によって行われる。脳出血がないと判定したときは、ステップ３０１に戻って次の画像の呼出を行う。脳出血が見あたらない間は、ステップ３０１，３０３の動作が繰り返される。
【００３１】
脳出血があると判定したときは、ステップ３０５で、平均値検出ＲＯＩ（Ｒｅａｓｉｏｎ  ｏｆ  Ｉｎｔｅｒｅｓｔ）設定を行う。平均値検出ＲＯＩ設定は、画像における正常な脳実質部のＣＴ値の平均値を求めるために行われる。
【００３２】
図４に、平均値検出ＲＯＩ設定の一例を示す。図４に示すように、脳実質部の脳出血を起こしていない部分に、円形のＲＯＩが設定される。ＲＯＩ設定は、例えば、ユーザーが画像の所望の部分を操作部３０のポインティングデバイス（ｐｏｉｎｔｉｎｇ  ｄｅｖｉｃｅ）でクリック（ｃｌｉｃｋ）すること等によって行われる。ＲＯＩの面積は予め定められている。なお、平均値検出ＲＯＩは、前もって、複数の画像に設定しておいてもよい。操作部３０は、本発明における設定手段の一例である。
【００３３】
次に、ステップ３０７で、脳出血部ＲＯＩ設定を行う。脳出血部ＲＯＩ設定は、画像における脳出血部を特定するために行われる。図５に、脳出血部ＲＯＩ設定の一例を示す。図５に示すように、脳出血部を包含するように楕円や長方形等の適宜の図形のＲＯＩが設定される。
【００３４】
ＲＯＩ設定は、例えば、ユーザーが画像の所望の部分を操作部３０のポインティングデバイスでクリックしてＲＯＩ図形を出現させ、その図形の位置、形状、大きさ、傾き等を、脳出血部を包含するように調整することによって行われる。操作部３０は、本発明における設定手段の一例である。
【００３５】
脳出血部ＲＯＩは脳出血部を包含するものであればよく、脳出血部の輪郭を精密になぞる必要はない。また、図５に示すように、同一の脳出血部を含む複数のＲＯＩを設定してよい。このため、ＲＯＩ設定は容易であり、かつ、労力と時間が少なくて済む。したがって、ユーザーにとってやりやすく、上中下いずれの脳領域にも何らの困難もなく適応できる。また、ユーザーの技量にも影響されない。
【００３６】
脳出血部ＲＯＩは骨部を含まないように設定するが、例えば図６に示すように、脳出血部が頭骨に近接しているときはそのような設定は困難なので、骨部を含んで設定せざるを得ない。その場合は、ＣＴ値の相違を利用して脳出血部と骨部の境界を検出し、境界の内側を脳出血部ＲＯＩとする。脳出血部と骨部の境界検出は、データ処理部１０により、既知のアルゴリズム（ａｌｇｏｒｉｔｈｍ）によって自動的に行われる。データ処理部１０は、本発明における骨部除外手段の一例である。
【００３７】
次に、ステップ３０９で、除外部の有無を判定する。除外部の有無の判定はユーザーの目視によって行われる。除外部とは、脳出血でないことが明らかな高輝度部であって、例えば、手術によって脳内に残置された人工物や自然発生の石灰化（Ｃａｌｃｉｆｉｃａｔｉｏｎ）部等である。
【００３８】
除外部があるときは、ステップ３１１で、除外部設定を行う。除外部設定の一例を図７に示す。図７における長方形部が除外部である。除外部設定は、除外すべき物体の輪郭を操作部３０のポインティングデバイス等でなぞることによって行われる。なお、除外部設定は除外物があるときだけの例外的な処理であり、頻度は小さい。操作部３０は、本発明における異物部除外手段の一例である。
【００３９】
特殊な場合として、例えば図８に示すように、血液が脳室内に存在する場合がある。これは、脳実質部の出血ではないので体積計算の対象外である。全自動方式では、このような場合は誤動作が発生するが、本装置では、ユーザーが目視によって判定するので、そのような事態は容易に回避可能である。
【００４０】
ステップ３１３で、設定登録を行う。これによって、平均値検出ＲＯＩ設定と脳出血部ＲＯＩ設定の内容が記憶部４０に記憶される。また、除外部があるときは除外部設定も登録される。これによって、１つの画像について、平均値検出ＲＯＩ設定と脳出血部ＲＯＩ設定（および除外部設定）が確定する。
【００４１】
ステップ３１５で、未呼出画像の有無を判定し、未呼出画像があるときはステップ３０１に戻って次の画像を呼び出し、その画像について、ステップ３０３−３１３の処理を行う。この動作を、未呼出画像がなくなるまで繰り返す。これによって、複数の画像の全てについて、平均値検出ＲＯＩ設定と脳出血部ＲＯＩ設定（および除外部設定）が画像ごとに確定する。ステップ３０３−３１５は、図２に示した選択プロセスＳ１に相当する。
【００４２】
ステップ３１７で、平均値算出を行う。平均値算出は、データ処理部１０によって行われる。データ処理部１０は、本発明における平均値算出手段の一例である。データ処理部１０は、平均値検出ＲＯＩ内のＣＴ値の平均値を画像ごとに求める。これによって、下記のような平均値が得られる。
【００４３】
【数６】
2007152105.tif 000006
【００４４】
ここで、
ｉ：画像番号
Ｉ：最大画像番号
ステップ３１９で、シフト量算出を行う。シフト量算出は、データ処理部１０によって行われる。データ処理部１０は、本発明におけるシフト量算出手段の一例である。データ処理部１０は、次式により、各画像の平均値を基準平均値とシフト量の和に分解する。
【００４５】
【数７】
2007152105.tif 000007
【００４６】
基準平均値としては、例えば１番目の画像の平均値が用いられる。なお、それに限らず、複数の画像のいずれか１つの平均値を基準平均値としてよい。
ステップ３２１で、ＣＴ値調整を行う。ＣＴ値調整は、データ処理部１０によって行われる。データ処理部１０は、本発明における調整手段の一例である。データ処理部１０は、シフト量を用いて脳出血部ＲＯＩの各画素のＣＴ値を調整する。ＣＴ値調整は次式により画像ごとに行われる。
【００４７】
【数８】
2007152105.tif 000008
【００４８】
ここで、
ｉ：画像番号
ｍ：同一画像における画素番号
このようなＣＴ値調整を行うことにより、複数画像間のＣＴ値シフトを調整して、平均値検出ＲＯＩ内の画素値条件を全画像を通じて同一にすることができる。
【００４９】
次に、ステップ３２３で、最大値特定を行う。最大値特定は、データ処理部１０によって行われる。データ処理部１０は、ＣＴ値調整済の画素値の最大値を脳出血部ＲＯＩごとに特定し、それら最大値のうちの最大のものを特定する。これによって、下記のように、ＲＯＩごとの最大値および複数のＲＯＩ全体を通じての最大値がそれぞれ得られる。
【００５０】
【数９】
2007152105.tif 000009
【００５１】
【数１０】
2007152105.tif 000010
【００５２】
ここで、
ｊ：ＲＯＩ番号
Ｊ：最大ＲＯＩ番号
ステップ３１７−３２３は、図２に示した解析／調整プロセスＳ２に相当する。
【００５３】
ステップ３２５で、体積計算を行う。ステップ３２５は、図２に示した計算プロセスＳ３に相当する。体積計算は、データ処理部１０によって行われる。データ処理部１０は、先ず画素の１個当たりの体積を求める。１個当たりの画素体積は次式によって求められる。
【００５４】
【数１１】
2007152105.tif 000011
【００５５】
ここで、
ＤＦＯＶ：デファインド（Ｄｅｆｉｎｅｄ）ＦＯＶの大きさ
ＳＦＯＶ：スキャン（ｓｃａｎ）ＦＯＶの大きさ
ＳｌｉｃｅＴｈｉｃｋｎｅｓｓ：スライス（ｓｌｉｃｅ）厚
Ｍａｇｎｉｆｉｃａｔｉｏｎ：倍率
データ処理部１０は、次に、ＣＴ値調整済の画素値を、複数の画像を通じてＣＴ値順にソート（ｓｏｒｔ）する。これによって、次式に示すような数列が得られる。
【００５６】
【数１２】
2007152105.tif 000012
【００５７】
ここで、
ｐｎ  ：ＣＴ値調整済の画素値
ｎ：複数の画像を通じての画素番号
Ｎ：複数の画像を通じての総画素数
データ処理部１０は、次に、画素値群における最大値を
【００５８】
【数１３】
2007152105.tif 000013
【００５９】
として特定する。データ処理部１０は、本発明における特定手段の一例である。
そして、データ処理部１０は、次式によって体積計算を行う。データ処理部１０は、本発明における計算手段の一例である。
【００６０】
【数１４】
2007152105.tif 000014
【００６１】
上式は、画素体積に重みを付けて全加算する式となっている。重み係数は、画素のＣＴ値と基準平均値との差を、最大値と基準平均値との差で正規化したものとしてある。ただし、基準平均値以下のＣＴ値を持つ画素の重み係数は０とする。
【００６２】
このため、最大のＣＴ値を持つ画素の重み係数は１となり、基準平均値以下のＣＴ値を持つ画素の重み係数０となり、最大値と基準平均値の間のＣＴ値を持つ画素の重み係数は１と０の間の値をとる。
【００６３】
このような重み付けにより、体積計算に寄与する画素はＣＴ値が基準平均値より大きいものだけとなる。そして、その寄与度は、最大のＣＴ値を持つ画素が最も大きく、以下、最大値からの隔たりに応じて寄与度が低くなる。
【００６４】
寄与度の低下は、パーシャルヴォリューム・エフェクトの程度を表す。上式は、パーシャルヴォリューム・エフェクトの大きい画素ほど寄与度を小さくして体積計算を行うものであり、脳出血部位の体積計算用の式として理にかなっている。また、基準平均値以下のＣＴ値を持つ画素の寄与度を０とすることも理にかなっている。
【００６５】
さらに、個々の画素値がノイズ（ｎｏｉｓｅ）を含んでいても、全加算により相殺されるのでノイズの影響も極めて少ない。このため、上式によって計算された体積は、ＣＴ値調整の効果とあいまって脳出血部位の体積を正確に示すものとなる。
【図面の簡単な説明】
【００６６】
【図１】本発明を実施するための最良の形態の一例の画像処理装置の構成を示す図である。
【図２】脳出血部位体積計算を行うためのプロセスを示す図である。
【図３】本発明を実施するための最良の形態の一例の画像処理装置の動作のフローチャートである。
【図４】平均値検出ＲＯＩ設定の一例を中間調の写真で示す図である。
【図５】脳出血部ＲＯＩ設定の一例を中間調の写真で示す図である。
【図６】脳出血部ＲＯＩ設定の一例を中間調の写真で示す図である。
【図７】除外部設定の一例を中間調の写真で示す図である。
【図８】脳の断層像を中間調の写真で示す図である。
【符号の説明】
【００６７】
１０  データ処理部
２０  表示部
３０  操作部
４０  記憶部
５０  入出力部
【図１】
2007152105.tif 000017
【図２】
2007152105.tif 000018
【図３】
2007152105.tif 000019
【図４】
2007152105.tif 000020
【図５】
2007152105.tif 000021
【図６】
2007152105.tif 000022
【図７】
2007152105.tif 000023
【図８】
2007152105.tif 000024