JP2937389B2 - 画像表示方法 - Google Patents

画像表示方法

Info

Publication number
JP2937389B2
JP2937389B2 JP2066695A JP6669590A JP2937389B2 JP 2937389 B2 JP2937389 B2 JP 2937389B2 JP 2066695 A JP2066695 A JP 2066695A JP 6669590 A JP6669590 A JP 6669590A JP 2937389 B2 JP2937389 B2 JP 2937389B2
Authority
JP
Japan
Prior art keywords
image
brightness
neural network
histogram
candidate
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
Application number
JP2066695A
Other languages
English (en)
Other versions
JPH03266188A (ja
Inventor
昭南 大橋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to JP2066695A priority Critical patent/JP2937389B2/ja
Priority to US07/552,520 priority patent/US5305204A/en
Priority to DE69031523T priority patent/DE69031523T2/de
Priority to EP90113780A priority patent/EP0409206B1/en
Publication of JPH03266188A publication Critical patent/JPH03266188A/ja
Application granted granted Critical
Publication of JP2937389B2 publication Critical patent/JP2937389B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Description

【発明の詳細な説明】 [発明の目的] (産業上の利用分野) 本発明は、ディジタル化された画像データの画素値を
用いてウィンドウに関するパラメータを設定し画像を表
示させるための画像表示装置に関する。
(従来の技術) 従来より例えば医用画像を表示する場合、画像メモリ
にディジタルで記録された画像データを読み出し、この
各画像データを階調変換して表示する装置、例えばモニ
タに表示するようにしていた。この場合、画像を見易く
表示するために、術者が表示された画像を見ながら表示
画像ごとにスイッチを操作して、階調変換に関するパラ
メータとしてのウィンドウレベル(以下WLという。),
ウィンドウワイズ(以下WWという。)(以下WL、WWを併
せてウィンドウという。)の値を調整し、表示画像を見
易く表示させていた。
また自動的にウィンドウを設定する方法としては、撮
影条件、例えばMRI(マグネティック・レゾナンス・イ
メージング)装置における撮影条件ごとに、予め前記ウ
ィンドウをプリセットしておき、これを読み出して前記
撮影条件に合わせてウィンドウをセットする方法もあ
る。
さらには、前記ウィンドウを自動的に設定して階調変
換する技術に関して、例えば画像表示装置(特開昭63-8
4526)が既に知られている。
この画像表示装置は、画像表示メモリへ記録する画像
データから抽出した画素値に対するヒストグラムを第13
図(a)に示すように作成し、このヒストグラムの最大
値H maxに対して定められたスレッシュホールドレベルT
hにより前記ヒストグラム上の最小値W1及び最大値Whを
求め、これにより第13図(b)に示すように画像表示の
ウィンドウワイズWとウィンドウレベルLを定めて階調
変換テーブルデータを作成するようにしたものである。
(発明が解決しようとする課題) しかしながら、上記従来の画像表示装置のウィンドウ
の自動設定にあっては、次のような問題がある。すなわ
ち公知技術では、単にヒストグラムのスレッシュホール
ドレベルThのみでウィンドウ・レベルLおよびウィンド
ウ・ワイズWを設定しているため、人間が設定したウィ
ンドウと一致しない場合が多々あった。このため観察に
適したウィンドウが設定できない場合には、操作者はこ
のウィンドウを観察に適合したウィンドウにいちいち再
設定しなければならず、操作者の負担になっていた。
さらには、前記ウィンドウを自動的に設定して階調変
換する技術に関して、例えば既に同一出願人が出願した
特願平1-186819がある。(この発明を第1の発明とす
る。)第1の発明はファジィー類似方式におけるメンバ
ーシップ関数類似方式によりウィンドウを自動的に設定
するものである。このためメンバーシップ関数などの設
定が必要であり、また1度設定したメンバーシップ関数
を装置ごとにあるいは病院ごとの特性に合わせて変更す
ることが困難であった。
そこで本発明の目的は、画像を表示する際のウィンド
ウの設定操作を簡単化して、操作者の負担を軽減でき、
しかも装置の特性に合った画像を容易に観察し得る画像
表示装置を提供することにある。
[発明の構成] (課題を解決するための手段) 本発明によれば、ディジタル画像データの各画素値
を、所定範囲内の画素値のみが連続的に変化する明るさ
値に変換されるように、明るさ値に変換する画像表示装
置において、ディジタル画像データから画素値のヒスト
グラムを求める手段と、前記ヒストグラム中の複数の所
定のパラメータを入力し、前記ディジタル画像データを
表示する際の表示画像の見易さを示す画像度を出力する
ニューラルネットワークと、前記所定範囲を変えたとき
の前記ニューラルネットワークの出力に応じて、前記所
定範囲の最適値を求める手段とを具備する画像表示装置
が提供される。
前記所定範囲の最適値を求める手段は、前記所定範囲
の候補を決める手段と、前記候補における画像度と前記
候補をずらした場合における画像度とを求める手段と、
前記候補における画像度と候補をずらした場合における
画像度とを比較し、候補における画像度が最大となるま
で前記候補をずらしながら、前記候補における画像度と
前記候補をずらした場合における画像度とを繰返し求め
る手段とを具備するものである。
前記ニューラルネットワークは学習精度が異なる2つ
の階層型のニューラルネットワークからなり、前記最適
値を求める手段は学習精度の低いニューラルネットワー
クから順に使って最適値を求め、各ニューラルネットワ
ークの出力が次のニューラルネットワークによる処理に
おける候補となるものである。
前記ニューラルネットワークはヒストグラム中の表示
画像の明るさに関する第1パラメータ、表示画像のコン
トラストに関する第2パラメータ、および明るさのバラ
ンスに関する第3パラメータを入力することを特徴とす
るものである。
前記第1パラメータは頻度が最も高い画素値の明るさ
データとヒストグラム上で体内と体外との境界にある画
素値の明るさデータを含み、前記第2パラメータはヒス
トグラム上の表示明るさ範囲の中間の明るさ部分の面積
データと最大明るさ部分の面積データとを含み、前記第
3パラメータはヒストグラム上の前記中間明るさより明
るい部分の面積と暗い部分の面積の比データを含むもの
である。
前記ニューラルネットワークは1枚の画像内の各部分
毎の見易さを示す画像度を求め、前記頻度が最も高い画
素値の明るさデータはヒストグラム内の各ピーク点の画
素値の明るさデータを含むものである。
(作用) このような手段を講じたことにより、次のような作用
を呈する。評価値をニューラルネットワークにより求め
るので、ファジィー類似方式におけるメンバーシップ関
数等を決定することなく、また装置ごとにあるいは病院
ごとの特性に合った画像を容易に観察することができ
る。
またニューラルネットワークにより評価値を求め、最
も大きい評価値をもつウィンドウでMRI画像を観察で
き、画像を表示する際にウィンドウ設定の操作が簡単に
なり、操作者の負担を軽減できる。
(実施例) 第1図は本発明に係る画像表示方法の全体の手順を示
すフロー図、第4図は前記画像表示方法を適用した画像
表示装置を示す概略ブロック図、第5図は典型的なMRI
画像のヒストグラムを示す概略図である。
第4図において、画像表示装置は、ディジタル化され
た画像データを記録する画像メモリ1,制御手段としての
計算機2,ウィンドウレベルおよびウィンドウワイズをマ
ニュアル入力操作するためのWL/WWスイッチ3,画像メモ
リ1からの画像データを入力して設定されたウィンドウ
で階調を持った画像として表示する表示装置4、ニュー
ラルネットワーク5から構成されている。このときの前
記画像データの画素値から階調への変換は、前記WL、WW
により決定されるものとなっている。この変換方法は、
例えばMRI装置で既に用いられている。
前記計算機2は、前記画像メモリ1から読み出される
画像データに基づき、最適のWL,WWを表示装置4に自動
的に設定するものとなっている。なおWL/WWスイッチ3
をマニュアル操作することにより、WL,WWを最適に設定
することもできる。
次に前記第1図、第4図、第5図を参照して、前記計
算機2により表示ウィンドウを求める方法を具体的に説
明する。ここでは、画像をMRI画像とする。第5図は典
型的なMRI画像のヒストグラムを示す概略図である。
1)計算機2は画像記憶装置(図示していない)から指
定された画像を読み出し、画像メモリ1に記憶する。
2)計算機2は、画像メモリ1の画像についてヒストグ
ラムを作成する(手順100)。作成された前記ヒストグ
ラムの典型的な例を第5図に示す。
3)次に計算機2は手順200において、背景を除去する
ために、このヒストグラムを解析し“A"を求める。“A"
の意味、および求め方は後に詳細に説明する。
4)最も観察したいと推定されるMRI値をMVPと名付け
る。次に手順300において、ヒストグラムを解析しMVPを
求める。MVPの定義、および求め方は後に詳細に説明す
る。
5)次に手順400において、その画像を表示するのに適
しているウィンドウを求める。この手順については後に
第2図、第3図を参照して詳細に説明する。
6)次に手順500において、求められたウィンドウによ
り、画像を表示する。
7)また自動設定されたウィンドウに対して観察者が必
要に応じ、第4図のWL、WWスイッチ3を操作してウィン
ドウを微調整することもできる。
したがって、画像を表示する場合には、表示ウィンド
ウの設定の操作が簡単になり、操作者の負担を軽減でき
る。
次に手順200において“A"の意味、および求め方を第
5図を参照して詳細に説明する。
第5図は手順100において作成されたMRI画像のヒスト
グラムの典型的な例である。人体のMRI画像の内、体外
に相当する部分は、一般的にMRI画素値(画像の各画素
のMRI値)が低くほぼ一定値である。これにより体外に
相当する部分は、第5図に示す斜線部分のような背景H
となる。この背景Hは観察を必要としない部分である。
そこで前記計算機2により前記背景Hと体内部分との
境界のMRI値を以下のごとく求める。このMRI値をAとす
る。第5図に示すようにMRI値Aを求めるには、まず最
低のMRI値以上、かつ一定以下のMRI値の範囲内であっ
て、ヒストグラムのピーク値H maxを求める。このピー
ク値H maxからMRI値の大きい方に数個のヒストグラム上
の点H1,H2…を求め、その点H1,H2…を通る直線Lnを最小
二乗法により求める。この直線Lnとヒストグラムの0を
表すX軸との交点におけるMRI値を求め、これをAとす
る。A以下のMRI値の部分を背景Hとする。
次に手順300におけるMVPの定義、および求め方を第5
図を参照して詳細に説明する。
MVPは、最も見たい部分と推定されるMRI値である。例
えば、その画像の中で背景を除き、最も多くの頻度を有
するMRI値が見易いように表示すれば、全体が見易い画
像になると考えて良い。したがって、第5図に示すよう
なヒストグラムを有する画像では、例えばA以上のMRI
値の中で最も頻度の高いMRI値をMVPと定めている。
したがって、MVPは、A以上のMRI値でヒストグラムの
ピークを探すことにより求められる。ここでピーク(Pe
aK)は、第8図に示すように1点の値が両端での数点の
値より互い高い値であると定義している。
前記ピークは1つのヒストグラムに複数存在するが、
それぞれのピークの重用度を求め、最も重用度の高いピ
ークをMVPとする。ピークの重用度は、第9図に示すよ
うに 1)ピークの間の距離 2)ピークの高さ 3)ピークの広がり などを考慮して評価する。
次に手順400におけるこの画像を表示するに適したウ
ィンドウの求め方について第2図の実施例を示すフロー
図を参照して詳細に説明する。
1)まず手順bにおいて、候補値として設定されたウィ
ンドウWL及びWWで画像を表示した場合の画像の見易さを
表現する評価値としての画像度をニューラルネットワー
クにより求める。この画像度の求め方については後に詳
細に説明する。
2)次に手順cにおいて、このWL及びWWを各々大きく変
化させて、その時の画像の見易さを評価し、その中で最
も高い画像度を持つWL,WWに設定する。
3)次に手順dにおいて、前記WL及びWWの変化量を順次
小さくすることにより最も大きい画像度をもつWL及びWW
を求める。
なおウィンドウに関するパラメータは、WLまたはWWの
うち少なくともいずれか一方であっても良い。
この方法の具体的な手順を第3図および第7図を参照
して説明する。第7図(a)はWW,WLを変化させる場合
を示す概略図、第7図(b)はWW,WLに対する画像度Q
を示す概略図、第3図は具体的な手順の実施例を示すフ
ロー図である。
1)手順4210において、前記計算機2によりウィンドウ
の初期値WL0、WW0を求める。初期値は、例えば第5図に
示すようにMVPを明るさの中心として、WL0=MVP,背景H
以外は全て見えるようにするために、WW0=2(MVP-A)
のように設定する。
2)手順4220において、WLc=WL0,WWc=WW0に設定す
る。ここでWLc、WWcは候補のウィンドウである。
3)手順4230において、ΔWW=WW0/2とし、ウィンドウ
の変化量に初期値を設定する。
4)手順4240において、WLc,WWcにおける画像度Qcを計
算する(前記手順bに対応)。この画像度の求め方につ
いては後に詳細に説明する。
5)次に手順4250において、 a)WLc+ΔWW,WWcにおける画像度Qd1を前記4)と同
様な方法で計算する。
b)WLc+ΔWW,WWc+ΔWWにおける画像度Qd2を計算す
る。
c)WLc,WWc+ΔWWにおける画像度Qd3を計算する。
d)WLc−ΔWW,WWc+ΔWWにおける画像度Qd4を計算す
る。
e)WLc−ΔWW,WWcにおける画像度Qd5を計算する。
f)WLc−ΔWW,WWc−ΔWWにおける画像度Qd6を計算す
る。
g)WLc,WWc−ΔWWにおける画像度Qd7を計算する。
h)WLc+ΔWW,WWc−ΔWWにおける画像度Qd8を計算す
る。
6)次に手順4260において、画像度Qd1〜Qd8の中で、
最大の画像度をQdを求め、その画像度QdにおけるWL,WW
をWLd,WWdとする。
7)次に第7図(b)に示すように、手順4270におい
て、画像度Qc<Qdならば、Qcからこれより高い画像度Qd
に移動する。すなわち、手順4280において、WLc=WLd,W
Wc=WWdとし、手順4240へ戻る(手順cに対応)。
8)Qc≧Qdならば、手順4290において、ΔWW=α×ΔWW
とする(手順dに対応)。ここでα<1なる定数であ
る。
9)手順4295において、ΔWW≧rならば、手順4240へ戻
る。ここでrは終了を決める定数である。ΔWW<rなら
ば、WLc,WWcが最終的に求めるWL,WWとなる。
このように、候補として設定されたWL,WWに対する画
像度と、WL,WWを各々大きく変化させたときにおける各
画像度とによりWL,WWが変更設定され、さらにWL,WWの変
化量を順次小さくすると、最大の画像度が自動的に求め
られるので、最適なウィンドウが設定できる。これによ
り観察者の操作負担を軽減でき、しかも人間が感ずる見
易い表示画像を得ることができる。
次に手順4240〜4250において用いられている画像度を
求める方式について詳細に説明する。画像度はある設定
されたウィンドウで画像を表示した場合に、そのウィン
ドウにより表示された画像が、どの程度見易いかを指数
化し、これを評価値として表示したものである。
本発明ではこの画像度をニューラルネットワークを用
いて求め、装置ごとに,病院ごとの特性に合った画像表
示方法を提供するものである。以下、具体的な実施例を
説明する。
前記第4図において、計算機2は前記画像度を求める
場合には、画像度の求めるための評価項目をニューラル
ネットワーク5に転送する。ニューラルネットワーク5
は画像度を求め、その画像度を計算機2に転送する。そ
の結果、計算機2は画像度を得ることができる。
以下、画像度をニューラルネットワーク5により求め
る方法について詳細に説明する。
第10図は前記ニューラルネットワーク5の構成を示す
図である。ニューラルネットワーク5は、例えば入力層
30,中間層32,出力層34の3層で構成されている。入力層
30は、入力層30a〜30eの5素子からなり、入力層30a〜3
0eに例えば、 1)MVPの明るさ 2)“A"の明るさ 3)明るさ12の面積 4)明るさ値Lの面積 5)明るさのバランス の評価項目(特徴)を入力する。本実施例では特徴を
5個としたが、これは本実施例に限定されるものではな
い。特徴の個数が増減した場合は、それに応じて入力層
30の素子数を増減させれば良い。これらの特徴の具体的
な求め方については後に説明する。
前記中間層32は例えば50素子からなり、前記入力層30
a〜30eともネットワークを構成している。
出力層34は1素子で構成され、前記出力層34から画像
度0〜1を出力する。
前記ニューラルネットワーク5の各素子は、前記実施
例に限定されないが、本実施例では以下に示すシグモイ
ド関数f(x)を用いている。
f(x)=1/{1+exp(−x+b)} ここで、bはバイアスであって、バイアス素子(出力
が1)のウェイトであり、学習により決定される。
ニューラルネットワーク5は予め学習しておく必要が
ある。この学習については後に詳細に説明する。ここで
は、ニューラルネットワーク5は学習が終了しているも
のとする。
第6図に示すようにMRI値はWL,WWに応じて、連続した
明るさとして表示される。ただしWL−WW/2以下のMRI値
は全て最も暗く、WL+WW/2以上のMRI値は、全て最も明
るく表示される。この明るさを便宜上第6図の縦軸に示
す数値、すなわち0.5〜16.5の連続値として表現してい
る。ここで画素値X,明るさYとし、直線部分だけに限定
すると、 Y=8.5+16×(X−WL)/WW …(1) X=WL+WW×(Y−8.5)/16 …(2) となる。
第11図は前記ニューラルネットワーク5により前記画
像度を求める手順を示すフロー図である。以下にこの手
順を説明する。まず最初に 1)手順a1では、設定されたウィンドウでのMVPの明る
さを求める。この明るさをMVPGとすると、(1)式のX
にMVPを、YにMVPGを代入して、 MVPG =8.5+16×(MVP−WL)/WW …(3) となる。(3)式を用いてMVPGを求め、それに正規化
係数1/16.5を乗算し、0〜1の値に正規化する。
2)手順b1では、手順a1と同様に設定されたウィンドウ
での“A"の明るさを求める。これをAGとすると AG=8.5+16×(A−WL)/WW …(4) となる。(4)式を用いてAGを求め、それに正規化係数
1/16.5を乗算し、0〜1の値に正規化する。
3)手順c1では、前記同様に設定されたウィンドウでの
明るさ12に相当するヒストグラムの面積を求める。
明るさ12は明るさ11.5〜12.5と考えて、(2)式に Y=11.5〜12.5を代入すると、 X12 =(3WW/16+WL)〜(WW/4+WL) …(5) となる。
したがって、明るさ12に相当するヒストグラムの面積
は(5)式の範囲のヒストグラムの面積を求めればよ
い。同様に正規化係数1/100を乗ずる。ただし、正規化
後、1を越えるものは1にする。
4)手順d1では、同様に設定されたウィンドウでの明る
さLに相当するヒストグラムの面積を求める。求められ
た明るさLは、第11図においても最も明るい明るさであ
るから、 Y=16.5に相当する。したがって、(2)式にY=1
6.5を代入して X16.5=WW/2+WL …(6) となる。したがって、X16.5以上のヒストグラムの面
積を求めれば、明るさLに相当するヒストグラムの面積
が求まる。同様に、正規化係数1/30を乗ずる。ただし、
正規化後、1を越えるものは1にする。
5)手順e1では、同様に設定されたウィンドウでの明る
さのバランスを求める。明るさのバランスは以下により
定義する。
(明るさ4〜8の面積)/(明るさ9〜13の面積) それぞれの明るさの面積は、前記手順c1と同様の方法
により求めることができる。同様に正規化の係数1/6を
乗ずる。ただし、正規化後、1を越えるものは1にす
る。
6)手順f1では、前記手順a1〜手順e1で求めたそれぞれ
の値をニューラルネットワーク5に転送する。
7)手順g1では、ニューラルネットワーク5よりニュー
ラルネットワーク5が計算した画像度を得る。
次に前述したニューラルネットワーク5の学習につい
て説明する。前述した如くニューラルネットワーク5は
予め学習しておく必要がある。学習方式は実施例に限定
されることはないが、本実施例ではバックプロパゲーシ
ョン法を用いる。ニューラルネットワーク5に一定数の
学習データを繰り返し提示し、バックプロパゲーション
法により学習を行ない、エラーが一定値以下になると、
終了する。本実施例では、以下に示す学習データを用い
て学習するが、この学習データの作成方法は実施例に限
定されることはない。
学習データは実際の複数枚の画像と、それを表示する
ために専門家が設定したウィンドウを基に作成する。
1)まず画像を1枚選択し、この画像に対して、専門家
(観察者)が設定したウィンドウをWLG,WWGとし、このW
LG,WWGに対して以下の25個の学習データを決定するため
の学習ウィンドウWLS,WWSをサンプリングする。
WLS=WLG−WWG/2 WLG−WWG/4 WLG WLG+WWG/4 WLG+WWG/2 WWS=WWG−WWG/2 WWG−WWG/4 WWG WWG+WWG/4 WWG+WWG/2 そして 2)これらのWLS,WWSに対して前記第11図に示す手順a1
〜e1によりMVP,“A"などの5つの特徴量を25組求め、学
習データの入力データを作成する。
3)またこれらのWLS,WWSに対する教師データを第12図
に示すように定め、前記2)の入力データと組み合わせ
て25組の学習データを作成する。
4)以下、適当な枚数の画像について、それぞれ25組の
学習データを作成する。
5)そして作成された全ての学習データを、エラーが一
定値以下になるまで、繰り返し学習する。
以上の学習により予め学習を行っておく。
したがって、ニューラルネットワーク5は学習が可能
なので、関数,重み係数を決定する煩わしさがなくな
る。また装置ごとに,あるいは病院ごとの特性に合った
画像を容易に観察することができる。
最後に本発明の変形例について説明する。
(1)先ず前記実施例はウィンドウレベルおよびウィン
ドウ幅の両方を変化させる方式であるが、ウィンドウレ
ベル,およびウィンドウ幅のいすれか一方を一定値に固
定し、他方を変化させるようにしても良い。
(2)さらに実施例は1枚数の画像ごとにヒストグラム
を作成し、1枚数の画像ごとに表示ウィンドウを探索し
ているが、同一グループに属する全部の画像、例えば1
回のボリュームスキャンで撮影した複数枚のスライスか
らヒストグラムを求め、それに関してウィンドウを決定
し、複数枚数の画像を同一の表示ウィンドウで観察する
ようにしても良い。この場合には1グループの画像がす
べて同一のウィンドウになる。また、ニューラルネット
ワークの教師データの作成にも、複数枚数の画像のヒス
トグラムを使用してもよい。
(3)さらに実施例では、表示ウィンドウは線形(リニ
アー)であるが、表示ウィンドウは非線形(ノンリニア
ー)であっても良い。
(4)さらに前記実施例では1枚の画像ごとにウィンド
ウを計算しているが、同一の条件で連続して撮影した画
像であれば、同一の表示ウィンドウで観察する方が良
い。このため1枚目で観察者が選択したウィンドウ(微
調整を行なった場合は微調整後のウィンドウ)とその画
像度を求めて記憶する。次に2枚目の画像の表示が指定
された場合、1枚目のウィンドウにおける2枚目の画像
の画像度を求め、それを前記記憶した1枚目の画像度と
比較し、その差が一定以下ならウィンドウを変更しない
で、1枚目のウィンドウにより2枚目の画像を表示す
る。その差が一定以上の場合には、本発明方法により求
めたウィンドウで表示するようにしても良い。
(5)さらに前記実施例では1枚の画像ごとにウィンド
ウを計算しているが、同一の条件で連続して撮影した画
像であれば、同一の表示ウィンドウで観察する方が良
い。このため同一グループの代表的な画像(例えば中心
位置の画像)でウィンドウを求め、同一グループの全部
の画像に対して同一のウィンドウで画像を表示しても良
い。
(6)さらに、前記実施例では観察目的をヒストグラム
のピークにより推定しているが、以下の情報、方法を用
いて観察部分の画素値を推定するようにしてもよい。
(1)簡単な画像認識を行ない比較的画素値の纒まって
いる部分を見つける。
(2)撮影条件、 (3)断面(アクシャル面、コロナル面、サジタル
面)、 (4)撮影目的、これは病院オーダーリングシステムか
ら得ることもできる。
また、これらの情報、ヒストグラムの特徴(全体の面
積、ピークの位置、形、高さなど)により、画像度を求
める項目、ニューラルネットワークのウエイトなどを異
なるものを使用するようにしても良い。
(7)まず前記実施例では、出力素子が1個であり、
「明るさ」,「コントラスト」の両方を統合した画像度
を表しているが、出力素子を「明るさ」,「コントラス
ト」に対して2個とし、これを「明るすぎる〜適当〜暗
すぎる」(0〜1),コントラストが「強すぎる〜適当
〜弱すぎる」(0〜1)のように表しても良い。この場
合、画像度は両者の重み付けの合計として求める。また
「明るさ」の教師データは、第12図においてWWS=WWG,
「コントラスト」の教師データはWLS=WLGとおけばよ
い。
(8)さらに実施例は、学習と実行を同一のニューラル
ネットワークで行なっているが、学習はワークステーシ
ョンなどの他のニューラルネットワークで行い、実行は
別の簡単なニューラルネットワークで行っても良い。ま
た学習したニューラルネットワークの重み係数のうち、
値の大きいものだけを取り出し、(ニューラルネットワ
ークを使用せずに)計算機のソフトウェアだけで計算し
ても良い。
(9)さらに実施例では、すべての画像について、同一
の重み係数を使用しているが、画像のタイプにより異な
る重み係数を使用しても良い。この場合は、計算機から
ニューラルネットワークに、画像タイプに応じて異なる
重み係数を転送するようにする。また学習も画像タイプ
ごとに別々の学習データを用いて行なう。
(10)さらに実施例では、学習データ用の画像の選択に
ついては説明していないが、例えば装置ごとに、10枚数
程度の画像(画像タイプごとに別々に)を学習画像とし
て記憶する。この画像によりニューラルネットワークを
学習させる。このようにすれば、装置に応じたニューラ
ルネットワーク(重み係数)が形成される。さらに新し
い学習画像を追加するか、あるいは古い画像を新しい画
像に更新して、再度全体の学習画像を用いて学習するよ
うにすれば、装置の変化にも対応可能である。
(11)さらに前記(9),(10)では、装置ごとに画像
タイプごとに別々の重み係数を設けるようにしている
が、これを装置ごと、観察者(医師)ごと,画像タイプ
ごとに別々の重み係数を設けるようにしても良い。
(12)実施例では、特徴量を線形で正規化した値を入力
としたが、正規化は例えば正弦関数や対数変換などで正
規化しても良い。
(13)実施例では、特徴の個数とニューラルネットワー
クの入力層の素子数を同一にしたが、1つの特徴に複数
の入力素子を割り当てて、特徴量を線形で正規化した値
と、正弦関数や対数変換などで正規化した値を入力して
も良い。
このほか本発明の要旨を逸脱しない範囲で種々変形実
施可能であるのは勿論である。
[発明の効果] 本発明によれば、評価値をニューラルネットワークに
より求めるので、ファジィー類似方式におけるメンバー
シップ関数等を決定することなく、装置ごとにあるいは
病院ごとの特性に合った画像を容易に観察することがで
きる。
またニューラルネットワークにより評価値を求め、最
も大きい評価値をもつウィンドウでMRI画像を観察で
き、画像を表示する際にウィンドウ設定の操作が簡単に
なり、操作者の負担を軽減できる画像表示装置を提供で
きる。
【図面の簡単な説明】
第1図は本発明に係る画像表示方法の全体の手順を示す
フロー図、第2図はウィンドウを求める方法を示すフロ
ー図、第3図は前記ウィンドウの求め方の具体的な手順
の実施例を示すフロー図、第4図は前記画像表示方法を
適用した画像表示装置を示す概略ブロック図、第5図は
典型的なMRI画像のヒストグラムを示す図、第6図はMRI
値をWL,WWに応じて連続した明るさとして表示した図、
第7図はWW,WLを変化させた場合のWW,WLに対する画像度
Qを示す概略図、第8図は前記MRI画像のピークの定義
を示す図、第9図はMVPの求め方を説明するための図、
第10図はニューラルネットワークの構成を示す図、第11
図は前記ニューラルネットワークで画像度を求める手順
を示すフロー図、第12図はWLS,WWSに対する教師データ
を示す図、第13図は従来のウィンドウの設定方法の一例
を示す図である。 1……画像メモリ、2……計算機、3……WL,WWスイッ
チ、4……表示装置、5……ニューラルネットワーク、
30……入力層、32……中間層、Q……画像度、WL,L……
ウィンドウレベル、WW,W……ウィンドウワイド、H……
背景、MVP……最も重要と推定されるMRI値。
フロントページの続き (58)調査した分野(Int.Cl.6,DB名) G06T 5/00

Claims (6)

    (57)【特許請求の範囲】
  1. 【請求項1】ディジタル画像データの各画素値を、所定
    範囲内の画素値のみが連続的に変化する明るさ値に変換
    されるように、明るさ値に変換する画像表示装置におい
    て、 ディジタル画像データから画素値のヒストグラムを求め
    る手段と、 前記ヒストグラム中の複数の所定のパラメータを入力
    し、前記ディジタル画像データを表示する際の表示画像
    の見易さを示す画像度を出力するニューラルネットワー
    クと、 前記所定範囲を変えたときの前記ニューラルネットワー
    クの出力に応じて、前記所定範囲の最適値を求める手段
    と、 を具備する画像表示装置。
  2. 【請求項2】前記所定範囲の最適値を求める手段は、前
    記所定範囲の候補を決める手段と、前記候補における画
    像度と前記候補をずらした場合における画像度とを求め
    る手段と、前記候補における画像度と候補をずらした場
    合における画像度とを比較し、候補における画像度が最
    大となるまで前記候補をずらしながら、前記候補におけ
    る画像度と前記候補をずらした場合における画像度とを
    繰返し求める手段とを具備することを特徴とする請求項
    1記載の画像表示装置。
  3. 【請求項3】前記ニューラルネットワークは学習精度が
    異なる2つの階層型のニューラルネットワークからな
    り、 前記最適値を求める手段は学習精度の低いニューラルネ
    ットワークから順に使って最適値を求め、各ニューラル
    ネットワークの出力が次のニューラルネットワークによ
    る処理における候補となることを特徴とする請求項2記
    載の画像表示装置。
  4. 【請求項4】前記ニューラルネットワークはヒストグラ
    ム中の表示画像の明るさに関する第1パラメータ、表示
    画像のコントラストに関する第2パラメータ、および明
    るさのバランスに関する第3パラメータを入力すること
    を特徴とする請求項1記載の画像表示装置。
  5. 【請求項5】前記第1パラメータは頻度が最も高い画素
    値の明るさデータとヒストグラム上で体内と体外との境
    界にある画素値の明るさデータを含み、前記第2パラメ
    ータはヒストグラム上の表示明るさ範囲の中間の明るさ
    部分の面積データと最大明るさ部分の面積データとを含
    み、前記第3パラメータはヒストグラム上の前記中間明
    るさより明るい部分の面積と暗い部分の面積の比データ
    を含むことを特徴とする請求項4記載の画像表示装置。
  6. 【請求項6】前記ニューラルネットワークは1枚の画像
    内の各部分毎の見易さを示す画像度を求め、前記頻度が
    最も高い画素値の明るさデータはヒストグラム内の各ピ
    ーク点の画素値の明るさデータを含むことを特徴とする
    請求項5記載の画像表示装置。
JP2066695A 1989-07-19 1990-03-16 画像表示方法 Expired - Lifetime JP2937389B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2066695A JP2937389B2 (ja) 1990-03-16 1990-03-16 画像表示方法
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
JP2066695A JP2937389B2 (ja) 1990-03-16 1990-03-16 画像表示方法

Publications (2)

Publication Number Publication Date
JPH03266188A JPH03266188A (ja) 1991-11-27
JP2937389B2 true JP2937389B2 (ja) 1999-08-23

Family

ID=13323333

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2066695A Expired - Lifetime JP2937389B2 (ja) 1989-07-19 1990-03-16 画像表示方法

Country Status (1)

Country Link
JP (1) JP2937389B2 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2616520B2 (ja) * 1991-08-30 1997-06-04 株式会社島津製作所 医療用画像表示装置
JPH05304622A (ja) * 1992-04-24 1993-11-16 Toshiba Corp ノイズ検出装置及びノイズリダクション装置
US6574373B1 (en) * 1999-02-12 2003-06-03 Eastman Kodak Company Method and apparatus for printing digital images
RU2431196C1 (ru) * 2010-03-31 2011-10-10 Закрытое Акционерное Общество "Импульс" Способ определения уровня яркости в зоне интереса цифрового медицинского рентгеновского изображения
CN109685866B (zh) * 2018-12-25 2023-07-18 上海联影医疗科技股份有限公司 晶体位置查找表构建方法、装置、计算机设备和存储介质

Also Published As

Publication number Publication date
JPH03266188A (ja) 1991-11-27

Similar Documents

Publication Publication Date Title
EP0409206B1 (en) Digital image display apparatus
CN100470560C (zh) 图像像素编码方法和图像处理方法
CN110399929B (zh) 眼底图像分类方法、装置以及计算机可读存储介质
CN1475969B (zh) 用于增强人像图像的方法和系统
Perner An architecture for a CBR image segmentation system
DE69629445T2 (de) Automatische Tonskalenabstimmung mittels Bildaktivitätsmessungen
US7869642B2 (en) Window leveling system and method
US7315638B2 (en) Image display system
KR100877457B1 (ko) 신호 처리 장치 및 방법과 신호 처리 프로그램을 기록한 기록매체
CN101959018A (zh) 图像处理设备、图像处理方法和程序
CN108334733B (zh) 医学图像显示方法、显示系统及计算机可读存储介质
EP1477929B1 (en) Signal processing
Moller-Jensen Classification of urban land cover based on expert systems, object models and texture
US20070127783A1 (en) Image processing apparatus, method and program for controlling flesh color of image
JP4303092B2 (ja) 年齢推定装置、年齢推定方法および年齢推定プログラム
Mouzon et al. Joint CNN and variational model for fully-automatic image colorization
JP7437738B2 (ja) 機械学習システム、学習用データセット生成システム、及び機械学習プログラム
JP2937389B2 (ja) 画像表示方法
JP2937390B2 (ja) 画像変換方法
Davies et al. Visual saliency as an aid to updating digital maps
JPH0464180A (ja) ディジタル画像表示装置
Yeganeh Cross dynamic range and cross resolution objective image quality assessment with applications
JP3198099B2 (ja) 画像処理方法及びその装置並びに記憶媒体
US20060098887A1 (en) Mehthod for image conversion
JPH0833722B2 (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: R313114

Free format text: JAPANESE INTERMEDIATE CODE: R313111

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: R313114

Free format text: JAPANESE INTERMEDIATE CODE: R313111

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: 20090611

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100611

Year of fee payment: 11

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100611

Year of fee payment: 11