(57)【要約】【課題】骨塩量等の骨に関する計測データを精度良く算出すること。【解決手段】この放射線画像検査方法は、基準体の像を含む画像データを用いて骨に関する計測を行う放射線画像検査方法であって、輝度値分布取得手段が、画像データに基づいて、長手方向に沿った基準体の像の輝度値の分布を得る工程（Ｓ１０１）と、第１及び第２画素取得手段が、前記長手方向における基準体の像の一端と他端との間の２つのエッジを抽出して、２つのエッジのそれぞれにおいて長手方向の連続画素数である第１及び第２画素数を得る工程（Ｓ１０３、Ｓ１０４）と、配置形態判定手段が、輝度値の分布、第１画素数及び第２画素数に基づいて、基準体の配置形態の正誤を判定する工程（Ｓ１０５）と、を含む。【選択図】図４
2012231996.tif 000002
【特許請求の範囲】
【請求項１】
基準体の像を含む画像データを用いて骨に関する計測を行う放射線画像検査方法であって、
輝度値分布取得手段が、前記画像データに基づいて、長手方向に沿った前記基準体の像の輝度値の分布を得る工程と、
第１及び第２画素数取得手段が、前記長手方向における前記基準体の像の一端と他端との間の２つのエッジを抽出して、前記２つのエッジのそれぞれにおいて前記長手方向の連続画素数である第１及び第２画素数を得る工程と、
配置形態判定手段が、前記輝度値の分布、前記第１画素数及び前記第２画素数に基づいて、前記基準体の配置形態の正誤を判定する工程と、
を含むことを特徴とする放射線画像検査方法。
【請求項２】
輝度値補正手段が、前記正誤を判定する工程により前記基準体の配置形態が誤と判定されたとき、前記輝度値の分布を補正する工程、
をさらに含むことを特徴とする請求項１記載の放射線画像検査方法。
【請求項３】
基準体の像を含む画像データを用いて骨に関する計測を行う放射線画像検査プログラムであって、
コンピュータを、
前記画像データに基づいて、長手方向に沿った前記基準体の像の輝度値の分布を得る輝度値分布取得手段、
前記長手方向における前記基準体の像の一端と他端との間の２つのエッジを抽出して、前記２つのエッジのそれぞれにおいて前記長手方向の連続画素数である第１及び第２画素数を得る第１及び第２画素数取得手段、
前記輝度値の分布、前記第１画素数及び前記第２画素数に基づいて、前記基準体の配置形態の正誤を判定する配置形態判定手段、
として機能させることを特徴とする放射線画像検査プログラム。
【請求項４】
基準体の像を含む画像データを用いて骨に関する計測を行う放射線画像検査装置であって、
前記画像データに基づいて、長手方向に沿った前記基準体の像の輝度値の分布を得る輝度値分布取得手段と、
前記長手方向における前記基準体の像の一端と他端との間の２つのエッジを抽出して、前記２つのエッジのそれぞれにおいて前記長手方向の連続画素数である第１及び第２画素数を得る第１及び第２画素数取得手段と、
前記輝度値の分布、前記第１画素数及び前記第２画素数に基づいて、前記基準体の配置形態の正誤を判定する配置形態判定手段と、
を含むことを特徴とする放射線画像検査装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線画像データを対象に骨に関する計測を行う放射線画像検査方法、放射線画像検査プログラム、及び放射線画像検査装置に関する。
【背景技術】
【０００２】
人間の骨の発育状態等を確認するために、骨塩量の計測が行われている。この骨塩量を計測する方法として、ＭＤ法やＤＩＰ法などが知られている。ＭＤ法では、人間の両手の中央にアルミスケールを置いてＸ線撮影を行う。このようにして撮影したＸ線フィルムを専用の解析装置により解析して骨塩量を得る。解析は、第２中手骨の中央をマイクロデンシトメータで走査し、得られた濃度をアルミスケール濃度に換算して、骨塩量に関する種々のデータを得る。
【０００３】
また、ＤＩＰ法では、アルミスケールと人間の両手とを同時に撮影した撮影済みＸ線フィルムを画像データに変換する。解析は、画像処理装置を用いてＭＤ法と同様の解析を自動的に行う。ＤＩＰ法は、ＭＤ法と比較して計測精度は落ちるものの解析時間を大幅に短縮することができる。
【０００４】
特許文献１には、被検骨の医用画像を利用して骨塩量計測を行うＭＤ法による骨計測装置が記載されている。この計測装置では、被検骨の骨管の方向を定量的に計測することにより、計測値の再現性が向上するように骨計測位置を決定することができる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００７−１５１６０９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１に記載されたＭＤ法による骨計測装置などでは、骨塩量に関する種々のデータに変換するための変換係数を画像データから生成するために、アルミスロープ等の基準体と共に人間の両手をＸ線撮影する。基準体の場所ごとの厚みは既知であるので、画像データの所定の場所における基準体の厚みと輝度値との関係が得られる。この関係に基づいて変換係数を生成する。しかしながら、Ｘ線撮影するとき誤った配置形態で基準体が配置されると、画像データの所定の場所における基準体の厚みと輝度値との関係を正確に得ることができない。これ故、骨塩量に関するデータを精度よく計測できないおそれがある。
【０００７】
本発明は、このような問題点に鑑みて為されたものである。本発明は、骨塩量等の骨に関する計測データを精度良く算出することができる放射線画像検査方法、放射線画像検査プログラム、及び放射線画像検査装置を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記問題を解決するため、本発明の放射線画像検査方法は、基準体の像を含む画像データを用いて骨に関する計測を行う放射線画像検査方法であって、輝度値分布取得手段が、画像データに基づいて、長手方向に沿った基準体の像の輝度値の分布を得る工程と、第１及び第２画素取得手段が、長手方向における基準体の像の一端と他端との間の２つのエッジを抽出して、２つのエッジのそれぞれにおいて長手方向の連続画素数である第１及び第２画素数を得る工程と、配置形態判定手段が、輝度値の分布、第１画素数及び第２画素数に基づいて、基準体の配置形態の正誤を判定する工程と、を含む。
【０００９】
或いは、本発明の放射線画像検査プログラムは、基準体の像を含む画像データを用いて骨に関する計測を行う放射線画像検査プログラムであって、コンピュータを、画像データに基づいて、長手方向に沿った基準体の像の輝度値の分布を得る輝度値分布取得手段、長手方向における基準体の像の一端と他端との間の２つのエッジを抽出して、２つのエッジのそれぞれにおいて長手方向の連続画素数である第１及び第２画素数を得る第１及び第２画素取得手段、輝度値の分布、第１画素数及び第２画素数に基づいて、基準体の配置形態の正誤を判定する配置形態判定手段、として機能させる。
【００１０】
或いは、本発明の放射線画像検査装置は、基準体の像を含む画像データを用いて骨に関する計測を行う放射線画像検査装置であって、画像データに基づいて、長手方向に沿った基準体の像の輝度値の分布を得る輝度値分布取得手段と、長手方向における基準体の像の一端と他端との間の２つのエッジを抽出して、２つのエッジのそれぞれにおいて長手方向の連続画素数である第１及び第２画素数を得る第１及び第２画素取得手段と、輝度値の分布、第１画素数及び第２画素数に基づいて、基準体の配置形態の正誤を判定する配置形態判定手段と、を含む。
【００１１】
このような放射線画像検査方法、放射線画像検査プログラム、或いは放射線画像検査装置によれば、画像データに基づいて長手方向に沿った基準体の像の輝度値の分布が得られる。さらに、長手方向における基準体の像の一端と他端との間の２つのエッジを抽出して、２つのエッジのそれぞれにおいて長手方向の連続画素数である第１及び第２画素数を得る。そして、輝度値の分布、第１画素数及び第２画素数に基づいて、基準体の配置形態の正誤を判定する。この判定結果に基づいて骨塩量等の骨に関する計測データを求めるための変換係数を得ることにより、計測データを精度良く算出できる変換係数を得ることができる。
【００１２】
本発明の放射線画像検査方法は、輝度値補正手段が、正誤を判定する工程により基準体の配置形態が誤と判定されたとき、輝度値の分布を補正する工程、をさらに含んでいてもよい。
【００１３】
このような構成によれば、誤った配置形態で撮影された基準体の像から得た輝度値の分布を補正する。この補正された輝度値の分布に基づいて骨塩量等の骨に関する計測データを求めるための変換係数を得ることにより、計測データを精度良く算出できる変換係数を得ることができる。
【発明の効果】
【００１４】
本発明によれば、骨塩量に関するデータを精度良く算出することができる変換係数を得ることができる。
【図面の簡単な説明】
【００１５】
【図１】本実施形態に用いられる基準体を説明するための図である。
【図２】本実施形態に用いられる放射線画像検査装置の概略構成を示す機能ブロック図である。
【図３】人間の両手と基準体とを撮影した画像データの一例を示す図である。
【図４】本実施形態に係る放射線画像検査方法の主要なステップを説明するための図である。
【図５】第１画素数を得る工程を説明するための図である。
【図６】第２画素数を得る工程を説明するための図である。
【図７】正しい配置形態の一例を示す図である。
【図８】正しい配置形態の別の例を示す図である。
【図９】誤った配置形態の一例を示す図である。
【図１０】誤った配置形態の別の例を示す図である。
【図１１】配置形態の正誤を判定するための主要なステップを説明するための図である。
【図１２】輝度値の分布を補正する工程を説明するための図である。
【図１３】記録媒体に記録されたプログラムを実行するためのコンピュータのハードウェア構成を示す図である。
【図１４】記録媒体に記録されたプログラムを実行するためのコンピュータの斜視図である。
【発明を実施するための形態】
【００１６】
以下、添付図面を参照しながら本発明による放射線画像検査方法、放射線画像検査プログラム、及び放射線画像検査装置の実施の形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。また、図１、図３、図５〜図１０、図１２に示すようにＸＹＺ座標系を設定し、以下の説明に用いる。以下の各実施形態において、ＸＹＺ座標系のＸ方向が画像データにおける水平方向であり、Ｙ方向が画像データにおける垂直方向であり、Ｚ方向が基準体の厚さ方向である。
【００１７】
本実施形態に用いられる基準体について説明する。図１は、本実施形態に用いられる基準体を説明するための斜視図である。本実施形態では、基準体はアルミニウム製のアルミスロープ１である。アルミスロープ１は、所定の傾きを有する斜面２を有している。本実施形態では、このアルミスロープ１を写した画像データを参照することによってどの輝度値がアルミ厚さに換算すると何ｍｍの厚さに相当するかが計算され、その計算結果により変換係数が生成される。
【００１８】
アルミスロープ１は、略楔形の形状を有している。すなわち、アルミスロープ１は、斜面２、先端面３、後端面４、基準側面５、平側面６、及び底面７を有している。斜面２、基準側面５、平側面６及び底面７を介して先端面３は後端面４に繋がっている。先端面３から後端面４までの距離Ｌ１は、例えば１５０ｍｍである。基準側面５から平側面６までの距離Ｌ２は、例えば１５ｍｍである。先端面３におけるアルミスロープ１の高さＬ３は、例えば０．７ｍｍである。後端面４におけるアルミスロープ１の高さＬ４は、例えば１５．４ｍｍである。斜面２は、＋Ｙ方向（先端面３から後端面４に向かう長手方向）に沿って厚さが増加する傾きを有する。例えば、Ｙ方向に１０ｍｍ進行すると、Ｚ方向に１ｍｍ進行するような傾きを有する。即ち、Ｙ方向の長さとＺ方向の最大厚との比が１０：１であるような傾きを有する。
【００１９】
基準側面５は、稜線５ａ及び複数の基準溝５ｂを有している。基準溝５ｂは、Ｙ方向に間隔を空けて設けられ、Ｚ方向に沿って斜面２から底面７まで延びている。先端面３からＹ方向に距離Ｌ５離間した位置に一の基準溝５ｂが設けられている。例えば、距離Ｌ５は１０ｍｍである。一の基準溝５ｂからＹ方向に１０ｍｍ離間した位置に隣接する別の基準溝５ｂが設けられている。即ち、隣り合う基準溝５ｂの間隔Ｌ６は１０ｍｍである。本実施形態では、基準溝５ｂはＹ方向に沿って１４個形成されている。基準溝５ｂのＸＹ平面に沿った断面は、略三角形の形状を有している。例えば、基準溝５ｂのＹ方向の幅Ｌ７は２ｍｍである。
【００２０】
平側面６は、稜線６ａを有しており、基準側面５の反対側に位置している。平側面６には、基準溝５ｂは設けられていない。また、底面７は、先端面３と直交するように形成されている。また、底面７は後端面４と直交するように形成されている。
【００２１】
図２は、本発明の好適な一実施形態に係る放射線画像検査装置である放射線画像検査装置５０の概略構成を示す機能ブロック図である。同図に示す放射線画像検査装置５０は、放射線透過画像を基に得られた画像データを取得し、その画像データを参照して骨塩量を測定するための情報処理装置である。この画像データは、上述したアルミスロープ１と検査対象者の手とを同時に放射線撮影した放射線透過画像をビットマップ形式やＪＰＥＧ形式等のデジタルデータに変換したものである。画像データは、通信ネットワークやＣＤ−ＲＯＭ、ＤＶＤ、半導体メモリ等の記録媒体を介して放射線画像検査装置５０に入力される。ここで、画像データは、放射線透過画像が写された放射線フィルムをデジタルカメラで撮影したり、専用の読取装置（図示せず）を用いて、放射線透過画像が記録されたイメージングプレートにレーザ照射をして輝尽性蛍光像を得ることによって生成される。
【００２２】
この放射線画像検査装置５０は、図２に示すように、機能的構成要素として、輝度値分布取得部（輝度値分布取得手段）５１、輝度値傾き算出部（輝度値傾き算出手段）５２、第１画素数取得部（第１画素数取得手段）５３ａ、第２画素数取得部（第２画素数取得手段）５３ｂ、配置形態判定部（配置形態判定手段）５４、輝度値補正部（輝度値補正手段）５５、測定処理部５６、及び計測値計算部５７を備える。
【００２３】
輝度値分布取得部５１は、アルミスロープ１の像を含む画像データに基づいて、アルミスロープ１の像における輝度値の分布を得る。具体的には、画像データを所定の閾値で二値化して、二値化された画像データを生成する。この画像データからアルミスロープ１の像を抽出する。そして、アルミスロープ１の像の領域と二値化する前の画像データとを用いて、アルミスロープ１の像における輝度値の分布を得る。
【００２４】
輝度値傾き算出部５２は、輝度値分布取得部５１によって取得された輝度値の分布に基づいて、Ｙ方向に沿った輝度値の傾きを算出する。
【００２５】
第１画素数取得部５３ａ及び第２画素数取得部５３ｂは、第１画素値及び第２画素値をカウントする。まず、画像データにおいて、アルミスロープ１の像の一端と他端とを繋ぐ一のエッジ及び他のエッジを抽出する。次に、画像データを所定の閾値で二値化して、Ｘ方向の所定幅の範囲おいて画素値“１”を有し、且つＹ方向に連続する画素数をカウントする。第１画素値は一のエッジにおいて連続する画素数である。第２画素値は他のエッジにおいて連続する画素数である。
【００２６】
配置形態判定部５４は、輝度値傾き算出部５２によって算出された輝度値の傾き、第１画素数取得部５３ａ及び第２画素数取得部５３ｂによってカウントされた第１画素数及び第２画素数に基づいて、アルミスロープ１の配置形態の正誤を判定する。
【００２７】
輝度値補正部５５は、配置形態判定部５４によってアルミスロープ１の配置形態が誤りであると判定された場合に、輝度値の分布を補正する。
【００２８】
測定処理部５６は、画像データに基づいて、所定の解析領域を設定するためのパラメータを算出する。具体的には、画像データに写された骨において、解析領域を設定するために必要な複数の基準位置を特定する。
【００２９】
計測値計算部５７は、画像データにおける所定の解析領域を設定する。次に、計測値計算部１９は、所定の解析領域における画素値を取得し、さらに上述した輝度値の分布に基づいて生成した変換係数を用いて該画素値を骨塩量評価に用いる骨密度パラメータとして算出する。そして、計測値計算部１９は、算出した骨密度パラメータをディスプレイやプリンタ等の外部出力デバイスや外部データに出力する。
【００３０】
図３は、放射線画像検査装置５０の処理対象となる画像データＤ１の一例を示す図である。画像データＤ１には、人間の左手８Ａ、右手８Ｂの像が表示されている。アルミスロープ１の像は、左手８Ａの像と右手８Ｂの像との間に表示されている。画像データＤ１において、明るく見える領域は輝度値が高い領域である。輝度値が高い領域は、例えばアルミスロープ１の厚みが大きいことを示す。暗く見える領域は輝度値が低い領域である。輝度値が低い領域は、例えばアルミスロープ１の厚みが小さいことを示す。
【００３１】
次に、本実施形態の放射線画像検査装置５０によって実行される骨塩量測定方法を説明すると共に、本実施形態の放射線画像検査方法について詳細に説明する。ここでは、図３に示された画像データＤ１から骨塩量を計算するための変換係数を生成するための輝度値の分布を得る工程を説明する。図４は、本実施形態に係る放射線画像検査方法の主要なステップを説明するための図である。
【００３２】
工程Ｓ１０１では、輝度値分布取得部５１によってアルミスロープ１の像における所定方向の輝度値の分布を得る。本実施形態では、所定方向は＋Ｙ方向である。まず、画像データＤ１を所定の閾値で二値化して、二値化画像を生成する。二値化された画像において、アルミスロープ１の像は例えば画素値“１”を有する。そして、二値化された画像データにおいて、Ｘ方向又はＹ方向に走査してアルミスロープ１の像を抽出して、アルミスロープ１の像の領域を示す座標データを得る。次に、アルミスロープ１の像の領域を示す座標データと画像データＤ１とを用いて、アルミスロープ１の像におけるＸ方向の輝度値の平均値を各Ｙ座標毎に算出する。上述した処理により、所定のＹ座標における輝度値が得られるので、Ｙ方向における輝度値の分布が得られる。
【００３３】
工程Ｓ１０２では、輝度値傾き算出部５２によって輝度値の分布に基づいて、Ｙ方向における輝度値の傾きを算出する。例えば、Ｙ方向に沿った距離を変数とした１次関数を用いて輝度値の分布を近似して、輝度値の分布の傾きを算出する。
【００３４】
工程Ｓ１０３では、第１画素数取得部５３ａによって第１画素値を得る。まず、画像データＤ１のアルミスロープ１を含む任意の範囲において、Ｙ方向におけるアルミスロープ１の像の一端１ａと他端１ｂとを繋ぐ一のエッジ１ｃを抽出する（図３参照）。図５の（ａ）部は、アルミスロープ１の像における左側のエッジ１ｃが抽出された画像データＤ２である。このエッジ１ｃは、基準側面５に含まれる稜線５ａ又は平側面６に含まれる稜線６ａの何れか一方であるが、この時点ではどちらであるか不明である。エッジ１ｃの抽出には、ソーベルフィルタを用いることができる。なお、エッジ１ｃの抽出には、ブリューイットフィルタ、ロバーツフィルタを用いても良い。
【００３５】
次に、画像データＤ２を所定の閾値で二値化して、図５の（ｂ）部に示される画像データＤ３を得る。画像データＤ３において、抽出されたエッジ１ｃでは、画素値“１”を有する複数の画素同士がＹ方向に沿って互いに隣接している。このように互いに隣接している画素の数をカウントする。即ち、画像データＤ３において、Ｘ方向の所定幅の範囲に画素値“１”を有する画素のＹ座標値がいくつ連続して存在しているかをカウントする。連続する画素値“１”を有する画素数の最大値を第１画素数として記録する。
【００３６】
工程Ｓ１０４では、第２画素数取得部５３ｂによって第２画素値を得る。まず、画像データＤ１のアルミスロープ１を含む任意の範囲において、Ｙ方向におけるアルミスロープ１の像の一端１ａと他端１ｂとを繋ぐ他のエッジ１ｄを抽出する。エッジ１ｄは、上述したエッジ１ｃの反対側にあるエッジである。図６の（ａ）部は、アルミスロープ１の像における右側のエッジ１ｄが抽出された画像データＤ４である。このエッジ１ｄは、基準側面５に含まれる稜線５ａ又は平側面６に含まれる稜線６ａの何れか一方であるが、この時点ではどちらであるか不明である。
【００３７】
次に、画像データＤ４を所定の閾値で二値化して、図６の（ｂ）部に示される画像データＤ５を得る。画像データＤ５において、抽出されたエッジ１ｄでは、画素値“１”を有する画素同士がＹ方向に沿って互いに隣接している。このように互いに隣接している画素の数をカウントする。即ち、画像データＤ５において、Ｘ方向の所定幅の範囲に画素値“１”を有する画素のＹ座標値が、いくつ連続して存在しているかをカウントする。連続する画素値“１”を有する画素数の最大値を第２画素数として記録する。
【００３８】
以上の工程により、輝度値の傾き、第１画素数及び第２画素数が得られる。次に、工程Ｓ１０５では、配置形態判定部５４によって、輝度値の傾き、第１画素数及び第２画素数に基づいて、アルミスロープ１の配置形態の正誤を判定する。
【００３９】
ここで、アルミスロープ１の配置形態について、正しい配置形態と誤った配置形態の例を説明する。図７及び図８は、アルミスロープ１の正しい配置形態を示す図である。アルミスロープ１は、アルミスロープ１が設置される撮像台９の主面９ａと底面７とが接触するように配置された形態が正しい配置形態である。図７の（ａ）部は正しく配置されたアルミスロープ１の側面図であり、図７の（ｂ）部はその正面図である。図７の（ａ）部を参照すると、アルミスロープ１は、撮像台９の主面９ａと底面７とが接触するように配置されている。後端面４の−Ｙ方向側（上側）に先端面３が位置しているので、Ｙ方向に沿ってアルミスロープ１の厚さは増す。即ち、Ｙ方向に沿った輝度値の傾きは正である。図７の（ｂ）部を参照すると、平側面６の−Ｘ方向側（左側）に基準側面５が位置している。従って、基準側面５が左側にあり且つＹ方向に沿った輝度値の傾きが正の場合は、アルミスロープ１の配置形態は正しい。
【００４０】
正しい配置形態の別の例を説明する。図８の（ａ）部は正しく配置されたアルミスロープ１の側面図であり、図８の（ｂ）部はその正面図である。図８の（ａ）部を参照すると、アルミスロープ１は、撮像台９の主面９ａと底面７とが接触するように配置されている。後端面４の＋Ｙ方向側（下側）に先端面３が位置しているので、Ｙ方向に沿ってアルミスロープ１の厚さは減少する。即ち、Ｙ方向に沿った輝度値の傾きは負である。図８の（ａ）部を参照すると、平側面６の＋Ｘ方向側（右側）に基準側面５が位置している。このように、基準側面５が右側にあり且つＹ方向に沿った輝度値の傾きが負の場合は、アルミスロープ１の配置形態は正しい。
【００４１】
図９及び図１０は、アルミスロープ１の誤った配置形態を示す図である。図７及び図８と異なる点は、撮像台９の主面９ａと斜面２とが接触するように配置されている点である。
【００４２】
図９の（ａ）部は誤って配置されたアルミスロープ１の側面図であり、図９の（ｂ）部はその正面図である。図９の（ａ）部を参照すると、アルミスロープ１は、撮像台９の主面９ａと斜面２とが接触するように配置されている。また、後端面４の＋Ｙ方向側（下側）に先端面３が位置しているので、Ｙ方向に沿ってアルミスロープ１の厚さは減少する。即ち、Ｙ方向に沿った輝度値の傾きは負である。図９の（ｂ）部を参照すると、平側面６の−Ｘ方向側（左側）に基準側面５が位置している。このように、基準側面５が左側にあり、Ｙ方向に沿った輝度値の傾きが負の場合は、アルミスロープ１の配置形態は誤りである。
【００４３】
図１０の（ａ）部は誤って配置されたアルミスロープ１の側面図であり、図１０の（ｂ）部はその正面図である。図１０の（ａ）部を参照すると、アルミスロープ１は、撮像台９の主面９ａと斜面２とが接触するように配置されている。後端面４の−Ｙ方向側（上側）に先端面３が位置しているので、Ｙ方向に沿ってアルミスロープ１の厚さは増加する。即ち、Ｙ方向に沿った輝度値の傾きは正である。図１０の（ｂ）部を参照すると、基準側面５は＋Ｘ方向側（右側）に位置している。このように、基準側面５が右側にあり、Ｙ方向に沿った輝度値の傾きが正の場合は、アルミスロープ１の配置形態は誤りである。
【００４４】
図１１は、工程Ｓ１０５において、配置形態判定部５４によって行われる配置形態の正誤を判定するための主要なステップを説明するための図である。工程Ｓ１１０では、基準側面５が左側にあるか否かを判定する。この判定は、第１画素数と第２画素数に基づいて行われる。基準側面５は、基準溝５ｂを備えているので、Ｙ方向に連続する画素値“１”の画素数は平側面６よりも少ない。これ故、第１画素数と第２画素数とを比較して、値が小さい方が基準側面５に含まれる稜線５ａである。即ち、アルミスロープ１の像における左側のエッジ１ｃにおける連続する画素の最大値である第１画素数が、第２画素数よりも小さい場合に、工程Ｓ１１０の判定は正となる。
【００４５】
工程Ｓ１１０において正と判定された場合（Ｓ１１０：ＹＥＳ）は、続いて工程Ｓ１１１が実行される。工程Ｓ１１１では、輝度値の傾きが正であるか否かを判定する。傾きが正の場合（Ｓ１１１：ＹＥＳ）は、図７に示される設置形態であるので、工程Ｓ１１２によりアルミスロープ１の設置形態は正しいと判定され、工程Ｓ１０５の判定処理を終了する。傾きが負の場合（Ｓ１１１：ＮＯ）は、図９に示される設置形態であるので、Ｓ１１３によりアルミスロープ１の設置形態は誤りである判定され、工程Ｓ１０５の処理を終了する。
【００４６】
一方、工程Ｓ１１０において基準側面５は左側ではないと判定された場合（Ｓ１１０：ＮＯ）は、工程Ｓ１１４が実行される。工程Ｓ１１４では、輝度値の傾きが負であるか否かを判定する。傾きが負である場合（Ｓ１１４：ＹＥＳ）は、図８に示される設置形態であるので、Ｓ１１２によりアルミスロープ１の設置形態は正しいと判定され、工程Ｓ１０５の判定処理を終了する。傾きが正の場合（Ｓ１１４：ＮＯ）は、図１０に示される設置形態であるので、Ｓ１１５によりアルミスロープ１の設置形態は誤りである判定され、工程Ｓ１０５の処理を終了する。
【００４７】
工程Ｓ１０６では、輝度値補正部５５によってアルミスロープ１の配置形態が誤りであると判定された場合に、輝度値の分布を補正する。図１２は、輝度値の分布を補正する工程を説明するための図である。図１２の（ａ）部は、アルミスロープ１が誤って配置された場合における輝度値の分布の一例を示す図であり、さらに後端面４付近における輝度値の分布を拡大した様子を示す。図１２の（ｂ）部は、アルミスロープ１が正しく配置された場合における輝度値の分布の一例を示す図であり、さらに後端面４付近における輝度値の分布を拡大した様子を示す。
【００４８】
図１２の（ａ）部の拡大図と図１２の（ｂ）部の拡大図とを比較すると、範囲Ｂ１において、輝度値の分布Ｇ１、Ｇ２が異なっている。即ち、図１２の（ａ）部の拡大図では、後端面４付近の範囲Ｂ１で、輝度値の傾きが負になる。そのため、範囲Ｂ１では輝度値とアルミスロープ１の厚さを変換する変換係数を精度よく得ることができない。即ち、範囲Ｂ１はアルミスロープ１が誤って配置された場合に正しく配置された場合と輝度値の分布が異なる範囲である。そこで、工程Ｓ１０２において算出された近似直線を用いて、図１２の（ａ）部における直線Ｇ３に示されるように、範囲Ｂ１における輝度値を外挿して、元の輝度値の分布と置き換える。輝度値の分布Ｇ２に外挿された直線Ｇ３を加えて補正した輝度値の分布によって、アルミスロープ１が正しく配置されたときの輝度値の分布Ｇ１が推計される。尚、この範囲Ｂ１の幅は、輝度値の分布の全体の幅が、アルミスロープ１の長さＬ１と同じになるように設定される。このような処理を実行することで、範囲Ｂ１における輝度値の分布を補正する。
【００４９】
以上の各工程を実行することで、画像データＤ１の輝度値をアルミスロープの厚さに精度良く変換できる変換係数を生成可能な輝度値の分布が得られる。
【００５０】
上述した工程Ｓ１０１〜工程Ｓ１０５の処理と並行して、骨塩量測定処理が実行される。骨塩量測定処理では、測定処理部５６によって放射線画像検査装置５０に入力された画像データＤ１において所定の解析領域が設定される。その後、計測値計算部５７によって、設定された解析領域の画素値と、変換係数とを用いて骨密度パラメータが算出される。最後に、計測値計算部５７によって算出された骨密度パラメータが外部に出力される。
【００５１】
以下、コンピュータを放射線画像検査装置５０として動作させる放射線画像検査プログラムについて説明する。
【００５２】
本発明の実施形態に係る放射線画像検査プログラムは、記録媒体に格納されて提供される。記録媒体としては、フロッピー（登録商標）ディスク、ＣＤ−ＲＯＭ、ＤＶＤ、あるいはＲＯＭ等の記録媒体、あるいは半導体メモリ等が例示される。
【００５３】
図１３は、記録媒体に記録されたプログラムを実行するためのコンピュータのハードウェア構成を示す図であり、図１４は、記録媒体に記憶されたプログラムを実行するためのコンピュータの斜視図である。コンピュータとして、ＣＰＵを具備しソフトウエアによる処理や制御を行なうサーバ装置、パーソナルコンピュータ等の各種データ処理装置を含む。
【００５４】
図１３に示すように、コンピュータ３０は、フロッピー（登録商標）ディスクドライブ装置、ＣＤ−ＲＯＭドライブ装置、ＤＶＤドライブ装置等の読取装置１２と、オペレーティングシステムを常駐させた作業用メモリ（ＲＡＭ）１４と、記録媒体１０に記憶されたプログラムを記憶するメモリ１６と、ディスプレイといった表示装置１８と、入力装置であるマウス２０及びキーボード２２と、データ等の送受を行うための通信装置２４と、プログラムの実行を制御するＣＰＵ２６とを備えている。コンピュータ３０は、記録媒体１０が読取装置１２に挿入されると、読取装置１２から記録媒体１０に格納された放射線画像検査プログラムにアクセス可能になり、当該放射線画像検査プログラムによって、本実施形態の放射線画像検査装置５０として動作することが可能になる。
【００５５】
図１４に示すように、放射線画像検査プログラムは、搬送波に重畳されたコンピュータデータ信号４１としてネットワークを介して提供されるものであってもよい。この場合、コンピュータ３０は、通信装置２４によって受信した放射線画像検査プログラムをメモリ１６に格納し、当該放射線画像検査プログラムを実行することができる。
【００５６】
以上説明した放射線画像検査装置５０及びこれを用いた放射線画像検査方法によれば、画像データＤ１に基づいてＹ方向に沿ったアルミスロープ１の像の輝度値の分布を得る。これ故、該輝度値の分布の傾きを算出することができる。さらに、Ｙ方向におけるアルミスロープ１の像の一端１ａと他端１ｂとの間の２つのエッジ１ｃ、１ｄを抽出して、２つのエッジ１ｃ、１ｄのそれぞれにおいてＹ方向の連続画素数である第１及び第２画素数を得る。そして、輝度値の分布、第１画素数及び第２画素数に基づいて、アルミスロープ１の配置形態の正誤を判定する。この判定結果に基づいて輝度値の分布について補正の要否を判断できる。この判断結果に基づいて骨塩量等の骨に関する計測データを求めるための変換係数を得ることにより、計測データを精度良く算出できる変換係数を得ることができる。
【００５７】
また、輝度値補正部５５による輝度値補正処理によれば、正誤の判定によりアルミスロープ１の配置形態が誤と判定されたとき、輝度値の分布を補正する工程をさらに含んでいてもよい。このような構成によれば、誤った配置形態で撮影されたアルミスロープ１の像から得た輝度値の分布を補正する。この補正された輝度値の分布に基づいて、骨塩量等の骨に関する計測データを求めるための変換係数を得ることにより、計測データを精度良く算出できる変換係数を得ることができる。
【００５８】
なお、本発明は、前述した実施形態に限定されるものではなく、他に様々な変形が可能である。上述した実施形態では、輝度値の分布を得る工程（Ｓ１０１）と輝度値の傾きを得る工程（Ｓ１０２）とを先に実施し、その後に第１画素を得る工程（Ｓ１０３）と第２画素数を得る工程（Ｓ１０４）とを実施した方法が例示されてる。例えば、第１画素を得る工程（Ｓ１０３）と第２画素数を得る工程（Ｓ１０４）を先に実施し、輝度値の分布を得る工程（Ｓ１０１）と輝度値の傾きを得る工程（Ｓ１０２）とを後から実施してもよい。
【００５９】
また、工程Ｓ１０５によりアルミスロープ１の配置形態が誤りであると判定されたとき、工程Ｓ１０６によって輝度値の分布を補正することなく、アルミスロープ１の配置形態が誤りである旨を、放射線画像検査装置５０がアラーム等により報知することとしてもよい。
【符号の説明】
【００６０】
１…アルミスロープ、２…斜面、３…先端面、４…後端面、５…基準側面、５ｂ…基準溝、６…平側面、７…底面、Ｄ１〜Ｄ５…画像データ。
【図１】
2012231996.tif 000003
【図２】
2012231996.tif 000004
【図４】
2012231996.tif 000005
【図７】
2012231996.tif 000006
【図８】
2012231996.tif 000007
【図９】
2012231996.tif 000008
【図１０】
2012231996.tif 000009
【図１１】
2012231996.tif 000010
【図１２】
2012231996.tif 000011
【図１３】
2012231996.tif 000012
【図１４】
2012231996.tif 000013
【図３】
2012231996.tif 000014
【図５】
2012231996.tif 000015
【図６】
2012231996.tif 000016