JP2003126057A

JP2003126057A - 画像処理装置

Info

Publication number: JP2003126057A
Application number: JP2001322060A
Authority: JP
Inventors: Ken Ishikawa; 謙石川
Original assignee: Hitachi Medical Corp
Current assignee: Hitachi Healthcare Manufacturing Ltd
Priority date: 2001-10-19
Filing date: 2001-10-19
Publication date: 2003-05-07
Anticipated expiration: 2021-10-19
Also published as: JP3967105B2

Abstract

(57)【要約】【課題】ツブレの影響を除去したダイナミックレンジ
圧縮を行える画像処理装置を提供する。【解決手段】上記課題は、各種物理量の局所的な計測
値である画素値の二次元分布を読み取って電気信号に変
換する画像読取手段１と、該変換された電気信号を処理
して原画像を生成し、その生成された原画像の画素値分
布、その画素値分布を平滑化して得られる平滑化画素値
分布、および前記画素値分布の局所的な変動成分の各種
パラメータを求めるともに、該求めた平滑化画素値分布
をパラメータとして与える第１の単調増加関数の演算結
果と前記各種パラメータの組み合わせを関連づけて与え
る第２の単調増加関数の演算結果とを加算して加算画素
値分布を演算する信号処理手段２と、該演算された加算
画素値分布を画像として出力する画像表示装置３とを備
えたことによって解決される。

Description

【発明の詳細な説明】【０００１】【発明が属する技術分野】本発明は、医用ディジタル画
像の階調特性を変更する画像処理装置に係り、特にコン
トラスト圧縮処理に好適な画像処理装置に関する。【０００２】【従来の技術】従来、画像情報の大局的な階調特性を変
更して、通常では可視不良域にあたる濃度域を可視域に
移動させ画像を見易くする画像処理方法として、例えば
Medical Imaging Technology誌、第16巻、第6号(1998年
11月発行)の第655頁から第665頁に掲載された、笠井、
藤田他による研究論文「腫瘤陰影自動検出アルゴリズム
における左右乳房画像の比較による偽陽性候補の削除」
の中で論じられているようなダイナミックレンジ圧縮処
理がよく用いられてきた。【０００３】この種のダイナミックレンジ圧縮処理で
は、原画像の画素値分布Ｏから該画素値分布Ｏをボケマ
スク法などで平滑化した画素値分布Ｂを求め、Ｂを濃度
変換したｋ(Ｂ)を原画像Ｏから差し引くことで、濃度域
の移動を行う。この処理は画像処理後の画素値分布をＰ
としたときに式(１)と表される。Ｐ＝Ｏ−ｋ(Ｂ)…(１) 式(１)を変形すれば、式(２)となる。Ｐ＝ｈ(Ｂ)＋Ｖ…(２) ここで、ｈ(Ｂ)は、式(３)に示されるように、画像情報
の大局的な階調特性の変換を表す関数である。ｈ(Ｂ)＝Ｂ−ｋ(Ｂ)…(３) 式(３)で示される階調特性は、単調増加している途中か
ら減少させるようなことはないので、ｈ(Ｂ)はＢの広義
の単調増加関数といえる。なお、白黒反転とダイナミッ
クレンジ圧縮処理と一緒に行うこともあるが、その場合
は演算式全体を符号反転した処理となる。また、式(２)
でＶは式(４)で示され、Ｏの局所的な変動成分に相当し
ている。Ｖ＝Ｏ−Ｂ…(４) 【０００４】こうした従来のダイナミックレンジ圧縮処
理で、比較的小さくて大きなコントラストを有する対象
を含む画像を処理すると、場合によっては処理後画像Ｐ
の画素値変動範囲が原画像Ｏの変動範囲を大きく逸脱し
てしまうことが生じる。そして、処理後画像を表示する
フィルムやＣＲＴの濃度や輝度の変動範囲は有限なの
で、画像情報の一部が、表示媒体の濃度や輝度の変動範
囲を超えたために失われてしまう現象(以下「ツブレ」
と称する)も生じるものであった。【０００５】【発明が解決しようとする課題】しかしながら、上記変
動範囲には比較的重要な診断情報が含まれることが多
く、上記ツブレによって貴重な診断情報が見えなくなる
おそれがあるという問題があった。【０００６】本発明の目的は、ツブレの影響を除去した
ダイナミックレンジ圧縮を行える画像処理装置を提供す
ることにある。【０００７】【課題を解決するための手段】上記目的は、各種物理量
の局所的な計測値である画素値の二次元分布を読み取っ
て電気信号に変換する画像読取手段と、該変換された電
気信号を処理して原画像を生成し、その生成された原画
像の画素値分布、その画素値分布を平滑化して得られる
平滑化画素値分布、および前記画素値分布の局所的な変
動成分の各種パラメータを求めるともに、該求めた平滑
化画素値分布をパラメータとして与える第１の単調増加
関数の演算結果と前記各種パラメータの組み合わせを関
連づけて与える第２の単調増加関数の演算結果とを加算
して加算画素値分布を演算する信号処理手段と、該演算
された加算画素値分布を画像として出力する画像表示装
置とを備えたことを特徴とする画像処理装置によって達
成される。【０００８】具体的には、各種物理量の局所的な計測値
である画素値の二次元分布を読み取って電気信号に変換
する画像読取手段と、この画像読取手段からの電気信号
を処理する信号処理手段と、この信号処理手段からの出
力信号を画像として出力する画像出力手段とから成る画
像処理装置において、上記信号処理手段が、上記原画像
の画素値分布Ｏから、該画素値分布Ｏを平滑化して得ら
れる画素値分布Ｂと、上記画素値分布Ｏの局所的な変動
成分Ｖとを求め、画像出力時の画素値分布をＰとしたと
きに演算式Ｐ＝ｈ(Ｂ)＋ｆ(Ｂ;Ｖ) (ただし、ｈ(Ｂ)はＢの単調増加関数である。又、ｆ
(Ｂ;Ｖ)はＢをパラメータとするＶの関数であり、Ｂが
一定であるとき、Ｖ₁＜Ｖ₂なるＶの組に対して常にｆ
(Ｂ;Ｖ₁)≦ｆ(Ｂ;Ｖ₂)となるＶに関する広義の単調増加
関数である。また、ｆ(Ｂ;０)＝０であり、Ｖのとり得
る最大値ＶmaxとＯのとり得る最大値Ｏmaxに対してｆ
(Ｂ;Ｖmax)≦Ｏmax−ｈ(Ｂ)となり、Ｖのとり得る最小
値ＶminとＯのとり得る最小値Ｏminに対してｆ(Ｂ;Ｖmi
n)≧Ｏmin−ｈ(Ｂ)となる性質を有する。)で表される演
算を行うようにしたものである。【０００９】また、関数ｆ(Ｂ;Ｖ)が、ｇ(0)=0 なる性質を持つ関数ｇ(ａ)とＶが正もしくは０であると
き、Ｏのとり得る最大値Omaxと定数αによりγ＝α(Ｏm
ax−ｈ(Ｂ)) で与えられ、Ｖが負であるとき、Ｏのとり得る最小値Ｏ
minと上式と同じ定数αによりγ＝α(ｈ(Ｂ)−Ｏmin) で与えられるパラメータγと原画像の画素値分布Ｏの関
数もしくはＯを平滑化して得られる計測値分布Ｂの関数
もしくは定数であるパラメータβとにより、演算式ｆ
(Ｂ;Ｖ)＝γｇ(βＶ/γ)によって規定されるものであ
る。つまり、上記組み合せはパラメータＢとＶで規定さ
れる。【００１０】【発明の実施の形態】まず、ダイナミックレンジ圧縮処
理によるツブレの現象について具体例を用いて明らかに
する。【００１１】例えば、図２(１)に断面が示されるよう
に、対象に図面上方からＸ線を照射して撮像する場合を
考える。この対象は、辺縁部分を肉厚に中央部分を肉薄
に形成された板状であって、辺縁部に近い部分には小さ
な穴が設けられ、中央部分にはＸ線吸収率が高い鉛のよ
うなＸ線高吸収体が置かれたものを使用する。ここでは
説明を簡単にするため、上記対象の一断面のＸ線透過量
分布を表わす一次元プロフィールで表わす。一次元プロ
フィール(図２(２))では、対象の辺縁部分を低い画素値
として対応づけ、中央部分は高い画素値として対応づけ
るが、前記穴の箇所は殆どＸ線吸収が起こらずに最大値
に近い高い画素値となり、前記Ｘ線高吸収体の箇所は殆
どＸ線吸収がされてしまい最小値に近い低い画素値とな
って、これらは小さい画像領域内で局所的に大きなコン
トラストを有する部分となっている。なお、対象の断面
(図２(１))と比較して一次元プロフィール(図２(２))の
変化が緩やかなのは散乱Ｘ線の影響を受けているためで
ある。【００１２】図２(２)の一次元プロフィールにダイナミ
ックレンジ圧縮処理を行うと、図２(３)の実線で示すプ
ロフィールに変化する。点線で示す処理前のプロフィー
ルと比較すると低画素値の肉厚の部分は画素値が上が
り、高画素値の肉薄の部分は画素値が下がり、原画像の
低空間周波数成分の階調特性が圧縮される方向に変化し
ていることがわかる。【００１３】しかし、その信号成分に重畳される高周波
数の変動成分の大きさは変わらないので、肉厚部に開け
られた穴の部分の画素値も押し上げられて、最大画素値
を超えツブレとなっている。又、肉薄部に置かれた金属
小片の部分の画素値も引き下げられ、こちらは最小画素
値を超えたツブレとなっている。【００１４】本発明による画像処理方法は、Ｘ線量、光
学濃度、輝度、Ｘ線吸収係数、水素原子濃度、放射性同
位元素濃度、超音波反射量、温度等の各種物理量の局所
的な計測値である画素値の二次元分布である原画像に対
して、画像情報を読み取って電気信号に変換しこの電気
信号を画像として出力するに際し、画像情報の大局的な
階調特性を変更して、通常では可視不良域にあたる濃度
域を可視域に移動させ画像を見易くするとともに、処理
に伴うツブレの発生を防止するものである。【００１５】このような画像処理を行うために、本発明
では、二次元画像上の位置座標をx,yとし、原画像の各
種物理量の局所計測値、例えば局所Ｘ線量である原画像
の画素値分布Ｏ(x,y)から、劾分布Ｏ(x,y)を平滑化して
得られる大局的な画素値分布Ｂ(x,y)と、上記画素値分
布Ｏ(x,y)の局所的な変動成分Ｖ(x,y)とを求め、画像出
力時すなわち画像処理後の画素値分布Ｐ(x,y)と次の式
(５)のように関係づける。Ｐ(x,y)＝ｈ(Ｂ(x,y))＋ｆ(Ｂ(x,y);Ｖ(x,y)) …(５) 【００１６】ここで、ｈ(Ｂ)はＢの広義の単調増加関数
であり、ダイナミックレンジ圧縮に相当する画像情報の
大局的な階調特性の変換を表す関数である。又、ｆ(Ｂ;
Ｖ)は画素値分布Ｂ(x,y)をパラメータとする変動成分Ｖ
(x,y)の関数である。ｆ(Ｂ;Ｖ)はＢが一定であるとき、
式(６)なる値の組に対して式(７)なる関係を満たす広義
の単調増加関数である。Ｖ₁＜Ｖ₂ …(６) ｆ(Ｂ;Ｖ₁)≦ｆ(Ｂ;Ｖ₂) …(７) 【００１７】ここで、Ｖは変動成分であるので正負の値
をとるが、Ｂの値によらず式(８)である。ｆ(Ｂ;０)＝０ …(８) そして、Ｂを一定にしたとき、Ｖのとり得る最大値Ｖma
xに対して式(９)を満たし、ｆ(Ｂ;Ｖmax)≦Ｏmax−ｈ(Ｂ) …(９) Ｖのとり得る最小値Ｖminに対して式(１０)を満たすよ
うにする。ｆ(Ｂ;Ｖmin)≧Ｏmin−ｈ(Ｂ) …(１０) 【００１８】ここで、式(９)のＯmaxは原画像がとり得
る最大画素値であり、式(１０)のＯminは原画像がとり
得る最小画素値である。こうすれば、式(５)、式(９)、
式(１０)より、式(１１)となり、Ｏmin≦Ｐ(x,y)≦Ｏmax …(１１) 画像として再生する際の処理後画像の濃度分布Ｐ(x,y)
は原画像がとり得る濃度範囲〔Ｏmin,Ｏmax〕を超える
ことがないのでツブレ現象は完全に除去できる。【００１９】上記の性質を満たす関数ｆ(Ｂ;Ｖ)とし
て、例えば、次のようなものが考えられる。即ち、式
(１２)とし、γをＶが正もしくは０の場合は式(１３)と
し、Ｖが負の場合は式(１４)とする。ｆ(Ｂ;Ｖ)＝(２γ／π)arctan〔(π/2)βＶ/γ〕…(１２) γ＝Ｏmax−ｈ(Ｂ) …(１３) γ＝Ｂ−Ｏmin …(１４) 式(１０)において、βは変動成分の強さを調整するパラ
メータであり、arctanはラジアン単位の逆正接関数であ
り、πは円周率である。【００２０】次に、このように構成された画像処理方法
の動作について、図１に示す上記の画像処理方法の実施
に使用する画像処理装置に即して説明する。【００２１】図１において、画像読取手段１は、各種物
理量の局所的な計測値である画素値の二次元分布である
原画像に対して、画像情報を読み取って電気信号に変換
するもので、例えば入射Ｘ線の線量分布を光電変換によ
って画素毎の信号電荷に変換し、マトリックス状に配置
されたスイッチング素子を駆動して順次信号電荷を読出
してディジタル信号に変換するＸ線平面検出器から成
る。信号処理手段２は、上記画像読取手段１から出力さ
れたディジタル信号を入力して所要の処理、例えばダイ
ナミックレンジ圧縮を行うものである。画像出力手段
は、上記信号処理手段２からの出力信号を画像として出
力するもので、信号処理手段のディジタル信号をアナロ
グの映像信号に変換してモニタ画面に表示する画像表示
装置３や、処理後の信号に従って輝度変調したレーザ光
をフィルムもしくは印画紙上に走査して画像のハードコ
ピーを作成するレーザプリンタ４などから成る。【００２２】このような画像処理装置を用いて、Ｘ線
量、光学濃度、輝度、Ｘ線吸収係数、水素原子濃度、放
射性同位元素濃度、超音波反射量、温度等の各種物理量
の局所的な計測値である画素値の二次元分布である原画
像に対して画像処理するには、まず、画像読取手段１で
原画像の画像情報を読み取りディジタル信号に変換す
る。次に、この画像読取手段１から出力されたディジタ
ル信号は、信号処理手段２へ入力し、その内部に設けら
れた第一の画像記憶装置5aに蓄積される。この画像記憶
装置５ａ内で画像情報は画素(i,j)ごとの濃度値に対応
したディジタル値Ｏ(i,j)として蓄えられる。ここで、
(i,j)は画素の番号をマトリックス形式で表わしたもの
であり、例えば原画像の左上隅から右にｉ番目、下にｊ
番目の画素を表わす。【００２３】次に、上記画像記憶装置5aの出力信号Ｏ
(i,j)を平滑化演算器６に入力し、矩形マスクによる平
滑化を行い、出力信号Ｂ(i,j)を第二の画像記憶装置５
ｂに蓄える。即ち、平滑化演算器６では式(１５)で表わ
される演算を行う。【００２４】ただし、上式でｎはマスクサイズであり、
除算n/2の小数点以下は切り捨てられる。マスクサイズ
ｎはｎ値群７の中から第一の切換器８ａによって選択す
る。平滑化された画像信号Ｂ(i,j)は更に画像情報の大
局的な階調特性を変換するためにｈテーブル群９中の第
二の切替器８ｂで選択したｈテーブルに入力してｈ(Ｂ
(i,j))に変換し、第三の画像記憶装置5cに蓄える。ｈテ
ーブル群９にはＩＣメモリで構成されたルックアップテ
ーブルであるｈテーブルが複数個含まれており、ＩＣメ
モリのアドレス入力に対応する入力値Ｂ(i,j)に応じ
て、ＩＣメモリのデータ出力に対応する出力値ｈ(Ｂ(i,
j))が関係づけられている。【００２５】以上の準備を行った後、実際の画像処理を
実行する。即ち、まず、第一の画像記憶装置５ａから読
み出された原画像の画素値Ｏ(i,j)と、第二の画像記憶
装置５ｂから読み出された平滑化した画像濃度Ｂ(i,j)
を第一の引き算器１０ａにＶ(i,j)＝Ｏ(i,j)−Ｓ(i,j) …(１６) を行う。この引き算器１０ａの出力Ｖ(i,j)は局所的変
動成分に相当する。【００２６】同時に、第三の画像記憶装置５ｃから読み
出された変換された平滑化画素値ｈ(Ｂ(i,j))を第二の
引き算器１０ｂと第三の引き算器１０ａ１０ｃにそれ
ぞれ入力し、二種類の減算を式(１７)、式(１８)のよう
に行うＯmax−ｈ(Ｂ(i,j)) …(１７) ｈ(Ｂ(i,j))−Ｏmin …(１８) ここで、Ｏmaxは画素値の変動範囲の上限、Ｏminは下限
である。【００２７】上記第一の引き算器１０ａから出力された
信号Ｖ(i,j)の符号ビットは選択器１１に入力し、第二
または第三の引き算器１０ｂ、１０ｃの出力のいずれか
を選択して出力させる。【００２８】即ち、選択器１１により、Ｖ(i,j)の符号
ビットが０である、即ちＶ(i,j)が正もしくは０の場合
は、選択器１１の出力が式(１７)の結果となるように、
符号ビットが1である、即ちＶ(i,j)が負の場合は式(１
８)の結果となるようにしている。ここで、選択器１１
の出力は前出の式(１３)、式(１４)で規定されるγと同
じもので、処理後の変動成分の上限、下限を与えるもの
である。即ちＶ(i,j)が正もしくは０の場合は式(１９)
となり、負の場合は式(２０)となる信号γ(i,j)を得
る。 γ(i,j)＝Ｏmax−Ｂ(i,j) …(１９) γ(i,j)＝Ｂ(i,j)−Ｏmin …(２０) 【００２９】一方、式(１２)において変動成分の強さを
調整するするパラメータであるβは、第三の切換器８ｃ
によって、β値群１２の中から選ばれる。逆数テーブル
１４は、前述のｈテーブルと同様にＩＣメモリで構成さ
れたルックアップテーブルであり、入力信号であるγ
(i,j)の逆数に対応した信号が得られるようになってい
る。パラメータβに対応する信号と局所的変動成分に対
応する信号Ｖ(i,j)及び逆数テーブル１４の出力は、第
一の乗算器１３ａ及び第二の乗算器１３ｂにより互いに
乗ぜられ、その結果式(２１)に対応した信号となり、逆
正接テーブル１５に出力する。 βＶ(i,j)／γ(i,j) …(２１) 【００３０】上記逆正接テーブル15は、前述のｈテーブ
ルと同様にＩＣメモリで構成されたルックアップテーブ
ルであり、入力をａとしたとき、式(２２)で規定される
出力ｂが得られるようになっている。ｂ＝(2/π)arctan(πa/2) …(２２) 【００３１】ここで、arctanはラジアン単位の逆正接関
数であり、πは円周率である。上記逆正接テーブル１５
の出力は、第三の乗算器１３ｃにより信号γ(i,j)と乗
ぜられ、その結果は、加算器１６により信号ｈ(Ｂ(i,
j))に加えられる。そして、この加算器１６の出力が本
発明の信号処理手段２の出力となる。【００３２】以上の説明から明らかなように、加算器１
６の出力は、ｈ(Ｂ(i,j))+(2/π)γ(i,j)arctan〔(π/2)βＶ(i,j)／γ(i,j)〕…(２３) となるが、これは、式(５)の右辺に式(１２)の右辺を代
入したものに対応している。従って、本実施例の信号処
理手段２は、ダイナミックレンジ圧縮のような画像情報
の大局的な階調特性の変換を行いながら、それに伴うツ
ブレ現象を完全に除去できる。【００３３】更に、本実施例の信号処理手段２において
は、処理に用いる演算式に含まれる全ての関数に広義の
単調性と連続性及び傾きの連続性が保証されているの
で、処理画像に不自然な階調の急激な変化や、急激な一
定値化や跳びが発生することがない。本実施例の信号処
理手段２は汎用の電子計算機によっても実現できる。し
かし、電子計算機では乗算、除算、逆正接関係計算など
に時間がかかる。それに対して、本実施例では乗算器１
３ａ〜１３ｃ、逆数テーブル１４、逆正接テーブル１５
を用いているので、高速の処理が可能である。【００３４】以上、本発明の画像処理方法及びその装置
の実施例について説明してきたが、本発明は上記の実施
例に限定されることなく各種の変形が可能である。【００３５】例えば、上記の実施例で用いられた(2/π)
arctan〔(π/2)ａ〕の逆正接関数を用いた関数は、次の
性質を備えた関数ｇ(ａ)もしくはそれを折線で近似した
関数に置き換えることができる。即ち、式(２４)〜式
(２９)のの諸性質である。ｇ(0)=0…(２４) 【００３６】これらの性質を満たす逆正接関数以外の関
数としては、例えば滴定曲線に類似した式(３０)などが
挙げられる。ｇ(ａ)＝｛2/〔１＋exp(−2a)〕｝＋１ …(３０) 上式でexpは指数関数を表わす。また、上記の実施例で
は画像情報はＸ線量として与えられているが、これを光
学濃度、輝度、Ｘ線吸収係数、水素原子濃度、放射性同
位元素濃度、超音波反射量、温度等で置き換えることも
可能である。【００３７】本実施形態により、ツブレを除去したダイ
ナミックレンジ圧縮を行うための画像処理装置を提供す
ることができる。これによって、例えば医用ディジタル
画像の鮮鋭化処理において処理画像の高画質化を達成で
き、良い診断情報を得ることができる。【００３８】【発明の効果】本発明は、ツブレの影響を除去したダイ
ナミックレンジ圧縮を行える画像処理装置を提供すると
いう効果を奏する。

【図面の簡単な説明】【図１】本発明による画像処理装置の実施形態を示すブ
ロック図。【図２】図１の画像処理装置で行うダイナミックレンジ
圧縮でのツブレの発生例と、本発明のツブレ抑制の効果
を説明する説明図。【符号の説明】１…画像読取手段、２…信号処理手段、３…画像表示装
置

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.⁷ 識別記号ＦＩテーマコート゛(参考）Ｇ０６Ｔ 1/00 ２９０Ｈ０４Ｎ 5/20 ５Ｃ０２１ 5/00 １００Ａ６１Ｂ 5/05 ３８０５Ｃ０７７Ｈ０４Ｎ 1/407 Ｇ０１Ｎ 24/02 ５２０Ｙ 5/20 Ａ６１Ｂ 6/00 ３５０Ｍ 5/325 Ｈ０４Ｎ 1/40 １０１ＥＦターム(参考） 4C093 AA01 AA22 AA26 CA50 FD01 FD11 FD13 FF06 FF08 FF34 FH02 4C096 AB50 AD14 DC09 DC12 DC40 4C301 EE20 JC07 JC12 JC20 KK03 LL02 4C601 EE30 JC04 JC15 JC18 JC19 KK03 LL01 LL02 5B057 AA07 CA08 CA12 CA16 CB08 CB12 CB16 CC01 CE05 CE11 DB02 DB09 DC23 5C021 PA66 PA76 XA31 5C077 LL04 MP01 PP02 PP10 PP47 PP48 PQ12 PQ23 SS02 SS06

Claims

【特許請求の範囲】【請求項１】各種物理量の局所的な計測値である画素
値の二次元分布を読み取って電気信号に変換する画像読
取手段と、該変換された電気信号を処理して原画像を生
成し、その生成された原画像の画素値分布、その画素値
分布を平滑化して得られる平滑化画素値分布、および前
記画素値分布の局所的な変動成分の各種パラメータを求
めるともに、該求めた平滑化画素値分布をパラメータと
して与える第１の単調増加関数の演算結果と前記各種パ
ラメータの組み合わせを関連づけて与える第２の単調増
加関数の演算結果とを加算して加算画素値分布を演算す
る信号処理手段と、該演算された加算画素値分布を画像
として出力する画像表示装置とを備えたことを特徴とす
る画像処理装置。