JP2004139464A - 画像処理装置、画像処理方法、プログラム及び記録媒体 - Google Patents
画像処理装置、画像処理方法、プログラム及び記録媒体 Download PDFInfo
- Publication number
- JP2004139464A JP2004139464A JP2002305024A JP2002305024A JP2004139464A JP 2004139464 A JP2004139464 A JP 2004139464A JP 2002305024 A JP2002305024 A JP 2002305024A JP 2002305024 A JP2002305024 A JP 2002305024A JP 2004139464 A JP2004139464 A JP 2004139464A
- Authority
- JP
- Japan
- Prior art keywords
- gradation conversion
- image
- image processing
- variable
- target image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012545 processing Methods 0.000 title claims abstract description 45
- 238000003672 processing method Methods 0.000 title claims description 25
- 238000006243 chemical reaction Methods 0.000 claims abstract description 148
- 238000000034 method Methods 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 16
- 230000006870 function Effects 0.000 description 63
- 238000011156 evaluation Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 210000000038 chest Anatomy 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010230 functional analysis Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 210000001370 mediastinum Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/92—Dynamic range modification of images or parts thereof based on global image properties
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
【課題】特定の階調変換関数を用いつつ、変換後画像の画素値ヒストグラムの平坦度が良好な階調変換を行うことのできる画像処理装置等を提供すること。
【解決手段】階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換する画像処理装置等とする。
【選択図】 図７
【解決手段】階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換する画像処理装置等とする。
【選択図】 図７
Description
【０００１】
【発明の属する技術分野】
本発明は、画像の階調変換処理を行うための画像処理装置、画像処理方法、プログラム及び記録媒体に関するものである。
【０００２】
【従来の技術】
画像全体のコントラストを強調するための代表的な手法としてヒストグラムイコライゼーションが知られている。これは階調変換を行なった後のヒストグラムの分布が平坦（均一）になるように階調変換関数を決める方法である（詳細なアルゴリズムは「ディジタル画像処理：Ｒｏｓｅｎｆｅｌｄ，Ｋａｋ著、長尾真監訳」参照）。ヒストグラムイコライゼーションを利用した階調変換自動化装置が特許第２８０８７７３号公報により知られている。尚、階調変換関数は階調変換曲線とも称され、実際の階調変換には階調変換関数に従ったルックアップテーブルが用いられる場合が多い。
【０００３】
また、Ｘ線画像も従来はフィルム画像が主であったが、近年では大画面エリアセンサを用いたディジタルＸ線撮影システムや、輝尽性蛍光体を用いたＣＲシステムなどのディジタル化が進んできている。このＸ線ディジタル画像においては、医師は従来のフィルム画像で診断してきているためフィルムと同等の画像を望む傾向がある。このためフィルムの特性曲線を階調曲線に用いて階調変換処理を行なうことによってフィルムと同等の画像を得ていた。
【０００４】
【発明が解決しようとする課題】
しかしながら、ヒストグラムイコライゼーションによる階調変換処理においては、変換後の画像のヒストグラムが平坦になることによってのみ階調曲線が決められる。したがって、その階調曲線の滑らかさは全く保障されていない。このため、Ｘ線画像の場合には従来のフィルム画像とかけ離れた画像となる場合も多く、従来のフィルム画像で診断してきた医師にとって診断しづらい画像となってしまうという問題がある。
【０００５】
本発明は、上記問題点に鑑みてなされたものであり、特定の階調変換関数を用いつつ、変換後画像の画素値ヒストグラムの平坦度が良好な階調変換を行うことのできる画像処理装置、画像処理方法、プログラム及び記録媒体を提供することを目的とする。
【０００６】
【課題を解決するための手段】
斯かる目的を達成するための本発明の第１の実施態様は、対象画像の階調変換処理を行う画像処理装置であって、階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを特徴とする画像処理装置である。
【０００７】
第２の実施態様は、対象画像の階調変換処理を行う画像処理装置に用いられる画像処理方法であって、階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを特徴とする画像処理方法である。
【０００８】
第３の実施態様は、対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムであって、前記所定の画像処理方法は、階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを含むことを特徴とする特徴とするプログラムである。
【０００９】
第４の実施態様は、対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムを記憶したコンピュータ読み取り可能な記録媒体であって、前記所定の画像処理方法は、階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを含むことを特徴とするコンピュータ読み取り可能な記録媒体である。
【００１０】
本発明の他の目的、特徴若しくは優位性等は、添付図面を参照してなされた後述の説明から明らかにされるが、以下に代表的な本発明の他の実施態様をいくつか示す。
第５の実施態様は、対象画像の階調変換処理を行う画像処理装置であって、階調変換特性を特定するための変数を用いて表わされる階調変換関数により対象画像を階調変換することにより、該変数を用いて表わされる画素値からなる画像を得る画像変換手段と、前記画像変換手段により得られた画像の画素値ヒストグラムを作成するヒストグラム作成手段と、前記ヒストグラム作成手段により得られた画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、前記階調変換関数を設定する階調変換関数設定手段と、前記階調変換関数設定手段により設定された階調変換関数により前記対象画像を階調変換する階調変換手段とを有することを特徴とする前記第１の実施態様に従った画像処理装置である。
【００１１】
第６の実施態様は、前記階調変換がルックアップテーブルを用いて行われることを特徴とする前記第１又は第５の実施態様に従った画像処理装置である。
【００１２】
第７の実施態様は、所定の評価基準に基づき前記画素値ヒストグラムの平坦度が最大となるように前記変数値が決定されることを特徴とする前記第１又は第５の実施態様に従った画像処理装置である。
【００１３】
第８の実施態様は、前記所定の評価基準は前記画素値ヒストグラムと平坦な画素値ヒストグラムとの平均二乗誤差であり、前記変数値は前記平均二乗誤差が最小となるように決定されることを特徴とする前記第７の実施態様に従った画像処理装置である。
【００１４】
第９の実施態様は、対象画像の階調変換処理を行う画像処理装置に用いられる画像処理方法であって、階調変換特性を特定するための変数を用いて表わされる階調変換関数により対象画像を階調変換することにより、該変数を用いて表わされる画素値からなる画像を得る画像変換工程と、前記画像変換工程において得られた画像の画素値ヒストグラムを作成するヒストグラム作成工程と、前記ヒストグラム作成工程において得られた画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、前記階調変換関数を設定する階調変換関数設定工程と、前記階調変換関数設定工程において設定された階調変換関数により前記対象画像を階調変換する階調変換工程とを含むことを特徴とする第２の実施態様に従った画像処理方法である。
【００１５】
第１０の実施態様は、前記階調変換がルックアップテーブルを用いて行われることを特徴とする前記第２又は第９の実施態様に従った画像処理方法である。
【００１６】
第１１の実施態様は、所定の評価基準に基づき前記画素値ヒストグラムの平坦度が最大となるように前記変数値が決定されることを特徴とする前記第２又は第９の実施態様に従った画像処理方法である。
【００１７】
第１２の実施態様は、前記所定の評価基準は前記画素値ヒストグラムと平坦な画素値ヒストグラムとの平均二乗誤差であり、前記変数値は前記平均二乗誤差が最小となるように決定されることを特徴とする前記第１１の実施態様に従った画像処理装置である。
【００１８】
第１３の実施態様は、対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムであって、前記所定の画像処理方法は、階調変換特性を特定するための変数を用いて表わされる階調変換関数により対象画像を階調変換することにより、該変数を用いて表わされる画素値からなる画像を得る画像変換工程と、前記画像変換工程において得られた画像の画素値ヒストグラムを作成するヒストグラム作成工程と、前記ヒストグラム作成工程において得られた画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、前記階調変換関数を設定する階調変換関数設定工程と、前記階調変換関数設定工程において設定された階調変換関数により前記対象画像を階調変換する階調変換工程とを含むことを特徴とする前記第３の実施態様に従ったプログラムである。
【００１９】
第１４の実施態様は、前記階調変換がルックアップテーブルを用いて行われることを特徴とする前記第３又は第１３の実施態様に従ったプログラムである。
【００２０】
第１５の実施態様は、所定の評価基準に基づき前記画素値ヒストグラムの平坦度が最大となるように前記変数値が決定されることを特徴とする前記第３又は第１３の実施態様に従ったプログラムである。
【００２１】
第１６の実施態様は、前記所定の評価基準は前記画素値ヒストグラムと平坦な画素値ヒストグラムとの平均二乗誤差であり、前記変数値は前記平均二乗誤差が最小となるように決定されることを特徴とする前記第１５の実施態様に従ったプログラムである。
【００２２】
第１７の実施態様は、対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムを記憶したコンピュータ読み取り可能な記録媒体であって、前記所定の画像処理方法は、階調変換特性を特定するための変数を用いて表わされる階調変換関数により対象画像を階調変換することにより、該変数を用いて表わされる画素値からなる画像を得る画像変換工程と、前記画像変換工程において得られた画像の画素値ヒストグラムを作成するヒストグラム作成工程と、前記ヒストグラム作成工程において得られた画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、前記階調変換関数を設定する階調変換関数設定工程と、前記階調変換関数設定工程において設定された階調変換関数により前記対象画像を階調変換する階調変換工程とを含むことを特徴とする前記第４の実施態様に従ったコンピュータ読み取り可能な記録媒体である。
【００２３】
第１８の実施態様は、前記階調変換がルックアップテーブルを用いて行われることを特徴とする前記第４又は第１７の実施態様に従ったコンピュータ読み取り可能な記録媒体である。
【００２４】
第１９の実施態様は、所定の評価基準に基づき前記画素値ヒストグラムの平坦度が最大となるように前記変数値が決定されることを特徴とする前記第４又は第１７の実施態様に従ったコンピュータ読み取り可能な記録媒体である。
【００２５】
第２０の実施態様は、前記所定の評価基準は前記画素値ヒストグラムと平坦な画素値ヒストグラムとの平均二乗誤差であり、前記変数値は前記平均二乗誤差が最小となるように決定されることを特徴とする前記第１９の実施態様に従ったコンピュータ読み取り可能な記録媒体である。
【００２６】
【発明の実施の形態】
以下、本発明の実施形態を添付図面を参照しながら詳細に説明する。
図１は本発明の一実施形態に係る画像処理装置を利用した放射線画像処理システムの全体構成を示す概略図である。
【００２７】
本実施形態に係る画像処理装置は、図１の一点鎖線内に示されるように、画像入力部５、Ｘ線撮影システム制御部６、画像処理部７及び画像保存部８で構成される。Ｘ線発生装置制御部４により制御されたＸ線源３より発生したＸ線は患者２を透過してＸ線センサ１により検知される。検知されたＸ線はディジタルＸ線画像として画像入力部５に入力される。入力されたディジタルＸ線画像は画像処理部７によってＸ線センサの補正処理、階調処理、空間フィルタリング処理等の画像処理がなされる。画像処理のなされたディジタルＸ線画像は診断モニタ９に表示されたり画像保存部８に保存されたりネットワーク１１を介してプリンタ１２、診断ワークステーション１３、画像データベース１４に出力されたりする。表示、出力された画像が満足のいくものでなかった場合には画像処理パラメータを変えるなどしながら画像処理、表示を繰り返し行う。以上の操作は操作部１０によって行われる。
【００２８】
図２は、このようなシステムで動作する階調変換処理の一例を概念的に示した図である。
ヒストグラムＨ（ｘ）１５は階調変換処理をする前の画像のヒストグラムである。このヒストグラムＨ（ｘ）１５を持つ画像に階調曲線ｆ（ｘ）１６による階調変換処理をかけて、変換された画像のヒストグラムがヒストグラムＨ’（ｘ）１７である。このヒストグラムＨ’（ｘ）１７が全く平坦（均一）である場合には、各画素値をとる画素の数は、画像全体の画素数を当該画像がとり得る画素値の数で割った平均画素数となる。本実施の形態では、このヒストグラムが全く均一である場合（平坦な画素値ヒストグラム）とＨ’（ｘ）１７との平均二乗誤差を求め、これを階調変換後画像の画素値ヒストグラムの平坦度の評価基準として用いる。尚、ここで、平坦度の評価基準は上述のような平均二乗誤差に限られず、階調変換後画像の画素値ヒストグラムと平坦な画素値ヒストグラムとのずれ量の尺度となるものであればよい。尚、このずれ量が小さいほど、平坦度が大きいものとする。
【００２９】
つまり、フィルムの特性曲線と同等になるような制約条件をかけながら、この平均二乗誤差の評価基準が最小となるような階調曲線ｆ（ｘ）１６を求めるものである。これを具体的に説明していく。
【００３０】
まず、ディジタル画像の画素値をχｉとする。画素値のとりうるレベル数をＭとすればχ１≦χｉ≦χＭである。例えば画素値として１２ビット持っている画像の場合にはＭ＝４０９６であり、χ１＝０，χＭ＝４０９５となる。１０ビットであれば、Ｍ＝１０２４であり、χ１＝０，χＭ＝１０２３となる。このディジタル画像のヒストグラムをＨ（χｉ）とする。
【００３１】
また、階調曲線（階調変換関数）をｆ（χｉ：ａ１，ａ２，・・・，ａＮ）とする。ここでｆ（χｉ：ａ１，ａ２，・・・，ａＮ）はＮ個のパラメータ（変数）ａ１，ａ２，・・・，ａＮによって規定される関数であるとする。即ち、階調変換関数ｆの階調変換特性は当該Ｎ個の変数の値により特定される。フィルムの特性曲線と同等の階調曲線の関数の一例として、シグモイド関数を用いた関数を次の式１に示す。
【００３２】
【数１】
【００３３】
この関数で（ａ１，ａ２，ａ３，ａ４，ａＮ）＝（１，５，３，０．２，１／１０２．４）とすれば、１０ビット画像を想定した図３のグラフになる。また、（ａ１，ａ２，ａ３，ａ４，ａＮ）＝（２，５，３，０．２，１／１０２．４）とすれば、図３のグラフのコントラストを倍にした図４のグラフの関数となる。
【００３４】
ただこの関数ではコントラストが左右対称のＳ字型曲線しかできないが、次の式２に示すような関数によって左右のコントラストの違うＳ字型曲線を実現することができる。
【００３５】
【数２】
【００３６】
この関数で（ａ１，ａ２，ａ３，ａ４，ａ５，ａ６，ａ７，ａ８，ａ９）＝（１，５，２，７，１／１０２．４，１／１０２．４，１，１．５，０．２）としたグラフを図５に示す。このような曲線は画素値の高い肺野のコントラストを大きくとって、画素値の低い縦隔でコントラストを小さくするために胸部用に用いられる。尚、ここで縦軸はフィルムの特性曲線に合わせて濃度とした。
【００３７】
このように階調曲線の関数を定義してパラメータをいろいろ変えれば、いろいろな種類のフィルムに対応した階調曲線をつくることができる。このような階調曲線ｆ（χ：ａ１，ａ２，・・・，ａＮ）によって変換した後のディジタル画像（各画素値はａ１，ａ２，・・・，ａＮの関数）のヒストグラムをＨ’（χｉ）とする。このＨ’（χｉ）はＨ’（ｆ（χｉ：ａ１，ａ２，・・・，ａＮ））と表わすことができる。
【００３８】
このＨ’（χｉ）と全画素値にわたって均一なヒストグラムＨ’’（χｉ）＝ｍ（ｍは平均画素数）との２乗誤差の和Ｉｅを求める。このＩｅは次の式３に示すような階調曲線関数のパラメータの関数となる。
【００３９】
【数３】
【００４０】
この２乗誤差の和Ｉｅが最小になるようにパラメータを決定する。このパラメータを制限のないパラメータとして扱ってもフィルムの特性曲線との同等性が保たれるのであれば、このＩｅ（ａ１，ａ２，・・・，ａＮ）を各々のパラメータａ１，ａ２，・・・，ａＮで偏微分して０として次の式４のような連立方程式を立てて、これを解けばパラメータａ１，ａ２，・・・，ａＮを求めることができる。
【００４１】
【数４】
【００４２】
この連立方程式を解いて求めたパラメータａ１，ａ２，・・・，ａＮを用いた階調曲線ｆ（χ：ａ１，ａ２，・・・，ａＮ）はフィルムとの同等性を保ちながら、変換後のヒストグラムの均一性を最大限にする階調曲線となっている。次に、上記連立方程式を解いて求めたパラメータａ１，ａ２，・・・，ａＮを用いた階調曲線ｆ（χ：ａ１，ａ２，・・・，ａＮ））に基づいて、画像入力部５に入力された画像情報の階調変換を行い、診断モニタ９等の情報出力手段にて出力する。
【００４３】
またパラメータに制限を与えないとフィルムの特性曲線との同等性が保たれない場合であるが、この場合で簡単な例として、特定のパラメータのみ制限のないパラメータとして扱い、その他のパラメータは定数で固定することが挙げられる。これは例えば式１でａ１，ａ２のみを変数とし、その他ａ３，ａ４，ａ５を固定する場合である。このようにすると図３、図４のようにＳ字形状を保ちながらコントラスト、横方向のシフトを行なうことができる。
【００４４】
これを用いてＩｅを求めればＩｅはａ１，ａ２の関数となり、上述のようにａ１，ａ２で偏微分して連立方程式を立てて解けばａ１，ａ２を求めることができる。さらに複雑なパラメータの制約条件がいろいろ考えられるが、そのような制約条件のもとで最適化したい場合には、例えば「関数解析による最適理論：Ｄａｖｉｄ Ｇ． Ｌｕｅｎｂｅｒｇｅｒ著、増淵正美、嘉納秀明共訳、コロナ社」ではラグランジュ乗数を用いた方法や不等式拘束条件の場合の方法など種々の方法が記載されており、これを使ってパラメータの最適化を行なえば良い。
【００４５】
尚、以上で説明した階調変換処理は、例えば、画像処理部７を図６のように構成することにより実現できる。図６において、６０は画像処理部７全体の制御を司るＣＰＵ、６１はＣＰＵ６０の動作に必要なプログラム、データ、対象画像等を記憶するための記憶手段としてのＲＡＭ、６２は階調変換関数ｆ（χ：ａ１，ａ２，・・・，ａＮ）により対象画像を階調変換する画像変換部、６３は画像変換部６２により得られた、変数で表わされる画像の画素値ヒストグラムＨ’（χｉ）を作成するヒストグラム作成部、６４はヒストグラム作成部６３により得られた画素値ヒストグラムＨ’（χｉ）と平坦ヒストグラムＨ’‘（χｉ）＝ｍとの２乗誤差を最小化する変数ａ１，ａ２，・・・，ａＮの値を求め、階調変換ルックアップテーブル（ＬＵＴ）を作成する階調変換ＬＵＴ作成部（階調変換関数設定部）、６５は階調変換ＬＵＴ作成部６４により得られた階調変換ＬＵＴを用いて対象画像を階調変換する階調変換部である。６６は上述の各部６０〜６５を通信可能に接続するＣＰＵバスである。
【００４６】
また、６２〜６５の各部は、例えば集積回路等を用いて構成することができる。また、６２〜６５の各部の機能は、ＣＰＵ６０がＲＡＭ６１等に格納された所定のプログラムを実行することによっても実現することができる。この場合、当該プログラムは、例えば図７のフローチャートに示される画像処理方法に対応したプログラムコードを含むものとすればよい。
【００４７】
次に、図７を参照して当該画像処理方法を説明する。まず、対象画像を階調変換関数ｆ（χ：ａ１，ａ２，・・・，ａＮ）により階調変換し、所定の変数で表わされた画素値からなる画像を得る（画像変換工程Ｓ７０）。次いで、前工程において得られた、変数で表わされる画像の画素値ヒストグラムＨ’（χｉ）を作成する（ヒストグラム作成工程Ｓ７１）。続いて、前工程において得られた画素値ヒストグラムＨ’（χｉ）と平坦ヒストグラムＨ’‘（χｉ）＝ｍとの２乗誤差を最小化する変数ａ１，ａ２，・・・，ａＮの値を求め、階調変換ＬＵＴを作成する（階調変換ＬＵＴ作成工程Ｓ７２）。前工程において得られた階調変換ＬＵＴを用いて対象画像を階調変換する（階調変換工程Ｓ７３）。
【００４８】
以上のように、本実施形態によれば、フィルムと同等な画質になるように変換することができる上に、変換後のヒストグラムがなるべく均一になるように効果的にコントラストの割り当てられた階調変換を実現できるため、この階調変換により得られた画像を診断医に提供することによって診断能を向上することができる。
【００４９】
各種デバイスと接続された装置或いはシステム内のコンピュータに対し、ソフトウェアのプログラムコードを供給し、そのシステム或いは装置のコンピュータ（ＣＰＵ或いはＭＰＵ）が、供給されたプログラムに従って上記各種デバイスを動作させることによって上記実施の形態の機能を実現したものも、本発明の範疇に含まれる。
【００５０】
また、この場合、上記ソフトウェアのプログラムコード自体が上述した実施の形態の機能を実現することになり、そのプログラムコード自体は本発明を構成する。そのプログラムコードの供給媒体としては、プログラム情報を所定の信号として伝搬させて供給するためのコンピュータネットワーク（ＬＡＮ、インターネット等のＷＡＮ、無線通信ネットワーク等）システムにおける通信媒体（光ファイバ等の有線回線や無線回線等）を用いることができる。
【００５１】
さらに、上記プログラムコードをコンピュータに供給するための手段、例えばかかるプログラムコードを格納した記憶媒体は本発明を構成する。かかるプログラムコードを記憶する記憶媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
【００５２】
また、供給されたプログラムコードをコンピュータが実行することにより、上述の実施の形態の機能が実現されるだけでなく、そのプログラムコードがコンピュータにおいて稼働しているＯＳ（オペレーティングシステム）或いは他のアプリケーションソフト等と共同して上述の実施の形態の機能が実現される場合にも、かかるプログラムコードが本発明の実施の態様に含まれることはいうまでもない。
【００５３】
さらに、供給されたプログラムコードがコンピュータの機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに格納された後、そのプログラムコードの指示に基づいてその機能拡張ボードや機能拡張ユニットに備わるＣＰＵ等が実際の処理の一部又は全部を行い、その処理によって上述した実施の形態の機能が実現される場合にも、かかるプログラムコードが本発明の実施の態様に含まれることはいうまでもない。
【００５４】
なお、上記実施の形態において示した各部の形状及び構造は、何れも本発明を実施するにあたっての具体化のほんの一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその精神、又はその主要な特徴から逸脱することなく、様々な形で実施することができる。
【００５５】
【発明の効果】
以上説明したように、本発明によれば、上述の目的を達成することができる。
【図面の簡単な説明】
【図１】本発明の一実施形態に係る画像処理装置を利用した放射線画像処理システムの全体構成を示す概念図である。
【図２】本発明の一実施形態に係る階調変換処理の一例を概念的に示した図である。
【図３】本発明の一実施形態に適用可能な階調曲線の一例を示したグラフである。
【図４】本発明の一実施形態に適用可能な階調曲線の一例を示したグラフである。
【図５】本発明の一実施形態に適用可能な階調曲線の一例を示したグラフである。
【図６】本発明の一実施形態に係る画像処理部の構成を示した図である。
【図７】本発明の一実施形態に係る画像処理方法を示したフローチャートである。
【符号の説明】
１ Ｘ線センサ
２ 患者
３ Ｘ線源
４ Ｘ線発生装置制御部
５ 画像入力部
６ Ｘ線撮影システム制御部
７ 画像入力部
８ 画像保存部
９ 診断モニタ
１０ 操作部
１１ ネットワーク
１２ プリンタ
１３ 診断ワークステーション
１４ 画像データベース
１５ 階調変換前ヒストグラム
１６ 階調曲線
１７ 階調変換後ヒストグラム
【発明の属する技術分野】
本発明は、画像の階調変換処理を行うための画像処理装置、画像処理方法、プログラム及び記録媒体に関するものである。
【０００２】
【従来の技術】
画像全体のコントラストを強調するための代表的な手法としてヒストグラムイコライゼーションが知られている。これは階調変換を行なった後のヒストグラムの分布が平坦（均一）になるように階調変換関数を決める方法である（詳細なアルゴリズムは「ディジタル画像処理：Ｒｏｓｅｎｆｅｌｄ，Ｋａｋ著、長尾真監訳」参照）。ヒストグラムイコライゼーションを利用した階調変換自動化装置が特許第２８０８７７３号公報により知られている。尚、階調変換関数は階調変換曲線とも称され、実際の階調変換には階調変換関数に従ったルックアップテーブルが用いられる場合が多い。
【０００３】
また、Ｘ線画像も従来はフィルム画像が主であったが、近年では大画面エリアセンサを用いたディジタルＸ線撮影システムや、輝尽性蛍光体を用いたＣＲシステムなどのディジタル化が進んできている。このＸ線ディジタル画像においては、医師は従来のフィルム画像で診断してきているためフィルムと同等の画像を望む傾向がある。このためフィルムの特性曲線を階調曲線に用いて階調変換処理を行なうことによってフィルムと同等の画像を得ていた。
【０００４】
【発明が解決しようとする課題】
しかしながら、ヒストグラムイコライゼーションによる階調変換処理においては、変換後の画像のヒストグラムが平坦になることによってのみ階調曲線が決められる。したがって、その階調曲線の滑らかさは全く保障されていない。このため、Ｘ線画像の場合には従来のフィルム画像とかけ離れた画像となる場合も多く、従来のフィルム画像で診断してきた医師にとって診断しづらい画像となってしまうという問題がある。
【０００５】
本発明は、上記問題点に鑑みてなされたものであり、特定の階調変換関数を用いつつ、変換後画像の画素値ヒストグラムの平坦度が良好な階調変換を行うことのできる画像処理装置、画像処理方法、プログラム及び記録媒体を提供することを目的とする。
【０００６】
【課題を解決するための手段】
斯かる目的を達成するための本発明の第１の実施態様は、対象画像の階調変換処理を行う画像処理装置であって、階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを特徴とする画像処理装置である。
【０００７】
第２の実施態様は、対象画像の階調変換処理を行う画像処理装置に用いられる画像処理方法であって、階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを特徴とする画像処理方法である。
【０００８】
第３の実施態様は、対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムであって、前記所定の画像処理方法は、階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを含むことを特徴とする特徴とするプログラムである。
【０００９】
第４の実施態様は、対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムを記憶したコンピュータ読み取り可能な記録媒体であって、前記所定の画像処理方法は、階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを含むことを特徴とするコンピュータ読み取り可能な記録媒体である。
【００１０】
本発明の他の目的、特徴若しくは優位性等は、添付図面を参照してなされた後述の説明から明らかにされるが、以下に代表的な本発明の他の実施態様をいくつか示す。
第５の実施態様は、対象画像の階調変換処理を行う画像処理装置であって、階調変換特性を特定するための変数を用いて表わされる階調変換関数により対象画像を階調変換することにより、該変数を用いて表わされる画素値からなる画像を得る画像変換手段と、前記画像変換手段により得られた画像の画素値ヒストグラムを作成するヒストグラム作成手段と、前記ヒストグラム作成手段により得られた画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、前記階調変換関数を設定する階調変換関数設定手段と、前記階調変換関数設定手段により設定された階調変換関数により前記対象画像を階調変換する階調変換手段とを有することを特徴とする前記第１の実施態様に従った画像処理装置である。
【００１１】
第６の実施態様は、前記階調変換がルックアップテーブルを用いて行われることを特徴とする前記第１又は第５の実施態様に従った画像処理装置である。
【００１２】
第７の実施態様は、所定の評価基準に基づき前記画素値ヒストグラムの平坦度が最大となるように前記変数値が決定されることを特徴とする前記第１又は第５の実施態様に従った画像処理装置である。
【００１３】
第８の実施態様は、前記所定の評価基準は前記画素値ヒストグラムと平坦な画素値ヒストグラムとの平均二乗誤差であり、前記変数値は前記平均二乗誤差が最小となるように決定されることを特徴とする前記第７の実施態様に従った画像処理装置である。
【００１４】
第９の実施態様は、対象画像の階調変換処理を行う画像処理装置に用いられる画像処理方法であって、階調変換特性を特定するための変数を用いて表わされる階調変換関数により対象画像を階調変換することにより、該変数を用いて表わされる画素値からなる画像を得る画像変換工程と、前記画像変換工程において得られた画像の画素値ヒストグラムを作成するヒストグラム作成工程と、前記ヒストグラム作成工程において得られた画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、前記階調変換関数を設定する階調変換関数設定工程と、前記階調変換関数設定工程において設定された階調変換関数により前記対象画像を階調変換する階調変換工程とを含むことを特徴とする第２の実施態様に従った画像処理方法である。
【００１５】
第１０の実施態様は、前記階調変換がルックアップテーブルを用いて行われることを特徴とする前記第２又は第９の実施態様に従った画像処理方法である。
【００１６】
第１１の実施態様は、所定の評価基準に基づき前記画素値ヒストグラムの平坦度が最大となるように前記変数値が決定されることを特徴とする前記第２又は第９の実施態様に従った画像処理方法である。
【００１７】
第１２の実施態様は、前記所定の評価基準は前記画素値ヒストグラムと平坦な画素値ヒストグラムとの平均二乗誤差であり、前記変数値は前記平均二乗誤差が最小となるように決定されることを特徴とする前記第１１の実施態様に従った画像処理装置である。
【００１８】
第１３の実施態様は、対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムであって、前記所定の画像処理方法は、階調変換特性を特定するための変数を用いて表わされる階調変換関数により対象画像を階調変換することにより、該変数を用いて表わされる画素値からなる画像を得る画像変換工程と、前記画像変換工程において得られた画像の画素値ヒストグラムを作成するヒストグラム作成工程と、前記ヒストグラム作成工程において得られた画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、前記階調変換関数を設定する階調変換関数設定工程と、前記階調変換関数設定工程において設定された階調変換関数により前記対象画像を階調変換する階調変換工程とを含むことを特徴とする前記第３の実施態様に従ったプログラムである。
【００１９】
第１４の実施態様は、前記階調変換がルックアップテーブルを用いて行われることを特徴とする前記第３又は第１３の実施態様に従ったプログラムである。
【００２０】
第１５の実施態様は、所定の評価基準に基づき前記画素値ヒストグラムの平坦度が最大となるように前記変数値が決定されることを特徴とする前記第３又は第１３の実施態様に従ったプログラムである。
【００２１】
第１６の実施態様は、前記所定の評価基準は前記画素値ヒストグラムと平坦な画素値ヒストグラムとの平均二乗誤差であり、前記変数値は前記平均二乗誤差が最小となるように決定されることを特徴とする前記第１５の実施態様に従ったプログラムである。
【００２２】
第１７の実施態様は、対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムを記憶したコンピュータ読み取り可能な記録媒体であって、前記所定の画像処理方法は、階調変換特性を特定するための変数を用いて表わされる階調変換関数により対象画像を階調変換することにより、該変数を用いて表わされる画素値からなる画像を得る画像変換工程と、前記画像変換工程において得られた画像の画素値ヒストグラムを作成するヒストグラム作成工程と、前記ヒストグラム作成工程において得られた画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、前記階調変換関数を設定する階調変換関数設定工程と、前記階調変換関数設定工程において設定された階調変換関数により前記対象画像を階調変換する階調変換工程とを含むことを特徴とする前記第４の実施態様に従ったコンピュータ読み取り可能な記録媒体である。
【００２３】
第１８の実施態様は、前記階調変換がルックアップテーブルを用いて行われることを特徴とする前記第４又は第１７の実施態様に従ったコンピュータ読み取り可能な記録媒体である。
【００２４】
第１９の実施態様は、所定の評価基準に基づき前記画素値ヒストグラムの平坦度が最大となるように前記変数値が決定されることを特徴とする前記第４又は第１７の実施態様に従ったコンピュータ読み取り可能な記録媒体である。
【００２５】
第２０の実施態様は、前記所定の評価基準は前記画素値ヒストグラムと平坦な画素値ヒストグラムとの平均二乗誤差であり、前記変数値は前記平均二乗誤差が最小となるように決定されることを特徴とする前記第１９の実施態様に従ったコンピュータ読み取り可能な記録媒体である。
【００２６】
【発明の実施の形態】
以下、本発明の実施形態を添付図面を参照しながら詳細に説明する。
図１は本発明の一実施形態に係る画像処理装置を利用した放射線画像処理システムの全体構成を示す概略図である。
【００２７】
本実施形態に係る画像処理装置は、図１の一点鎖線内に示されるように、画像入力部５、Ｘ線撮影システム制御部６、画像処理部７及び画像保存部８で構成される。Ｘ線発生装置制御部４により制御されたＸ線源３より発生したＸ線は患者２を透過してＸ線センサ１により検知される。検知されたＸ線はディジタルＸ線画像として画像入力部５に入力される。入力されたディジタルＸ線画像は画像処理部７によってＸ線センサの補正処理、階調処理、空間フィルタリング処理等の画像処理がなされる。画像処理のなされたディジタルＸ線画像は診断モニタ９に表示されたり画像保存部８に保存されたりネットワーク１１を介してプリンタ１２、診断ワークステーション１３、画像データベース１４に出力されたりする。表示、出力された画像が満足のいくものでなかった場合には画像処理パラメータを変えるなどしながら画像処理、表示を繰り返し行う。以上の操作は操作部１０によって行われる。
【００２８】
図２は、このようなシステムで動作する階調変換処理の一例を概念的に示した図である。
ヒストグラムＨ（ｘ）１５は階調変換処理をする前の画像のヒストグラムである。このヒストグラムＨ（ｘ）１５を持つ画像に階調曲線ｆ（ｘ）１６による階調変換処理をかけて、変換された画像のヒストグラムがヒストグラムＨ’（ｘ）１７である。このヒストグラムＨ’（ｘ）１７が全く平坦（均一）である場合には、各画素値をとる画素の数は、画像全体の画素数を当該画像がとり得る画素値の数で割った平均画素数となる。本実施の形態では、このヒストグラムが全く均一である場合（平坦な画素値ヒストグラム）とＨ’（ｘ）１７との平均二乗誤差を求め、これを階調変換後画像の画素値ヒストグラムの平坦度の評価基準として用いる。尚、ここで、平坦度の評価基準は上述のような平均二乗誤差に限られず、階調変換後画像の画素値ヒストグラムと平坦な画素値ヒストグラムとのずれ量の尺度となるものであればよい。尚、このずれ量が小さいほど、平坦度が大きいものとする。
【００２９】
つまり、フィルムの特性曲線と同等になるような制約条件をかけながら、この平均二乗誤差の評価基準が最小となるような階調曲線ｆ（ｘ）１６を求めるものである。これを具体的に説明していく。
【００３０】
まず、ディジタル画像の画素値をχｉとする。画素値のとりうるレベル数をＭとすればχ１≦χｉ≦χＭである。例えば画素値として１２ビット持っている画像の場合にはＭ＝４０９６であり、χ１＝０，χＭ＝４０９５となる。１０ビットであれば、Ｍ＝１０２４であり、χ１＝０，χＭ＝１０２３となる。このディジタル画像のヒストグラムをＨ（χｉ）とする。
【００３１】
また、階調曲線（階調変換関数）をｆ（χｉ：ａ１，ａ２，・・・，ａＮ）とする。ここでｆ（χｉ：ａ１，ａ２，・・・，ａＮ）はＮ個のパラメータ（変数）ａ１，ａ２，・・・，ａＮによって規定される関数であるとする。即ち、階調変換関数ｆの階調変換特性は当該Ｎ個の変数の値により特定される。フィルムの特性曲線と同等の階調曲線の関数の一例として、シグモイド関数を用いた関数を次の式１に示す。
【００３２】
【数１】
この関数で（ａ１，ａ２，ａ３，ａ４，ａＮ）＝（１，５，３，０．２，１／１０２．４）とすれば、１０ビット画像を想定した図３のグラフになる。また、（ａ１，ａ２，ａ３，ａ４，ａＮ）＝（２，５，３，０．２，１／１０２．４）とすれば、図３のグラフのコントラストを倍にした図４のグラフの関数となる。
【００３４】
ただこの関数ではコントラストが左右対称のＳ字型曲線しかできないが、次の式２に示すような関数によって左右のコントラストの違うＳ字型曲線を実現することができる。
【００３５】
【数２】
この関数で（ａ１，ａ２，ａ３，ａ４，ａ５，ａ６，ａ７，ａ８，ａ９）＝（１，５，２，７，１／１０２．４，１／１０２．４，１，１．５，０．２）としたグラフを図５に示す。このような曲線は画素値の高い肺野のコントラストを大きくとって、画素値の低い縦隔でコントラストを小さくするために胸部用に用いられる。尚、ここで縦軸はフィルムの特性曲線に合わせて濃度とした。
【００３７】
このように階調曲線の関数を定義してパラメータをいろいろ変えれば、いろいろな種類のフィルムに対応した階調曲線をつくることができる。このような階調曲線ｆ（χ：ａ１，ａ２，・・・，ａＮ）によって変換した後のディジタル画像（各画素値はａ１，ａ２，・・・，ａＮの関数）のヒストグラムをＨ’（χｉ）とする。このＨ’（χｉ）はＨ’（ｆ（χｉ：ａ１，ａ２，・・・，ａＮ））と表わすことができる。
【００３８】
このＨ’（χｉ）と全画素値にわたって均一なヒストグラムＨ’’（χｉ）＝ｍ（ｍは平均画素数）との２乗誤差の和Ｉｅを求める。このＩｅは次の式３に示すような階調曲線関数のパラメータの関数となる。
【００３９】
【数３】
この２乗誤差の和Ｉｅが最小になるようにパラメータを決定する。このパラメータを制限のないパラメータとして扱ってもフィルムの特性曲線との同等性が保たれるのであれば、このＩｅ（ａ１，ａ２，・・・，ａＮ）を各々のパラメータａ１，ａ２，・・・，ａＮで偏微分して０として次の式４のような連立方程式を立てて、これを解けばパラメータａ１，ａ２，・・・，ａＮを求めることができる。
【００４１】
【数４】
この連立方程式を解いて求めたパラメータａ１，ａ２，・・・，ａＮを用いた階調曲線ｆ（χ：ａ１，ａ２，・・・，ａＮ）はフィルムとの同等性を保ちながら、変換後のヒストグラムの均一性を最大限にする階調曲線となっている。次に、上記連立方程式を解いて求めたパラメータａ１，ａ２，・・・，ａＮを用いた階調曲線ｆ（χ：ａ１，ａ２，・・・，ａＮ））に基づいて、画像入力部５に入力された画像情報の階調変換を行い、診断モニタ９等の情報出力手段にて出力する。
【００４３】
またパラメータに制限を与えないとフィルムの特性曲線との同等性が保たれない場合であるが、この場合で簡単な例として、特定のパラメータのみ制限のないパラメータとして扱い、その他のパラメータは定数で固定することが挙げられる。これは例えば式１でａ１，ａ２のみを変数とし、その他ａ３，ａ４，ａ５を固定する場合である。このようにすると図３、図４のようにＳ字形状を保ちながらコントラスト、横方向のシフトを行なうことができる。
【００４４】
これを用いてＩｅを求めればＩｅはａ１，ａ２の関数となり、上述のようにａ１，ａ２で偏微分して連立方程式を立てて解けばａ１，ａ２を求めることができる。さらに複雑なパラメータの制約条件がいろいろ考えられるが、そのような制約条件のもとで最適化したい場合には、例えば「関数解析による最適理論：Ｄａｖｉｄ Ｇ． Ｌｕｅｎｂｅｒｇｅｒ著、増淵正美、嘉納秀明共訳、コロナ社」ではラグランジュ乗数を用いた方法や不等式拘束条件の場合の方法など種々の方法が記載されており、これを使ってパラメータの最適化を行なえば良い。
【００４５】
尚、以上で説明した階調変換処理は、例えば、画像処理部７を図６のように構成することにより実現できる。図６において、６０は画像処理部７全体の制御を司るＣＰＵ、６１はＣＰＵ６０の動作に必要なプログラム、データ、対象画像等を記憶するための記憶手段としてのＲＡＭ、６２は階調変換関数ｆ（χ：ａ１，ａ２，・・・，ａＮ）により対象画像を階調変換する画像変換部、６３は画像変換部６２により得られた、変数で表わされる画像の画素値ヒストグラムＨ’（χｉ）を作成するヒストグラム作成部、６４はヒストグラム作成部６３により得られた画素値ヒストグラムＨ’（χｉ）と平坦ヒストグラムＨ’‘（χｉ）＝ｍとの２乗誤差を最小化する変数ａ１，ａ２，・・・，ａＮの値を求め、階調変換ルックアップテーブル（ＬＵＴ）を作成する階調変換ＬＵＴ作成部（階調変換関数設定部）、６５は階調変換ＬＵＴ作成部６４により得られた階調変換ＬＵＴを用いて対象画像を階調変換する階調変換部である。６６は上述の各部６０〜６５を通信可能に接続するＣＰＵバスである。
【００４６】
また、６２〜６５の各部は、例えば集積回路等を用いて構成することができる。また、６２〜６５の各部の機能は、ＣＰＵ６０がＲＡＭ６１等に格納された所定のプログラムを実行することによっても実現することができる。この場合、当該プログラムは、例えば図７のフローチャートに示される画像処理方法に対応したプログラムコードを含むものとすればよい。
【００４７】
次に、図７を参照して当該画像処理方法を説明する。まず、対象画像を階調変換関数ｆ（χ：ａ１，ａ２，・・・，ａＮ）により階調変換し、所定の変数で表わされた画素値からなる画像を得る（画像変換工程Ｓ７０）。次いで、前工程において得られた、変数で表わされる画像の画素値ヒストグラムＨ’（χｉ）を作成する（ヒストグラム作成工程Ｓ７１）。続いて、前工程において得られた画素値ヒストグラムＨ’（χｉ）と平坦ヒストグラムＨ’‘（χｉ）＝ｍとの２乗誤差を最小化する変数ａ１，ａ２，・・・，ａＮの値を求め、階調変換ＬＵＴを作成する（階調変換ＬＵＴ作成工程Ｓ７２）。前工程において得られた階調変換ＬＵＴを用いて対象画像を階調変換する（階調変換工程Ｓ７３）。
【００４８】
以上のように、本実施形態によれば、フィルムと同等な画質になるように変換することができる上に、変換後のヒストグラムがなるべく均一になるように効果的にコントラストの割り当てられた階調変換を実現できるため、この階調変換により得られた画像を診断医に提供することによって診断能を向上することができる。
【００４９】
各種デバイスと接続された装置或いはシステム内のコンピュータに対し、ソフトウェアのプログラムコードを供給し、そのシステム或いは装置のコンピュータ（ＣＰＵ或いはＭＰＵ）が、供給されたプログラムに従って上記各種デバイスを動作させることによって上記実施の形態の機能を実現したものも、本発明の範疇に含まれる。
【００５０】
また、この場合、上記ソフトウェアのプログラムコード自体が上述した実施の形態の機能を実現することになり、そのプログラムコード自体は本発明を構成する。そのプログラムコードの供給媒体としては、プログラム情報を所定の信号として伝搬させて供給するためのコンピュータネットワーク（ＬＡＮ、インターネット等のＷＡＮ、無線通信ネットワーク等）システムにおける通信媒体（光ファイバ等の有線回線や無線回線等）を用いることができる。
【００５１】
さらに、上記プログラムコードをコンピュータに供給するための手段、例えばかかるプログラムコードを格納した記憶媒体は本発明を構成する。かかるプログラムコードを記憶する記憶媒体としては、例えばフレキシブルディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、ＲＯＭ等を用いることができる。
【００５２】
また、供給されたプログラムコードをコンピュータが実行することにより、上述の実施の形態の機能が実現されるだけでなく、そのプログラムコードがコンピュータにおいて稼働しているＯＳ（オペレーティングシステム）或いは他のアプリケーションソフト等と共同して上述の実施の形態の機能が実現される場合にも、かかるプログラムコードが本発明の実施の態様に含まれることはいうまでもない。
【００５３】
さらに、供給されたプログラムコードがコンピュータの機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに格納された後、そのプログラムコードの指示に基づいてその機能拡張ボードや機能拡張ユニットに備わるＣＰＵ等が実際の処理の一部又は全部を行い、その処理によって上述した実施の形態の機能が実現される場合にも、かかるプログラムコードが本発明の実施の態様に含まれることはいうまでもない。
【００５４】
なお、上記実施の形態において示した各部の形状及び構造は、何れも本発明を実施するにあたっての具体化のほんの一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその精神、又はその主要な特徴から逸脱することなく、様々な形で実施することができる。
【００５５】
【発明の効果】
以上説明したように、本発明によれば、上述の目的を達成することができる。
【図面の簡単な説明】
【図１】本発明の一実施形態に係る画像処理装置を利用した放射線画像処理システムの全体構成を示す概念図である。
【図２】本発明の一実施形態に係る階調変換処理の一例を概念的に示した図である。
【図３】本発明の一実施形態に適用可能な階調曲線の一例を示したグラフである。
【図４】本発明の一実施形態に適用可能な階調曲線の一例を示したグラフである。
【図５】本発明の一実施形態に適用可能な階調曲線の一例を示したグラフである。
【図６】本発明の一実施形態に係る画像処理部の構成を示した図である。
【図７】本発明の一実施形態に係る画像処理方法を示したフローチャートである。
【符号の説明】
１ Ｘ線センサ
２ 患者
３ Ｘ線源
４ Ｘ線発生装置制御部
５ 画像入力部
６ Ｘ線撮影システム制御部
７ 画像入力部
８ 画像保存部
９ 診断モニタ
１０ 操作部
１１ ネットワーク
１２ プリンタ
１３ 診断ワークステーション
１４ 画像データベース
１５ 階調変換前ヒストグラム
１６ 階調曲線
１７ 階調変換後ヒストグラム
Claims (4)
- 対象画像の階調変換処理を行う画像処理装置であって、
階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを特徴とする画像処理装置。 - 対象画像の階調変換処理を行う画像処理装置に用いられる画像処理方法であって、
階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを特徴とする画像処理方法。 - 対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムであって、
前記所定の画像処理方法は、
階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを含むことを特徴とする特徴とするプログラム。 - 対象画像の階調変換処理を行う所定の画像処理方法をコンピュータに実行させるためのプログラムを記憶したコンピュータ読み取り可能な記録媒体であって、
前記所定の画像処理方法は、
階調変換特性を特定するための変数を用いて表わされる階調変換関数を使用し、該階調変換関数により対象画像を変換して得られる画像の画素値ヒストグラムの平坦度に基づいて前記変数の値を決定し、決定された該変数値により定まる階調変換関数により前記対象画像を階調変換することを含むことを特徴とするコンピュータ読み取り可能な記録媒体。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002305024A JP4323770B2 (ja) | 2002-10-18 | 2002-10-18 | 画像処理装置、画像処理方法、プログラム及び記録媒体 |
| US10/673,292 US7433515B2 (en) | 2002-10-18 | 2003-09-30 | Image processing apparatus and method for performing gradation conversion processing for an object image |
| EP03256300A EP1411472A3 (en) | 2002-10-18 | 2003-10-07 | Gradation conversion of an image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002305024A JP4323770B2 (ja) | 2002-10-18 | 2002-10-18 | 画像処理装置、画像処理方法、プログラム及び記録媒体 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2004139464A true JP2004139464A (ja) | 2004-05-13 |
| JP2004139464A5 JP2004139464A5 (ja) | 2005-12-02 |
| JP4323770B2 JP4323770B2 (ja) | 2009-09-02 |
Family
ID=32040870
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002305024A Expired - Fee Related JP4323770B2 (ja) | 2002-10-18 | 2002-10-18 | 画像処理装置、画像処理方法、プログラム及び記録媒体 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US7433515B2 (ja) |
| EP (1) | EP1411472A3 (ja) |
| JP (1) | JP4323770B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017168101A (ja) * | 2016-03-16 | 2017-09-21 | トムソン ライセンシングＴｈｏｍｓｏｎ Ｌｉｃｅｎｓｉｎｇ | 拡張された高ダイナミックレンジ（“ｈｄｒ”）−ｈｄｒトーンマッピングのための方法、装置、及びシステム |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4803982B2 (ja) * | 2004-09-14 | 2011-10-26 | オリンパス株式会社 | 画像処理装置、画像記録装置、及び画像処理方法 |
| US8285041B2 (en) * | 2004-09-14 | 2012-10-09 | Olympus Corporation | Image processing apparatus, image recording apparatus, and image processing method |
| RU2407221C1 (ru) * | 2006-10-10 | 2010-12-20 | Ниппон Телеграф Энд Телефон Корпорейшн | Способ и устройство управления кодированием с внутренним предсказанием, программа для них и запоминающий носитель, который хранит программу |
| JP5665393B2 (ja) * | 2010-07-05 | 2015-02-04 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
| CN103324823A (zh) * | 2012-03-05 | 2013-09-25 | 株式会社东芝 | 阅片报告制作支援装置 |
| IT202000027669A1 (it) | 2020-11-18 | 2022-05-18 | Paolo Gobbi Frattini S R L | “insieme ad elementi connettibili e disconnettibili per la ricostituzione di farmaci e nutrimenti fluidi somministrabili a pazienti con sostanze attive in polvere o gel.” |
Family Cites Families (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5983149A (ja) * | 1982-11-04 | 1984-05-14 | Fuji Photo Film Co Ltd | 階調補正曲線の作成方法 |
| JP2808773B2 (ja) | 1990-01-09 | 1998-10-08 | 株式会社日立製作所 | 階調変換自動化装置 |
| US5046118A (en) * | 1990-02-06 | 1991-09-03 | Eastman Kodak Company | Tone-scale generation method and apparatus for digital x-ray images |
| US5224177A (en) * | 1991-10-31 | 1993-06-29 | The University Of Chicago | High quality film image correction and duplication method and system |
| US5333065A (en) | 1991-12-09 | 1994-07-26 | Agfa-Gevaert | Signal to density mapping with controlled contrast |
| IL106691A (en) * | 1993-08-13 | 1998-02-08 | Sophis View Tech Ltd | System and method for diagnosis of living tissue diseases |
| US5581370A (en) | 1995-06-05 | 1996-12-03 | Xerox Corporation | Image-dependent automatic area of interest enhancement |
| WO1998037738A1 (en) | 1997-02-21 | 1998-08-27 | Direct Radiography Corp. | Image data processing for digital x-ray detectors |
| US6449390B1 (en) | 1997-09-24 | 2002-09-10 | Canon Kabushiki Kaisha | Image processing apparatus and method therefor |
| US6714674B1 (en) | 1998-04-30 | 2004-03-30 | General Electric Company | Method for converting digital image pixel values |
| JP4328399B2 (ja) * | 1998-08-31 | 2009-09-09 | キヤノン株式会社 | 画像処理装置、画像処理方法、及び記憶媒体 |
| US6370265B1 (en) | 1999-11-24 | 2002-04-09 | Direct Radiography Corp. | Method for generating gray scale transfer functions for use in displaying a digital radiogram |
| JP2001243464A (ja) | 2000-02-29 | 2001-09-07 | Canon Inc | 画像処理装置、画像処理システム、画像処理方法、及び記憶媒体 |
| JP2001299733A (ja) * | 2000-04-27 | 2001-10-30 | Konica Corp | Ｐｃｉ放射線画像処理装置、ｐｃｉ放射線画像検出処理装置、ｐｃｉ放射線画像出力装置及びｐｃｉ画像診断支援装置 |
| US6574307B1 (en) * | 2001-11-21 | 2003-06-03 | Ge Medical Systems Global Technology Company Llc | Method and apparatus for enhancing the contrast of a medical diagnostic image that includes foreign objects |
| US6614877B2 (en) * | 2001-11-21 | 2003-09-02 | Ge Medical Systems Global Technology Company Llc | Method and apparatus for enhancing the contrast of a medical diagnostic image acquired using collimation |
-
2002
- 2002-10-18 JP JP2002305024A patent/JP4323770B2/ja not_active Expired - Fee Related
-
2003
- 2003-09-30 US US10/673,292 patent/US7433515B2/en not_active Expired - Fee Related
- 2003-10-07 EP EP03256300A patent/EP1411472A3/en not_active Withdrawn
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017168101A (ja) * | 2016-03-16 | 2017-09-21 | トムソン ライセンシングＴｈｏｍｓｏｎ Ｌｉｃｅｎｓｉｎｇ | 拡張された高ダイナミックレンジ（“ｈｄｒ”）−ｈｄｒトーンマッピングのための方法、装置、及びシステム |
| JP7179443B2 (ja) | 2016-03-16 | 2022-11-29 | インターデジタル マディソン パテント ホールディングス， エスアーエス | 拡張された高ダイナミックレンジ（“ｈｄｒ”）－ｈｄｒトーンマッピングのための方法、装置、及びシステム |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4323770B2 (ja) | 2009-09-02 |
| EP1411472A3 (en) | 2005-02-16 |
| US20040091139A1 (en) | 2004-05-13 |
| EP1411472A2 (en) | 2004-04-21 |
| US7433515B2 (en) | 2008-10-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5665903B2 (ja) | 画像処理装置及び方法、画像処理システム、並びにプログラム | |
| WO2020070834A1 (ja) | 学習済みモデルの製造方法、輝度調整方法および画像処理装置 | |
| JP2004041694A (ja) | 画像生成装置およびプログラム、画像選択装置、画像出力装置、画像提供サービスシステム | |
| US11406340B2 (en) | Method for converting tone of chest X-ray image, storage medium, image tone conversion apparatus, server apparatus, and conversion method | |
| JP2002010139A (ja) | デジタル画像のトーンスケールパラメータを自動決定する方法 | |
| JP2005109790A (ja) | 医用画像処理装置 | |
| JPWO2007013300A1 (ja) | 異常陰影候補検出方法及び異常陰影候補検出装置 | |
| JP2008520346A (ja) | Ｘ線撮影画像処理向けビジュアルプレファランス用語の変換方法 | |
| US7680352B2 (en) | Processing method, image processing system and computer program | |
| JP4323770B2 (ja) | 画像処理装置、画像処理方法、プログラム及び記録媒体 | |
| JP2002017685A (ja) | 医用画像表示方法および医用画像表示装置 | |
| JP4770839B2 (ja) | 画像処理方法、画像処理装置及び画像処理プログラム | |
| JP2004030596A (ja) | 画像階調変換方法、画像階調変換装置、システム、プログラム及び記憶媒体 | |
| JP2001344601A (ja) | 画像処理装置及び画像処理プログラム | |
| US10776906B2 (en) | Radiographic image processing device, storage medium, and radiographic image processing method | |
| JP2004147084A (ja) | 画像処理装置、画像処理方法、プログラム、および記憶媒体 | |
| JP2018149166A (ja) | 放射線画像処理装置 | |
| JP2006026396A (ja) | 画像処理装置及び方法、制御プログラム並びに記憶媒体 | |
| CN113936068A (zh) | 伪影校正方法、装置和存储介质 | |
| JP2005020337A (ja) | 画像処理方法および装置並びにプログラム | |
| CN102802530B (zh) | 图像处理装置、x线图像诊断装置及图像处理方法 | |
| JP2002330948A (ja) | 画像処理制御方法および装置 | |
| JP2005027168A (ja) | 画像処理装置及び方法 | |
| JP2006167187A (ja) | 医用画像表示装置 | |
| JP2002086705A (ja) | 画像形成装置及び階調特性評価方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20051018 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051018 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080715 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080916 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081021 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20081222 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20090224 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090326 |
|
| A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20090508 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20090602 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090605 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120612 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120612 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130612 Year of fee payment: 4 |
|
| LAPS | Cancellation because of no payment of annual fees |