JPH07129751A - Medical picture processor - Google Patents
Medical picture processorInfo
- Publication number
- JPH07129751A JPH07129751A JP5292365A JP29236593A JPH07129751A JP H07129751 A JPH07129751 A JP H07129751A JP 5292365 A JP5292365 A JP 5292365A JP 29236593 A JP29236593 A JP 29236593A JP H07129751 A JPH07129751 A JP H07129751A
- Authority
- JP
- Japan
- Prior art keywords
- circuit
- region
- image
- interest
- area
- 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
- 238000000605 extraction Methods 0.000 claims abstract description 30
- 239000000284 extract Substances 0.000 abstract description 4
- 210000001835 viscera Anatomy 0.000 abstract 2
- 230000006870 function Effects 0.000 description 18
- 238000010586 diagram Methods 0.000 description 16
- 210000000056 organ Anatomy 0.000 description 15
- 238000002595 magnetic resonance imaging Methods 0.000 description 8
- 238000000034 method Methods 0.000 description 8
- 238000005259 measurement Methods 0.000 description 6
- 206010028980 Neoplasm Diseases 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000012447 hatching Effects 0.000 description 2
- 239000000470 constituent Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
Description
【０００１】[0001]
【産業上の利用分野】本発明は、Ｘ線ＣＴ画像やＭＲＩ
（磁気共鳴イメージング）画像等の医用画像について各
種の画像処理を施す医用画像処理装置に関し、特に表示
された医用画像から臓器抽出や任意の領域抽出を行う際
にその抽出のための画素値のしきい値を最適に選択する
ことができる医用画像処理装置に関する。The present invention relates to X-ray CT images and MRI.
(Magnetic resonance imaging) A medical image processing apparatus for performing various image processing on a medical image such as an image, particularly when performing organ extraction or arbitrary region extraction from the displayed medical image The present invention relates to a medical image processing apparatus capable of optimally selecting a threshold value.
【０００２】[0002]
【従来の技術】従来の医用画像診断装置は、複数の医用
画像のデータを記録しておく記憶装置と、この記憶装置
から読み出した画像データを表示するための処理を行う
画像表示回路と、この画像表示回路から出力された画像
データを画像として表示する表示装置と、この表示装置
に表示された画像上で関心領域を決定する関心領域決定
回路と、この関心領域決定回路で決定された関心領域内
の画素値に対してかける重み関数を生成するフィルタ発
生回路と、このフィルタ発生回路で生成された重み関数
を上記関心領域に対してかけ合わせるフィルタ回路と、
上記関心領域内での画素値ごとの度数を計算してグラフ
を作成して上記画像表示回路へ送出するグラフ作成回路
と、上記各構成要素に操作指令を入力する入力装置とを
備えて成っていた。そして、表示画像から特定の臓器や
領域の抽出を行うには、画像の画素値をしきい値で２値
化して抽出していた。このとき、上記しきい値の決定
は、表示のレベル、ウインドウ幅を絞り込んで表示され
る領域を観察しながら行い、その決定したしきい値は、
キーボードなどの入力装置から画素値の値を直接入力し
ていた。また、関心領域決定回路で決定された関心領域
内での画素値ごとの度数を計算してグラフを作成し、こ
のグラフ化されたヒストグラムからしきい値を決定し
て、特定の臓器や領域の抽出を行う方法もあった。2. Description of the Related Art A conventional medical image diagnostic apparatus includes a storage device for recording data of a plurality of medical images, an image display circuit for performing processing for displaying image data read from the storage device, A display device that displays the image data output from the image display circuit as an image, a region of interest determination circuit that determines a region of interest on the image displayed on the display device, and a region of interest determined by the region of interest determination circuit A filter generating circuit for generating a weighting function to be applied to pixel values in the inside, and a filter circuit for multiplying the weighting function generated by the filter generating circuit to the region of interest,
It comprises a graph creating circuit for calculating the frequency for each pixel value in the region of interest to create a graph and sending it to the image display circuit, and an input device for inputting an operation command to each component. It was Then, in order to extract a specific organ or region from the display image, the pixel value of the image is binarized with a threshold value and then extracted. At this time, the threshold value is determined while observing the display level and the area displayed by narrowing the window width, and the determined threshold value is
The pixel value was directly input from an input device such as a keyboard. In addition, a graph is created by calculating the frequency for each pixel value in the region of interest determined by the region of interest determination circuit, and the threshold value is determined from this graphed histogram to determine the specific organ or region. There was also a way to do the extraction.
【０００３】[0003]
【発明が解決しようとする課題】しかし、このような従
来の医用画像処理装置における表示画像からの特定領域
の抽出においては、特にしきい値の値のとり方が重要で
あった。とりわけＭＲＩ画像については、計測条件、静
磁場強度などによる画像のコントラストのつき方や画素
値の変動が大きく、しきい値も画像によって変動させる
必要があった。従来の医用画像処理装置では、表示画像
からのしきい値の決定は半自動化されているが、最適な
領域抽出を行うためには、表示画像を観察しながら試行
錯誤的にしきい値を上げ下げして微調整を行う必要があ
った。また、しきい値により抽出される画像上の領域の
面積も大きく変化するが、従来は、上記しきい値を決定
する際には、そのしきい値と抽出される画像上の領域の
面積との関係については考慮されてはいなかった。従っ
て、最適なしきい値を容易に選択して決定することはで
きなかった。However, in extracting a specific area from a display image in such a conventional medical image processing apparatus, how to set a threshold value is particularly important. In particular, with regard to MRI images, there are large variations in the contrast of images and pixel values due to measurement conditions, static magnetic field strength, etc., and it was necessary to vary the threshold value depending on the image. In the conventional medical image processing apparatus, the determination of the threshold value from the display image is semi-automatic, but in order to perform the optimum area extraction, the threshold value is raised and lowered by trial and error while observing the display image. It was necessary to make fine adjustments. Further, although the area of the region on the image extracted by the threshold changes greatly, conventionally, when determining the threshold, the threshold and the area of the region on the extracted image are Was not considered. Therefore, the optimum threshold value cannot be easily selected and determined.
【０００４】そこで、本発明は、このような問題点に対
処し、表示された医用画像から臓器抽出や任意の領域抽
出を行う際にその抽出のための画素値のしきい値を最適
に選択することができる医用画像処理装置を提供するこ
とを目的とする。Therefore, the present invention addresses such a problem and optimally selects the threshold value of the pixel value for the extraction of organs or arbitrary regions from the displayed medical image. It is an object of the present invention to provide a medical image processing apparatus that can be used.
【０００５】[0005]
【課題を解決するための手段】上記目的を達成するため
に、本発明による医用画像処理装置は、複数の医用画像
のデータを記録しておく記憶装置と、この記憶装置から
読み出した画像データを表示するための処理を行う画像
表示回路と、この画像表示回路から出力された画像デー
タを画像として表示する表示装置と、この表示装置に表
示された画像上で関心領域を決定する関心領域決定回路
と、この関心領域決定回路で決定された関心領域内の画
素値に対してかける重み関数を生成するフィルタ発生回
路と、このフィルタ発生回路で生成された重み関数を上
記関心領域に対してかけ合わせるフィルタ回路と、上記
関心領域内での画素値ごとの度数を計算してグラフを作
成して上記画像表示回路へ送出するグラフ作成回路と、
上記各構成要素に操作指令を入力する入力装置とを備え
て成る医用画像処理装置において、上記グラフ作成回路
からのグラフのデータを入力し指定された画素値の上限
と下限のしきい値内で画素度数の面積を計算する画素度
数面積計算回路と、上記入力装置で入力されたしきい値
の上限及び下限をグラフ上に標示線として標示するしき
い値選択回路と、このしきい値選択回路で選択された範
囲の画素値のみの領域を抽出して表示する領域抽出回路
とを設けたものである。In order to achieve the above object, a medical image processing apparatus according to the present invention stores a storage device in which data of a plurality of medical images is recorded and image data read from the storage device. An image display circuit that performs a process for displaying, a display device that displays the image data output from the image display circuit as an image, and a region-of-interest determination circuit that determines a region of interest on the image displayed on the display device. And a filter generating circuit for generating a weighting function for multiplying the pixel values in the ROI determined by the ROI, and the weighting function generated by the filter generating circuit for multiplying the ROI. A filter circuit, a graph creation circuit for calculating a frequency for each pixel value in the region of interest to create a graph and sending it to the image display circuit,
In a medical image processing apparatus comprising an input device for inputting an operation command to each of the constituent elements, within the upper and lower threshold values of the pixel value specified by inputting the data of the graph from the graph creating circuit A pixel frequency area calculation circuit for calculating the area of the pixel frequency, a threshold value selection circuit for indicating the upper and lower limits of the threshold value input by the input device as a marking line on the graph, and this threshold value selection circuit. An area extraction circuit for extracting and displaying an area of only the pixel value in the range selected in step 1 is provided.
【０００６】[0006]
【作用】このように構成された医用画像処理装置は、画
素度数面積計算回路でグラフ作成回路から出力されたグ
ラフのデータを入力して指定された画素値の上限と下限
のしきい値内で画素度数の面積を計算し、しきい値選択
回路により入力装置で入力されたしきい値の上限及び下
限をグラフ上に標示線として表示し、領域抽出回路で上
記しきい値選択回路により選択された範囲の画素値のみ
の領域を抽出して表示するように動作する。これによ
り、表示された医用画像から臓器抽出や任意の領域抽出
を行う際にその抽出のための画素値のしきい値を最適に
選択して、特定領域を抽出することができる。The medical image processing apparatus configured as described above receives the data of the graph output from the graph creation circuit in the pixel frequency area calculation circuit, and within the upper and lower threshold values of the specified pixel value. The area of the pixel frequency is calculated, the upper and lower limits of the threshold value input by the input device are displayed as a marking line on the graph by the threshold selection circuit, and selected by the threshold selection circuit by the area extraction circuit. It operates so as to extract and display an area having only pixel values in the selected range. Thus, when performing organ extraction or arbitrary area extraction from the displayed medical image, the threshold value of the pixel value for the extraction can be optimally selected to extract the specific area.
【０００７】[0007]
【実施例】以下、本発明の実施例を添付図面に基づいて
詳細に説明する。図１は本発明による医用画像処理装置
の実施例を示すブロック図である。この医用画像処理装
置は、予め記録された医用画像について各種の画像処理
を施すもので、図１に示すように、記憶装置１と、画像
表示回路２と、表示装置３と、関心領域決定回路４と、
フィルタ発生回路５と、フィルタ回路６と、グラフ作成
回路７と、入力装置８とを備え、さらに画素度数面積計
算回路９と、しきい値選択回路１０と、領域抽出回路１
１とを設けて成る。Embodiments of the present invention will now be described in detail with reference to the accompanying drawings. FIG. 1 is a block diagram showing an embodiment of a medical image processing apparatus according to the present invention. This medical image processing apparatus performs various kinds of image processing on medical images recorded in advance, and as shown in FIG. 1, a storage device 1, an image display circuit 2, a display device 3, and a region of interest determination circuit. 4 and
The filter generation circuit 5, the filter circuit 6, the graph creation circuit 7, and the input device 8 are provided, and the pixel frequency area calculation circuit 9, the threshold value selection circuit 10, and the region extraction circuit 1 are provided.
1 and are provided.
【０００８】上記記憶装置１は、図示外のＸ線ＣＴ装置
やＭＲＩ装置などの医用画像診断装置で計測された複数
の医用画像のデータを記録しておくもので、例えば磁気
テープ、磁気ディスク、光ディスクなどから成る。ま
た、画像表示回路２は、上記記憶装置１から読み出した
画像データを表示するための処理を行うもので、その内
部には、文字情報や記号等を表示するためのキャラクタ
メモリ１２と、後述のヒストグラム等の図形情報を表示
するためのグラフィックメモリ１３と、画像を表示する
ためのディスプレイイメージメモリ１４とを有してい
る。さらに、表示装置３は、上記画像表示回路２から出
力された画像データを画像として表示するもので、例え
ばテレビモニタから成る。なお、図１において、画像メ
モリ１５は、上記記憶装置１から読み出した画像データ
を表示装置３に表示するために一時的に記憶するもので
ある。The storage device 1 stores data of a plurality of medical images measured by a medical image diagnostic device such as an X-ray CT device or an MRI device (not shown). For example, a magnetic tape, a magnetic disk, It consists of an optical disk. Further, the image display circuit 2 performs a process for displaying the image data read from the storage device 1, and inside thereof, a character memory 12 for displaying character information, symbols and the like, and a later-described It has a graphic memory 13 for displaying graphic information such as a histogram and a display image memory 14 for displaying an image. Further, the display device 3 displays the image data output from the image display circuit 2 as an image, and includes, for example, a television monitor. In FIG. 1, the image memory 15 temporarily stores the image data read from the storage device 1 for displaying on the display device 3.
【０００９】関心領域決定回路４は、上記表示装置３に
表示された画像上で関心領域を決定するもので、後述の
入力装置８を操作者が操作して関心領域の範囲が入力さ
れるようになっている。フィルタ発生回路５は、上記関
心領域決定回路４で決定された関心領域内の画素値に対
してかける重み関数を生成するもので、上記関心領域内
の画素値について線形あるいは非線形の変換を行う関数
を生成して、変化の小さいところは大きくなるように、
また不要な範囲は狭帯域に圧縮するようになっている。
また、フィルタ回路６は、上記フィルタ発生回路５で生
成された重み関数を前記関心領域決定回路４で決定され
た関心領域に対してかけ合わせるものである。The region-of-interest determining circuit 4 determines the region of interest on the image displayed on the display device 3, so that the operator operates the input device 8 described later to input the range of the region of interest. It has become. The filter generating circuit 5 generates a weighting function to be applied to the pixel value in the ROI determined by the ROI determining circuit 4, and is a function for performing linear or non-linear conversion on the pixel value in the ROI. , So that small changes are large,
The unnecessary range is compressed to a narrow band.
The filter circuit 6 multiplies the weighting function generated by the filter generating circuit 5 with the region of interest determined by the region of interest determining circuit 4.
【００１０】グラフ作成回路７は、前記関心領域決定回
路４で決定された関心領域内での画素値ごとの度数を計
算してヒストグラムを作成して前記画像表示回路２へ送
出するものである。また、入力装置８は、上記各構成要
素に対し操作指令を入力するもので、例えばマウス、ト
ラックボール、ジョイスティック等から成る。The graph creating circuit 7 calculates the frequency for each pixel value in the region of interest determined by the region of interest determining circuit 4, creates a histogram, and sends it to the image display circuit 2. The input device 8 is for inputting an operation command to each of the above-mentioned components, and is composed of, for example, a mouse, a trackball, a joystick, or the like.
【００１１】ここで、本発明においては、上記グラフ作
成回路７に画素度数面積計算回路９及びしきい値選択回
路１０が接続され、このしきい値選択回路１０に領域抽
出回路１１が接続されている。上記画素度数面積計算回
路９は、グラフ作成回路７で作成されて出力されたヒス
トグラムのデータを入力し指定された画素値の上限と下
限のしきい値内で画素度数の面積を計算するもので、そ
の計算結果の面積データは前記画像表示回路２へ送出さ
れるようになっている。また、しきい値選択回路１０
は、上記グラフ作成回路７で作成されたヒストグラム上
に前記入力装置８で入力されたしきい値の上限及び下限
を標示線として表示するもので、その標示線のデータは
前記画像表示回路２へ送出されるようになっている。さ
らに、領域抽出回路１１は、前記決定された関心領域内
について上記しきい値選択回路１０で選択された範囲の
画素値のみの領域を抽出して標示するもので、その抽出
領域のデータは画像表示回路２へ送出されるようになっ
ている。In the present invention, the pixel frequency area calculating circuit 9 and the threshold value selecting circuit 10 are connected to the graph creating circuit 7, and the area extracting circuit 11 is connected to the threshold value selecting circuit 10. There is. The pixel frequency area calculation circuit 9 inputs the data of the histogram created and output by the graph creation circuit 7 and calculates the area of the pixel frequency within the upper and lower threshold values of the specified pixel value. The area data as the calculation result is sent to the image display circuit 2. In addition, the threshold selection circuit 10
Displays the upper and lower limits of the threshold value input by the input device 8 as marking lines on the histogram created by the graph creating circuit 7, and the data of the marking lines are sent to the image display circuit 2. It will be sent out. Further, the area extracting circuit 11 extracts and marks an area having only pixel values in the range selected by the threshold selecting circuit 10 within the determined area of interest, and the data of the extracted area is an image. The signal is sent to the display circuit 2.
【００１２】次に、このように構成された医用画像処理
装置におけるしきい値の決定及び領域抽出の動作につい
て、図２に示すフローチャートを参照して説明する。ま
ず、図１に示す記憶装置１から処理対象の画像を読み出
し、一旦画像メモリ１５に書き込む。そして、上記画像
メモリ１５から画像データを読み出し、画像表示回路２
を介して表示装置３に対象画像を表示する（図２のステ
ップＡ）。このときの対象画像の表示例は、図３に示す
ように、表示装置３の画面の一部に対象画像Ｉ0が表示
される。この対象画像Ｉ0は、表示すべき画像の全体を
表示してもよいし、一部分を表示してもよいし、或いは
拡大又は縮小して表示してもよい。Next, the operation of determining the threshold value and extracting the area in the medical image processing apparatus configured as described above will be described with reference to the flowchart shown in FIG. First, the image to be processed is read from the storage device 1 shown in FIG. 1 and is once written in the image memory 15. Then, the image data is read from the image memory 15, and the image display circuit 2
The target image is displayed on the display device 3 via (step A in FIG. 2). In the display example of the target image at this time, as shown in FIG. 3, the target image I 0 is displayed on a part of the screen of the display device 3. The target image I 0 may be displayed as a whole, a part of the image to be displayed, or may be displayed in an enlarged or reduced size.
【００１３】次に、操作者は、上記表示された対象画像
Ｉ0を見ながら図１に示す入力装置８を操作して、関心
領域決定回路４の動作により上記対象画像Ｉ0上で関心
領域を決定する（ステップＢ）。この関心領域の表示例
は、図３に示すように、対象画像Ｉ0上で診断部位の断
層像１６の特定領域を例えば四角形で囲んで関心領域Ｒ
が表示される。そして、この関心領域Ｒの設定は、表示
装置３の画面上にグラフィックの点列で上記特定領域を
囲むことにより行うが、その関心領域Ｒは、断層像１６
の一部を円や四角形などでおおまかに指定するだけでも
よいし、断層像１６の全体でもよいし、或いは断層像１
６上の複数の領域を指定してもよい。Next, the operator operates the input device 8 shown in FIG. 1 while looking at the displayed target image I 0 , and the operation of the ROI determining circuit 4 causes the ROI on the target image I 0. Is determined (step B). In the display example of the region of interest, as shown in FIG. 3, the region of interest R is formed by enclosing a specific region of the tomographic image 16 of the diagnosis region on the target image I 0 with, for example, a rectangle.
Is displayed. Then, the region of interest R is set by enclosing the specific region on the screen of the display device 3 with a graphic sequence of points.
A part of the tomographic image 16 may be roughly specified by a circle or a square, or the entire tomographic image 16 may be specified, or the tomographic image 1
A plurality of areas on 6 may be designated.
【００１４】次に、上記のように決定した関心領域Ｒで
よいかどうか判定する（ステップＣ）。もし、上記の関
心領域Ｒでは不適切の場合は、“NO”側へ進んでステッ
プＢに戻り、再度関心領域Ｒを決定し直す。また、上記
の関心領域Ｒでよい場合は、これで関心領域を確定し、
“YES”側へ進む。そして、上記確定された関心領域Ｒ
内の画素値に対して重み関数（フィルタ）をかけるかど
うか判定する（ステップＤ）。いま、フィルタをかけな
いとすると、ステップＤは“NO”側へ進んでステップＧ
へ跳ぶ。Next, it is judged whether or not the region of interest R determined as described above is acceptable (step C). If the region of interest R is inappropriate, the process proceeds to "NO", returns to step B, and the region of interest R is determined again. If the region of interest R described above is acceptable, the region of interest is fixed with this,
Proceed to the "YES" side. Then, the determined region of interest R
It is determined whether or not to apply a weighting function (filter) to the pixel values in (step D). Now, assuming no filtering is done, step D goes to the "NO" side and step G
Jump to.
【００１５】すると、ステップＧでは、グラフ作成回路
７により、図３に示す関心領域Ｒ内での画素値ごとの度
数を計算してヒストグラムが作成され、画像表示回路２
を介して表示装置３に表示する。また、これと同時に、
画素度数面積計算回路９により、上記グラフ作成回路７
からのヒストグラムのデータを入力し指定された画素値
の上限と下限のしきい値内で画素度数の面積が計算さ
れ、画像表示回路２を介して表示装置３に表示する。こ
のときのヒストグラムムの表示例は、図３に示すよう
に、表示装置３の画面の一部に横軸が画素値で縦軸が画
素度数のヒストグラム１７が表示される。このヒストグ
ラム１７では、関心領域Ｒ内の例えば臓器ごとの画素値
の偏りがそれぞれピークとなって表示される。また、画
素度数の面積の表示例は、例えば上記ヒストグラム１７
の表示位置の上方に面積Ａとして数値表示される。その
後、上記のように表示されたヒストグラム１７及び面積
Ａでよいかどうか判定する（ステップＨ）。これでよい
場合は、“YES”側へ進みステップＪに入る。Then, in step G, the graph creating circuit 7 calculates the frequency for each pixel value in the region of interest R shown in FIG. 3 to create a histogram, and the image display circuit 2
Is displayed on the display device 3 via. At the same time,
The pixel frequency area calculation circuit 9 causes the graph creation circuit 7
The area of the pixel frequency is calculated within the specified upper and lower threshold values of the pixel value by inputting the data of the histogram from, and displayed on the display device 3 via the image display circuit 2. In the display example of the histogram at this time, as shown in FIG. 3, a histogram 17 in which the horizontal axis is the pixel value and the vertical axis is the pixel frequency is displayed on a part of the screen of the display device 3. In the histogram 17, for example, the deviation of the pixel value for each organ in the region of interest R is displayed as a peak. An example of displaying the area of the pixel frequency is, for example, the histogram 17 described above.
The area A is numerically displayed above the display position of. Then, it is determined whether the histogram 17 and the area A displayed as described above are sufficient (step H). If this is acceptable, proceed to the "YES" side and enter step J.
【００１６】操作者は、入力装置８を操作して、表示装
置３に表示されたヒストグラム１７上でしきい値の上限
及び下限を入力し、しきい値を決定する（ステップ
Ｊ）。このとき、上記入力装置８からのしきい値の上限
及び下限の入力値は、しきい値選択回路１０を介して画
像表示回路２へ送られ、表示装置３に表示されたヒスト
グラム１７上に２本の垂直なバー（標示線）１８ａ，１
８ｂとしてリアルタイムに表示される。図３において
は、一方の破線バー１８ａでしきい値の下限が設定さ
れ、他方の破線バー１８ｂでしきい値の上限が設定され
た状態を示している。The operator operates the input device 8 to input the upper and lower limits of the threshold value on the histogram 17 displayed on the display device 3 to determine the threshold value (step J). At this time, the input values of the upper limit and the lower limit of the threshold value from the input device 8 are sent to the image display circuit 2 via the threshold value selection circuit 10 and are displayed on the histogram 17 displayed on the display device 3 as 2 values. Vertical bars (marking lines) 18a, 1 of the book
8b is displayed in real time. In FIG. 3, the lower limit of the threshold value is set by one broken line bar 18a, and the upper limit of the threshold value is set by the other broken line bar 18b.
【００１７】次に、図１に示す領域抽出回路１１の動作
により、上記のように決定されて表示されたしきい値の
上限及び下限の範囲の画素値のみの領域を抽出して表示
する（ステップＫ）。このときの抽出領域の表示例は、
図４（ａ）に示すように、例えば特定の臓器１９ａ，１
９ｂのようになる。前記しきい値選択回路１０は、上記
臓器１９ａ，１９ｂの部分を例えば高輝度で対象画像Ｉ
0上に表示して、抽出された部分が明瞭にわかるように
表示する。又は、表示装置３の画面上で上記対象画像Ｉ
0が表示された位置とは別の場所に、図４（ｂ）に示す
ように、画像Ｉ0′内に抽出された臓器１９ａ，１９ｂ
の部分のみを画像あるいはグラフィックで表示してもよ
いし、さらに、図４（ｃ）に示すように、関心領域Ｒの
部分を切り出しその中に抽出された臓器１９ａ，１９ｂ
の部分のみを表示してもよい。なお、上述の図４
（ａ），（ｂ）においては、関心領域Ｒの外部にも上記
しきい値の上限及び下限の範囲の画素値を有する部分が
存在する場合は、その部分についても抽出して表示して
もよい。Next, the area extracting circuit 11 shown in FIG. 1 operates to extract and display an area having only pixel values in the upper and lower threshold ranges determined and displayed as described above ( Step K). The display example of the extraction area at this time is
As shown in FIG. 4A, for example, specific organs 19a, 1
It looks like 9b. The threshold selection circuit 10 uses the target image I with high brightness, for example, for the parts of the organs 19a and 19b.
It is displayed above 0 so that the extracted part can be clearly seen. Alternatively, the target image I is displayed on the screen of the display device 3.
As shown in FIG. 4B, the organs 19a and 19b extracted in the image I 0 ′ are placed at a position different from the position where 0 is displayed.
May be displayed as an image or a graphic, and as shown in FIG. 4 (c), the region of interest R is cut out and the organs 19a and 19b extracted therein are extracted.
You may display only the part of. Note that FIG.
In (a) and (b), if there is a portion having pixel values in the upper and lower limits of the threshold value outside the region of interest R, that portion may also be extracted and displayed. Good.
【００１８】次に、操作者は、図４に示すように抽出し
て表示された領域を見てこれでよいかどうか判定する
（ステップＬ）。これでよい場合は、“YES”側へ進ん
で処理を終了する。もし、この抽出領域では不適切であ
れば、“NO”側へ進みステップＨは“YES”側へ進んで
ステップＪに戻り、しきい値の決定をやり直す。また、
上記関心領域Ｒの外部にも抽出された部分が存在し、そ
の部分も抽出領域に含めたい場合は、ステップＬで“N
O”側へ進み、その後ステップＨ及びステップＣを介し
てステップＢまで戻り、再度関心領域Ｒを決定し直す。
その後、ステップＣを介してステップＤへ進み、この状
態でフィルタをかけるかどうか判定する。Next, the operator looks at the area extracted and displayed as shown in FIG. 4 and determines whether or not this is enough (step L). If this is all right, proceed to the “YES” side and end the processing. If this extraction area is inappropriate, the process proceeds to "NO" side, the step H proceeds to "YES" side, returns to step J, and the threshold value is decided again. Also,
If there is an extracted part outside the region of interest R, and if that part is also to be included in the extracted region, in step L, "N
The process proceeds to the O ″ side, then returns to step B via step H and step C, and the region of interest R is determined again.
After that, the process proceeds to Step D through Step C, and it is determined whether or not the filter is applied in this state.
【００１９】ここで、上記決定し直した関心領域Ｒ内に
ついてその画像のヒストグラムから容易に臓器等の境界
を決定することが困難と思われる場合には、画素値に対
して重み関数（フィルタ）をかけることとし、ステップ
Ｄは“YES”側へ進む。そして、上記関心領域決定回路
４で決定された関心領域Ｒ内の画素値に対してかける重
み関数を、フィルタ発生回路５で生成する（ステップ
Ｅ）。この重み関数は、画素値の変化の小さいところは
大きくなるように、また不要な範囲は狭帯域に圧縮する
ように与えればよい。その生成の仕方は、ヒストグラム
均等化法を行って自動的に発生させたり、予め何種類か
の重み関数を用意しておいたり、或いはトラックボール
などの入力装置８を操作して操作者が手動で入力しても
よい。いずれの場合にも、操作者はヒストグラムを見な
がら図５に示すように、重み関数２０を表示装置３に表
示されたヒストグラム１７に対応付けて与えることがで
きる。If it is difficult to easily determine the boundary of an organ or the like from the histogram of the image in the re-determined region of interest R, a weighting function (filter) is applied to the pixel value. And step D proceeds to the “YES” side. Then, the filter generating circuit 5 generates a weighting function to be applied to the pixel values in the region of interest R determined by the region of interest determining circuit 4 (step E). This weighting function may be given so that the portion where the change of the pixel value is small becomes large and the unnecessary range is compressed to a narrow band. The method of generating the histogram may be automatically generated by performing a histogram equalization method, some kinds of weighting functions may be prepared in advance, or an operator may manually operate the input device 8 such as a trackball. You may enter with. In any case, the operator can give the weighting function 20 in association with the histogram 17 displayed on the display device 3 while looking at the histogram, as shown in FIG.
【００２０】次に、このように生成された重み関数（フ
ィルタ）２０を画像表示回路２を介して表示装置３の画
面に表示する（ステップＦ）。この重み関数２０の表示
例は、図５（ａ）に示すように、同図（ｂ）に示すヒス
トグラム１７に対応付けてグラフで表示されるため、画
素値と重み関数２０との対応が明確になり、例えば臓器
１９ａ，１９ｂの画素境界を選択し易いヒストグラム１
７を得ることができるようになる。そして、画素値ごと
に重み関数をかけ合わせ、画素値の線形あるいは非線形
の変換を行えばよい。その後、前述のようにステップＧ
→Ｈ→Ｊ→Ｋ→Ｌのように進んで、しきい値を決定する
と共に、特定領域を抽出して表示する。これにより、表
示装置３に表示されたヒストグラム１７上で簡便に、し
かも抽出すべき領域の画素度数の面積の増減割合を意識
しながら、最適なしきい値を設定することができる。Next, the weighting function (filter) 20 thus generated is displayed on the screen of the display device 3 via the image display circuit 2 (step F). As shown in FIG. 5A, the display example of the weighting function 20 is displayed in a graph in association with the histogram 17 shown in FIG. 5B, so that the correspondence between the pixel value and the weighting function 20 is clear. Histogram 1 that makes it easy to select pixel boundaries of organs 19a and 19b, for example.
You will be able to get 7. Then, the weighting function may be multiplied for each pixel value to perform linear or non-linear conversion of the pixel value. Then step G as described above.
→ H → J → K → L, the threshold value is determined, and the specific area is extracted and displayed. As a result, the optimum threshold value can be set easily on the histogram 17 displayed on the display device 3 while being aware of the increase / decrease rate of the pixel frequency of the region to be extracted.
【００２１】図６は本発明における画像処理の第二の実
施例を説明するための表示画面の一例を示した図であ
る。この第二の実施例は、対象画像Ｉ0の中に複数の関
心領域Ｒ1，Ｒ2を設定し、それぞれの関心領域Ｒ1，Ｒ2
に対しヒストグラム１７1，１７2を表示して、しきい値
の設定を行うものである。すなわち、表示装置３の画面
に表示された対象画像Ｉ0上の断層像１６には、それぞ
れ特定領域を囲んだ二つの関心領域Ｒ1，Ｒ2が設定され
ており、第一の関心領域Ｒ1に対応して第一のヒストグ
ラム１７1が作成表示され、第二の関心領域Ｒ2に対応し
て第二のヒストグラム１７2が作成表示されている。な
お、第二の関心領域Ｒ2内に存在する符号２１は、例え
ば腫瘍部を示しており、その周囲より高い画素値を有し
高輝度に表示されている。FIG. 6 is a diagram showing an example of a display screen for explaining the second embodiment of the image processing in the present invention. The second embodiment sets a plurality of regions of interest R 1, R 2 in the target image I 0, each of region of interest R 1, R 2
On the other hand, the histograms 17 1 and 17 2 are displayed to set the threshold value. That is, in the tomographic image 16 on the target image I 0 displayed on the screen of the display device 3, two regions of interest R 1 and R 2 surrounding the specific region are set, and the first region of interest R 1 the first histogram 17 1 appears created corresponding to 1, the second histogram 17 2 are displayed created corresponding to the second region of interest R 2. Note that the reference numeral 21 existing in the second region of interest R 2 indicates, for example, a tumorous part, which has a pixel value higher than the surroundings and is displayed with high brightness.
【００２２】この状態で、まず、第一のヒストグラム１
７1上でしきい値の下限２２ａ及び上限２２ｂを設定す
る。このようにしきい値が設定されると、次に画素度数
面積計算回路９により前記第一の関心領域Ｒ1における
しきい値の下限２２ａと上限２２ｂの範囲内の画素度数
の面積が計算され、表示装置３に面積Ａ1として数値表
示される。これと同時に、第二のヒストグラム１７2上
においても、第一のヒストグラム１７1上で設定したし
きい値の範囲が高輝度などで区別して表示される。ま
た、上記と同様にして、第二のヒストグラム１７2上で
しきい値の下限２３ａ及び上限２３ｂを設定する。する
と、画素度数面積計算回路９により前記第二の関心領域
Ｒ2におけけるしきい値の下限２３ａと上限２３ｂの範
囲内の画素度数の面積が計算され、表示装置３に面積Ａ
2として数値表示される。このとき、第二の関心領域Ｒ2
内には高輝度に表示される腫瘍部２１が存在し、この腫
瘍部２１は図７の第二のヒストグラム１７2上でハッチ
ンングを施した領域に対応している。In this state, first, the first histogram 1
The lower limit 22a and the upper limit 22b of the threshold value are set on 7 1 . When the threshold value is set in this way, the pixel frequency area calculation circuit 9 next calculates the area of the pixel frequency within the range of the lower limit 22a and the upper limit 22b of the threshold in the first region of interest R 1 , The display device 3 numerically displays the area A 1 . At the same time, the range of the threshold value set on the first histogram 17 1 is also displayed on the second histogram 17 2 so as to be distinguished by high brightness. Further, similarly to the above, the lower limit 23a and the upper limit 23b of the threshold value are set on the second histogram 17 2 . Then, the pixel frequency area calculation circuit 9 calculates the area of the pixel frequency within the range of the lower limit 23a and the upper limit 23b of the threshold in the second region of interest R 2 , and the area A is displayed on the display device 3.
Numerically displayed as 2 . At this time, the second region of interest R 2
There is a tumor part 21 displayed in high brightness inside, and this tumor part 21 corresponds to the hatched region on the second histogram 17 2 in FIG. 7.
【００２３】ここで、図６に示す対象画像Ｉ0におい
て、第一及び第二の関心領域Ｒ1，Ｒ2で囲んだ臓器１９
ａ，１９ｂが略対称に存在するものであれば、抽出され
る面積は略等しくなるはずである。しかし、第二の関心
領域Ｒ2内には腫瘍部２１が存在するために、抽出され
る面積は異なることとなる。このことは、面積Ａ1と面
積Ａ2との数値表示が異なることから理解できる。この
とき、第二の関心領域Ｒ2内で腫瘍部２１を選択しない
場合は、第二のヒストグラム１７2においてしきい値を
多少変更して、しきい値の下限２３ａ及び上限２３ｂを
決定する。しかしながら、第二の関心領域Ｒ2内に腫瘍
部２１を含める場合は、図７に示すように、第二のヒス
トグラム１７2において一ヶ所のしきい値２３ａ，２３
ｂ以外に、もう一ヶ所のしきい値の下限２４ａ及び上限
２４ｂも選択する必要がある。このようにすると、図７
において、第一及び第二の関心領域Ｒ1，Ｒ2で抽出され
る面積を略等しくすることができる。このことは、面積
Ａ1と面積Ａ2との数値表示が略等しくなることからも理
解できる。Here, in the target image I 0 shown in FIG. 6, an organ 19 surrounded by the first and second regions of interest R 1 and R 2.
If a and 19b exist substantially symmetrically, the extracted areas should be substantially equal. However, since the tumor part 21 exists in the second region of interest R 2 , the extracted areas will be different. This can be understood from the fact that the numerical representations of the area A 1 and the area A 2 are different. At this time, when the tumor part 21 is not selected in the second region of interest R 2 , the threshold value is slightly changed in the second histogram 17 2 to determine the lower limit 23a and the upper limit 23b of the threshold value. However, when the tumorous part 21 is included in the second region of interest R 2 , as shown in FIG. 7, one threshold value 23 a, 23 in the second histogram 17 2 is included.
In addition to b, it is also necessary to select another threshold lower limit 24a and upper limit 24b. In this way, FIG.
In, the areas extracted by the first and second regions of interest R 1 and R 2 can be made substantially equal. This can be understood from the fact that the numerical representations of the area A 1 and the area A 2 are substantially equal.
【００２４】図８は本発明における画像処理の第三の実
施例を説明するための表示画面の一例を示した図であ
る。この第三の実施例は、２枚の対象画像Ｉ1，Ｉ2から
より精度の良いしきい値を決定しようとするものであ
る。すなわち、計測条件等によりコントラストの異なる
２枚の画像から、特定領域をより正確かつ容易に抽出す
ることを可能とするものである。計測条件の異なるもの
としては、ＭＲＩ装置におけるエコー時間や計測シーケ
ンスの異なる画像がある。また、計測装置の異なるもの
としては、Ｘ線ＣＴ装置による計測画像やＭＲＩ装置に
よる計測画像などがある。図８においては、第一の対象
画像Ｉ1は、例えばＭＲＩ装置で計測された第一エコー
画像であり、第二の対象画像Ｉ2は、同じくＭＲＩ装置
で計測された第二エコー画像である。FIG. 8 is a diagram showing an example of a display screen for explaining the third embodiment of the image processing in the present invention. The third embodiment is intended to determine a more accurate threshold value from the two target images I 1 and I 2 . That is, the specific area can be more accurately and easily extracted from the two images having different contrasts depending on the measurement conditions and the like. Images with different echo times and measurement sequences in the MRI apparatus are examples of different measurement conditions. Further, as a different measuring device, there are a measurement image by an X-ray CT device, a measurement image by an MRI device, and the like. In FIG. 8, the first target image I 1 is, for example, the first echo image measured by the MRI apparatus, and the second target image I 2 is the second echo image also measured by the MRI apparatus. .
【００２５】このように、表示装置３の画面に２枚の対
象画像Ｉ1，Ｉ2が表示された状態で、それぞれの断層像
１６上にて各々関心領域Ｒ1，Ｒ2を同一の位置に決定す
る。そして、各々の関心領域Ｒ1，Ｒ2内の画像に対して
フィルタをかける場合は、同一のフィルタをかけてもよ
いし、別々に異なるフィルタをかけてもよい。その後、
しきい値の決定は、どちらの対象画像Ｉ1，Ｉ2上で行っ
てもよい。図８においては、第一の対象画像Ｉ1上の関
心領域Ｒ1についてのヒストグラム２５1上でしきい値の
下限２６ａと上限２６ｂとが設定され、この範囲でしき
い値が決定されている。このとき、前述と同様に、上記
しきい値の範囲に対応した抽出領域も表示される。ま
た、第二の対象画像Ｉ2上においても、その関心領域Ｒ2
についてのヒストグラム２５2上で対応する同じ抽出領
域が表示され、第一及び第二の関心領域Ｒ1，Ｒ2で抽出
される面積Ａ1，Ａ2は、等しい。さらに、第二のヒスト
グラム２５2上において、上記しきい値の範囲と同じ画
素値を有しハッチングの施された範囲２７が高輝度で表
示され、当該第二のヒストグラム２５2上でも抽出領域
を確認することができる。このようにして、１枚の対象
画像Ｉ1からでは例えば臓器の境界が不明確な場合で
も、異種の対象画像Ｉ2から高精度にしきい値を決定し
て、特定領域を抽出することができる。Thus, in the state where the two target images I 1 and I 2 are displayed on the screen of the display device 3, the regions of interest R 1 and R 2 are located at the same position on the respective tomographic images 16. To decide. When the images in the regions of interest R 1 and R 2 are filtered, the same filter may be applied or different filters may be applied separately. afterwards,
The threshold value may be determined on either of the target images I 1 and I 2 . In Figure 8, a lower 26a and upper 26b threshold on the histogram 25 1 of region of interest R 1 on the first target image I 1 is set, the threshold in this range is determined . At this time, similarly to the above, the extraction area corresponding to the range of the threshold value is also displayed. Further, also on the second target image I 2 , the region of interest R 2
The same extraction area corresponding is displayed on the histogram 25 2 for the area A 1, A 2, which is extracted by the first and second region of interest R 1, R 2 are equal. Further, in the second histogram 25 on 2, range 27 subjected to the hatching have the same pixel value range of the threshold value is displayed with high luminance, the extraction region on the second histogram 25 2 You can check. In this way, even if the boundary of an organ is unclear from one target image I 1 , for example, it is possible to accurately determine the threshold value from the different target images I 2 and extract the specific region. .
【００２６】図９は本発明における画像処理の第四の実
施例を説明するための表示画面の一例を示した図であ
る。この第四の実施例は、例えば診断部位の多数枚のス
ライス画像から成る対象画像からそれぞれ特定領域を抽
出することにより、立体画像の抽出を可能とするもので
ある。すなわち、一連のスライス画像から例えば代表的
な４枚を選択して対象画像Ｉ1〜Ｉ4とし、それぞれにつ
いて関心領域Ｒ1〜Ｒ4を決定し、さらに各関心領域Ｒ1
〜Ｒ4内の画素値についてヒストグラムを作成する。す
ると、図９に示すように、４枚の対象画像Ｉ1〜Ｉ4の側
方に、それぞれの関心領域Ｒ1〜Ｒ4内について作成され
たヒストグラム２７1〜２７4が表示される。FIG. 9 is a diagram showing an example of a display screen for explaining the fourth embodiment of the image processing in the present invention. In the fourth embodiment, for example, a stereoscopic image can be extracted by extracting a specific region from a target image composed of a large number of slice images of a diagnosis site. That is, for example, four representative images are selected from a series of slice images as target images I 1 to I 4 , regions of interest R 1 to R 4 are determined for each, and each region of interest R 1 is further determined.
Creating a histogram for pixel values ~R 4. Then, as shown in FIG. 9, on the side of the four target image I 1 ~I 4, the histogram 27 1-27 4 created for each region of interest R 1 to R 4 are displayed.
【００２７】次に、上記いずれかのヒストグラム２７1
〜２７4上でしきい値を決定する。図９では、第一のヒ
ストグラム２７1上でしきい値の下限２８ａと上限２８
ｂの範囲を設定した状態を示している。これにより、各
関心領域Ｒ1〜Ｒ4内において、図１０でハッチングを施
して示すように、上記しきい値の下限２８ａ及び上限２
８ｂの範囲の画素値のみの領域が抽出され、例えば高輝
度に表示される。同時に、各抽出領域の面積が計算さ
れ、それぞれ面積Ａ1〜Ａ4として数値表示される。この
場合、各ヒストグラム２７1〜２７4上でも各対象画像Ｉ
1〜Ｉ4における抽出領域を確認することができる。な
お、上記しきい値の設定は、第一のヒストグラム２７1
上に限られず、他のヒストグラム２７2〜２７4のいずれ
を用いて行ってもよい。また、対象画像Ｉ1〜Ｉ4ごとに
しきい値を変化させてもよい。このように、多数枚の対
象画像Ｉ1〜Ｉ4より画素度数の面積を意識してそれぞれ
の特定領域を抽出することにより、例えば診断部位のス
ライス方向のスライス画像の面積変化から抽出しようと
する立体全体をおおまかにとらえた上で、立体画像の抽
出が可能となる。Next, one of the above histograms 27 1
To determine the threshold on to 27 4. In FIG. 9, the lower limit 28a and the upper limit 28 of the threshold value are displayed on the first histogram 27 1.
The state where the range of b is set is shown. As a result, in each of the regions of interest R 1 to R 4 , as shown by hatching in FIG. 10, the lower limit 28 a and the upper limit 2 of the threshold value are set.
A region having only pixel values in the range of 8b is extracted and displayed with high brightness, for example. At the same time, the area of each extraction region is calculated and numerically displayed as the areas A 1 to A 4 . In this case, each target image I is also displayed on each of the histograms 27 1 to 27 4.
The extraction area in 1 to I 4 can be confirmed. In addition, the above threshold value is set by the first histogram 27 1
The histogram is not limited to the above, and any of the other histograms 27 2 to 27 4 may be used. Further, the threshold value may be changed for each of the target images I 1 to I 4 . In this way, by extracting the respective specific regions from the large number of target images I 1 to I 4 by considering the area of the pixel frequency, it is attempted to extract from the area change of the slice image in the slice direction of the diagnosis region, for example. It is possible to extract a stereoscopic image after roughly capturing the entire stereoscopic image.
【００２８】図１１はヒストグラム表示の他の例を示す
説明図である。以上の説明では、ヒストグラムの表示
は、図１１（ａ）に示すように、横軸に画素値をとると
共に縦軸に画素度数をとって棒グラフのように表示する
ものとしたが、これに限らず、図１１（ｂ）に示すよう
に、横軸に画素度数をとり、画素値の範囲を示したもの
とし画素値の範囲ごとに例えば異なる輝度で表示するよ
うにしてもよい。FIG. 11 is an explanatory diagram showing another example of the histogram display. In the above description, the histogram is displayed as a bar graph with the pixel value on the horizontal axis and the pixel frequency on the vertical axis as shown in FIG. 11A, but the present invention is not limited to this. Alternatively, as shown in FIG. 11B, the pixel frequency may be taken on the horizontal axis, and the range of pixel values may be shown, and for example, different brightness may be displayed for each range of pixel values.
【００２９】図１２はヒストグラム表示の更に他の例を
示す説明図である。この表示例は、縦軸の画素度数を積
分値で示したものである。図１２（ａ）の折れ線グラフ
表示においては、しきい値の下限２９ａ及び上限２９ｂ
の範囲が領域抽出をするために設定した部分であり、矢
印で示した範囲３０が画素度数すなわち面積に対応す
る。また、図１２（ｂ）の棒グラフ表示においては、し
きい値の下限３１ａ及び上限３１ｂの範囲が領域抽出を
するために設定した部分であり、ハッチングを施して示
した範囲３２が画素度数すなわち面積に対応する。FIG. 12 is an explanatory view showing still another example of the histogram display. In this display example, the pixel frequency on the vertical axis is shown by an integrated value. In the line graph display of FIG. 12A, the lower limit 29a and the upper limit 29b of the threshold value are displayed.
The range is a part set for area extraction, and the range 30 indicated by an arrow corresponds to the pixel frequency, that is, the area. Further, in the bar graph display of FIG. 12B, the range of the lower limit 31a and the upper limit 31b of the threshold value is a portion set for area extraction, and the hatched range 32 is the pixel frequency, that is, the area. Corresponding to.
【００３０】[0030]
【発明の効果】本発明は以上のように構成されたので、
画素度数面積計算回路でグラフ作成回路から出力された
ヒストグラムのデータを入力して指定された画素値の上
限と下限のしきい値内で画素度数の面積を計算し、しき
い値選択回路により入力装置で入力されたしきい値の上
限及び下限をヒストグラム上に標示線として表示し、領
域抽出回路で上記しきい値選択回路により選択された範
囲の画素値のみの領域を抽出して表示することができ
る。これにより、表示された医用画像から臓器抽出や任
意の領域抽出を行う際にその抽出のための画素値のしき
い値を容易かつ最適に選択して、特定領域を抽出するこ
とができる。Since the present invention is constructed as described above,
The pixel frequency area calculation circuit inputs the histogram data output from the graph creation circuit to calculate the pixel frequency area within the upper and lower thresholds of the specified pixel value, and the threshold selection circuit inputs it. Display the upper and lower limits of the threshold value input by the device as marking lines on the histogram, and extract and display the area of only the pixel values in the range selected by the threshold value selection circuit by the area extraction circuit. You can Accordingly, when performing organ extraction or arbitrary area extraction from the displayed medical image, the threshold value of the pixel value for the extraction can be easily and optimally selected to extract the specific area.
【図１】本発明による医用画像処理装置の実施例を示す
ブロック図である。FIG. 1 is a block diagram showing an embodiment of a medical image processing apparatus according to the present invention.
【図２】上記医用画像処理装置におけるしきい値の決定
及び領域抽出の動作を説明するためのフローチャートで
ある。FIG. 2 is a flowchart for explaining the operations of threshold determination and area extraction in the medical image processing apparatus.
【図３】対象画像及びヒストグラムの表示例を示す説明
図である。FIG. 3 is an explanatory diagram showing a display example of a target image and a histogram.
【図４】対象画像からの特定領域の抽出の表示例を示す
説明図である。FIG. 4 is an explanatory diagram showing a display example of extraction of a specific area from a target image.
【図５】フィルタをかける際の重み関数の表示例を示す
説明図である。FIG. 5 is an explanatory diagram showing a display example of a weighting function when a filter is applied.
【図６】本発明における画像処理の第二の実施例を説明
するための表示画面の一例を示した図である。FIG. 6 is a diagram showing an example of a display screen for explaining a second embodiment of image processing according to the present invention.
【図７】上記第二の実施例においてしきい値を二ケ所に
設定した例を示す説明図である。FIG. 7 is an explanatory diagram showing an example in which threshold values are set in two places in the second embodiment.
【図８】本発明における画像処理の第三の実施例を説明
するための表示画面の一例を示した図である。FIG. 8 is a diagram showing an example of a display screen for explaining a third embodiment of image processing in the present invention.
【図９】本発明における画像処理の第四の実施例を説明
するための表示画面の一例を示した図である。FIG. 9 is a diagram showing an example of a display screen for explaining a fourth embodiment of image processing according to the present invention.
【図１０】図９において特定領域が抽出された状態を示
す説明図である。10 is an explanatory diagram showing a state in which a specific area is extracted in FIG.
【図１１】ヒストグラム表示の他の例を示す説明図であ
る。FIG. 11 is an explanatory diagram showing another example of histogram display.
【図１２】ヒストグラム表示の更に他の例を示す説明図
である。FIG. 12 is an explanatory diagram showing still another example of histogram display.
１…記憶装置 ２…画像表示回路 ３…表示装置 ４…関心領域決定回路 ５…フィルタ発生回路 ６…フィルタ回路 ７…グラフ作成回路 ８…入力装置 ９…画素度数面積計算回路 １０…しきい値選択回路 １１…領域抽出回路 DESCRIPTION OF SYMBOLS 1 ... Storage device 2 ... Image display circuit 3 ... Display device 4 ... Region of interest determination circuit 5 ... Filter generation circuit 6 ... Filter circuit 7 ... Graph creation circuit 8 ... Input device 9 ... Pixel frequency area calculation circuit 10 ... Threshold selection Circuit 11 ... Region extraction circuit
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｇ０６Ｔ 5/20 9/20 9191−5Ｌ Ｇ０６Ｆ 15/68 ４００ Ａ 7459−5Ｌ 15/70 ３３５ Ｚ ─────────────────────────────────────────────────── ─── Continuation of the front page (51) Int.Cl. 6 Identification code Internal reference number FI Technical display location G06T 5/20 9/20 9191-5L G06F 15/68 400 A 7459-5L 15/70 335 Z
Claims (1)
記憶装置と、この記憶装置から読み出した画像データを
表示するための処理を行う画像表示回路と、この画像表
示回路から出力された画像データを画像として表示する
表示装置と、この表示装置に表示された画像上で関心領
域を決定する関心領域決定回路と、この関心領域決定回
路で決定された関心領域内の画素値に対してかける重み
関数を生成するフィルタ発生回路と、このフィルタ発生
回路で生成された重み関数を上記関心領域に対してかけ
合わせるフィルタ回路と、上記関心領域内での画素値ご
との度数を計算してグラフを作成して上記画像表示回路
へ送出するグラフ作成回路と、上記各構成要素に操作指
令を入力する入力装置とを備えて成る医用画像処理装置
において、上記グラフ作成回路からのグラフのデータを
入力し指定された画素値の上限と下限のしきい値内で画
素度数の面積を計算する画素度数面積計算回路と、上記
入力装置で入力されたしきい値の上限及び下限をグラフ
上に標示線として標示するしきい値選択回路と、このし
きい値選択回路で選択された範囲の画素値のみの領域を
抽出して表示する領域抽出回路とを設けたことを特徴と
する医用画像処理装置。1. A storage device for recording data of a plurality of medical images, an image display circuit for performing processing for displaying image data read from the storage device, and an image output from the image display circuit. A display device that displays data as an image, a region of interest determination circuit that determines a region of interest on the image displayed on the display device, and a pixel value in the region of interest determined by the region of interest determination circuit A filter generating circuit that generates a weighting function, a filter circuit that multiplies the weighting function generated by this filter generating circuit to the region of interest, and a frequency is calculated for each pixel value in the region of interest to calculate a graph. In a medical image processing apparatus comprising a graph creating circuit for creating and sending the graph to the image display circuit, and an input device for inputting an operation command to each of the components, Pixel frequency area calculation circuit that inputs the graph data from the drawing circuit and calculates the pixel frequency area within the specified upper and lower threshold values of the pixel value, and the threshold value input by the above input device A threshold selection circuit for indicating the upper and lower limits of the above as a marking line on the graph, and an area extraction circuit for extracting and displaying an area of only pixel values in the range selected by this threshold selection circuit are provided. A medical image processing apparatus characterized by the above.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5292365A JPH07129751A (en) | 1993-10-29 | 1993-10-29 | Medical picture processor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP5292365A JPH07129751A (en) | 1993-10-29 | 1993-10-29 | Medical picture processor |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH07129751A true JPH07129751A (en) | 1995-05-19 |
Family
ID=17780866
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP5292365A Pending JPH07129751A (en) | 1993-10-29 | 1993-10-29 | Medical picture processor |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH07129751A (en) |
Cited By (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09212646A (en) * | 1996-02-02 | 1997-08-15 | Fuji Photo Film Co Ltd | Image analysis device |
| JP2001056862A (en) * | 1999-05-07 | 2001-02-27 | Monorisu:Kk | Method and device for tracking image area |
| JP2002042110A (en) * | 2000-07-27 | 2002-02-08 | Ateeru:Kk | Image measuring system, image measuring method, and computer-readable recording medium recording image measuring program |
| JP2003010171A (en) * | 2001-06-27 | 2003-01-14 | Hitachi Medical Corp | Method and device for supporting medical image diagnosis |
| JP2004093510A (en) * | 2002-09-03 | 2004-03-25 | Horiba Ltd | Analysis result display method, x-ray analyzer, and computer program |
| JP2004358001A (en) * | 2003-06-05 | 2004-12-24 | Ziosoft Inc | Medical image processor, medical image processing method and program |
| JP2005278786A (en) * | 2004-03-29 | 2005-10-13 | Toshiba Corp | Medical image processor |
| JP2005342516A (en) * | 2004-06-02 | 2005-12-15 | Ge Medical Systems Global Technology Co Llc | Method and apparatus for simultaneous displaying inversion mode ultrasonic image and histogram information |
| JP2006167187A (en) * | 2004-12-16 | 2006-06-29 | Hitachi Medical Corp | Medical image display apparatus |
| JP2006180987A (en) * | 2004-12-27 | 2006-07-13 | Hitachi Medical Corp | Medical image display device |
| JP2006325629A (en) * | 2005-05-23 | 2006-12-07 | Ge Medical Systems Global Technology Co Llc | Three-dimensional interest region setting method, image acquisition apparatus, and program |
| JP2008009627A (en) * | 2006-06-28 | 2008-01-17 | Canon Machinery Inc | Recognized image creation method and positioning device for recognition target |
| JP2008289916A (en) * | 2000-06-30 | 2008-12-04 | Hitachi Medical Corp | Image diagnosis supporting device |
| JP2009219610A (en) * | 2008-03-14 | 2009-10-01 | Yokohama Tlo Co Ltd | Organ area specifying method and organ area specifying instrument |
| JP2009545800A (en) * | 2006-08-02 | 2009-12-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | How to reconstruct a voxel cluster map in an image |
| US7639855B2 (en) | 2003-04-02 | 2009-12-29 | Ziosoft, Inc. | Medical image processing apparatus, and medical image processing method |
| JP2010119850A (en) * | 2008-11-22 | 2010-06-03 | General Electric Co <Ge> | System, apparatus, and process for automated medical image segmentation using statistical model |
| JP2011083636A (en) * | 2003-05-08 | 2011-04-28 | Hitachi Medical Corp | Ultrasonic diagnostic apparatus |
| JP2013042767A (en) * | 2011-08-22 | 2013-03-04 | Asahikawa Medical College | Image processor, image processing method, and program |
| JP2013056040A (en) * | 2011-09-08 | 2013-03-28 | Olympus Corp | Fluorescence observation device |
| JP2013169433A (en) * | 2012-02-23 | 2013-09-02 | Aze Ltd | Medical image generation processing device and program |
| JP2013233470A (en) * | 2013-08-30 | 2013-11-21 | Canon Inc | Processing apparatus, processing method, and program |
| JP2015156894A (en) * | 2014-02-21 | 2015-09-03 | 学校法人日本大学 | Medical image processor, medical object area extraction method thereof and medical object area extraction processing program |
| JP2020151450A (en) * | 2019-03-19 | 2020-09-24 | 緯創資通股▲ふん▼有限公司Ｗｉｓｔｒｏｎ Ｃｏｒｐｏｒａｔｉｏｎ | Image identification method and image identification device |
| WO2020188970A1 (en) * | 2019-03-15 | 2020-09-24 | コニカミノルタ株式会社 | Display control device, radiographic image capturing system, and program |
-
1993
- 1993-10-29 JP JP5292365A patent/JPH07129751A/en active Pending
Cited By (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09212646A (en) * | 1996-02-02 | 1997-08-15 | Fuji Photo Film Co Ltd | Image analysis device |
| JP2001056862A (en) * | 1999-05-07 | 2001-02-27 | Monorisu:Kk | Method and device for tracking image area |
| JP2008289916A (en) * | 2000-06-30 | 2008-12-04 | Hitachi Medical Corp | Image diagnosis supporting device |
| JP2002042110A (en) * | 2000-07-27 | 2002-02-08 | Ateeru:Kk | Image measuring system, image measuring method, and computer-readable recording medium recording image measuring program |
| JP2003010171A (en) * | 2001-06-27 | 2003-01-14 | Hitachi Medical Corp | Method and device for supporting medical image diagnosis |
| JP2004093510A (en) * | 2002-09-03 | 2004-03-25 | Horiba Ltd | Analysis result display method, x-ray analyzer, and computer program |
| US7639855B2 (en) | 2003-04-02 | 2009-12-29 | Ziosoft, Inc. | Medical image processing apparatus, and medical image processing method |
| US9471981B2 (en) | 2003-05-08 | 2016-10-18 | Hitachi Medical Corporation | Reference image display method for ultrasonography and ultrasonic diagnosis apparatus |
| JP2011083636A (en) * | 2003-05-08 | 2011-04-28 | Hitachi Medical Corp | Ultrasonic diagnostic apparatus |
| JP2004358001A (en) * | 2003-06-05 | 2004-12-24 | Ziosoft Inc | Medical image processor, medical image processing method and program |
| JP2005278786A (en) * | 2004-03-29 | 2005-10-13 | Toshiba Corp | Medical image processor |
| JP2005342516A (en) * | 2004-06-02 | 2005-12-15 | Ge Medical Systems Global Technology Co Llc | Method and apparatus for simultaneous displaying inversion mode ultrasonic image and histogram information |
| JP2006167187A (en) * | 2004-12-16 | 2006-06-29 | Hitachi Medical Corp | Medical image display apparatus |
| JP2006180987A (en) * | 2004-12-27 | 2006-07-13 | Hitachi Medical Corp | Medical image display device |
| JP2006325629A (en) * | 2005-05-23 | 2006-12-07 | Ge Medical Systems Global Technology Co Llc | Three-dimensional interest region setting method, image acquisition apparatus, and program |
| JP2008009627A (en) * | 2006-06-28 | 2008-01-17 | Canon Machinery Inc | Recognized image creation method and positioning device for recognition target |
| JP2009545800A (en) * | 2006-08-02 | 2009-12-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | How to reconstruct a voxel cluster map in an image |
| JP2009219610A (en) * | 2008-03-14 | 2009-10-01 | Yokohama Tlo Co Ltd | Organ area specifying method and organ area specifying instrument |
| JP2010119850A (en) * | 2008-11-22 | 2010-06-03 | General Electric Co <Ge> | System, apparatus, and process for automated medical image segmentation using statistical model |
| JP2013042767A (en) * | 2011-08-22 | 2013-03-04 | Asahikawa Medical College | Image processor, image processing method, and program |
| JP2013056040A (en) * | 2011-09-08 | 2013-03-28 | Olympus Corp | Fluorescence observation device |
| US9417188B2 (en) | 2011-09-08 | 2016-08-16 | Olympus Corporation | Fluorescence observation device |
| JP2013169433A (en) * | 2012-02-23 | 2013-09-02 | Aze Ltd | Medical image generation processing device and program |
| JP2013233470A (en) * | 2013-08-30 | 2013-11-21 | Canon Inc | Processing apparatus, processing method, and program |
| JP2015156894A (en) * | 2014-02-21 | 2015-09-03 | 学校法人日本大学 | Medical image processor, medical object area extraction method thereof and medical object area extraction processing program |
| WO2020188970A1 (en) * | 2019-03-15 | 2020-09-24 | コニカミノルタ株式会社 | Display control device, radiographic image capturing system, and program |
| JP2020151450A (en) * | 2019-03-19 | 2020-09-24 | 緯創資通股▲ふん▼有限公司Ｗｉｓｔｒｏｎ Ｃｏｒｐｏｒａｔｉｏｎ | Image identification method and image identification device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH07129751A (en) | Medical picture processor | |
| US8270688B2 (en) | Method for intelligent qualitative and quantitative analysis assisting digital or digitized radiography softcopy reading | |
| US7324673B1 (en) | Image processing apparatus, image processing method, and image processing program storage medium | |
| US6058218A (en) | Enhanced visualization of weak image sources in the vicinity of dominant sources | |
| JP5252804B2 (en) | Method for forming a four-dimensional medical image | |
| US8107700B2 (en) | System and method for efficient workflow in reading medical image data | |
| DE69629445T2 (en) | Automatic tone scale adjustment using image activity measurements | |
| JP7194143B2 (en) | Systems and methods to facilitate review of liver tumor cases | |
| JP5367704B2 (en) | Image processing apparatus, image processing method, and image processing program | |
| US20050251014A1 (en) | Lesion marking and characterization quality assurance method and system | |
| US8442282B2 (en) | Computer-aided imaging diagnostic processing apparatus and computer-aided imaging diagnostic processing method | |
| JP2008515466A (en) | Method and system for identifying an image representation of a class of objects | |
| JP2013126575A (en) | Image processing device and method, and program | |
| US5509084A (en) | Method and apparatus for postprocessing medical images | |
| JP6775331B2 (en) | Processing to create a transmitted image without artificial noise | |
| US20060198552A1 (en) | Image processing method for a digital medical examination image | |
| JP6102236B2 (en) | Radiation image processing device | |
| CN107787506A (en) | Select the transmission function for showing medical image | |
| US20090190810A1 (en) | Semiautomatic myocardial region of interest definition | |
| JPH01147619A (en) | Graph display device and improvement in visibility thereof | |
| JP2001087228A (en) | Image reading support device | |
| JP3460495B2 (en) | Electronic image display method and electronic image display device | |
| JP7301724B2 (en) | MEDICAL IMAGE PROCESSING APPARATUS, METHOD AND PROGRAM | |
| JP2006026396A (en) | Image processing system and method, control program, and storage medium | |
| JP5998903B2 (en) | Radiation image processing apparatus, radiation image processing method, and program |