(57)【要約】【課題】  元の画像と同じ階調に感じられるようにしたまま、ヒストグラムエクスパンドによる効果的な拡大表示を行う。【解決手段】  拡大表示する領域の部分画像をモニターに表示した時の輝度範囲に対応するＪＮＤの範囲をモニターの全輝度範囲に対応するＪＮＤの範囲に変換することによって人間の目に同じように感じる階調変換テーブルを作成し、この階調変換テーブルによって部分画像を階調変換して拡大表示する。【選択図】  図１
2005295417.tif 000002
【特許請求の範囲】
【請求項１】
拡大前と拡大後で画像のコントラストが人間の目に同じように感じるように階調変換を行って拡大表示する画像処理装置。
【請求項２】
請求項１においてＪＮＤを介して人間の目に同じように感じる階調変換テーブルを作成して階調変換を行って拡大表示する画像処理装置。
【請求項３】
請求項２において拡大表示する領域の部分画像をモニターに表示した時の輝度範囲に対応するＪＮＤの範囲をモニターの全輝度範囲に対応するＪＮＤの範囲に変換することによって、
人間の目に同じように感じる階調変換テーブルを作成する画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はディジタル画像の拡大表示に係るもので、特にＸ線ディジタル画像のようにダイナミックレンジの広いディジタル画像の拡大表示に関するものである。
【背景技術】
【０００２】
ディジタル画像処理において、対象となる画像のヒストグラムをとり、そのヒストグラムにおける画素値の最小値から最大値までの範囲を広げることによって画像の強調を行うヒストグラムエクスパンドの手法が知られている。
【０００３】
（例えば、「画像処理ハンドブック」画像処理ハンドブック編集委員会編  昭晃堂  ｐ２６３〜ｐ２６４参照）。
【０００４】
このヒストグラムエクスパンドの手法を拡大表示に用いた方法が例えば特開平５−２２８１３３に開示されている。
【０００５】
これは拡大表示する範囲の部分画像のヒストグラムをとり、そのヒストグラムにおける画素値の最小値から最大値までの範囲を広げることによって拡大画像の強調を行い、より効果的な拡大表示を行うものである。
【発明の開示】
【発明が解決しようとする課題】
【０００６】
人間の視覚は、輝度によりコントラストの感じ方が変わることが知られている。
【０００７】
このような感度の違いをＢａｒｔｅｎモデルを基に定義したものとしてＪＮＤ（Ｊｕｓｔ  Ｎｏｔｉｃｅａｂｌｅ  Ｄｉｆｆｅｒｅｎｃｅ）が知られており、医療の標準的な規格であるＤＩＣＯＭ規格で用いられている。
【０００８】
このＪＮＤは平均的な人が識別可能な輝度の差を表わすものである。
【０００９】
ＪＮＤは輝度範囲０．０５ｃｄ／ｍ２〜４０００ｃｄ／ｍ２において定義され、この輝度範囲においてＪＮＤは１０２３となる。つまり、０．０５ｃｄ／ｍ２〜４０００ｃｄ／ｍ２の輝度範囲においては、平均的な人の目で１０２３段階の輝度差を識別できるということになる。
【００１０】
このＪＮＤのインデックスをｊ、輝度をＬとした時、ｊからＬへの関数Ｌ（ｊ）は次式によって表わすことができる。
【００１１】
【数１】
2005295417.tif 000003
【００１２】
ここで、Ｉｎは自然対数を表し、
ａ＝−１．３０１１８７７，ｂ＝−０．０２５８４０１９，ｃ＝０．０８０２４２６３６，ｄ＝−０．１０３２０２２９，ｅ＝０．１３６４６６９９，ｆ＝０．０２８７４５６２０，ｇ＝−０．０２５４６８４０４，ｈ＝−０．００３１９７８９７７，ｋ＝０．０００１２９９２６３４，ｍ＝０．００１３６３５３３４
である。
【００１３】
この関数のグラフを図７に示す。
【００１４】
逆にＬからｊへの関数ｊ（Ｌ）は次式で表わすことができる。
【００１５】
ｊ（Ｌ）＝Ａ＋Ｂ・ｌｏｇ１０（Ｌ）＋Ｃ・（ｌｏｇ１０（Ｌ））２＋
Ｄ・（ｌｏｇ１０（Ｌ））３＋Ｅ・（ｌｏｇ１０（Ｌ））４＋Ｆ・（ｌｏｇ１０（Ｌ））５＋Ｇ・（ｌｏｇ１０（Ｌ））６＋Ｈ・（ｌｏｇ１０（Ｌ））７＋Ｉ・（ｌｏｇ１０（Ｌ））８
ここで、ｌｏｇ１０は常用対数を表し、
Ａ＝７１．４９８０６８，Ｂ＝９４．５９３０５３，Ｃ＝４１．９１２０５３，Ｄ＝９．８２４７００４，Ｅ＝０．２８１７５４０７，Ｆ＝−１．１８７８８４５５，Ｇ＝−０．１８０１４３４９，Ｈ＝０．１４７１０８９９，Ｉ＝−０．０１７０４６８４５
である。
【００１６】
この関数のグラフを図８に示す。
【００１７】
以上のように人の目は輝度によってコントラストの感じ方が変わってしまうことがわかる。したがって、拡大する時にヒストグラムエクスパンドを行えば、画像全体の輝度が変わってしまい拡大前のコントラストと拡大後のコントラストが人の目には違って感じられてしまうという問題点がある。
【課題を解決するための手段】
【００１８】
このような問題点を解決するために、拡大前と拡大後で画像のコントラストが人間の目に同じように感じるように階調変換を行って拡大表示する。
【００１９】
このための処理の概念図を図２に、フローチャートを図３に示す。
【００２０】
この処理手順は、
（１）拡大表示する領域を設定する。
【００２１】
（２）切り出された部分画像のヒストグラムを求める。
【００２２】
（３）（２）で求めたヒストグラムの画素値範囲（最小値、最大値）を求める。
【００２３】
（４）モニターの表示特性より（３）の画素値範囲に対応する輝度範囲を求める。
【００２４】
（５）輝度−ＪＮＤ（Ｊｕｓｔ  Ｎｏｔｉｃｅａｂｌｅ  Ｄｉｆｆｅｒｅｎｃｅ）テーブルを用いて、（４）の輝度範囲を表示装置の最大輝度範囲に広げたときに同一に見えるための画素値変換テーブルを作成する。
【００２５】
（６）（５）の画素値変換テーブルを用いて部分画像を変換して拡大表示する。
である。
【００２６】
この拡大表示を部分画像の切り出し範囲のＲＯＩをマウスでリアルタイムに動かしながら行えばより効果的である。以上のようにして元の画像とコントラストは同じに感じるようにしたまま、ヒストグラムエクスパンドによる強調された効果的な拡大表示を実現することができる。
【発明の効果】
【００２７】
常に元の画像と同じ階調に感じられるようにしたまま、ヒストグラムエクスパンドにより表示装置の最大限のコントラストを用いながら拡大表示を行うことができる。このため、違和感のない効果的な拡大表示となり、画像の正確で精密な観察が可能となる。
【発明を実施するための最良の形態】
【００２８】
［実施例］
図２は本発明による放射線画像処理システムの一例の全体構成を示す概略図である。
【００２９】
４．Ｘ線発生装置制御部により制御された３．Ｘ線源より発生したＸ線は２．患者を透過して１．Ｘ線センサーにより検知される。検知されたＸ線はディジタルＸ線画像として５．画像入力部に入力される。入力されたディジタルＸ線画像は７．画像処理部によってＸ線センサーの補正処理、階調処理、空間フィルタリング処理等の画像処理がなされる。画像処理のなされたディジタルＸ線画像は９．診断モニターに表示されたり８．画像保存部に保存されたり１１．ネットワークを介して１２．プリンター、１３．診断ワークステーション、１４．画像データベースに出力されたりする。表示、出力された画像が満足のいくものでなかった場合には画像処理パラメータを変えるなどしながら画像処理、表示を繰り返し行う。以上の操作は１０．操作部によって行われる。
【００３０】
このようなシステムで動作する実施例について説明する。
【００３１】
まず、拡大表示をするための範囲を設定する。設定はあらかじめ表示されている原画像上に矩形のＲＯＩ（関心領域）をマウスなどのポインティングデバイスにより指定する。拡大範囲が指定されたら原画像からこの範囲を切り出す。そして切り出した部分画像のヒストグラムを求める。このヒストグラムより部分画像の画素値範囲を決定する。画素値範囲を決定する最も単純な方法は、ヒストグラムにおける画素値の最小値から最大値までを画素値範囲とする方法である。ただし、ノイズや異常値が部分画像内に含まれていた場合、単純に最小値から最大値とするとこのノイズや異常値のために実質的な画素値範囲よりも広めの画素値範囲としてしまうことがある。このような事態を避けるためには図４に示すように画素値の大きい方と小さい方から例えば全体の画素数の１％ずつを除外した範囲を画素値範囲とすれば良い。
【００３２】
部分画像の画素値範囲が求まったら、この画素値範囲の画像をモニターに表示した時の輝度範囲を求める。これは図５に示すようにモニター表示特性の画素値範囲に対応する輝度範囲によって求めることができる。
【００３３】
輝度範囲が求まったら、輝度−ＪＮＤ変換関数ｊ（Ｌ）及びＪＮＤ−輝度変換関数Ｌ（ｊ）を用いて、部分画像の輝度範囲を表示装置の最大輝度範囲に広げたときに同一に見えるための画素値変換テーブルを作成する。この画素値変換テーブルを作成する方法の概念図を図６に示す。６０１．ではまず拡大する部分画像の画素値をモニターの表示特性である画素値→輝度変換テーブルによって輝度に変換する。このとき、６０６．部分画像の画素値範囲は原画像の全体画素値範囲よりも狭いから当然その輝度範囲（６０７．）もモニター全体の輝度範囲よりも狭くなっている。次にこの部分画像の輝度をＪＮＤに変換する。このＪＮＤへの変換には輝度−ＪＮＤ変換関数ｊ（Ｌ）を用いる。このＪＮＤをモニター全体の輝度範囲に対応するＪＮＤに変換する（６０３．）。そして今度は逆にＪＮＤ−輝度変換関数Ｌ（ｊ）を用いて、モニター全体の輝度範囲に対応するＪＮＤから輝度に変換する（６０４．）。この輝度に対応する画素値を求めれば部分画像の輝度範囲を表示装置の最大輝度範囲に広げたときに同一に見えるための画素値の対応が決定する。これをすべての画素値について対応関係を求めれば画素値変換テーブルを作成することができる。
【００３４】
このようにして画素値変換テーブルができたら、この画素値変換テーブルを用いて部分画像を変換して拡大表示する。
【００３５】
以上のようにして原画像のコントラストが変わらなく感じるように画像を拡大して表示することができる。
【図面の簡単な説明】
【００３６】
【図１】システム構成例を表す図
【図２】処理の流れを説明する概念図
【図３】処理の流れを示すフローチャート
【図４】画素値範囲を説明する概念図
【図５】モニター表示特性を説明する概念図
【図６】画素値の変換の流れを説明する概念図
【図７】ＪＮＤ−輝度変換テーブルを示すグラフ
【図８】輝度−ＪＮＤ変換テーブルを示すグラフ
【図１】
2005295417.tif 000004
【図２】
2005295417.tif 000005
【図３】
2005295417.tif 000006
【図４】
2005295417.tif 000007
【図５】
2005295417.tif 000008
【図６】
2005295417.tif 000009
【図７】
2005295417.tif 000010
【図８】
2005295417.tif 000011