(57)【要約】【課題】  骨密度を正確に評価する。【解決手段】  下顎骨とリファレンスバーとが並んで撮影されたレントゲン撮影画像２１４のうち、リファレンスバー領域２４８が選択されると、この領域２４８の輝度の平均値および偏差が、文字列２２４および２２６によって表示される。一方、評価対象領域２５４が選択されると、この領域２５４の輝度が、リファレンスバー領域２４８の輝度の平均値および偏差、並びに基準平均値、基準偏差に基づいて、つまり一定の基準の下、補正される。そして、補正後の輝度は、ヒストグラム（曲線２６４および２６６）で表示される。さらに、補正後の輝度の平均値および偏差が求められ、これらは文字列２３０および２３２によって表示される。そして、これらヒストグラム、補正後の輝度の平均値および偏差から、骨密度が定量的に評価される。【選択図】        図１６
2005058751.tif 000002
【特許請求の範囲】
【請求項１】
下顎骨のレントゲン撮影画像を基に骨密度を評価する骨密度評価装置であって、
上記レントゲン撮影画像は上記下顎骨と並んで配置された標本体の画像を含み、
上記標本体の画像の濃淡度を検出する検出手段と、
上記検出手段による検出結果が基準値と一致するように上記レントゲン撮影画像の濃淡度を補正する補正手段と、
上記補正手段によって補正された補正後濃淡度に基づいて上記骨密度を評価する評価手段とを具備することを特徴とする、骨密度評価装置。
【請求項２】
上記評価手段は上記レントゲン撮影画像のうち特定領域の上記補正後濃淡度に基づいて評価を行う、請求項１に記載の骨密度評価装置。
【請求項３】
上記特定領域は第１小臼歯の周囲の歯槽骨部分に対応する領域を含む、請求項２に記載の骨密度評価装置。
【請求項４】
上記検出手段は上記標本体の画像の特定部分の濃淡度を検出する、請求項１ないし３のいずれかに記載の骨密度評価装置。
【請求項５】
上記標本体の画像の濃淡度は該画像の場所によって異なり、
上記検出手段は上記標本体の画像の濃淡度の平均および偏差の一方または両方を検出する、請求項１ないし３のいずれかに記載の骨密度評価装置。
【請求項６】
上記基準値を任意に設定する設定手段をさらに備える、請求項１ないし５のいずれかに記載の骨密度評価装置。
【請求項７】
上記基準値は特定の上記レントゲン撮影画像についての上記検出手段による検出結果に基づいて設定された、請求項１ないし５のいずれかに記載の骨密度評価装置。
【請求項８】
上記評価手段は上記補正後濃淡度をヒストグラムで表示する表示手段を含む、請求項１ないし７のいずれかに記載の骨密度評価装置。
【請求項９】
上記評価手段は上記補正後濃淡度に基づいて上記骨密度の程度を判定する判定手段を含む、請求項１ないし８のいずれかに記載の骨密度評価装置。
【請求項１０】
上記評価手段による複数の上記レントゲン撮影画像のそれぞれに係る複数の評価結果を同時に出力する出力手段をさらに備える、請求項１ないし９のいずれかに記載の骨密度評価装置。
【請求項１１】
下顎骨のレントゲン撮影画像を基に骨密度を評価する骨密度評価方法であって、
上記レントゲン撮影画像は上記下顎骨と並んで配置された標本体の画像を含み、
上記標本体の画像の濃淡度を検出する検出過程と、
上記検出過程における検出結果が基準値と一致するように上記レントゲン撮影画像の濃淡度を補正する補正過程と、
上記補正過程において補正された補正後濃淡度に基づいて上記骨密度を評価する評価過程とを具備することを特徴とする、骨密度評価方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、骨密度評価装置および骨密度評価方法に関し、特に例えば下顎骨のレントゲン撮影画像を基に骨密度を評価する、骨密度評価装置および骨密度評価方法に関する。
【背景技術】
【０００２】
近年、下顎骨の骨密度（骨塩量）から身体全体の骨密度を評価する研究が注目されている。特に、非特許文献１には、下顎骨のうち歯槽骨の骨密度、具体的には小臼歯間（第１小臼歯と第２小臼歯との間）で、かつ歯頚部（エナメル質とセメント質との結合部分）から歯根尖側に向けて約６［ｍｍ］離れた位置における骨密度が、身体全体の骨密度と密接に関係することが、紹介されている。そして、歯槽骨のレントゲン撮影画像の濃淡度から当該歯槽骨の骨密度を推測し、その推測結果から身体全体の骨密度を評価する実験例が、記載されている。
【非特許文献１】Grethe Jonasson, GudrunBankvall, Stavros Kiliaridis : Estimation of skeletal bone mineral density bymeans of the trabecular pattern of the alveolar bone, its interdentalthickness, and the bone mass of the mandible. “ORAL SURGERY ORAL MEDICINE ORAL PATHOLOGY” September 2001, Volume 92, Number 3, p.346-p.352
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかし、上述の非特許文献１に開示された従来技術は、飽くまで研究レベルの技術であり、単に歯槽骨のレントゲン撮影画像の濃淡度から身体全体の骨密度を感覚的に推測できることを紹介したものに過ぎない。つまり、この従来技術では、身体全体の骨密度を正確に評価することができない。なぜなら、レントゲン撮影画像の撮影時や現像時、保管時等の諸条件（撮影時間の長さや現像液の温度、レントゲンフィルムの保管環境等）によって当該レントゲン撮影画像の濃淡度に差異が生じ、換言すれば個々のレントゲン撮影画像毎に骨密度を評価する際の基準が異なるからである。
【０００４】
そこで、この発明は、骨密度を正確に評価することができる実用的な骨密度評価装置および骨密度評価方法を提供すること、を目的とする。
【課題を解決するための手段】
【０００５】
第１の発明は、下顎骨のレントゲン撮影画像を基に骨密度を評価する骨密度評価装置であって、当該レントゲン撮影画像は下顎骨と並んで配置された標本体の画像を含む。そして、標本体の画像の濃淡度を検出する検出手段と、この検出手段による検出結果が基準値と一致するようにレントゲン撮影画像の濃淡度を補正する補正手段と、この補正手段によって補正された補正後濃淡度に基づいて骨密度を評価する評価手段と、を具備することを特徴とするものである。
【０００６】
即ち、第１の発明では、１枚のレントゲン撮影画像に下顎骨の画像と標本体の画像とが並んで表示されている。このうち標本体の画像の濃淡度が、検出手段によって検出される。そして、この検出手段による検出結果が基準値と一致するように、レントゲン撮影画像の濃淡度、詳しくは下顎骨の画像の濃淡度が、補正手段によって補正される。さらに、この補正手段による補正後濃淡度に基づいて、下顎骨の骨密度、ひいては身体全体の骨密度が、評価手段によって評価される。つまり、特定の標本体の画像の濃淡度が基準値と一致するように補正された、換言すれば一定の基準の下で補正された、レントゲン撮影画像の濃淡度に基づいて、骨密度が評価される。
【０００７】
なお、評価手段は、レントゲン撮影画像のうち特定領域の補正後濃淡度に基づいて評価を行うようにしてもよい。
【０００８】
この場合、特定領域は、第１小臼歯の周囲の歯槽骨部分に対応する領域を含むのが、望ましい。
【０００９】
また、検出手段は、標本体の画像の特定部分、例えば最も濃淡度の高い部分、或いは低い部分の濃淡度を検出するものであってもよい。つまり、当該特定部分の濃淡度が基準値と一致するように、レントゲン撮影画像の濃淡度が補正されるようにしてもよい。
【００１０】
さらに、例えば標本体の厚さ寸法が一様でない場合のように、当該標本体の画像の濃淡度が場所によって異なることがある。このような場合は、検出手段によって、当該濃淡度の平均または偏差を検出するのが、望ましい。つまり、標本体画像の濃淡度の平均または偏差が基準値と一致するように、レントゲン撮影画像の濃淡度が補正されるようにしてもよい。
【００１１】
そしてさらに、上述の基準値を任意に設定する設定手段を設けてもよい。
【００１２】
また、基準値は、特定のレントゲン撮影画像についての検出手段による検出結果に基づいて設定されたものであってもよい。つまり、当該特定のレントゲン撮影画像の濃淡度を基準として、他のレントゲン撮影画像の濃淡度が補正されるようにしてもよい。この場合、他のレントゲン撮影画像の濃淡度は、特定のレントゲン撮影画像の濃淡度に合わせられる。
【００１３】
そして、評価手段は、補正後濃淡度をヒストグラムで表示する表示手段を含むものとしてもよい。
【００１４】
さらに、評価手段は、補正後濃淡度に基づいて骨密度の程度、例えば当該骨密度が正常であるか否か、或いは改善を要するか否かを判定する判定手段を含むものとしてもよい。
【００１５】
また、複数のレントゲン撮影画像のそれぞれについての評価手段による複数の評価結果を、同時に出力する出力手段をさらに設けてもよい。このようにすれば、例えば骨密度の経時的変化を把握するのに、極めて有用である。
【００１６】
第２の発明は、第１の発明に対応するいわゆる方法発明である。即ち、下顎骨のレントゲン撮影画像を基に骨密度を評価する骨密度評価方法であって、当該レントゲン撮影画像は下顎骨と並んで配置された標本体の画像を含む。そして、標本体の画像の濃淡度を検出する検出過程と、この検出過程における検出結果に基づいてレントゲン撮影画像の濃淡度を補正する補正過程と、この補正過程において補正された補正後濃淡度に基づいて骨密度を評価する評価過程と、を具備することを特徴とするものである。
【発明の効果】
【００１７】
この発明によれば、一定の基準の下で補正された下顎骨のレントゲン撮影画像の濃淡度に基づいて、骨密度が評価される。従って、上述した従来技術とは異なり、レントゲン撮影画像毎の濃淡度の差異に関係なく、骨密度を定量的に、つまり正確に評価することができる。即ち、十分に実用的な骨密度評価装置または骨密度評価方法を実現できる。
【発明を実施するための最良の形態】
【００１８】
この発明に係る歯科医療用の骨密度評価システムの第１実施形態について、図１〜図２１を参照して説明する。
【００１９】
図１に示すように、この第１実施形態の骨密度評価システム１０は、パーソナルコンピュータ（以下、ＰＣと言う。）１２を備えている。このＰＣ１２には、ＣＤ−ＲＯＭ（Compact Disc ROM）１８から、骨密度評価プログラムがインストールされる。ＰＣ１２は、この骨密度評価プログラムが起動されることで、骨密度評価装置として機能する。さらに、ＰＣ１２には、画像入力手段としてのフィルムスキャナ１４と、印刷手段としてのレーザプリンタ１６とが、接続されている。また、ＰＣ１２は、命令入力手段としてのキーボード２０およびマウス２２と、表示手段としてのディスプレイ２４とを、備えている。
【００２０】
このような構成の骨密度評価システム１０では、レントゲンフィルム２６に撮影された下顎骨の撮影画像が、フィルムスキャナ１４によって読み取られ、ディジタル画像データに変換される。変換された画像データは、ＰＣ１２に入力され、当該ＰＣ１２内の図示しないハードディスクに、例えば８ビットまたは２４ビットのビットマップ形式で記録される。この要領で、ハードディスクには、同一患者の同一箇所（下顎骨）についてそれぞれ異なる日に撮影された複数の画像データが、順次記録される。
【００２１】
なお、レントゲンフィルム２６には、下顎骨のうち主に第１小臼歯付近の画像、より具体的には、図２に示すように、第１小臼歯２８と、これと隣り合う犬歯３０および第２小臼歯３２と、これらの歯２８〜３２を支える歯槽骨３４との画像が、撮影されている。また、これらの歯２８〜３２（厳密には各歯２８〜３２の画像）と並んで、例えば当該各歯２８〜３２の上方に、後述するリファレンスバー３６も撮影されている。
【００２２】
ＰＣ１２は、かかる下顎骨の画像データ、具体的には図２に点線３８で示すように第１小臼歯２８の歯根４０の下側半分とその周囲の歯槽骨３４との画像の濃淡度から、当該歯槽骨３４の骨密度（骨塩量）を評価し、ひいては身体全体の骨密度を評価する。例えば、当該濃淡度が高い（画像が明るい）場合は骨密度が高く、当該濃淡度が低い（画像が暗い）場合は骨密度が低い、と評価する。そして、この評価結果をディスプレイ２４に表示する。この評価結果はまた、必要に応じてハードディスクに記録され、或いはレーザプリンタ２４によってカルテ（またはカルテに添付される用紙）４２に印刷される。
【００２３】
このように第１小臼歯２８の歯根４０の下側半分とその周囲の歯槽骨３４とを評価の対象とするのは、上述した従来技術で紹介されているように、当該第１小臼歯２８の周囲（第１小臼歯と第２小臼歯との間）の歯槽骨３４の骨密度が、身体全体の骨密度と密接に関係するからである。また、歯槽骨３４の上側は歯周病による影響を受け易く、例えば歯周病に罹ると歯槽骨３４の上側の骨密度が低下するため、かかる歯周病による影響を排除するために当該歯槽骨３４（歯根４０）の上側半分を評価対象から外してある。なお、歯根４０自体の成分（画像の濃淡度）は、身体全体の骨密度との間で特別な関係はないので、当該歯根４０については評価対象から外してもよい。ただし、当該歯根４０の成分は殆ど不変であるので、これを評価対象に入れるか否かは骨密度の評価結果には特に影響しない。よって、この第１実施形態では、歯根４０の下側半分についても評価対象に入れてある。
【００２４】
ここで、改めて評価対象について厳密に説明すると、当該評価対象とされる上述の点線３８で囲まれた矩形領域の上側端縁は、第１小臼歯２８の歯頚部から歯根尖までの距離Ｄの半分の位置に設定される。そして、当該矩形領域３８の下側端縁は、第１小臼歯２８の歯根尖の位置に設定される。さらに、矩形領域３８の右側端縁は、各歯２８〜３２の歯頚部を結ぶ直線４４上における第１小臼歯２８と犬歯３０との中間点の位置に合わせられる。これと同様に、左側端縁は、当該直線４４上における第１小臼歯２８と第２小臼歯との中間点の位置に合わせられる。
【００２５】
また、上述のリファレンスバー３６は、図３に示すような標本体としてのアルミニウム製ブロック（以下、このアルミニウム製ブロックについても同じ符号３６を付して説明する。）の撮影画像である。即ち、アルミニウム製ブロック３６は、その長手方向において厚さ寸法が段階的に変化する階段状のものであり、当該階段の段差ΔＴおよび間隔（階段の踏板に相当する部分の長さ）Ｐのそれぞれは、一定とされている。また、段数は、７段〜９段程度とされている。さらに、最上段の上面には、後述するＸ線４６（図５参照）を遮断するべく、アルミニウム箔４８が貼着されている。なお、このアルミニウム製ブロック３６の長さ寸法Ｌは、２０［ｍｍ］程度であり、幅寸法Ｗは、１０［ｍｍ］程度である。また、最下段部分での厚さ寸法（底面５０から最下段の上面までの高さ寸法）Ｔａは、１［ｍｍ］前後であり、最上段部分での厚さ寸法（底面５０から最上段の上面までの高さ寸法）Ｔｂは、例えば６［ｍｍ］〜８［ｍｍ］程度である。
【００２６】
このアルミニウム製ブロック３６は、撮影前に、図４に示すように縦長長方形のレントゲンフィルム２６の撮影面（露光面）５２上に貼着される。具体的には、アルミニウム製ブロック３６は、その長手方向をレントゲンフィルム２６の上端縁に沿う方向に延伸させ、かつ底面５０をレントゲンフィルム２６の撮影面５２に対向させた状態で、当該撮影面５２の上端寄りの部分に貼着される。そして、かかるアルミニウム製ブロック３６が貼着されたレントゲンフィルム２６は、撮影時において、図５に示すように、上端縁（アルミニウム製ブロック３８が設けられている側の端縁）を図示しない患者の上顎側（図５の紙面の表面側）に向け、かつ撮影面５２を外方に向けた状態で、撮影対象である各歯２８〜３２および歯槽骨３４の内側に設置される。このとき、レントゲンフィルム２６は、図示しないフィルムホルダによって支持される。そして、このレントゲンフィルム２６の撮影面５２に向けて、測定対象（各歯２８〜３２および歯槽骨３４）の外側に設置されたレントゲン撮影装置（Ｘ線照射器）５４から、Ｘ線４６が照射される。これによって、上述した図２の如く、各歯２８〜３２の撮影画像と、歯槽骨３４の撮影画像と、アルミニウム製ブロック３６の撮影画像（リファレンスバー３６）とが、１枚のレントゲンフィルム２６に並んで写し出される。
【００２７】
なお、アルミニウム製ブロック３６は上述の如く階段状の形状をしているので、その撮影画像であるリファレンスバー３６は、長手方向において、当該アルミニウム製ブロック３６の厚さ寸法（高さ寸法）に応じた濃淡度を示す。即ち、厚みの小さいところほど、Ｘ線に対する透過率が高いので、暗く（黒っぽく）示され、厚みの大きいところほど、Ｘ線に対する透過率が低いので、明るく（白っぽく）示される。そして、アルミニウム箔４８が貼着されたところは、Ｘ線を殆ど通さないので、最も明るく（例えば真っ白に）示される。
【００２８】
さて、今、上述の如く同一患者についてそれぞれ異なる日に撮影された複数（ここでは４枚以上）の画像データが、既にＰＣ１２内のハードディスクに記録されているとする。この状態で、上述した骨密度評価プログラムが起動されると、ＰＣ１２のディスプレイ２４には、図６に示すようなメイン画面１００が表示される。
【００２９】
即ち、このメイン画面１００の一番上には、横長のタイトルバー１０２が表示される。このタイトルバー１０２内には、例えば患者のカルテの番号を表す横書きの文字列１０４が左詰で表示される。そして、タイトルバー１０２の下方に、当該タイトルバー１０２と同様の横長のメニューバー１０６が表示される。このメニューバー１０６内には、メイン画面１００において操作可能なメニューの内容を表す複数の文字列１０８，１０８，…が、横一列に表示される。さらに、メニューバー１０６の下方には、ツールバー１１０が表示され、このツールバー１１０内には、上述のメニューを図案化した複数のツールボタン１１２，１１２，…が、横一列に表示される。そして、ツールバー１１０の下方に、矩形のフレーム領域１１４が表示される。
【００３０】
このフレーム領域１１４の上側には、概ね正方形の４つのピクチャボックス１１６，１１６，…が横一列に表示される。そして、これらピクチャボックス１１６，１１６，…の下方には、“画像１”，“画像２”，“画像３”および“画像４”という当該各ピクチャボックス１１６，１１６，…に表示される画像の見出しを表す４つの文字列１１８，１１８，…が、この順番で左から表示される。さらに、それぞれの文字列１１８の下方には、後述するサブ画面（ダイアログボックス）２００における補正後の輝度Ｙｏ’［ｉ，ｊ］の平均値Ｍｏ’および偏差Ｄｏ’を表す２つの文字列１２０および１２２が、上下２段で表示される。
【００３１】
また、各ピクチャボックス１１６，１１６，…の並びのさらに右横方には、当該各ピクチャボックス１１６，１１６，…のいずれを有効とするのかを選択するための４つのラジオボタン１２４，１２４，…が、縦一列に表示される。そして、それぞれのラジオボタン１２４の右横方には、当該ラジオボタン１２４に対応する“画像Ａ”（Ａは１〜４のいずれかの数値）を示す文字列１２６が表示されている。さらに、ラジオボタン１２４，１２４，…の並びの上方には、当該ラジオボタン１２４，１２４，…の機能を表す“画像選択”という文字列１２８が表示される。一方、ラジオボタン１２４，１２４，…の並びの下方には、後述する基準平均値ＳＭｂおよび基準偏差ＳＤｂを表す２つの文字列１３０および１３２が、上下２段で表示される。
【００３２】
そして、フレーム領域１１４の下部には、ヒストグラム表示エリア１３４が表示される。このヒストグラム表示エリア１３４については後で詳しく説明するが、当該ヒストグラム表示エリア１３４の横軸には、一定の間隔で垂直の目盛線１３６，１３６，…が記されている。そして、それぞれの目盛線１３６の下方に、当該目盛線１３６の指示値を表す文字列１３８が表示されている。
【００３３】
この図６に示すメイン画面１００において、上述したマウス２２の操作によって“画像１”に対応するラジオボタン１２４がＯＮ（クリック）される、つまり当該“画像１”に対応するピクチャボックス１１６が有効化されるとする。そして、この状態で、マウス２２の操作によってハードディスクから任意の画像データが読み込まれるとする。すると、図７に示すように、“画像１”のピクチャボックス１１６内に、当該画像データに従うレントゲン撮影画像１４０が表示される。なお、このレントゲン撮影画像１４０は、グレースケールで表示される。また、元の画像データがカラーデータの場合には、当該画像データはハードディスクからの読み込み時にグレースケールデータに変換される。この画像データを読み込むためのコマンドは、メニューバー１０６内の“ファイル”メニューに格納されている。
【００３４】
これと同様の要領で、他の“画像２”，“画像３”および“画像４”についても、それぞれに対応するラジオボタン１２４がＯＮされた状態でハードディスクから画像データが読み込まれると、図８に示すように、読み込まれた画像データに従うレントゲン撮影画像１４０がそれぞれに対応するピクチャボックス１１６内に表示される。なお、これらのピクチャボックス１１６，１１６，…に表示されるレントゲン撮影画像１４０，１４０，…は、同一患者のものであり、“画像１”，“画像２”，“画像３”および“画像４”の順番で撮影日が古い（“画像１”の撮影日が一番古い）ものとする。
【００３５】
このようにして各ピクチャボックス１１６，１１６，…内にレントゲン撮影画像１４０，１４０，…が表示された後、“画像１”に対応するラジオボタン１２４がＯＮされ、この状態で、当該“画像１”のレントゲン撮影画像１４０について評価対象とする領域を選択する旨の操作がマウス２２によって成されると、ディスプレイ２４に、図９に示すようなサブ画面２００が表示される。なお、この領域選択を命令するためのコマンドは、メニューバー１０６内にある“編集”メニューに格納されている。
【００３６】
図９において、サブ画面２００の一番上には、横長のタイトルバー２０２が表示される。このタイトルバー２０２には、例えばカルテの番号と、現在有効化されている“画像Ａ”（ここでは“画像１”）という情報と、を表す横書きの文字列２０４が、左詰で表示される。そして、タイトルバー２０２の下方に、矩形のフレーム領域２０６が表示される。
【００３７】
このフレーム領域２０６の左上部分には、後述する基準平均値ＳＭｂおよび基準偏差ＳＤｂを表す２つの文字列２０８および２１０が、上下２段で表示される。そして、これらの文字列２０８および２１０の並びの下方には、概ね正方形の領域選択エリア２１２が表示される。この領域選択エリア２１２には、現在有効化されている“画像Ａ”のレントゲン撮影画像１４０を拡大した拡大画像２１４が表示される。さらに、この領域選択エリア２１２の右横方には、当該領域選択エリア２１２内に表示されている拡大画像２１４を回転させるためのスライダ２１６が表示される。例えば、マウス２２によってこのスライダ２１６が上方側へスライドされると、拡大画像２１４は反時計回りに回転する。これとは反対に、スライダ２１６が下方側へスライドされると、拡大画像２１４は時計回りに回転する。このスライダ２１６による拡大画像２１４の回転角度は、当該スライダ２１６の上方に表示された文字列２１８によって表される。
【００３８】
一方、フレーム領域２０６の右上部分には、横長のリファレンスバー表示エリア２２０が表示される。そして、このリファレンスバー表示エリア２２０の下方の左寄りの部分には、概ね正方形の評価対象表示エリア２２２が表示される。これらのエリア２２０および２２２については、後で詳しく説明する。
【００３９】
さらに、評価対象表示エリア２２２の右横方には、後述する平均値Ｍｂおよび偏差Ｄｂを表す２つの文字列２２４および２２６が、上下２段で表示される。そして、これらの文字列２２４および２２６の並びの上方には、当該各文字列２２４および２２６の見出しを表す“リファレンスバー”という文字列２２８が表示される。そしてさらに、文字列２２４および２２６の並びの下方には、後述する補正後平均値Ｍｏ’および補正後偏差Ｄｏ’を表す２つの文字列２３０および２３２が、上下２段で表示される。そして、これらの文字列２３０および２３２の並びと上述の文字列２２４および２２６の並びとの間には、文字列２３０および２３２の見出しを表す“評価対象（補正値）”という文字列２３４が表示される。
【００４０】
また、評価対象表示エリア２２２の下方には、２つのラジオボタン２３６および２３６が、横一列に表示される。そして、それぞれのラジオボタン２３６の右横方には、当該ラジオボタン２３６の見出しを表す文字列２３８が表示される。さらに、ラジオボタン２３６および２３６の並びの上方には、当該ラジオボタン２３６および２３６の機能を表す“選択範囲”という文字列２４０が表示される。
【００４１】
そして、フレーム領域２０６の下部には、ヒストグラム表示エリア２４２が表示される。このヒストグラム表示エリア２４２については後で詳しく説明するが、当該ヒストグラム表示エリア２４２の横軸には、一定の間隔で垂直の目盛線２４４．２４４．…が記されている。そして、それぞれの目盛線２４４の下方には、当該目盛線２４４の指示値を表す文字列２４６が表示される。
【００４２】
この図９において、マウス２２の操作によって“リファレンスバー”と記されたラジオボタン２３６がＯＮされるとする。そして、領域選択エリア２１２に表示されたレントゲン撮影画像２１４上において、図１０に点線２４８で示すように、マウス２２の操作によって上述したリファレンスバー３６に対応する部分２５０が選択されたとする。具体的には、当該点線２４８で示される領域選択用の矩形の描画開始点Ｓがマウス２２のクリック操作によって指定された後、当該矩形２４８の描画終了点Ｅがマウス２２のドラッグ・アンド・ドロップ操作によって指定されたとする。なお、このとき、リファレンスバー２５０（３６）の長手方向（図１０の左右方向）において、矩形２４８の一辺が当該リファレンスバー２５０の一端縁と一致し、かつ当該一辺と対向する矩形２４８の他辺がリファレンスバー２５０の他端縁と一致するようにするのが、肝要である。一方、リファレンスバー２５０の幅方向（図１０の上下方向）においては、矩形２４８の各辺（上辺および下辺）が当該リファレンスバー２５０の周縁よりも内側にあればよい。
【００４３】
このようにしてマウス２２の操作によってリファレンスバー２５０としての領域２４８が選択されると、図１１に示すように、当該リファレンスバー領域２４８の拡大画像２５２が、リファレンスバー表示エリア２２０に表示される。そして、このリファレンスバー領域２４８を構成する各画素の輝度Ｙｂ［ｉ，ｊ］（［ｉ，ｊ］：各画素の座標）が２５６階調に振り分けられる（変換される）と共に、次の数１に基づいて、それぞれの階調ｘ（ｘ＝０〜２５５）毎の度数Ｈｂ［ｘ］が求められる。
【００４４】
【数１】
2005058751.tif 000003
【００４５】
この数１において、ｎｂ［ｘ］は、階調ｘに振り分けられた画素数であり、Ｎｂは、リファレンスバー領域２４８を構成する全ての画素数である。このように階調ｘ毎の画素数ｎｂ［ｘ］を全画素数Ｎｂで除することによって、いわゆる正規化された度数Ｈｂ［ｘ］が求められる。
【００４６】
この数１の演算結果Ｈｂは、図１２に示すテーブル３００に記録される。このテーブル３００は、骨密度評価プログラムが起動されたときに、ハードディスク内に形成される。
【００４７】
そして、この度数Ｈｂ［ｘ］を用いて、リファレンスバー領域２４８における輝度Ｙｂ［ｉ，ｊ］の平均値Ｍｂおよび偏差Ｄｂが求められる。具体的には、平均値Ｍｂは、数２によって求められ、偏差Ｄｂは、数３によって求められる。
【００４８】
【数２】
2005058751.tif 000004
【００４９】
【数３】
2005058751.tif 000005
【００５０】
これら数２および数３によって求められた平均値Ｍｂおよび偏差Ｄｂもまた、上述のテーブル３００に記録される。そして、これら平均値Ｍｂおよび偏差Ｄｂは、図１１に示すように、文字列２２４および２２６によって表示される。
【００５１】
さらに、現在有効化されている画像が“画像１”である場合には、平均値Ｍｂおよび偏差Ｄｂは、それぞれ上述した基準平均値ＳＭｂおよび基準偏差ＳＤｂとして設定される。換言すれば、図１２のテーブル３００において、“画像１”の欄に記録された平均値Ｍｂおよび偏差Ｄｂは、それぞれ基準平均値ＳＭｂおよび基準偏差Ｄｂとして取り扱われる。そして、これらの基準平均値ＳＭｂおよび基準偏差Ｄｂは、図１１に示すように、文字列２０８および２１０によって表示される。
【００５２】
このようにしてリファレンスバー領域２４８の選択が行われた後、次に、“評価対象”と記されたラジオボタン２３６がＯＮされたとする。そして、領域選択エリア２１２に表示されたレントゲン撮影画像２１４上において、図１３に点線２５４で示すように、上述した図２における矩形領域３８に対応する部分が選択されたとする。即ち、第１小臼歯２８に対応する画像２５６のうち歯根４０の下側半分に対応する部分２５８と、その周囲の歯槽骨３４に対応する部分２６０とが、領域選択用の矩形２５４で囲まれたとする。
【００５３】
すると、図１４に示すように、矩形２５４で囲まれた評価対象としての領域の拡大画像２６２が、評価対象表示エリア２２２に表示される。そして、当該評価対象領域２５４を構成する各画素の輝度Ｙｏ［ｉ，ｊ］が２５６階調に振り分けられると共に、次の数４に基づいて、それぞれの階調ｘ毎の度数Ｈｏ［ｘ］が求められる。
【００５４】
【数４】
2005058751.tif 000006
【００５５】
この数４において、ｎｏ［ｘ］は、階調ｘに振り分けられた画素数であり、Ｎｏは、評価対象領域２５４を構成する全ての画素数である。この数４の演算結果Ｈｏは、上述した図１２のテーブル３００に記録される。
【００５６】
そして、この数４の演算結果Ｈｏ［ｘ］を用いて、評価対象領域２５４における輝度Ｙｏ［ｉ，ｊ］の平均値Ｍｏおよび偏差Ｄｏが求められる。具体的には、平均値Ｍｏは、数５によって求められ、偏差Ｄｏは、数６によって求められる。
【００５７】
【数５】
2005058751.tif 000007
【００５８】
【数６】
2005058751.tif 000008
【００５９】
これら数５および数６によって求められた平均値Ｍｏおよび偏差Ｄｏもまた、上述のテーブル３００に記録される。
【００６０】
さらに、上述した基準平均値ＳＭｂおよび基準偏差ＳＤｂと、数２および数３によって求められた平均値Ｍｂおよび偏差Ｄｂとに基づいて、評価対象領域２５４の画像の輝度Ｙｏ［ｉ，ｊ］が補正される。具体的には、次の数７によって補正後輝度Ｙｏ’［ｉ，ｊ］が求められる。
【００６１】
【数７】
2005058751.tif 000009
【００６２】
そして、上述した数４〜数６と同様の要領で、当該補正後輝度Ｙｏ’［ｉ，ｊ］に基づく階調毎度数Ｈｏ’［ｘ］，平均値Ｍｏ’および偏差Ｄｏ’が求められる。これら補正後の階調毎度数Ｈｏ’［ｘ］，平均値Ｍｏ’および偏差Ｄｏ’もまた、上述のテーブル３００に記録される。なお、有効化されている“画像Ａ”が“画像１”である場合は、補正後輝度Ｙｏ’［ｉ，ｊ］は補正前の輝度Ｙｏ［ｉ，ｊ］と等価になる。従って、補正後の階調毎度数Ｈｏ’［ｘ］，平均値Ｍｏ’および偏差Ｄｏ’もまた、補正前の階調毎度数Ｈｏ［ｘ］，平均値Ｍｏおよび偏差Ｄｏと等価になる。そして、補正後の平均値Ｍｏ’および偏差Ｄｏ’は、図１４に示すように、文字列２３０および２３２によって表される。
【００６３】
さらに、ヒストグラム表示エリア２４２には、上述の階調毎度数Ｈｏ［ｘ］および補正後の階調毎度数Ｈｏ’［ｘ］のそれぞれに従う２つのヒストグラムが表示される。即ち、横軸に階調ｘを表し、縦軸に当該階調ｘ毎の度数Ｈｏ［ｘ］またはＨｏ’［ｘ］を表した２本の曲線２６４および２６６が表示される。これらの曲線２６４および２６６は、互いに異なる色で表示され、例えば曲線２６４は青色で表示され、曲線２６６は赤色で表示される。なお、有効化されている“画像Ａ”が“画像１”の場合には、各曲線２６４および２６６は互いに重複して表示され、例えば曲線２６４の上に曲線２６６が重ねて表示される。従って、この場合、外観上は１本の曲線２６６のみが表示されているように見える。また、上述した目盛線２２４，２２４，…およびこれらの指示値を表す文字列２４６，２４６，…は、６４階調毎に表示されている。
【００６４】
このように図１４に示すようなサブ画像２００が表示された後、タイトルバー２０２の右端にある“×”印のボタン（いわゆる閉じるボタン）２６８がクリックされると、ディスプレイ２４の表示は、図１５に示すようなメイン画面１００に戻る。
【００６５】
即ち、補正後の平均値Ｍｏ’および偏差Ｄｏ’が、“画像１”のピクチャボックス１１６の下方の文字列１２０および１２２によって表される。そして、上述した基準平均値ＳＭｂおよび基準偏差ＳＤｂが、メイン画面１００の右端部分に表示された文字列１３０および１３２によって表される。さらに、ヒストグラム表示エリア１３４には、横軸に階調ｘを表し、縦軸に当該階調ｘ毎の補正後の度数Ｈｏ’［ｘ］を表したヒストグラム、つまり図１４における曲線２６６と同様の曲線１４２が表示される。なお、このヒストグラム表示エリア１３４においても、目盛線１３６，１３６，…およびこれらの目盛値を表す文字列１３８，１３８，…は、６４階調毎に表示される。
【００６６】
これと同じ要領で、“画像２”についてもサブ画像２００が表示され、当該サブ画像２００上でリファレンスバー領域２４８と評価対象領域２５４とが選択されると、当該サブ画像２００は図１６に示すようになる。
【００６７】
即ち、サブ画面２００の右端部分に表示された文字列２２４および２２６によって、リファレンスバー領域２４８の輝度Ｙｂ［ｉ，ｊ］の平均値Ｍｂおよび偏差Ｄｂが表される。換言すれば、上述した数２および数３の演算結果ＭｂおよびＤｂが表示される。そして、文字列２３０および２３２によって、評価対象領域２５４の補正後の輝度Ｙｏ’［ｉ，ｊ］の平均値Ｍｏ’および偏差Ｄｏ’が表される。つまり、数５および数６の算出結果Ｍｏ’およびＤｏ’が表示される。さらに、ヒストグラム表示エリア２４２に、補正前の階調毎度数Ｈｏ［ｘ］および補正後の階調毎度数Ｈｏ’［ｘ］のそれぞれに従う２つの曲線２６４および２６６が、互いに異なる色で表示される（図１６においてはこの色の違いを線種の違いで表現している）。
【００６８】
なお、ここで言う補正、つまり上述の数７に基づく補正とは、次のような意味を持つ。即ち、数７に基づいて補正されるということは、“画像Ａ”と“画像１”との間で互いのリファレンスバー領域２４８（図１０参照）の輝度Ｙｂ［ｉ，ｊ］が等価となるように、当該“画像Ａ”の評価対象領域２５４の輝度Ｙｏ［ｉ，ｊ］が補正される、ということである。換言すれば、“画像Ａ”の輝度Ｙｏ［ｉ，ｊ］の基準が、“画像１”の輝度Ｙｏ［ｉ，ｊ］の基準に合わせられる、ということである。
【００６９】
また、図１６に示すように、領域選択エリア２１２内に表示されているレントゲン撮影画像２１４が傾いている場合は、スライダ２１６によって当該レントゲン撮影画像２１４の傾きを矯正する。具体的には、例えばレントゲン撮影画像２１４上のリファレンスバー２５０が水平方向に沿って延伸していないときは、当該リファレンスバー２５０が水平方向に沿って延伸するように当該レントゲン撮影画像２１４の傾きを矯正する。そして、この矯正後、リファレンスバー領域２４８を選択する。一方、各歯２８〜３２に対応する画像２５６，２７０および２７２の並びが水平方向に沿っていない場合には、この並びが水平方向に沿うようにレントゲン撮影画像２１４の傾きを矯正する。そして、この矯正後、評価対象領域２５４を選択する。
【００７０】
このようにして“画像２”についてリファレンスバー領域２４８と評価対象領域２５６とが選択された後、これと同様の要領で、“画像３”および“画像４”のそれぞれについて当該各領域２４８および２５６が選択され、ディスプレイ２４の表示がメイン画面１００に戻されると、当該メイン画面１００は図１７に示すようになる。
【００７１】
即ち、各“画像Ａ”の補正後の平均値Ｍｏ’および偏差Ｄｏ’が、それぞれに対応するピクチャボックス１１６，１１６，…の下方の文字列１２０および１２２によって表示される。そして、上述したように、メイン画面１００の右端部分に表示された文字列１３０および１３２によって、基準平均値ＳＭｂおよび基準偏差ＳＤｂが表される。さらに、ヒストグラム表示エリア１３４に、各“画像Ａ”の補正後の階調毎度数Ｈｏ’［ｘ］に従う４つのヒストグラム（曲線１４２，１４２，…）が、互いに異なる色で表示される。
【００７２】
このように、この第１実施形態によれば、評価対象領域２５４の位置が一様に規定されると共に、一定の基準の下で評価対象領域２５４の輝度Ｙｏ［ｉ，ｊ］が補正される。そして、この補正後輝度Ｙｏ’［ｉ，ｊ］に基づいて、それぞれ撮影日の異なる各“画像Ａ”間で濃淡度が比較される。従って、この比較結果から時間の経過に伴う骨密度の変化を容易に把握することができる。即ち、骨密度を感覚的にしか評価できなかった従来技術とは異なり、当該骨密度を定量的に、つまり正確に評価することができる。
【００７３】
なお、図１７に示すメイン画面１００の内容、言わば骨密度の評価結果は、上述したようにハードディスクに記録（保存）することができる。また、当該評価結果は、レーザプリンタ１６によってカルテ４２に印刷することもできる。このように評価結果を記録したり印刷したりするためのコマンドは、メイン画面１００のメニューバー１０６内にある“ファイル”メニューに格納されている。そして、メイン画面１００のタイトルバー１０２の右端にある“×”印のボタン１４４がクリックされると、骨密度評価プログラムが終了される。
【００７４】
このように骨密度評価プログラムが起動されているとき、ＰＣ１２（厳密にはＰＣ１２内のＣＰＵ）は、図１８〜図２１のフローチャートに示される各処理を実行する。
【００７５】
即ち、図１８に示すように、骨密度評価プログラムが起動されると、ＰＣ１２は、ステップＳ１に進み、上述した図１２のテーブル３００の内容を参照する。そして、ステップＳ３に進み、当該ステップＳ１における参照結果に基づいて、ディスプレイ２４にメイン画面１００を表示する。
【００７６】
メイン画面１００の表示後、ＰＣ１２は、ステップＳ５に進み何らかの操作が成されるのを待つ。つまり、コマンド待ち状態となる。そして、何らかの操作が成されると、ＰＣ１２は、ステップＳ７に進み、その操作内容を判定する。
【００７７】
ここで、骨密度評価プログラムを終了する旨の操作が成される、つまり上述の“×”印のボタン１４４がマウス２２によってクリックされると、ＰＣ１２は、当該骨密度評価プログラムを終了する。
【００７８】
一方、いずれかのラジオボタン１２４がクリック（ＯＮ）されると、ＰＣ１２は、ステップＳ９に進み、クリックされたラジオボタン１２４のみにチェックマーク（黒丸）が表示されるように、各ラジオボタン１２４，１２４，…の表示を制御する。そして、ステップＳ１１に進み、当該クリックされたラジオボタン１２４に対応する“画像Ａ（Ａは１〜４のいずれかの数字）”を有効化した後、ステップＳ５に戻る。
【００７９】
また、ステップＳ７においてハードディスクから任意の画像を読み込む旨の操作が成されたと判定すると、ＰＣ１２は、ステップＳ１３に進む。そして、このステップＳ１３において画像を読み出した後、ステップＳ１５において当該読み出した画像を“画像Ａ”に対応するピクチャボックス１１６に表示する。さらに、ステップＳ１７において、“画像Ａ”の補正後の平均値Ｍｏ’および偏差Ｄｏ’を表示するための文字列１２０および１２２をクリアする（ゼロにする）。そして、ステップＳ５に戻る。
【００８０】
さらに、ステップＳ７において領域選択の旨の操作が成されたと判断すると、ＰＣ１２は、ステップＳ１９に進み、“画像Ａ”が既に読み込まれているか否かを判断する。ここで、“画像Ａ”が未だ読み込まれていない場合は、ステップＳ２１に進み、図には示さないが、“画像Ａ”が未だ読み込まれていないことを表すメッセージをディスプレイ２４に表示する。そして、このメッセージを一定時間表示した後、ステップＳ１０５に戻る。
【００８１】
一方、“画像Ａ”が読み込まれている場合には、ステップＳ１９からステップＳ２３に進み、領域選択処理を行う。即ち、上述したサブ画面２００の機能を実現するための処理を実行する。この領域選択処理の詳細については、後で説明する。そして、このステップＳ２３の領域選択処理を終えると、ＰＣ１２は、ステップＳ１に戻る。
【００８２】
そして、ステップＳ７において評価結果を保存する旨の操作が成されたと判断すると、ＰＣ１２は、ステップＳ２５に進む。そして、このステップＳ２５において評価結果、つまりメイン画面１００の内容をハードディスクに記録した後、ステップＳ１０５に戻る。
【００８３】
また、ステップＳ７において評価結果を印刷する旨の操作が成されたと判断すると、ＰＣ１２は、ステップＳ２７に進み、当該評価結果をレーザプリンタ１６に印刷させる。そして、印刷終了後、ステップＳ５に戻る。
【００８４】
さらに、ステップＳ７において上述以外の操作が成されたと判断すると、ＰＣ１２は、ステップＳ２９に進み、その操作に応じた処理を実行した後、ステップＳ５に戻る。
【００８５】
さて、上述のステップＳ２３の領域選択処理において、ＰＣ１２は、具体的に次のような処理を実行する。即ち、図１９に示すように、領域選択処理に入ると、ＰＣ１２は、まず、ステップＳ１０１においてディスプレイ２４にサブ画面２００を表示する。そして、ステップＳ１０３において何らかの操作が成されるのを待ち、何らかの操作が成されると、ステップＳ１０５に進んで当該操作の内容を判定する。
【００８６】
ここで、サブ画面２００の表示を終了する旨の操作が成されると、つまり“×”印のボタン２６８がクリックされると、ＰＣ１２は、ステップＳ１０７に進み、当該サブ画面２００の表示を終了する。そして、改めてメイン画面１００を表示するべく、図１８のステップＳ１に戻る。
【００８７】
一方、ステップＳ１０５においてスライダ２１６が操作されたと判断すると、ＰＣ１２は、ステップＳ１０９に進む。そして、このステップＳ１０９において、マウス２２の操作に応じてサブ画面２００上のスライダ２１６の表示状態を制御し、さらにステップＳ１１１において当該スライダ２１６の位置に応じて文字列２１８（レントゲン撮影画像２１４の回転角度）の表示を制御する。そしてさらに、ステップＳ１１３において、レントゲン撮影画像２１４の傾きを制御した後、ステップＳ１１５に進み、マウス２２による操作が継続されているか否かを判断する。ここで、操作が継続されている場合は、ステップＳ１０９に戻り、継続されていない場合には、ステップＳ１０３に戻る。
【００８８】
そして、ステップＳ１０５においていずれかのラジオボタン２３６がクリックされたと判断すると、ＰＣ１２は、ステップＳ１１７に進み、クリックされたラジオボタン２３６にのみチェックマークが表示されるように、各ラジオボタン２３６および２３６の表示を制御する。そして、ステップＳ１１９において、当該クリックされたラジオボタン２３６に対応する領域がリファレンスバー領域２４８および評価対象領域２５４のいずれであるのかを認識した後、ステップＳ１０３に戻る。
【００８９】
さらに、ステップＳ１０５において上述のリファレンスバー３６に対応する領域２４８が選択された、つまり“リファレンスバー”と記されたラジオボタン２３６がＯＮされている状態で領域選択エリア２１２内のレントゲン撮影画像２１４上においてマウス２２の操作が成されたと判断すると、ＰＣ１２は、図２０のステップＳ１２１に進む。そして、このステップＳ１２１において、マウス２２の操作に応じてリファレンスバー領域２４８を設定した後、ステップＳ１２３に進み、当該リファレンスバー領域２４８の画像をリファレンスバー表示エリア２２０に拡大表示する。
【００９０】
そして、ＰＣ１２は、ステップＳ１２５において、リファレンスバー領域２４８の画像の輝度Ｙｂ[ｉ，ｊ]を取得する。このとき、当該輝度Ｙｂ[ｉ，ｊ]を２５６階調に振り分ける。そして、ステップＳ１２７において、上述した数１に基づいて階調毎度数Ｈｂ[ｘ]を演算し、ステップＳ１２９において、当該演算結果Ｈｂ[ｘ]を上述した図１２のテーブル３００に記録する。
【００９１】
さらに、ＰＣ１２は、ステップＳ１３１において、上述の数２に基づいて平均値Ｍｂを求めると共に、数３に基づいて偏差Ｄｂを求め、ステップＳ１３３において、これらの演算結果ＭｂおよびＤｂをテーブル３００に記録する。そして、ステップＳ１３５において、これら平均値Ｍｂおよび偏差Ｄｂを文字列２２４およびお２２６で表した後、ステップＳ１３７において、現在有効化されているのは“画像１”であるか否かを判断する。
【００９２】
ここで、“画像１”が有効化されていると判断すると、ＰＣ１２は、ステップＳ１３９に進み、上述の平均値Ｍｂおよび偏差Ｄｂを、それぞれ基準平均値ＳＭｂおよび基準偏差Ｄｂとして、文字列２０８および２１０によって表示する。そして、この表示後、図１９のステップＳ１０３に戻る。一方、現在有効化されているのが“画像１”ではない場合は、ステップＳ１３９をスキップして、直接ステップＳ１０３に戻る。
【００９３】
また、図１９のステップＳ１０５において評価対象領域２５４が選択された、つまり“評価対象”と記されたラジオボタン２３６がＯＮされている状態で領域選択エリア２１２内のレントゲン撮影画像２１４上においてマウス２２の操作が成されたと判断すると、ＰＣ１２は、図２１のステップＳ１４１に進む。そして、このステップＳ１４１において、マウス２２の操作に応じて評価対象領域２５４を設定した後、ステップＳ１４３に進み、当該評価対象領域２５４内の画像を評価対象表示エリア２２２に拡大表示する。
【００９４】
そして、ステップＳ１４５において、評価対象領域２５４内の画像の輝度Ｙｏ[ｉ，ｊ]を取得する。このとき、当該輝度Ｙｏ[ｉ，ｊ]を２５６階調に振り分ける。そして、ステップＳ１４７において、上述した数４に基づいて階調毎度数Ｈｏ[ｘ]を演算し、ステップＳ１４９において、当該演算結果Ｈｏ[ｘ]を図１２のテーブル３００に記録する。さらに、ＰＣ１２は、ステップＳ１５１において、階調毎度数Ｈｏ[ｘ]に従うヒストグラム、つまり曲線２６４を、ヒストグラム表示エリア２４２に表示する。
【００９５】
このステップＳ１５１の実行後、ＰＣ１２は、ステップＳ１５３に進む。そして、このステップＳ１５３において、上述した数５に基づいて平均値Ｍｏを求めると共に、数６に基づいて偏差Ｄｏを求め、ステップＳ１５５において、これらの演算結果ＭｏおよびＤｏをテーブル３００に記録する。そして、ステップＳ１５７において、基準平均値ＳＭｂおよび基準偏差ＳＤｂが設定されているか否かを判断する。
【００９６】
ここで、基準平均値ＳＭｂおよび基準偏差ＳＤｂが設定されていない場合、ＰＣ１２は、図１９のステップＳ１０３に戻る。一方、基準平均値ＳＭｂおよび基準偏差ＳＤｂが設定されている場合は、ステップＳ１５９に進み、上述した数７に基づいて、補正後輝度Ｙｏ’［ｉ，ｊ］を求める。このときも、当該補正後輝度’Ｙｏ[ｉ，ｊ]を２５６階調に振り分ける。
【００９７】
そして、ステップＳ１６１に進み、補正後輝度Ｙｏ’［ｉ，ｊ］に基づく階調毎度数Ｈｏ’[ｘ]を求めた後、ステップＳ１６３において当該演算結果Ｈｏ’[ｘ]を上述のテーブル３００に記録する。そして、ステップＳ１６５において、補正後の階調毎度数Ｈｏ’[ｘ]に従うヒストグラム、つまり曲線２６６を、ヒストグラム表示エリア２４２に表示する。
【００９８】
さらに、ＰＣ１２は、ステップＳ１６７において、補正後の階調毎度数Ｈｏ’[ｘ]に基づいて平均値Ｍｏ’を求めると共に、これら補正後の階調毎度数Ｈｏ’[ｘ]および平均値Ｍｏ’に基づいて偏差Ｄｏ’求める。そして、ステップＳ１６９において、これらの演算結果Ｍｏ’およびＤｏ’をテーブル３００に記録する。さらに、ステップＳ１７１において、これら補正後の平均値Ｍｏ’および偏差Ｄｏ’を文字列２３０および２３２によって表示した後、図１９のステップＳ１０３に戻る。
【００９９】
以上のように、第１実施形態の骨密度評価システム１０は、図１に示すような比較的に簡単かつ安価な構成であるので、開業医のような比較的に小規模な医療施設においても容易に導入することができる。しかも、上述した従来技術とは異なり、定量的に骨密度を評価できるので、歯科治療過程において前置的なスクリーニング（骨粗鬆症であるか否かの判断）が可能となる。このことは、後々の治療計画を立てるのに極めて有用であり、例えば下顎骨の強健さが要求されるインプラント治療の事前検査に特に貢献する。
【０１００】
なお、この第１実施形態においては、ＰＣ１２に骨密度評価プログラムを実行させることで当該ＰＣ１２を骨密度評価装置として機能させるようにしたが、これに限らない。例えば、この第１実施形態と同様の機能を有する専用の骨密度評価装置を構成してもよい。
【０１０１】
また、画像入力手段として、フィルムスキャナ１４を用いたが、これに限らない。例えばＣＣＤ（Charge Coupled Device）カメラ等の撮影手段を用いてもよいし、レントゲン撮影時に直接画像を取り込むいわゆるディジタルカセッテを用いてもよい。そして、印刷手段としてのレーザプリンタ１６に代えて、インクジェットプリンタやドットインパクトプリンタ等の他の印刷方式を採用するプリンタを用いてもよい。また、磁気カードなどの記録媒体に、評価結果を記録してもよい。さらに、ＰＣ１２に入力される画像や当該ＰＣ１２から出力される結果を、インターネット等の通信手段によって遠隔地との間で送受信できるようにしてもよい。
【０１０２】
そして、標本体として階段状のアルミニウム製ブロック３６を用いたが、これに限らない。例えば、厚みが一定のアルミニウム製ブロックを用いてもよいし、アルミニウム以外の材料で構成されたブロックを用いてもよい。また、アルミニウム製ブロック３６に代えて、骨密度が既値である擬似骨を用いてもよい。かかる擬似骨を用いれば、当該擬似骨の画像の濃淡度を基に、下顎骨の骨密度を推測することができ、ひいては身体全体の骨蜜度を推測することができる。
【０１０３】
また、上述したメイン画面１００において４つのレントゲン撮影画像１４０，１４０，…を同時に表示すると共に、これら４つのレントゲン撮影画像１４０，１４０，…に係る評価結果（補正後平均値Ｍｏ’，補正後偏差Ｄｏ’およびヒストグラム（曲線１４２，１４２，…））を同時に表示するようにしたが、これに限らない。即ち、４以外の数のレントゲン撮影画像１４０およびその評価結果を表示するようにしてもよい。
【０１０４】
そして、上述した数７に基づいて評価対象領域２５４の輝度Ｙｏ［ｉ，ｊ］を補正したが、当該数７以外の演算式（アルゴリズム）に基づいて補正を行ってもよい。
【０１０５】
さらに、評価対象領域２５４の形状は矩形に限らない。例えば、歯根４０（図２参照）に対応する部分２５８を当該評価対象領域２５４から排除してもよい。また、当該評価対象領域２５４の形状を任意に設定できるようにしてもよい。このようにすれば、患者毎に微妙に異なる下顎骨の形状（第１小臼歯２８と、これと隣り合う犬歯３０および第２小臼歯３２との並び）の違いに、柔軟かつ的確に対応することができる。
【０１０６】
また、評価対象領域２５４は、この第１実施形態で説明した部分（第１小臼歯２８の歯根４０の下側半分とその周囲の歯槽骨３４との部分）以外の部分としてもよい。ただし、この場合も、各“画像Ａ”間で当該評価対象領域２５４を統一させなければならないことは、言うまでもない。
【０１０７】
次に、この発明の第２実施形態について、図２２〜図４０を参照して説明する。
【０１０８】
上述したように、評価対象領域２５４（図１３参照）の輝度Ｙｏ［ｉ，ｊ］は、骨密度に相関する。即ち、補正後の輝度Ｙｏ’［ｉ，ｊ］の平均値Ｍｏ’が大きいほど、骨密度は高く、当該平均値Ｍｏ’が小さいほど、骨密度は低い。このたび、５０歳〜６９歳（平均年齢約６０歳）の計３５人の人間（サンプル）ついて、実際に、公知のＤＸＡ（Dual energy X-ray Absorptiometry）法による骨密度測定を実施すると共に、その測定値（ＢＭＤ（Bone
Mineral Density）値）と当該補正後平均値Ｍｏ’とを比較してみた。併せて、Ｔ値（２０歳〜４４歳までのＢＭＤ値の平均値を１００としたときの実際のＢＭＤ値の比率）についても、補正後平均値Ｍｏ’と比較してみた。その結果を、図２２に示す。
【０１０９】
この図２２に示すように、補正後平均値Ｍｏ’は、ＢＭＤ値が大きいほど大きく、ＢＭＤ値が小さいほど小さい、という傾向にある。その証拠に、ＢＭＤ値と補正後平均値Ｍｏ’との間の相関係数は“０．６４３９”であり、つまり、これら両者間に正の相関がある、と言える。一方、Ｔ値との関係についても同様である。即ち、Ｔ値と補正後平均値Ｍｏ’との間の相関関数は“０．６４８１”であり、これら両者間にも正の相関がある。このことからも、上述の如く補正後平均値Ｍｏ’が骨密度に相関することは、明らかである。
【０１１０】
ところで、一般（国際的）には、骨密度はＴ値に基づいて評価されることが多い。具体的には、Ｔ値が７０％以下の場合は、骨粗鬆症であると診断される。そして、Ｔ値が７１％〜８０の場合は、骨減少症であると診断される。さらに、Ｔ値が８１％以上の場合には、正常であると診断され、特に当該Ｔ値が１００％以上の場合には、極めて正常（安心）であると診断される。ここで、これらのＴ値の区分毎に、図２２における補正後平均値Ｍｏ’の平均を求めると、図２３に示すようになる。そして、この図２３の関係をグラフで表すと、図２４のようになる。
【０１１１】
これら図２３および図２４においても、当然に、補正後平均値Ｍｏ’はＴ値と相関し、ひいては骨密度と相関する。換言すれば、補正後平均値Ｍｏ’に基づいて骨密度の程度（高低）を判定することができる、と推測される。例えば、補正後平均値Ｍｏ’の平均値が７１．４１％以下である場合は、骨粗鬆症である、厳密には骨粗鬆症の疑いがある、と判定することができる。そして、補正後平均値Ｍｏ’の平均値が８４．８９％以下である場合には、骨減少症である、厳密には骨減少症の疑いがある、と判定することができる。そして、これ以外の場合は、骨密度が正常である、と判定することができる。
【０１１２】
そこで、この第２実施形態では、補正後平均値Ｍｏ’に基づいて骨密度の程度を判定する機能を、上述した骨密度評価プログラムに組み込む。また、第１実施形態においては、“画像１”の輝度Ｙ［ｉ，ｊ］が言わば不変とされたが、この第２実施形態においては、当該“画像１”の輝度Ｙ［ｉ，ｊ］を或る一定の基準に基づいて補正する。これ以外については、第１実施形態と同様であるので、当該第１実施形態と同様の部分についての詳細な説明は省略する。
【０１１３】
即ち、第２実施形態における骨密度評価プログラムが起動されると、ＰＣ１２のディスプレイ２４には、図２５に示すようなメイン画面４００が表示される。
【０１１４】
このメイン画面４００は、上述した第１実施形態におけるメイン画面１００（図６等参照）に対し、エディットボックス４０２および４０４と、文字列４０６，４０８および４１０，４１０，…とを付加したものである。このうち、エディットボックス４０２および４０４は、フレーム領域１１４の右側略中央に表示され、詳しくは基準偏差ＳＤｂを表す文字列１３２の下方に上下２段で表示される。そして、これらのエディットボックス４０２および４０４の左横方に、“要注意”および“注意”という当該各エディットボックス４０２および４０４の見出しを表す２つの文字列４０６および４０８が表示される。なお、各エディットボックス４０２および４０４の一方、例えば“要注意”という見出し（文字列）４０６が付されたエディットボックス４０２には、上述の如く補正後平均値Ｍｏ’の平均値に基づいて骨粗鬆症であるか否かを判定する際の基準となる値、例えば７１．４１％の近似値である“７０”という値が、デフォルト値として入力される。そして、“注意”という見出し４０８が付されたエディットボックス４０４には、当該補正後平均値Ｍｏ’の平均値に基づいて骨減少症であるか否かを判定する際の基準となる値、例えば８４．８９％の近似値である“８５”という値が、デフォルト値として入力される。
【０１１５】
そして、それぞれのピクチャボックス１１６の下方に表示された３つの文字列１１８，１２０および１２２のさらに下方に、当該ピクチャボックス１１６内に表示されるレントゲン撮影画像１４０に基づく判定結果を表す文字列４１０が、表示される。この文字列４１０（判定結果）の表示内容については、後で詳しく説明するが、骨密度評価プログラムの起動直後は、当該文字列４１０として“チェック未”という文字列が表示される。
【０１１６】
さて、今、メイン画面４００において、第１実施形態と同様の要領で“画像１”〜“画像４”に対応する画像データがハードディスクから順次読み込まれ、これによって、図２６に示すように、各ピクチャボックス１１６，１１６，…内に当該“画像１”〜“画像４”に対応するレントゲン撮影画像１４０，１４０，…が表示されたとする。そして、“画像１”に対応するラジオボタン１２４がＯＮされた状態、つまり“画像１”が有効化された状態で、これから輝度Ｙ［ｉ，ｊ］の補正を開始する旨の操作が成されたとする。すると、ディスプレイ２４に、図２７に示すようなサブ画面５００が表示される。なお、かかる補正開始を指示するためのコマンドは、メニューバー１０６内の“編集”メニューに格納されている。
【０１１７】
図２７に示すように、サブ画面５００の一番上には、第１実施形態と同様の横長のタイトルバー５０２が表示される。そして、このタイトルバー５０２には、カルテの番号と、現在有効化されている“画像Ａ”（ここでは“画像１”）と、を表す横書きの文字列５０４が、例えば左詰で表示される。さらに、タイトルバー５０２の下方に、矩形のフレーム領域５０６が表示される。
【０１１８】
このフレーム領域２０６の左側には、第１実施形態のサブ画面２００における領域選択エリア２１２よりも大きめの領域選択エリア５０８が、表示される。そして、この領域選択エリア５０８には、現在有効化されている“画像Ａ”のレントゲン撮影画像１４０を拡大した拡大画像５０９が表示される。
【０１１９】
さらに、領域選択エリア５０８の下方であって、フレーム領域５０６の左下隅の部分には、“上”，“下”，“左”および“右”という個別の文字列５１０，５１０，…が付された４つのラジオボタン５１２，５１２，…が、縦２行横２列に表示されている。そして、これらのラジオボタン５１２，５１２，…（文字列５１０，５１０，…）の上方には、当該ラジオボタン５１２，５１２，…の機能を表す“領域抽出”という文字列５１４が表示される。さらに、ラジオボタン５１２，５１２，…（文字列５１０，５１０，…）の一群の右横方には、“参照チェック”という文字列５１４が付されたチェックボックス５１６が、表示される。
【０１２０】
一方、領域選択エリア５０８の右横方であって、フレーム領域２０６の上部には、４つのラジオボタン５１８，５１８，…が、縦一列に表示される。そして、各ラジオボタン５１８，５１８，…の右横方には、当該各ラジオボタン５１８，５１８，…の見出しを表す個別の文字列、例えば“リファレンスバー”，“最高値（白）”，“最低値（黒）”および“評価対象”という文字列５２０，５２０，…が表示される。さらに、一群のラジオボタン５１８，５１８，…（文字列５２０，５２０，…）の上方には、当該ラジオボタン５１８，５１８，…の機能を表す“選択範囲”という文字列５２２が表示される。
【０１２１】
そして、一群のラジオボタン５１８，５１８，…の右横方であって、フレーム領域５０６の右上隅の部分には、注記５２４を表示するための注記エリア５２６が設けられている。なお、サブ画面５００が表示された直後、換言すれば当該サブ画面２００によって後述する補正処理が実行される前は、当該注記５２４として“未処理”という文字列が表示される。
【０１２２】
さらに、注記エリア５２６の下方には、２つのエディットボックス５２８および５３０が、横並びに表示される。そして、これらのエディットボックス５２８および５３０の上方に、“基準最高値”および“基準最低値”という当該各エディットボックス５２８および５３０の見出しを表す２つの文字列５３２および５３４が表示される。なお、各エディットボックス５２８および５３０の一方、例えば“基準最高値”という見出し（文字列）５３２が付されたエディットボックス５２８には、後述する基準最高値Ｙｓｍａｘが入力される。ただし、サブ画面５００の表示直後は、当該基準最高値Ｙｓｍａｘとして“１８０”というデフォルト値が、入力される。そして、他方の“基準最低値”という見出し５３４が付されたエディットボックス５３０には、後述する基準最低値Ｙｓｍｉｎが入力される。サブ画面５００の表示直後は、このエディットボックス５３０に“３０”というデフォルト値が入力される。
【０１２３】
そしてさらに、各エディットボックス５２８および５３０の下方には、横長のリファレンスバー表示エリア５３６が表示される。そして、このリファレンスバー表示エリア５３６の左側上方には、“リファレンスバー”という当該リファレンスバー表示エリア５３６の見出しを表す文字列５３８が、表示される。
【０１２４】
また、リファレンスバー表示エリア５３６の左側下方には、概ね正方形の評価対象表示エリア５４０が表示される。そして、この評価対象表示エリア５４０の左側上方には、“評価対象”という当該評価対象表示エリア５４０の見出しを表す文字列５４２が、表示される。
【０１２５】
さらに、評価対象表示エリア５４０の右横方には、概ね正方形のヒストグラム表示エリア５４４が表示される。このヒストグラム表示エリア５４４の横軸には、一定の間隔で垂直の目盛線５４６．５４６．…が記されている。また、このヒストグラム表示エリア５４４の上方には、３つのラジオボタン５４８，５４８，…が、横一列に表示される共に、各ラジオボタン５４８，５４８，…の右横方には、“原画”，“バー”および“対象”という当該各ラジオボタン５４８，５４８，…の見出しを表す個別の文字列５５０，５５０，…が表示される。さらに、各ラジオボタン５４８，５４８，…（文字列５５０，５５０，…）の上方には、当該各ラジオボタン５４８，５４８，…の機能を表す“ヒストグラム”という文字列５５２が表示される。
【０１２６】
かかるサブ画面５００において、例えば、今、図２８に示すように、“選択範囲”という文字列５２２に付属する４つのラジオボタン５１８，５１８，…のうち、“リファレンスバー”という見出し５２０が付されたラジオボタン５１８がＯＮされるとする。そして、領域選択エリア５０８内のレントゲン撮影画像５０９上において、同図に点線５５４で示すように、マウス２２（図１参照）の操作によって上述したリファレンスバー３６（図２参照）に対応する部分５５６が選択されたとする。なお、この選択手順は、上述の第１実施形態（図１０）においてリファレンスバー３６に対応する部分２５０を選択するのと同様である。
【０１２７】
このように点線５５４で示される言わばリファレンスバー領域が選択されると、このリファレンスバー領域５５４内における輝度Ｙ［ｉ，ｊ］の最高値および最低値が、検出される。ここで、最高値としては、基本的に、リファレンスバー領域５５４のうち上述したアルミニウム箔４８が貼着された部分（図３参照）に対応する部分の輝度Ｙ［ｉ，ｊ］が検出される。一方、最低値としては、リファレンスバー領域５５４のうちリファレンスバー３６の最下段部分（厚さ寸法がＴａの部分）に対応する部分の輝度Ｙ［ｉ，ｊ］が検出される。
【０１２８】
これら検出された最高値および最低値は、それぞれ輝度最高値Ｙｂｍａｘおよび輝度最低値Ｙｂｍｉｎとして、図２９に示すテーブル６００に記録される。なお、このテーブル６００は、骨密度評価プログラムが起動されたときに、ハードディスク内に形成される。また、このテーブル６００には、上述した基準最高値Ｙｓｍａｘおよび基準最低値Ｙｓｍｉｎも記録されている。
【０１２９】
そして、このテーブル６００に記録された輝度最高値Ｙｂｍａｘ，輝度最低値Ｙｂｍｉｎ，基準最高値Ｙｓｍａｘおよび基準最低値Ｙｓｍｉｎに基づいて、領域選択エリア５０８に表示されている拡大画像５０９の輝度Ｙ［ｉ，ｊ］が補正される。つまり、次の数８に基づいて、補正後輝度Ｙ’［ｉ，ｊ］が求められる。
【０１３０】
【数８】
2005058751.tif 000010
【０１３１】
そして、この数８によって求められた補正後輝度Ｙ’［ｉ，ｊ］に基づいて、領域選択エリア５０８内の拡大画像５０９が改めて表示し直される。つまり、補正後輝度Ｙ’［ｉ，ｊ］に基づくレントゲン撮影画像１４０の拡大画像５０９が、当該領域選択エリア５０８に表示される。
【０１３２】
さらに、このレントゲン撮影画像１４０（拡大画像５０９）の全画素の補正後輝度Ｙ’［ｉ，ｊ］が、２５６階調に振り分けられると共に、上述した数１と同様の次の数９に基づいて、それぞれの階調ｘ毎の度数Ｈａ’［ｘ］が求められる。
【０１３３】
【数９】
2005058751.tif 000011
【０１３４】
なお、この数９において、Ｎａは、レントゲン撮影画像１４０の全画素数であり、ｎａ’［ｘ］は、当該全画素数Ｎａのうち階調ｘに振り分けられた画素数である。この数９の演算結果Ｈａ’［ｘ］もまた、上述のテーブル６００に記録される。
【０１３５】
そしてさらに、補正後輝度Ｙ’［ｉ，ｊ］に従うリファレンスバー領域５５４の拡大画像５５８が、リファレンスバー表示エリア５３６に表示される。そして、このリファレンスバー領域５５４についても、当該リファレンスバー領域５５４を構成する各画素の輝度Ｙｂ’［ｉ，ｊ］が２５６階調に振り分けられると共に、次の数１０に基づいて、それぞれの階調ｘ毎の補正後度数Ｈｂ’［ｘ］が求められる。
【０１３６】
【数１０】
2005058751.tif 000012
【０１３７】
この数１０において、Ｎｂは、リファレンスバー領域５５４を構成する全ての画素数であり、ｎｂ’［ｘ］は、当該全画素数Ｎｂのうち階調ｘに振り分けられた画素数である。この数１０の演算結果Ｈｂ’［ｘ］もまた、テーブル６００に記録される。
【０１３８】
この数１０によって求められた補正後度数Ｈｂ’［ｘ］に基づいて、さらにリファレンスバー領域５５４における補正後輝度Ｙｂ’［ｉ，ｊ］の平均値Ｍｂ’が求められる。具体的には、上述した数２と同様の次の数１１によって補正後平均値Ｍｂ’が求められる。
【０１３９】
【数１１】
2005058751.tif 000013
【０１４０】
この数１１の演算結果Ｍｂ’もまた、テーブル６００に記録される。そして、現在有効化されている画像が“画像１”である場合は、この補正後平均値Ｍｂ’は、後述する基準平均値ＳＭｂ’として設定される。
【０１４１】
さらに、補正後度数Ｈｂ’［ｘ］および補正後平均値Ｍｂ’に基づいて、リファレンスバー領域５５４における補正後輝度Ｙｂ’［ｉ，ｊ］の偏差Ｄｂ’が求められる。具体的には、上述した数３と同様の数１２によって補正後偏差Ｄｂ’が求められる。
【０１４２】
【数１２】
2005058751.tif 000014
【０１４３】
そして、この数１１の演算結果Ｄｂ’も、テーブル６００に記録される。また、現在有効化されている画像が“画像１”である場合は、この補正後偏差Ｄｂ’は、後述する基準偏差ＳＤｂ’として設定される。
【０１４４】
これらの補正後平均値Ｍｂ’および補正後偏差Ｄｂ’は、注記５２４として注記エリア５２６に表示される。また、この注記エリア５２６の上部には、当該注記５２４として、“リファレンスバー”という見出しを表す文字列も表示される。
【０１４５】
ここで、上述した“ヒストグラム”という文字列５５２に付属する３つのラジオボタン５４８，５４８，…のうち、“バー”という見出し５５０が付されたラジオボタン５４８がＯＮされるとする。すると、ヒストグラム表示エリア５４４に、リファレンスバー領域５５２の補正後度数Ｈｂ’［ｘ］に従う曲線（ヒストグラム）５６０が表示される。従って、この曲線５６０から、リファレンスバー領域５５２の補正後輝度Ｙｂ’［ｉ，ｊ］の分布を認識することができる。
【０１４６】
また、図には示さないが、“原画”という見出し５５０が付されたラジオボタン５４８がＯＮされると、ヒストグラム表示エリア５４４には、レントゲン撮影画像１４０全体の補正後度数Ｈａ’［ｘ］に従うヒストグラムが表示される。従って、この補正後度数Ｈａ’［ｘ］に従うヒストグラムから、レントゲン撮影画像１４０全体の補正後輝度Ｙ’［ｉ，ｊ］の分布を認識することができる。なお、この時点で、“対象”という見出し５５０が付されたラジオボタン５４８がＯＮされても、ヒストグラム表示エリア５４４には、何ら表示されない。
【０１４７】
次に、図３０に示すように、“選択範囲”という文字列５２２に付属する４つのラジオボタン５１８，５１８，…のうち、“最高値（白）”という見出し５２０が付されたラジオボタン５１８がＯＮされるとする。なお、当該“最高値（白）”というラジオボタン５１８がＯＮされることによって、上述の“リファレンスバー”というラジオボタン５１８がＯＦＦされる。これと同時に、リファレンスバー領域（点線）５５４の表示も消える。
【０１４８】
この図３０に示す状態において、領域選択エリア５０８内のレントゲン撮影画像５０９上の任意の位置でクリックされると、そのクリック位置の輝度Ｙ［ｉ，ｊ］が、輝度最高値Ｙｂｍａｘとして設定し直される。つまり、テーブル６００内の輝度最高値Ｙｂｍａｘが更新される。そして、この更新後の輝度最高値Ｙｂｍａｘは、上述の注記エリア５２６に注記５２４として表示される。
【０１４９】
さらに、更新後の輝度最高値Ｙｂｍａｘに基づいて、補正後輝度Ｙ’［ｉ，ｊ］が改めて補正される。そして、この言わば再補正後の輝度Ｙ’［ｉ，ｊ］に基づいて、レントゲン撮影画像１１４（５０９）全体の階調毎度数Ｈａ’［ｘ］，リファレンスバー領域５５４（拡大画像５５８）の階調毎度数Ｈｂ’［ｘ］，輝度平均値Ｍｂ’および輝度偏差Ｄｂ’が、改めて補正される。これらの補正結果は、領域選択エリア５０８に表示されたレントゲン撮影画像５０９，リファレンスバー表示エリア５３６に表示された拡大画像５５８，およびヒストグラム表示エリア５４４に表示された曲線５６０にも、反映される。
【０１５０】
なお、輝度最高値Ｙｂｍａｘとしては、上述したようにリファレンスバー領域５５４のうちアルミニウム箔４８が貼着された部分に対応する部分、つまり最も白く（明るく）撮影されている部分の輝度Ｙ’［ｉ，ｊ］が設定されるのが、望ましい。なぜなら、これ以外の部分が設定された場合は、正確な補正後輝度Ｙ’［ｉ，ｊ］を得られない可能性があるからである。ただし、例えばレントゲン撮影画像５０９（１１４）上においてアルミニウム箔４８に対応する部分が良好に撮影されておらず、かついずれかの歯２８，３０または３２に金属を用いた治療が施されている場合（例えば金冠または銀冠が被覆されている場合）には、当該金属に対応する部分の輝度Ｙ’［ｉ，ｊ］が輝度最高値Ｙｂｍａｘとして設定されるようにする。このようにすれば、アルミニウム箔４８に対応する部分が良好に撮影されていなくても、比較的に望ましい補正後輝度Ｙ’［ｉ，ｊ］を得ることができる。
【０１５１】
続いて、図３１に示すように、“選択範囲”に付属する４つのラジオボタン５１８，５１８，…のうち、“最低値（黒）”という見出し５２０が付されたラジオボタン５１８がＯＮされるとする。そして、この状態で、レントゲン撮影画像５０９上の任意の位置においてクリックされると、そのクリック位置の輝度Ｙ［ｉ，ｊ］が、輝度最低値Ｙｂｍｉｎとして設定し直される。つまり、テーブル６００内の輝度最低値Ｙｂｍｉｎが更新される。そして、この更新後の輝度最低値Ｙｂｍｉｎは、注記エリア５２６に注記５２４として表示される。
【０１５２】
さらに、上述の輝度最高値Ｙｂｍａｘが更新されたときと同様に、更新後の輝度最低値Ｙｂｍｉｎに基づいて、補正後輝度Ｙ’［ｉ，ｊ］が改めて補正される。そして、この再補正後輝度Ｙ’［ｉ，ｊ］に基づいて、階調毎度数Ｈａ’［ｘ］およびＨｂ’［ｘ］，輝度平均値Ｍｂ’並びに輝度偏差Ｄｂ’が、改めて補正されると共に、この補正結果は、領域選択エリア５０８内のレントゲン撮影画像５０９，リファレンスバー表示エリア５３６内の拡大画像５５８，およびヒストグラム表示エリア５４４内の曲線５６０に、反映される。
【０１５３】
なお、輝度最低値Ｙｂｍａｘとしては、上述したようにリファレンスバー領域５５４のうちリファレンスバー３６の最下段部分に対応する部分の輝度Ｙ［ｉ，ｊ］が設定されるのが、望ましい。ただし、例えば当該リファレンスバー３６の最下段部分が良好に撮影されていない場合には、レントゲン撮影画像５０９上で最も黒く（暗く）撮影されている部分の輝度Ｙ’［ｉ，ｊ］が設定されるようにする。このようにすれば、リファレンスバー３６の最下段部分が良好に撮影されていない場合でも、比較的に望ましい補正後輝度Ｙ’［ｉ，ｊ］を得ることができる。
【０１５４】
さらに続いて、図３２に示すように、上述した４つのラジオボタン５１８，５１８，…のうち、“評価対象”という見出し５２０が付されたラジオボタン５１８がＯＮされたとする。これによって、評価対象領域５６２の選択が開始される。
【０１５５】
即ち、領域選択エリア５０８内のレントゲン撮影画像５０９上において、第１実施形態（図１３参照）のときと同様に、マウス２２の操作によって評価対象領域５６２が選択される。ただし、第１実施形態では、評価対象領域２５４の形状が矩形に限定されたが、この第２実施形態においては、当該評価対象領域５６２を多角形状に選択することができる。具体的には、マウス２２による最初のクリック操作によって、評価対象領域５６２（枠）の始点が設定される。そして、マウス２２によってクリック操作が成される度に、評価対象領域５６２の頂点が設定される。そして、評価対象領域３６２の始点と同じ位置でマウス２２がダブルクリックされると、当該評価対象領域３６２の終点が設定される。従って、例えば矩形状の評価対象領域５６２を形成することもできるが、図３２に示すように、第１小臼歯２８の歯根４０に対応する部分５６４を当該評価対象領域５６２から外すこともできる。つまり、歯根４０の周囲の歯槽骨３４に対応する部分５６６のみを評価対象領域５６２として選択することができる。
【０１５６】
このようにして評価対象領域５６２が選択されると、サブ画面５００の上に、図３３に示すようなメッセージボックス７００がオーバレイ表示される。このメッセージボックス７００には、選択された評価対象領域５６２の確定の可否を問うメッセージ、例えば“この領域でキャプチャしてもよろしいですか？”という文字列７０２が表示される。そして、この文字列７０２の下方に、“はい”，“いいえ”および“キャンセル”という個別の文字列が付された３つのボタン７０４，７０６および７０８が、横一列に表示される。
【０１５７】
このメッセージボックス７００において、例えば“はい”というボタン７０４がクリックされると、評価対象領域５６２が確定される。そして、上述した数４と同様の次の数１３に基づいて、確定された評価対象領域５６２の補正後度数Ｈｏ’［ｘ］が求められる。
【０１５８】
【数１３】
2005058751.tif 000015
【０１５９】
なお、この数１３において、Ｎｏは、評価対象領域５６２を構成する全ての画素数であり、ｎｏ’［ｘ］は、当該全画素数Ｎｏのうち階調ｘの画素数である。この数１３の演算結果Ｈｏ’［ｘ］は、上述したテーブル６００に記録される。
【０１６０】
そして、この数１３によって求められた補正後度数Ｈｏ’［ｘ］に基づいて、評価対象領域５６２における補正後輝度Ｙｏ’［ｉ，ｊ］の平均値、つまり上述した補正後平均値Ｍｏ’が求められる。具体的には、上述した数５と同様の次の数１４によって、補正後平均値Ｍｏ’が求められる。
【０１６１】
【数１４】
2005058751.tif 000016
【０１６２】
この数１４の演算結果Ｍｏ’もまた、テーブル６００に記録される。この数１４によって求められた補正後平均値Ｍｏ’および上述の補正後度数Ｈｏ’［ｘ］に基づいて、評価対象領域５６２における補正後輝度Ｙｏ’［ｉ，ｊ］の偏差Ｄｏ’が求められる。具体的には、上述した数６と同様の数１５によって補正後偏差Ｄｏ’が求められる。なお、この数１５の演算結果Ｄｏ’もまた、テーブル６００に記録される。
【０１６３】
【数１５】
2005058751.tif 000017
【０１６４】
このようにして評価対象領域５６２の補正後度数Ｈｏ’［ｘ］，補正後平均値Ｍｏ’および補正後偏差Ｄｏ’が求められた後、メッセージボックス７００の表示が消える。そして、サブ画面５００が、図３４に示すような状態に遷移する。
【０１６５】
即ち、評価対象表示エリア５４０内に、評価対象領域５６２の拡大画像５６８が表示されると共に、注記エリア５２６に、補正後平均値Ｍｏ’および補正後偏差Ｍｏ’を表す文字列が、注記５２４として表示される。また、この注記エリア５２６の上部には、当該注記５２４として、“評価対象”という見出しを表す文字列も表示される。
【０１６６】
ここで、上述した“ヒストグラム”に付属する３つのラジオボタン５４８，５４８，…のうち、“対象”という見出し５５０が付された右側にあるラジオボタン５４８がＯＮされるとする。すると、ヒストグラム表示エリア５４４に、評価対象領域５６２の補正後度数Ｈｏ’［ｘ］に従う曲線５６０が表示される。従って、この曲線５６０から、評価対象領域５６２の補正後輝度Ｙｏ’［ｉ，ｊ］の分布を認識することができる。
【０１６７】
なお、“原画”というラジオボタン５４８がＯＮされると、上述したようにヒストグラム表示エリア５４４には、レントゲン撮影画像１４０全体の補正後度数Ｈａ’［ｘ］に従うヒストグラムが表示される。そして、“バー”というラジオボタン５４８がＯＮされると、リファレンスバー領域５５２の補正後度数Ｈｂ’［ｘ］に従うヒストグラムが表示される。
【０１６８】
また、評価対象領域５６２が確定された後で、上述の如く輝度最高値Ｙｂｍａｘが更新されると、この更新後の輝度最高値Ｙｂｍａｘに基づいて、補正後輝度Ｙ’［ｉ，ｊ］が改めて補正され、ひいては当該評価対象領域５６２の補正後度数Ｈｏ’［ｘ］，補正後平均値Ｍｏ’および補正後偏差Ｄｏ’も改めて補正される。そして、これらの補正結果は、評価対象表示エリア５４０の拡大画像５６８、およびヒストグラム表示エリア５４４の曲線５６０にも、反映される。このことは、輝度最低値Ｙｂｍｉｎが更新されたときも、同様である。
【０１６９】
図３３に戻って、メッセージボックス７００において、“いいえ”というボタン７０６がクリックされると、当該メッセージボックス７００の表示が消える。この場合、図３２のサブ画面５００が再度現れる。そして、このサブ５００において、評価対象領域５６２の形状および寸法を任意に変更することができる。この変更手順については、一般の画像編集ソフトウェアによる手順と同様であるので、詳しい説明を省略する。
【０１７０】
また、メッセージボックス７００において“キャンセル”というボタン７０８がクリックされると、当該メッセージボックス７００の表示が消えると共に、サブ画面５００が最初の状態、つまり評価対象領域５６２の始点が設定される前の状態に戻る。なお、メッセージボックス７００の右上にある“×”印の閉じるボタン７１０がクリックされた場合も、同様である。
【０１７１】
さらに、サブ画面５００において、上述した“基準最高値”という見出し５３２が付されたエディットボックス５２８の数値が変更されると、これに従ってテーブル６００内の基準最高値Ｙｓｍａｘが更新される。そして、この更新後の基準最高値Ｙｓｍａｘに基づいて、補正後輝度Ｙ’［ｉ，ｊ］が改めて補正され、ひいてはテーブル６００内の各パラメータ、つまり各補正後度数Ｈａ’、Ｈｂ’，Ｈｏ’，補正後平均値Ｍｂ’，Ｍｏ’，補正後偏差Ｄｂ’およびＤｏ’も改めて補正される。そして、これらの補正結果は、領域選択エリア５０８のレントゲン撮影画像５０９，評価対象表示エリア５４０の拡大画像５６８、およびヒストグラム表示エリア５４４の曲線５６０にも、反映される。このことは、“基準最低値”という見出し５３４が付されたエディットボックス５３０の数値が変更されたときも、同様である。ただし、これら基準最高値Ｙｓｍａｘおよび基準最低値Ｙｓｍｉｎは、飽くまで基準値であるので、不用意に変更されるのは好ましくない。これらが変更されるのは、例えば上述の如くリファレンスバー領域５５４が良好に撮影されていない場合等のように輝度最高値Ｙｂｍａｘまたは輝度最低値Ｙｂｍｉｎを正確に捉えることができない場合に留めておくのが、望ましい。
【０１７２】
なお、この第２実施形態では、上述したように評価対象領域５６２を多角形状に選択することができるが、その反面、当該評価対象領域５６２として一定の形状でかつ一定寸法の領域を設定し難くなる場合がある。そこで、この第２実施形態においては、一定形状かつ一定寸法の評価対象領域５６２を設定できるようにするべく、当該評価対象領域５６２の設定を支援する領域設定支援機能を備えている。
【０１７３】
即ち、“参照チェック”という見出し５１４が付されたチェックボックス５１６がクリックされると、図３５に示すように、当該チェックボックス５１６にチェックマーク５７０が入力される。これによって、上述の領域設定支援機能が有効化される。そして、“領域抽出”という文字列５１４の下方にある４つのラジオボタン５１２，５１２，…のいずれか、例えば“上”という文字列５１０が付されたラジオボタン５１２がＯＮされるとする。この状態で、レントゲン撮影画像５０９上の第１小臼歯２８の最上端に対応する部分がクリックされると、当該部分に所定のマーク、例えば“×”というマーク５７２が付される。
【０１７４】
これと同様に、例えば“下”，“左”および“右”のラジオボタン５１２，５１２，…がこの順番でＯＮされると共に、レントゲン撮影画像５０９上の第１小臼歯２８の最下端，左端および右端のそれぞれに対応する部分がこの順番でクリックされると、図３６に示すように、当該第１小臼歯２８の最下端，左端および右端に対応する部分のそれぞれに、マーク５７２が付される。そして、このように“上”，“下”，“左”および“右”に対応する４つのマーク５７２，５７２，…が付されると、これらのマーク５７２，５７２，…を各辺に１つずつ置く矩形の目印線５７４が、自動的に描画される。さらに、この目印線５７４の中央に、当該目印線５７４で表される矩形を半分に縮小した別の目印線５７６が描画されると共に、互いに対向するマーク５７２，５７２，…を結ぶ直線状の目印線５７８および５８０が描画される。
【０１７５】
従って、図３７に示すように、これらの目印線５７４〜５８０を参考にして、評価対象領域５６２を描画する、つまり始点，各頂点および終点を指定すれば、当該評価対象領域５６２を設定し易くなり、ひいては一定形状かつ一定寸法の評価対象領域５６２を設定できるようになる。なお、ラジオボタン５１２，５１２，…をＯＮする順番は、任意である。また、それぞれのマーク５７２の位置は、任意に変更することができる。そして、上述したチェックマーク５７０が外されると（チェックボックス５１６が再度クリックされると）、各マーク５７２，５７２，…および目印線５７４〜５８０の表示が消える。
【０１７６】
このようにサブ画面５００による補正を終えた後、タイトルバー５０２の右端にある閉じるボタン５８２がクリックされると、ディスプレイ２４の表示は、図３８に示すようなメイン画面１００に戻る。
【０１７７】
即ち、補正後平均値Ｍｏ’および補正後偏差Ｄｏ’が、“画像１”のピクチャボックス１１６の下方の文字列１２０および１２２によって表される。そして、上述した基準平均値ＳＭｂ’および基準偏差ＳＤｂ’が、メイン画面１００の右端部分に表示された文字列１３０および１３２によって表される。さらに、“チェック”という文字列４１０の横方に、補正後平均値Ｍｏ’に基づく判定結果が表示される。ここでは、補正後平均値Ｍｏ’（＝８３．３３５１）が、“要注意”というエディットボックス４０２に入力された値（＝７０）よりも大きく、かつ“注意”というエディットボックス４０４に入力された値以下であるので、当該判定結果として骨減少症の可能性があることを示す“注意”という結果が表示される。なお、補正後平均値Ｍｏ’が、“要注意”というエディットボックス４０２に入力された値以下である場合には、当該判定結果として骨粗鬆症の可能性があることを示す“要注意”という結果が表示される。また、補正後平均値Ｍｏ’が、“注意”というエディットボックス４０４に入力された値よりも大きい場合には、当該判定結果として“正常”という結果が表示される。
【０１７８】
そしてさらに、ヒストグラム表示エリア１３４に、補正後度数Ｈｏ’［ｘ］に従う曲線４１２が表示される。この曲線４１２は、上述した図３４における曲線５６０と相似形のものである。なお、図３８からは判らないが、“画像１”のピクチャボックス１１６には、補正後輝度Ｙ’［ｉ，ｊ］に基づく（つまり補正後の）レントゲン撮影画像１４０が表示される。
【０１７９】
これと同じ要領で、他の“画像Ａ”、つまり“画像２”，“画像３”および“画像４”についてもサブ画像５００が表示され、当該サブ画像５００上でリファレンスバー領域５５４と評価対象領域５６２とが選択されると、当該サブ画像５００は図３９に示すようになる。なお、図３９は、“画像２”についてのサブ画像５００である。
【０１８０】
この“画像１”以外のサブ画面５００では、当該“画像１”とは一部異なる方法によって、補正処理が行われる。即ち、リファレンスバー領域５５４が選択されると、上述した数１，数２および数３に基づいて当該リファレンスバー領域５５４の補正前の度数Ｈｂ［ｘ］，輝度平均値Ｍｂおよび輝度偏差Ｄｂが求められる。これらの演算結果Ｈｂ［ｘ］，ＭｂおよびＤｂは、上述したテーブル６００に記録される。そして、これら補正前の度数Ｈｂ［ｘ］，輝度平均値Ｍｂおよび輝度偏差Ｄｂに基づいて、レントゲン撮影画像１４０（拡大画像５０９）の輝度Ｙ［ｉ，ｊ］が補正される。具体的には、上述した数７と同様の次の数１６に基づいて、補正後輝度Ｙ’［ｉ，ｊ］が求められる。
【０１８１】
【数１６】
2005058751.tif 000018
【０１８２】
つまり、“画像１”以外の“画像Ａ”については、当該“画像１”のリファレンスバー領域５５４の補正後輝度Ｙｂ’［ｉ，ｊ］に基づいて輝度補正が行われる。なお、これ以降は、“画像１”の場合と同様の手順により補正処理が行われる。
【０１８３】
このようにして全ての“画像Ａ”についてサブ画面５００による補正処理が行われた後、メイン画面４００に戻されると、当該メイン画面４００は図４０に示すようになる。
【０１８４】
即ち、各“画像Ａ”の補正後平均値Ｍｏ’および補正後偏差Ｄｏ’が、それぞれに対応するピクチャボックス１１６の下方の文字列１２０および１２２によって表示される。そして、補正後平均値Ｍｏ’に基づく判定結果が、文字列４１０によって表される。また、ピクチャボックス１１６内のレントゲン撮影画像１４０は、補正後輝度Ｙ’［ｉ，ｊ］に基づいて表示し直される。さらに、ヒストグラム表示エリア１３４に、各“画像Ａ”の補正後度数Ｈｏ’［ｘ］に従う４つの曲線４１２，４１２，…が互いに異なる色で表示される。
【０１８５】
このように、第２実施形態によれば、基準最高値Ｙｓｍａｘおよび基準最低値Ｙｓｍｉｎという一定の基準に基づいて、“画像１”の輝度Ｙ［ｉ，ｊ］が補正される。そして、この一定の基準に基づいて補正された“画像１”のリファレンスバー領域５５４の補正後輝度Ｙｂ’［ｉ，ｊ］に基づいて、他の“画像Ａ”、つまり“画像２”，“画像３”および“画像４”の輝度補正が行われる。換言すれば、これらの“画像２”，“画像３”および“画像４”についても、一定の基準に基づいて輝度補正が行われることになる。従って、各“画像Ａ”の補正後輝度Ｙ’［ｉ，ｊ］、厳密には補正後平均値Ｍｏ’から、骨密度の程度を判定することができる。しかも、この判定結果がメイン画面４００上で（文字列４１０によって）表示されるので、患者側にとっても自身の骨密度を把握し易くなる。勿論、この判定結果は、印刷したり、記録媒体に記録したりすることもできる。
【図面の簡単な説明】
【０１８６】
【図１】この発明の第１実施形態を示す骨密度評価システムの全体構成図である。
【図２】同第１実施形態における撮影後のレントゲンフィルムの一態様を示す図解図である。
【図３】同第１実施形態におけるアルミニウム製ブロックの概観斜視図である。
【図４】同第１実施形態における撮影前のレントゲンフィルムの一態様を示す図解図である。
【図５】同第１実施形態においてレントゲン撮影時の状態を示す図解図である。
【図６】同第１実施形態においてＰＣのディスプレイに表示されるメイン画面の一態様を示す図解図である。
【図７】図６とは異なるメイン画面の態様を示す図解図である。
【図８】図７とはさらに異なるメイン画面の態様を示す図解図である。
【図９】同第１実施形態においてＰＣのディスプレイに表示されるサブ画面の一態様を示す図解図である。
【図１０】同第１実施形態におけるサブ画面の一部分を拡大して示す図解図である。
【図１１】図９とは異なるサブ画面の態様を示す図解図である。
【図１２】同第１実施形態におけるテーブルの態様を概念的に示す図解図である。
【図１３】図１０とは異なる態様を示す拡大部分の図解図である。
【図１４】図１１とはさらに異なるサブ画面の態様を示す図解図である。
【図１５】図８とはさらに異なるメイン画面の態様を示す図解図である。
【図１６】図１４とはさらに異なるサブ画面の態様を示す図解図である。
【図１７】図１５とはさらに異なるメイン画面の態様を示す図解図である。
【図１８】同第１実施形態におけるＰＣの動作を示すフローチャートである。
【図１９】図１８における領域選択処理の詳細を示すフローチャートである。
【図２０】図１９に続くフローチャートである。
【図２１】図１９に続く図２０とは別のフローチャートである。
【図２２】同第１実施形態によって得られた補正後輝度と実際のＢＭＤ値およびそのＴ値との比較を示す一覧表である。
【図２３】図２２におけるＴ値の所定区分と当該区分毎の補正後輝度の平均値との比較を示す表である。
【図２４】図２３の関係を表すグラフである。
【図２５】この発明の第２実施形態においてＰＣのディスプレイに表示されるメイン画面の一態様を示す図解図である。
【図２６】図２５とは異なるメイン画面の態様を示す図解図である。
【図２７】同第２実施形態においてＰＣのディスプレイに表示されるサブ画面の一態様を示す図解図である。
【図２８】図２７とは異なるサブ画面の態様を示す図解図である。
【図２９】同第２実施形態におけるテーブルの態様を概念的に示す図解図である。
【図３０】図２８とはさらに異なるサブ画面の態様を示す図解図である。
【図３１】図３０とはさらに異なるサブ画面の態様を示す図解図である。
【図３２】図３１とはさらに異なるサブ画面の態様を示す図解図である。
【図３３】同第２実施形態においてＰＣのディスプレイに表示されるメッセージボックスの一態様を示す図解図である。
【図３４】図３２とはさらに異なるサブ画面の態様を示す図解図である。
【図３５】図３４とはさらに異なるサブ画面の態様を示す図解図である。
【図３６】図３５とはさらに異なるサブ画面の態様を示す図解図である。
【図３７】図３６とはさらに異なるサブ画面の態様を示す図解図である。
【図３８】図２６とはさらに異なるメイン画面の態様を示す図解図である。
【図３９】図３７とはさらに異なるサブ画面の態様を示す図解図である。
【図４０】図３８とはさらに異なるメイン画面の態様を示す図解図である。
【符号の説明】
【０１８７】
１０  骨密度評価システム
１２  ＰＣ
１４  フィルムスキャナ
１６  レーザプリンタ
１００  メイン画面
２００  サブ画面
２１４  レントゲン撮影画像
２４２  ヒストグラム表示エリア
２４８  リファレンスバー領域
２５０  リファレンスバー
２５４  評価対象領域
【図１】
2005058751.tif 000019
【図２】
2005058751.tif 000020
【図３】
2005058751.tif 000021
【図４】
2005058751.tif 000022
【図５】
2005058751.tif 000023
【図６】
2005058751.tif 000024
【図７】
2005058751.tif 000025
【図８】
2005058751.tif 000026
【図９】
2005058751.tif 000027
【図１０】
2005058751.tif 000028
【図１１】
2005058751.tif 000029
【図１２】
2005058751.tif 000030
【図１３】
2005058751.tif 000031
【図１４】
2005058751.tif 000032
【図１５】
2005058751.tif 000033
【図１６】
2005058751.tif 000034
【図１７】
2005058751.tif 000035
【図１８】
2005058751.tif 000036
【図１９】
2005058751.tif 000037
【図２０】
2005058751.tif 000038
【図２１】
2005058751.tif 000039
【図２２】
2005058751.tif 000040
【図２３】
2005058751.tif 000041
【図２４】
2005058751.tif 000042
【図２５】
2005058751.tif 000043
【図２６】
2005058751.tif 000044
【図２７】
2005058751.tif 000045
【図２８】
2005058751.tif 000046
【図２９】
2005058751.tif 000047
【図３０】
2005058751.tif 000048
【図３１】
2005058751.tif 000049
【図３２】
2005058751.tif 000050
【図３３】
2005058751.tif 000051
【図３４】
2005058751.tif 000052
【図３５】
2005058751.tif 000053
【図３６】
2005058751.tif 000054
【図３７】
2005058751.tif 000055
【図３８】
2005058751.tif 000056
【図３９】
2005058751.tif 000057
【図４０】
2005058751.tif 000058