(57)【要約】【課題】  適切な条件あるいは適切に近い条件で画像処理することが可能な画像処理方法および画像処理装置ならびに画像処理プログラムを実現する。【解決手段】  被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理装置であって、前記放射線画像の特徴量を算出する特徴量算出手段と、前記特徴量算出手段で算出された特徴量を特徴量評価関数により評価する特徴量評価手段と、前記特徴量評価手段で評価された前記特徴量の評価結果に基づいて画像処理のパラメータを決定するパラメータ決定手段と、前記パラメータ決定手段で決定された前記パラメータにより画像処理を実行する画像処理手段と、を有することを特徴とする。【選択図】      図１
2005109867.tif 000002
【特許請求の範囲】
【請求項１】
被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理方法であって、
前記放射線画像の特徴量を算出する特徴量算出ステップと、
前記特徴量算出ステップで算出された特徴量を特徴量評価関数により評価する特徴量評価ステップと、
前記特徴量評価ステップで評価された前記特徴量の評価結果に基づいて画像処理のパラメータを決定するパラメータ決定ステップと、
前記パラメータ決定ステップで決定された前記パラメータにより画像処理を実行する画像処理ステップと、
を有することを特徴とする画像処理方法。
【請求項２】
前記画像処理ステップで実行される画像処理は階調処理であり、
前記特徴量評価ステップは、画像処理としての階調処理における階調処理条件を参照して、前記特徴量評価関数により前記特徴量の評価を行う、
ことを特徴とする請求項１記載の画像処理方法。
【請求項３】
前記特徴量評価ステップで評価された前記特徴量の評価結果としての特徴量評価関数の出力が最大あるいは最小になる状態を用いて、前記パラメータ決定ステップで画像処理のパラメータを決定する、
ことを特徴とする請求項１または請求項２のいずれかに記載の画像処理方法。
【請求項４】
前記特徴量算出ステップで得られる放射線画像の特徴量は、画像全体あるいは画像の一部の統計量に基づいた特徴量であり、
前記特徴量評価関数は、１以上の変数によりパラメータライズされた変換方法により前記放射線画像を変換した際の特徴量の変化に基づいて決定される、
ことを特徴とする請求項３記載の画像処理方法。
【請求項５】
前記特徴量は、前記画像全体あるいは前記画像の一部の領域あるいは診断上重要となる関心領域における領域の最大値もしくは最小値のいずれかである、
ことを特徴とする請求項４記載の画像処理方法。
【請求項６】
特徴量は、領域のメジアンである、
ことを特徴とする請求項４記載の画像処理方法。
【請求項７】
特徴量は、領域の最頻値である、
ことを特徴とする請求項４記載の画像処理方法。
【請求項８】
特徴量は、領域の分散もしくは標準偏差のいずれか一方である、
ことを特徴とする請求項４記載の画像処理方法。
【請求項９】
統計量は画素値あるいは前記領域によって決定される請求項５から請求項７のいずれかに記載の統計量と、統計量から一定値を加算もしくは減算の少なくとも一方を実行した値との組から決定される量であって、
前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される、
ことを特徴とする請求項４記載の画像処理方法。
【請求項１０】
特徴量は、隣接あるいは所定間隔の画素の画素値差分であり、前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される特徴量が、前記特徴量変換後の差分値の変化率から決定される、
ことを特徴とする請求項４記載の画像処理方法。
【請求項１１】
前記統計量から一定値を加算もしくは減算の少なくとも一方を実行した値と前記特徴量との組から決定される特徴量は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つであり、
前記特徴量の変換後の値から決定される値は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つの変化率である、
ことを特徴とする請求項９記載の画像処理方法。
【請求項１２】
前記パラメータライズされた変換方法は、階調処理のためのルックアップテーブルによる画素値変換である、
ことを特徴とする請求項９乃至請求項１１のいずれかに記載の画像処理方法。
【請求項１３】
前記パラメータライズされた変換方法は、周波数強調処理である、
ことを特徴とする請求項９乃至請求項１１のいずれかに記載の画像処理方法。
【請求項１４】
前記パラメータライズされた変換方法は、イコライゼーション処理である、
ことを特徴とする請求項９乃至請求項１１のいずれかに記載の画像処理方法。
【請求項１５】
前記領域がエッジ検出処理により決定される、
ことを特徴とする請求項５乃至請求項９のいずれかに記載の画像処理方法。
【請求項１６】
前記領域が関心領域設定処理により決定される、
ことを特徴とする請求項５乃至請求項９のいずれか記載の画像処理方法。
【請求項１７】
前記エッジ検出処理は、ウェーブレット変換が使用される、
ことを特徴とする請求項１５記載の画像処理方法。
【請求項１８】
前記エッジ検出処理は、ガウス−ラプラシアンピラミッド変換が使用される、
ことを特徴とする請求項１５記載の画像処理方法。
【請求項１９】
被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理装置であって、
前記放射線画像の特徴量を算出する特徴量算出手段と、
前記特徴量算出手段で算出された特徴量を特徴量評価関数により評価する特徴量評価手段と、
前記特徴量評価手段で評価された前記特徴量の評価結果に基づいて画像処理のパラメータを決定するパラメータ決定手段と、
前記パラメータ決定手段で決定された前記パラメータにより画像処理を実行する画像処理手段と、
を有することを特徴とする画像処理装置。
【請求項２０】
前記画像処理手段で実行される画像処理は階調処理であり、
前記特徴量評価手段は、画像処理としての階調処理における階調処理条件を参照して、前記特徴量評価関数により前記特徴量の評価を行う、
ことを特徴とする請求項１９記載の画像処理装置。
【請求項２１】
前記特徴量評価手段で評価された前記特徴量の評価結果としての特徴量評価関数の出力が最大あるいは最小になる状態を用いて、前記パラメータ決定手段で画像処理のパラメータを決定する、
ことを特徴とする請求項１９または請求項２０のいずれかに記載の画像処理装置。
【請求項２２】
前記特徴量算出手段で得られる放射線画像の特徴量は、画像全体あるいは画像の一部の統計量に基づいた特徴量であり、
前記特徴量評価関数は、１以上の変数によりパラメータライズされた変換方法により前記放射線画像を変換した際の特徴量の変化に基づいて決定される、
ことを特徴とする請求項２１記載の画像処理装置。
【請求項２３】
特徴量は、前記特徴量は、前記画像全体あるいは前記画像の一部の領域あるいは診断上重要となる関心領域における領域の最大値もしくは最小値のいずれかである、
ことを特徴とする請求項２２記載の画像処理装置。
【請求項２４】
特徴量は、領域のメジアンである、
ことを特徴とする請求項２２記載の画像処理装置。
【請求項２５】
特徴量は、領域の最頻値である、
ことを特徴とする請求項２２記載の画像処理装置。
【請求項２６】
特徴量は、領域の分散もしくは標準偏差のいずれか一方である、
ことを特徴とする請求項２２記載の画像処理装置。
【請求項２７】
統計量は画素値あるいは前記領域によって決定される請求項５から請求項７のいずれかに記載の統計量と、統計量から一定値を加算もしくは減算の少なくとも一方を実行した値との組から決定される量であって、
前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される、
ことを特徴とする請求項２２記載の画像処理装置。
【請求項２８】
特徴量は、隣接あるいは所定間隔の画素の画素値差分であり、前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される特徴量が、前記特徴量変換後の差分値の変化率から決定される、
ことを特徴とする請求項２２記載の画像処理装置。
【請求項２９】
前記統計量から一定値を加算もしくは減算の少なくとも一方を実行した値と前記特徴量との組から決定される特徴量は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つであり、
前記特徴量の変換後の値から決定される値は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つの変化率である、
ことを特徴とする請求項２７記載の画像処理装置。
【請求項３０】
前記パラメータライズされた変換方法は、階調処理のためのルックアップテーブルによる画素値変換である、
ことを特徴とする請求項２７乃至請求項２９のいずれかに記載の画像処理装置。
【請求項３１】
前記パラメータライズされた変換方法は、周波数強調処理である、
ことを特徴とする請求項２７乃至請求項２９のいずれかに記載の画像処理装置。
【請求項３２】
前記パラメータライズされた変換方法は、イコライゼーション処理である、
ことを特徴とする請求項２７乃至請求項２９のいずれかに記載の画像処理装置。
【請求項３３】
前記領域がエッジ検出処理により決定される、
ことを特徴とする請求項２３乃至請求項２７記載の画像処理装置。
【請求項３４】
前記領域が関心領域設定処理により決定される、
ことを特徴とする請求項２３乃至請求項２７記載の画像処理装置。
【請求項３５】
前記エッジ検出処理は、ウェーブレット変換が使用される、
ことを特徴とする請求項３３記載の画像処理装置。
【請求項３６】
前記エッジ検出処理は、ガウス−ラプラシアンピラミッド変換が使用される、
ことを特徴とする請求項３３記載の画像処理装置。
【請求項３７】
被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理プログラムであって、
前記放射線画像の特徴量を算出する特徴量算出ルーチンと、
前記特徴量算出ルーチンで算出された特徴量を特徴量評価関数により評価する特徴量評価ルーチンと、
前記特徴量評価ルーチンで評価された前記特徴量の評価結果に基づいて画像処理のパラメータを決定するパラメータ決定ルーチンと、
前記パラメータ決定ルーチンで決定された前記パラメータにより画像処理を実行する画像処理ルーチンと、
を有することを特徴とする画像処理プログラム。
【請求項３８】
前記画像処理ルーチンで実行される画像処理は階調処理であり、
前記特徴量評価ルーチンは、画像処理としての階調処理における階調処理条件を参照して、前記特徴量評価関数により前記特徴量の評価を行う、
ことを特徴とする請求項３７記載の画像処理プログラム。
【請求項３９】
前記特徴量評価ルーチンで評価された前記特徴量の評価結果としての特徴量評価関数の出力が最大あるいは最小になる状態を用いて、前記パラメータ決定ルーチンで画像処理のパラメータを決定する、
ことを特徴とする請求項３７または請求項３８のいずれかに記載の画像処理プログラム。
【請求項４０】
前記特徴量算出ルーチンで得られる放射線画像の特徴量は、画像全体あるいは画像の一部の統計量に基づいた特徴量であり、
前記特徴量評価関数は、１以上の変数によりパラメータライズされた変換方法により前記放射線画像を変換した際の特徴量の変化に基づいて決定される、
ことを特徴とする請求項３９記載の画像処理プログラム。
【請求項４１】
特徴量は、前記特徴量は、前記画像全体あるいは前記画像の一部の領域あるいは診断上重要となる関心領域における領域の最大値もしくは最小値のいずれかである、
ことを特徴とする請求項４０記載の画像処理プログラム。
【請求項４２】
特徴量は、領域のメジアンである、
ことを特徴とする請求項４０記載の画像処理プログラム。
【請求項４３】
特徴量は、領域の最頻値である、
ことを特徴とする請求項４０記載の画像処理プログラム。
【請求項４４】
特徴量は、領域の分散もしくは標準偏差のいずれか一方である、
ことを特徴とする請求項４０記載の画像処理プログラム。
【請求項４５】
統計量は画素値あるいは前記領域によって決定される請求項５から請求項７のいずれかに記載の統計量と、統計量から一定値を加算もしくは減算の少なくとも一方を実行した値との組から決定される量であって、
前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される、
ことを特徴とする請求項４０記載の画像処理プログラム。
【請求項４６】
特徴量は、隣接あるいは所定間隔の画素の画素値差分であり、前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される特徴量が、前記特徴量変換後の差分値の変化率から決定される、
ことを特徴とする請求項４０記載の画像処理プログラム。
【請求項４７】
前記統計量から一定値を加算もしくは減算の少なくとも一方を実行した値と前記特徴量との組から決定される特徴量は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つであり、
前記特徴量の変換後の値から決定される値は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つの変化率である、
ことを特徴とする請求項４５記載の画像処理プログラム。
【請求項４８】
前記パラメータライズされた変換方法は、階調処理のためのルックアップテーブルによる画素値変換である、
ことを特徴とする請求項４５乃至請求項４７のいずれかに記載の画像処理プログラム。
【請求項４９】
前記パラメータライズされた変換方法は、周波数強調処理である、
ことを特徴とする請求項４５乃至請求項４７のいずれかに記載の画像処理プログラム。
【請求項５０】
前記パラメータライズされた変換方法は、イコライゼーション処理である、
ことを特徴とする請求項４５乃至請求項４７のいずれかに記載の画像処理プログラム。
【請求項５１】
前記領域がエッジ検出処理により決定される、
ことを特徴とする請求項４１乃至請求項４５記載の画像処理プログラム。
【請求項５２】
前記領域が関心領域設定処理により決定される、
ことを特徴とする請求項４１乃至請求項４５記載の画像処理プログラム。
【請求項５３】
前記エッジ検出処理は、ウェーブレット変換が使用される、
ことを特徴とする請求項５１記載の画像処理プログラム。
【請求項５４】
前記エッジ検出処理は、ガウス−ラプラシアンピラミッド変換が使用される、
ことを特徴とする請求項５１記載の画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は放射線画像を処理する際の画像処理方法および画像処理装置ならびに画像処理プログラムに関し、さらに詳しくは、診断等に適した放射線画像を得ることができる画像処理方法および画像処理装置ならびに画像処理プログラムに関する。
【背景技術】
【０００２】
近年、放射線画像を直接デジタル画像として撮影できる装置が開発されている。たとえば、被写体に照射された放射線量を検出し、その検出量に対応して形成される放射線画像を電気信号として得る装置としては、輝尽性蛍光体を用いたディテクタを用いる方法が特開昭55-12429号公報、特開昭63-189853号公報など、多数開示されている。
【０００３】
このような装置では、シート状の基板に輝尽性蛍光体を塗布、あるいは蒸着等によって固着したディテクタに、いったん被写体を透過した放射線を照射して輝尽性蛍光体に放射線を吸収させる。
【０００４】
その後、この輝尽性蛍光体を光または熱エネルギ−で励起することにより、この輝尽性蛍光体が上記吸収によって蓄積している放射線エネルギ−を蛍光として放射させ、この蛍光を光電変換して画像信号を得るようにしている。
【０００５】
一方、照射された放射線の強度に応じた電荷を光導電層に生成し、生成された電荷を二次元的に配列された複数のコンデンサに蓄積し、それら蓄積された電荷を取り出すことにより得られる放射線画像検出装置が提案されている。
【０００６】
このような放射線画像検出装置では、フラットパネルディテクタ（ＦＰＤ）と呼ばれるものを使用している。この種のＦＰＤは、特開平9-90048号公報に記載されているように、照射された放射線強度に応じた蛍光を発する蛍光体と、蛍光体から発する蛍光を直接または縮小光学系を介して受光して光電変換を行うフォトダイオードやＣＣＤのような光電変換素子の組み合わせによって実現されるものが知られている。
【０００７】
また特開平6-342098号公報に記載されているように、照射された放射を直接電荷に変換するものも知られている。
【０００８】
これらの放射線画像検出装置では、取得した画像に対して、診断に適した画像となるよう、階調変換処理やエッジ強調処理などの画像処理を行うのが一般的である。
【０００９】
なお、このようにして得られた画像データに基づく放射線画像を表示したり出力する場合、撮影条件の変動を受けることなく見やすい放射線画像となるように画像処理が行われている。
【００１０】
このため、放射線画像の所望の領域内の画像データの累積ヒストグラムを求めたのち、この累積ヒストグラムが所定の割合となるデータレベルを基準信号値と設定して画像処理を行うことが、以下の特許文献１に記載されている。
【００１１】
また、高信号値領域と低信号値領域の分布状態に基づいて、画像処理条件を定めて、適切な画像処理を行うことが、以下の特許文献２に記載されている。
【特許文献１】特開平６−６１３２５号公報（第１頁、図１）
【特許文献２】特開２０００−１５７５１８７号公報（第１〜５頁、図４）
【発明の開示】
【発明が解決しようとする課題】
【００１２】
ところで、放射線画像では、被写体の撮影部位によって高濃度領域（放射線透過量の大きい領域）と低濃度領域（放射線透過量の小さい領域）の割合が大きく異なる。また、被写体の状態、例えば患者の呼吸状態によって肺野内の濃度が著しく異なるものとされる。      
【００１３】
このため、放射線画像の所望の領域内の画像データの累積ヒストグラムを求めたのち、この累積ヒストグラムが所定の割合となるデータレベルを基準信号値と設定して画像処理を行う場合（特開平６−６１３２５号公報）には、高濃度領域と低濃度領域の割合の違いによっては、診断等に適した見やすい放射線画像を得ることができない恐れがある。
【００１４】
例えば、基準信号値が所定の濃度となるように階調処理を行った場合、低濃度領域が支配的であるときには基準信号値が低い値となるので全体的に高濃度な画像となり、高濃度領域が支配的であるときには基準信号値が高い値となるので全体的に低濃度な画像となる。
【００１５】
しかし、例えば肋骨のような部位では上部肋骨部（高濃度領域の多い肺野部分が含まれる）や下部肋骨部（高濃度領域の少ない腹部部分が含まれる）を撮影して診断を行う場合でも、診断で必要とされる部分は高濃度領域と低濃度領域の両方にまたがっており、画像全体の濃度が低濃度あるいは高濃度に偏ることは、診断を行う上で好ましくない。
【００１６】
すなわち、従来は、画像処理のパラメータを決定するために様々な手法が提案されているが、その手法は画像内の特定領域（関心領域など）の信号値やヒストグラムから決定される信号値あるいはパラメータとして一意的に与えるなどのように決定され、得られる領域やパラメータが不適切な場合、不適切な画像処理が得られる場合があった。
【００１７】
本発明は以上のような課題に鑑みてなされたものであって、画像処理条件決定の失敗により診断が困難となる状況の発生を防止し、いかなる場合であっても適切な条件あるいは適切に近い条件で画像処理することが可能な画像処理方法および画像処理装置ならびに画像処理プログラムを実現することを目的とする。
【課題を解決するための手段】
【００１８】
すなわち、前記した課題を解決は、以下に列記する発明により解決される。
【００１９】
（１）請求項１記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理方法であって、前記放射線画像の特徴量を算出する特徴量算出ステップと、前記特徴量算出ステップで算出された特徴量を特徴量評価関数により評価する特徴量評価ステップと、前記特徴量評価ステップで評価された前記特徴量の評価結果に基づいて画像処理のパラメータを決定するパラメータ決定ステップと、前記パラメータ決定ステップで決定された前記パラメータにより画像処理を実行する画像処理ステップと、を有することを特徴とする画像処理方法である。
【００２０】
（２）請求項２記載の発明は、前記画像処理ステップで実行される画像処理は階調処理であり、前記特徴量評価ステップは、画像処理としての階調処理における階調処理条件を参照して、前記特徴量評価関数により前記特徴量の評価を行う、ことを特徴とする請求項１記載の画像処理方法である。
【００２１】
（３）請求項３記載の発明は、前記特徴量評価ステップで評価された前記特徴量の評価結果としての特徴量評価関数の出力が最大あるいは最小になる状態を用いて、前記パラメータ決定ステップで画像処理のパラメータを決定する、ことを特徴とする請求項１または請求項２のいずれかに記載の画像処理方法である。
【００２２】
（４）請求項４記載の発明は、前記特徴量算出ステップで得られる放射線画像の特徴量は、画像全体あるいは画像の一部の統計量に基づいた特徴量であり、前記特徴量評価関数は、１以上の変数によりパラメータライズされた変換方法により前記放射線画像を変換した際の特徴量の変化に基づいて決定される、ことを特徴とする請求項３記載の画像処理方法である。
【００２３】
（５）請求項５記載の発明は、前記特徴量は、前記画像全体あるいは前記画像の一部の領域あるいは診断上重要となる関心領域における領域の最大値もしくは最小値のいずれかである、ことを特徴とする請求項４記載の画像処理方法である。
【００２４】
（６）請求項６記載の発明は、特徴量は、領域のメジアンである、ことを特徴とする請求項４記載の画像処理方法である。
【００２５】
（７）請求項７記載の発明は、特徴量は、領域の最頻値である、ことを特徴とする請求項４記載の画像処理方法である。
【００２６】
（８）請求項８記載の発明は、特徴量は、領域の分散もしくは標準偏差のいずれか一方である、ことを特徴とする請求項４記載の画像処理方法である。
【００２７】
（９）請求項９記載の発明は、統計量は画素値あるいは前記領域によって決定される請求項５から請求項７のいずれかに記載の統計量と、統計量から一定値を加算もしくは減算の少なくとも一方を実行した値との組から決定される量であって、前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される、ことを特徴とする請求項４記載の画像処理方法である。
【００２８】
（１０）請求項１０記載の発明は、特徴量は、隣接あるいは所定間隔の画素の画素値差分であり、前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される特徴量が、前記特徴量変換後の差分値の変化率から決定される、ことを特徴とする請求項４記載の画像処理方法である。
【００２９】
（１１）請求項１１記載の発明は、前記統計量から一定値を加算もしくは減算の少なくとも一方を実行した値と前記特徴量との組から決定される特徴量は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つであり、前記特徴量の変換後の値から決定される値は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つの変化率である、ことを特徴とする請求項９記載の画像処理方法である。
【００３０】
（１２）請求項１２記載の発明は、前記パラメータライズされた変換方法は、階調処理のためのルックアップテーブルによる画素値変換である、ことを特徴とする請求項９乃至請求項１１のいずれかに記載の画像処理方法である。
【００３１】
（１３）請求項１３記載の発明は、前記パラメータライズされた変換方法は、周波数強調処理である、ことを特徴とする請求項９乃至請求項１１のいずれかに記載の画像処理方法である。
【００３２】
（１４）請求項１４記載の発明は、前記パラメータライズされた変換方法は、イコライゼーション処理である、ことを特徴とする請求項９乃至請求項１１のいずれかに記載の画像処理方法である。
【００３３】
（１５）請求項１５記載の発明は、前記領域がエッジ検出処理により決定される、ことを特徴とする請求項５乃至請求項９のいずれかに記載の画像処理方法である。
【００３４】
（１６）請求項１６記載の発明は、前記領域が関心領域設定処理により決定される、ことを特徴とする請求項５乃至請求項９のいずれかに記載の画像処理方法である。
【００３５】
（１７）請求項１７記載の発明は、前記エッジ検出処理は、ウェーブレット変換が使用される、ことを特徴とする請求項１５記載の画像処理方法である。
【００３６】
（１８）請求項１８記載の発明は、前記エッジ検出処理は、ガウス−ラプラシアンピラミッド変換が使用される、ことを特徴とする請求項１５記載の画像処理方法である。
【００３７】
（１９）請求項１９記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理装置であって、前記放射線画像の特徴量を算出する特徴量算出手段と、前記特徴量算出手段で算出された特徴量を特徴量評価関数により評価する特徴量評価手段と、前記特徴量評価手段で評価された前記特徴量の評価結果に基づいて画像処理のパラメータを決定するパラメータ決定手段と、前記パラメータ決定手段で決定された前記パラメータにより画像処理を実行する画像処理手段と、を有することを特徴とする画像処理装置である。
【００３８】
（２０）請求項２０記載の発明は、前記画像処理手段で実行される画像処理は階調処理であり、前記特徴量評価手段は、画像処理としての階調処理における階調処理条件を参照して、前記特徴量評価関数により前記特徴量の評価を行う、ことを特徴とする請求項１９記載の画像処理装置である。
【００３９】
（２１）請求項２１記載の発明は、前記特徴量評価手段で評価された前記特徴量の評価結果としての特徴量評価関数の出力が最大あるいは最小になる状態を用いて、前記パラメータ決定手段で画像処理のパラメータを決定する、ことを特徴とする請求項１９または請求項２０のいずれかに記載の画像処理装置である。
【００４０】
（２２）請求項２２記載の発明は、前記特徴量算出手段で得られる放射線画像の特徴量は、画像全体あるいは画像の一部の統計量に基づいた特徴量であり、前記特徴量評価関数は、１以上の変数によりパラメータライズされた変換方法により前記放射線画像を変換した際の特徴量の変化に基づいて決定される、ことを特徴とする請求項２１記載の画像処理装置である。
【００４１】
（２３）請求項２３記載の発明は、前記特徴量は、前記画像全体あるいは前記画像の一部の領域あるいは診断上重要となる関心領域における領域の最大値もしくは最小値のいずれかである、ことを特徴とする請求項２２記載の画像処理装置である。
【００４２】
（２４）請求項２４記載の発明は、特徴量は、領域のメジアンである、ことを特徴とする請求項２２記載の画像処理装置である。
【００４３】
（２５）請求項２５記載の発明は、特徴量は、領域の最頻値である、ことを特徴とする請求項２２記載の画像処理装置である。
【００４４】
（２６）請求項２６記載の発明は、特徴量は、領域の分散もしくは標準偏差のいずれか一方である、ことを特徴とする請求項２２記載の画像処理装置である。
【００４５】
（２７）請求項２７記載の発明は、統計量は画素値あるいは前記領域によって決定される請求項５から請求項７のいずれかに記載の統計量と、統計量から一定値を加算もしくは減算の少なくとも一方を実行した値との組から決定される量であって、前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される、ことを特徴とする請求項２２記載の画像処理装置である。
【００４６】
（２８）請求項２８記載の発明は、特徴量は、隣接あるいは所定間隔の画素の画素値差分であり、前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される特徴量が、前記特徴量変換後の差分値の変化率から決定される、ことを特徴とする請求項２２記載の画像処理装置である。
【００４７】
（２９）請求項２９記載の発明は、前記統計量から一定値を加算もしくは減算の少なくとも一方を実行した値と前記特徴量との組から決定される特徴量は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つであり、前記特徴量の変換後の値から決定される値は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つの変化率である、ことを特徴とする請求項２７記載の画像処理装置である。
【００４８】
（３０）請求項３０記載の発明は、前記パラメータライズされた変換方法は、階調処理のためのルックアップテーブルによる画素値変換である、ことを特徴とする請求項２７乃至請求項２９のいずれかに記載の画像処理装置である。
【００４９】
（３１）請求項３１記載の発明は、前記パラメータライズされた変換方法は、周波数強調処理である、ことを特徴とする請求項２７乃至請求項２９のいずれかに記載の画像処理装置である。
【００５０】
（３２）請求項３２記載の発明は、前記パラメータライズされた変換方法は、イコライゼーション処理である、ことを特徴とする請求項２７乃至請求項２９のいずれかに記載の画像処理装置である。
【００５１】
（３３）請求項３３記載の発明は、前記領域がエッジ検出処理により決定される、ことを特徴とする請求項２３乃至請求項２７記載の画像処理装置である。
【００５２】
（３４）請求項３４記載の発明は、前記領域が関心領域設定処理により決定される、ことを特徴とする請求項２３乃至請求項２７記載の画像処理装置である。
【００５３】
（３５）請求項３５記載の発明は、前記エッジ検出処理は、ウェーブレット変換が使用される、ことを特徴とする請求項３３記載の画像処理装置である。
【００５４】
（３６）請求項３６記載の発明は、前記エッジ検出処理は、ガウス−ラプラシアンピラミッド変換が使用される、ことを特徴とする請求項３３記載の画像処理装置である。
【００５５】
（３７）請求項３７記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理プログラムであって、前記放射線画像の特徴量を算出する特徴量算出ルーチンと、前記特徴量算出ルーチンで算出された特徴量を特徴量評価関数により評価する特徴量評価ルーチンと、前記特徴量評価ルーチンで評価された前記特徴量の評価結果に基づいて画像処理のパラメータを決定するパラメータ決定ルーチンと、前記パラメータ決定ルーチンで決定された前記パラメータにより画像処理を実行する画像処理ルーチンと、を有することを特徴とする画像処理プログラムである。
【００５６】
（３８）請求項３８記載の発明は、前記画像処理ルーチンで実行される画像処理は階調処理であり、前記特徴量評価ルーチンは、画像処理としての階調処理における階調処理条件を参照して、前記特徴量評価関数により前記特徴量の評価を行う、ことを特徴とする請求項３７記載の画像処理プログラムである。
【００５７】
（３９）請求項３９記載の発明は、前記特徴量評価ルーチンで評価された前記特徴量の評価結果としての特徴量評価関数の出力が最大あるいは最小になる状態を用いて、前記パラメータ決定ルーチンで画像処理のパラメータを決定する、ことを特徴とする請求項３７または請求項３８のいずれかに記載の画像処理プログラムである。
【００５８】
（４０）請求項４０記載の発明は、前記特徴量算出ルーチンで得られる放射線画像の特徴量は、画像全体あるいは画像の一部の統計量に基づいた特徴量であり、前記特徴量評価関数は、１以上の変数によりパラメータライズされた変換方法により前記放射線画像を変換した際の特徴量の変化に基づいて決定される、ことを特徴とする請求項３９記載の画像処理プログラムである。
【００５９】
（４１）請求項４１記載の発明は、特徴量は、前記特徴量は、前記画像全体あるいは前記画像の一部の領域あるいは診断上重要となる関心領域における領域の最大値もしくは最小値のいずれかである、ことを特徴とする請求項４０記載の画像処理プログラムである。
【００６０】
（４２）請求項４２記載の発明は、特徴量は、領域のメジアンである、ことを特徴とする請求項４０記載の画像処理プログラムである。
【００６１】
（４３）請求項４３記載の発明は、特徴量は、領域の最頻値である、ことを特徴とする請求項４０記載の画像処理プログラムである。
【００６２】
（４４）請求項４４記載の発明は、特徴量は、領域の分散もしくは標準偏差のいずれか一方である、ことを特徴とする請求項４０記載の画像処理プログラムである。
【００６３】
（４５）請求項４５記載の発明は、統計量は画素値あるいは前記領域によって決定される請求項５から請求項７のいずれかに記載の統計量と、統計量から一定値を加算もしくは減算の少なくとも一方を実行した値との組から決定される量であって、前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される、ことを特徴とする請求項４０記載の画像処理プログラムである。
【００６４】
（４６）請求項４６記載の発明は、特徴量は、隣接あるいは所定間隔の画素の画素値差分であり、前記放射線画像を前記パラメータライズされた変換方法により変換された際の特徴量の変化に基づいて決定される特徴量が、前記特徴量変換後の差分値の変化率から決定される、ことを特徴とする請求項４０記載の画像処理プログラムである。
【００６５】
（４７）請求項４７記載の発明は、前記統計量から一定値を加算もしくは減算の少なくとも一方を実行した値と前記特徴量との組から決定される特徴量は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つであり、前記特徴量の変換後の値から決定される値は、それらの差分の平均、最大、最小、分散、標準偏差の少なくとも一つの変化率である、ことを特徴とする請求項４５記載の画像処理プログラムである。
【００６６】
（４８）請求項４８記載の発明は、前記パラメータライズされた変換方法は、階調処理のためのルックアップテーブルによる画素値変換である、ことを特徴とする請求項４５乃至請求項４７のいずれかに記載の画像処理プログラムである。
【００６７】
（４９）請求項４９記載の発明は、前記パラメータライズされた変換方法は、周波数強調処理である、ことを特徴とする請求項４５乃至請求項４７のいずれかに記載の画像処理プログラムである。
【００６８】
（５０）請求項５０記載の発明は、前記パラメータライズされた変換方法は、イコライゼーション処理である、ことを特徴とする請求項４５乃至請求項４７のいずれかに記載の画像処理プログラムである。
【００６９】
（５１）請求項５１記載の発明は、前記領域がエッジ検出処理により決定される、ことを特徴とする請求項４１乃至請求項４５記載の画像処理プログラムである。
【００７０】
（５２）請求項５２記載の発明は、前記領域が関心領域設定処理により決定される、ことを特徴とする請求項４１乃至請求項４５記載の画像処理プログラムである。
【００７１】
（５３）請求項５３記載の発明は、前記エッジ検出処理は、ウェーブレット変換が使用される、ことを特徴とする請求項５１記載の画像処理プログラムである。
【００７２】
（５４）請求項５４記載の発明は、前記エッジ検出処理は、ガウス−ラプラシアンピラミッド変換が使用される、ことを特徴とする請求項５１記載の画像処理プログラムである。
【発明の効果】
【００７３】
本発明では、放射線画像の特徴量を算出し、算出された特徴量を特徴量評価関数により評価し、評価された前記特徴量の評価結果に基づいて画像処理のパラメータを決定し、決定された前記パラメータにより画像処理を実行する。
【００７４】
したがって、放射線画像毎の特徴量に基づいて評価して決定したパラメータにより画像処理を実行するので、例えば、各信号値に対して特徴量に基づいて、重要度に応じた重み付けをして、各画素の階調処理用ＬＵＴにより変換した際の信号の増幅率を算出し、増幅率と重みを掛け合わせた値のすべての画素での総和が最大になるようなＬＵＴを決定することで適切な画像処理が実行できる。ここで特徴量に基づいた重要度とは、たとえば一般に診断上重要なポイントは画像中央に位置することが多く、さらにエッジ成分をもつことなどから、特徴量として画像中心との距離を測る画像中央度と微分フィルタなどを用いて検出した各画素におけるエッジ強度を使用し、それらの値に基づいて決定される量である。また、増幅率としてはたとえば各画素値におけるＬＵＴの傾きやそれに相当する量を考えることができる。
【００７５】
この結果、画像処理条件決定の失敗により診断が困難となる状況の発生を防止し、いかなる場合であっても適切な条件あるいは適切に近い条件で画像処理することが可能になる。
【発明を実施するための最良の形態】
【００７６】
以下、図面を参照して本発明の実施をするための最良の形態を詳細に説明する。
【００７７】
本発明の実施をするための最良の形態の画像処理方法および画像処理装置ならびに画像処理プログラムの好適な実施の形態について説明する。なお、これにより本発明が限定されるものではない。
【００７８】
なお、本実施例の形態例の各手段は、ハードウェアやファームウェア、またはソフトウェアで構成することが可能である。このため、画像処理方法の各ステップ，画像処理装置の各手段，画像処理プログラムの各ルーチンの処理手順に沿った機能ブロック図としての図１を示す。
【００７９】
以下、本実施の形態例の構成および動作について、図１のブロック図、図２のフローチャート、および、その他の説明図を参照しつつ詳細に説明する。なお、図１における各手段は、画像形成装置の各手段を示すだけでなく、画像処理方法の各ステップ、画像処理プログラムの各ルーチンを示すものである。
【００８０】
〈全体構成および処理の流れ〉
1)全体構成：
放射線発生装置３０、放射線画像読取器４０、画像処理装置１００が、図１に示されるように構成されている。
【００８１】
また、画像処理装置１００内部には、制御手段１０１、画像データ生成手段１１０、特徴量算出手段１２０、特徴量統合手段１３０、特徴量評価手段１４０、パラメータ決定手段１５０が、図１に示すように構成される。
【００８２】
2)処理の流れ：
・制御手段１０１は放射線画像撮影・読み取り、および、画像処理のパラメータ決定に伴う各種処理についての制御を行う。
・放射線発生装置３０からの放射線は被写体５を透過し、この被写体５を透過した放射線は放射線画像読取器４０により読み取られる。
・放射線画像読取器４０により読み取られた信号は、画像データ生成手段１１０で画像データに変換される（図２Ｓ１）。
・特徴量算出手段１２０は、放射線画像データの特徴量を算出する（図２Ｓ２）。
・特徴量統合手段１３０は、放射線画像データ特徴量を統合する（図２Ｓ３）。
・特徴量評価手段１４０は、算出あるいは統合された特徴量を特徴量評価関数により評価する（図２Ｓ４）。
・パラメータ決定手段１５０は、評価された前記特徴量の評価結果に基づいて、画像処理のパラメータ（画像処理条件）を決定する（図２Ｓ５）。
・画像処理手段１６０は、画像データ生成手段１１０からの画像データについて、パラメータ決定手段１５０で決定されたパラメータに従って画像処理を実行する（図２Ｓ６）。
【００８３】
〈各手段、各ステップの詳細内容〉
（１）各部の動作・制御：
制御手段１０１では、まず、撮影部位または撮影方向等の情報をユーザインターフェースなどから取得する。これらの情報は、ユーザが撮影部位等を特定することで行う。例えば、表示部とタッチパネルとを兼ね備えた当該画像処理装置のユーザインタフェース（図示せず）から、撮影部位が表示されたボタンを押すこと等により入力される。そのほか、磁気カード、バーコード、ＨＩＳ（病院内情報システム：ネットワークによる情報管理）等を利用しても行われる。
【００８４】
（２）放射線画像入力：
放射線発生器３０は制御手段１０１によって制御されており、放射線発生器３０から放射された放射線は、被写体５を通して放射線画像読取器４０の前面に装着されている撮像パネルに照射される。そして、放射線画像読取器４０では、被写体５を透過した放射線を検知して、画像信号として取得する。
【００８５】
具体的な構成例としては、輝尽性蛍光体プレートを用いたものとして、特開平１１-１４２９９８号公報や特開２００２-１５６７１６号公報に記載されたものがある。また、フラットパネルディテクタ（FPD）を入力装置として用いるものには、特開平６-３４２０９８号公報に記載された、検出したX線を直接電荷に変換し、画像信号として取得するものや、特開平９-９００４８号公報に記載された、検出したＸ線を一旦光に変換した後、その光を受光して電荷に変換する、間接方式のものがある。
【００８６】
なお、放射線画像読取器４０は、放射線画像を記録した銀塩フィルムにレーザ、蛍光灯等の光源からの光を照射し、この銀塩フィルムの透過光を光電変換して画像データを生成してもよい。また、放射線量子計数型検出器を用いて放射線エネルギーを直接電気信号に変換して画像データを生成する構成であってもよい。
【００８７】
被写体５の放射線画像を得る際には、放射線発生器３０と放射線画像読取器４０の撮像パネルの間に被写体５が位置するものとされて、放射線発生器３０から放射された放射線が被写体５に照射されると共に、被写体５を透過した放射線が撮像パネルに入射される。
【００８８】
（３）関心領域設定：
ところで、放射線画像の撮影に際しては、例えば診断に必要とされない部分に放射線が照射されないようにするため、あるいは診断に必要とされない部分に放射線が照射されて、この部分で散乱された放射線が診断に必要とされる部分に入射されて分解能が低下することを防止するため、被写体５の一部や放射線発生器３０に鉛板等の放射線非透過物質を設置して、被写体５に対する放射線の照射野を制限する照射野絞りが行われる。
【００８９】
この照射野絞りが行われた場合、照射野内領域と照射野外領域の画像データを用いてレベルの変換処理やその後の階調処理を行うものとすると、照射野外領域の画像データによって、照射野内の診断に必要とされる部分の画像処理が適正に行われなくなってしまう。このため、画像データ生成手段１１０では、照射野内領域と照射野外領域を判別する照射野認識が行われる。
【００９０】
照射野認識では、例えば特開昭６３−２５９５３８号で示される方法が用いられて、図３Ａに示すように撮像面上の所定の位置Ｐから撮像面の端部側に向かう線分上の画像データを用いて例えば微分処理が行われる。この微分処理によって得られた微分信号Ｓdは、図３Ｂに示すように照射野エッジ部で信号レベルが大きくなるため、微分信号Ｓdの信号レベルを判別して１つの照射野エッジ候補点ＥＰ1が求められる。この照射野エッジ候補点を求める処理を、撮像面上の所定の位置を中心として放射状に行うことにより複数の照射野エッジ候補点ＥＰ1〜ＥＰkが求められる。このようにして得られた複数の照射野エッジ候補点ＥＰ1〜ＥＰkの隣接するエッジ候補点を直線あるいは曲線で結ぶことにより照射野エッジ部が求められる。
【００９１】
また、特開平５−７５７９号で示される方法を用いることもできる。この方法では、撮像面を複数の小領域に分割したとき、照射野絞りによって放射線の照射が遮られた照射野外の小領域では、略一様に放射線の放射線量が小さくなり画像データの分散値が小さくなる。また、照射野内の小領域では、被写体によって放射線量が変調されることから照射野外に比べて分散値が高くなる。さらに、照射野エッジ部を含む小領域では最も放射線量が小さい部分と被写体によって変調された放射線量の部分が混在することから分散値は最も高くなる。このことから、分散値によって照射野エッジ部を含む小領域が判別される。
【００９２】
また、特開平７−１８１６０９号で示される方法を用いることもできる。この方法では、画像データを所定の回転中心に関して回転移動させて、平行状態検出手段によって照射野の境界線が画像上に設定された直交座標の座標軸と平行となるまで回転を行うものとし、平行状態が検出されると、直線方程式算出手段によって回転角度と回転中心から境界線までの距離によって回転前の境界の直線方程式が算出される。その後、複数の境界線に囲まれる領域を直線方程式から決定することで、照射野の領域を判別することができる。また照射野エッジ部が曲線である場合には、境界点抽出手段で画像データに基づき例えば１つの境界点を抽出し、この境界点の周辺の境界候補点群から次の境界点を抽出する。以下同様に、境界点の周辺の境界候補点群から境界点を順次抽出することにより、照射野エッジ部が曲線であっても判別することができる。
【００９３】
照射野認識が行われると、放射線画像読取器からの画像データＤＴの分布を所望のレベルの分布に変換する際に、放射線画像読取器からの画像データＤＴのレベルの分布を決定するための領域（以下「関心領域」という）の設定を行う。
【００９４】
この関心領域の設定では、たとえば胸部正面画像では診断上重要な領域がすべて含まれるように肺野全体に設定される。この関心領域の設定は例えば図７のように画像内の1/3〜2/3の領域での水平方向の平均プロファイルを作成し、はじめに画像中央の最小値Pを決定する。次に画像の両端に向かって、最大値を検出する（図７中のA、B）。次にそれぞれの最大値から画像端に向かって、最大値と最小値間の所定の割合になる点を探すことで関心領域の左右端を決定することができる。この所定の割合については例えば経験的に4/5程度にすることで、適切に左右端を決定することができる。さらに得られた左右端の範囲での垂直方向の平均プロファイルを作成し、水平方向と同様な処理を行い、図中のX,Y,Qを検出することで上下端を決定できる。このほか診断上重要な部分は画像中央に撮影されることことを考慮し関心領域画像中央に設定する、頭部などのように診断上重要な部分が広範囲に渡る場合は照射野領域全体を関心領域に設定する方法などが考えられる。
【００９５】
（４）特徴量算出：
ここで、特徴量算出手段１３０は、放射線画像内の特徴量を算出する。この特徴量として、たとえば、画像に微分フィルタあるいはラプラシアンフィルタなどの高周波領域を抽出するフィルタ処理を行うことによって得られるエッジ成分の絶対値を画像のエッジ領域のコントラストとしてとることができる。
【００９６】
図４（ａ）は頚椎についてのオリジナルの画像、図４(ｂ)は図４（ａ）の頚椎画像について抽出されたエッジを示した様子である。
【００９７】
この場合、コントラストの大きさ、すなわちエッジ成分の絶対値の大きさにしたがって点数をつけ、点数がどの画素値にどのように分布しているかのヒストグラムを作成することで、コントラスト（エッジ）の大きさを各画素への重みとして用いることができる。ここで点数のつけ方については、たとえば絶対値がある経験的に定められた閾値A以下の場合は絶対値の値を点数とし、A以上の場合は０とすることで人体から構成されると予想されるエッジのみにエッジに応じた点数を与えることができる。この他、閾値Aを用いて、図８のような関数を用いて滑らかに変化させてもよい。
【００９８】
また、たとえば、エッジ領域の抽出にはラプラシアンピラミッド法やウェーブレットなどの多重解像度法を用いることもできる。これらの多重解像度法を用いた高周波成分抽出については、たとえば特開平５−２４４５０８号公報や特開平６−２７４６１４号公報に詳しい記載がある。
【００９９】
上記コントラストにより算出された重み付けの値を画像内のすべてのエッジ領域について総和した値が特徴量に相当する。このほか、特徴量としては、たとえば、画像平均濃度や、画像内の画素信号値の分散値や中央度（どれだけ中央（画像中心）に近いかの度合い）、なども考えられ、このほか
(a) 領域内の画素信号値の最大値もしくは最小値のいずれか、
(b) 領域内の画素信号値のメジアン、
(c) 領域内の画素信号値の最頻値（モード）、
(d) 領域内の画素信号値の分散もしくは標準偏差のいずれか一方、
(e) 隣接あるいは所定間隔の画素の画素値差分、
なども考えられる。
【０１００】
（５）特徴量統合：
ここで特徴量統合手段１３０では、特徴量算出手段１２０で算出された複数の特徴量のなかで、より診断に重要な組み合わせが重要視されるように結合される。これにより、診断に必要な領域を重視した処理が可能となる。
【０１０１】
より具体的には例えば、たとえば、(x,y)を画素、P(x,y)を画素値とし、画像エッジ部のコントラストW_1(x,y)、画像中央度W_2(x,y)とした場合、結合重みをA(P,W_1,W_2)とすれば特徴量Cを、
C=ΣA(P,W_1,W_2)P(x,y)+A(P,W_1,W_2)W_1(x,y)+A(P,W_1,W_2)W_2(x,y)
で与えることができる。ここでΣはエッジ領域すべてでの総和を表すものとする。AをたとえばP1<P<P2,W1<W_1<W2,W3>W_2を満たす領域でのみ２を取り、その他の領域では１をとるようにすれば、信号値が人体領域と推定されるP1〜P2の範囲かつ、照射野のエッジなどを除いた人体により形成されるエッジの範囲W1〜W2にエッジがあり、ある一定値W3よりが画像中央にある画素にのみ、２倍の得点を配点することができる。これらのP1,P2,W1〜W3は判別分析法などを使用し、経験的に求められる値である。この結合には上記のような方法のほかファジイ積分を用いた手法も考えられる。
【０１０２】
（６）特徴量評価：
ここで、特徴量評価手段１４０は、特徴量算出手段１２０で算出され、特徴量統合手段１３０で統合された特徴量を、特徴量評価関数により評価する。
【０１０３】
以下、ルックアップテーブルの使用効率を最大にする階調処理の例について説明する。この特徴量評価関数Eとしてのは、基本式として、
Ｅ＝Σｆ（Δｘ′／Δｘ）Ｗ（ｘ）、
と示すことができる。
【０１０４】
ここで、ｆは特徴量評価関数、Δｘはｘにおける画素値の微小変化であり、具体的にはある一定値δをｘを中心に可算、減算したもので区間(ｘ−δ,ｘ＋δ)の長さを表す。またΔｘ'は階調特性を決定するルックアップテーブルであるＬＵＴ通過後のXの値ｘ'の微小変化で具体的にはＬＵＴによるXの変換後の値をＬＵＴ(ｘ)で表せば区間(ＬＵＴ(ｘ−δ),ＬＵＴ(ｘ＋δ))の長さを表す。また、Ｗ（ｘ）はｘの重みである。図５（ａ）では、ｘ１におけるΔｘ１とΔｘ１′、ｘ２におけるΔｘ２とΔｘ２′、の例を示している。このような特徴量評価関数ｆを、エッジと認められる全ての点で行ってΣを算出する。
【０１０５】
なお、特徴量評価関数ｆとしては、
Δｘ′／Δｘ＝１０では１０、
Δｘ′／Δｘ＝５では５、
Δｘ′／Δｘ＝１では０、
Δｘ′／Δｘ＝１／５では−５、
Δｘ′／Δｘ＝１／１０では−１０、
などとすることができる。このような特徴量評価関数ｆとしては、たとえば、図６（ａ）あるいは図６（ｂ）のような特性になる。
【０１０６】
また、たとえば、階調特性を決定するルックアップテーブル（図５参照）のシフト値Ｓ、グラデーション値Ｇを変数として、あるＳＧ値のルックアップテーブルで変換した後の画像からエッジ検出処理を行い、その比を取る方法が考えられる。このＳＧ値は具体的には基本となるＬＵＴに対して、ＬＵＴ(Ａｘ−Ｂ)により画像を変換した場合、グラデーション値ＧをＡ、シフト値ＳをＢと定義する。
【０１０７】
すなわち、評価関数をＥとすると、
Ｅ（Ｓ，Ｇ）＝ＥＤＧＥ（Ｓ，Ｇ）／ＯＲＧ＿Ｅ（Ｓ，Ｇ）
とすればよい。
【０１０８】
ここで、ＥＤＧＥ（Ｓ，Ｇ）はルックアップテーブル変換後の画像の平均コントラスト、ＯＲＧ＿Ｅ（Ｓ，Ｇ）はオリジナルの画像の平均コントラストである。このコントラストはたとえば、直接画像から求めるものではなく、入力は常に一定値とし、あるしきい値を超えるエッジ点でのルックアップテーブル通過後、この一致値の振幅がどう変わったかを評価するものでも良い。すなわち、以下のようである。
Ｅ（Ｓ，Ｇ）＝Σ｛（ＬＵＴ（Ｓ，Ｇ，ｘ＋Ａ）−ＬＵＴ（Ｓ，Ｇ，ｘ−Ａ））／（２＊Ａ）＊Ｗ（ｘ）｝…（Ａ）、で与えられる。ここでＡは画素値ｘにおける微小区間を表す定数で例えば１２ビット階調の画像でに対して５などの非常に小さい値が選ばれる。またＷ(ｘ)は各画素に対する重要度をあらわした重みであり、ＬＵＴ(Ｓ，Ｇ，ｘ)は与えられたＳＧ値でのルックアップテーブルで画像内の点ｘを変換した際の出力値を示し、Σは検出されたエッジ点の全ての領域での和を表す。ここでは、コントラストの評価を例にしたが、同様に変換後の画像と変換前の画像の分散の比を評価関数とすることもできるし、同様に変換後の画像の平均濃度を評価関数とすることもできる。
【０１０９】
たとえば、分散の比の場合、
Ｅ（Ｓ）＝Σ｛ＶＡＲ（Ｓ）／ＯＲＧ＿Ｖ（Ｓ）｝
となる。ここで、ＶＡＲ（Ｓ）は変換後の平均分散、ＯＲＧ＿Ｖ（Ｓ）はオリジナル画像の分散とし、Σは画像のある領域内の点における和を表す。
【０１１０】
さらに、たとえば、平均の比の場合、
Ｅ（Ｓ）＝ＡＶＥ（Ｓ，Ｇ）−Ａ
となる。ここで、ＡＶＥ（Ｓ，Ｇ）は変換後の平均信号値、Ａは定数である。
【０１１１】
（７）パラメータ決定：
ここで、パラメータ決定手段１５０は、特徴量評価手段１４０で特徴量評価関数により評価された結果であるＥ（Ｓ，Ｇ）が最大あるいは最小になるようなＳ値、Ｇ値を求める。
【０１１２】
このようにすることで、たとえば、特徴量をコントラストとすれば、Ｅ（Ｓ，Ｇ）が最大になるようにすることで、変換後の画像の全体としての重要領域のコントラスト増幅率を最大になるようにすることができる。
【０１１３】
また、たとえば、特徴量を分散とすれば、評価関数が最小になるようにすることで、ある領域内の画像のばらつきが最も小さくなるようにＳ値を決定できる。
【０１１４】
また、たとえば、  特徴量が平均値の場合には、評価関数が最小になるようにすることで、ある領域内の画像画像平均値が特定の定数Ａと等しくなるようにできる。
【０１１５】
さらに、これらの関数を組み合わせることで、評価手段１４０において別の評価関数を作成することも可能である。
【０１１６】
また、以上の特徴量評価結果を単独に用いるのではなく、複数を組み合わせて、変換後の画像の全体として重要領域のコントラスト増幅率、平均値を所望の状態にすることもできる。具体的にはコントラスト増幅率の評価関数をＥＣ（Ｓ，Ｇ）、平均値をＥＡ（Ｓ）とすれば、
Ｅ（Ｓ，Ｇ）＝ＥＡ（Ｓ）＊ＥＣ（Ｓ，Ｇ）、
を評価関数とすることで、コントラストが大きく濃度が所定値Ａに近いほど高い得点を算出することができ、したがって、平均画素値がＡに近くコントラストの増幅率が高いＬＵＴを決定することができる。
【０１１７】
なお、以上の説明で、パラメータライズされた変換方法は、階調処理のためのルックアップテーブルによる画素値変換の他に、周波数強調処理やイコライゼーション処理などが考えられる。これらの処理については特開２００１−１２０５２４号公報、特開昭６２−６２３７３号公報、特開平９−４４６４５号公報や特許公報２６６３１８号などに詳しい記載がある。評価関数として、処理前と処理後の分散の比やエッジサイズの変化率などをとることで、これらの処理に対して適切なパラメータを決定することが可能となる。例えば、周波数強調処理は強調係数β、イコライゼーション処理における信号補正関数の傾きの絶対値Kなどによって処理が決定されるが、周波数強調処理の評価関数として
ＥＦ(β)={ＥＤＧＥ(β)／ＥＤＧＥ_ORG(β)}／Σ{ＶＡＲ(β)／ＯＲＧ＿Ｖ(β)}-Ａ、
選べば、エッジの増幅度と画素分散値の比が特定値Ａにつねに近くなるように強調度を決定することができ、エッジの強調と画像のざらつきをコントロールすることが可能となる。ここでＥＤＧＥ(β)はオリジナル画像の平均コントラスト、ＥＤＧＥ_ORG(β)は周波数強調処理後の平均コントラスト、ＶＡＲ(β)は周波数強調後の画像平均分散値の ＯＲＧ＿Ｖ(β)はオリジナル画像の平均信号値、Ａは経験的に定まる定数である。
【０１１８】
さらにイコラゼーション処理の場合、Ｄ_ORG(Ｋ)をオリジナル画像のダイナミックレンジ（画像内の最大信号値-最小信号値）、Ｄ(Ｋ)をイコライゼーション処理後のダイナミックレンジ、Ｌ(Ｋ,ｘ)を画素ｘのイコライゼーション処理後の値をＬＵＴにより変換した場合の値としたとき、
ＥＦ(Ｋ)=Σ{Ｌ(Ｋ,ｘ＋Ａ)-Ｌ(Ｋ,ｘ−Ａ)／２＊Ａ}/α{Ｄ(Ｋ)/Ｄ_ORG(Ｋ)}、
で与える。ここでΣはたとえば検出されたエッジ領域すべてにおける総和を表すものとする。このようにすることで画像全体のコントラストを考慮し、かつ診断上重要な領域のイコライゼーション処理後のコントラスト増幅率が適切になるように補正の強度を決定することができる。上式のαは適切な値が算出されるようにするための補正係数で経験的に求められる値である。
【０１１９】
（８）画像処理：
画像処理手段１６０は、画像データ生成手段１１０からの画像データについて、以上の説明のようにしてパラメータ決定手段１５０で決定されたパラメータにしたがって、画像処理を実行する。
【０１２０】
具体的には例えば階調処理においては、上記の評価関数（Ａ）により決定されたＳＧ値とＬＵＴを使用して、画像データの各画素XをＬＵＴ（Ｇｘ−Ａ）によって変換する。このようにすることで、処理済み画像は与えられた評価関数、すなわち上記例では変換後の重要画素のコントラスト増大率が最大になるように変換された画像となる。また、周波数強調処理やイコライゼーション処理においても同様に上記のような評価関数により決定した強調係数あるいは補正係数により、処理が行われる。この処理済み画像データは、図示されない画像表示手段、各種外部機器、画像出力装置などに供給される。
【０１２１】
以上のように、放射線画像毎の特徴量に基づいて評価して決定したパラメータにより画像処理を実行するので、例えば、各信号値に対して特徴量に基づいて、重要度に応じた重み付けをして、各画素の階調処理用ＬＵＴにより変換した際の信号の増幅率を算出し、増幅率と重みを掛け合わせた値のすべての画素での総和が最大になるようなＬＵＴを決定することで適切な画像処理が実行できる。ここで特徴量に基づいた重要度とは、たとえば一般に診断上重要なポイントは画像中央に位置することが多く、さらにエッジ成分をもつことなどから、特徴量として画像中心との距離を測る画像中央度と微分フィルタなどを用いて検出した各画素におけるエッジ強度を使用し、それらの値に基づいて決定される量である。また、増幅率としてはたとえば各画素値におけるＬＵＴの傾きやそれに相当する量を考えることができる。
【０１２２】
この結果、画像処理条件決定の失敗により診断が困難となる状況の発生は防止され、画像データ毎に適切なパラメータが決定されるようになるので、いかなる場合であっても適切な条件あるいは適切に近い条件で画像処理することが可能になる。
【図面の簡単な説明】
【０１２３】
【図１】本発明の一実施の形態例の全体構成を機能的に示す機能ブロック図である。
【図２】本発明の一実施の形態例の全体処理の流れを示すフローチャート図である。
【図３】本発明の実施の形態例における処理の様子を示す説明図である。
【図４】本発明の実施の形態例における処理の様子を示す説明図である。
【図５】本発明の実施の形態例における処理の様子を示す説明図である。
【図６】本発明の実施の形態例における処理の様子を示す説明図である。
【図７】本発明の実施の形態例における処理の様子を示す説明図である。
【図８】本発明の実施の形態例における処理の様子を示す説明図である。
【符号の説明】
【０１２４】
５  被写体
３０  放射線発生装置
４０  放射線画像読取器
１００  画像処理装置
１０１  制御手段
１１０  画像データ生成手段
１２０  特徴量算出手段
１３０  特徴量統合手段
１４０  特徴量評価手段
１５０  パラメータ決定手段
１６０  画像処理手段
【図１】
2005109867.tif 000003
【図２】
2005109867.tif 000004
【図３】
2005109867.tif 000005
【図４】
2005109867.tif 000006
【図５】
2005109867.tif 000007
【図６】
2005109867.tif 000008
【図７】
2005109867.tif 000009
【図８】
2005109867.tif 000010