JPH01130287A - Gray level histogram arithmetic unit - Google Patents
Gray level histogram arithmetic unitInfo
- Publication number
- JPH01130287A JPH01130287A JP62288371A JP28837187A JPH01130287A JP H01130287 A JPH01130287 A JP H01130287A JP 62288371 A JP62288371 A JP 62288371A JP 28837187 A JP28837187 A JP 28837187A JP H01130287 A JPH01130287 A JP H01130287A
- Authority
- JP
- Japan
- Prior art keywords
- histogram
- address
- image
- window
- density
- 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.)
- Pending
Links
- 238000006243 chemical reaction Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000010586 diagram Methods 0.000 description 5
- 238000000034 method Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
Description
【発明の詳細な説明】
[発明の目的]
(産業上の利用分野)
本発明は、医用画像等の濃度ヒストグラムを得る濃度ヒ
ストグラム演算装置に関する。DETAILED DESCRIPTION OF THE INVENTION [Object of the Invention] (Industrial Application Field) The present invention relates to a density histogram calculation device for obtaining density histograms of medical images and the like.
(従来の技術)
例えばX線CT装置で搬影されたCT像のような医用画
像をデイスプレィに表示して診断に供する場合、画像を
構成する多数の画素(ピクセル)の各画素濃度を求めて
第9図に示すような濃度CCT値の度合)ヒストグラム
を作成することが行われている。第9図で縦軸は頻度、
横軸は濃度を示しており、求めた各画素濃度を同一面上
に分布させるごとにより、その画像の最大濃度、最小濃
度、平均濃度等を知ることができる。このような結果を
基にして診断に供することができる。(Prior Art) For example, when displaying a medical image such as a CT image transmitted by an X-ray CT device on a display for diagnosis, it is necessary to calculate the density of each pixel of a large number of pixels that make up the image. A histogram of density CCT values as shown in FIG. 9 is created. In Figure 9, the vertical axis is frequency;
The horizontal axis indicates the density, and by distributing the determined pixel densities on the same plane, it is possible to know the maximum density, minimum density, average density, etc. of the image. Diagnosis can be made based on such results.
ところで従来においてこのような濃度ヒストグラムを得
るには、プロセッサを用意しこれによって1画素ずつの
濃度を読込んで、予め設定された濃度グループごとにグ
ループ分けをしてその頻度をカウントすることが行われ
ている。しかしこのようJプロセッサによって1画素ず
つ濃度を求める演算方法では処理に時間がかかり、画像
の任意位置にROI(関心領域)を設定しこのROIを
任意に移動してROIの濃度ヒストグラムを求めるとい
う実際に行われる処理方法では実用に耐えられるだけの
演算速度が得られない場合が多い。Conventionally, in order to obtain such a density histogram, a processor was prepared, which read the density of each pixel, divided into preset density groups, and counted the frequency. ing. However, this calculation method of calculating the density of each pixel using the J processor takes a long time, and in practice, it involves setting an ROI (region of interest) at an arbitrary position in the image and moving this ROI arbitrarily to obtain the density histogram of the ROI. In many cases, the processing method used in this paper cannot achieve a calculation speed sufficient for practical use.
(発明が解決しようとする問題点)
このように従来の濃度ヒストグラム演算装置では処理速
度が遅いので、リアルタイムで濃度ヒストグラムを得る
のは不可能であるという問題がある。(Problems to be Solved by the Invention) As described above, since the processing speed of the conventional density histogram calculation device is slow, there is a problem that it is impossible to obtain a density histogram in real time.
本発明は以上のような問題に対処してなされたもので、
リアルタイムで濃度ヒストグラムを得ることができる濃
度ヒストグラム演算装置を提供することを目的とするも
のである。The present invention has been made in response to the above-mentioned problems.
It is an object of the present invention to provide a density histogram calculation device that can obtain a density histogram in real time.
[発明の構成]
(問題点を解決するための手段)
上記目的を達成するために本発明は、多数の画素のデー
タをウィンド変換をうまく利用することによって頻度を
グレイコードとして変換し、このグレイコードをアドレ
スとしてヒストグラムメモリに書き込むと共にこのヒス
トグラムのアドレスごとのデータを読み出してインクリ
メンタに加え、このインクリメンタの出力をヒストグラ
ムメモリの同一アドレスに再杏込みすることを特徴とす
るものである。[Structure of the Invention] (Means for Solving the Problems) In order to achieve the above object, the present invention converts the frequency of data of a large number of pixels into a Gray code by making good use of window transformation, and converts the frequency into a Gray code. The code is written as an address in the histogram memory, the data for each address of this histogram is read out and added to an incrementer, and the output of this incrementer is re-embedded at the same address in the histogram memory.
(作 用)
ウィンド変換の出力でおるグレイコードを利用すること
によりこれがアドレスとしてヒストグラムメモリに書き
込まれる。ヒストグラムメモリの各アドレスのデータは
インクリメンタに加えられ、インクリメンタの各出ノ〕
はヒストグラムメモリの対応した同一アドレスに再書込
まれる。これによってグレイコードに対応した各画素濃
度ごとの頻度がヒストグラムメモリによってカウントさ
れるので、濃度ヒストグラムが得られる。このような演
Wffi理はウィンド変換の出力を利用して行われるの
で、ウィンド変換速度と比例した速度で行われる。従っ
て高速で演算処理を行うことができ、リアルタイムで濃
度ヒストグラムを得ることができる。(Operation) By using the Gray code output from the window transformation, this is written into the histogram memory as an address. The data at each address in the histogram memory is added to the incrementer, and the data at each address in the incrementer is added to the incrementer.
is rewritten to the same corresponding address in the histogram memory. As a result, the frequency of each pixel density corresponding to the Gray code is counted by the histogram memory, so that a density histogram is obtained. Since such a Wffi operation is performed using the output of the window transform, it is performed at a speed proportional to the window transform speed. Therefore, calculation processing can be performed at high speed, and a density histogram can be obtained in real time.
(実施例)
第1図は本発明の痩度ヒストグラム演算装置の実施例を
示すブロック図で、ウィンドプロセッサ1は画像メモリ
4に格納されている原画像(例えばCT像)を最適のウ
ィンド変換条件でウィンド変換してこの結果をフレーム
メモリ3に格納するためのものである。画像メモリ4に
はCT値のデータが例えば12ビツトで格納されており
、このCT値データはウィンド変換されて例えば8ビツ
トのデータに圧縮されてフレームメモリ3に格納される
。ウィンドテーブル2にはウィンド変換条件に対応して
、表示されるCT値の範囲を決めるウィンド幅(W、W
)及びその範囲の中心のCT値を決めるウィンドレベル
(W、L)から成るパラメータが格納されている。(Embodiment) FIG. 1 is a block diagram showing an embodiment of the leanness histogram calculation device of the present invention, in which a wind processor 1 converts an original image (for example, a CT image) stored in an image memory 4 under optimal window conversion conditions. This is for performing window conversion and storing the result in the frame memory 3. The image memory 4 stores CT value data in, for example, 12 bits, and this CT value data is window-converted and compressed into, for example, 8-bit data, which is then stored in the frame memory 3. Wind table 2 contains window widths (W, W) that determine the range of displayed CT values in accordance with the window conversion conditions.
) and the wind level (W, L) that determines the CT value at the center of the range.
フレームメモリ3には画像データが表示イメージで格・
納されており、この画像データはD/Aコンバータ8に
よってD/A変換された後CRTデイスプレィ9に画像
として表示される。すなわちX線CT装置で測定された
被検体のCT、ifiの分イ「が、観察したい部位に応
じてこの部位が最適のコントラストで表示されるように
設定されたウィンド変換条件でウィンド変換され濃淡画
像がデイスプレィ9に表示されることになる。Image data is stored in frame memory 3 as a display image.
This image data is D/A converted by a D/A converter 8 and then displayed as an image on a CRT display 9. In other words, the CT of the subject measured by an X-ray CT device, i.e., the portion of ifi, is converted into shading using window conversion conditions set so that the area to be observed is displayed with the optimum contrast. The image will be displayed on the display 9.
ビットプレーン7はビットマツプデイスプレィの機能を
有しており、グラフィックコントローラ6によって高速
にグラフィック表示が行われる。The bit plane 7 has a bitmap display function, and graphics controller 6 performs graphic display at high speed.
M P LJ 5は専用のプロセッサから成っており、
装置全体の制御動作を行う。ROI入力手段12は71
クス、トラッカーボール等から構成され、デイスプレィ
9に画像の表示が行われているとき、任意位置にROI
を設定するとその画像上にはグラフィックコントローラ
6の機能によりその設定に基づいた四辺形1円形、楕円
形等の任意の形状のROIが高速に表示される。またこ
のROIは任意位置に移動可能に構成されると共に拡大
、縮小も自在となるように構成されている。10.、.
11はデータバスである。M P LJ 5 consists of a dedicated processor,
Performs control operations for the entire device. ROI input means 12 is 71
When an image is displayed on the display 9, ROI can be placed at any position.
When set, an ROI of arbitrary shape such as a quadrilateral, a circle, an ellipse, etc. based on the setting is displayed on the image at high speed by the function of the graphic controller 6. Further, this ROI is configured to be movable to any position and also configured to be enlarged and reduced. 10. ,.
11 is a data bus.
ウィンドプロセッサ1はW、W及びW、Lの2つのパラ
メータから決定されるウィンド変換条件によって、CT
値から表示のためのグレイコードを決定して出力する。The wind processor 1 performs CT according to the window conversion conditions determined from two parameters W, W and W, L.
Determines the gray code for display from the value and outputs it.
第2図はこのようなウィンド変換の様子を示すもので、
縦軸はグレイコードで例えば2B=256段階からから
成り、横軸はCT値で例えば2 +2 = 4096段
階から成っている。Figure 2 shows this kind of window conversion.
The vertical axis represents the Gray code, which is made up of, for example, 2B=256 steps, and the horizontal axis is the CT value, which is made up of, for example, 2+2=4096 steps.
12ビツトのCT値は8ビツトのグレイコードに対応す
るように圧縮され、このために複数のグループ分けされ
各グループは1つのCT値D1゜D2 、D3・・・で
代表されるような圧縮方法がとられる。例えばCT値D
1で代表されるグループには成る範囲のCT値が含まれ
ており、これらのCT値はすべてグレイコード1に対応
される。このようなウィンド変換条件はウィンドテーブ
ル2に格納しておくことができる。 ′
ウィンドプロセッサ1から出力されたグレイコードはヒ
ストグラムメモリ13にアドレスとして書込まれる。ヒ
ストグラムメモリ13は第3図に示すようにOから25
5に至る256段階のアドレスを有しており、各アドレ
スはグレイコードに対応している。例えばウィンドプロ
セッサ1から出力されたグレイコード2はヒストグラム
メモリ13のアドレス2として書込まれる。ヒストグラ
ムメモリ13の内容は初期状態ですべてOにリセットさ
れている。このヒストグラムメモリ13の各アドレスの
データはインクリメンタ14に加えられ、このインクリ
メンタ14はこれをインクリメント(1を加える)とし
てヒストグラムメモリ13の同一アドレスに再書込みを
行う。The 12-bit CT value is compressed to correspond to the 8-bit Gray code, and for this reason, it is divided into multiple groups, and each group is represented by one CT value D1, D2, D3, etc. is taken. For example, CT value D
The group represented by 1 includes a range of CT values, all of which correspond to Gray code 1. Such window conversion conditions can be stored in the window table 2. ' The Gray code output from the window processor 1 is written into the histogram memory 13 as an address. As shown in FIG.
It has 256 levels of addresses up to 5, and each address corresponds to a Gray code. For example, Gray code 2 output from the window processor 1 is written as address 2 in the histogram memory 13. The contents of the histogram memory 13 are all reset to O in the initial state. The data at each address in the histogram memory 13 is added to the incrementer 14, and the incrementer 14 increments the data (adds 1) and rewrites the data to the same address in the histogram memory 13.
次に本実施例の作用を説明する。Next, the operation of this embodiment will be explained.
ヒストグラムメモリ13の内容がすべてOにリレットさ
れている状態で、ウィンドプロセッサ1から例えばグレ
イコード2がヒストグラムメモリ13にアドレス2とし
て書込まれたとすると、このアドレス2のデータO(O
にリセットされている)がインクリメンタ14に加えら
れる。インクリメンタ14はこのデータOをインクリメ
ントして(すなわち1を加えて)この結果の1をヒスト
グラムメモリ13の同一アドレスでおるアドレス2に再
書込みを行う。これによってアドレス2にはデータ1が
格納されることになる。第4図(a)はこの様子を示し
ている。次にウィンドプロセッサ1から例えばグレイコ
ード4がヒストグラムメモリ13のアドレス4として書
込まれたとすると、同様な原理ににつで第4図(b)、
のようにアドレス4にはデータ1が格納されることにな
る。If, for example, Gray code 2 is written from the window processor 1 to the histogram memory 13 as address 2 while the contents of the histogram memory 13 are all retted to O, then the data O(O
) is added to the incrementer 14. The incrementer 14 increments this data O (that is, adds 1) and rewrites the resulting 1 to address 2 of the histogram memory 13, which is the same address. As a result, data 1 will be stored at address 2. FIG. 4(a) shows this situation. Next, suppose that, for example, Gray code 4 is written as address 4 in the histogram memory 13 from the window processor 1, the same principle can be used as shown in FIG. 4(b).
Data 1 will be stored at address 4 as shown in FIG.
次にグレイコード2が再びヒストグラムメモリ13のア
ドレス2に書込まれたとすると、このアドレス2からは
既に格納されているデータ1がインクリメンタ14に加
えられる。するとインクリメンタ14はインクリメント
したデータ2をアドレス2に再書込みする。これによっ
てアドレス2には第4図(C)のようにデータ2が格納
されることになる。Next, if Gray code 2 is written again to address 2 of the histogram memory 13, the already stored data 1 is added to the incrementer 14 from this address 2. Then, the incrementer 14 rewrites the incremented data 2 to address 2. As a result, data 2 is stored at address 2 as shown in FIG. 4(C).
このようにしてすべてのアドレスに対してグレイコード
の書き込みが終了すると、MPU5がヒストグラムメモ
リ13のアドレス1から順々にアドレス255に至るま
でそのデータを読み取ることにより各アドレスに対応し
た各グレイコードの数(頻度)をカウントすることがで
きる。ざらに第2図の関係に基づいて各グレイコードに
対応したCT値Di 、D2 、D3・・・に置換えて
この頻度を示せば第5図のような濃度ヒストグラムを得
ることができる。これらの動作はウィンドブ[]セッサ
1及σグラフィックコントローラ6の機能を利用するこ
とにより容易に行うことができる。ウィンド変換処理は
通常10乃至30回/秒の処理速度が可能なので、この
逆数の1710乃至1730秒で濃度ヒストグラムの演
n処理を行うことができることになる。従って高速で演
l’Jffi理を行うことができるので、リアルタイム
で濃度ヒストグラムを得ることができる。After writing the Gray code to all addresses in this way, the MPU 5 reads the data from address 1 to address 255 in the histogram memory 13, thereby writing each Gray code corresponding to each address. The number (frequency) can be counted. Roughly based on the relationship shown in FIG. 2, a density histogram as shown in FIG. 5 can be obtained by replacing the CT values Di, D2, D3, . . . corresponding to each gray code and indicating the frequency. These operations can be easily performed by using the functions of the window processor 1 and the σ graphics controller 6. Since window conversion processing can normally be performed at a processing speed of 10 to 30 times/second, the density histogram calculation process can be performed in the reciprocal of 1710 to 1730 seconds. Therefore, since the calculation can be performed at high speed, the density histogram can be obtained in real time.
また、デイスプレィ9に画像の表示が行われているとき
、ROI入力手段12によってその画像上の任意位置に
ROI 15を設定することにより、第6図のようにそ
の画像上のROI 15内にRO115のみの濃度ヒス
トグラムを表示することができる。また第7図のように
ROI 15外にこの濃度ヒストグラムを表示させるこ
ともできる。Further, when an image is displayed on the display 9, by setting the ROI 15 at an arbitrary position on the image using the ROI input means 12, an ROI 115 is placed within the ROI 15 on the image as shown in FIG. Only the concentration histogram can be displayed. Further, this density histogram can also be displayed outside the ROI 15 as shown in FIG.
ざらにROI 15を移動したり、拡大、縮小させるこ
とによりその都度濃度ヒストグラムの表示を行わせるこ
とができる。ざらにまた第8図のように画像上の任意の
特定領域16のみを対象としてこの濃度ヒストグラムの
表示を行わせることも可能である。By roughly moving, enlarging, or reducing the ROI 15, the density histogram can be displayed each time. Furthermore, as shown in FIG. 8, it is also possible to display the density histogram only for an arbitrary specific area 16 on the image.
画像上の必要な位置のみにROIを設定してこのROI
内の濃度ヒストグラムを求める方法は実際上堰も要望が
高いので、このような処理が早い演算速度で実現できる
ことは極めて有効である。Set the ROI only at the required position on the image and use this ROI.
Since a method for determining the concentration histogram within the range is in high demand in practice, it is extremely effective to be able to perform such processing at a high calculation speed.
本分実施例中で示した画像の圧縮方法は一例を示したも
のであり、必要に応じてその内容は任意に変更すること
ができる。The image compression method shown in the main embodiment is merely an example, and its contents can be arbitrarily changed as necessary.
[発明の効果]
以上述べたように本発明によれば、ウィンド変換の出力
を利用して濃度ヒストグラムを求めるようにしたので、
高速な演算処理が可能となりリアルタイムで濃度ヒスト
グラムを得ることができるt[Effects of the Invention] As described above, according to the present invention, since the density histogram is obtained using the output of the window transformation,
It enables high-speed calculation processing and allows concentration histograms to be obtained in real time.
第1図は本発明の濃度ヒス1〜グラム演痺装置の実施例
を示すブロック図、第2図は本実施例装置によるウィン
ド変換方法の説明図、第3図は本実施例装置のヒストグ
ラムメモリの配置図、第4図(a)乃至(C)は本発明
の濃度ヒストグラム演算の原理の説明図、第5図は本実
施例にJ:つて1qられた濃度ヒストグラム、第6図乃
至第8図は本発明の他の実施例の濃度ヒストグラムの表
示例、第9図は濃度ヒストグラムの説明図である。
1・・・ウィンドプロセッサ、
2・・・ウィンドテーブル、3・・・フレームメモリ、
4・・・画像メモリ、12・・・ROI入力手段、13
・・・ヒストグラムメモリ、
14・・・インクリメンタ。
代理人 弁理士 則 近 憲 化量
近 藤 猛第3図
(C)
弔8図 カ9図FIG. 1 is a block diagram showing an embodiment of the concentration hiss 1 to gram numbing device of the present invention, FIG. 2 is an explanatory diagram of a window conversion method by the device of this embodiment, and FIG. 3 is a histogram memory of the device of this embodiment. 4(a) to 8(C) are explanatory diagrams of the principle of density histogram calculation according to the present invention. FIG. The figure is a display example of a density histogram according to another embodiment of the present invention, and FIG. 9 is an explanatory diagram of the density histogram. 1... Wind processor, 2... Wind table, 3... Frame memory,
4... Image memory, 12... ROI input means, 13
...Histogram memory, 14...Incrementer. Agent Patent Attorney Rules
Takeshi Kondo Figure 3 (C) Funeral Figure 8 Figure C9
Claims (2)
ことにより画像の濃度ヒストグラムを得る濃度ヒストグ
ラム演算装置において、多数の画素のデータをグレイコ
ードに変換するウィンド変換手段と、ウインド変換手段
から出力されたグレイコードをアドレスとして入力し各
アドレスごとの頻度をカウントするヒストグラム演算手
段と、ヒストグラム演算手段から前記アドレスに対応し
たデータを読み出しこのデータをインクリメントして前
記ヒストグラム演算手段の同一アドレスにカウントすべ
きデータとして再書込みを行うインクリメンタとを備え
たことを特徴とする濃度ヒストグラム演算装置。(1) In a density histogram calculation device that obtains a density histogram of an image by determining the density of each pixel of a large number of pixels constituting the image, a window conversion means that converts data of a large number of pixels into a gray code, and a window conversion means. a histogram calculation means that inputs the output gray code as an address and counts the frequency of each address; and a histogram calculation means that reads data corresponding to the address from the histogram calculation means, increments this data, and counts it to the same address of the histogram calculation means. 1. A density histogram calculation device comprising: an incrementer that rewrites data as data to be processed.
段を備え、このROIのみをヒストグラム演算の対象と
する特許請求の範囲第1項記載の濃度ヒストグラム演算
装置。(2) The density histogram calculation device according to claim 1, further comprising ROI input means for setting an ROI at an arbitrary position in an image, and in which only this ROI is subjected to histogram calculation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP62288371A JPH01130287A (en) | 1987-11-17 | 1987-11-17 | Gray level histogram arithmetic unit |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP62288371A JPH01130287A (en) | 1987-11-17 | 1987-11-17 | Gray level histogram arithmetic unit |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH01130287A true JPH01130287A (en) | 1989-05-23 |
Family
ID=17729332
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP62288371A Pending JPH01130287A (en) | 1987-11-17 | 1987-11-17 | Gray level histogram arithmetic unit |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH01130287A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009302043A (en) * | 2008-06-13 | 2009-12-24 | Furutech Co Ltd | Power plug |
JP2012198139A (en) * | 2011-03-22 | 2012-10-18 | Olympus Corp | Image processing program, image processing device, measurement analysis device and image processing method |
-
1987
- 1987-11-17 JP JP62288371A patent/JPH01130287A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009302043A (en) * | 2008-06-13 | 2009-12-24 | Furutech Co Ltd | Power plug |
JP2012198139A (en) * | 2011-03-22 | 2012-10-18 | Olympus Corp | Image processing program, image processing device, measurement analysis device and image processing method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE69028195T2 (en) | IMAGE TREATMENT UNIT | |
JPH01130287A (en) | Gray level histogram arithmetic unit | |
JPH07131653A (en) | Image gradation compression display circuit | |
JPH063607B2 (en) | Image processing device | |
KR920003702B1 (en) | Image processing method and apparatus for the search of volume and its center of 3-d image | |
JPS61172541A (en) | Double window display apparatus | |
JPH06178115A (en) | Image processor | |
JP3455027B2 (en) | X-ray image display device | |
CN112767508A (en) | Pseudo color display method, system, device, and medium for respiratory gated image | |
JPS62117080A (en) | Picture processing device | |
JPS6364217B2 (en) | ||
JPH02210393A (en) | Image display device | |
JP2814769B2 (en) | Image processing device | |
JPH032534B2 (en) | ||
JPH0333328B2 (en) | ||
JP2002330948A (en) | Method and device for controlling image processing | |
JPH01137375A (en) | Picture processor | |
JPH0546744A (en) | Image processor | |
JP3015227B2 (en) | Image processing device | |
JP3077265B2 (en) | Shading image processing device | |
JP2601059B2 (en) | Image processing device | |
JPH0410669B2 (en) | ||
JPH06100900B2 (en) | Method and apparatus for high-speed color change display of color image | |
JP2021175111A (en) | Image processing device, image processing method, and control program | |
JP2696997B2 (en) | Digital image data converter |