JPH0833722B2 - 画像表示方法 - Google Patents
画像表示方法Info
- Publication number
- JPH0833722B2 JPH0833722B2 JP1186819A JP18681989A JPH0833722B2 JP H0833722 B2 JPH0833722 B2 JP H0833722B2 JP 1186819 A JP1186819 A JP 1186819A JP 18681989 A JP18681989 A JP 18681989A JP H0833722 B2 JPH0833722 B2 JP H0833722B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- value
- procedure
- window
- image display
- 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.)
- Expired - Lifetime
Links
- 238000000034 method Methods 0.000 title claims description 60
- 238000011156 evaluation Methods 0.000 claims description 19
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 8
- 241000282412 Homo Species 0.000 description 6
- 230000003247 decreasing effect Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 238000005520 cutting process Methods 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004313 glare Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Controls And Circuits For Display Device (AREA)
Description
【発明の詳細な説明】 ［発明の目的］ （産業上の利用分野） 本発明は、ディジタル画像データを画像表示するため
の画像表示方法に関し、特に表示ウィンドウ（ウィンド
ウレベル，ウィンドウ幅）を自動設定できる画像表示方
法に関する。
の画像表示方法に関し、特に表示ウィンドウ（ウィンド
ウレベル，ウィンドウ幅）を自動設定できる画像表示方
法に関する。
（従来の技術） 従来より例えば医用画像を表示する画像表示装置にお
いては、画像メモリにディジタルで記録された画像デー
タを読み出し、この各画像データを階調変換して表示装
置、例えばモニタに表示するようにしていた。また表示
画像の細かい濃度変化を知るため、術者が表示された画
像を見ながら表示画像ごとにスイッチを操作して、ウィ
ンドウレベル，ウィンドウ幅（以下ウィンドウワイズと
いう。）の値を調整しこれらの値を最適な値に設定して
いた。
いては、画像メモリにディジタルで記録された画像デー
タを読み出し、この各画像データを階調変換して表示装
置、例えばモニタに表示するようにしていた。また表示
画像の細かい濃度変化を知るため、術者が表示された画
像を見ながら表示画像ごとにスイッチを操作して、ウィ
ンドウレベル，ウィンドウ幅（以下ウィンドウワイズと
いう。）の値を調整しこれらの値を最適な値に設定して
いた。
また自動的に設定する方法としては、撮影条件、例え
ばMRI装置における撮影条件ごとに予め前記ウィンドウ
レベル，ウィンドウワイズをプリセットしておき、この
ウィンドウレベル，ウィンドウワイズを読み出して前記
撮影条件に合わせてウィンドウレベル，ウィンドウワイ
ズをセットする方法もある。
ばMRI装置における撮影条件ごとに予め前記ウィンドウ
レベル，ウィンドウワイズをプリセットしておき、この
ウィンドウレベル，ウィンドウワイズを読み出して前記
撮影条件に合わせてウィンドウレベル，ウィンドウワイ
ズをセットする方法もある。
さらには、前記ウィンドウを自動的に設定して階調変
換する技術に関して、例えば画像表示装置（特開昭63−
84526）が既に知られている。
換する技術に関して、例えば画像表示装置（特開昭63−
84526）が既に知られている。
この画像表示装置は、画像表示メモリへ記録する画像
データから抽出した画素値に対するヒストグラムを第７
図（ａ）に示すように作成し、このヒストグラムの最大
値Hmaxに対して定められたスレッシュホールドレベルTh
により前記ヒストグラム上の最小値Wl及び最大値Whを求
め、これにより第７図（ｂ）に示すように画像表示のウ
ィンドウＷとレベルＬを定めて階調変換テーブルデータ
を作成するようにしたものである。
データから抽出した画素値に対するヒストグラムを第７
図（ａ）に示すように作成し、このヒストグラムの最大
値Hmaxに対して定められたスレッシュホールドレベルTh
により前記ヒストグラム上の最小値Wl及び最大値Whを求
め、これにより第７図（ｂ）に示すように画像表示のウ
ィンドウＷとレベルＬを定めて階調変換テーブルデータ
を作成するようにしたものである。
（発明が解決しようとする課題） しかしながら、上記従来の画像表示装置のウィンドウ
の自動設定にあっては、次のような問題がある。すなわ
ち公知技術では、単にヒストグラムのスレッシュホール
ドレベルThのみでウィンドウ・レベルＬ及びウィンドウ
・ワイドＷを設定しているが、人間は明るさ、ぎらつき
などの複数の項目により判断しているため、このウィン
ドウ・レベルＬ及びウィンドウ・ワイドＷは人間が感じ
る画像の見易さに対応するウィンドウ・レベルと一致し
ない場合が多々あった。このため観察に適したウィンド
ウが設定できない場合には、術者はこのウィンドウを観
察に適合したウィンドウにいちいち再設定しなければな
らず、術者の操作負担になっていた。
の自動設定にあっては、次のような問題がある。すなわ
ち公知技術では、単にヒストグラムのスレッシュホール
ドレベルThのみでウィンドウ・レベルＬ及びウィンドウ
・ワイドＷを設定しているが、人間は明るさ、ぎらつき
などの複数の項目により判断しているため、このウィン
ドウ・レベルＬ及びウィンドウ・ワイドＷは人間が感じ
る画像の見易さに対応するウィンドウ・レベルと一致し
ない場合が多々あった。このため観察に適したウィンド
ウが設定できない場合には、術者はこのウィンドウを観
察に適合したウィンドウにいちいち再設定しなければな
らず、術者の操作負担になっていた。
そこで本発明の目的は、見易い表示画像のウィンドウ
を自動的に設定して術者の操作負担を軽減できる画像表
示方法を提供することにある。
を自動的に設定して術者の操作負担を軽減できる画像表
示方法を提供することにある。
［発明の構成］ （課題を解決する為の手段） 本発明の上記の課題を解決し目的を達成する為に次の
ような手段を講じた。本発明は、ディジタル化された画
像データを用いてウィンドウに関するパラメータを設定
し画像を表示させるための画像表示方法において、前記
画像データに基づき画素値に対するヒストグラムを作成
する手順と、候補値として前記パラメータを複数種類設
定し、設定された前記パラメータに対してそれぞれ画像
の見易さを表現する画像度を求め、求められたそれぞれ
の画像度を比較し、最も大きい画像度をもつ前記パラメ
ータを最終設定とする手順とからなることを特徴とす
る。
ような手段を講じた。本発明は、ディジタル化された画
像データを用いてウィンドウに関するパラメータを設定
し画像を表示させるための画像表示方法において、前記
画像データに基づき画素値に対するヒストグラムを作成
する手順と、候補値として前記パラメータを複数種類設
定し、設定された前記パラメータに対してそれぞれ画像
の見易さを表現する画像度を求め、求められたそれぞれ
の画像度を比較し、最も大きい画像度をもつ前記パラメ
ータを最終設定とする手順とからなることを特徴とす
る。
また、前記パラメータを最終設定する手順は、最初に
候補値として前記パラメータを設定し、設定されたパラ
メータに対して前記画像度を求める手順と、このパラメ
ータを各々大きく変化させこれらに対する画像度を求め
る手順と、これらの画像度に基づき前記候補値のパラメ
ータを変更設定する手順と、前記パラメータの変化量を
順次小さくすることにより最も大きい画像度をもつ前記
パラメータを最終設定とする手順とを含むことを特徴と
する。
候補値として前記パラメータを設定し、設定されたパラ
メータに対して前記画像度を求める手順と、このパラメ
ータを各々大きく変化させこれらに対する画像度を求め
る手順と、これらの画像度に基づき前記候補値のパラメ
ータを変更設定する手順と、前記パラメータの変化量を
順次小さくすることにより最も大きい画像度をもつ前記
パラメータを最終設定とする手順とを含むことを特徴と
する。
また前記パラメータはウィンドウレベル，ウィンドウ
幅のうち少なくともどちらか一方であることを特徴とす
る。
幅のうち少なくともどちらか一方であることを特徴とす
る。
また画像度を決定するための定数を所定のプログラム
により決定することを特徴とする。
により決定することを特徴とする。
また画像度をルールを用いて求めることを特徴とす
る。
る。
また前記画像データは核磁気共鳴診断装置によって得
られたものであることを特徴とする。
られたものであることを特徴とする。
（作用） このような手段を講じたことにより、次のような作用
を呈する。画像データに基づき画素値に対するヒストグ
ラムを作成し、候補値としてパラメータを複数種類設定
し、設定されたパラメータに対して求められたそれぞれ
の画像度を比較して、最も大きい画像度をもつパラメー
タを最終設定とするので、最適なウィンドウ設定の操作
が簡単になる。
を呈する。画像データに基づき画素値に対するヒストグ
ラムを作成し、候補値としてパラメータを複数種類設定
し、設定されたパラメータに対して求められたそれぞれ
の画像度を比較して、最も大きい画像度をもつパラメー
タを最終設定とするので、最適なウィンドウ設定の操作
が簡単になる。
また候補として設定されたパラメータに対する画像度
と、このパラメータを各々大きく変化させたときにおけ
る各画像度とによりパラメータが変更設定され、さらに
パラメータの変化量を順次小さくすると、最大の画像度
をもつウィンドウが自動的に求められるので、最適なウ
ィンドウ設定の操作が簡単になる。これにより術者の操
作負担を軽減でき、しかも人間が感ずる見易い表示画像
を得ることができる。また定数を所定のプログラムによ
り決定すれば、術者の所望とするウィンドを設定でき
る。さらに、画像度をルールを用いて求めれば、エキス
パートシステムとして動作させることもできる。
と、このパラメータを各々大きく変化させたときにおけ
る各画像度とによりパラメータが変更設定され、さらに
パラメータの変化量を順次小さくすると、最大の画像度
をもつウィンドウが自動的に求められるので、最適なウ
ィンドウ設定の操作が簡単になる。これにより術者の操
作負担を軽減でき、しかも人間が感ずる見易い表示画像
を得ることができる。また定数を所定のプログラムによ
り決定すれば、術者の所望とするウィンドを設定でき
る。さらに、画像度をルールを用いて求めれば、エキス
パートシステムとして動作させることもできる。
（実施例） 第１図は本発明に係る画像表示方法の手順の一実施例
を示す概略図、第２図は前記画像表示方法を適用した画
像表示装置の一実施例を示す概略ブロック図である。
を示す概略図、第２図は前記画像表示方法を適用した画
像表示装置の一実施例を示す概略ブロック図である。
第１図における画像表示方法は、画像データに基づき
画素値に対するヒストグラムを作成する手順a,候補値と
して設定されたパラメータとしてのウィンドウレベル及
びウィンドウワイズに対する画像度を求める手順b,この
ウィンドウレベル及びウィンドウワイズを各々大きく変
化させこれらに対する各画像度に基づき前記候補値のウ
ィンドウレベル及びウィンドウワイズを変更設定する手
順c,前記ウィンドウレベル及びウィンドウワイズの変化
量を順次小さくすることにより最も大きい画像度をもつ
ウィンドウレベル及びウィンドウワイズを最終設定する
手順ｄとからなる。
画素値に対するヒストグラムを作成する手順a,候補値と
して設定されたパラメータとしてのウィンドウレベル及
びウィンドウワイズに対する画像度を求める手順b,この
ウィンドウレベル及びウィンドウワイズを各々大きく変
化させこれらに対する各画像度に基づき前記候補値のウ
ィンドウレベル及びウィンドウワイズを変更設定する手
順c,前記ウィンドウレベル及びウィンドウワイズの変化
量を順次小さくすることにより最も大きい画像度をもつ
ウィンドウレベル及びウィンドウワイズを最終設定する
手順ｄとからなる。
第２図において、画像表示装置は、ディジタル化され
た画像データを記録する画像メモリ1,制御手段としての
計算機2,ウィンドウレベルおよびウィンドウワイズをマ
ニュアル入力操作するためのWL/WWスイッチ3,画像メモ
リ１からの画像データを入力してWL/WWスイッチ３の操
作入力により所定のウィンドウ設定で階調を持った画像
として表示する表示装置４から構成されている。このと
きの前記画像データの画素値から階調への変換は、前記
ウィンドウレベル（以下WLという。），ウィンドウワイ
ズ（以下WWという。）により決定されるものとなってい
る。この変換方法は、例えばＸ線CT装置,MRI装置で既に
用いられている。
た画像データを記録する画像メモリ1,制御手段としての
計算機2,ウィンドウレベルおよびウィンドウワイズをマ
ニュアル入力操作するためのWL/WWスイッチ3,画像メモ
リ１からの画像データを入力してWL/WWスイッチ３の操
作入力により所定のウィンドウ設定で階調を持った画像
として表示する表示装置４から構成されている。このと
きの前記画像データの画素値から階調への変換は、前記
ウィンドウレベル（以下WLという。），ウィンドウワイ
ズ（以下WWという。）により決定されるものとなってい
る。この変換方法は、例えばＸ線CT装置,MRI装置で既に
用いられている。
前記計算機２は、前記画像メモリ１から読み出される
画像データに基づき、最適のWL,WWを表示装置４に自動
的に設定するものとなっている。なおWL/WWスイッチ３
をマニュアル操作することにより、WL,WWを最適に設定
することもできる。
画像データに基づき、最適のWL,WWを表示装置４に自動
的に設定するものとなっている。なおWL/WWスイッチ３
をマニュアル操作することにより、WL,WWを最適に設定
することもできる。
次に前記第１図及び第２図を参照して前記計算機２に
より最適のWL,WWを自動設定する画像表示方法を具体的
に説明する。ここでは、表示画像をMRI画像とする。第
３図は典型的なMRI画像のヒストグラムを示す概略図で
ある。
より最適のWL,WWを自動設定する画像表示方法を具体的
に説明する。ここでは、表示画像をMRI画像とする。第
３図は典型的なMRI画像のヒストグラムを示す概略図で
ある。
１）まず計算機２により画像メモリ１から画像データ
が読み出され、この画像データに対して第３図に示すMR
I画像のヒストグラムが作成される（第１図手順ａ）。
人体のMRI画像の内、体外に相当する部分は、一般的にM
RI値（画素値）が低くほぼ一定値である。これにより体
外に相当する部分は、第３図に示す斜線部分のような背
景Ｈとなる。この背景Ｈは観察しない不要部分である。
が読み出され、この画像データに対して第３図に示すMR
I画像のヒストグラムが作成される（第１図手順ａ）。
人体のMRI画像の内、体外に相当する部分は、一般的にM
RI値（画素値）が低くほぼ一定値である。これにより体
外に相当する部分は、第３図に示す斜線部分のような背
景Ｈとなる。この背景Ｈは観察しない不要部分である。
２）そこで前記計算機２により背景Ｈを暗くするカッ
ト処理が以下のごとく行なわれる。最初に前記背景Ｈと
体内部との境界のMRI値をもとめ、第３図に示すように
このMRI値をＡとする。すなわちMRI値Ａを求めるには、
まず最低のMRI値以上、かつこのMRI値の一定範囲内であ
って、ヒストグラムのピーク値Hmaxを求める。このピー
ク値HmaxからMRI値の大きい方に数個のヒストグラム上
の点H1，H2…を求め、その点H1，H2…を通る直線Lnを最
小二乗法により求める。この直線LnとMRI値を表すＸ軸
との交点におけるMRI値をＡを求め、Ａ以下のMRI値をカ
ット処理することにより背景Ｈを暗くする。
ト処理が以下のごとく行なわれる。最初に前記背景Ｈと
体内部との境界のMRI値をもとめ、第３図に示すように
このMRI値をＡとする。すなわちMRI値Ａを求めるには、
まず最低のMRI値以上、かつこのMRI値の一定範囲内であ
って、ヒストグラムのピーク値Hmaxを求める。このピー
ク値HmaxからMRI値の大きい方に数個のヒストグラム上
の点H1，H2…を求め、その点H1，H2…を通る直線Lnを最
小二乗法により求める。この直線LnとMRI値を表すＸ軸
との交点におけるMRI値をＡを求め、Ａ以下のMRI値をカ
ット処理することにより背景Ｈを暗くする。
次に画像度について説明する。画像度は、あるウィン
ドウの場合に、そのウィンドウにより表示された画像
が、どの程度見易いかを指数化し、これを数値として表
示したものである。この画像度は次のようにして求めら
れる。画像度をＱとし、画像度を決定する項目の数をｋ
個とすると、 Ｑは により計算される。ここでWnは項目ｎの重み係数であ
り、Fnは項目ｎの評価値である。重み係数Wnはその項目
が画像度に関連する割合により経験的に決定したもので
あり、例えばW1＝5,W2＝4,W3＝3,W4＝2,W5＝２のような
数値である。また項目ごとの評価値Vnは、項目ごとに異
なる評価関数Fnにより求められる。Wn,Fnを最初から最
適に設定することが困難な場合もあるが、その場合には
ウィンドウの設定結果を見ながら最適になるようにWn,F
nを調整する。
ドウの場合に、そのウィンドウにより表示された画像
が、どの程度見易いかを指数化し、これを数値として表
示したものである。この画像度は次のようにして求めら
れる。画像度をＱとし、画像度を決定する項目の数をｋ
個とすると、 Ｑは により計算される。ここでWnは項目ｎの重み係数であ
り、Fnは項目ｎの評価値である。重み係数Wnはその項目
が画像度に関連する割合により経験的に決定したもので
あり、例えばW1＝5,W2＝4,W3＝3,W4＝2,W5＝２のような
数値である。また項目ごとの評価値Vnは、項目ごとに異
なる評価関数Fnにより求められる。Wn,Fnを最初から最
適に設定することが困難な場合もあるが、その場合には
ウィンドウの設定結果を見ながら最適になるようにWn,F
nを調整する。
ここでは例えば項目の１つであるMVPに関する評価値
を求めるための評価関数F1について説明する。V1＝F
1（MVP）である。
を求めるための評価関数F1について説明する。V1＝F
1（MVP）である。
第４図に示すようにMRI値はWL,WWに応じて、連続した
明るさとして表示される。ただしWL−WW/2以下のMRI値
は全て最も暗く、WL＋WW/2以上のMRI値は、全て最も明
るく表示される。この明るさを便宜上第４図の縦軸に示
す数値、すなわち0.5〜16.5の連続値として表現してい
る。
明るさとして表示される。ただしWL−WW/2以下のMRI値
は全て最も暗く、WL＋WW/2以上のMRI値は、全て最も明
るく表示される。この明るさを便宜上第４図の縦軸に示
す数値、すなわち0.5〜16.5の連続値として表現してい
る。
前記第３図に示すようにMVPは、最も見たい部分と推
定されるMRI値である。このMVPは第３図に示すような単
峰性のヒストグラムを有する画像では、例えばＡ以上の
MRI値の中で最も頻度の高いMRI値をMVPと定めている。
このMVPが見易い明るさになるようにWL,WWを設定するこ
とが望ましい。見易い明るさとは経験的に真ん中の明る
さより少し明るいと考えられる。従って、ここではMVP
が“10"の明るさになるとき、最も（MVPに関する。）評
価値が高くなるように、“10"から離れるに従って評価
値が低くなるようにF1を設定する。F1は厳密に定義され
るものではなく、ファジィ理論におけるメンバーシップ
関数のように、人間の感覚により定義（設定）する。例
えばF1は第５図に示すような非線形関数に設定する。
定されるMRI値である。このMVPは第３図に示すような単
峰性のヒストグラムを有する画像では、例えばＡ以上の
MRI値の中で最も頻度の高いMRI値をMVPと定めている。
このMVPが見易い明るさになるようにWL,WWを設定するこ
とが望ましい。見易い明るさとは経験的に真ん中の明る
さより少し明るいと考えられる。従って、ここではMVP
が“10"の明るさになるとき、最も（MVPに関する。）評
価値が高くなるように、“10"から離れるに従って評価
値が低くなるようにF1を設定する。F1は厳密に定義され
るものではなく、ファジィ理論におけるメンバーシップ
関数のように、人間の感覚により定義（設定）する。例
えばF1は第５図に示すような非線形関数に設定する。
同様にF2〜FKの関数を設定することができるが、例え
ば第２の項目は“A"点の明るさに関するものであり、
“A"点以下は背景であるから、暗いほうが望ましい。従
って、“A"点の明るさが0.5であるときには、評価値が
高くなり、明るくなるに従い、評価値が下がる関数を設
定する。
ば第２の項目は“A"点の明るさに関するものであり、
“A"点以下は背景であるから、暗いほうが望ましい。従
って、“A"点の明るさが0.5であるときには、評価値が
高くなり、明るくなるに従い、評価値が下がる関数を設
定する。
以下同様に、「明るさのバランス」に関する項目な
ど、これ以外の項目に関する関数を設定する。
ど、これ以外の項目に関する関数を設定する。
３）次に前記計算機２によりWLc＝WL0,WWc＝WW0に設
定する。WL0，WW0は初期値であり、例えば第３図に示す
ようにMVPを明るさの中心として、WL0＝MVP,背景Ｈ以外
は全て見えるようにするために、WW0＝２（MVP−Ａ）の
ように設定する。
定する。WL0，WW0は初期値であり、例えば第３図に示す
ようにMVPを明るさの中心として、WL0＝MVP,背景Ｈ以外
は全て見えるようにするために、WW0＝２（MVP−Ａ）の
ように設定する。
ここで設定したWLc,WWcはウィンドの候補値であり、
以下のような手順により最適なWL,WWを求める。
以下のような手順により最適なWL,WWを求める。
まず、設定されたWL,WWにおける画像の見易さを画像
度により評価する（第１図手順ｂ）。次にWL,WWをそれ
ぞれ大きく変化させて、同様にその時の画像の見易さを
評価し、その中で最も高い画像度を持つWL,WWに設定す
る（手順ｃ）。そして順次WL,WWの変化量を少なくして
いき、最終的に見易いWL,WWに設定する（手順ｄ）。
度により評価する（第１図手順ｂ）。次にWL,WWをそれ
ぞれ大きく変化させて、同様にその時の画像の見易さを
評価し、その中で最も高い画像度を持つWL,WWに設定す
る（手順ｃ）。そして順次WL,WWの変化量を少なくして
いき、最終的に見易いWL,WWに設定する（手順ｄ）。
この方法の具体的な手順を第６図を参照して説明す
る。第６図（ａ）はWW,WLを変化させる場合を示す概略
図、第６図（ｂ）はWW,WLに対する画像度Ｑを示す概略
図である。
る。第６図（ａ）はWW,WLを変化させる場合を示す概略
図、第６図（ｂ）はWW,WLに対する画像度Ｑを示す概略
図である。
ΔWW＝WW0/2とする。
WLc,WWcにおける画像度Qcを計算する（前記手順ｂ
に対応）。この画像度は前述した方法により求められ
る。
に対応）。この画像度は前述した方法により求められ
る。
ａ）次にWLc＋ΔWW,WWcにおける画像度Qd1を前記
と同様な方法で計算する。
と同様な方法で計算する。
ｂ）WLc＋ΔWW,WWc＋ΔWWにおける画像度Qd2を計算
する。
する。
ｃ）WLc,WWc＋ΔWWにおける画像度Qd3を計算する。
ｄ）WLc−ΔWW,WWc＋ΔWWにおける画像度Qd4を計算
する。
する。
ｅ）WLc−ΔWW,WWcにおける画像度Qd5を計算する。
ｆ）WLc−ΔWW,WWc−ΔWWにおける画像度Qd6を計算
する。
する。
ｇ）WLc,WWc−ΔWWにおける画像度Qd7を計算する。
ｈ）WLc＋ΔWW,WWc−ΔWWにおける画像度Qd8を計算
する。
する。
そして画像度Qd1〜Qd8の中で、最大の画像度をQdを
求め、またその画像度をQdにおけるWL,WWをWLd,WWdとす
る。
求め、またその画像度をQdにおけるWL,WWをWLd,WWdとす
る。
第５図（ｂ）に示すように画像度Qc＜Qdならば、Qcか
らこれより高い画像度Qdに移動する。すなわち、WLc＝W
Ld,WWc＝WWdとし、に戻る（手順ｃに対応）。
らこれより高い画像度Qdに移動する。すなわち、WLc＝W
Ld,WWc＝WWdとし、に戻る（手順ｃに対応）。
Qc≧Qdならば、ΔWW＝α×ΔWWとする（手順ｄに対
応）。ここでα＜１なる定数である。
応）。ここでα＜１なる定数である。
ΔWW≧ｒならば、へ戻る。ここでｒは終了を決め
る定数である。
る定数である。
ΔWW＜ｒならば、WLc,WWcが最終的に求めるWL,WWとな
る。
る。
このように本実施例によれば、候補として設定された
WL,WWに対する画像度と、WL,WWを各々大きく変化させた
ときにおける各画像度とによりWL,WWが変更設定され、
さらにWL,WWの変化量を順次小さくすると、最大の画像
度が自動的に求められるので、最適なウィンドウ設定の
操作が簡単になる。これにより術者の操作負担を軽減で
き、しかも人間が感ずる見易い表示画像を得ることがで
きる。
WL,WWに対する画像度と、WL,WWを各々大きく変化させた
ときにおける各画像度とによりWL,WWが変更設定され、
さらにWL,WWの変化量を順次小さくすると、最大の画像
度が自動的に求められるので、最適なウィンドウ設定の
操作が簡単になる。これにより術者の操作負担を軽減で
き、しかも人間が感ずる見易い表示画像を得ることがで
きる。
さらに、ここではWn,Fnは最初から最適に設定するこ
とが困難であるため、ウィンドウの設定結果を見ながら
最適になるように、人間がWn,Fnを調整したが、これを
プログラムで自動的に設定する方法もある。例えば、簡
単のためにWnだけを自動的に設定する場合について説明
する。ここで項目数を５個、すなわちWnをW1〜W5の場合
である。例えばWnは１〜５の整数値とする。この場合W1
〜W5の組み合わせは3125組になる。
とが困難であるため、ウィンドウの設定結果を見ながら
最適になるように、人間がWn,Fnを調整したが、これを
プログラムで自動的に設定する方法もある。例えば、簡
単のためにWnだけを自動的に設定する場合について説明
する。ここで項目数を５個、すなわちWnをW1〜W5の場合
である。例えばWnは１〜５の整数値とする。この場合W1
〜W5の組み合わせは3125組になる。
また、基準となる画像を例えば10枚選択する。この画
像に対して専門家である放射線技師が設定したウィンド
ウをプログラムに入力する。プログラムはW1〜W5の3125
組のそれぞれの場合について、10枚の画像について最適
なWL,WWを求め放射線技師が設定したウィンドと比較
し、最も差の少ないW1〜W5を見付け出す。こうすれば、
専門家である放射線技師に最も近いWL,WWを設定できるW
1〜W5を決定できる。同様に、Fnに関する定数につてい
も最適の値が決定できる。またこうすれば、画像表示装
置ごとの特性も定数の中に組み込むことができる。定数
は１度設定すればほとんど変える必要がないので、計算
時間はそれ程問題にならない。
像に対して専門家である放射線技師が設定したウィンド
ウをプログラムに入力する。プログラムはW1〜W5の3125
組のそれぞれの場合について、10枚の画像について最適
なWL,WWを求め放射線技師が設定したウィンドと比較
し、最も差の少ないW1〜W5を見付け出す。こうすれば、
専門家である放射線技師に最も近いWL,WWを設定できるW
1〜W5を決定できる。同様に、Fnに関する定数につてい
も最適の値が決定できる。またこうすれば、画像表示装
置ごとの特性も定数の中に組み込むことができる。定数
は１度設定すればほとんど変える必要がないので、計算
時間はそれ程問題にならない。
また以上の方式は手続き型のプログラムで実現するこ
ともできるが、ルールを用いてエキスパートシステムと
して実現することもできる。画像の見易さは人間の感覚
に依存するところが大きいのでエキスパートシステムと
することも有効である、 なお、上述した実施例では、MRI画像について説明し
たが、画像はディジタル画像であれば、MRI画像に限定
されることなく、例えばＸ線CT画像，核医学等の画像で
あっても良い。また前述したヒストグラムの形状が異な
る場合には、ヒストグラムの形状に応じて場合分けを行
ない、適用項目を選択する方法にすればよい。また２個
以上の峰を持つヒストグラムの場合には、さらに複雑な
方法によりMVPを求める。どの様な方法によりMVPを求め
ても、本発明を適用することができる。本発明の要旨を
逸脱しない範囲で種々変形実施可能であるのは勿論であ
る。
ともできるが、ルールを用いてエキスパートシステムと
して実現することもできる。画像の見易さは人間の感覚
に依存するところが大きいのでエキスパートシステムと
することも有効である、 なお、上述した実施例では、MRI画像について説明し
たが、画像はディジタル画像であれば、MRI画像に限定
されることなく、例えばＸ線CT画像，核医学等の画像で
あっても良い。また前述したヒストグラムの形状が異な
る場合には、ヒストグラムの形状に応じて場合分けを行
ない、適用項目を選択する方法にすればよい。また２個
以上の峰を持つヒストグラムの場合には、さらに複雑な
方法によりMVPを求める。どの様な方法によりMVPを求め
ても、本発明を適用することができる。本発明の要旨を
逸脱しない範囲で種々変形実施可能であるのは勿論であ
る。
［発明の効果］ 本発明は画像データに基づき画素値に対するヒストグ
ラムを作成し、候補値としてパラメータを複数種類設定
し、設定されたパラメータに対して求められたそれぞれ
の画像度を比較して、最も大きい画像度をもつパラメー
タを最終設定とするので、最適なウィンドウ設定の操作
が簡単になる。
ラムを作成し、候補値としてパラメータを複数種類設定
し、設定されたパラメータに対して求められたそれぞれ
の画像度を比較して、最も大きい画像度をもつパラメー
タを最終設定とするので、最適なウィンドウ設定の操作
が簡単になる。
また候補として設定されたパラメータに対する画像度
と、このパラメータを各々大きく変化させたときにおけ
る各画像度とによりパラメータが変更設定され、さらに
パラメータの変化量を順次小さくすると、最大の画像度
をもつウィンドウが自動的に求められるので、最適なウ
ィンドウ設定の操作が簡単になる。これにより術者の操
作負担を軽減でき、しかも人間が感ずる見易い表示画像
を得ることができる。また定数を所定のプログラムによ
り決定すれば、術者の所望とするウィンドを設定でき
る。さらに、画像度をルールを用いて求めれば、エキス
パートシステムとして動作させることもできる画像表示
方法を提供できる。
と、このパラメータを各々大きく変化させたときにおけ
る各画像度とによりパラメータが変更設定され、さらに
パラメータの変化量を順次小さくすると、最大の画像度
をもつウィンドウが自動的に求められるので、最適なウ
ィンドウ設定の操作が簡単になる。これにより術者の操
作負担を軽減でき、しかも人間が感ずる見易い表示画像
を得ることができる。また定数を所定のプログラムによ
り決定すれば、術者の所望とするウィンドを設定でき
る。さらに、画像度をルールを用いて求めれば、エキス
パートシステムとして動作させることもできる画像表示
方法を提供できる。
第１図は本発明に係る画像表示方法の手順の一実施例を
示す概略図、第２図は画像表示方法を適用した画像表示
装置の一実施例を示す概略ブロック図、第３図は典型的
なMRI画像のヒストグラムを示す概略図、第４図はMRI値
をWL,WWに応じて変換した明るさを示す概略図、第５図
は評価関数を示す概略図、第６図は画像度を示す概略
図、第７図は従来の画像表示装置における画像表示方法
の一例を示す概略図である。 １……画像メモリ、２……計算機、３……WL/WWスイッ
チ、４……表示装置、Ｑ……画像度、WL,L……ウィンド
ウレベル、WW,W……ウィンドウワイド、Ｈ……背景、MV
P……MRI値の最大値、Ｆ……評価関数。
示す概略図、第２図は画像表示方法を適用した画像表示
装置の一実施例を示す概略ブロック図、第３図は典型的
なMRI画像のヒストグラムを示す概略図、第４図はMRI値
をWL,WWに応じて変換した明るさを示す概略図、第５図
は評価関数を示す概略図、第６図は画像度を示す概略
図、第７図は従来の画像表示装置における画像表示方法
の一例を示す概略図である。 １……画像メモリ、２……計算機、３……WL/WWスイッ
チ、４……表示装置、Ｑ……画像度、WL,L……ウィンド
ウレベル、WW,W……ウィンドウワイド、Ｈ……背景、MV
P……MRI値の最大値、Ｆ……評価関数。
Claims (6)
- 【請求項１】画像データをウインドウ幅及びウインドウ
レベルに基づいて輝度変換し、表示するための画像表示
方法において、 前記輝度変換された画像データの中の最高頻度を示す画
素値に割り当てられる輝度に基づいて第１の評価値を求
める第１の手順と、 前記輝度変換された画像データの中の背景画素を除く最
低画素値に割り当てられる輝度に基づいて第２の評価値
を求める第２の手順と、 前記第１の評価値と前記第２の評価値とを加算すること
により画像値を求める第３の手順と、 前記ウインドウ幅と前記ウインドウレベルの少なくとも
一方を変化させて、前記画像値が最大になるウインドウ
幅とウインドウレベルを求める第４の手順と、 前記第４の手順で求められたウインドウ幅とウインドウ
レベルに基づいて前記画像データを輝度変換する第５の
手順と、 前記第５の手順で輝度変換された画像データを表示する
第６の手順とからなり、 前記第１の評価値は、前記最高頻度を示す画素値が中間
輝度より明るく表示されるように、前記最高頻度を示す
画素値に割り当てられる輝度が前記中間輝度より明るい
所定輝度に近いほど高い値を示すことを特徴とする画像
表示方法。 - 【請求項２】前記第４の手順は、前記画像データの画素
値毎の頻度を表す第１のヒストグラムに基づいて、前記
ウインドウ幅と前記ウインドウレベル各々の初期値を求
めるサブ手順を含むことを特徴とする請求項１記載の画
像表示方法。 - 【請求項３】前記サブ手順は、前記第１のヒストグラム
から背景領域を除去した第２のヒストグラムを作成し、
前記第２のヒストグラムの最高頻度に対応するピーク画
素値MVPを求め、前記第２のヒストグラムの最低画素値
Ａを求め、前記ピーク画素値MVPを前記ウインドウレベ
ルの初期値とし、前記最低画素値Ａから前記ピーク画素
値MVPまでの２倍を前記ウインドウ幅の初期値とするこ
とを特徴とする請求項２記載の画像表示方法。 - 【請求項４】前記第２の評価値は前記最低画素値に割り
当てられる輝度が低いほど高くなり、前記割り当てられ
る輝度が高いほど低くなる値を示すことを特徴とする請
求項１記載の画像表示方法。 - 【請求項５】前記第４の手順は、前記第１の手順乃至前
記第３の手順により前記画像値が求められる毎に、前記
ウインドウ幅と前記ウインドウレベルの変化量を減少さ
せていくことを特徴とする請求項１記載の画像表示方
法。 - 【請求項６】前記第３の手順は前記第１の評価値と前記
第２の評価値とを重ね付け加算することにより画像値を
求めることを特徴とする請求項１記載の画像表示方法。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1186819A JPH0833722B2 (ja) | 1989-07-19 | 1989-07-19 | 画像表示方法 |
| US07/552,520 US5305204A (en) | 1989-07-19 | 1990-07-16 | Digital image display apparatus with automatic window level and window width adjustment |
| DE69031523T DE69031523T2 (de) | 1989-07-19 | 1990-07-18 | Digitales Bildanzeigegerät |
| EP90113780A EP0409206B1 (en) | 1989-07-19 | 1990-07-18 | Digital image display apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1186819A JPH0833722B2 (ja) | 1989-07-19 | 1989-07-19 | 画像表示方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0351891A JPH0351891A (ja) | 1991-03-06 |
| JPH0833722B2 true JPH0833722B2 (ja) | 1996-03-29 |
Family
ID=16195156
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1186819A Expired - Lifetime JPH0833722B2 (ja) | 1989-07-19 | 1989-07-19 | 画像表示方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0833722B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE602004028099D1 (de) | 2003-02-07 | 2010-08-26 | Panasonic Corp | Beleuchtungseinrichtung, einen sockel verwendend, um ein flaches led-modul auf einen kühlkörper zu montieren |
| JP5530626B2 (ja) * | 2008-12-19 | 2014-06-25 | 株式会社日立メディコ | 磁気共鳴イメージング装置、及び画像処理装置 |
| JP5513749B2 (ja) * | 2009-01-30 | 2014-06-04 | 株式会社東芝 | Ｘ線撮影装置及びｘ線画像処理方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61194576A (ja) * | 1985-02-23 | 1986-08-28 | Toshiba Corp | 画像処理装置 |
| JP2721152B2 (ja) * | 1986-06-17 | 1998-03-04 | 株式会社東芝 | 画像データのウインドウ設定装置 |
| JPS6384526A (ja) * | 1986-09-30 | 1988-04-15 | 株式会社 日立メデイコ | 画像表示装置 |
| FI880818A (fi) * | 1987-03-18 | 1988-09-19 | Gen Electric | Bildskaerm med automatisk bildjustering. |
-
1989
- 1989-07-19 JP JP1186819A patent/JPH0833722B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0351891A (ja) | 1991-03-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| DE69629445T2 (de) | Automatische Tonskalenabstimmung mittels Bildaktivitätsmessungen | |
| EP0409206B1 (en) | Digital image display apparatus | |
| US5150421A (en) | System for automated transformation of gray level of image | |
| JP4130025B2 (ja) | ２次データ・セットの情報をマージする方法及びシステム、強調された画像データを再生する方法及びシステム、強調されたディジタル・データを再生するシステム、１次データ・セット・ワードを強調する方法、並びに強調されたディジタル・データを表示する方法 | |
| US7315638B2 (en) | Image display system | |
| KR100327745B1 (ko) | 엑스선 컴퓨터 단층 촬영 장치 및 영상 진단 장치 | |
| JP5039310B2 (ja) | 脳内出血部位セグメンテーション装置 | |
| US20020183607A1 (en) | Method and system for visualizing a body volume and computer program product | |
| US20040170247A1 (en) | Displaying image data using automatic presets | |
| US5218623A (en) | Method and apparatus for specifying slice planes in x-ray computed tomography | |
| US5748801A (en) | Method of setting threshold values for extracting image data | |
| DE10239646A1 (de) | Interaktive Echtzeit-Segmentierung von pulmonalen Knoten mit Steuerparametern | |
| JPH07129751A (ja) | 医用画像処理装置 | |
| JPS6311146A (ja) | 画像診断装置 | |
| US7079140B2 (en) | Diagnostic device having means for setting transfer functions | |
| EP1741058B1 (en) | Method and system for automatically improving the usability of a medical picture | |
| JP5042533B2 (ja) | 医用画像表示装置 | |
| JPH0833722B2 (ja) | 画像表示方法 | |
| US9454811B2 (en) | Image processing apparatus, image processing method and computer-readable storage medium storing program | |
| JP2937389B2 (ja) | 画像表示方法 | |
| JPH043969B2 (ja) | ||
| JP3748457B2 (ja) | ウインドウ操作装置 | |
| US7945112B2 (en) | Method for displaying digital grey scale images at desired tonal value on a screen of a display device and a display device | |
| DE10224756B4 (de) | Verfahren und Vorrichtung zur Erzeugung einer kombinierten Parameterkarte | |
| JP2937390B2 (ja) | 画像変換方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313114 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R360 | Written notification for declining of transfer of rights |
Free format text: JAPANESE INTERMEDIATE CODE: R360 |
|
| R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313114 |
|
| R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313114 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080329 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090329 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100329 Year of fee payment: 14 |
|
| EXPY | Cancellation because of completion of term | ||
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100329 Year of fee payment: 14 |