(57)【要約】【課題】経過観察において断層画像群に含まれる病変あるいは病変候補が悪性病変であるか否かを判断するための適切な情報を出力することが可能な医用画像診断支援装置、その作動方法、およびプログラムを提供する。【解決手段】画像取得部４１は今回撮影した断層画像群２１Ｎと前回撮影した断層画像群２１Ｐを取得する。座標取得部４２は断層画像群２１Ｎおよび２１Ｐに含まれる病変２２の３次元的な中心の座標を取得する。画素値分布取得部４３は中心の座標から半径Ｒの球面上の断層画像群２１Ｎおよび２１Ｐの対象画素６０Ａの画素値分布を半径Ｒを変数として複数取得する。画素値分布画像生成部４４は画素値分布を半径Ｒに関して並べた画素値分布画像６３ＡＮおよび６３ＡＰを生成する。入出力制御部４０は画素値分布画像６３ＡＮおよび６３ＡＰを出力する。【選択図】図５
2015221141.tif 000002
【特許請求の範囲】
【請求項１】
被写体を撮影した複数枚の断層画像で構成される断層画像群を取得する画像取得部と、
前記画像取得部で取得した前記断層画像群に含まれる病変あるいは病変候補の３次元的な中心または重心の座標を取得する座標取得部と、
前記座標取得部で取得した前記座標から半径Ｒの球面上の前記断層画像群の対象画素の画素値分布、または前記座標取得部で取得した前記座標から半径Ｒの円周上の前記断層画像の対象画素の画素値分布を、前記半径Ｒを変数として複数取得する画素値分布取得部と、
前記画素値分布取得部で取得した前記画素値分布を前記半径Ｒに関して並べた画素値分布画像を生成する画素値分布画像生成部と、
前記画素値分布画像生成部で生成した前記画素値分布画像を出力する出力制御部とを備えることを特徴とする医用画像診断支援装置。
【請求項２】
前記出力制御部は、撮影時期の異なる複数の前記断層画像群に基づき前記画素値分布画像生成部で生成した複数の前記画素値分布画像を出力する請求項１に記載の医用画像診断支援装置。
【請求項３】
撮影時期の異なる複数の前記断層画像群に基づき前記画素値分布画像生成部で生成した複数の前記画素値分布画像の差分画像を生成する差分画像生成部を備え、
前記出力制御部は、前記差分画像生成部で生成した前記差分画像を出力する請求項１または２に記載の医用画像診断支援装置。
【請求項４】
前記画素値分布取得部は、前記画素値の分布量を度数で表し、かつ前記画素値の範囲を階級で表すヒストグラムを生成するヒストグラム生成部と、
前記ヒストグラム生成部で生成した前記ヒストグラムに基づき前記画素値分布を生成する画素値分布生成部とを有する請求項１ないし３のいずれか１項に記載の医用画像診断支援装置。
【請求項５】
前記画素値分布生成部は、前記画素値の分布量を長さ方向にとり、前記画素値の分布量の大小を長短で表現し、前記画素値の範囲の高低を色の違いで表現した前記画素値分布を、前記ヒストグラムを元に生成する請求項４に記載の医用画像診断支援装置。
【請求項６】
前記画素値分布生成部は、長さが前記半径Ｒの２乗に比例する直線状の前記画素値分布を生成し、
前記画素値分布画像生成部は、前記直線状の前記画素値分布を前記半径Ｒに関して並べ、階段状の前記画素値分布画像を生成する請求項５に記載の医用画像診断支援装置。
【請求項７】
前記画素値分布生成部は、長さが前記半径Ｒの２乗に比例する円弧状の前記画素値分布を生成し、
前記画素値分布画像生成部は、前記円弧状の前記画素値分布を前記半径Ｒに関して並べ、扇形状の前記画素値分布画像を生成する請求項５に記載の医用画像診断支援装置。
【請求項８】
前記画素値分布生成部は、長さを正規化して揃えた直線状の前記画素値分布を生成し、
前記画素値分布画像生成部は、前記直線状の前記画素値分布を前記半径Ｒに関して並べ、矩形状の前記画素値分布画像を生成する請求項５に記載の医用画像診断支援装置。
【請求項９】
前記画素値分布生成部は、前記画素値の範囲を長さ方向にとり、前記画素値の分布量の大小を色の違いで表現した前記画素値分布を、前記ヒストグラムを元に生成する請求項４に記載の医用画像診断支援装置。
【請求項１０】
前記画素値の分布量は、前記画素値の範囲毎の前記対象画素の個数、
前記画素値の範囲毎の前記対象画素の個数を全対象画素の個数で除算した、前記画素値の範囲毎の前記対象画素の個数の割合、
前記球面を複数に分割した微小領域を前記対象画素とした場合の前記画素値の範囲毎の前記対象画素の面積、
または前記球面を複数に分割した微小領域を前記対象画素とした場合の前記画素値の範囲毎の前記対象画素の面積を、前記球面の面積で除算した、前記対象画素の面積の割合のいずれかである請求項４ないし９のいずれか１項に記載の医用画像診断支援装置。
【請求項１１】
画像取得部により、被写体を撮影した複数枚の断層画像で構成される断層画像群を取得する画像取得ステップと、
座標取得部により、前記画像取得ステップで取得した前記断層画像群に含まれる病変あるいは病変候補の３次元的な中心または重心の座標を取得する座標取得ステップと、
画素値分布取得部により、前記座標取得ステップで取得した前記座標から半径Ｒの球面上の前記断層画像群の対象画素の画素値分布、または前記座標取得ステップで取得した前記座標から半径Ｒの円周上の前記断層画像の対象画素の画素値分布を、前記半径Ｒを変数として複数取得する画素値分布取得ステップと、
画素値分布画像生成部により、前記画素値分布取得ステップで取得した前記画素値分布を前記半径Ｒに関して並べた画素値分布画像を生成する画素値分布画像生成ステップと、
出力制御部により、前記画素値分布画像生成ステップで生成した前記画素値分布画像を出力する出力制御ステップとを備えることを特徴とする医用画像診断支援装置の作動方法。
【請求項１２】
被写体を撮影した複数枚の断層画像で構成される断層画像群を取得する画像取得機能と、
前記画像取得機能で取得した前記断層画像群に含まれる病変あるいは病変候補の３次元的な中心または重心の座標を取得する座標取得機能と、
前記座標取得機能で取得した前記座標から半径Ｒの球面上の前記断層画像群の対象画素の画素値分布、または前記座標取得機能で取得した前記座標から半径Ｒの円周上の前記断層画像の対象画素の画素値分布を、前記半径Ｒを変数として複数取得する画素値分布取得機能と、
前記画素値分布取得機能で取得した前記画素値分布を前記半径Ｒに関して並べた画素値分布画像を生成する画素値分布画像生成機能と、
前記画素値分布画像生成機能で生成した前記画素値分布画像を出力する出力制御機能とを、コンピュータに実行させることを特徴とする医用画像診断支援プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医用画像診断支援装置、医用画像診断支援装置の作動方法、および医用画像診断支援プログラムに関する。
【背景技術】
【０００２】
従来、医療分野において、ＣＴ（Computed tomography）装置やＭＲＩ（Magnetic Resonance Imaging）装置等のモダリティで被写体（患者）を撮影した複数枚の断層画像で構成される断層画像群に基づく画像診断が行われている。こうした画像診断において、一定期間を開けて撮影を行い、これにより得られた撮影時期の異なる複数の断層画像群を比較読影することで、例えば、断層画像群に含まれる病変が癌等の悪性病変であるか否かを判断する経過観察が行われる。病変が悪性病変の疑いありと判断した場合は精密検査や治療等の適切な処置が講じられる。
【０００３】
特許文献１には、経過観察における病変が悪性病変であるか否かの判断をアシストするため、病変の成長度を算出し、算出した成長度を表示する医用画像診断支援装置が記載されている。この医用画像診断支援装置は、まず、気管支分岐部のスライス位置情報等の生体臓器の構造物の位置に基づき、撮影時期の異なる複数の断層画像群、例えば今回撮影した断層画像群と前回撮影した断層画像群の３次元的な位置合わせをする。次いで、今回の断層画像群から病変を抽出し、抽出した病変と、抽出した病変と同一位置に存在する前回の断層画像群の病変とを解析して、それぞれの病変の特徴量を算出する。そして、算出した特徴量に基づき病変の成長度を算出している。特徴量の例として、病変のうちの高い画素値の領域と低い画素値の領域のランレングスの割合である粗雑度や、今回の断層画像群の病変中心部の画素のうち、前回の断層画像群の病変の画素の画素値から変化した個数、病変の長径、短径、体積等が挙げられている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００８−１７３２１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
特許文献１では、今回と前回の断層画像群の３次元的な位置合わせをした後、この３次元的な位置合わせが正確に行われていることを前提に、今回の断層画像群から抽出した病変と、抽出した病変と同一位置に存在する前回の断層画像群の病変とを解析して、それぞれの病変の特徴量を算出している。このため、もし３次元的な位置合わせが正確に行われていないと、特徴量、ひいては特徴量に基づき算出する成長度の信頼性が低くなってしまう。また、今回の断層画像群からの病変の抽出が正確に行われない場合も、成長度の信頼性が低くなってしまう。
【０００６】
ここで、撮影時期の異なる複数の断層画像群に含まれる病変は、その病変がたとえ全く経時変化がないものであっても、撮影時の被写体の位置、被写体の向きの変化や、被写体の呼吸、被写体の心拍等の周辺臓器の動きの変化等によって、病変の位置および病変の向きが違ってくる。被写体内における病変の位置は、例えば病変の中心や重心等の１点の基準点を決定し、この基準点に基づき撮影時期の異なる複数の断層画像群の間で位置合わせを行えば、ある程度位置合わせの精度の低下による影響を抑制することができるが、病変の向きは、断層画像の解像度（スライス厚、画素数等）や病変自体の多少の撓みも考慮すると、完全に一致させることは難しい。ここで病変の向きとは、２次元の断層画像においては、病変の中心を通り、かつ２次元平面に垂直な軸を回転軸とする２次元の回転方向、３次元の断層画像群においては、病変の中心を回転中心とする３次元の回転方向である。特に形状および内部構造が均一でない病変程、向きを一致させることができないことによる成長度の信頼性低下への影響が大きくなる。また、病変は、隣接する正常領域と峻別可能な明瞭な境界があるとは限らず、この場合の病変を正確に抽出することの難しさが、向きを合わせることの精度にも影響する。
【０００７】
このように、特許文献１に記載の成長度は、位置合わせの精度および病変の抽出精度によって信頼性が左右される。したがって、特許文献１に記載の成長度は、経過観察において断層画像群に含まれる病変が悪性病変であるか否かを判断するための情報としては不適切である。
【０００８】
さらに、病変には様々な大きさがあり、大きさが１ｃｍ程度と比較的小さいものもある。こうした比較的小さな病変は、早期癌の初期状態である可能性があり、重要な経過観察対象の１つであるため、病変の大きさや濃度（画素値の高低）の微妙な経時変化を正しく認識して、その病変が悪性病変であるか否かを判断する必要がある。しかしながら、比較的小さな病変は、比較的大きな病変と比べて画像の情報量が少ないため、位置合わせ、および病変の抽出を正確に行うことがより難しくなる。このため、比較的小さな病変に対しては、正確な位置合わせを前提とした成長度の信頼性はより低くなり、比較的小さな病変が悪性病変であるか否かを成長度で判断することはリスクがある。また、成長度という定量的な数値だけでは、病変の大きさや濃度の微妙な経時変化は視覚的には分からないため、定性的な経時変化を正しく認識することができない。
【０００９】
本発明は、上記課題に鑑みてなされたもので、位置合わせの精度および病変の抽出精度による信頼性への影響が小さく、経過観察において断層画像群に含まれる病変あるいは病変候補が悪性病変であるか否かを判断するための適切な情報を出力することが可能な医用画像診断支援装置、医用画像診断支援装置の作動方法、および医用画像診断支援プログラムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
上記目的を達成するために、本発明の医用画像診断支援装置は、被写体を撮影した複数枚の断層画像で構成される断層画像群を取得する画像取得部と、画像取得部で取得した断層画像群に含まれる病変あるいは病変候補の３次元的な中心または重心の座標を取得する座標取得部と、座標取得部で取得した座標から半径Ｒの球面上の断層画像群の対象画素の画素値分布、または座標取得部で取得した座標から半径Ｒの円周上の断層画像の対象画素の画素値分布を、半径Ｒを変数として複数取得する画素値分布取得部と、画素値分布取得部で取得した画素値分布を半径Ｒに関して並べた画素値分布画像を生成する画素値分布画像生成部と、画素値分布画像生成部で生成した画素値分布画像を出力する出力制御部とを備えている。
【００１１】
出力制御部は、撮影時期の異なる複数の断層画像群に基づき画素値分布画像生成部で生成した複数の画素値分布画像を出力することが好ましい。
【００１２】
撮影時期の異なる複数の断層画像群に基づき画素値分布画像生成部で生成した複数の画素値分布画像の差分画像を生成する差分画像生成部を備え、出力制御部は、差分画像生成部で生成した差分画像を出力することが好ましい。
【００１３】
画素値分布取得部は、画素値の分布量を度数で表し、かつ画素値の範囲を階級で表すヒストグラムを生成するヒストグラム生成部と、ヒストグラム生成部で生成したヒストグラムに基づき画素値分布を生成する画素値分布生成部とを有することが好ましい。
【００１４】
画素値分布生成部は、例えば、画素値の分布量を長さ方向にとり、画素値の分布量の大小を長短で表現し、画素値の範囲の高低を色の違いで表現した画素値分布を、ヒストグラムを元に生成する。
【００１５】
画素値分布生成部は、例えば、長さが半径Ｒの２乗に比例する直線状の画素値分布を生成し、画素値分布画像生成部は、直線状の画素値分布を半径Ｒに関して並べ、階段状の画素値分布画像を生成する。また、画素値分布生成部は、長さが半径Ｒの２乗に比例する円弧状の画素値分布を生成し、画素値分布画像生成部は、円弧状の画素値分布を半径Ｒに関して並べ、扇形状の画素値分布画像を生成してもよい。あるいは、画素値分布生成部は、長さを正規化して揃えた直線状の画素値分布を生成し、画素値分布画像生成部は、直線状の画素値分布を半径Ｒに関して並べ、矩形状の画素値分布画像を生成してもよい。
【００１６】
画素値分布生成部は、画素値の範囲を長さ方向にとり、画素値の分布量の大小を色の違いで表現した画素値分布を、ヒストグラムを元に生成してもよい。
【００１７】
画素値の分布量は、例えば、画素値の範囲毎の対象画素の個数、画素値の範囲毎の対象画素の個数を全対象画素の個数で除算した、画素値の範囲毎の対象画素の個数の割合、球面を複数に分割した微小領域を対象画素とした場合の画素値の範囲毎の対象画素の面積、または球面を複数に分割した微小領域を対象画素とした場合の画素値の範囲毎の対象画素の面積を、球面の面積で除算した、対象画素の面積の割合のいずれかである。
【００１８】
また、本発明の医用画像診断支援装置の作動方法は、画像取得部により、被写体を撮影した複数枚の断層画像で構成される断層画像群を取得する画像取得ステップと、座標取得部により、画像取得ステップで取得した断層画像群に含まれる病変あるいは病変候補の３次元的な中心または重心の座標を取得する座標取得ステップと、画素値分布取得部により、座標取得ステップで取得した座標から半径Ｒの球面上の断層画像群の対象画素の画素値分布、または座標取得ステップで取得した座標から半径Ｒの円周上の断層画像の対象画素の画素値分布を、半径Ｒを変数として複数取得する画素値分布取得ステップと、画素値分布画像生成部により、画素値分布取得ステップで取得した画素値分布を半径Ｒに関して並べた画素値分布画像を生成する画素値分布画像生成ステップと、出力制御部により、画素値分布画像生成ステップで生成した画素値分布画像を出力する出力制御ステップとを備えている。
【００１９】
さらに、本発明の医用画像診断支援プログラムは、被写体を撮影した複数枚の断層画像で構成される断層画像群を取得する画像取得機能と、画像取得機能で取得した断層画像群に含まれる病変あるいは病変候補の３次元的な中心または重心の座標を取得する座標取得機能と、座標取得機能で取得した座標から半径Ｒの球面上の断層画像群の対象画素の画素値分布、または座標取得機能で取得した座標から半径Ｒの円周上の断層画像の対象画素の画素値分布を、半径Ｒを変数として複数取得する画素値分布取得機能と、画素値分布取得機能で取得した画素値分布を半径Ｒに関して並べた画素値分布画像を生成する画素値分布画像生成機能と、画素値分布画像生成機能で生成した画素値分布画像を出力する出力制御機能とを、コンピュータに実行させる。
【発明の効果】
【００２０】
本発明によれば、断層画像群に含まれる病変あるいは病変候補の３次元的な中心または重心の座標から半径Ｒの球面上の断層画像群の対象画素の画素値分布、または断層画像群に含まれる病変あるいは病変候補の３次元的な中心または重心の座標から半径Ｒの円周上の断層画像の対象画素の画素値分布を半径Ｒに関して並べた画素値分布画像を生成し、生成した画素値分布画像を出力するので、病変の３次元的な中心または重心の座標の１点の基準点による位置合わせだけで済む。また、病変の向きの違いは、半径Ｒの球面上の断層画像群の対象画素の画素値分布、または半径Ｒの円周上の断層画像の対象画素の画素値分布を取得することで相殺することができ、病変の向きを正確に合わせる必要がない。さらに、画素値分布画像で画素値分布を２次元で表現することによって、病変の大きさや濃度の微妙な経時変化を正しく認識することができる。したがって、位置合わせの精度および病変の抽出精度による信頼性への影響が小さく、経過観察において断層画像群に含まれる病変あるいは病変候補が悪性病変であるか否かを判断するための適切な情報を出力することが可能な医用画像診断支援装置、医用画像診断支援装置の作動方法、および医用画像診断支援プログラムを提供することができる。
【図面の簡単な説明】
【００２１】
【図１】医療情報ネットワークシステムを示す図である。
【図２】断層画像と断層画像群を示す説明図である。
【図３】画像ファイルの内容を示す図である。
【図４】支援装置を構成するコンピュータを示すブロック図である。
【図５】支援装置のＣＰＵの機能を示すブロック図である。
【図６】病変抽出部による病変の抽出、および座標算出部による病変の中心または重心の座標の算出の様子を示す説明図である。
【図７】病変の中心から半径Ｒ（Ｋ）の球面ＳＰ（Ｋ）上の対象画素の画素値分布を取得する様子を示す説明図である。
【図８】病変の中心から半径Ｒ（Ｋ）の球面ＳＰ（Ｋ）上の対象画素の画素値分布を取得する様子を示す説明図である。
【図９Ａ】病変の中心から半径Ｒ（Ｋ）の球面ＳＰ（Ｋ）上の対象画素を示す説明図である。
【図９Ｂ】図９Ａの平面図である。
【図１０Ａ】病変の中心から半径Ｒ（Ｋ）の球面ＳＰ（Ｋ）上の対象画素を示す説明図である。
【図１０Ｂ】図１０Ａの平面図である。
【図１１】ヒストグラムを示す図である。
【図１２】ヒストグラムを元に画素値分布を生成する様子を示す説明図である。
【図１３】画素値分布画像を示す図である。
【図１４】今回と前回撮影した断層画像を表示する比較読影ウィンドウを示す図である。
【図１５】断層画像と画素値分布画像を表示する比較読影ウィンドウを示す図である。
【図１６】画素値分布画像を拡大表示した比較読影ウィンドウを示す図である。
【図１７】支援装置と診療科端末、医用画像ＤＢサーバ間の各種情報の遣り取りを示す説明図である。
【図１８】支援装置の作動手順を示すフローチャートである。
【図１９】差分画像生成部を設けた支援装置のＣＰＵの機能を示すブロック図である。
【図２０】差分画像を示す図である。
【図２１】直線状の画素値分布から円弧状の画素値分布を生成し、円弧状の画素値分布から扇形状の画素値分布画像を生成する様子を示す説明図である。
【図２２】長さを正規化して揃えた直線状の画素値分布から生成した矩形状の画素値分布画像を示す図である。
【図２３】画素値の範囲を長さ方向にとり、画素値の分布量の大小を色の違いで表現した画素値分布を、ヒストグラムを元に生成する様子を示す説明図である。
【図２４】図２３の画素値分布から生成した矩形状の画素値分布画像を示す図である。
【図２５】球面ＳＰ（Ｋ）を複数の微小領域に分割し、微小領域を球面上の対象画素とする様子を示す説明図である。
【発明を実施するための形態】
【００２２】
[第１実施形態]
図１において、医療情報ネットワークシステム１０は、病院などの医療施設に構築され、ＣＴ装置１１、医用画像データベース（以下、ＤＢと略す）サーバ１２、診療科端末１３、医用画像診断支援装置（以下、単に支援装置という）１４、およびこれらを相互通信可能に接続するネットワーク１５等で構成される。ネットワーク１５は、例えば、医療施設内に敷設されたＬＡＮ（Local Area Network）や移動体通信ネットワークである。
【００２３】
医療情報ネットワークシステム１０は、上記の他にも、病院情報システム（ＨＩＳ；Hospital Information System）や放射線情報システム（ＲＩＳ；Radiology Information System）を備えている。ＨＩＳは、電子カルテ、会計情報、検査予約情報、投薬処方情報等、医療施設の医療行為に関わる各種情報を、診療科、放射線科等の医療施設の各部門から受け付け、各種情報を保管および管理する。検査予約情報は、診療科の医師（診断医）により診療科端末１３の電子カルテ上で発行される。検査予約情報は、血液検査、内視鏡検査等の種々の検査予約の他、ＣＴ装置１１による断層画像撮影の検査予約である撮影オーダを含む。撮影オーダは、断層画像撮影を医師や放射線技師等のＣＴ装置１１の使用者に指示するためのものである。
【００２４】
ＨＩＳには、被写体情報が記憶されている。被写体情報は電子カルテに登録されている。被写体情報は、例えば、個々の被写体を識別するための被写体ＩＤ（IDentification）、被写体氏名、性別、生年月日、年齢、身長、体重等の項目を有する（図３参照）。
【００２５】
ＲＩＳは、診療科端末１３からの撮影オーダを受け付けて、撮影オーダを保管および管理する。撮影オーダは、例えば、個々のオーダを識別するためのオーダＩＤ、当該撮影オーダを発行した診断医の医師ＩＤ、当該撮影オーダによる断層画像撮影の対象の被写体ＩＤ、当該撮影オーダによる断層画像撮影を担当する使用者ＩＤ、撮影部位および／または向き等の項目を有する（図３参照）。撮影部位には、頭部、胸部、腹部といった人体の各部位がある。また向きには仰向け、うつ伏せ、横臥等がある。撮影オーダには、これらの項目の他に、ＲＩＳで撮影オーダを受け付けた日時や、経過観察等の断層画像撮影の目的、診断医から使用者への申し渡し事項等の項目が設けられている。
【００２６】
ＲＩＳは、撮影オーダをＣＴ装置１１に送信する。使用者は、ＣＴ装置１１のコンソールのディスプレイで撮影オーダを確認し、確認した撮影オーダに応じた撮影条件をＣＴ装置１１に設定して断層画像撮影を行う。撮影条件は、放射線、例えばＸ線を被写体に向けて照射するＸ線源に与える管電圧、管電流、Ｘ線の照射時間、各断層画像の撮影間隔であるスライス厚等の項目を有する（図３参照）。
【００２７】
診療科端末１３は診断医によって操作される。診療科端末１３は、ディスプレイ１６と入力デバイス１７を有する。ディスプレイ１６は、マウスやキーボード等の入力デバイス１７の操作に応じた各種操作画面を表示する。操作画面にはＧＵＩ(Graphical User Interface)が配され、診療科端末１３はこのＧＵＩを通じて入力デバイス１７からの操作指示の入力を受け付ける。
【００２８】
ＣＴ装置１１は、設定された撮影条件に応じて断層画像撮影を行い、図２に示すように複数枚の断層画像２０で構成される断層画像群２１を得る。図２では、撮影部位として胸部が指定され、設定されたスライス厚で、断層画像２０で表されるＸＹ平面に垂直なＺ方向に順次スキャンして断層画像撮影を行った結果得られる複数枚の断層画像２０を示している。また、図２では、左の肺野下部に病変２２が存在し、この病変２２が、連続する３枚の断層画像２０に映り込んでいる状態を示している。ＣＴ装置１１は、得られた断層画像２０の画像ファイル２３を作成し、作成した画像ファイル２３を医用画像ＤＢサーバ１２に送信する。医用画像ＤＢサーバ１２は、いわゆるＰＡＣＳ（Picture Archiving and Communication System）サーバであり、ＣＴ装置１１から送信された画像ファイル２３を受信して、画像ファイル２３を保管および管理する。
【００２９】
図３において、画像ファイル２３は、例えばＤＩＣＯＭ（Digital Imaging and Communication in Medicine）規格に準拠したファイル形式であり、タグ記憶領域２５と画像記憶領域２６とを有する。画像記憶領域２６には断層画像２０が記憶される。タグ記憶領域２５には、被写体ＩＤ、被写体氏名等の被写体情報、オーダＩＤ、撮影部位および／または向き等の撮影オーダ、および管電圧（単位ｋＶ）、管電流（単位ｍＡ）、Ｘ線の照射時間（単位ｍｓｅｃ）、スライス厚（単位ｍｍ）等の撮影条件の各項目が設けられている。
【００３０】
また、タグ記憶領域２５には、ファイルＩＤ、撮影日時、および病変ＩＤの項目が設けられている。ファイルＩＤは、個々の画像ファイル２３を識別するための番号であり、画像ファイル２３の作成時にＣＴ装置１１により自動的に付される。図２で説明したように、断層画像撮影においては、１件の撮影オーダに対して、複数枚の断層画像２０で構成される断層画像群２１が得られるため、ファイルＩＤは、「Ｆ２１６」等の撮影オーダ毎に割り振られる符号と、「１」等の複数枚の断層画像２０のうちの何枚目かを表す符号をハイフンで繋げた構成である。撮影オーダ毎に割り振られる符号が同じ画像ファイル２３の断層画像２０が、１件の撮影オーダで得られる断層画像群２１を構成する。
【００３１】
病変ＩＤは、個々の病変２２を識別するための番号であり、画像記憶領域２６に記憶された断層画像２０から病変２２が抽出された場合に付される。
【００３２】
医用画像ＤＢサーバ１２は、付帯情報の各項目を検索キーとする外部装置からの断層画像群２１の取得要求に応じて、検索キーに該当する断層画像群２１を検索し、検索した断層画像群２１を取得要求元の外部装置に送信する。取得要求は、検索キーの他に取得要求元の外部装置の識別情報（ＩＰ（Internet Protocol）アドレス）等を含む。医用画像ＤＢサーバ１２は、この識別情報を元に断層画像群２１の送信先を特定する。医用画像ＤＢサーバ１２は、断層画像群２１だけでなく、特定の断層画像２０の取得要求に応じて、特定の断層画像２０を検索して送信することも可能である。なお、断層画像２０や断層画像群２１を外部装置に送信すると表現しているが、実際には画像ファイル２３を外部装置に送信する。
【００３３】
医用画像ＤＢサーバ１２は、同一の被写体に対して一定期間を開けて撮影を行って得られた撮影時期の異なる複数の断層画像群、例えば、今回撮影した断層画像群２１Ｎと前回撮影した断層画像群２１Ｐの取得要求５６（図５参照）を支援装置１４から受けて、断層画像群２１Ｎおよび２１Ｐを検索し、検索した断層画像群２１Ｎおよび２１Ｐを支援装置１４に送信する。なお、以下の説明では、今回の撮影に関する画像に添字の「Ｎ」を、前回の撮影に関する画像に添字の「Ｐ」を付して区別する。
【００３４】
取得要求に含まれる検索キーとしては、被写体ＩＤや撮影日時、オーダＩＤ、使用者ＩＤ等、画像ファイル２３のタグ記憶領域２５に記憶された任意の項目を指定することが可能である。また、例えば撮影日時が６か月前の全ての断層画像群２１を指定する等、検索範囲を指定することも可能であり、被写体ＩＤがＰ０５００で、かつ使用者ＩＤがＲ０００１の断層画像群２１等、複数の項目を用いていわゆるアンド検索やオア検索をすることも可能である。
【００３５】
診療科端末１３は、断層画像群２１の取得要求を医用画像ＤＢサーバ１２に送信し、この取得要求に応じて医用画像ＤＢサーバ１２から送信された断層画像群２１を受信する。また、診療科端末１３は、断層画像群２１の解析要求５５（図５参照）を支援装置１４に送信し、この解析要求５５に応じて支援装置１４から送信された画素値分布画像６３Ａ（図１３参照）を受信する。断層画像群２１の解析要求５５には、画素値分布画像６３Ａを生成する元となる断層画像群２１を指定するための検索キーや、解析要求５５を発した診療科端末１３の識別情報（ＩＰ（Internet Protocol）アドレス）等が含まれている。
【００３６】
ディスプレイ１６は、医用画像ＤＢサーバ１２からの断層画像２０や支援装置１４からの画素値分布画像６３Ａ等を表示する（図１４、図１５、および図１６参照）。診療科端末１３は、ディスプレイ１６に表示された断層画像２０に関する診断医の診断結果の入力を、入力デバイス１７を介して受け付ける。診断結果には、断層画像群２１Ｎと断層画像群２１Ｐを比較読影し、病変２２が癌等の悪性病変であるか否かを判断した経過観察の結果が含まれる。病変が悪性病変の疑いありの判断結果の場合は、精密検査や治療等の適切な処置が講じられる。
【００３７】
支援装置１４は、診療科端末１３で診断医が経過観察を行う際に、病変２２が癌等の悪性病変であるか否かを判断するための適切な情報として、断層画像群２１Ｎおよび２１Ｐに基づき生成した画素値分布画像６３ＡＮおよび６３ＡＰを診療科端末１３に出力する。
【００３８】
図４において、支援装置１４を構成するコンピュータは、ストレージデバイス３０、メモリ３１、ＣＰＵ（Central Processing Unit）３２、および通信部３３を備えている。これらはデータバス３４を介して相互接続されている。
【００３９】
ストレージデバイス３０は、支援装置１４を構成するコンピュータに内蔵、またはケーブルやネットワークを通じて接続されたハードディスクドライブである。ストレージデバイス３０には、オペレーティングシステム等の制御プログラムや、医用画像診断支援プログラム（以下、単に支援プログラムという）３５を含む各種アプリケーションプログラムが記憶されている。支援プログラム３５は、コンピュータを、支援装置１４として機能させるためのプログラムである。
【００４０】
メモリ３１は、ＣＰＵ３２が処理を実行するためのワークメモリである。ＣＰＵ３２は、ストレージデバイス３０に記憶されたプログラムをメモリ３１へロードして、プログラムにしたがった処理を実行することにより、コンピュータの各部を統括的に制御する。
【００４１】
通信部３３は、ネットワーク１５を介した、外部装置との各種情報の伝送を行うネットワークインターフェースである。通信部３３は、診療科端末１３から断層画像群２１Ｎおよび２１Ｐの解析要求５５を、医用画像ＤＢサーバ１２から取得要求５６に応じた断層画像群２１Ｎおよび２１Ｐをそれぞれ受信し、かつ医用画像ＤＢサーバ１２に取得要求５６を、診療科端末１３に画素値分布画像６３ＡＮおよび６３ＡＰをそれぞれ送信する。
【００４２】
図５において、ＣＰＵ３２は、支援プログラム３５を起動すると、メモリ３１と協働して、入出力制御部（出力制御部に相当）４０、画像取得部４１、座標取得部４２、画素値分布取得部４３、および画素値分布画像生成部４４として機能する。座標取得部４２は、病変抽出部５０と座標算出部５１とを有する。また、画素値分布取得部４３は、ヒストグラム生成部５２と画素値分布生成部５３とを有する。
【００４３】
入出力制御部４０は、通信部３３で受信した各種情報を受け付ける入力制御機能、および通信部３３に各種情報を出力する出力制御機能を担う。具体的には、入出力制御部４０は、診療科端末１３からの断層画像群２１Ｎおよび２１Ｐの解析要求５５、および取得要求５６に応じて医用画像ＤＢサーバ１２から送信されて通信部３３で受信した断層画像群２１Ｎおよび２１Ｐを受け付ける。また、入出力制御部４０は、受け付けた解析要求５５に対応する取得要求５６を通信部３３に出力し、通信部３３を介して取得要求５６を医用画像ＤＢサーバ１２に送信させる。さらに、入出力制御部４０は、画素値分布画像６３ＡＮおよび６３ＡＰを通信部３３に出力し、通信部３３を介して解析要求元の診療科端末１３に送信する。
【００４４】
画像取得部４１は、入出力制御部４０で受け付けた医用画像ＤＢサーバ１２からの断層画像群２１Ｎおよび２１Ｐを取得する画像取得機能を担う。画像取得部４１は、取得した断層画像群２１Ｎおよび２１Ｐを座標取得部４２および画素値分布取得部４３に出力する。
【００４５】
座標取得部４２は、画像取得部４１で取得した断層画像群２１Ｎおよび２１Ｐに含まれる病変２２の３次元的な中心の座標を取得する座標取得機能を担う。座標取得部４２の病変抽出部５０は、断層画像群２１Ｎおよび２１Ｐを構成する複数枚の断層画像２０Ｎおよび２０Ｐの全てに対して病変２２の抽出を行う。病変抽出部５０は、例えば、断層画像２０Ｎおよび２０Ｐの画素の画素値（ＣＴ値）と、病変２２とその他の領域を識別するための予め設定された閾値との高低を比較し、閾値以上の画素値の画素の領域を病変２２として抽出する。病変抽出部５０は、病変２２の抽出結果を座標算出部５１に出力する。
【００４６】
座標算出部５１は、病変抽出部５０からの病変２２の抽出結果に基づき、病変２２の３次元的な中心の座標を算出する。座標算出部５１は、断層画像２０Ｎおよび２０Ｐのそれぞれに基づき算出した座標が一致または類似する病変２２を同一の病変２２として認識し、当該断層画像２０Ｎおよび２０Ｐのそれぞれの画像ファイル２３の病変ＩＤの項目に同じ番号を付す。
【００４７】
病変２２の中心とは、文字通り断層画像２０Ｎおよび２０Ｐにおいて病変２２として抽出した画素の領域の幾何学的な中心であり、例えば断層画像２０Ｎおよび２０Ｐで病変２２として抽出した画素の領域が円形の場合はその円の中心である。
【００４８】
連続する複数枚の断層画像２０に同一の病変２２が映り込んでいる場合は、座標算出部５１は、連続する複数枚の断層画像２０から抽出した同一の病変２２のうち、大きさが最も大きいものの中心の座標を算出する。座標算出部５１は、連続する複数枚の断層画像２０のうち、病変２２が最も大きく中心の座標を算出した断層画像２０の画像ファイル２３に病変ＩＤを付す。座標取得部４２は、座標算出部５１で算出した座標を画素値分布取得部４３に出力する。
【００４９】
図６は、病変抽出部５０による病変２２の抽出、および座標算出部５１による病変２２の中心の座標の算出の様子を模式的に示したものである。断層画像２０の画素６０の画素値の高低をハッチングの濃淡で示しており、ハッチングの濃さが薄い程、画素値が高いことを示している。下段の実線で囲む領域５７が、画素値が閾値以上の画素６０の領域であり、病変抽出部５０により病変２２として抽出された領域である。この領域５７は平面形状が円形ないしは正方形であり、領域５７内での画素値の高い部分が右上側に偏在している。符号Ｃ１は病変２２の中心である。中心Ｃ１の座標は、中心Ｃ１上の画素６０の断層画像群２１における位置を示す（Ｘ、Ｙ、Ｚ）の３次元座標で表される。
【００５０】
画素値分布取得部４３は、座標取得部４２で取得した病変２２の中心Ｃ１の座標から半径Ｒの球面ＳＰ上の断層画像群２１の画素６０（以下、対象画素６０Ａと表現する）の画素値分布６２Ａ（図１２参照）を、半径Ｒを変数として複数取得する画素値分布取得機能を担う。
【００５１】
図７、図８は、中心Ｃ１の座標から半径Ｒ（Ｋ）の球面ＳＰ（Ｋ）を示す図であり、図７はＸＹの２次元、図８はスキャン方向のＺ方向を含むＸＹＺの３次元で表現したものである。Ｋは、半径Ｒ、球面ＳＰの番号を示す０以上の自然数である（Ｋ＝０、１、２、３、・・・、Ｎ−２、Ｎ−１、Ｎ）。半径Ｒ（０）の球面ＳＰ（０）は中心Ｃ１であり、この場合の対象画素６０Ａは中心Ｃ１上の１個の画素６０である。半径Ｒ（Ｋ）と半径Ｒ（Ｋ＋１）の差分ΔＲは、例えば画素６０の大きさ以下の値である。半径Ｒ（Ｎ）は、半径Ｒの最大値であり、病変２２とその周辺近傍の領域を包含する大きさである。半径Ｒ（Ｎ）は、例えば数十〜数百個分の画素６０の大きさである。半径Ｒ（Ｎ）は、早期癌の初期状態である可能性がある大きさが１ｃｍ程度の病変２２を対象とする場合は、メートル換算すると約１．５ｃｍ〜２．０ｃｍに設定される。
【００５２】
図９Ａ、図１０Ａにおいてハッチングで示すように、本実施形態では、断層画像群２１の画素６０のうち、球面ＳＰ（Ｋ）内に画素６０の中心が含まれ、かつ球面ＳＰ（Ｋ）の球面上の位置にある断層画像群２１の画素６０を対象画素６０Ａとする。図９Ａでは、中心Ｃ１を有する断層画像２０Ａの対象画素６０Ａを示し、図１０Ａでは、断層画像２０Ａに前後し、かつ球面ＳＰ（Ｋ）と交わる断層画像２０Ｂおよび２０Ｃの対象画素６０Ａを示している。図９Ｂ、図１０Ｂは、それぞれ図９Ａ、図１０Ａの平面図である。
【００５３】
図９Ａにおいて、断層画像２０Ａの対象画素６０Ａは計２４個存在する。また、図１０Ａにおいて、断層画像２０Ｂおよび２０Ｃの対象画素６０Ａはそれぞれ１６個存在する。つまり、この３枚の断層画像２０Ａ、２０Ｂ、および２０Ｃの場合の対象画素６０Ａは計５６個存在することになる。対象画素６０Ａの個数は、半径Ｒ（０）の場合が１個と最も少なく、半径Ｒ（Ｋ）が大きくなる程多くなる。
【００５４】
画素値分布取得部４３のヒストグラム生成部５２は、対象画素６０Ａに関して、その画素値が−１０００〜−９０１、−９００〜−８０１、・・・、−１００〜−１、０〜９９、１００〜１９９等の画素値の範囲のいずれにあたるかを判定し、画素値の範囲毎の対象画素６０Ａの個数を計数する。本実施形態では、この対象画素６０Ａの個数を画素値の分布量として用いる。
【００５５】
図１１に示すように、ヒストグラム生成部５２は、画素値の範囲毎に計数した対象画素６０Ａの個数を縦軸の度数で表し、かつ画素値の範囲を横軸の階級で表すヒストグラム６１を半径Ｒ（Ｋ）毎に生成する。病変２２は中心Ｃ１に近付く程、画素値が高くなる傾向があるので、半径Ｒ（Ｋ）が小さい程、ヒストグラム６１は画素値の範囲が高い横軸の右側の度数が多いものとなる。反対に半径Ｒ（Ｋ）が大きい程、ヒストグラム６１は画素値の範囲が低い横軸の左側の度数が多いものとなる。ヒストグラム生成部５２は、生成したヒストグラム６１を画素値分布生成部５３に出力する。
【００５６】
図１２に示すように、画素値分布生成部５３は、ヒストグラム生成部５２で生成したヒストグラム６１を元に、対象画素６０Ａの個数を長さ方向にとり、対象画素６０Ａの個数の大小を長短で表現し、画素値の範囲の高低を色の違いで表現した画素値分布６２Ａを生成する。より具体的には、画素値分布生成部５３は、ヒストグラム６１の各画素値の範囲の対象画素６０Ａの個数を示す縦棒を、画素値の範囲の高い横軸の右側から順に積み重ねて、各画素値の範囲の対象画素６０Ａの個数の大小を示す直線状の帯を作り、この帯に画素値の範囲の高低に応じた色をつけて画素値分布６２Ａとする。
【００５７】
図１２では、画素値の範囲の高低をハッチングの濃淡で示しており、ハッチングの濃さが薄い程、画素値の範囲が高いことを示している。画素値の範囲の高低の表現は、例えば、画素値の範囲の高い部分を白、中間の部分をグレー、低い部分を黒とし、画素値の範囲の高低を白、グレー、および黒のグラデーションで表す。あるいは、白、グレー、および黒のモノクロではなく、画素値の範囲の高い部分を赤、中間の部分を緑、低い部分を青で表現する等、多色で表現してもよい。
【００５８】
画素値分布６２Ａは、半径Ｒ（Ｋ）が小さい程、対象画素６０Ａの個数が少ないので長さが短く、半径Ｒ（Ｋ）が小さい程、病変２２の中心Ｃ１近くのため画素値の範囲の高い部分の占める割合が多いものとなる。また、球面ＳＰ（Ｋ）の面積は半径Ｒ（Ｋ）の２乗に比例し、対象画素６０Ａの個数も半径Ｒ（Ｋ）の２乗に比例するので、画素値分布６２Ａの長さも、半径Ｒ（Ｋ）の２乗に比例するものとなる。画素値分布取得部４３は、画素値分布生成部５３で生成した画素値分布６２Ａを画素値分布画像生成部４４に出力する。
【００５９】
図１３に示すように、画素値分布画像生成部４４は、画素値分布取得部４３で取得した画素値分布６２Ａを半径Ｒ（Ｋ）に関して並べることで、画素値分布６２Ａを２次元で表現した画素値分布画像６３Ａを生成する画素値分布画像生成機能を担う。より具体的には、画素値分布画像生成部４４は、半径Ｒ（０）の画素値分布６２Ａの右側に、半径Ｒ（１）の画素値分布６２Ａの下端を揃えて並べ、これを半径Ｒ（Ｎ）の画素値分布６２Ａまで順に繰り返すことで、図に示す右上がりの階段状の画素値分布画像６３Ａを生成する。
【００６０】
画素値分布画像６３Ａは、半径Ｒ（Ｋ）が小さい程、画素値の範囲の高い部分の占める割合が多く、半径Ｒ（Ｋ）が大きくなるにつれ、画素値の範囲の低い部分の占める割合が多くなる。画素値分布画像６３Ａの半径Ｒ（Ｎ）付近の色が変化する境目が病変２２と正常領域の境界を示し、ひいては病変２２の大きさを示す。画素値分布画像生成部４４は、生成した画素値分布画像６３Ａを入出力制御部４０に出力する。
【００６１】
画素値分布取得部４３は、断層画像２０Ｎおよび２０Ｐのそれぞれに基づき算出した病変２２の中心Ｃ１の座標に基づいて、画素値分布６２ＡＮおよび６２ＡＰを生成する。画素値分布画像生成部４４は、画素値分布取得部４３で生成された画素値分布６２ＡＮおよび６２ＡＰに基づき、断層画像２０Ｎおよび２０Ｐに対応する画素値分布画像６３ＡＮおよび６３ＡＰを生成する。画素値分布画像生成部４４は、画素値分布画像６３ＡＮおよび６３ＡＰに、断層画像２０Ｎおよび２０Ｐに付した病変ＩＤを付す。これにより、断層画像２０Ｎおよび２０Ｐと、これに対応する画素値分布画像６３ＡＮおよび６３ＡＰが病変ＩＤによって関連付けられる。
【００６２】
病変抽出部５０で病変２２が複数個抽出された場合は、抽出された複数個の病変２２に対して、座標の算出、画素値分布６２Ａの生成、画素値分布画像６３Ａの生成をそれぞれ行う。したがって病変２２が複数個抽出された場合、画素値分布画像６３Ａは抽出された病変２２の個数分生成される。画素値分布画像６３Ａが複数個の病変２２のうちのいずれの病変２２に基づき生成したものであるかは、病変ＩＤにより区別する。
【００６３】
経過観察において、診療科端末１３は、図１４に示す比較読影ウィンドウ７０Ａをディスプレイ１６に表示させる。比較読影ウィンドウ７０Ａは、断層画像２０Ｎおよび２０Ｐを比較読影するための画面であり、取得要求に応じて医用画像ＤＢサーバ１２から送信された断層画像２０Ｎおよび２０Ｐを左右に並べて表示する表示枠７１Ａおよび７１Ｂと、被写体情報等の断層画像２０Ｎおよび２０Ｐに共通する付帯情報を表示する表示枠７２とを有する。表示枠７１Ａおよび７１Ｂの上方には、ファイルＩＤを表示する表示枠７３、撮影日時を表示する表示枠７４が設けられている。
【００６４】
比較読影ウィンドウ７０Ａには、画素値分布画像６３Ａの表示を指示する指示ボタン７５が配置されている。この指示ボタン７５が選択されると、診療科端末１３は、表示枠７１Ａおよび７１Ｂに表示されている断層画像２０Ｎおよび２０Ｐが属する断層画像群２１Ｎおよび２１Ｐの解析要求５５を支援装置１４に送信する。
【００６５】
指示ボタン７５が選択された場合、比較読影ウィンドウ７０Ａは、図１５に示す比較読影ウィンドウ７０Ｂに表示が遷移する。比較読影ウィンドウ７０Ｂには、表示枠７１Ａおよび７１Ｂの下部に、画素値分布画像６３ＡＮおよび６３ＡＰを表示する表示枠７６Ａおよび７６Ｂが追加される。表示枠７１Ａおよび７１Ｂに表示される断層画像２０Ｎおよび２０Ｐは、表示枠７６Ａおよび７６Ｂに表示された画素値分布画像６３ＡＮおよび６３ＡＰと病変ＩＤで関連付けられたものである。
【００６６】
図１６に示すように、表示枠７６Ａおよび７６Ｂは、画素値分布画像６３ＡＮおよび６３ＡＰをより細部まで比較読影することを可能とするため、比較読影ウィンドウ７０Ｂから独立して比較読影ウィンドウ７０Ｃで拡大表示することが可能となっている。画素値分布画像６３ＡＮおよび６３ＡＰを拡大表示することによって、病変２２と正常領域の境界の僅かな濃度の変化等も見逃さずに認識することができる。
【００６７】
以下、上記構成による作用について、図１７および図１８を参照して説明する。まず、診断医により、診療科端末１３からＲＩＳに撮影オーダが発行される。撮影オーダは、ＲＩＳからＣＴ装置１１に送信され、ＣＴ装置１１のコンソールのディスプレイに表示される。使用者は、撮影オーダの内容をディスプレイで確認し、撮影オーダに応じた撮影条件をＣＴ装置１１に設定して断層画像撮影を行う。これにより得られた複数枚の断層画像２０で構成される断層画像群２１は、同じ符号のファイルＩＤが割り振られた画像ファイル２３として医用画像ＤＢサーバ１２に送信され、医用画像ＤＢサーバ１２にて保管および管理される。
【００６８】
診断医は、診療科端末１３を操作して経過観察を行う。診断医は、被写体ＩＤや撮影日時を検索キーとする断層画像群２１Ｎおよび２１Ｐの取得要求を、診療科端末１３を通じて医用画像ＤＢサーバ１２に送信する。医用画像ＤＢサーバ１２は、診療科端末１３からの取得要求に応じた断層画像群２１Ｎおよび２１Ｐを検索し、検索した断層画像群２１Ｎおよび２１Ｐを診療科端末１３に送信する。
【００６９】
診療科端末１３では、医用画像ＤＢサーバ１２からの断層画像群２１Ｎおよび２１Ｐが受信され、これに基づき比較読影ウィンドウ７０Ａがディスプレイ１６に表示される。診断医は、この比較読影ウィンドウ７０Ａに映る断層画像２０Ｎおよび２０Ｐを比較読影し、経過観察を行う。
【００７０】
比較読影中、診断医は、画素値分布画像６３ＡＮおよび６３ＡＰを表示させるために指示ボタン７５を選択する。図１７（Ａ）に示すように、指示ボタン７５が選択されると、その際に表示枠７１Ａおよび７１Ｂに表示されている断層画像２０Ｎおよび２０Ｐが属する断層画像群２１Ｎおよび２１Ｐの解析要求５５が診療科端末１３から支援装置１４に送信される。
【００７１】
支援装置１４では、診療科端末１３からの解析要求５５が通信部３３で受信される。図１８のステップＳ１０に示すように、通信部３３で受信された解析要求５５は、入出力制御部４０で受け付けられる。図１７（Ｂ）および図１８のステップＳ１１に示すように、入出力制御部４０で受け付けられた解析要求５５に対応する取得要求５６が、通信部３３を介して医用画像ＤＢサーバ１２に送信される。医用画像ＤＢサーバ１２は、支援装置１４からの取得要求５６に応じた断層画像群２１Ｎおよび２１Ｐを検索し、図１７（Ｃ）に示すように、検索した断層画像群２１Ｎおよび２１Ｐを支援装置１４に送信する。
【００７２】
支援装置１４では、医用画像ＤＢサーバ１２からの断層画像群２１Ｎおよび２１Ｐが通信部３３で受信される。通信部３３で受信された断層画像群２１Ｎおよび２１Ｐは、入出力制御部４０で受け付けられる。入出力制御部４０で受け付けられた断層画像群２１Ｎおよび２１Ｐは、画像取得部４１に出力される。これにより断層画像群２１Ｎおよび２１Ｐが画像取得部４１で取得される（図１８のステップＳ１２）。断層画像群２１Ｎおよび２１Ｐは、座標取得部４２および画素値分布取得部４３に出力される。
【００７３】
座標取得部４２では、病変抽出部５０により断層画像群２１Ｎおよび２１Ｐに含まれる病変２２が抽出され（図１８のステップＳ１３）、その抽出結果が座標算出部５１に出力される。座標算出部５１では、病変２２の中心Ｃ１の座標が算出される（図１８のステップＳ１４）。これにより病変２２の中心Ｃ１の座標が取得される（図１８のステップＳ１５）。病変２２の中心Ｃ１の座標は、画素値分布取得部４３に出力される。
【００７４】
病変２２は、隣接する正常領域と峻別可能な明瞭な境界があるとは限らないため、病変２２をより正確に抽出しようとすると複雑な処理が必要で処理に時間が掛かる。しかし、病変抽出部５０による病変２２の抽出は、後段の座標算出部５１による病変２２の中心Ｃ１の座標を算出するための前処理であるため、病変２２を正確に抽出する必要はなく、病変２２として抽出した領域５７に正常領域が一部含まれていても一向に構わない。また、より正確に病変２２を抽出したとしても、抽出精度が低い方法で病変２２を抽出した場合と比べて、病変２２の中心Ｃ１の座標に大きな誤差は生じないため、病変抽出部５０による病変２２の抽出は、高度な正確性は求められない。したがって、病変２２の抽出を簡単な処理で短時間に済ませることができる。
【００７５】
画素値分布取得部４３では、ヒストグラム生成部５２により、画素値の範囲毎に計数した対象画素６０Ａの個数を度数で表し、かつ画素値の範囲を階級で表すヒストグラム６１が半径Ｒ（Ｋ）毎に生成される（図１８のステップＳ１６）。そして、このヒストグラム６１に基づき、画素値分布生成部５３で画素値分布６２ＡＮおよび６２ＡＰが生成される（図１８のステップＳ１７）。これにより画素値分布６２ＡＮおよび６２ＡＰが取得される（図１８のステップＳ１８）。画素値分布６２ＡＮおよび６２ＡＰは、画素値分布画像生成部４４に出力される。
【００７６】
画素値分布画像生成部４４では、画素値分布６２ＡＮおよび６２ＡＰを半径Ｒ（Ｋ）に関して並べた画素値分布画像６３ＡＮおよび６３ＡＰが生成される（図１８のステップＳ１９）。画素値分布画像６３ＡＮおよび６３ＡＰは入出力制御部４０に出力される。入出力制御部４０に出力された画素値分布画像６３ＡＮおよび６３ＡＰは、通信部３３に出力され、図１７（Ｄ）および図１８のステップＳ２０に示すように、通信部３３を介して診療科端末１３に送信される。
【００７７】
診療科端末１３では、比較読影ウィンドウ７０Ｂがディスプレイ１６に表示され、断層画像２０Ｎおよび２０Ｐに加えて、画素値分布画像６３ＡＮおよび６３ＡＰが閲覧可能となる。診断医は、必要に応じて比較読影ウィンドウ７０Ｃを表示し、画素値分布画像６３ＡＮおよび６３ＡＰを拡大して詳細に観察する。
【００７８】
断層画像群２１Ｎおよび２１Ｐに含まれる病変２２は、今回と前回の各撮影時の被写体の位置、被写体の向きの変化や、被写体の呼吸、被写体の心拍等の周辺臓器の動きの変化等によって、病変２２の位置および病変２２の向きが違ってくるが、病変２２の位置の違いは、断層画像２０Ｎおよび２０Ｐの病変２２の３次元的な中心Ｃ１の座標の１点の基準点による位置合わせだけで解消することができる。病変２２の向きの違いは、病変２２に対して異方性がない半径Ｒ（Ｋ）の球面ＳＰ（Ｋ）上の対象画素６０Ａの画素値の分布状況を取得することで相殺することができ、病変２２の向きは正確に合わせる必要がない。したがって、断層画像２０の解像度（スライス厚、画素数）や病変２２自体の撓みを考慮する必要もない。病変２２自体に撓みがあった場合でも、画素値分布画像６３Ａ全体としてみれば撓みがない場合と比べて僅かな差であるため、病変２２が悪性病変であるか否かの判断への病変２２自体の撓みの影響を相殺することができる。
【００７９】
早期癌の初期状態である可能性がある大きさが１ｃｍ程度の病変２２に対しては、位置合わせ、および病変の抽出を正確に行うことがより難しくなるが、本発明では上述のように正確な位置合わせ、および正確な病変の抽出は必要ないため、より効果を発揮する。このため、癌の早期発見および迅速な処置に寄与することができる。
【００８０】
画素値分布画像６３Ａは、それ自体で、特許文献１の成長度のような定量的な数値では分からない病変２２の大きさや濃度を診断医に視覚的に認識させることができる。このため、画素値分布画像６３ＡＮおよび６３ＡＰを診療科端末１３に出力して閲覧可能とすることで、病変２２の大きさの拡大傾向、濃度の増加傾向といった、病変２２が悪性病変であるか否かを判断する際の重要な判断材料を診断医に提供することができる。
【００８１】
また、従来は、１つの病変２２に対して、その病変２２が悪性病変であるか否かを判断する場合に、複数枚の断層画像２１Ｎおよび２１Ｐを比較読影する必要があり時間が掛かっていたが、本実施形態では、１つの病変２２に対して１枚の画素値分布画像６３ＡＮおよび６３ＡＰを比較読影するだけで病変２２が悪性病変であるか否かを判断することができ、経過観察の作業を効率化することができる。
【００８２】
画素値分布画像６３Ａは、断層画像２０の画素値ではなく、その分布を示すものであるため、撮影時の被写体の位置、被写体の向きの変化や、被写体の呼吸、被写体の心拍等の周辺臓器の動きの変化等によって、断層画像群２０Ｎと断層画像２０Ｐとで病変２２およびその周辺近傍の正常領域の画素値に変化があったとしても、病変２２の大きさの拡大傾向、濃度の増加傾向の認識に与える影響は少ない。
【００８３】
なお、病変２２の中心の座標の代わりに、病変２２の重心の座標を取得してもよい。病変２２の重心は、病変２２として抽出した領域５７の画素値の高低を考慮して決定する。このため病変２２の重心は、図６において符号Ｃ２で示すように、領域５７の幾何学的な中心Ｃ１とは必ずしも一致しない。例えば領域５７の一方に画素値が高い部分が偏っていた場合、幾何学的な中心を病変２２の画素値が高い部分に向けてシフトさせて、病変２２の重心とする。図６に示す重心Ｃ２は、右上側に偏在した画素値の高い部分に応じて中心Ｃ１から右上側にシフトしている。
【００８４】
幾何学的な中心Ｃ１をどれだけシフトさせるかは、領域５７内での画素値が高い部分と他の部分との画素値の差や、領域５７における画素値が高い部分が占める割合等に応じて決定する。
【００８５】
また、座標取得部４２で取得した座標から半径Ｒ（Ｋ）の球面ＳＰ（Ｋ）上の断層画像群２１の画素６０を対象画素６０Ａとする代わりに、座標取得部４２で取得した座標から半径Ｒ（Ｋ）の円の円周上の断層画像２０の画素６０を対象画素６０Ａとしてもよい。この場合、対象画素６０Ａは、例えば図９Ａに示す断層画像２０Ａの２４個の対象画素６０Ａのみとなり、図１０Ａに示す断層画像２０Ｂおよび２０Ｃの１６個の対象画素６０Ａは除外される。
【００８６】
［第２実施形態］
上記第１実施形態では、支援装置１４で、断層画像群２１Ｎおよび２１Ｐに基づき画素値分布画像６３ＡＮおよび６３ＡＰを生成し、生成した画素値分布画像６３ＡＮおよび６３ＡＰを診療科端末１３に出力しているが、画素値分布画像６３ＡＮおよび６３ＡＰの差分画像を生成し、生成した差分画像と画素値分布画像６３ＡＮおよび６３ＡＰとを診療科端末１３に出力してもよい。
【００８７】
図１９において、本実施形態の支援装置のＣＰＵ８０は、上記第１実施形態の各機能部４０〜４４の他に、差分画像生成部８１として機能する。差分画像生成部８１は、画素値分布画像６３ＡＮおよび６３ＡＰを画素値分布画像生成部４４から受け取る。そして、図２０に示す画素値分布画像６３ＡＮおよび６３ＡＰの差分画像８２を生成する。
【００８８】
画素値分布画像６３ＡＮおよび６３ＡＰは、画素値の範囲の高低を色の違いで表現したものであるため、差分画像生成部８１は、例えば、画素値分布画像６３ＡＮおよび６３ＡＰで画素値の範囲が同じ部分については黒（図２０で濃いハッチングで示す）とし、画素値の範囲が何段階異なるかでグレーから白（図２０で薄いハッチングで示す）に変化するよう色付けした差分画像８２を生成する。例えば、画素値分布画像６３ＡＮで１００〜１９９、画素値分布画像６３ＡＰで０〜９９で、画素値の範囲が１段階異なる部分については黒に近いグレーとし、画素値分布画像６３ＡＮで１００〜１９９、画素値分布画像６３ＡＰで−５００〜−４０１で、画素値の範囲が６段階異なる部分については白に近いグレーとする。差分画像生成部８１は、生成した差分画像８２を入出力制御部４０に出力する。
【００８９】
入出力制御部４０は、画素値分布画像６３ＡＮおよび６３ＡＰに加えて、差分画像８２を通信部３３に出力し、通信部３３を介して診療科端末１３に送信させる。診療科端末１３は、画素値分布画像６３ＡＮおよび６３ＡＰに加えて、差分画像８２を比較読影ウィンドウ７０Ｂ、７０Ｃに表示する。
【００９０】
差分画像８２の全体に変化がない場合は、画素値分布画像６３ＡＮおよび６３ＡＰに差がない、すなわち病変２２に経時変化がないことが分かり、反対に差分画像８２に変化がある場合は、病変２２が成長もしくは減退しており、病変２２に経時変化があることが分かる。このため、画素値分布画像６３ＡＮおよび６３ＡＰの表示に加えて差分画像８２を表示することで、より病変２２の大きさの拡大傾向、濃度の増加傾向を正しく認識することができる。
【００９１】
［第３実施形態］
上記各実施形態では、長さが半径Ｒ（Ｋ）の２乗に比例する直線状の画素値分布６２Ａを半径Ｒ（Ｋ）に関して単純に並べて、階段状の画素値分布画像６３Ａを生成しているが、画素値分布６２および画素値分布画像６３の形状は直線状および階段状に限らない。例えば図２１に示すように、画素値分布生成部５３で、ヒストグラム６１を元に生成した直線状の画素値分布６２Ａから、円弧状の画素値分布６２Ｂを生成し、この円弧状の画素値分布６２Ｂを、画素値分布画像生成部４４で半径Ｒ（Ｋ）に関して並べ、扇形状の画素値分布画像６３Ｂを生成してもよい。円弧状の画素値分布６２Ｂは、上記第１実施形態の画素値分布６２Ａと同じく、半径Ｒ（Ｋ）の２乗に比例するものとなる。
【００９２】
階段状の画素値分布画像６３Ａは、矩形状の断層画像２０に対して形状が著しく異なり、また、半径Ｒ（Ｋ）が大きい程、画素値分布６２Ａの長さが長くなり、半径Ｒ（Ｋ）が小さい場合の長さが短い画素値分布６２Ａと比べて目立つため、観察にある程度の慣れが必要である。対して、円弧状の画素値分布６２Ｂは、階段状の画素値分布画像６３Ａと比べて断層画像２０の形状に似ており、半径Ｒ（Ｋ）が大きい場合の長さが長い画素値分布６２Ｂも短い画素値分布６２Ｂと比べて目立たないため、断層画像２０に近い感覚で観察することができる。
【００９３】
［第４実施形態］
上記各実施形態では、長さが異なる画素値分布６２Ａ、６２Ｂを半径Ｒ（Ｋ）に関して並べて、階段状の画素値分布画像６３Ａまたは扇形状の画素値分布６２Ｂを生成しているが、本発明はこれに限定されない。画素値分布６２の長さを正規化して揃え、この長さを揃えた画素値分布６２を半径Ｒ（Ｋ）に関して並べ、矩形状の画素値分布画像６３を生成してもよい。
【００９４】
この場合、ヒストグラム生成部５２は、画素値の分布量として、上記各実施形態の画素値の範囲毎の対象画素６０Ａの個数ではなく、画素値の範囲毎の対象画素６０Ａの個数を全対象画素６０Ａの個数で除算した、画素値の範囲毎の対象画素６０Ａの個数の割合を算出する。全対象画素６０Ａの個数とは、ある球面ＳＰ（Ｋ）に対する対象画素６０Ａの個数であり、図９および図１０の例でいえば、全対象画素６０Ａの個数は５６個である。
【００９５】
例えば、−１０００〜−９０１の画素値の範囲の対象画素６０Ａの個数が４個で、全対象画素６０Ａの個数が１００個であった場合、−１０００〜−９０１の画素値の範囲の画素値の分布量は、４／１００＝０．０４となる。各画素値の範囲の対象画素６０Ａの個数の割合を加算すると１になる。ヒストグラム生成部５２は、画素値の範囲毎の対象画素６０Ａの個数の割合を度数で表し、画素値の範囲を階級で表すヒストグラムを生成する。
【００９６】
画素値分布生成部５３は、上記第１実施形態と同様に、ヒストグラムの画素値の範囲毎の対象画素６０Ａの個数の割合を示す縦棒を、画素値の範囲の高いものから順に積み重ねて、対象画素６０Ａの個数の割合の大小を示す直線状の帯を作り、この帯に画素値の範囲の高低に応じた色をつけた画素値分布６２Ｃ（図２２参照）を生成する。上記第１、第３実施形態の画素値分布６２Ａ、６２Ｂは、半径Ｒ（Ｋ）が小さい程、対象画素６０Ａの個数が少ないので長さが短くなるが、本実施形態の画素値分布６２Ｃは、加算すると１になる対象画素６０Ａの個数の割合を画像化したものであるため、半径Ｒ（Ｋ）に関わらず同じ長さとなる。したがって、画素値分布６２Ｃは、画素値分布６２Ａの長さを正規化して揃えたものとなる。
【００９７】
図２２に示すように、画素値分布画像生成部４４は、長さを正規化して揃えた直線状の画素値分布６２Ｃを半径Ｒ（Ｋ）に関して並べ、矩形状の画素値分布画像６３Ｃを生成する。この矩形状の画素値分布画像６３Ｃは、階段状の画素値分布画像６３Ａおよび円弧状の画素値分布６２Ｂのような特徴的な形状でなく、断層画像２０と同じ形状であるため、より断層画像２０に近い感覚で観察することができる。
【００９８】
［第５実施形態］
上記各実施形態では、対象画素６０Ａの個数や割合といった画素値の分布量を長さ方向にとり、画素値の分布量の大小を長短で表現し、画素値の範囲の高低を色の違いで表現した画素値分布６２Ａ〜６２Ｃを生成しているが、本発明はこれに限定されない。画素値の範囲を長さ方向にとり、画素値の分布量の大小を色の違いで表現した画素値分布を生成してもよい。
【００９９】
図２３において、画素値分布生成部５３は、ヒストグラム生成部５２で生成したヒストグラム６１を元に、画素値の範囲を長さ方向にとり、画素値の分布量の大小を色の違いで表現した画素値分布６２Ｄを生成する。より具体的には、画素値分布生成部５３は、直線状の帯に、ヒストグラム６１の各画素値の範囲に応じて等間隔で区切った複数の区画を作り、この区画に画素値の分布量の大小に応じた色をつけて画素値分布６２Ｄとする。
【０１００】
図２３では、画素値の分布量の大小をハッチングの濃淡で示しており、ハッチングの濃さが薄い程、画素値の分布量が大きいことを示している。画素値の分布量の大小の表現は、例えば、画素値の分布量が大きい部分を白、中間の部分をグレー、小さい部分（画素値の分布量が０の場合も含む）を黒とし、画素値の分布量の大小を白、グレー、および黒のグラデーションで表す。画素値分布６２Ａの場合と同じく、多色で表現してもよい。
【０１０１】
図２４に示すように、画素値分布画像生成部４４は、画素値分布６２Ｄを半径Ｒ（Ｋ）に関して並べ、矩形状の画素値分布画像６３Ｄを生成する。この矩形状の画素値分布画像６３Ｄは、上記第４実施形態の画素値分布画像６３Ｃと同様、断層画像２０と同じ形状であるため、より断層画像２０に近い感覚で観察することができる。
【０１０２】
上記各実施形態では、病変２２の中心Ｃ１の座標から半径Ｒ（Ｋ）の球面ＳＰ（Ｋ）内に中心が含まれ、かつ球面ＳＰ（Ｋ）の球面上の位置にある断層画像群２１の画素６０を対象画素６０Ａとしているが、球面ＳＰ（Ｋ）内に中心が含まれているか否かに関わらず、球面ＳＰ（Ｋ）の球面上の位置にある断層画像群２１の画素６０を対象画素６０Ａとしてもよい。
【０１０３】
また、図２５に示すように、球面ＳＰ（Ｋ）を複数の微小領域に分割し、この微小領域を対象画素６０Ａとしてもよい。この場合の対象画素６０Ａは、球面ＳＰ（Ｋ）上に作成したいわば仮想的な画素であるため、このままでは画素値をもたない。そこで、対象画素６０Ａの画素値を、対象画素６０Ａの近傍の断層画像２０の画素６０の画素値から補間により算出する。例えば、ある対象画素６０Ａの近傍の断層画像２０の画素６０が４個あり、４個の画素６０の画素値がそれぞれ１５、２５、３８、４６であった場合、その対象画素６０Ａの画素値は、４個の画素６０の画素値の平均値を補間値とすると、（１５＋２５＋３８＋４６）／４＝３１となる。
【０１０４】
この場合、ヒストグラム生成部５２は、画素値の分布量として、画素値の範囲毎の対象画素６０Ａの面積（微小領域の面積）を積算し、対象画素６０Ａの面積を度数で表し、かつ画素値の範囲を階級で表すヒストグラムを半径Ｒ（Ｋ）毎に生成する。画素値分布生成部５３は、画素値の範囲毎の対象画素６０Ａの面積を長さ方向にとり、対象画素６０Ａの面積の大小を長短で表現し、画素値の範囲の高低を色の違いで表現した画素値分布６２を生成する。画素値分布６２は、上記第１実施形態の画素値分布６２Ａのように直線状であってもよいし、上記第３実施形態の画素値分布６２Ｂのように円弧状であってもよい。あるいは、上記第５実施形態の画素値分布６２Ｃのように、画素値の範囲を長さ方向にとり、画素値の範囲毎の対象画素６０Ａの面積の大小を色の違いで表現してもよい。
【０１０５】
画素値の分布量としては、画素値の範囲毎の対象画素６０Ａの面積に代えて、画素値の範囲毎の対象画素６０Ａの面積を、球面ＳＰ（Ｋ）の面積で除算した対象画素６０Ａの面積の割合を採用してもよい。画素値の分布量を対象画素６０Ａの面積の割合とすれば、第４実施形態の画素値分布６２Ｄと同様に、長さを正規化して揃えた画素値分布６２とすることができる。
【０１０６】
上記第１実施形態では、断層画像２０Ｎおよび２０Ｐのそれぞれに基づき算出した座標が一致または類似する病変２２を同一の病変２２として認識しているが、本発明はこれに限定されない。例えば、比較読影ウィンドウ７０Ａにて、断層画像２０Ｎおよび２０Ｐ上で同一と思しき病変２２の領域を診断医に指定させ、指定された領域に対して病変の抽出、座標の算出、画素値分布６２および画素値分布画像６３の生成を行ってもよい。この場合は、解析要求５５に、病変２２と思しき領域が指定された断層画像２０Ｎおよび２０Ｐの画像ファイル２３のファイルＩＤと、指定された領域の座標情報とを付す。支援装置１４では指定された領域にのみ病変の抽出等を行う。このため、断層画像群２１Ｎおよび２１Ｐを構成する複数枚の断層画像２０Ｎおよび２０Ｐの全てに対して病変２２の抽出等を行うよりも効率的である。
【０１０７】
病変２２の中心Ｃ１または重心Ｃ２の座標を手動で指定してもよい。例えば、比較読影ウィンドウ７０Ａにて、断層画像２０Ｎおよび２０Ｐ上で診断医が病変２２の中心Ｃ１または重心Ｃ２と思しき画素６０を指定する。この場合は、解析要求５５に、指定された画素６０の座標を付す。座標取得部４２は、入出力制御部４０から解析要求５５を受け取り、解析要求５５に付された座標を取得する。このようにした場合、病変の抽出および座標の算出が不要であるため、処理を簡素化することができる。ただし、病変２２の中心Ｃ１または重心Ｃ２が診断医の感覚に委ねられるため、より正確性を重視する場合は上記各実施形態のように病変を抽出して座標を算出したほうが好ましい。
【０１０８】
気管支分岐部のスライス位置情報等の生体臓器の構造物の位置に基づき、断層画像群２１Ｎおよび２１Ｐの３次元的な位置合わせをしたうえで病変の抽出等を行ってもよい。この場合、断層画像群２１Ｎのみに対して病変の抽出、座標の算出を行い、断層画像群２１Ｎに対する病変の抽出結果、座標の算出結果を、断層画像群２１Ｐに適用してもよい。
【０１０９】
診療科端末１３から支援装置１４に解析要求５５を送信する際には、診療科端末１３は医用画像ＤＢサーバ１２からの断層画像群２１Ｎおよび２１Ｐを受信しているので、解析要求５５に断層画像群２１Ｎおよび２１Ｐを添付してもよい。この場合は、支援装置１４から医用画像ＤＢサーバ１２に、解析要求５５に対応する取得要求５６を送信する必要がなくなる。
【０１１０】
上記各実施形態では、医用画像ＤＢサーバ１２および診療科端末１３とは別に支援装置１４を設けているが、医用画像ＤＢサーバ１２または診療科端末１３に支援装置１４の機能を担わせてもよい。また、診療科端末１３と支援装置１４の機能を一体化し、支援装置１４のディスプレイに比較読影ウィンドウ７０Ａ〜７０Ｃを表示してもよい。この場合、出力制御部は、ディスプレイへの画素値分布画像６３の出力制御機能を担う。
【０１１１】
画素値分布画像６３の送信先としては、上記各実施形態の診療科端末１３に限らず、断層画像群２１の読影業務を担う読影医が操作する読影端末でもよい。また、断層画像群２１を医用画像ＤＢサーバ１２から取得するのではなく、断層画像群２１が記憶されたコンパクトディスク等の記憶媒体から取得してもよい。さらに、ＣＴ装置１１や診療科端末１３は１台でも複数台でもよい。また、ＣＴ装置１１に加えて、あるいは代えて、ＭＲＩ装置等の他のモダリティを設けてもよい。
【０１１２】
上記各実施形態では、今回と前回撮影した２つの断層画像群を元に画素値分布画像を生成する例を説明しているが、今回、前回、前々回の３つの断層画像群等、２つ以上の断層画像群を元に画素値分布画像を生成してもよい。
【０１１３】
なお、病変２２には、最終的に診断医によって例えば癌であると確定診断が下された悪性病変等の狭義の病変（請求項１、１１、１２における「病変」に相当する）の他に、診断の初期において悪性が疑われて悪性病変の候補とされたが、最終的には良性病変と確定診断された病変候補（請求項１、１１、１２における「病変候補」に相当する）が含まれる。
【符号の説明】
【０１１４】
１４  医用画像診断支援装置（支援装置）
２０  断層画像
２１、２１Ｎおよび２１Ｐ  断層画像群
２２  病変
３２、８０  ＣＰＵ
３５  医用画像診断支援プログラム（支援プログラム）
４０  入出力制御部（出力制御部）
４１  画像取得部
４２  座標取得部
４３  画素値分布取得部
４４  画素値分布画像生成部
５２  ヒストグラム生成部
５３  画素値分布生成部
６０Ａ  対象画素
６１  ヒストグラム
６２Ａ、６２Ｂ、６２Ｃ、６２Ｄ  画素値分布
６３Ａ、６３ＡＮおよび６３ＡＰ、６３Ｂ、６３Ｃ、６３Ｄ  画素値分布画像
８１  差分画像生成部
８２  差分画像
【図１】
2015221141.tif 000003
【図２】
2015221141.tif 000004
【図３】
2015221141.tif 000005
【図４】
2015221141.tif 000006
【図５】
2015221141.tif 000007
【図６】
2015221141.tif 000008
【図７】
2015221141.tif 000009
【図８】
2015221141.tif 000010
【図９Ａ】
2015221141.tif 000011
【図９Ｂ】
2015221141.tif 000012
【図１０Ａ】
2015221141.tif 000013
【図１０Ｂ】
2015221141.tif 000014
【図１１】
2015221141.tif 000015
【図１２】
2015221141.tif 000016
【図１３】
2015221141.tif 000017
【図１４】
2015221141.tif 000018
【図１５】
2015221141.tif 000019
【図１６】
2015221141.tif 000020
【図１７】
2015221141.tif 000021
【図１８】
2015221141.tif 000022
【図１９】
2015221141.tif 000023
【図２０】
2015221141.tif 000024
【図２１】
2015221141.tif 000025
【図２２】
2015221141.tif 000026
【図２３】
2015221141.tif 000027
【図２４】
2015221141.tif 000028
【図２５】
2015221141.tif 000029