密度分析方法、密度分析方法を実装したコンピュータプログラムおよび密度分析システム

【課題】Ｘ線ＣＴ装置で被検体を撮像して得られる撮像データに基づいて被検体の密度を分析する密度分析について分析精度を高める。
【解決手段】Ｘ線ＣＴ装置で被検体を撮像して得られる撮像データに基づいて被検体の密度を分析する密度分析方法において、撮像データに設定される複数の密度分析対象部位について、密度分析対象部位ごとにそれを中心にして複数の画素からなる密度検査領域を設定し、その密度検査領域における複数の画素の各密度値の平均値を求め、求めた平均値を密度分析対象部位の密度値として密度分析を行うようにしている。

【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線ＣＴ装置で得られる撮像データを用いて被検体の密度を分析する密度分析方法、密度分析方法を実装したコンピュータプログラムおよび密度分析システムに関する。
【背景技術】
【０００２】
Ｘ線ＣＴ装置では、加速器やＸ線管などのＸ線源から照射されるＸ線が被検体を透過する際に吸収・減衰される量を検出することで得られる撮像データに画像再構成処理を施すことにより被検体の２次元断面画像や３次元画像を生成する。被検体を透過する際のＸ線の吸収・減衰はそのＸ線が透過した部分における被検体の密度に相関する。すなわちＸ線ＣＴ装置による被検体の画像を矩形の小領域などとして構成する画素におけるＣＴ値は、各画素に対応する被検体の部位の密度に相関した値を持つ。このことからＸ線ＣＴ装置による撮像データを用いて被検体を密度に関して分析することが行われている。例えばＸ線ＣＴ装置で撮像した断層像から電池内の活性物質の密度分布を求めることにより電池特性を評価し、製作時や使用時における電池の品質をチェックする技術がその例である（特許文献１）。この他に例えば特許文献２〜特許文献４などにもＸ線ＣＴ装置による撮像データから被検体の密度分析を行う技術が開示されている。
【０００３】
【特許文献１】特許第３１２５１２０号公報
【特許文献２】特開平７−１３４０８８号公報
【特許文献３】特開平８−３２７５２５号公報
【特許文献４】特開平１０−２２７７３４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
Ｘ線ＣＴ装置による密度分析の有用性はますます高まる方向にあり、それに伴って、より微小な密度差まで分析する高精度な密度分析に対する要求が高まってきている。例えば、電池内の活性物質の密度分布による電池特性の評価では、放電または充電後における活性物質の密度分布の差を分析することが有用であるが、このような場合には、より高精度な密度分析が求められる。
【０００５】
しかし、Ｘ線ＣＴ装置には分解能がある。現在のＸ線ＣＴ装置における分解能は、被検体を構成する物質の密度の２％程度であるのが一般的であるが、その分解能の近辺ではノイズの影響が大きくなる。そのため分解能の近辺では有意な密度分析が行えなくなり、高精度な密度分析の要求に応えることができていないのが実情である。
【０００６】
本発明は、このような事情を背景になされたものであり、Ｘ線ＣＴ装置における分解能の近辺でも有意な密度分析を可能とする密度分析方法の提供を第１の目的とし、またそのような密度分析方法を実装したコンピュータプログラムの提供を第２の目的とし、さらにそのような密度分析方法を実行するための密度分析システムの提供を第３の目的としている。
【課題を解決するための手段】
【０００７】
上記第１の目的のために本発明では、Ｘ線ＣＴ装置で被検体を撮像して得られる撮像データに基づいて前記被検体の密度を分析する密度分析方法において、前記撮像データに設定される複数の密度分析対象部位について、前記密度分析対象部位ごとにそれを中心にして複数の画素からなる密度検査領域を設定し、前記密度検査領域における複数の画素の各密度値の平均値を求め、求めた平均値を前記密度分析対象部位の密度値として密度分析を行うようにしたことを特徴としている。
【０００８】
また上記第１の目的のために本発明では、Ｘ線ＣＴ装置で被検体を撮像して得られる撮像データに基づいて前記被検体の密度を分析する密度分析方法において、前記Ｘ線ＣＴ装置により同一の被検体について少なくとも２つ以上の撮像データをそれぞれ異なる時刻に取得するか、または形状、材質が同一の複数の被検体のそれぞれについて撮像データを取得する処理過程、前記撮像データ取得過程で得られた複数の撮像データについて位置合わせを行う処理過程、前記複数の撮像データのそれぞれに設定される複数の密度分析対象部位について、前記密度分析対象部位ごとにそれを中心にして複数の画素からなる密度検査領域を設定する処理過程、前記処理過程で設定された密度検査領域における複数の画素の各密度値の平均値を求める処理過程、および前記処理過程で求められた密度値の平均値を前記密度分析対象部位の密度値とし、前記複数の撮像データ間での対応密度分析対象部位間について密度差を求める処理過程を含むことを特徴としている。
【０００９】
上記第２の目的のために本発明では、上記のような密度分析方法を実装したコンピュータプログラムについて、前記Ｘ線ＣＴ装置により取得された撮像データを読み込む撮像データ読み込み機能、前記撮像データ読み込み機能で読み込まれた複数の撮像データについて位置合わせを行うための位置合わせ機能、前記複数の撮像データのそれぞれに設定される複数の密度分析対象部位について、前記密度分析対象部位ごとにそれを中心にして複数の画素からなる密度検査領域を設定するための密度検査領域設定機能、前記密度検査領域設定機能で設定された密度検査領域における複数の画素の各密度値の平均値を求めるための密度算出機能、および前記密度算出機能で求められた密度値の平均値を前記密度分析対象部位の密度値とし、前記複数の撮像データ間での対応密度分析対象部位間について密度差を求めるための密度差算出機能を含むものとしている。
【００１０】
上記第３の目的のために本発明では、Ｘ線ＣＴ装置で被検体を撮像して得られる撮像データに基づいて前記被検体の密度に関して分析する密度分析方法に用いられる密度分析システムにおいて、上記のような密度分析方法を実行するための密度分析装置を備えていることを特徴としている。
【００１１】
また本発明では、上記のような密度分析システムについて、前記密度分析装置は、上記のようなコンピュータプログラムを記憶するプログラム記憶装置と、前記プログラム記憶装置から読み出した前記コンピュータプログラムを実行する演算装置とで構成するものとしている。
【発明の効果】
【００１２】
本発明では、密度分析対象部位を中心にして複数の画素からなる密度検査領域を設定し、この密度検査領域について、そこに含まれる複数の画素の各密度値の平均値を求め、これを密度分析対象部位における代表密度値とするようにしている。このようにすることにより、画素ごとにランダムに生じるノイズの影響をキャンセルして大幅に低減することができる。この結果、密度分析の精度を大幅に高めることができ、Ｘ線ＣＴ装置における分解能の近辺でも有意な密度分析が可能となる。
【発明を実施するための最良の形態】
【００１３】
以下、本発明を実施する上で好ましい形態について説明する。図１に第１の実施形態による密度分析方法における処理の流れを示し、図２にその処理におけるＸ線ＣＴ装置による被検体の撮像の様子を模式化して示す。本実施形態は、同一の被検体における密度の時間的な変化、例えば被検体に対して何らかの処理を施した場合における処理前後の密度差を分析する場合の例である。そのような例としては、例えば放充電を繰り返す電池における放電または充電の前後での密度差分析がある。
【００１４】
被検体１を撮像するＸ線ＣＴ装置２は、加速器やＸ線管などで構成されＸ線２ｘを照射するＸ線源２ａ、Ｘ線源２ａからのＸ線２ｘを検出するＸ線検出器２ｂ、および被検体１を載せて回転させるターンテーブル２ｃを備えている。そしてターンテーブル２ｃを回転させながら被検体１にＸ線２ｘを照射しつつ、被検体１を透過したり被検体１の周囲を通過したりしたＸ線２ｘをＸ線検出器２ｂで検出して撮像データを取得する。
【００１５】
このようなＸ線ＣＴ装置２により時刻Ｔ１において被検体１を撮像し、被検体１の撮像データ１０１を取得する（処理１００１）。撮像データ１０１は、何らかの処理が施される前の被検体１の撮像データである。それから被検体１に所定の処理を施す（処理１００２）。この処理は、被検体１が放充電を繰り返す電池の場合であれば、放電や充電である。処理後の時刻Ｔ２で新たに被検体１をＸ線ＣＴ装置２により撮像して撮像データ１０２を取得する（処理１００３）。撮像データ１０１、１０２は、被検体１を異なる時刻Ｔ１、Ｔ２で別個に撮像して得られたものである。このため撮像データ１０１、１０２それぞれを再構成して生成される画像中における被検体１の位置は、図２に示すように、ずれを生じているのが通常である。そこで撮像データ１０１、１０２について位置合わせ、より具体的には撮像データ１０１、１０２における被検体１の位置を合わせる位置合わせを行う（処理１００４）。
【００１６】
撮像データの位置合わせは、被検体１の幾何学的形状に基づいて行うことができるし、また被検体の重心や慣性モーメント主軸に基づいて行うことができ、さらに撮像データ１０１、１０２の形状誤差を最小化する方法などで行うこともできる。以下では被検体１の幾何学的形状に基づく位置合せについて説明する。幾何学的形状に基づく位置合わせの例のイメージを図３に示す。幾何学的形状に基づいて撮像データ１０１と１０２の位置合わせを行うには、被検体１における基準となる点（基準原点）と、被検体１が有する特徴的な幾何学的形状の画像座標系４に対する傾き角度（回転角度）を求める必要がある。これらを求める一つの方法では、被検体１が交叉する直線を伴う輪郭を有するものである場合であれば、それを特徴的な幾何学的形状として、撮像データ１０１と撮像データ１０２の対応位置に、交叉する２本の直線８、８と直線８´、８´のそれぞれを定義する。撮像データ１０１に交叉する２本の直線８、８を定義するには、撮像データ１０１において、被検体１を表すデータ部と被検体１の周囲の空気を表すデータ部との境界近傍を画像上で複数点計測し、それで求まる複数の点に最善に適合する直線を計算から求める方法を用いることができる。直線８、８が定義できたらその交点を求め、それを基準原点５とする。また直線８の方向ベクトルから被検体１が画像座標系４に対してなす角度として回転角度を求める。撮像データ１０２についても同様に、直線８´、８´を定義してその交点を基準原点５´とし、直線８´の方向ベクトルから被検体１の回転角度を求める。
【００１７】
基準原点を求める他の方法としては、図３に併せて示す例のように、撮像データ１０１と撮像データ１０２の対応位置に円７、７´を定義し、その中心点を６、６´を基準原点とする方法が可能である。以上のようにして基準原点と回転角度が求まったら、撮像データ１０１と撮像データ１０２に対してアフィン変換（一次変換）を施して位置合わせをする。図４にアフィン変換後の撮像データ１０１´と撮像データ１０２´の例を示す。
【００１８】
ここで、被検体が位置合わせの指標となる幾何学的形状を有しない場合もある。そのような被検体の場合には、位置合わせの指標とし易い直線や平面あるいは円筒面などの幾何学的形状を有する位置合せ用指標体を被検体に取り付けて撮像を行うようにする。図５にその例を示す。図５の例では被検体１の輪郭が円弧形状の部分とぎざぎざな形状部分からなり、これらを回転角度の指標とすることが困難である。そこで、断面形状が四角形である位置合せ用指標体９を取り付け、これを被検体１と一緒に撮像する。このような被検体１については、その円弧形状輪郭部分を利用して基準原点を求める。すなわち、被検体１の円弧形状の輪郭部分に沿って複数の点を計測し、その複数の点に最善に適合する円弧１０、１０´の中心点として基準原点５、５´を求める。一方、被検体１の回転角度は位置合せ用指標体９から求める。すなわち、位置合せ用指標体９について図３に関して説明したのと同様にして直線８、８´を求め、この直線８、８´により被検体１の回転角度を求める。
【００１９】
Ｘ線ＣＴ装置２は、被検体１の２次元断面画像をそれに垂直な方向に積み重ねることで被検体１の３次元形状画像を得ることもできる。またＸ線ＣＴ装置２が図２の例のような扇状のＸ線２ｘ（ファンビーム）に代えて円錐状のＸ線（コーンビーム）を照射する場合には、直接に被検体１の３次元形状画像を得ることができる。被検体１について３次元形状の撮像データを得る場合には、その３次元形状撮像データについて位置合わせを行う。図６にそのような例を示す。図６の例のような３次元形状の被検体１については、まず撮像データ１０１における被検体１の上端面上の点を画像上で複数計測し、それらの点に最善に適合する平面１２を求める。それから、被検体１の内側面あるいは外側面上の点を画像上で複数計測し、それらの点に最善に適合する円筒面１１を求める。そして円筒面１１の中心軸１１ａと平面１２の交点１１ｐとして基準原点を求める。また撮像データ１０１の１つの軸を円筒面１１の中心軸１１ａとし、残りの２つの軸は、被検体１に取り付けた位置合せ用指標体９を利用して決定し、これら３つの軸から回転角度を求める。撮像データ１０２についても同様にして平面１２´と円筒面１１´を求め、それから円筒面１１´の中心軸１１ａ´と平面１２´の交点１１ｐ´として基準原点を求め、また円筒面１１´の中心軸１１ａ´として求めた１つの軸と被検体１に取り付けた位置合せ用指標体９を利用して求めた残りの２つの軸から回転角度を求める。
【００２０】
以上のようにして被検体１の位置合わせを終えたら、その撮像データ１０１´、１０２´のそれぞれについて密度検査領域を設定・配列する（処理１００５）。密度検査領域とは、被検体１を構成する物質の密度に対応するＣＴ値を撮像データから取得するについてＣＴ値をサンプリングする広がりの範囲を規定する領域である。この密度検査領域の設定は、撮像データに設定される複数のＣＴ値取得対象部位（密度分析対象部位）のそれぞれについて、その密度分析対象部位を中心にして適切な数の画素を密度分析対象部位ごとに設定することでなされる。適切な画素数つまり密度検査領域のサイズは、撮像データにおけるノイズの大きさに相関し、そのノイズの大きさは、被検体を構成物質の種類（被検体の材質）や大きさ、さらにはＸ線のエネルギー条件により異なる。したがって密度検査領域のサイズは、密度分析時における上記のようなパラメータを考慮して決定することになる。一の例では、撮像データが２次元の場合、密度分析対象部位のサイズを２×２の４画素ないし３×３の９画素（画素のサイズは例えば２００μ程度）とし、これを中心にして周囲に隣接する状態で数百〜数千画素を密度検査領域として設定する。このような密度検査領域は、密度分析対象部位の設定間隔によっては互いに重なり合う状態になることもある。
【００２１】
図４の一次変換後の撮像データ１０１´、１０２´に密度検査領域を設定・配列した例を図７に示す。図の例では複数の密度分析対象部位それぞれの密度検査領域３を同じ大きさにしてある。このように各密度検査領域３を同じ大きさに設定してあれば、撮像データ１０１´、１０２´それぞれの同一座標位置に密度検査領域３が配列されることになり、密度検査領域３内の各画素におけるＣＴ値を平均化した値の撮像データ１０１´、１０２´間での差から撮像データ１０１と撮像データ１０２の間での密度変化量を２次元断面画像または３次元形状画像上での分布として求めることができる。
【００２２】
ここで、以上の例では被検体について位置合わせを行った後に密度検査領域３を設定・配列するようにしていたが、密度検査領域３の設定・配列を予め決めておき、その設定・配列を一次変換により撮像データ１０１と撮像データ１０２のそれぞれに適合性させる方法とすることも可能である。図８に、予め設定・配列を決めておいた密度検査領域３を撮像データ１０１と撮像データ１０２のそれぞれに適合させる場合のイメージを示す。ＣＴ値は、デジタルデータとして表現されているため、画像を一次変換する場合には、有限桁数で表現されていることによる打切り誤差を生じ、撮像データ１０１と撮像データ１０１´あるいは撮像データ１０２と撮像データ１０２´の間で若干の誤差を含む。この誤差は微小な密度変化に対して影響を与える可能性がある。こうしたことを考慮すると、予め決めてある設定・配列の密度検査領域を撮像データに適合させる手法は、密度分析の高精度化に対し、より有効であるといえる。
【００２３】
撮像データが３次元の場合には、密度検査領域も３次元となる。図９に３次元撮像データにおける３次元の密度検査領域３の設定・配置の例を示す。このように３次元撮像データに３次元の密度検査領域３を設定・配置することで、密度変化の３次元分布（空間分布）を得ることができる。
【００２４】
図１０には密度検査領域の形状についての他の例を示す。（ａ）は円弧形状にした密度検査領域３ａの例であり、（ｂ）は部分扇形状にした密度検査領域３ｂの例であり、（ｃ）は球形にした密度検査領域３ｃの例である。このように密度検査領域は、任意の形状で設定することが可能である。
【００２５】
密度検査領域の設定を終えたら各密度検査領域３について平均密度を算出する（処理１００６）。平均密度は以下の式により計算する。
【００２６】
【数１】

ここで、ρ_ｉバーは撮像データ１０１´または撮像データ１０２´に設定した個々の密度検査領域における平均密度、ｉ（＝１，２…，Ｍ）は撮像データ１０１´、１０２´それぞれに設定した密度検査領域の数、ρ_ｉｊは各密度検査領域に含まれる個々の画素におけるＣＴ値を密度に変換した値、ｊ（＝１，２…，Ｎ）は各密度検査領域に含まれる画素の数である。
【００２７】
ＣＴ値から密度値への変換は以下の式による。
【００２８】
【数２】

ここで、ρは密度、ｘはＣＴ値であり、ａとｂは以下の式により求められる。
【００２９】
【数３】

【００３０】
【数４】

（３）式と（４）式において〈〉で囲まれた値は平均値を表す。なお必要に応じて分散あるいは標準偏差を計算することもできる。分散と標準偏差は下記の（５）、（６）式で計算する。
【００３１】
【数５】

【００３２】
【数６】

【００３３】
平均密度の算出を終えたら、撮像データ１０１´と撮像データ１０２´それぞれの対応する密度検査領域間で密度差を求める（処理１００７）。密度差は以下の式により計算する。
【００３４】
【数７】

ここで、（ρ_ｉバー）_Ａは撮像データ１０１´に設定した密度検査領域の内のｉ番目の密度検査領域における平均密度、（ρ_ｉバー）_Ｂは撮像データ１０２´に設定した密度検査領域の内のｉ番目の密度検査領域における平均密度、Δρ_ｉバーはこれらの密度差である。
【００３５】
以上の処理を撮像データ１０１´、１０２´のそれぞれに設定した複数の密度検査領域の全てについて行うことで、被検体１に対する局所的な密度の時間変化量（処理前後の変化量）を求めることができる。得られた密度の変化量を各密度検査領域の中心座標における密度値（代表密度値）として、その密度値を等高線で結ぶことにより、密度変化量の２次元あるいは３次元分布図を作成することができる。密度変化量の２次元分布図の例を図１１に示す。また、各密度検査領域における密度変化量を複数の密度検査領域全体にわたって積分することにより、被検体の重量の時間的変化量を算出することもできる。
【００３６】
以上の例では、位置合わせした後の撮像データ１０１´、１０２´に設定した密度検査領域３について平均密度を求め、それについて差分を求めるようにしていたが、これに代えて、位置合わせした後の撮像データ１０１´、１０２´それぞれの対応する画素ごとに密度差を求め、その密度差について密度検査領域３ごとに平均値を求めるようにすることも可能である。また以上の例は、撮像データが２つの場合であったが、これに限られず、撮像データを３つ以上取得する場合にも同様な処理が可能である。
【００３７】
以上のように本発明では、密度分析対象部位を中心にして適切な数の画素を含む密度検査領域を設定し、この密度検査領域について、そこに含まれる各画素のＣＴ値から求まる各密度値から平均値を求め、これを密度分析対象部位における代表密度値とするようにしている。このようにすることにより、画素ごとにランダムに生じるノイズの影響をキャンセルして大幅に低減することができる。この結果、密度分析の精度を大幅に高めることができ、Ｘ線ＣＴ装置における分解能の近辺でも有意な密度分析が可能となる。
【００３８】
こうした本発明による密度分析方法は、以上の実施形態で説明したような同一の被検体における処理前後での密度変化の分析のように分析対象の密度差が微小である場合に特に有効であるが、一般的な密度分析についてもその分析の高精度性は有用である。
【００３９】
図１２に第２の実施形態による密度分析方法における処理の流れを示す。本実施形態は、量産される製品を被検体として、各製品における密度分布の差異を分析することで製品のばらつきを評価する場合などに適用することを前提にした密度分析に関するものである。処理２００１で、図２における構成においてＸ線ＣＴ装置２により被検体１の撮像データ２０１を取得する。続く処理２００２では被検体１と同一形状の被検体１ａ（被検体１と同一の製品）について撮像データ２０２を取得する。以下の処理２００３〜２００６は図１における処理１００４〜１００７に対応しており、それぞれにおける処理は図１に関して説明したのと同様である。
【００４０】
以下では本発明による密度分析方法を実装したコンピュータプログラムの実施形態について説明する。一実施形態によるコンピュータプログラムを実行させた場合にコンピュータの表示装置に表示される画面の例を図１３〜図１９に示す。本実施形態は、上記第１の実施形態による密度分析方法を実装したコンピュータプログラムの場合である。コンピュータプログラムは、撮像データ読み込み機能を含むとともに、図１に関して上で説明した各処理についての機能、すなわち位置合わせのための機能（位置合わせ機能）、密度検査領域を設定するための機能（密度検査領域設定機能）、密度検査領域ごとに平均密度を求める機能（密度算出機能）、および複数の撮像データ間で密度差を求める機能（密度差算出機能）を含んでいる。以下ではこれらの各機能について説明する。
【００４１】
図１３は、図１における処理１００１と処理１００３で得られた撮像データを撮像データ読み込み機能により読み込む処理過程で表示される画面１４の例であり、Ｘ線ＣＴ装置２による撮像で得られた撮像データ１０１と撮像データ１０２をコンピュータに読み込んでそれぞれの画像を表示している状態を示している。
【００４２】
図１４は、位置合わせ機能を呼び出して図１における処理１００４を行っている状態での画面の例を示している。すなわち位置合わせのための設定、および位置合わせ後に撮像データを一次変換した場合の画面の遷移を示している。位置合わせ機能には、被検体または被検体に取り付けた位置合せ用指標体が有する特徴的な幾何学的形状を計測する手段およびそれに基づいた位置合わせ処理を行う手段が含まれている。図の例ではそれらの手段をツール１５としてまとめて画面に表示している。また図の例では、直線、円、円弧、平面、円筒面などの幾何学的形状を計測可能としている。この他に、円錐や球なども計測できるようにするようにする場合もある。位置合わせ機能による処理の内容は上で説明したのと同様なので、ここでは説明を省略する。
【００４３】
図１５は、密度検査領域設定機能（密度検査領域配置機能）により撮像データ１０１´と撮像データ１０２´に密度検査領域を設定する処理を行っている状態での画面の遷移を示す。密度検査領域設定機能は、密度検査領域の形状を決定する手段、密度検査領域のサイズを決定する手段、および密度検査領域の配列を決定する手段を含んでおり、それらはツール１６としてまとめて画面に表示される。図の例では、密度検査領域の形状として、直線、円弧、任意の曲線、矩形、円、多角形、直方体、球などを指定可能としてある。密度検査領域設定機能における密度検査領域配列決定手段には、画面１７を参照しながら配列を決定する方式や数式により配列を決定する方式を適用することができる。また密度検査領域配列決定手段には、画面１７上で配列した密度検査領域を撮像データ１０１´、１０２´に反映させて配列状態を確認するためのプレビュー機能を含ませてある。密度検査領域設定機能における密度検査領域サイズ決定手段には、画面１７上でマウスなどの入力手段により密度検査領域のサイズを決定する方式や密度検査領域のサイズを数値で入力する方式を適用することができる。一旦決定した密度検査領域サイズを変更することもできる。図１６に示すのは、図１５の状態における密度検査領域のサイズを拡大修正した場合の画面の例である。
【００４４】
位置合わせ時に密度検査領域の一次変換を選択することもでき、それを選択した場合にも同様の機能を利用できる。ツール１６における画面１７に密度検査領域のサイズや配列を指定すると、撮像データ１０１、１０２上に、位置合わせ時に取得した結果に応じた一次変換を密度検査領域に対して施した状態が表示される。この処理の様子を図１７に示す。
【００４５】
図１８は、密度算出機能（密度計測機能）により密度の算出処理を行っている状態での画面の遷移を示している。密度算出機能は、密度算出対象の撮像データを選択する手段、密度の算出を実行する手段、および算出結果を表示する手段を含んでおり、それらはツール１８としてまとめて画面に表示される。
【００４６】
図１９は、密度差算出機能（密度差計算機能）により密度差の算出処理を行っている状態での画面の遷移を示している。密度差算出機能は、密度差算出対象の撮像データを選択する手段、密度差の算出を実行する手段、および算出結果を表示する手段を含んでおり、それらはツール１９としてまとめて画面に表示される。図の例では、密度差算出結果として密度差分布図２０も表示するようにしてある。
【００４７】
ここで、以上の実施形態では、撮像データ１０１と撮像データ１０２に対して同等に行う処理を両撮像データについて同時的に施すようにしているが、各撮像データに順次的に施すような形態とすることもできる。
【００４８】
以上のようなコンピュータプログラムを用いることにより、本発明による密度分析方法をより簡易にかつ効率的に進めることが可能となる。
【００４９】
次に、本発明による密度分析方法を実施する上で好ましい密度分析システムの形態例について説明する。図２０に第１の実施形態による密度分析システムの構成を示す。本実施形態における密度分析システムは、被検体を撮像するＸ線ＣＴ装置２、Ｘ線ＣＴ装置２で撮像された画像（撮像データ）を受信する画像受信装置（撮像データ受信装置）２１、受信した画像を記憶しておく画像記憶装置（撮像データ記憶装置）２２、画像記憶装置２２に記憶されている画像を呼び出して密度分析を実行する密度分析装置２３、密度分析装置２３による分析結果を記憶しておく密度分析結果記憶装置２４、画像や密度分析過程を表示する表示装置２５を備えて構成される。
【００５０】
図２１に示すのは、第２の実施形態による密度分析システムの構成である。本実施形態における密度分析システムは、第１の実施形態におけるそれと基本的には同様である。相違しているのは、密度分析装置２３に相当する要素として、密度分析プログラム記憶装置２６と演算装置（コンピュータ）２７からなる装置を備えていることである。すなわち本実施形態の密度分析システムでは、密度分析プログラム記憶装置２６と演算装置２７が「密度分析装置」を構成し、密度分析プログラム記憶装置２６に記憶させてある密度分析プログラムを演算装置２７で実行させることで密度分析を行うことになる。その密度分析プログラムには、上で説明したコンピュータプログラムが用いられる。
【産業上の利用可能性】
【００５１】
本発明は、Ｘ線ＣＴ装置による撮像データを用いる密度分析について、分析精度を大幅に高め、Ｘ線ＣＴ装置における分解能の近辺でも有意な密度分析が可能とするものであり、さまざまな分野で求められる密度分析に広く適用することができる。
【図面の簡単な説明】
【００５２】
【図１】第１の実施形態による密度分析方法における処理の流れを示す図である。
【図２】Ｘ線ＣＴ装置による被検体の撮像の様子を模式化して示す図である。
【図３】幾何学的形状に基づく位置合わせの例のイメージを示す図である。
【図４】一次変換後の撮像データの例を示す図である。
【図５】位置合せ用指標体を被検体に取り付けて撮像した撮像データの例を示す図である。
【図６】３次元形状撮像データについて位置合わせの例のイメージを示す図である。
【図７】図４の一次変換後の撮像データに密度検査領域を設定・配列した例を示す図である。
【図８】予め設定・配列を決めたおいた密度検査領域を撮像データに適合させる場合のイメージを示す図である。
【図９】３次元撮像データにおける３次元の密度検査領域の設定・配置の例を示す図である。
【図１０】密度検査領域の形状についての他の例を示す図である。
【図１１】密度変化量の２次元分布図の例を示す図である。
【図１２】第２の実施形態による密度分析方法における処理の流れを示す図である。
【図１３】撮像データを撮像データ読み込み機能により読み込む処理過程で表示される画面の例を示す図である。
【図１４】位置合わせ機能により位置合わせ処理を行っている状態での画面の遷移を示す図である。
【図１５】密度検査領域設定機能により密度検査領域を設定する処理を行っている状態での画面の遷移を示す図である。
【図１６】図１５の状態における密度検査領域のサイズを拡大修正した場合の画面の例を示す図である。
【図１７】位置合わせ時に密度検査領域の一次変換を選択した場合の処理の様子を示す図である。
【図１８】密度算出機能により密度の算出処理を行っている状態での画面の遷移を示す図である。
【図１９】密度差算出機能により密度差の算出処理を行っている状態での画面の遷移を示す図である。
【図２０】第１の実施形態による密度分析システムの構成を示す図である。
【図２１】第２の実施形態による密度分析システムの構成を示す図である。
【符号の説明】
【００５３】
１被検体
２Ｘ線ＣＴ装置
３密度検査領域
１０１撮像データ
１０２撮像データ

【特許請求の範囲】
【請求項１】
Ｘ線ＣＴ装置で被検体を撮像して得られる撮像データに基づいて前記被検体の密度を分析する密度分析方法において、
前記撮像データに設定される複数の密度分析対象部位について、前記密度分析対象部位ごとにそれを中心にして複数の画素からなる密度検査領域を設定し、前記密度検査領域における複数の画素の各密度値の平均値を求め、求めた平均値を前記密度分析対象部位の密度値として密度分析を行うようにしたことを特徴とする密度分析方法。
【請求項２】
Ｘ線ＣＴ装置で被検体を撮像して得られる撮像データに基づいて前記被検体の密度を分析する密度分析方法において、
前記Ｘ線ＣＴ装置により同一の被検体について少なくとも２つ以上の撮像データをそれぞれ異なる時刻に取得するか、または形状、材質が同一の複数の被検体のそれぞれについて撮像データを取得する処理過程、前記撮像データ取得過程で得られた複数の撮像データについて位置合わせを行う処理過程、前記複数の撮像データのそれぞれに設定される複数の密度分析対象部位について、前記密度分析対象部位ごとにそれを中心にして複数の画素からなる密度検査領域を設定する処理過程、前記処理過程で設定された密度検査領域における複数の画素の各密度値の平均値を求める処理過程、および前記処理過程で求められた密度値の平均値を前記密度分析対象部位の密度値とし、前記複数の撮像データ間での対応密度分析対象部位間について密度差を求める処理過程を含むことを特徴とする密度分析方法。
【請求項３】
請求項１または請求項２に記載の密度分析方法を実装したコンピュータプログラムであって、
前記Ｘ線ＣＴ装置により取得された撮像データを読み込む撮像データ読み込み機能、前記撮像データ読み込み機能で読み込まれた複数の撮像データについて位置合わせを行うための位置合わせ機能、前記複数の撮像データのそれぞれに設定される複数の密度分析対象部位について、前記密度分析対象部位ごとにそれを中心にして複数の画素からなる密度検査領域を設定するための密度検査領域設定機能、前記密度検査領域設定機能で設定された密度検査領域における複数の画素の各密度値の平均値を求めるための密度算出機能、および前記密度算出機能で求められた密度値の平均値を前記密度分析対象部位の密度値とし、前記複数の撮像データ間での対応密度分析対象部位間について密度差を求めるための密度差算出機能を含んでいることを特徴とするコンピュータプログラム。
【請求項４】
Ｘ線ＣＴ装置で被検体を撮像して得られる撮像データに基づいて前記被検体の密度に関して分析する密度分析方法に用いられる密度分析システムにおいて、
請求項１または請求項２に記載の密度分析方法を実行するための密度分析装置を備えていることを特徴とする密度分析システム。
【請求項５】
前記密度分析装置は、請求項３に記載のコンピュータプログラムを記憶するプログラム記憶装置と、前記プログラム記憶装置から読み出した前記コンピュータプログラムを実行する演算装置とで構成される請求項４に記載の密度分析システム。

【図１】