(57)【特許請求の範囲】
【請求項１】
骨部組織の構造の状態を表す指標値を取得する骨計測方法であって、
少なくとも前記骨部組織を含む被検体を被写体とする放射線画像に対してモーフォロジー演算に基づく骨部組織の強調処理を施し、
前記強調処理された骨部組織の構造の状態を表す骨部組織画像にスターボリューム法に基づいた処理を施して前記指標値を取得することを特徴とする骨計測方法。
【請求項２】
骨部組織の構造の状態を表す指標値を取得する骨計測方法であって、
少なくとも前記骨部組織を含む被検体を被写体とする放射線画像に対してモーフォロジー演算に基づく骨部組織の強調処理を施し、
前記強調処理された骨部組織の構造の状態を表す骨部組織画像にノードストラット解析法に基づいた処理を施して前記指標値を取得することを特徴とする骨計測方法。
【請求項３】
前記強調処理が、スケルトン処理であることを特徴とする請求項１または２記載の骨計測方法。
【請求項４】
骨部組織の構造の状態を表す指標値を算出する骨構造指標値算出手段を備えた骨計測装置であって、
少なくとも前記骨部組織を含む被検体を被写体とする放射線画像に対してモーフォロジー演算に基づく骨部組織の強調処理を施すモーフォロジー演算手段をさらに備え、
前記骨構造指標値算出手段が、前記強調処理された骨部組織の構造を表す骨部組織画像にスターボリューム法に基づいた処理を施して前記指標値を取得するものであることを特徴とする骨計測装置。
【請求項５】
骨部組織の構造の状態を表す指標値を算出する骨構造指標値算出手段を備えた骨計測装置であって、
少なくとも前記骨部組織を含む被検体を被写体とする放射線画像に対してモーフォロジー演算に基づく骨部組織の強調処理を施すモーフォロジー演算手段をさらに備え、
前記骨構造指標値算出手段が、前記強調処理された骨部組織の構造を表す骨部組織画像にノードストラット解析法に基づいた処理を施して前記指標値を取得するものであることを特徴とする骨計測装置。
【請求項６】
前記モーフォロジー演算手段による強調処理が、スケルトン処理であることを特徴とする請求項４または５記載の骨計測装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は骨計測方法および装置に関し、詳しくは骨粗鬆症等の診断に有用な人体等の骨部組織の構造を、放射線画像に基づいて、定量的な指標値として求める方法および装置に関するものである。
【０００２】
【従来の技術】
骨塩定量、すなわち骨の中のカルシウムの量を定量的に測定することは骨折予防の診断のために有用である。
【０００３】
つまり、骨塩量は、骨の内部を構成する海綿質である骨梁の疎密、すなわち骨密度によってに定まるため、骨密度が疎であれば骨部画像における骨部の陰影の濃度が高くなり、骨密度が密であれば骨部画像における骨部の陰影の濃度は低くなる。
【０００４】
したがって、骨中のカルシウムの微量変化を知ることは骨粗鬆症の早期発見を可能にし、骨折予防のうえでも重要な効果がある。
【０００５】
そこで従来、ＭＤ法(Microdensitometry) 、ＳＰＡ法(Single Photon Absorptiometry)、ＤＰＡ法(Dual Photon Absorptiometry)、ＱＤＲ法(Quantitative Digited Radiography)、ＱＣＴ法(Quantitative Computer Tomography)、ＤＱＣＴ法(Dual energy Quantitative Computer Tomography)等、数々の骨塩定量の方法が提案され、実施されている。
【０００６】
ところで、これらの各手法はいわゆる骨密度を計量するものであり、確かに骨粗鬆症の診断の指標値として従来から広く一般に用いられているが、最近は、この骨密度値の他に、骨構造状態を指標値として定量的に表すことが提案されている。
【０００７】
すなわち、骨密度値は、単純には体積に対する質量の比の値であるから、荷重を支持する骨としての実体である骨梁の分布状態（骨構造）の差異が反映されるものではないが、骨粗鬆症は本来は骨の強度の問題であって、この強度は骨梁の分布状態に大きく依存するものであるから、骨構造状態を指標値化することは、骨粗鬆症の診断においては極めて有用な技術である。
【０００８】
このような骨構造状態を表す指標値としては、例えば、スターボリューム（Star volume ）、ノードストラット解析（Node-strut analysis ）により得られる指標置などが知られている。
【０００９】
なお、本明細書においては、説明の冗長を避けるため、指標値であるスターボリュームを算出する処理方法をスターボリューム法と、ノードストラット解析により指標値を算出する処理方法をノードストラット解析法と、それぞれ称するものとする。
【００１０】
スターボリューム法では、図３に示すように、骨髄腔内の所定の点ａから全ての方向について骨梁に遮られることなく見渡せる範囲の骨髄腔体積の平均を意味するＶｍ（marrow space star volume）、およびこれと同様に骨梁の中の点ｂから全ての方向に骨梁の端に至る範囲の骨梁体積の平均を意味するＶｔ（trabecular star volume）の２種類の指標値が定義されている。なお、図３においては斜線部が骨梁を、その他の部分が骨髄腔を示す。Star volume は、標本のサンプリング法の工夫により骨髄腔の骨梁の大きさをｍｍ3 あるいはμｍ3 といった３次元の値として表す隔たりのない立体学的指標とされている。Ｖｍは骨梁の連続性が高い場合は小さく、骨梁の消滅や穿孔が多いときは大きくなり、逆にＶｔは骨梁の連続性が高い場合は大きく、骨梁の消滅や穿孔が多いときは小さくなる。
【００１１】
骨髄腔内の任意の点ｉにおけるＶｍｉは以下の式（１）により定義される。
【００１２】
Ｖｍｉ＝（π／３）×ｌ0m3           （１）
但し、点ｉを中心とする任意の方向において、骨髄腔が連続する長さをｌ0 としたとき、ｌ0m3 は点ｉを中心として放射状に全ての方向においてｌ03を求めた場合の平均値を表す。
【００１３】
また、骨梁上の任意の点ｊにおけるＶｔｊは以下の式（２）により定義される。
【００１４】
Ｖｔｊ＝（π／３）×Σｌ14／Σｌ1   （２）
但し、ｌ1 は点ｊを中心とする任意の方向において、骨梁が連続する長さを示す。また、ｌ1 は点ｊを中心として放射状に全ての方向において求められ、Σはその全ての方向における和を求めることを表す。そして、ＶｍｉおよびＶｔｊをサンプリング点ごとに算出し、その平均値をＶｍおよびＶｔとすればよい。
【００１５】
なお、骨梁構造の指標値としては、ＶｍとＶｔのいずれを用いてもよい。
【００１６】
そして例えば図４に示すように、被検体の年齢、骨密度値およびＶｔを軸とする三次元座標系が形成され、入力された結果がこの三次元座標内のいずれの位置に存在するかに応じて、骨部組織の状態を判定する。すなわち、被検体の年齢、骨密度値およびＶｔが点Ｐ１で示す位置にあれば状態Ａ（例えば、骨粗鬆症の疑いがある状態）にあり、点Ｐ２で示す位置にあれば状態Ｂ（注意を要する状態）にあり、点Ｐ３で示す位置にあれば状態Ｃ（正常な状態）にある、と判定するものである。
【００１７】
一方、ノードストラット解析法は、骨梁の連続性を二次元的に評価する方法であり、３個以上の骨梁の結合点をＮｄ（Node；ノード、図５における白点）、他の骨梁と結合していない終末点をＴｍ（Terminus、図５における黒点）と定義する。そして、図５に示すように、これらをつなぐ骨梁の中心線（strut ；ストラット）を、ＮｄＮｄ（結合点同士をつなぐストラット）、ＮｄＴｍ（結合点と終末点とをつなぐストラット）、ＴｍＴｍ（終末点同士をつなぐストラット）、ＣｔＮｄ（皮質骨（Ｃｔ）と結合点とをつなぐストラット）、およびＣｔＴｍ（皮質骨と終末点とをつなぐストラット）に分類し、各ストラットの長さを測定する。
【００１８】
次いでこれらの各ストラットの長さと結合点Ｎｄの数ＮNdおよびＴｍの数ＮTmに基づいて以下の指標値を定義する。
【００１９】
（１）長さの指標値
すべてのストラットの長さの和ＴSLに対する各ストラットの長さの割合（％）
ＮｄＮｄ／ＴSL
ＮｄＴｍ／ＴSL
ＴｍＴｍ／ＴSL
ＣｔＮｄ／ＴSL
ＣｔＴｍ／ＴSL
（２）面積に対する長さの指標値
骨部組織（骨梁および皮質骨の和、または骨梁のみ）の面積ＴＶに対するストラットの長さの割合（ｍｍ／ｍｍ2 ）
ＴSL／ＴＶ
ＮｄＮｄ／ＴＶ
ＮｄＴｍ／ＴＶ
ＴｍＴｍ／ＴＶ
ＣｔＮｄ／ＴＶ
ＣｔＴｍ／ＴＶ
（３）数の指標値
ＮNd／ＴＶ；骨部組織の面積当たりの結合点の数（／ｍｍ2 ）
ＮTm／ＴＶ；骨部組織の面積当たりの終末点の数（／ｍｍ2 ）
ＮNd／ＮTm；結合点の数と終末点の数との比
以上の指標値は、結合点Ｎｄ関連の指標値が大きいほど骨梁の連続性が高く、終末点Ｔｍの関連のパラメータが大きいほど骨梁の連続性が低いことを示す。
【００２０】
このノードストラット解析法は、骨梁の連続性を直接かつ簡便に計測できる方法として有用性が高い。
【００２１】
【発明が解決しようとする課題】
しかし、上述した骨構造の状態を表す指標値を求める方法は、骨のスライス試片という実体の標本を対象とするものであるため、被検体から非侵襲的に求めることができないという、実用上の問題があった。
【００２２】
一方、上記骨部組織を含む放射線画像に基づいて上記骨部組織の構造状態を指標値化する試みもあるが、放射線画像における、軟部組織との重なり、濃度・コントラストの変動等の要因のため、実用上の計測精度を確保することができないという問題がある。
【００２３】
本発明は、上記事情に鑑みなされたものであって、骨部組織の構造状態を、放射線画像に基づいて精度よく指標値化することを可能にした骨計測方法および装置を提供することを目的とするものである。
【００２４】
【課題を解決するための手段】
本発明の骨計測方法は、放射線画像を表す画像情報に対してモーフォロジー演算に基づく骨部組織の強調処理を施すことにより、骨部組織の構造を抽出または強調処理し、その抽出または強調処理された骨部組織の構造を表す画像に基づいて、当該構造状態を表す指標値を算出するものである。
【００２５】
すなわち本発明の骨計測方法は、骨部組織の構造の状態を表す指標値を取得する骨計測方法であって、
少なくとも前記骨部組織を含む被検体を被写体とする放射線画像に対してモーフォロジー演算に基づく骨部組織の強調処理を施し、
前記強調処理された骨部組織の構造を表す骨部組織画像に対して、前記指標値を取得する処理を施すことを特徴とするものである。
【００２６】
ここで、骨部組織とは、骨梁とそれ以外の空洞部分である骨髄腔、およびこれらの外周を覆う皮質骨等を意味し、骨部組織の構造の状態とは主として骨梁と骨髄腔との分布状態を意味する。
【００２７】
骨部組織の構造の状態を表す指標値としては、前述したスターボリューム法によるＶｍ若しくはＶｔ、またはノードストラット解析法によるＮｄＮｄ／ＴSL、ＮｄＴｍ／ＴSL、ＴｍＴｍ／ＴSL、ＣｔＮｄ／ＴSL、ＣｔＴｍ／ＴSL、ＴSL／ＴＶ、ＮｄＮｄ／ＴＶ、ＮｄＴｍ／ＴＶ、ＴｍＴｍ／ＴＶ、ＣｔＮｄ／ＴＶ、ＣｔＴｍ／ＴＶ、ＮNd／ＴＶ、ＮTm／ＴＶ、ＮNd／ＮTm等を適用することができる。
【００２８】
また、モーフォロジー演算に基づく骨部組織の強調処理としては、スケルトン処理を適用することができる。ここでモーフォロジー演算について説明する。
【００２９】
モーフォロジー演算（以下、モーフォロジー処理ともいう）とは、原画像のうち異常な陰影等の特定の画像部分だけを選択的に抽出する、モーフォロジー（Morphology；モフォロジーまたはモルフォロジーとも称する）のアルゴリズムに基づく処理であり、特に乳癌における特徴的形態である微小石灰化像を検出するのに有効な手法として研究されているが、対象画像としてはこのようなマンモグラムにおける微小石灰化像に限るものではない。
【００３０】
そしてこのモーフォロジー処理は、抽出しようとする画像部分の大きさ、形状に対応した構造要素Ｂを用いた処理を行なうものであり、複雑なバックグラウンド情報に影響されにくい、抽出した画像が歪まない、などの特徴がある。
【００３１】
すなわち、この手法は一般の微分処理に比べて、石灰化像のサイズ・形状・濃度分布などの幾何学的情報をよりよく保って検出することができる。
【００３２】
以下、このモーフォロジー処理の概要を、マンモグラムにおける微小石灰化像の検出に適用した例について説明する。
【００３３】
（モーフォロジーの基本演算）
モーフォロジー処理は一般的にはＮ次元空間における集合論として展開されるが、直感的な理解のために２次元の濃淡画像を対象として説明する。
【００３４】
濃淡画像を座標（ｘ，ｙ）の点が濃度値ｆ（ｘ，ｙ）に相当する高さをもつ空間とみなす。ここで、濃度値ｆ（ｘ，ｙ）は、濃度が低い（ＣＲＴに表示した場合には輝度が高い）程大きな画像信号値となる高輝度高信号レベルの信号とする。
【００３５】
まず、簡単のため上記２次元の濃淡画像の断面に相当する１次元の関数ｆ（ｘ）を考える。モーフォロジー演算に用いる構造要素ｇは次式（３）に示すように、原点について対称な対称関数
【００３６】
【数１】
0003656695.tif 000002
【００３７】
であり、定義域内で値が０で、その定義域Ｇが下記式（４）であるとする。
【００３８】
【数２】
0003656695.tif 000003
【００３９】
このとき、モーフォロジー演算の基本形は式（５）〜（８）に示すように、非常に簡単な演算となる。
【００４０】
【数３】
0003656695.tif 000004
【００４１】
すなわち、ダイレーション（dilation）処理は、注目画素を中心とした、±ｍ（構造要素Ｂに応じて決定される値であって、図６中のマスクサイズに相当）の幅の範囲内の最大値を探索する処理であり（同図（Ａ）参照）、一方、エロージョン（erosion ）処理は、注目画素を中心とした、±ｍの幅の範囲内の最小値を探索する処理である（同図（Ｂ）参照）。また、オープニング（opening ）処理はエロージョン処理後にダイレーション処理を行なう処理、すなわち最小値の探索の後に最大値を探索する処理であり、クロージング（closing ）処理は、ダイレーション処理後にエロージョン処理を行なう処理、すなわち最大値の探索の後に最小値を探索する処理に相当する。
【００４２】
つまりオープニング処理は、低輝度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより空間的に狭い範囲で変動する凸状の濃度変動部分（周囲部分よりも輝度が高い部分）を取り除くことに相当する（同図（Ｃ）参照）。
【００４３】
一方、クロージング処理は、高輝度側から濃度曲線ｆ（ｘ）を滑らかにし、マスクサイズ２ｍより空間的に狭い範囲で変動する凹状の濃度変動部分（周囲部分よりも輝度が低い部分）を取り除くことに相当する（同図（Ｄ）参照）。
【００４４】
なお、構造要素ｇが原点に対して対称ではない場合の、式（５）に示すダイレーション演算をミンコフスキー（Minkowski ）和、式（６）に示すエロージョン演算をミンコフスキー差という。
【００４５】
ここで、濃度の高いもの程大きな値となる高濃度高信号レベルの信号の場合においては、濃度値ｆ（ｘ）の画像信号値が高輝度高信号レベルの場合に対して大小関係が逆転するため、高濃度高信号レベルの信号に対するダイレーション処理と高輝度高信号レベルに対するエロージョン処理（同図（Ｂ））とは一致し、高濃度高信号レベルの信号に対するエロージョン処理と高輝度高信号レベルに対するダイレーション処理（同図（Ａ））とは一致し、高濃度高信号レベルの信号に対するオープニング処理と高輝度高信号レベルに対するクロージング処理（同図（Ｄ））とは一致し、高濃度高信号レベルの信号に対するクロージング処理と高輝度高信号レベルに対するオープニング処理（同図（Ｃ））とは一致する。
【００４６】
なお、本項では高輝度高信号レベルの画像信号（輝度値）の場合について説明する。
【００４７】
（石灰化陰影検出への応用）
石灰化陰影の検出には、原画像から平滑化した画像を引き去る差分法が考えられる。単純な平滑化法では石灰化陰影と細長い形状の非石灰化陰影（乳腺、血管および乳腺支持組織等）との識別が困難であるため、東京農工大の小畑らは、多重構造要素を用いたオープニング処理に基づく下記式（９）で表されるモーフォロジーフィルターを提案している（「多重構造要素を用いたモルフォロジーフィルタによる微小石灰化像の抽出」電子情報通信学会論文誌 D-II Vol.J75-D-II No.7 P1170 〜1176  1992年７月、「モルフォロジーの基礎とそのマンモグラム処理への応用」MEDICAL IMAGING TECHNOLOGY Vol.12 No.1 January 1994 等）。
【００４８】
【数４】
0003656695.tif 000005
【００４９】
ここでＢi （ｉ＝１，２，…，ｎ）は、直線状の大きさがｍ画素でｎ個（例えば図７に示すものでは、９画素４方向でありｍ＝９，ｎ＝４）の構造要素（これらを全体として以下、ｍ画素ｎ方向の多重構造要素という）である。構造要素Ｂi を検出対象である石灰化陰影よりも大きく設定すれば、上記オープニング演算による処理で、構造要素Ｂi よりも細かな信号変化部分（空間的に狭い範囲で信号が変動する画像部分）であって周囲よりも輝度値の大きい凸状の部分である石灰化陰影は取り除かれる。一方、細長い形状の乳腺の陰影等の非石灰化陰影はその長さが構造要素Ｂi よりも長く、その傾きが４つの構造要素Ｂi のいずれかに一致すればオープニング処理（式（９）の第２項の演算）をしてもそのまま残る。したがってオープニング処理によって得られた平滑化画像（石灰化陰影のみが取り除かれた画像）を原画像ｆから引き去ることで、小さな石灰化陰影のみが含まれる画像が得られる。これが式（９）の考え方である。
【００５０】
なお、前述したように、高濃度高信号レベルの信号の場合においては、石灰化陰影は周囲の画像部分よりも濃度値が低くなり、石灰化陰影は周囲部分に対して濃度値の小さい凹状の信号変化部分となるため、オープニング処理に代えてクロージング処理を適用し、式（９）に代えて式（10）を適用する。
【００５１】
【数５】
0003656695.tif 000006
【００５２】
なお、モーフォロジー演算の一例である式（10）のクロージング処理を具体的に説明する。
【００５３】
すなわち、高濃度高信号レベルの画像信号である濃度値Ｓorg についてのモーフォロジー演算によれば、例えば図８（１）の実線に示すような濃度値Ｓorg の分布を有する画像データに対して、同図（２）に示すような直線状の３画素の構造要素Ｂで、最大値処理（ダイレーション処理）を行うことにより、ある注目画素の濃度値Ｓi は、その注目画素を中心として互いに隣接する３画素（構造要素Ｂにより決定される）の中の最大値Ｓi+1 を採用したＳi ′に変換される。この演算を全画素について行うことにより、濃度値Ｓorg ′の分布を有する同図（１）の破線で示す最大値信号に変換される。
【００５４】
次に、この最大値処理により得られた最大値信号に対してさらに構造要素Ｂによる最小値処理（エロージョン処理）を考えると、同図（１）の破線で示された注目画素の最大値信号Ｓi ′は、その注目画素を中心として互いに隣接する３画素の中の最小値Ｓi-1 ′を採用したＳi ″（＝Ｓi ）に変換される。この演算を全画素について行うことにより、最大値処理後の最小値信号Ｓorg ″の分布は同図（１）の一点鎖線で示すものとされる。この一点鎖線で示された画像信号は、もとの実線のオリジナルの画像データに対して、構造要素Ｂよりも空間的に狭い範囲で信号が変動する画像部分が消え、構造要素Ｂよりも空間的に広い範囲で変動する信号値の変化部分である画像部分や変動のない画像部分はもとの形状のまま残っていることを示している。すなわち、以上の処理（クロージング処理）は、画像濃度の分布を高濃度側から平滑化する処理として作用する。
【００５５】
このようにクロージング処理で得られた値（Ｓorg に対して最大値処理を行なった後にさらに最小値処理を行なった値）を原画像信号Ｓorg から差し引くことにより得られた値Ｓmor は、上記クロージング処理で消された空間的に狭い範囲で変動する信号値の変化部分である画像部分を表す。
【００５６】
ここで、本来、画像信号は２次元の要素である位置（ｘ，ｙ）と、３次元目の要素である信号値ｆ（ｘ，ｙ）を有するが、上記説明においては、理解の容易化のために、この２次元上に展開された画像の所定の断面に現れた、１次元状の画像信号分布曲線について説明した。
【００５７】
したがって実際には、以上の説明を２次元画像に適用する必要があり、多重構造要素を用いるのも２次元画像に対応させるためである。
【００５８】
次にモーフォロジー演算に基づくスケルトン処理について説明する。
【００５９】
スケルトン処理とは一般に図形の骨格（skeleton）を抽出する処理であり、骨格とは図形に内接する円盤の中心の集合としてとらえることができる。すなわち例えば図９（Ａ）〜（Ｅ）に示す各図形（中太の実線で表記）の骨格はそれぞれ太い実線で示すものとなる。
【００６０】
以下、このスケルトン処理を、上述したモーフォロジー演算により行なう場合について説明する。この場合、スケルトン処理は下記式（11）または（12）で表すことができる。
【００６１】
【数６】
0003656695.tif 000007
【００６２】
【数７】
0003656695.tif 000008
【００６３】
ここで式（11）と式（12）とは前述した、画像を高濃度高信号レベルの画像信号として表すか、または高輝度高信号レベルの画像信号として表すかの差異によるものであり、高濃度高信号レベルの画像信号で表された画像から低濃度（高輝度）の画像部分の骨格を抽出する場合は式（11）を適用し、一方、高輝度高信号レベルの画像信号で表された画像から低輝度（高濃度）の画像部分の骨格を抽出する場合は式（12）を適用するものであり、作用自体の実質的な差はない。
【００６４】
例えばネガフイルム（高濃度高信号レベル）上においては骨部は他の画像部分に比べてその濃度は低いものとなり、骨梁の存在する部分は濃度が低く、存在しない部分は濃度が高くなる。したがって、周囲よりも濃度の低い部分となる骨梁を対象としてスケルトン処理を行なうことに該当するため式（11）を適用すればよい。
【００６５】
ここで式（11）における構造要素Ｂを半径ｒの円とし、図10に示す図形に対してスケルトン処理を施した状態を示す。図10に示す図形は、その輪郭よりも外側領域が濃度の高い部分であり、内側が濃度の低い部分である。
【００６６】
この図形に対してまず図示上段では構造要素Ｂによるエロージョン処理が施される。λ＝０（構造要素Ｂによる０回のエロージョン処理）では図形は何ら変化がない。
【００６７】
λ＝１（構造要素Ｂによる１回のエロージョン処理）では図形が構造要素Ｂの半径ｒ分だけ内側に埋め込められる。
【００６８】
λ＝２（構造要素Ｂによる２回のエロージョン処理）では図形の円から突出した部分が完全に消失する。
【００６９】
同様の操作を繰り返すことにより、λ＝Ｎ−１（構造要素ＢによるＮ−１回のエロージョン処理）で図形は半径ｒ以下の円のみとなる。
【００７０】
一方、図示中段は、この構造要素Ｂによる各回（λ＝０，１，２，…，Ｎ−１，Ｎ）のエロージョン処理を施した画像に対して、さらにそれぞれ構造要素Ｂによるオープニング処理を施した図形である。
【００７１】
この図示上段の図形から図示中段の図形を、処理の回数を対応させて差し引いた図形が図示下段である。
【００７２】
λ＝１において、元の図形の円から突出した部分の骨格要素が抽出され、λ＝Ｎ−１において、元の図形の円の骨格要素が抽出されていることが分かる。
【００７３】
このように元の図形に対してエロージョン処理を施し、さらにオープニング処理を施し、処理回数を対応させて差し引き、この結果の和集合を求めたものが式（11）の意味するところである。
【００７４】
式（12）では濃度の高低が式（11）とは逆転した図形から骨格要素を抽出するのに有効であり、元の図形に対してダイレーション処理を施し、さらにクロージング処理を施し、処理回数を対応させて差し引き、この結果の和集合を求めることにより骨格要素を抽出することを意味する。
【００７５】
なお上記式（11）、（12）中の和集合として、ｎ1 の値が比較的大きいもののみの和集合、すなわち例えばｎ1 ＝０，１の場合を除いたλ＝２，３，４，５等の和集合を採用し、この場合に得られる骨格要素を表示することにより、骨梁の変化をより見易くすることができ、好ましい。これは、λが０や１のときは特に、非常に空間周波数の高いノイズ成分も抽出されるため、和集合からこれらのノイズ成分を除外することによって、より読影性能の高い画像を得ることができるからである。
【００７６】
以上がモーフォロジー演算およびスケルトン処理の内容である。
【００７７】
なお、モーフォロジー演算に基づく強調処理の対象としては、骨部組織を含む被検体を被写体とした通常の放射線画像に代えて、エネルギーサブトラクション処理用の、放射線エネルギー分布が互いに異なる２つ以上の放射線画像に基づいて骨部組織が抽出又は強調処理されたエネルギーサブトラクション画像を適用することもできる。
【００７８】
本発明の骨計測装置は、放射線画像を表す画像情報に対してモーフォロジー演算に基づく骨部組織の強調処理を施すことにより、骨部組織の構造を抽出し、その抽出（強調処理）された骨部組織の構造を表す画像に基づいて、当該構造状態を表す指標値を算出するものである。
【００７９】
すなわち本発明の骨計測装置は、骨部組織の構造の状態を表す指標値を算出する骨構造指標値算出手段を備えた骨計測装置であって、
少なくとも前記骨部組織を含む被検体を被写体とする放射線画像に対してモーフォロジー演算に基づく骨部組織の強調処理を施すモーフォロジー演算手段をさらに備え、
前記骨構造指標値算出手段が、前記強調処理された骨部組織の構造を表す骨部組織画像に対して、前記指標値を取得する処理を施すものであることを特徴とするものである。
【００８０】
ここで、骨部組織、骨部組織の構造の状態、指標値、モーフォロジー演算に基づく骨部組織の強調処理、およびスケルトン処理については、上述した本発明の骨計測方法におけるものと同義である。
【００８１】
【発明の効果】
本発明の骨計測方法および装置によれば、放射線画像を表す画像情報に対して、スケルトン処理等の、モーフォロジー演算に基づく骨部組織の強調処理を施すことにより、被検体に非侵襲で骨部組織の構造状態を精度よく抽出または強調処理することができ、その抽出または強調処理された骨部組織の構造を表す画像に基づいて、スターボリューム法による指標値やノードストラット解析法による指標値等の、当該骨部組織の構造状態を表す指標値を算出することにより、特に骨粗鬆症の診断に有用な情報である骨部組織の構造状態を、定量的な指標値として提供することができる。
【００８２】
【発明の実施の形態】
以下、本発明の骨計測装置の具体的な実施の形態について図面を用いて説明する。
【００８３】
図１は本発明の骨計測装置の一実施形態の構成を示すブロック図である。
【００８４】
図示の骨計測装置100 は、入力された、軟部組織Paと骨部組織Pbとを含む人体等を被写体とした放射線画像Ｐを表す画像信号（以下、放射線画像情報という）Ｓに対して、モーフォロジー演算に基づいたスケルトン処理を施すことにより当該骨部組織Pbに相当する画像（スケルトン画像）を表す画像情報Sbを抽出するモーフォロジー演算手段10と、このモーフォロジー演算手段10により抽出された画像情報Sbに対して、当該骨部組織Paの構造状態を表す指標値Ｖｔを求める骨構造指標値算出手段20とを備えた構成である。
【００８５】
ここで、上記本実施形態の骨計測装置 100に入力される画像情報Ｓは、例えば、放射線を照射するとこの放射線のエネルギーの一部が蓄積され、その後、可視光やレーザー光等の励起光を照射すると蓄積された放射線エネルギーに応じた輝尽発光を示す蓄積性蛍光体（輝尽性蛍光体）を利用して、人体等の少なくとも骨部組織を含む被写体の放射線画像情報を一旦、シート状の蓄積性蛍光体（蓄積性蛍光体シート）に記録し、この蓄積性蛍光体シートをレーザー光等の励起光で走査して輝尽発光光を生じせしめ、この輝尽発光光を光電的に読み取って得られた画像信号である。なお、この画像情報としては、単一の放射線画像に基づいて得られた画像情報だけでなく、同一被写体についての、エネルギー分布状態が互いに異なる２つのエネルギーサブトラクション処理用画像情報に基づいて得られる、骨部組織をある程度抽出した画像情報を適用することもでき、また、骨部組織の構造状態の指標値を求めようとする画像部分が、全体画像のうちの一部である場合は、当該画像部分にＲＯＩ（region of interest；関心領域）設定を行ったうえで、当該ＲＯＩ内の画像部分に対応する画像情報のみを適用してもよい。
【００８６】
モーフォロジー演算手段10によるスケルトン処理は、具体的には図２に示すように、処理対象である画像情報Ｓに対して、ダイレーション処理を８回（ｎ＝０〜７）行い、各ダイレーション処理後の各画像情報Ａ（０）、Ａ（１）、…、Ａ（７）に対してそれぞれエロージョン処理を施して画像情報Ｂ（０）、Ｂ（１）、…、Ｂ（７）を得、これらの各結果Ｂ（ｎ）をそれぞれのダイレーション処理前の画像情報Ａ（ｎ−１）（ただし、Ａ（−１）はＳに等しい）から減じる演算を行って差分画像情報Ｓ（ｎ）＝｛Ａ（ｎ−１）−Ｂ（ｎ）｝を求める。
【００８７】
そして得られた各差分画像情報Ｓ（ｎ）をそれぞれ予め設定された閾値により２値化してスケルトン画像情報Ｓkl（ｎ）＝｛Ｓkl（０），Ｓkl（１），Ｓkl（２），Ｓkl（３），Ｓkl（４），Ｓkl（５），Ｓkl（６），Ｓkl（７）｝を得る。
【００８８】
さらにこれらのスケルトン画像情報Ｓkl（ｎ）のうち、ｉ≦ｎ≦ｊなる（ｉ，ｊ）を用いて、論理和ＳＵＭ（ｉ，ｊ）＝Ｓkl（ｉ）＋Ｓkl（ｉ＋１）＋…＋Ｓkl（ｊ）が表すサムセット画像情報を得る。一例としては、ＳＵＭ（０，１）、ＳＵＭ（２，５）、ＳＵＭ（３，６）、ＳＵＭ（４，７）の４つが適当である。
【００８９】
そしてこれらの４つのサムセット画像情報ＳＵＭ（０，１）、ＳＵＭ（２，５）、ＳＵＭ（３，６）、ＳＵＭ（４，７）と８つのスケルトン画像情報Ｓkl（０），Ｓkl（１），Ｓkl（２），Ｓkl（３），Ｓkl（４），Ｓkl（５），Ｓkl（６），Ｓkl（７）が、スケルトン処理の結果である骨部画像情報Sbとして、モーフォロジー演算手段10から出力される。
【００９０】
また、骨構造指標値算出手段20による、骨部画像情報Sbが表す骨部組織Paの構造状態を表す指標値Ｖｔを求める処理としては、本実施形態においては、前述したスターボリューム法が適用される。そして本実施形態のスターボリューム法では、骨梁の中の点から延ばす放射状の全方向として、例えば９度ごとの40方向を設定し、骨梁の中の点からこれらの全ての方向（40方向）について、骨梁の端に至る範囲（図３参照）の体積の平均を意味する、骨部組織の構造状態を表す１つの指標値Ｖｔを求める。
【００９１】
ここでまず、任意のサンプリング点ｊについての指標値Ｖｔｊを次式（２）にしたがって求める。
【００９２】
Ｖｔｊ＝（π／３）×Σｌ14／Σｌ1   （２）
但し、ｌ1 は点ｊを中心とする任意の方向において、骨梁が連続する長さを示す。また、ｌ1 は点ｊを中心として放射状に全ての方向において求められ、Σはその全ての方向における和を求めることを表す。
【００９３】
そして、Ｖｔｊを各サンプリング点ごとに算出し、それらの平均値を指標値Ｖｔとすればよい。
【００９４】
次に本実施形態の骨計測装置100 の作用について説明する。
【００９５】
まず、前述した骨部画像情報Sbを含む放射線画像情報Ｓがモーフォロジー演算手段10に入力され、このモーフォロジー演算手段10は入力された画像情報Ｓに対してスケルトン処理（図２）を施して骨部画像情報Sbを抽出する。
【００９６】
得られた骨部画像情報Sbは骨構造指標値算出手段20に入力され、骨構造指標値算出手段20は、入力された骨部画像情報Sbに基づき、骨部組織Paの構造状態を表す１つの指標値であるＶｔ値を上記式（２）にしたがって算出する。
【００９７】
このように算出された指標値Ｖｔは骨計測装置100 から出力され、骨粗鬆症の診断等に供される。
【００９８】
このように本実施形態の骨計測装置100 によれば、モーフォロジー演算に基づくスケルトン処理により、放射線画像中における骨部組織を軟部組織から精度よく分離した画像情報を得ることができ、そしてこのように精度よく抽出された骨部組織を表す画像情報に基づいて当該骨部組織の構造を表す指標値を求めることにより、骨構造の診断に有用な情報を定量的なものとして、かつ放射線画像から精度よく得ることができる。
【００９９】
なお、上記実施形態における骨構造指標値算出手段20は、骨部組織の構造を表す指標値として骨梁体積の平均を意味するＶｔ値を算出するものとしたが、本発明の骨計測装置はこの実施形態に限るものではなく、骨髄腔体積の平均を意味するＶｍ値（式（１））や、ノードストラット解析法による、すべてのストラットの長さの和ＴSLに対する各ストラットの長さの割合（％）、骨部組織（ここでは、骨梁および皮質骨の和、または骨梁のみ）の面積に対するストラットの長さ、結合点または終末点の数等を骨部組織の構造を表す指標値として求める骨構造指標値算出手段を適用することができる。
【図面の簡単な説明】
【図１】本発明の骨計測装置の一実施形態を示すブロック図
【図２】図１に示した骨計測装置におけるモーフォロジー演算手段10によるスケルトン処の流れを示すフローチャート
【図３】スターボリューム法による指標値を説明するための図
【図４】骨部組織の状態を判定する方法を示す図
【図５】ノードストラット解析法による指標値を説明するための図
【図６】モーフォロジー処理の基本的な作用を示す図、（Ａ）ダイレーション（dilation）処理、（Ｂ）エロージョン（erosion ）処理、（Ｃ）オープニング（opening ）処理、（Ｄ）クロージング（closing ）処理
【図７】モーフォロジー処理に用いる構造要素Ｂi の一例を示す図
【図８】クロージング処理を具体的に説明する図
【図９】種々の図形とその骨格（skeleton；スケルトン）を示す図
【図１０】スケルトン処理を具体的に説明する図
【符号の説明】
10    モーフォロジー演算手段
20    骨構造指標値算出手段
100   骨計測装置
【図１】
0003656695.tif 000009
【図２】
0003656695.tif 000010
【図３】
0003656695.tif 000011
【図４】
0003656695.tif 000012
【図５】
0003656695.tif 000013
【図６】
0003656695.tif 000014
【図７】
0003656695.tif 000015
【図８】
0003656695.tif 000016
【図９】
0003656695.tif 000017
【図１０】
0003656695.tif 000018