JP3748457B2 - Window operation device - Google Patents

Window operation device Download PDF

Info

Publication number
JP3748457B2
JP3748457B2 JP32040394A JP32040394A JP3748457B2 JP 3748457 B2 JP3748457 B2 JP 3748457B2 JP 32040394 A JP32040394 A JP 32040394A JP 32040394 A JP32040394 A JP 32040394A JP 3748457 B2 JP3748457 B2 JP 3748457B2
Authority
JP
Japan
Prior art keywords
window
amount
value
δwl
histogram
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 - Fee Related
Application number
JP32040394A
Other languages
Japanese (ja)
Other versions
JPH08180181A (en
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 JP32040394A priority Critical patent/JP3748457B2/en
Publication of JPH08180181A publication Critical patent/JPH08180181A/en
Application granted granted Critical
Publication of JP3748457B2 publication Critical patent/JP3748457B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

【0001】
【産業上の利用分野】
本発明は、例えば医用画像撮影装置にて撮影された画像データのウインドウ幅、ウインドウレベルが好適となるように制御するウインドウ操作装置に関する。
【0002】
【従来の技術】
医用画像診断装置にて撮影される例えばX線透視画像等の医用画像は、診察する際に濃淡値諧調変換(以下、ウインドウ変換という)を加えて、診察対象となる部位の諧調を際立たせる。ウインドウ変換には、ウインドウレベル(WL)の調整と、このウインドウレベルを中心としてどの画素値まで諧調をつけるかを決めるウインドウ幅(WW)の調整とがあり、通常はこの2つを調整する。つまり、撮影されたX線画像内には、骨部、臓器部等それぞれ画素値の異なる部位が存在するので、当該画像上で所望の臓器を観察したい場合にはこの画素値近傍にウインドウレベル、ウインドウ幅を設定してこの部位の濃淡を際立たせれば良い。
【0003】
このようなウインドウ変換の従来方法としては以下に示す4通りの方法があり、以下に説明する。
【0004】
(a) 回転ダイヤルを使用する方式
(b) スナップスイッチを使用する方式
(c) 置数キーを設ける方式
(d) ソフトウエアによる制御方式
<(a)回転ダイヤルを使用する方式>
この方式は、回転ダイヤルの回転量に比例した分だけウインドウ幅、ウインドウレベルを変化させるものであり、ダイヤルの回転をロータリーエンコーダにて検出し、この検出結果に定数を乗じてウインドウ幅、ウインドウレベルの調整量を決定する。また、回転ダイヤルの回転速度を検出し、この検出結果に基づき回転速度が速いときにはウインドウ幅、ウインドウレベルの変化量が大きくなるように制御することもできる。
【0005】
<(b)スナップスイッチを使用する方法>
この方式は、スナップスイッチの切換えによりウインドウ幅、ウインドウレベルの変化量を変換するものである。即ち、ウインドウ幅、ウインドウレベルの切換えスイッチと、変化速度の速い・遅いを切換えるスイッチと、増加・減少を切換えるスイッチを配設し、各スイッチ操作により8通りの操作モードを設定することたできる。例えば、ウインドウ幅を速い速度で増加させる方向に変化させる等の設定が可能である。
【0006】
<(c)置数キーを設ける方式>
この方式は、多数の置数キーを設定し、各置数キーに従ってウインドウ幅、ウインドウレベルを決定するものである。
【0007】
<(d) ソフトウエアによる制御方式>
画面上にウインドウ幅、ウインドウレベルを文字やグラフィックで表示し、マウスやキーボードで操作する方式であり、データ処理的には前記した(b)、(c)の方式と同一であるが、ウインドウ変換をハードウエアを使用せずにコンピュータ画面を見ながらマウスやキーボードを用いて操作する方法である。
【0008】
【発明が解決しようとする課題】
上記した各従来方法によれば以下に示す如くの欠点がある。
【0009】
(a)回転ダイヤルを使用する方式においては、X線CT画像等の医用画像データは一般に、−2048〜2047迄の広範囲な帯域のデータを有しており、しかし、濃淡値として表示するのは、ほんの一部である。例えば、CT画像では空気は−1000、水は0、内臓組織は一般に+40程度、骨組織は+500程度から+1000程度の範囲であり、通常は主に内臓組織に濃淡値表示帯域を設定する。
【0010】
通常は画像再構成時に関心濃淡値を決定し画像情報に設定するが、骨組織/脳組織/内臓のいずれかを診断するかは画像再構成時には厳密には決定しておらず、内臓組織と骨組織、脳組織と骨組組織等へのウインドウ設定を繰り返す事がある。
【0011】
また、回転ダイヤル方式のウインドウ操作(操作速度によってウインドウ変化倍率を変化させる方式)では、回転ダイヤルをいきおいよく回す事によって所望のウインドウ値に速やかに近付ける事が出来る。しかし、表示帯域は一般に狭いため簡単に行き過ぎてしまい、回転ダイヤルを逆回転させたり、ゆっくり操作するなどを行う事により所望のウインドウ値に設定しなければならない。この様に、大きく離れた濃淡値帯域間の設定を繰り返す際の操作の制度が悪い。
【0012】
(b)スナップ・スイッチを使用する方式は、操作速度のバリエーションが予め用意した種類のみであり操作に制限がある。そこで、スイッチを増やすと操作速度のバリエーションを増やす事が出来るが、操作が煩雑になる。
【0013】
この様にスイッチ数に関して、操作の煩雑/操作性から増やしても、減らしても操作性が低下する。実際の装置では、本方式をウインドウ値の微調整のみに使用し大きく離れたウインドウ値の設定には(c)の置数キーを設ける方式を併用している。
【0014】
(c)多数の置数キーを設ける方式では、WL/WWの必要な置数キーを揃えるだけで操作パネル上が置数キーで一杯になってしまうという問題があり、実際の装置では(b)方式と併用しそれぞれバランスをとり実用的な装置を開発する事が出来たが、複数種のスイッチを多数操作するため煩雑で操作性が悪い。
【0015】
(d)ソフトウエア制御方式では、キーボードからのウインドウ値入力、変化量入力、或いは、マウスによる操作等操作が煩雑であるという問題がある。
【0016】
この発明はこのような従来の課題を解決するためになされたもので、その目的とするところは、容易に所望のウインドウ幅、ウインドウレベルに設定することのできるウインドウ操作装置を提供することにある。
【0017】
【課題を解決するための手段】
上記目的を達成するために、請求項1記載の発明は、入力される画像データのウインドウ幅及びウインドウレベルを操作手段を調整して所望の値に設定するウインドウ操作装置において、
表示されている画像のウインドウ幅の大きさに応じて前記操作手段における操作量に対する前記ウインドウ幅の変化量を変化させる制御手段を有し、
前記制御手段は、前記画像データ全体の濃淡発生頻度分布又は予め設定された濃淡値発生頻度分布に基づき、データが多く発生している濃淡値近傍にて前記操作手段における操作量を緩やかとすべく制御することを特徴とするものである。
【0023】
上述の如く構成された請求項1の発明によれば、画像データのヒストグラムに応じてウインドウ幅、ウインドウレベルが設定される。例えば、ヒストグラムの値の小さい画素値においては、この画素値にウインドウレベルを設定する可能性は少ないので操作量に対する変化量を大きくする。また、反対に、ヒストグラムの値が大きい画素値においては、この画素値にウインドウレベルを設定する可能性は大きく、微調整が容易なように、操作量に対する変化量を小さくする。これによれば、微調整が必要なウインドウ幅が小さい領域にて変化量が小さくなるので微調整が容易となり、反対に、ウインドウ幅が大きい領域では変化量がおおきくなるのでウインドウ幅の変更が容易となる。
【0024】
請求項2の発明によれば、ヒストグラムに対応した操作感度代表値が求められる。ここで、操作感度代表値とは、ヒストグラムの大きさに応じて決められるウインドウレベルである。そして、操作スイッチの操作速度に基づき一次変換量が求められ、更に、この一次変換量と操作感度代表値から二次変換量が求められる。従って、操作者の操作速度に応じてウインドウレベルの変化量が切替わり、目標のウインドウレベルに容易に到達させることができる。
【0029】
【実施例】
以下、本発明の実施例を図面に基づいて説明する。図1は本発明の第1実施例に係るウインドウ操作装置の構成を示すブロック図である。同図に示すように、このウインドウ操作装置は、X線撮影装置等の医用画像診断装置にて撮影された医用画像を取り込むインターフェース1と、取り込まれた医用画像データを格納する画像メモリ2と、該画像メモリ2に格納された画像データのヒストグラムを計算するヒストグラム計算部3と、計算されたヒストグラムデータを格納するヒストグラムメモリ7と、予め設定された複数のヒストグラムデータが格納される既定義ヒストグラムメモリ8と、前記ヒストグラムメモリ7及び既定義ヒストグラムメモリ8のいづれか一方のデータを選択するヒストグラム選択部9とを有している。また、画像メモリ2の出力は図16に示す表示ユニット側に画像データを出力する。
【0030】
ヒストグラム計算部3は、例えば、図2(a)に示すように、画素値毎のデータ発生頻度をヒストグラムとして求める処理を行う。
【0031】
既定義ヒストグラムメモリ8は、予め定義されたヒストグラムデータを記憶させるものであり、例えば、被検体頭部や胸部の画像のヒストグラムデータが格納されている。
【0032】
また、このウインドウ操作装置は、前記したヒストグラム選択手段にて選択されたヒストグラムデータに基づきウインドウ幅及びウインドウレベルの操作感度(操作量に対する変化量の割合)の代表値を決定する操作感度代表値計算部10と、この代表値を格納する操作感度代表値メモリ11と、ウインドウレベルの操作を行うWL回転スイッチ4と、該スイッチ4の操作量ΔWLから速度一次変換値ΔWL′(後述)を求める速度一次変換部5と、更に速度二次変換値ΔWL″を求める速度二次変換部6と、加算部12と、ウインドウレベルを決定して出力するウインドウレベル決定部13とを有している。
【0033】
操作感度代表値計算部10は、ヒストグラムデータhist(val) を基に、以下に示す(1)式により操作感度代表値oratio(val) を求める。
【0034】
【数1】
oratio(val) =max (hist(data i ) )*k/hist(val) …(1)
iは画像データ数、kは定数、val は画素値であり例えば−211〜211−1までの諧調を有する。
【0035】
つまり、この(1)式では、画像データの大きさに対する操作感度を示しており、ヒストグラムデータhist(val) が小さいところでは操作感度代表値oratio(val) が大きく、少ない操作量で大きく変化する。また、反対にhist(val) が大きいときにはoratio(val) は小さくなるので、微調整が容易となる。(1)式を図に示すと図3の如くとなり、データ頻度の小さい画素値近傍では操作感度が極端に大きくなる。
【0036】
そこで、これを防止する為に操作感度代表値oratio(val) にしきい値を設ける。その結果、例えば、図4に示す如くの特性曲線が得られる。このようにすれば、操作感度の極端な変化を防止することができるようになる。また、これを図2(a)に示したヒストグラムに当てはめると、同図(b)に示す如くとなり、画素値の頻度の高いところで操作感度代表値F(x) が小さくなっていることが理解できる。
【0037】
また、ヒストグラムメモリのデータサイズを減らし、またウインドウ音化をスムースにするために、ヒストグラムデータを薮点で平均化し、使用時に補間あるいは、そのまま使用する。
【0038】
速度一次変換部5は、回転スイッチの回転速度を強調したウインドウレベルの変化量ΔWL′を求めて出力するものであり、例えば、単位時間当たりの変化量ΔWLの二乗の値を出力する。即ち、以下に示す(2)式である。
【0039】
ΔWL′=ΔWL*|ΔWL| …(2)
速度二次変換部6は、速度一次変換部5にて求められた変化量ΔWL′と、操作感度代表値関数oratioから、次の(3)式にて二次変化量ΔWL″を求める。
ΔWL″=ΔWL′*oratio(WL-1)…(3)
ここで、oratio(WL-1)は、時系列的に変化する操作感度代表値の前回の数値である。
【0040】
加算部12は、前回のウインドウレベルWL-1にΔWL″を加算して今回のウインドウレベルWLを求める。
【0041】
ウインドウレベルメモリ13は、加算部12にて求められたウインドウレベルWLを記憶し出力するものである。
【0042】
また、図16に示す表示ユニットの画像読み出し部23は、画像メモリ2に画像データのアドレス信号を供給して、このアドレスに対応した画像データを読み出す。ウインドウ関数発生部24は、本願発明のウインドウ操作装置にて決定されたウインドウ幅、ウインドウレベルに基づいて、ウインドウ関数をルックアップテーブルデータとして発生する。濃淡値変換部25は入力の12bit 画像をルックアップテーブルに基づき濃淡値への変換処理を行う。画像処理部26は、濃淡値が変換された後の画像信号をビデオ信号に変換する等の処理を加えてモニタ27に出力するものである。
【0043】
このように構成された本発明の第1実施例の動作を以下に説明する。まず、インターフェース1から医用画像が入力されると、この画像データは一旦画像メモリ2に格納され、表示ユニットに送られてモニタ27にて表示される。また、この画像メモリ2内に格納された画像データはヒストグラム計算部3にてヒストグラムが計算され、ヒストグラムメモリ7内に格納される。
【0044】
そして、ヒストグラム選択部9にて、計算されたヒストグラムまたは既定義のヒストグラムのいずれかが選択され、操作感度代表値計算部10に送られる。操作感度代表値計算部10では、前記した(1)式に従って操作感度代表値を演算し、この演算結果を操作感度代表値メモリ11内に格納する。
【0045】
いま、WL回転スイッチ4が操作されると速度一次変換部5にてスイッチ4の変化量ΔWLに基づき(2)式に従って一次変化量ΔWL′が求められ、更に、速度二次変換部6にて、二次変換量ΔWL″が(3)式よって求められる。このときの(3)式における「oratio(WL-1)」の値は、操作感度代表値計算部10に記憶された値である。
【0046】
その後、加算部12において次の(3)式に示す加算が行われ今回のウインドウレベルWLが求められる。
WL=WL-1+ΔWL″ …(4)
そして、このウインドウレベル値WLはウインドウレベル決定部13より出力される。その結果、図5に示す如くの操作速度、操作感度、画素値との関係が得られる。
【0047】
このようにして、本実施例によれば、画像のヒストグラムデータ(実際の画像データから求めたもの又は既定義のもの)、及びWL回転スイッチの操作量に基づいてウインドウレベルの変化量を可変としている。即ち、ヒストグラムデータから、画素値の多い部位については変化量を小さくし、画素値の少ない部位については変化量が大きくなるように設定される。また、回転スイッチの回転操作速度が速いときにはウインドウレベルの変化量が大きくなるように設定される。従って、操作者の操作に対して感応的にウインドウレベルが変化するようになり、操作性が著しく向上する。
【0048】
次に本発明の第2実施例について説明する。図6は本発明の第2実施例に係るウインドウ操作装置の構成を示すブロック図である。図示のように、この実施例では、指定組織帯域中心算出部21、指定組織帯域中心地メモリ22が、それぞれ操作感度代表値計算部10、操作感度代表値メモリ11の代わりに配置されている点で前記した第1実施例と異なる。
【0049】
指定組織帯域中心算出部21は、ヒストグラムデータhist(val) 、及び操作者による設定帯域から次の(5)式に従って指定組織帯中心値i-centerを計算し、
i-center=max(hist(val)) …(5)
val は帯域指定範囲内(内臓ならば0〜50、量ならば500〜1500等)そして、速度二次変換部からのWL-1を元に感度を求め、速度二次変換部へ感度を設定する。
感度=|i-center−WL-1
【0050】
つまり、この第2実施例では、関心のある所望部位を指定するとこの指定された部位についてのヒストグラムに対して操作感度が変化するよう設定される、例えば図7(a)に示す如くのヒストグラムが与えられた際に、操作者が特に内臓組織について観察したい場合には、この内臓組織について帯域を設定すると、この部位についてのヒストグラムに対して操作感度の変化量の設定が行われ、同図(b)に示す如くの特性曲線が得られる。また、速度二次変換部の処理によって図8に示す如くの操作速度、操作感度、画素値の関係を得ることができる。
【0051】
このようにして、第2実施例によれば、所望の帯域を設定するとこの部位についてのヒストグラムに対して操作感度の設定が行われるので、観察したい組織が決まっている際には特に有用である。
【0052】
図9は本発明の第3実施例に係るウインドウ操作装置の構成を示すブロック図である。この実施例では、現在のウインドウ幅WWと、ウインドウ幅の操作量に基づいて操作感度を設定するものである。同図におけるWW回転スイッチ31を回転させると、この変化量WWからこの回転速度ΔWWが求められ、速度一次変換部32にて次の(6)式によりこの回転速度ΔWWを強調したWW値の変化速度ΔWW′が求められる。
ΔWW′=ΔWW*|ΔWW| …(6)
【0053】
その結果、図10に示す如くの一次関数的な変化を示す感度特性が、図12に示すように二次関数的となり、ウインドウ幅WWが大きくなるにつれて操作感度の増加率が上昇する。
【0054】
その後、速度二次変換部33ではWWメモリ35に記憶されている前回のウインドウ幅WW-1に基づき、二次変換量ΔWW″を次の(7)式にて求める。
ΔWW″=ΔWW′*(1+WW-1)/B …(7)
Bは定数
【0055】
その後、加算部34にて二次変換量ΔWW″を前回のウインドウ幅WW-1に加算して今回のウインドウ幅を求め出力する。その結果、WW回転スイッチ31の回転速度に応じた操作感度を得ることができるようになる。この際の特性曲線は、図10に示した操作感度に対しては図11の如くとなり、図12に示した操作感度に対しては図13のようになる。同図から明らかなように、回転スイッチの回転速度が速いほど操作感度が大きくなっている。
【0056】
次に本発明の第4実施例について説明する。図14は、第4実施例に係る構成図であり、前記した図9の実施例と比べて感度補正部36を具備する点で異なる。
【0057】
感度補正部36は、前回のウインドウ幅WW-1からΔWWだけ増加して今回のウインドウ幅となる際に例えば1000、100等のきりの良い数値を通過するときには、このきりの良い数値付近にて操作量に対する変化量が低下するように、以下に示す如くの処理を行う。
【0058】
今、変化量を低下させたい帯域の中心値をkvaln (nは1〜M、Mはきりのよい値として設定する総点数)とし、Cn をkvaln 前後の低下させたい帯域の半値とし、ウインドウ幅の変化量をΔWW''' 、処理後のウインドウ幅をWW''' とすると、以下の如くとなる。
【数2】
(kvaln −Cn )<WW-1<(kvaln +Cn )の範囲において、
(kvaln −Cn )<WW-1 ならば、
ΔWW''' =ΔWW″*(WW-1−kvaln )/Cn
(kvaln −Cn )>WW-1 ならば、
ΔWW''' =ΔWW″*(kvaln − WW-1)/Cn
(WW''' /WW″)<0.5であれば、
ΔWW''' =ΔWW″*0.5
上記以外のとき;ΔWW''' =ΔWW″ …(8)
【0059】
こうして、上記の方法を用いてウインドウ幅WW''' が求められ、kvaln 近傍にてウインドウ幅の操作感を穏やかにすることができるようになる。例えば、kvaln 、Cn を、
n =int(kvaln 0.5 )
とすると、図15(a)に示す如くのデータが得られ、操作感度特性は同図(b)のように、各きりの良いウインドウ幅の地点で操作感度が小さくなるように設定され、この付近における操作が容易になる。
【0060】
このように、ウインドウ幅の設定は通常、100や1000等きりの良い値に設定することが多く、第4実施例では、ウインドウ幅がきりの良い値の近傍にて操作感度が低下するので、操作性が向上する。
【0061】
図17は、本発明の第5実施例を示すブロック図である。この実施例の構成は、前記した第1実施例とほぼ同一であり、感度補正部41が新たに配設されている。
【0062】
感度補正部41は、前回のウインドウレベルWL-1にΔWLが加えられ今回のウインレベルに変化する際に、例えば、1000、100等のきりの良い数値を通過する際には、このきりの良い数値付近にて操作量に対する変化量が低下するように、以下に示す如くの処理を行う。
【数3】
(kvaln −Cn )<WL-1<(kvaln +Cn )の範囲において、
(kvaln −Cn )<WL-1 ならば、
ΔWL''' =ΔWL″*(WL-1−kvaln )/Cn
(kvaln −Cn )>WL-1 ならば、
ΔWL''' =ΔWL″*(kvaln − WL-1)/Cn
(WL''' /WL″)<0.5であれば、
ΔWL''' =ΔWL″*0.5
上記以外のとき;ΔWL''' =ΔWL″ …(9)
【0063】
こうして、上記の方法を用いてウインドウレベルWL''' が求められ、kvaln 近傍にてウインドウ幅の操作感を穏やかにすることができるようになる。
【0064】
このようにして、第5実施例では100や1000のようにきりの良いレベルにて操作感度が低下するように設定されるので、きりの良い所望の値にウインドウレベルを容易に設定することができるようになる。
【0065】
また、図2に示した第2実施例について感度補正部を設け、前記第5実施例と同様の方法を用いたウインドウレベルの設定をしても良い。この場合の構成は、図18のようになり、動作は前述した通りである。
【0066】
また、通常、肺、内臓、骨等の組織部分にウインドウレベルを設定することが多いので、ウインドウレベルの変化量を低下させたい値kvaln を予め−1000(肺)、40(内臓)、1000(骨)等の値に設定しても良い。このようにすれば、図22に示すような感度特性を得ることができる。
【0067】
次に本発明の第6実施例について説明する。この実施例では、ウインドウレベルを操作する回転スイッチを連続して操作すると徐々に変化速度が増加するように制御するものである。そして、この実施例では、図1、図6に示した速度一次変換部5における処理が異なる。以下、図1を用いて具体的に説明する。
【0068】
まず、WL回転スイッチ4の入力データΔWLの今回値ΔWLn と前回値ΔWLn-1 との差Δn ΔWLを次の(10)式を用いて求める。
Δn ΔWL=ΔWLn −ΔWLn-1 …(10)
【0069】
そして、このΔn ΔWLの値が一定幅で連続している際には(例えば、m回)速度一次変換部5の出力を次の(11)式に示すように増加させる。
【数4】
ΔWLn ′=ΔWLn *|ΔWLn |*f(m)…(11)
例えば、f(m)=mk
kは定数である
【0070】
このような設定によれば、操作者が連続にWL回転スイッチ4を操作すると、この連続している時間に対応して操作量に対する変化量が大きくなる。従って、例えばウインドウレベルを1000変更する際にもスイッチ4の操作量を大きくする必要はなく、緩やかな回転を連続しているだけで容易に所定のウインドウレベルに到達させることができるようになる。
【0071】
次に、本発明の第7実施例について説明する。図19は、第7実施例に係る構成図であり、図1に示した第1実施例の速度一次変換部5の構成が帯域判定部41、連続性計数部42、演算部43から構成されている点で第1実施例とは異なっている。
【0072】
帯域判定部41は、回転ダイヤル4の操作量ΔWLと所定の下限値TL との比較を行い、現在どの帯域に属するかを判定する。そして、|ΔWL|が連続して下限値TL を越えた際にはこの回数mが連続性計数部42にてカウントされる。また、下限値TL 以下となった際にはmを初期化(即ち、「1」)する。
【0073】
そして、この回数mを用いて、演算部43では次の(12)式にて一次変換量ΔWL′を演算する。
【数5】
ΔWLn ′=ΔWLn *|ΔWLn |*f1 (m) …(12)
例えば、f1 (m)=mk
kは定数である
ここで、(12)式による処理を加えないときのΔWLの特性図、ΔWL′の特性図を図23に示し、(12)式による処理を加えたときのΔWLの特性図、ΔWL′の特性図を図24に示す。
【0074】
そして、このように構成された第7実施例では、前記第6実施例と同様の効果を得ることができると共に、操作量ΔWLが下限値TL を連続して越えているときに一次変換量ΔWL′が大きくなるように制御され、反対に操作量ΔWLが下限値TL よりも小さくなったときには初期値に戻るので、操作者の操作意志に適応したウインドウレベルの調整を行うことができるようになる。
【0075】
次に、第7実施例の変形例について説明する。図20は該変形例の構成を示すブロック図であり、図19に示した帯域判定部41が帯域ゼロ判定部44となっている点で異なっている。
【0076】
帯域ゼロ判定部44は、回転スイッチ4の操作量がゼロであるかどうかを判定し、ゼロでない場合が複数回連続したときは連続性計数部42にてこの回数mを計数する。その他については前記した第7実施例と同様に、演算部43では以下の(13)式にて一次変換量ΔWL′を演算する。
【数6】
ΔWLn ′=ΔWLn *|ΔWLn |*f2 (m) …(13)
2 (m)=mk
kは定数である
ここで、(13)式による処理を加えたときのΔWLの特性図、ΔWL′の特性図を図25に示す。
【0077】
そして、この変形例によれば、前記第7実施例が所定の下限値TL を用いて、ΔWLがこの下限値TL よりも大きいかどうかを比較したのに対して、ΔWLがゼロであるかどうかを見て連続に回転スイッチ4を操作しているかどうかを判定している。従って、この変形例は、第7実施例の下限値TL をゼロとした場合である。
【0078】
また、第7実施例の方法と変形例の方法とを併用して一次変換量ΔWLを求めることも可能である。この場合は図21に示す如くの構成となり、演算部43での演算式は次の(14)式となる。
【0079】
【数7】
ΔWLn ′=ΔWLn *|ΔWLn |*f1 (m)*f2 (n) …(14)
このような構成によれば、より一層感応的な操作が可能となる。
【0080】
【発明の効果】
以上説明したように、本発明によれば、微調整が必要なウインドウ幅が小さい領域にて変化量が小さくなるので微調整が容易となり、反対に、ウインドウ幅が大きい領域では変化量が大きくなるのでウインドウ幅の変更が容易となる。
【図面の簡単な説明】
【図1】本発明の第1実施例の構成を示すブロック図。
【図2】第1実施例に係り、腹部画像を例とした濃淡帯域と操作感度の関係を示す特性図。
【図3】画素値と操作感度の関係を示す特性図。
【図4】しきい値を設定したときの画素値と操作感度の関係を示す特性図。
【図5】第1実施例に係る操作速度、操作感度、画素値の関係を示す特性図。
【図6】本発明の第2実施例の構成を示すブロック図。
【図7】第2実施例に係り、腹部画像を例とした濃淡帯域と操作感度の関係を示す特性図。
【図8】第2実施例に係る操作速度、操作感度、画素値の関係を示す特性図。
【図9】本発明の第3実施例の構成を示すブロック図。
【図10】第3実施例に係り、ウインドウ幅と操作感度の関係を示す第1の特性図である。
【図11】第3実施例に係り、ウインドウ幅、操作速度、操作感度の関係を示す第1の特性図。
【図12】第3実施例に係り、ウインドウ幅と操作感度の関係を示す第2の特性図である。
【図13】第3実施例に係り、ウインドウ幅、操作速度、操作感度の関係を示す第2の特性図。
【図14】本発明の第4実施例の構成を示すブロック図。
【図15】第4実施例に係り、ウインドウ幅と感度の関係を示す特性図である。
【図16】表示ユニットの構成を示すブロック図である。
【図17】本発明の第5実施例の構成を示すブロック図。
【図18】第5実施例の変形例の構成を示すブロック図。
【図19】本発明の第7実施例の構成を示すブロック図。
【図20】第7実施例の変形例の構成を示すブロック図。
【図21】第7実施例の他の変形例の構成を示すブロック図。
【図22】第5実施例に係る処理を加えたときのウインドウレベルと操作感度の関係を示す特性図。
【図23】第7実施例に係る処理を加えないときのΔWLおよびΔWL′の変化を示す特性図。
【図24】第7実施例に係る処理を加えたときのΔWLおよびΔWL′の変化を示す特性図。
【図25】第7実施例の変形例に係る処理を加えたときのΔWLおよびΔWL′の変化を示す特性図。
【符号の説明】
1 インターフェース 2 画像メモリ 3 ヒストグラム計算部
4 ウインドウレベル回転スイッチ 5,32 速度一次変換部
6,33 速度二次変換部 7 ヒストグラムメモリ
8 既定義ヒストグラムメモリ 9 ヒストグラム選択部
10 操作感度代表値計算部 11 操作感度代表値メモリ
12,34 加算部 13 ウインドウレベルメモリ
21 指定組織帯域中心算出部 22 指定組織帯域中心値メモリ
31 ウインドウ幅回転スイッチ 35 ウインドウ幅メモリ
36,41 感度補正部
[0001]
[Industrial application fields]
The present invention relates to a window operation device that controls a window width and a window level of image data captured by, for example, a medical image capturing device.
[0002]
[Prior art]
A medical image such as an X-ray fluoroscopic image taken by a medical image diagnostic apparatus, for example, is subjected to grayscale gradation conversion (hereinafter referred to as window conversion) at the time of medical examination to highlight the gradation of a part to be examined. The window conversion includes an adjustment of the window level (WL) and an adjustment of the window width (WW) that determines to which pixel value the gradation is to be centered around the window level. Usually, these two are adjusted. That is, since there are parts with different pixel values, such as bone parts and organ parts, in the captured X-ray image, if you want to observe a desired organ on the image, the window level near the pixel value, What is necessary is just to set the window width and make the shade of this part stand out.
[0003]
As a conventional method of such window conversion, there are the following four methods, which will be described below.
[0004]
(A) Method using a rotary dial
(B) Method using a snap switch
(C) Method of providing numeric keys
(D) Control method by software
<(A) Method using a rotary dial>
This method changes the window width and window level by an amount proportional to the amount of rotation of the rotary dial. The rotation of the dial is detected by a rotary encoder, and the detection result is multiplied by a constant to obtain the window width and window level. Determine the adjustment amount. It is also possible to detect the rotational speed of the rotary dial, and to control the amount of change in the window width and window level to be large when the rotational speed is high based on the detection result.
[0005]
<(B) Method of using snap switch>
In this method, the change amount of the window width and the window level is converted by switching the snap switch. That is, a switch for switching the window width and window level, a switch for switching fast / slow change speed, and a switch for switching increase / decrease can be provided, and eight operation modes can be set by operating each switch. For example, setting such as changing the window width in a direction of increasing at a high speed is possible.
[0006]
<(C) Method of providing a numeric key>
In this method, a large number of numeric keys are set, and the window width and window level are determined according to each numeric key.
[0007]
<(D) Control method by software>
This is a method of displaying the window width and window level on the screen as characters and graphics, and operating with a mouse or keyboard. The data processing is the same as the methods (b) and (c) described above. Is operated using a mouse or keyboard while looking at the computer screen without using hardware.
[0008]
[Problems to be solved by the invention]
Each of the above conventional methods has the following drawbacks.
[0009]
(A) In a method using a rotary dial, medical image data such as an X-ray CT image generally has data in a wide band from -2048 to 2047, but it is displayed as a gray value. , Just a part. For example, in the CT image, air is -1000, water is 0, visceral tissue is generally in the range of about +40, and bone tissue is in the range of about +500 to about +1000. Usually, a gray value display band is mainly set in the visceral tissue.
[0010]
Normally, the gray value of interest is determined and set in the image information at the time of image reconstruction. However, whether to diagnose bone tissue / brain tissue / internal organ is not strictly determined at the time of image reconstruction. The window setting for bone tissue, brain tissue and frame tissue may be repeated.
[0011]
Further, in the case of a rotary dial type window operation (a method in which the window change magnification is changed according to the operation speed), a desired window value can be quickly approached by rotating the rotary dial. However, since the display band is generally narrow, it easily goes too far, and the desired window value must be set by rotating the rotary dial in reverse or by slowly operating it. Thus, the system of operation when repeating the setting between the gray value bands that are far apart is poor.
[0012]
(B) The method using the snap switch is limited to the operation speed variations, and the operation is limited. Therefore, increasing the number of switches can increase the variation in operation speed, but the operation becomes complicated.
[0013]
As described above, regarding the number of switches, the operability deteriorates even if the number of switches is increased or decreased due to the complicated operation / operability. In an actual apparatus, this method is used only for fine adjustment of a window value, and a method of providing a numeric key of (c) is used in combination for setting a window value far apart.
[0014]
(C) In the system in which a large number of numeric keys are provided, there is a problem that the numeric keys on the operation panel are filled with the numeric keys required for WL / WW. ) And a practical device that can be balanced with each other and developed. However, since a large number of types of switches are operated, it is cumbersome and poor in operability.
[0015]
(D) In the software control method, there is a problem that operations such as window value input from the keyboard, change amount input, or mouse operation are complicated.
[0016]
The present invention has been made to solve such a conventional problem, and an object of the present invention is to provide a window operating device that can easily set a desired window width and window level. .
[0017]
[Means for Solving the Problems]
In order to achieve the above object, the invention according to claim 1 is a window operating device for adjusting a window width and a window level of input image data to desired values by adjusting an operation means.
Control means for changing the amount of change in the window width with respect to the operation amount in the operation means according to the size of the window width of the displayed image;
The control means is configured to moderate the operation amount in the operation means in the vicinity of the gray value where a lot of data is generated based on the gray frequency occurrence frequency distribution of the entire image data or a preset gray value occurrence frequency distribution. It is characterized by controlling.
[0023]
According to the invention of claim 1 configured as described above, the window width and the window level are set according to the histogram of the image data. For example, in a pixel value having a small histogram value, since there is little possibility of setting a window level for this pixel value, the change amount with respect to the operation amount is increased. On the other hand, for a pixel value having a large histogram value, the window level is likely to be set for this pixel value, and the change amount with respect to the operation amount is reduced so that fine adjustment is easy. According to this, since the amount of change is small in the region where the window width that needs fine adjustment is small, the fine adjustment is easy. On the other hand, the amount of change is large in the region where the window width is large, so the window width can be easily changed. It becomes.
[0024]
According to the invention of claim 2, the operation sensitivity representative value corresponding to the histogram is obtained. Here, the operation sensitivity representative value is a window level determined according to the size of the histogram. Then, the primary conversion amount is obtained based on the operation speed of the operation switch, and the secondary conversion amount is obtained from the primary conversion amount and the operation sensitivity representative value. Therefore, the amount of change in the window level is switched according to the operation speed of the operator, and the target window level can be easily reached.
[0029]
【Example】
Embodiments of the present invention will be described below with reference to the drawings. FIG. 1 is a block diagram showing the configuration of the window operating device according to the first embodiment of the present invention. As shown in the figure, the window operating device includes an interface 1 for capturing a medical image captured by a medical image diagnostic apparatus such as an X-ray imaging apparatus, an image memory 2 for storing the captured medical image data, A histogram calculation unit 3 for calculating a histogram of image data stored in the image memory 2, a histogram memory 7 for storing calculated histogram data, and a predefined histogram memory for storing a plurality of preset histogram data 8 and a histogram selection unit 9 for selecting one of the data of the histogram memory 7 and the predefined histogram memory 8. Further, the image memory 2 outputs image data to the display unit shown in FIG.
[0030]
For example, as shown in FIG. 2A, the histogram calculation unit 3 performs a process of obtaining the data occurrence frequency for each pixel value as a histogram.
[0031]
The predefined histogram memory 8 stores predefined histogram data, and stores, for example, histogram data of images of the subject's head and chest.
[0032]
Further, the window operating device calculates an operation sensitivity representative value for determining a representative value of the window width and window level operation sensitivity (the ratio of the change amount to the operation amount) based on the histogram data selected by the histogram selecting means. Unit 10, an operation sensitivity representative value memory 11 for storing the representative value, a WL rotation switch 4 for performing a window level operation, and a speed for obtaining a speed primary conversion value ΔWL ′ (described later) from an operation amount ΔWL of the switch 4. A primary conversion unit 5, a speed secondary conversion unit 6 for obtaining a speed secondary conversion value ΔWL ″, an addition unit 12, and a window level determination unit 13 that determines and outputs a window level are provided.
[0033]
The operation sensitivity representative value calculation unit 10 obtains the operation sensitivity representative value ratio (val) based on the histogram data hist (val) according to the following equation (1).
[0034]
[Expression 1]
oratio (val) = max (hist (datai)) * K / hist (val) (1)
i is the number of image data, k is a constant, and val is a pixel value.11~ 211It has a tone of up to -1.
[0035]
In other words, the equation (1) indicates the operation sensitivity with respect to the size of the image data. When the histogram data hist (val) is small, the operation sensitivity representative value ratio (val) is large and changes greatly with a small operation amount. . On the other hand, when hist (val) is large, oratio (val) is small, and fine adjustment is easy. The expression (1) is shown in FIG. 3, and the operation sensitivity becomes extremely large in the vicinity of a pixel value having a low data frequency.
[0036]
Therefore, in order to prevent this, a threshold is provided for the operation sensitivity representative value ratio (val). As a result, for example, a characteristic curve as shown in FIG. 4 is obtained. In this way, an extreme change in operation sensitivity can be prevented. Further, when this is applied to the histogram shown in FIG. 2A, it becomes as shown in FIG. 2B, and it is understood that the operation sensitivity representative value F (x) becomes small at a high frequency of pixel values. it can.
[0037]
Further, in order to reduce the data size of the histogram memory and to make the window sound smooth, the histogram data is averaged at saddle points and is interpolated or used as it is.
[0038]
The speed primary conversion unit 5 calculates and outputs a window level change amount ΔWL ′ in which the rotation speed of the rotary switch is emphasized, and outputs, for example, a square value of the change amount ΔWL per unit time. That is, the following equation (2).
[0039]
ΔWL ′ = ΔWL * | ΔWL | (2)
The speed secondary conversion unit 6 obtains the secondary change amount ΔWL ″ by the following equation (3) from the change amount ΔWL ′ obtained by the speed primary conversion unit 5 and the operation sensitivity representative value function oratio.
ΔWL ″ = ΔWL ′ * oratio (WL-1) ... (3)
Where oratio (WL-1) Is the previous value of the operational sensitivity representative value that changes in time series.
[0040]
The adding unit 12 calculates the previous window level WL-1ΔWL ″ is added to the current window level WL.
[0041]
The window level memory 13 stores and outputs the window level WL obtained by the adding unit 12.
[0042]
Further, the image reading unit 23 of the display unit shown in FIG. 16 supplies an address signal of image data to the image memory 2 and reads the image data corresponding to this address. The window function generator 24 generates a window function as look-up table data based on the window width and window level determined by the window operating device of the present invention. The gray value conversion unit 25 converts the input 12-bit image into a gray value based on a lookup table. The image processing unit 26 performs processing such as converting the image signal after the gray value is converted into a video signal and outputs the video signal to the monitor 27.
[0043]
The operation of the first embodiment of the present invention thus configured will be described below. First, when a medical image is input from the interface 1, the image data is temporarily stored in the image memory 2, sent to the display unit, and displayed on the monitor 27. The image data stored in the image memory 2 is calculated by the histogram calculation unit 3 and stored in the histogram memory 7.
[0044]
Then, either the calculated histogram or the predefined histogram is selected by the histogram selection unit 9 and sent to the operation sensitivity representative value calculation unit 10. The operation sensitivity representative value calculation unit 10 calculates an operation sensitivity representative value according to the above-described equation (1), and stores the calculation result in the operation sensitivity representative value memory 11.
[0045]
Now, when the WL rotation switch 4 is operated, the primary speed change unit 5 obtains the primary change amount ΔWL ′ in accordance with the equation (2) based on the change amount ΔWL of the switch 4. , The secondary conversion amount ΔWL ″ is obtained by the equation (3). “Oratio (WL in the equation (3) at this time”-1The value “)” is a value stored in the operation sensitivity representative value calculation unit 10.
[0046]
Thereafter, the addition shown in the following equation (3) is performed in the adding unit 12 to obtain the current window level WL.
WL = WL-1+ ΔWL ″ (4)
The window level value WL is output from the window level determination unit 13. As a result, the relationship between the operation speed, the operation sensitivity, and the pixel value as shown in FIG. 5 is obtained.
[0047]
In this way, according to the present embodiment, the amount of change in the window level is made variable based on the histogram data of the image (obtained from actual image data or already defined) and the operation amount of the WL rotation switch. Yes. That is, from the histogram data, the change amount is set to be small for a portion with a large pixel value, and the change amount is set to be large for a portion with a small pixel value. Further, when the rotational operation speed of the rotary switch is fast, the window level change amount is set to be large. Accordingly, the window level changes sensitively to the operation of the operator, and the operability is remarkably improved.
[0048]
Next, a second embodiment of the present invention will be described. FIG. 6 is a block diagram showing the configuration of the window operating device according to the second embodiment of the present invention. As shown in the figure, in this embodiment, the designated tissue band center calculation unit 21 and the designated tissue band center memory 22 are arranged in place of the operation sensitivity representative value calculation unit 10 and the operation sensitivity representative value memory 11, respectively. This differs from the first embodiment described above.
[0049]
The designated tissue band center calculation unit 21 calculates the designated tissue band center value i-center according to the following equation (5) from the histogram data hist (val) and the band set by the operator:
i-center = max (hist (val)) (5)
val is within the band specification range (0 to 50 for internal organs, 500 to 1500 for quantities, etc.) and WL from the speed secondary conversion unit-1Based on the above, obtain the sensitivity and set the sensitivity to the speed secondary conversion unit.
Sensitivity = | i-center-WL-1
[0050]
That is, in the second embodiment, when a desired part of interest is designated, the operation sensitivity is set to change with respect to the histogram for the designated part. For example, a histogram as shown in FIG. When given, if the operator wants to observe the internal organ tissue in particular, setting the band for the internal organ tissue sets the amount of change in the operation sensitivity for the histogram for this region. A characteristic curve as shown in b) is obtained. Further, the relationship between the operation speed, the operation sensitivity, and the pixel value as shown in FIG. 8 can be obtained by the processing of the speed secondary conversion unit.
[0051]
Thus, according to the second embodiment, when a desired band is set, the operation sensitivity is set for the histogram for this part, which is particularly useful when the tissue to be observed is determined. .
[0052]
FIG. 9 is a block diagram showing the configuration of the window operating device according to the third embodiment of the present invention. In this embodiment, the operation sensitivity is set based on the current window width WW and the operation amount of the window width. When the WW rotation switch 31 in the figure is rotated, the rotation speed ΔWW is obtained from the change amount WW, and the change of the WW value in which the rotation speed ΔWW is emphasized by the following equation (6) in the speed primary conversion unit 32. A speed ΔWW ′ is determined.
ΔWW ′ = ΔWW * | ΔWW | (6)
[0053]
As a result, the sensitivity characteristic showing a linear function change as shown in FIG. 10 becomes a quadratic function as shown in FIG. 12, and the increase rate of the operation sensitivity increases as the window width WW increases.
[0054]
Thereafter, in the speed secondary conversion unit 33, the previous window width WW stored in the WW memory 35 is displayed.-1Based on the above, the secondary conversion amount ΔWW ″ is obtained by the following equation (7).
ΔWW ″ = ΔWW ′ * (1 + WW-1) / B (7)
B is a constant
[0055]
Thereafter, the addition unit 34 sets the secondary conversion amount ΔWW ″ to the previous window width WW.-1Is added to and the current window width is obtained and output. As a result, it is possible to obtain operation sensitivity corresponding to the rotation speed of the WW rotation switch 31. The characteristic curve at this time is as shown in FIG. 11 for the operation sensitivity shown in FIG. 10, and as shown in FIG. 13 for the operation sensitivity shown in FIG. As is clear from the figure, the operational sensitivity increases as the rotational speed of the rotary switch increases.
[0056]
Next, a fourth embodiment of the present invention will be described. FIG. 14 is a configuration diagram according to the fourth embodiment, which differs from the above-described embodiment of FIG. 9 in that a sensitivity correction unit 36 is provided.
[0057]
The sensitivity correction unit 36 uses the previous window width WW.-1In order to decrease the amount of change with respect to the manipulated value in the vicinity of this good numerical value, for example, when passing a good numerical value such as 1000, 100, etc. Processing as shown is performed.
[0058]
Now, kval the center value of the band you want to reduce the amount of changen(Where n is 1 to M, M is the total number of points set as a clear value), and CnKvalnAssuming that the bandwidth to be reduced before and after is half, the amount of change in the window width is ΔWW ′ ″, and the window width after processing is WW ′ ″, the result is as follows.
[Expression 2]
(Kvaln-Cn) <WW-1<(Kvaln+ Cn)
(Kvaln-Cn) <WW-1  Then
ΔWW ′ ″ = ΔWW ″ * (WW-1−kvaln) / Cn
(Kvaln-Cn)> WW-1  Then
ΔWW ′ ″ = ΔWW ″ * (kvaln-WW-1) / Cn
If (WW ′ ″ / WW ″) <0.5,
ΔWW ′ ″ = ΔWW ″ * 0.5
Other than above; ΔWW ′ ″ = ΔWW ″ (8)
[0059]
Thus, the window width WW '' 'is obtained using the above method, and kvalnThe operational feeling of the window width can be made gentle in the vicinity. For example, kvaln, CnThe
Cn= Int (kvaln 0.5)
Then, data as shown in FIG. 15 (a) is obtained, and the operation sensitivity characteristic is set so that the operation sensitivity becomes small at each point of a good window width as shown in FIG. 15 (b). Operation in the vicinity becomes easy.
[0060]
As described above, the setting of the window width is usually set to a good value such as 100 or 1000, and in the fourth embodiment, the operation sensitivity decreases in the vicinity of the good value of the window width. Operability is improved.
[0061]
FIG. 17 is a block diagram showing a fifth embodiment of the present invention. The configuration of this embodiment is almost the same as that of the first embodiment described above, and a sensitivity correction unit 41 is newly provided.
[0062]
Sensitivity correction unit 41 determines the previous window level WL.-1When ΔWL is added to the current WIN level, for example, when a clear numerical value such as 1000 or 100 is passed, the amount of change with respect to the manipulated variable is decreased in the vicinity of the clear numerical value. In addition, the following processing is performed.
[Equation 3]
(Kvaln-Cn) <WL-1<(Kvaln+ Cn)
(Kvaln-Cn) <WL-1  Then
ΔWL ′ ″ = ΔWL ″ * (WL-1−kvaln) / Cn
(Kvaln-Cn) > WL-1  Then
ΔWL '' '= ΔWL ″ * (kvaln− WL-1) / Cn
If (WL ″ ′ / WL ″) <0.5,
ΔWL ′ ″ = ΔWL ″ * 0.5
Other than above; ΔWL ′ ″ = ΔWL ″ (9)
[0063]
Thus, using the above method, the window level WL '' 'is obtained and kvalnThe operational feeling of the window width can be made gentle in the vicinity.
[0064]
In this way, in the fifth embodiment, since the operation sensitivity is set to decrease at a sharp level such as 100 or 1000, the window level can be easily set to a desired desired value. become able to.
[0065]
Further, a sensitivity correction unit may be provided for the second embodiment shown in FIG. 2, and the window level may be set using the same method as in the fifth embodiment. The configuration in this case is as shown in FIG. 18, and the operation is as described above.
[0066]
Also, since the window level is usually set for tissue parts such as the lungs, internal organs, and bones, the value kval for which the amount of change in the window level should be reducednMay be set in advance to a value such as −1000 (lung), 40 (internal organs), 1000 (bone), or the like. In this way, sensitivity characteristics as shown in FIG. 22 can be obtained.
[0067]
Next, a sixth embodiment of the present invention will be described. In this embodiment, control is performed so that the rate of change gradually increases when the rotary switch for operating the window level is continuously operated. In this embodiment, the processing in the speed primary conversion unit 5 shown in FIGS. 1 and 6 is different. This will be specifically described below with reference to FIG.
[0068]
First, the current value ΔWL of the input data ΔWL of the WL rotation switch 4nAnd previous value ΔWLn-1Difference fromnΔWL is obtained using the following equation (10).
ΔnΔWL = ΔWLn-ΔWLn-1  (10)
[0069]
And this ΔnWhen the value of ΔWL is continuous with a constant width (for example, m times), the output of the speed primary conversion unit 5 is increased as shown in the following equation (11).
[Expression 4]
ΔWLn'= ΔWLn* | ΔWLn| * F (m) (11)
For example, f (m) = mk
k is a constant
[0070]
According to such setting, when the operator continuously operates the WL rotation switch 4, the amount of change with respect to the operation amount increases corresponding to the continuous time. Therefore, for example, even when the window level is changed to 1000, it is not necessary to increase the operation amount of the switch 4, and it is possible to easily reach the predetermined window level only by continuing the gentle rotation.
[0071]
Next, a seventh embodiment of the present invention will be described. FIG. 19 is a configuration diagram according to the seventh embodiment. The configuration of the speed primary conversion unit 5 of the first embodiment shown in FIG. 1 includes a band determination unit 41, a continuity counting unit 42, and a calculation unit 43. This is different from the first embodiment.
[0072]
The band determination unit 41 operates the rotation amount 4 WL of the rotary dial 4 and a predetermined lower limit TLTo determine which band the current band belongs to. And | ΔWL | is continuously lower limit TLThe number m is counted by the continuity counting unit 42 when the number exceeds. Also, the lower limit value TLWhen it becomes below, m is initialized (namely, "1").
[0073]
Then, using this number m, the calculation unit 43 calculates the primary conversion amount ΔWL ′ by the following equation (12).
[Equation 5]
ΔWLn'= ΔWLn* | ΔWLn| * F1(M) ... (12)
For example, f1(M) = mk
k is a constant
Here, FIG. 23 shows a characteristic diagram of ΔWL and a characteristic diagram of ΔWL ′ when the processing according to the equation (12) is not applied, and FIG. 23 shows a characteristic diagram of ΔWL and a characteristic of ΔWL ′ when the processing according to the equation (12) is applied. The figure is shown in FIG.
[0074]
In the seventh embodiment configured as described above, the same effect as that of the sixth embodiment can be obtained, and the operation amount ΔWL is lower than the lower limit value T.LIs controlled so as to increase the primary conversion amount ΔWL ′, while the manipulated variable ΔWL is lower than the lower limit value T.LWhen the value becomes smaller than the initial value, the initial value is restored, so that the window level can be adjusted in accordance with the operation intention of the operator.
[0075]
Next, a modification of the seventh embodiment will be described. FIG. 20 is a block diagram showing the configuration of the modified example, which is different in that the band determining unit 41 shown in FIG.
[0076]
The band zero determination unit 44 determines whether or not the operation amount of the rotary switch 4 is zero. When the operation amount is not zero, the continuity counting unit 42 counts the number m. In other respects, as in the seventh embodiment, the calculation unit 43 calculates the primary conversion amount ΔWL ′ by the following equation (13).
[Formula 6]
ΔWLn'= ΔWLn* | ΔWLn| * F2(M) ... (13)
f2(M) = mk
k is a constant
Here, FIG. 25 shows a characteristic diagram of ΔWL and a characteristic diagram of ΔWL ′ when the processing according to the equation (13) is applied.
[0077]
According to this modification, the seventh embodiment has a predetermined lower limit TLAnd ΔWL is the lower limit TLIn contrast, whether or not the rotation switch 4 is continuously operated is determined by checking whether or not ΔWL is zero. Therefore, this modification is the lower limit T of the seventh embodiment.LIs zero.
[0078]
It is also possible to obtain the primary conversion amount ΔWL by using the method of the seventh embodiment and the method of the modified example together. In this case, the configuration is as shown in FIG. 21, and the calculation formula in the calculation unit 43 is the following formula (14).
[0079]
[Expression 7]
ΔWLn'= ΔWLn* | ΔWLn| * F1(M) * f2(N) (14)
According to such a configuration, a more sensitive operation is possible.
[0080]
【The invention's effect】
As described above, according to the present invention, the amount of change becomes small because the amount of change is small in the region where the window width that needs fine adjustment is small, and conversely, the amount of change becomes large in the region where the window width is large. Therefore, it becomes easy to change the window width.
[Brief description of the drawings]
FIG. 1 is a block diagram showing the configuration of a first embodiment of the present invention.
FIG. 2 is a characteristic diagram showing the relationship between the shading band and the operation sensitivity according to the first embodiment, taking an abdominal image as an example;
FIG. 3 is a characteristic diagram showing a relationship between a pixel value and operation sensitivity.
FIG. 4 is a characteristic diagram showing a relationship between a pixel value and operation sensitivity when a threshold is set.
FIG. 5 is a characteristic diagram showing a relationship among operation speed, operation sensitivity, and pixel value according to the first embodiment.
FIG. 6 is a block diagram showing the configuration of a second embodiment of the present invention.
FIG. 7 is a characteristic diagram showing the relationship between the shading band and the operation sensitivity in the second example, using an abdominal image as an example;
FIG. 8 is a characteristic diagram showing a relationship among operation speed, operation sensitivity, and pixel value according to the second embodiment.
FIG. 9 is a block diagram showing the configuration of a third embodiment of the present invention.
FIG. 10 is a first characteristic diagram showing a relationship between window width and operation sensitivity according to the third embodiment.
FIG. 11 is a first characteristic diagram showing the relationship between window width, operation speed, and operation sensitivity according to the third embodiment.
FIG. 12 is a second characteristic diagram showing the relationship between the window width and the operation sensitivity according to the third embodiment.
FIG. 13 is a second characteristic diagram showing the relationship between the window width, the operation speed, and the operation sensitivity according to the third embodiment.
FIG. 14 is a block diagram showing a configuration of a fourth embodiment of the present invention.
FIG. 15 is a characteristic diagram showing a relationship between window width and sensitivity according to the fourth embodiment.
FIG. 16 is a block diagram illustrating a configuration of a display unit.
FIG. 17 is a block diagram showing a configuration of a fifth embodiment of the present invention.
FIG. 18 is a block diagram showing a configuration of a modified example of the fifth embodiment.
FIG. 19 is a block diagram showing a configuration of a seventh embodiment of the present invention.
FIG. 20 is a block diagram showing a configuration of a modified example of the seventh embodiment.
FIG. 21 is a block diagram showing a configuration of another modified example of the seventh embodiment.
FIG. 22 is a characteristic diagram showing the relationship between the window level and the operation sensitivity when processing according to the fifth embodiment is applied.
FIG. 23 is a characteristic diagram showing changes in ΔWL and ΔWL ′ when the process according to the seventh example is not applied.
FIG. 24 is a characteristic diagram showing changes in ΔWL and ΔWL ′ when a process according to the seventh example is applied.
FIG. 25 is a characteristic diagram showing changes in ΔWL and ΔWL ′ when a process according to a modification of the seventh embodiment is applied.
[Explanation of symbols]
1 Interface 2 Image memory 3 Histogram calculator
4 Window level rotary switch 5,32 Speed primary converter
6,33 Speed secondary conversion unit 7 Histogram memory
8 Predefined histogram memory 9 Histogram selector
10 Operation Sensitivity Representative Value Calculation Unit 11 Operation Sensitivity Representative Value Memory
12, 34 Adder 13 Window level memory
21 Designated Organization Band Center Calculation Unit 22 Designated Organization Band Center Value Memory
31 Window width rotation switch 35 Window width memory
36, 41 Sensitivity correction unit

Claims (3)

入力される画像データのウインドウ幅及びウインドウレベルを操作手段を調整して所望の値に設定するウインドウ操作装置において、
表示されている画像のウインドウ幅の大きさに応じて前記操作手段における操作量に対する前記ウインドウ幅の変化量を変化させる制御手段を有し、
前記制御手段は、前記画像データ全体の濃淡発生頻度分布又は予め設定された濃淡値発生頻度分布に基づき、データが多く発生している濃淡値近傍にて前記操作手段における操作量を緩やかとすべく制御することを特徴とするウインドウ操作装置。
In a window operation device that adjusts the window width and window level of input image data to a desired value by adjusting the operation means,
Control means for changing the amount of change in the window width with respect to the operation amount in the operation means according to the size of the window width of the displayed image;
The control unit is configured to moderate the operation amount in the operation unit in the vicinity of the gray value where a lot of data is generated based on the gray frequency occurrence frequency distribution of the whole image data or a preset gray value occurrence frequency distribution. A window operating device characterized by controlling.
入力される画像データのウインドウ幅及びウインドウレベルを操作手段を調整して所望の値に設定するウインドウ操作装置において、
表示されている画像のウインドウ幅の大きさに応じて前記操作手段における操作量に対する前記ウインドウ幅の変化量を変化させる制御手段と、
前記画像データを格納する画像メモリと、
前記画像データから画素値のヒストグラムを計算する第一計算手段と、
前記ヒストグラムのデータから各画素値に対応する操作感度代表値を求める第二計算手段と、
前記画像データのウインドウレベルを調整するウインドウレベルスイッチと、
前記ウインドウレベルスイッチの操作速度に基づき一次変換量を求める第一変換手段と、
前記一次変換量と前記操作感度代表値に基づき二次変換量を求める第二変換手段と、
前記二次変換量をウインドウレベルに加算して新たなウインドウレベルとする加算手段とを備えることを特徴とするウインドウ操作装置。
In a window operation device that adjusts the window width and window level of input image data to a desired value by adjusting the operation means,
Control means for changing the amount of change in the window width with respect to the operation amount in the operation means in accordance with the size of the window width of the displayed image;
An image memory for storing the image data;
First calculation means for calculating a histogram of pixel values from the image data;
Second calculation means for obtaining an operation sensitivity representative value corresponding to each pixel value from the histogram data;
A window level switch for adjusting a window level of the image data;
First conversion means for obtaining a primary conversion amount based on the operation speed of the window level switch;
Second conversion means for obtaining a secondary conversion amount based on the primary conversion amount and the operation sensitivity representative value;
An addition means for adding the secondary conversion amount to a window level to obtain a new window level.
前記第一計算手段は、前記画像データの所望帯域の画素値についてのみヒストグラムを計算し、
前記第二計算手段は、前記ヒストグラムについて操作感度代表値を求めることを特徴とする請求項2記載のウインドウ操作装置。
The first calculation means calculates a histogram only for pixel values in a desired band of the image data,
3. The window operating device according to claim 2, wherein the second calculating means obtains an operation sensitivity representative value for the histogram.
JP32040394A 1994-12-22 1994-12-22 Window operation device Expired - Fee Related JP3748457B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP32040394A JP3748457B2 (en) 1994-12-22 1994-12-22 Window operation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP32040394A JP3748457B2 (en) 1994-12-22 1994-12-22 Window operation device

Publications (2)

Publication Number Publication Date
JPH08180181A JPH08180181A (en) 1996-07-12
JP3748457B2 true JP3748457B2 (en) 2006-02-22

Family

ID=18121080

Family Applications (1)

Application Number Title Priority Date Filing Date
JP32040394A Expired - Fee Related JP3748457B2 (en) 1994-12-22 1994-12-22 Window operation device

Country Status (1)

Country Link
JP (1) JP3748457B2 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7869642B2 (en) 2006-01-31 2011-01-11 Agfa Inc. Window leveling system and method
JP2009125140A (en) * 2007-11-20 2009-06-11 Konica Minolta Medical & Graphic Inc Image processing apparatus, image processing method, and program
CN102485179B (en) * 2010-12-02 2016-03-09 深圳迈瑞生物医疗电子股份有限公司 Image window width/window position control method, device and medical imaging device
JP7177611B2 (en) * 2017-07-10 2022-11-24 キヤノン株式会社 Image display device, image display method and program
US20190012067A1 (en) * 2017-07-10 2019-01-10 Canon Kabushiki Kaisha Image displaying apparatus, image display method, and storage medium
JP7473313B2 (en) * 2018-10-05 2024-04-23 キヤノンメディカルシステムズ株式会社 Medical image processing device, medical image processing method, and medical image processing program

Also Published As

Publication number Publication date
JPH08180181A (en) 1996-07-12

Similar Documents

Publication Publication Date Title
JP4545291B2 (en) Method and apparatus for controlling the dynamic range of digital diagnostic images
EP0429191A2 (en) Imaging apparatus and methods
JPH0614193A (en) Automatic gradation-scale forming method in digital radiation image
JPS6111615B2 (en)
US20040051710A1 (en) Image display system
US4079417A (en) Digital video window control
US4030119A (en) Video window control
JP4807824B2 (en) Medical diagnostic imaging system
US6215900B1 (en) Image data processing for digital x-ray detectors
JP3748457B2 (en) Window operation device
JPH1014908A (en) X-ray ct apparatus and image display method and apparatus therefor
KR101798083B1 (en) Apparatus and method for processing medical image, and computer readable recording medium related to the method
US6101234A (en) Apparatus and method for displaying computed tomography fluoroscopy images
JP2003198938A (en) Method and apparatus for enhancing contrast of medical diagnostic image that includes foreign object
JP2011120668A (en) Method for optimizing condition for observing local area of digital medical radiation image
JP4314001B2 (en) Image display device
JP4626138B2 (en) Image processing apparatus, image processing method, and program
US20060098006A1 (en) Method of displaying images
JP2012070814A (en) Image processor, program, and x-ray ct apparatus
CN114098769A (en) Medical image display method and module and medical image equipment
JP2005027168A (en) Image processor and image processing method
JP2715841B2 (en) Image processing device
JP7238185B1 (en) Image generating device, medical device, and storage medium
JPH0351891A (en) Image display method
JPH1084504A (en) X-ray image display device

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20040820

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040831

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20041029

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050415

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20050427

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050614

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20050620

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20051128

LAPS Cancellation because of no payment of annual fees