【発明の詳細な説明】 （産業上の利用分野） 本発明は、画像のパターン認識方法および装置に関す
るものであり、詳細には、被写体の放射線画像の各画素
の濃度を表わす画像データに基づいて、放射線画像内の
所定の画素P0が、該放射線画像を構成する円形パターン
内の画素であるか否かを判定する円形パターン判定方法
および装置に関するものである。Description: TECHNICAL FIELD The present invention relates to a method and an apparatus for recognizing an image pattern. More specifically, the present invention relates to an image pattern recognition method and apparatus based on image data representing the density of each pixel of a radiation image of a subject. The present invention relates to a circular pattern determining method and apparatus for determining whether a predetermined pixel P0 in a radiation image is a pixel in a circular pattern forming the radiation image.
（従来の技術） 記録された放射線画像を読み取って各画素の濃度を表
わす画像データを得、この画像データに適切な画像処理
を施した後、画像を再生記録することは種々の分野で行
なわれている。たとえば、後の画像処理に適合するよう
に設計されたガンマ値の低いＸ線フイルムを用いてＸ縁
画像を記録し、このＸ線画像が記録されたフイルムから
Ｘ線画像を読み取って電気信号（画像データ）に変換
し、この画像データに画像処理を施した後コピー写真等
に可視像として再生することにより、コントラスト，シ
ャープネス，粒状性等の画質性能の良好な再生画像を得
ることが行なわれている（特公昭61−5193号公報参
照）。2. Description of the Related Art In various fields, a recorded radiation image is read, image data representing the density of each pixel is obtained, and after appropriate image processing is performed on the image data, the image is reproduced and recorded. ing. For example, an X-edge image is recorded using an X-ray film having a low gamma value designed to be suitable for later image processing, and the X-ray image is read from the film on which the X-ray image is recorded, and an electric signal ( After converting the image data into image data and performing image processing on the image data and reproducing the image data as a visible image on a copied photograph or the like, a reproduced image having good image quality performance such as contrast, sharpness, and graininess can be obtained. (See Japanese Patent Publication No. 61-5193).
また本願出願人により、放射線（Ｘ線，α線，β線，
γ線，電子線，紫外線等）を照射するとこの放射線エネ
ルギーの一部が蓄積され、その後可視光等の励起光を照
射すると蓄積されたエネルギーに応じて輝尽発光を示す
蓄積性蛍光体（輝尽性蛍光体）を利用して、人体等の被
写体の放射線画像情報をシート状の蓄積性蛍光体に一旦
記録し、この蓄積性蛍光体シートをレーザー光等の励起
光で走査して輝尽発光光を生ぜしめ、得られた輝尽発光
光を光電的に読み取って画像データを得、この画像デー
タに基づき被写体の放射線画像を写真感光材料等の記録
材料、CRT等に可視像として出力させる放射線画像再生
システムがすでに提案されている（特開昭55−12429
号，同56−11395号，同55−163472号，同56−104645
号，同55−116340号等）。In addition, the applicant (X-ray, α-ray, β-ray,
Irradiation with gamma rays, electron beams, ultraviolet rays, etc. accumulates a part of this radiation energy, and then irradiation with excitation light, such as visible light, causes a stimulable phosphor (luminous) to emit stimulated emission according to the accumulated energy. Radiation image information of a subject such as a human body is temporarily recorded in a sheet-shaped stimulable phosphor using a stimulable phosphor, and the stimulable phosphor sheet is scanned with excitation light such as laser light to radiate the light. Generates luminescent light, obtains image data by photoelectrically reading the obtained stimulating luminescent light, and outputs a radiation image of the subject as a visible image to a recording material such as a photographic photosensitive material or a CRT based on this image data. There has already been proposed a radiation image reproduction system for causing the
Nos. 56-11395, 55-163472, 56-104645
No. 55-116340, etc.).
このシステムは、従来の銀塩写真を用いる放射線写真
システムと比較して極めて広い放射線露出域にわたって
画像を記録しうるという実用的な利点を有している。す
なわち、蓄積性蛍光体においては、放射線露光量に対し
て蓄積後に励起によって輝尽発光する発光光の光量が極
めて広い範囲にわたって比例することが認められてお
り、従って種々の撮影条件により放射線露光量がかなり
大幅に変動しても、蓄積性蛍光体シートより放射される
輝尽発光光の光量を読取ゲインを適当な値に設定して光
電変換手段により読み取って電気信号に変換し、この電
気信号を用いて写真感光材料等の記録材料、CRT等の表
示装置に放射線画像を可視像として出力させることによ
って、放射線露光量の変動に影響されない放射線画像を
得ることができる。This system has the practical advantage of being able to record images over a very large radiation exposure area compared to conventional radiographic systems using silver halide photography. That is, in the case of the stimulable phosphor, it has been recognized that the amount of emitted light that is stimulated by excitation after accumulation is proportional to the radiation exposure amount over an extremely wide range. Even if fluctuates considerably, the amount of the stimulating light emitted from the stimulable phosphor sheet is read by the photoelectric conversion means with the reading gain set to an appropriate value and converted into an electric signal. By outputting a radiation image as a visible image on a recording material such as a photographic light-sensitive material or a display device such as a CRT using, a radiation image which is not affected by a change in radiation exposure can be obtained.
上記Ｘ線フイルムや蓄積性蛍光体シート等を用いたシ
ステム、特に人体の医療診断用として構成されたシステ
ムにおいて、近年、単に観察（診断）に適した良好な画
質性能を備えた再生画像を得ることに加えて、画像のパ
ターン認識が行なわれてきている（たとえば特開昭62−
125481号公報参照）。In recent years, in a system using the X-ray film or the stimulable phosphor sheet, particularly a system configured for medical diagnosis of a human body, a reproduced image having good image quality suitable for observation (diagnosis) is simply obtained. In addition, pattern recognition of images has been performed (for example,
No. 125481).
ここで画像のパターン認識とは、画像データに種々の
処理を施すことにより、複雑な放射線画像から目的とす
るパターンを抽出する操作をいい、たとえば人体の胸部
Ｘ線画像のような種々の線状，円形状のパターンの入り
混じった非常に複雑な画像から、たとえば腫瘍に対応す
る略円形の陰影を抽出する操作等をいう。Here, image pattern recognition refers to an operation of extracting a target pattern from a complex radiation image by performing various processes on image data. For example, various linear patterns such as a chest X-ray image of a human body are used. , An operation of extracting, for example, a substantially circular shadow corresponding to a tumor from an extremely complicated image mixed with a circular pattern.
このように複雑な放射線画像（たとえば人体の胸部Ｘ
線画像または乳房Ｘ線画像）において目的とするパター
ン（たとえば腫瘍影）を抽出し、その抽出したパターン
を明示した可視画像を再生表示することにより、観察者
の観察の補助（たとえば医師の診断の補助）を行なわせ
ることができる。Such a complicated radiographic image (for example, the chest X of the human body)
A target pattern (for example, a tumor shadow) is extracted from an X-ray image or a breast X-ray image, and a visible image in which the extracted pattern is specified is reproduced and displayed, thereby assisting the observer in observing (for example, a doctor's diagnosis). Assistance).
（発明が解決しようとする課題） 上記特開昭62−125481号公報には、たとえば人体のＸ
線画像上の腫瘍影（円形パターン）の認識を行なうため
に円形パターンに感度の高い第一のフィルタと線形パタ
ーンに感度の高い第二のフィルタとを併用して、上記Ｘ
線画像上を走査することにより、円形パターン，線形パ
ターンを認識する方法が記載されている。(Problems to be Solved by the Invention) Japanese Patent Application Laid-Open No. 62-125481 discloses, for example,
In order to recognize a tumor shadow (circular pattern) on a line image, a first filter having high sensitivity to a circular pattern and a second filter having high sensitivity to a linear pattern are used in combination with the above X.
A method of recognizing a circular pattern and a linear pattern by scanning a line image is described.
上記第一および第二のフィルタは、それぞれ円形パタ
ーン，線形パターンを洩れなく抽出するとともにそれぞ
れ円形パターン，線形パターンでないものを確実に除外
するという抽出精度が必ずしも十分でないという問題点
がある。これを解決する方法として上記２つのフィルタ
を併用して抽出精度をあげることも考えられるが、２つ
のフィルタを使用するとその分処理が複雑化することに
なる。The first and second filters have a problem that the extraction accuracy of extracting a circular pattern and a linear pattern without omission, respectively, and reliably excluding a non-circular pattern and a linear pattern, respectively, is not always sufficient. As a method of solving this, it is conceivable to increase the extraction accuracy by using the two filters in combination, but the use of two filters complicates the processing.
本発明は、上記事情に鑑み、単一のフィルタで円形パ
ターンの線形パターンとを精度良く分離し、放射線画像
上に指定された所定の画素P0が円形パターン内の画素で
あるか否かを精度良く認識することのできる円形パター
ン判定方法および装置を提供することを目的とするもの
である。In view of the above circumstances, and a linear pattern of the circular pattern accurately separated by a single filter, a given pixel P0 that is specified on the radiographic image whether the pixel in the circular pattern It is an object of the present invention to provide a method and an apparatus for determining a circular pattern that can be accurately recognized.
（課題を解決するための手段） 本発明の円形パターン判定方法および装置は、前述し
たように、被写体の放射線画像の各画素の濃度を表わす
画像データに基づいて、前記放線画像内の所定の画素P0
が、該放射線画像を構成する円形パターン内の画素であ
るか否かを判定する方法および装置に関するものであ
る。(Means for Solving the Problems) As described above, the method and the apparatus for determining a circular pattern according to the present invention use the predetermined pixel in the radiation image based on the image data representing the density of each pixel of the radiation image of the subject. P0
Relates to a method and an apparatus for determining whether or not the pixel is a pixel in a circular pattern constituting the radiation image.
本発明の第一の円形パターン判定方法は、上記方法に
おいて、 前記所定の画素P0を含む中央領域内の多数の画素に対
応する前記画素データの平均的な値Q0と、該中央領域を
取り巻く多数の各周辺領域毎に、該各周辺領域内の多数
の画素に対応する前記画像データの平均的な値Qi（ｉ＝
1,2,…,n）とを求め、 前記中央領域に対応する前記平均的な値Q0と、前記各
周辺領域に対応する前記平均的な値Qi（ｉ＝1,2,…,n）
のそれぞれとの各差分Δｉ（ｉ＝1,2,…,n）を求め、 多数の前記差分Δｉを代表する第一の特性値と多数の
前記差分Δｉのばらつきを表わす第二の特性値とを求
め、 これら第一および第二の特性値との比率を求め、 この比率を所定のしきい値と比較することにより、前
記所定の画素P0が前記円形パターン内の画素であるか否
かを判定することを特徴とするものである。First circular pattern determination method of the present invention, in the above method, the average value Q0 of the pixel data corresponding to the number of pixels in the central area including the predetermined pixel P0, the central region For each of a number of surrounding areas, the average value Qi (i =i) of the image data corresponding to a number of pixels in each of the surrounding areas.
, N), and the average value Q0 corresponding to the central region and the average value Qi (i = 1, 2,. n)
Each differenceΔ i (i = 1,2, ... , n) of the respective determined, the second representing the variation of the first characteristic value and a number of the difference deltai representing the number of the difference deltai of obtains a characteristic value, determine the ratio between these first and second characteristic value, by comparing the ratio with a predetermined threshold, the predetermined pixel P0 is a pixel in the circular pattern Or not.
本発明の円形パターン判定方法は、種々のサイズの円
形パターンや、たとえば楕円形パターン等の種々に変形
した円形パターン等に対処し得るように、以下に示すよ
うに代表値Qiを求めるステップを備えていてもよい。即
ち、本発明の第二の円形パターン判定方法は、 被写体の放射線画像の各画素の濃度を表わす画像デー
タに基づいて、前記放線画像内の所定の画素P0が、該放
射線画像を構成する円形パターン内の画素であるか否か
を判定する円形パターン判定方法において、 前記所定の画素P0を含む中央領域内の多数の画素に対
応する前記画素データの平均的な値Q0と、前記所定の画
素P0から前記放射線画像の周囲に延びる複数の各線分Li
（ｉ＝1,2,…,n）上の、前記所定の画素P0から複数の所
定距離rij（ｊ＝1,2,…,m）離れた各画素Pijを含む各周
辺領域毎に、該各周辺領域内の多数の画素に対応する前
記画像データの平均的な値Qijとを求め、 前記各線分Li毎に、前記平均的な値の代表値Qiを求
め、 前記中央領域に対応する前記平均的な値Q0と、各前記
代表値Qiとの各差分Δｉ（ｉ＝1,2,…,n）を求め、 多数の前記差分Δｉを代表する第一の特性値と多数の
前記差分Δｉのばらつきを表わす第二の特性値とを求
め、 これら第一および第二の特性値との比率を求め、 この比率を所定のしきい値と比較することにより、前
記所定の画素P0が前記円形パターン内の画素であるか否
かを判定することを特徴とする。Circular pattern determination method of the present invention, and a circular pattern of various sizes, for example as may be addressed in a circular pattern or the like various modifications, such as elliptical pattern, determining a representative value Qi as shown below You may have. That is, the second circular pattern determination method according to the present invention is configured such that, based on image data representing the density of each pixel of the radiation image of the subject, a predetermined pixel P0 in the radiation image is a circle forming the radiation image. in the circular pattern determination method of determining whether a pixel in the pattern, the average value Q0 of the pixel data corresponding to the number of pixels in the central area including the predetermined pixel P0, the predetermined A plurality of line segments Li extending from the pixel P0 to the periphery of the radiation image.
(I = 1, 2,..., N), for each peripheral region including each pixel Pij separated from the predetermined pixel P0 by a plurality of predetermined distances rij (j = 1, 2,..., M) to obtains the average value Qij of the image data corresponding to the number of pixels of the respective peripheral region, for each of the line segment Li, obtains a representative value Qi of the average value, the A difference Δi (i = 1, 2,..., N) between the average value Q0 corresponding to the central region and each of the representative values Qi is obtained, and a difference representing the number of the differences Δi is obtained. obtains a second characteristic value representing the variation of one characteristic value and number of the difference deltai, determine the ratio between these first and second characteristic value, comparing this ratio with a predetermined threshold value it makes wherein said predetermined pixels P0 to determine whether the pixels in the circular pattern.
また、本発明の第一の円形パターン判定装置は、前記
第一の円形パターン判定方法を実施する装置であって、 前記所定の画素P0を含む中央領域内の多数の画素に対
応する前記画素データの平均的な値Q0と、該中央領域を
取り巻く多数の各周辺領域毎に、該各周辺領域内の多数
の画素に対応する前記画像データの平均的な値Qi（ｉ＝
1,2,…,n）とを求める平均演算手段、 前記中央領域に対応する前記平均的な値Q0と、前記各
周辺領域に対応する前記平均的な値Qi（ｉ＝1,2,…,n）
のそれぞれとの各差分Δｉ（ｉ＝1,2,…,n）を求める差
分演算手段、 多数の前記差分Δｉを代表する第一の特性値と多数の
前記差分Δｉのばらつきを表わす第二の特性値とを求め
る特性値演算手段、 これら第一および第二の特性値との比率を求める比率
演算手段、および この比率を所定のしきい値と比較することにより、前
記所定の画素P0が前記円形パターン内の画素であるか否
かを判定する判定手段を備えたことを特徴とするもので
ある。Further, the pixel first circular pattern determination device of the present invention is an apparatus for implementing the first circular pattern judging method, which corresponds to the number of pixels in the central area including the predetermined pixel P0 An average value Q0 of the data and, for each of a number of peripheral regions surrounding the central region, an average value Qi of the image data corresponding to a number of pixels in each of the peripheral regions (i =
, N), the average value Q0 corresponding to the central region, and the average value Qi (i = 1, 2) corresponding to each peripheral region. ,…, N)
Each differenceΔ i (i = 1,2, ... , n) of the respective differential calculating means for obtaining a represents the variation of the first characteristic value and a number of the difference deltai representing the number of the differential deltai A characteristic value calculating means for obtaining a second characteristic value; a ratio calculating means for obtaining a ratio between the first and second characteristic values; and comparing the ratio with a predetermined threshold value to obtain the predetermined pixel. is characterized in that P0 is provided with a judging means for judging whether or not pixels in the circular pattern.
さらに、本発明の第二の円形パターン判定装置は、前
記第二の円形パターン判定方法を実施する装置であっ
て、 前記所定の画素P0を含む中央領域内の多数の画素に対
応する前記画像データの平均的な値Q0と、前記所定の画
素P0から前記放射線画像の周囲に延びる複数の各線分Li
（ｉ＝1,2,…,n）上の、前記所定の画素P0から複数の所
定距離rij（ｊ＝1,2,…,m）離れた各画素Pijを含む各周
辺領域毎に、該各周辺領域内の多数の画素に対応する前
記画像データの平均的な値Qijとを求める平均演算手
段、 前記各線分Li毎に、前記平均的な値の代表値Qiを求め
る代表値演算手段、 前記中央領域に対応する前記平均的な値Q0と、各前記
代表値Qiとの各差分Δｉ（ｉ＝1,2,…,n）を求める差分
演算手段、 多数の前記差分Δｉを代表する第一の特性値と多数の
前記差分Δｉのばらつきを表わす第二の特性値とを求め
る特性演算手段、 これら第一および第二の特性値との比率を求める比率
演算手段、および この比率を所定のしきい値と比較することにより、前
記所定の画素P0が前記円形パターン内の画素であるか否
かを判定する判定手段を備えたことを特徴とするもので
ある。Further, the image second circular pattern determination device of the present invention is an apparatus for implementing the second circular pattern judging method, which corresponds to the number of pixels in the central area including the predetermined pixel P0 An average value Q0 of data and a plurality of line segments Li extending from the predetermined pixel P0 to the periphery of the radiation image
(I = 1, 2,..., N), for each peripheral region including each pixel Pij separated from the predetermined pixel P0 by a plurality of predetermined distances rij (j = 1, 2,..., M) Average calculation means for obtaining an average value Qij of the image data corresponding to a large number of pixels in each of the peripheral regions; and for each of the line segments Li , a representative value Qi of the average value A representative value calculating means for calculating; a difference calculating means for calculating each difference Δi (i = 1, 2,..., N) between the average value Q0 corresponding to the central region and each of the representative values Qi ; first characteristic value and the number of the difference deltai second characteristic value and the characteristic calculating means for calculating a representative of the variation of the representative of a number of the difference deltai, these first and ratio of the second characteristic value ratio calculation means for obtaining, and by comparing the ratio with a predetermined threshold, the predetermined pixel P0 is determined manually determines whether the pixels in the circular pattern It is characterized by having a step.
ここで、前記「放射線画像の各画素の濃度」における
「濃度」とは、放射線量に対応した信号値をいい、具体
的には、例えばＸ線写真フィルムの濃度や、蓄積性蛍光
体からの輝尽発光光の強度等をいう。Here, the “density” in the “density of each pixel of the radiation image” refers to a signal value corresponding to a radiation dose, and specifically, for example, the density of an X-ray photographic film, It refers to the intensity of stimulated emission light.
また、前記「平均的な値」とは、典型的には対象とす
る両域内の多数の画像データの相加平均演算により求め
た平均値をいうが、これのみでなく、たとえば相乗平均
値，メジアン値等実質的に平均値を置き換え得る値を含
むものである。ただし、以後簡単のため、これらを代表
する用語として「平均値」を用いる。The “average value” typically refers to an average value obtained by arithmetic averaging of a large number of image data in both regions of interest, but is not limited to this. For example, a geometric mean value, It includes values that can substantially replace the average value, such as the median value. However, for the sake of simplicity, "average value" will be used as a representative term.
また、前記「第一の特性値」とは、前記多数の差分Δ
ｉを代表する値をいい、具体的には上記広義の意味にお
ける多数の差分Δｉの平均値，最大値，最小値等をい
う。また、前記「第二の特性値」とは、前記多数の差分
Δｉのばらつきを表す値をいい、具体的には、該多数の
差分Δｉの分散値，最大値−最小値等をいう。Further, the “first characteristic value” refers to the number of differences Δ
refers to a value representing ai, specifically refers average of a number of difference deltai in the sense of the broad, maximum value, minimum value, and the like. Further, the "second characteristic value" refers to a value representing the variation of the number of difference deltai, specifically, the dispersion value of the difference deltai of said multiple, the maximum value - refers to the minimum value, etc. .
また、前記「複数の各線分Li（ｉ＝1,2,…,n）」の
「複数」（即ちｎの値）は特定の数に制限されるもので
はなく、たとえば必要な判定精度、演算時間等を考慮し
て、任意に定められるものである。The “plurality” (ie, the value of n) of the “plurality of each line segment Li (i = 1, 2,..., N)” is not limited to a specific number. It is arbitrarily determined in consideration of the calculation time and the like.
また、前記「複数の所定距離rij（ｊ＝1,2,…,m）」
のｍの値も上記ｎの値と同様に特定の値に制限されるも
のではなく、また、たとえば放射線画像上の所定の方向
に直径を有する楕円パターンを本発明にいう円形パター
ンに含ませる場合に該所定の方向にはそれと直交する方
向と比較してｍの値を大きくとること等、ｍの値は前記
各線分Li（ｉ＝1,2,…,n）毎に異なっていてもよい。ま
た、各線分Liにおける最小距離riminから最大距離rimax
まで連続的に距離を変えることも、rimin〜rimaxの間に
多数の所定距離rij（ｊ＝1,2,…,m）が含まれると観念
することができ、このように連続的に距離を変える態様
も前記「複数の所定距離rij（ｊ＝1,2,…,m）」に含ま
れる。Further, the “plurality of predetermined distances rij (j = 1, 2,..., M)”
The value of m is not limited to a specific value similarly to the value of n. For example, when an elliptical pattern having a diameter in a predetermined direction on a radiographic image is included in the circular pattern according to the present invention. Even if the value of m is different for each of the line segments Li (i = 1, 2,..., N), for example, the value of m is larger in the predetermined direction than in the direction orthogonal thereto. Good. In addition, the maximum distance rimax from the minimum distance rimin in each line segment Li
It can be thought that the distance is continuously changed up to rimin to rimax , and that a large number of predetermined distances rij (j = 1, 2,...,M ) are included. The plurality of predetermined distances rij (j = 1, 2,..., M) are also included in the above-described embodiment.
さらに前記「代表値」とは、中央領域の平均的なQ0と
各線分Li上の複数の平均的な値Qij（ｊ＝1,2,…,m）と
の差異の有無や程度を表現するに適する値をいい、典型
的には、中央の画像データがその周囲の画像データより
も小さな値を有する円形パターンを対象とする場合はた
とえば各線分Li毎の複数の平均的な値Qijの最大値をい
い、中央の画像データがその周囲の画像データよりも大
きな値を有する円形パターンを対象とする場合はたとえ
ば最小値をいう。ただし、たとえば前記代表値として最
大値が採用される場合において、ある線分Lkに関して前
記各平均的な値Qijが全て負となった場合等その線分Lk
に関して有効な最大値を求め得ない場合に、単純計算に
よる最大値を代表値とすることだけでなく、たとえば0.
0を代表値としてもよいものである。このように本発明
においては、中央領域の平均的な値Q0との差異を有効に
表わしていない代表値が得られた場合に該代表値に代え
て他の値を採用した場合の該他の値も前記代表値に含ま
れる。また、（最大値）−（一定値），（最小値）＋
（一定値）等、実質的に最大値や最小値を表わしている
値も上記代表値に含まれる。このように他の値を採用し
得ることは、本発明の第一の円形パターン判定方法およ
び第一の円形パターン判定装置についても同様であり、
求められた差分Δｉが中央領域の平均的な値Q0との差異
を有効に表わしていない場合、この差分Δｉに代えてた
とえば0.0等を採用してもよい。Further, the "representative value" refers to the presence or absence and degree of a difference between the average Q0 in the central region and a plurality of average values Qij (j = 1, 2,..., M) on each linesegment Li. refers to a value suitable to represent, typically, a central image data is a plurality of the average case, for example every line segment Li to target circular pattern having a smaller value than the image data of the surrounding of It refers to the maximum value of the value Qij , and for example, it refers to the minimum value when the central image data targets a circular pattern having a larger value than the surrounding image data. However, for example, when the maximum value is adopted as the representative value, when the average value Qij is all negative with respect to a certain line segment Lk , the line segment Lk
When a valid maximum value cannot be obtained for, not only the maximum value obtained by simple calculation is used as the representative value but also, for example, 0.
0 may be set as the representative value. As described above, in the present invention, when a representative value that does not effectively represent the difference from the average value Q0 of the central region is obtained, the other value when another value is adopted instead of the representative value is obtained. Is also included in the representative value. Also, (maximum value)-(constant value), (minimum value) +
A value substantially representing the maximum value or the minimum value, such as (constant value), is also included in the representative value. It is the same for the first circular pattern determination method and the first circular pattern determination device of the present invention that the other values can be adopted as described above,
If the obtained difference deltai is not effectively represent the difference between the average value Q0 of the central region may be employed, for example 0.0, such as in place of the difference deltai.
また、実質的に上記演算を行なう方法，装置であれば
本発明に含まれるものであり、たとえば前記代表値とし
て最大値が採用され、差分演算式として Δｉ＝Qi−Q0 …（１） が採用される場合において、各線分Li毎の複数の平均的
な値Qijの最大値Qiを求めた後（１）式に従って差分演
算を行なうことに代えて、複数の平均的な値Qijそれぞ
れの差分演算Δij＝Qij−Q0を行ない、しかる後に各線
分Li毎の差分演算Δijの最大値を求めても、その最大値
は上記（１）と同一となるが、たとえばこのような演算
も本発明の内包されるものである。Further, any method or apparatus for performing the above-described calculation is included in the present invention. For example, a maximum value is adopted as the representative value, and Δi = Qi −Q0 (1) in the case where) is employed, instead of performing the difference operation according to each line segment after the maximum value Qi of a plurality of average values Qij for each Li (1) formula, a plurality of average performs value Qij each difference calculation Δij = Qij -Q0, whereafter even the maximum value of the difference calculation deltaij for each line segment Li, the maximum value is the same as above (1) However, for example, such an operation is also included in the present invention.
（作用） 本発明の円形パターン判定方法および装置は、上記の
構成を有し、第一の特性値と第二の特性値の比較により
所定の画素P0が円形パターン内の画素であるか否かを判
定するようにしたため、上記構成の単一のフィルタのみ
で円形パターンと線形パターンとが精度良く分離され、
従って精度の良い判定が行なわれる。The circular pattern judging method and apparatus (action) the present invention has the configuration described above, whether a predetermined pixel P0 by comparing the first characteristic value and the second characteristic value is a pixel in a circular pattern In order to determine whether the circular pattern and the linear pattern are accurately separated by only a single filter of the above configuration,
Therefore, an accurate determination is made.
また、本発明の第二の円形パターン判定方法および第
二の円形パターン判定装置は、所定の画素P0から周囲に
延びる複数の各線分Li毎に該各線分Liに対応する複数の
周辺領域の平均的な値Qijの代表値Qiを求めてこれらの
代表値Qiを用いるようにしたため、選定された複数の所
定距離rij（ｊ＝1,2,…ｍ）に応じて小さな円形パター
ンから大きな円形パターンまでの種々の円形パターンが
混在していても、またたとえば楕円パターンのような変
形した円形パターンが種々存在していてもこれらに対し
て十分な判定能力を備えたフィルタが実現される。The second circular pattern judging method and the second circular pattern determination device of the present invention, a plurality of peripheral corresponding to respective segment Li in each of a plurality of line segments Li extending around a predetermined pixel P0 Since the representative values Qi of the average values Qij of the area are obtained and these representative values Qi are used, according to the plurality of selected predetermined distances rij (j = 1, 2,... M). Even if various circular patterns from a small circular pattern to a large circular pattern are mixed, and even if there are various deformed circular patterns such as an elliptical pattern, a filter having a sufficient judgment ability for them. Is realized.
（実 施 例） 以下、本発明の実施例について、図面を参照して説明
する。尚、ここでは前述した蓄積性蛍光体シートを用
い、人体の乳房内に典型的には略球形として生じる腫瘍
の陰影を抽出する例について説明する。この腫瘍影は典
型的には可視画像上では周囲と比べ白っぽい（濃度が低
い）略円形パターンとして現われる。(Example) Hereinafter, an example of the present invention will be described with reference to the drawings. Here, an example will be described in which the above-described stimulable phosphor sheet is used to extract a shadow of a tumor which typically occurs as a substantially spherical shape in the breast of a human body. The tumor shadow typically appears as a whitish (low density) substantially circular pattern on the visible image compared to the surroundings.
第３図は、Ｘ線撮影装置の一例の概略図である。 FIG. 3 is a schematic view of an example of an X-ray imaging apparatus.
このＸ線撮影装置10のＸ線源11からＸ線12が人体13の
乳房13aに向けて照射され、人体13を透過したＸ線12aが
蓄積性蛍光体シート14に照射されることにより、人体の
乳房13aの透過Ｘ線画像がシート14に蓄積記録される。X-rays 12 are radiated from an X-ray source 11 of the X-ray imaging apparatus 10 toward a breast 13a of a human body 13, and the X-rays 12a transmitted through the human body 13 are radiated to a stimulable phosphor sheet 14, whereby a human body is irradiated. The transmission X-ray image of the breast 13a is stored and recorded on the sheet 14.
第４図は、Ｘ線画像読取装置の一例と、本発明の円形
パターン判定装置の一実施例を内包するコンピュータシ
ステムとを表わした斜視図である。FIG. 4 is a perspective view showing an example of an X-ray image reading device and a computer system including an embodiment of the circular pattern determination device of the present invention.
Ｘ線画像が記録された蓄積性蛍光体シート14がＸ線画
像読取装置20の所定位置にセットされる。この所定位置
にセットされた蓄積性蛍光体シート14は、モータ21によ
り駆動されるエンドレスベルト等のシート搬送手段22に
より、矢印Ｙ方向に搬送（副走査）される。一方、レー
ザー光源23から発せられた光ビーム24はモータ25により
駆動され矢印方向に高速回転する回転多面鏡26によって
反射偏向され、ｆθレンズ等の集束レンズ27を通過した
後、ミラー28により光路を変えて前記シート14に入射し
副走査の方向（矢印Ｙ方向）と略垂直な矢印Ｘ方向に主
走査する。シート14の励起光24が照射された箇所から
は、蓄積記録されているＸ線画像情報に応じた光量の輝
尽発光光29が発散され、この輝尽発光光29は光ガイド30
によって導かれ、フォトマルチプライヤ（光電子増倍
管）31によって光電的に検出される。上記光ガイド30は
アクリル板等の導光性材料を成形して作られたものであ
り、直線状をなす入射端面30aが蓄積性蛍光体シート14
上の主走査線に沿って延びるように配され、円環状に形
成された射出端面30bにフォトマルチプライヤ31の受光
面が結合されている。入射端面30aから光ガイド30内に
入射した輝尽発光光29は、該光ガイド30の内部を全反射
を繰り返して進み、射出端面30bから射出してフォトマ
ルチプライヤ31に受光され、Ｘ線画像を表わす輝尽発光
光29がフォトマルチプライヤ31によって電気信号に変換
される。The stimulable phosphor sheet 14 on which the X-ray image is recorded is set at a predetermined position of the X-ray image reading device 20. The stimulable phosphor sheet 14 set at this predetermined position is conveyed (sub-scanning) in the direction of arrow Y by sheet conveying means 22 such as an endless belt driven by a motor 21. On the other hand, a light beam 24 emitted from a laser light source 23 is reflected and deflected by a rotating polygon mirror 26 driven by a motor 25 and rotated at a high speed in the direction of the arrow, and after passing through a focusing lens 27 such as an fθ lens, the optical path is changed by a mirror 28. Instead, the light is incident on the sheet 14 and is main-scanned in an arrow X direction substantially perpendicular to the sub-scanning direction (arrow Y direction). From the portion of the sheet 14 irradiated with the excitation light 24, stimulated emission light 29 of an amount corresponding to the accumulated and recorded X-ray image information is diverged, and the stimulated emission light 29 is
And is photoelectrically detected by a photomultiplier (photomultiplier tube) 31. The light guide 30 is formed by molding a light-guiding material such as an acrylic plate, and a linear incident end face 30a is formed on the stimulable phosphor sheet 14.
The light receiving surface of the photomultiplier 31 is coupled to the emission end face 30b which is arranged so as to extend along the upper main scanning line and is formed in an annular shape. The stimulated emission light 29 that has entered the light guide 30 from the incident end face 30a travels through the inside of the light guide 30 by repeating total reflection, exits from the emission end face 30b, is received by the photomultiplier 31, and is received by the X-ray image. Is converted into an electric signal by the photomultiplier 31.
フォトマルチプライヤ31から出力されたアナログ出力
信号S0は対数増幅器32で対数的に増幅され、A/D変換器3
3でディジタル化され、電気信号としての画像データ
（各画素の濃度を表わす画像データ）S1が得られる。The analog output signal S0 output from the photomultiplier 31 is logarithmically amplified by a logarithmic amplifier 32, and the A / D converter 3
Image data (image data representing the density of each pixel) S1 which is digitized in step 3 as an electric signal is obtained.
得られた画像データS1は、コンピュータシステム40に
入力される。このコンピュータシステム40は、本発明の
円形パターン判定装置の一例を内包する構成を有するも
のであり、CPUおよび内部メモリが内蔵された本体部41,
補助メモリとしてのフロッピィディスクが挿入されドラ
イブされるドライブ部42,オペレータがこのコンピュー
タシステム40に必要な指示等を入力するためのキーボー
ド43および必要な情報を表示するためのCRTディスプレ
イ44から構成されている。The obtained image data S1 is input to the computer system 40. The computer system 40 has a configuration including an example of the circular pattern determination device of the present invention, and includes a main body 41 having a built-in CPU and an internal memory.
It comprises a drive unit 42 into which a floppy disk as an auxiliary memory is inserted and driven, a keyboard 43 for an operator to input necessary instructions and the like to the computer system 40, and a CRT display 44 for displaying necessary information. I have.
コンピュータシステム40に入力された画像データS1に
基づいて、Ｘ線画像の各画素を所定の画素P0として、該
所定の画素P0が腫瘍影（円形パターン）内の画素である
か否かが判定され、これをＸ線画像全体にわたって繰り
返すことにより、該Ｘ線画像上の腫瘍影が抽出される。Based on the image data S1 input to the computer system 40, as the predetermined pixels P0 to each pixel of the X-ray image, whether the predetermined pixel P0 is a pixel in the tumor shadow (circular pattern) It is determined, and this is repeated over the entire X-ray image, whereby a tumor shadow on the X-ray image is extracted.
第１図は、上記腫瘍影を抽出するための本発明の一実
施例としての実空間フィルタを説明するために、Ｘ線画
像上の所定の画素P0を中心に該画像上に仮想的に描いた
図である。該所定の画素P0が上記腫瘍影内の画素である
か否かが判定される。ここに示すようなフィルタを用い
てＸ線画像上を走査することにより、Ｘ線画像上の腫瘍
影が抽出される。Figure 1, in order to explain the actual spatial filter as an embodiment of the present invention for extracting the tumor shadow, virtually around a predetermined pixel P0 on the X-ray image on the image FIG. The predetermined pixel P0 is whether a pixel of said tumor Kagenai is determined. By scanning the X-ray image using the filter as shown here, a tumor shadow on the X-ray image is extracted.
第１図に示すように、Ｘ線画像内の所定の画素P0から
該Ｘ線画像の周囲に延びる複数（ここで８本）の線分Li
（ｉ＝1,2,…,8）を想定し、さらに所定の画素P0を中心
とした半径ｒの円Ｒを想定する。また所定の画素P0を含
む中心領域Q0と、線分Li（ｉ＝1,2,…,8）のそれぞれと
円Ｒとの交点の各画素Pi（ｉ＝1,2,…,8）を含む各周辺
領域Qiを考える。尚、上記半径ｒ、中央領域Q0,各周辺
領域Qiの面積、および想定する周辺領域の数等は、対象
とする円形パターンの大きさ、判定精度、演算速度等を
考慮して適切に定められる。また、本実施例では、所定
の画素Q0から等距離ｒだけ離れた各画素Piを想定してい
るが、たとえば第１図のＸ方向に直径をもつ楕円パター
ンと対象とする場合、画素P1,P5として画素P0から遠距
離にある画素を選択する等、各画素Pi毎に所定の画素P0
からの距離が異なっていてもよい。As shown in FIG. 1, a line segment of a plurality (here eight) from predetermined pixels P0 in the X-ray image extending around the X-ray image Li
(I = 1,2, ..., 8 ) assume, further assuming a circle R of radius r centered at the given pixel P0. Also, each pixel Pi (i = 1, 2,...) At the intersection of the circle R with the center area Q0 including the predetermined pixel P0 and each of the line segments Li (i = 1, 2,..., 8). consider each peripheral area Qi, including 8). Incidentally, the radius r, the center region Q0, Suto of each peripheral area of the region Qi, and assuming that the peripheral region, the size of the circular pattern of interest, determination accuracy, to properly consider the operation speed, etc. Determined. In this embodiment, each pixelPi is assumed to be equidistant from the predetermined pixel Q0 by an equal distance r. For example, when an elliptical pattern having a diameter in the X direction in FIG. P1, P equal to the pixel P0 as5 selects a pixel in the far, each pixel Pi predetermined pixel P per0
May be different.
上記のようにして想定した中央領域Q0および各周辺領
域Qi内の多数の画素に対応する多数の画像データの平均
値Q0,Qi（ｉ＝1,2,…,8）が求められる。尚、ここでは
簡単のため、各領域Q0,Qi（ｉ＝1,2,…,8）を指す記号
と該各領域内の画像データの平均値を指す記号とで同一
の番号を用いている。Average values Q0 , Qi (i = 1, 2,..., 8) of a large number of image data corresponding to a large number of pixels in the central region Q0 and the peripheral regions Qi assumed as described above are obtained. Can be Here, for the sake of simplicity, the same number is used for a symbol indicating each area Q0 , Qi (i = 1, 2,..., 8) and a symbol indicating the average value of image data in each area. ing.
次に中央領域の平均値Q0と各周辺領域の平均値Qiのそ
れぞれと各差分Δｉ（ｉ＝1,2,…,8）が Δｉ＝Qi−Q0 …（１） として求められる。Next, the average value Q0 of the central region and the average value Qi of each peripheral region and each difference Δi (i = 1, 2,..., 8) are defined as Δi = Qi −Q0 (1). Desired.
本実施例では、これらの差分Δｉ（ｉ＝1,2,…,8）を
代表する第一の特性値、およびこれらの差分Δｉ（ｉ＝
1,2,…,8）のばらつきを表わす第二の特性値としてそれ
ぞれこれらの差分Δｉの平均値が採用され、これら平均値，分散σ２が求められる。In this embodiment, the first characteristic value representing the difference Δi (i = 1, 2,..., 8) and the difference Δi (i =
1,2, ..., a second average value of the differential deltai respectively as the characteristic value representing the variation of 8) Are adopted, and the average value and the variance σ2 are obtained.
次に、これら平均値，分散σ２の比率Ｃが、として求められ、この比率C1が所定のしきい値Th1と比
較され、C1≧Th1の場合は比較的平均値が大きく分散
σ２が小さいため所定の画素P0が腫瘍影内の画素である
と判定され、C1＜Th1の場合は所定の画素P0は腫瘍影外
の画素であると判定される。Next, the ratio C of these average value and variance σ2 is Obtained as this ratio C1 is compared with a predetermined threshold value Th1, a given pixel P0 for a relatively average value large variance sigma2 is small in the case of C1 ≧ Th1 is the pixel tumor Kagenai is determined to be the predetermined pixel P0 in the case of C1 <Th1 is determined as a pixel of tumor Kagegai.
第２図は、本発明の他の実施例としての空間フィルタ
を説明するために、Ｘ線画像上の所定の画素P0を中心に
該画像上に仮想的に描いた図である。FIG. 2, for explaining a spatial filter as another embodiment of the present invention, is a diagram depicting virtually on the image around a predetermined pixel P0 on the X-ray image.
第２図に示すように、Ｘ線画像内の所定の画素P0から
該Ｘ線画像の周囲に延びる複数（ここでは８本）の線分
Li（ｉ＝1,2,…,8）を想定し、さらに所定の画素P0を中
心とした、それぞれ半径r1,r2,r3の３つの円Rj（ｊ＝1,
2,3）を想定する。所定の画素P0を含む中央領域をQ0と
し、各線分Liと各円Rjとの各交点に位置する各画素Pij
（第２図にはP11,P12,P13,P51,P52,P53について記号を
示してある。）を含む各周辺領域をPij（ｉ＝1,2,…,8;
j＝1,2,3）（ただし第２図には、明示的にはQ0およびQ
11,Q12,Q13,Q51,Q52,Q53のみ示してある。）とする。As shown in FIG. 2, a plurality of (eight in this case) line segments extending from a predetermined pixel P0 in the X-ray image to the periphery of the X-ray image
Assuming Li (i = 1, 2,..., 8), three circles Rj (j = 1,2 ) with respective radii r1 , r2 , r3 centered on a predetermined pixel P0
2,3) is assumed. The central area including the predetermined pixel P0 is defined as Q0, and each pixel Pij located at each intersection of each line segmentLi and each circle Rj
(In FIG. 2, symbols are shown for P11 , P12 , P13 , P51 , P52 , and P53. ) Each peripheral region is represented by Pij (i = 1, 2,..., 8;
j = 1,2,3) (However, in FIG. 2, Q0 and Q
Only11 , Q12 , Q13 , Q51 , Q52 and Q53 are shown. ).
この各領域Q0およびQij（ｉ＝1,2,…,8;j＝1,2,3）毎
に、該各領域Q0,Qij内の多数の各画素に対応する多数の
各画像データの平均値Q0,Qij（ｉ＝1,2,…,8;j＝1,2,
3）が求められる。尚、ここでも簡単のため、各領域Q0,
Qij（ｉ＝1,2,…,8;j＝1,2,3）を指す記号と該各領域内
の画像データの平均値を指す記号とで同一の記号を用い
ている。For each of the regions Q0 and Qij (i = 1, 2,..., 8; j = 1, 2, 3), a large number of pixels corresponding to a large number of pixels in each of the regions Q0 , Qij Average values Q0 , Qij (i = 1, 2,..., 8; j = 1, 2,
3) is required. Here, for simplicity, each region Q0 ,
The same symbol is used for the symbol indicating Qij (i = 1, 2,..., 8; j = 1, 2, 3) and the symbol indicating the average value of the image data in each area.
次に中央領域の平均値Q0と各周辺領域の平均値Qijの
それぞれとの各差分Δij（ｉ＝1,2,…,8;j＝1,2,3）が Δij＝Qij−Q0 …（２） として求められ、さらに各線分Li毎に、差分Δijの最大
値Δｉが求められる。本実施例ではこの最大値Δｉが本
発明にいう代表値である。Next, each difference Δij (i = 1, 2,..., 8; j = 1, 2, 3) between the average value Q0 of the central region and the average value Qij of each peripheral region is Δij = Qij -Q0 ... determined as (2), further each line segment Li, the maximum value deltai of the difference deltaij is calculated. In this embodiment the maximum value deltai is the representative value according to the present invention.
次に、最大値Δｉ（ｉ＝１〜８）を代表する第一の特
性値Ｕと最大値Δｉ（ｉ＝１〜８）のばらつきを表わす
第二の特性値Ｖとが求められる。このために、まず以下
の演算式に従って各特性値U1〜U4,V1〜V4が求められ
る。Next, a second characteristic value V representing the variation of the maximum value Δi (i = 1~8) first characteristic value representative of the U and the maximum value Δi (i = 1~8) is obtained. For this purpose, first, the characteristic values U1 to U4 and V1 to V4 are obtained according to the following arithmetic expressions.
U1＝（Δ１＋Δ２＋Δ５＋Δ６）/4 …（３） U2＝（Δ２＋Δ３＋Δ６＋Δ７）/4 …（４） U3＝（Δ３＋Δ４＋Δ７＋Δ８）/4 …（５） U4＝（Δ４＋Δ５＋Δ８＋Δ１）/4 …（６） V1＝U1/U3 …（７） V2＝U2/U4 …（８） V3＝U3/U1 …（９） V4＝U4/U2 …（10） ここで、たとえば（３）式に従って特性値U1を求める
場合について説明すると、隣接する２つの領域（Δ１と
Δ２、またはΔ５とΔ６）について加算することは平滑
化を意味し、画素P0を挾んだ互いに反対側の領域（Δ１
＋Δ２とΔ５＋Δ６）について加算することはたとえば
濃度勾配のある領域に腫瘍影があってもその腫瘍影を検
出することができるようにするためである。U 1 = (Δ 1 + Δ 2 + Δ 5 + Δ 6) / 4 ... (3) U 2 = (Δ 2 + Δ 3 + Δ 6 + Δ 7) / 4 ... (4) U 3 = (Δ 3 + Δ 4 + Δ 7 + Δ 8 ) / 4 ... (5) U 4 = (Δ 4 + Δ 5 + Δ 8 + Δ 1) / 4 ... (6) V 1 = U 1 / U 3 ... (7) V 2 = U 2 / U 4 ... (8) V3 = U3 / U1 (9) V4 = U4 / U2 (10) Here, for example, a case where the characteristic value U1 is obtained according to the equation (3) will be described. delta1 and delta2 or delta adding about5 and delta6) means smoothing, region of opposite sides sandwiching the pixel P0 (delta1,
The addition of + Δ2 and Δ5 + Δ6 ) is performed so that even if there is a tumor shadow in an area having a concentration gradient, the tumor shadow can be detected.
また、たとえば（７）式に従って特性値V1を求める場
合について説明すると、特性値U1と特性値U3とは互いに
直交する方向について求めた特性値であり、したがって
第２図に示す腫瘍影７が円形であればV1≒1.0となり円
形から外れる場合、即ち画素P0が肋骨影のように直線状
の陰影内にある場合はV1は1.0から外れることになる。Further, for example, (7) will be described for obtaining the characteristic value V1 in accordance with formula, is a characteristic value determined for directions perpendicular to each other and the characteristic value U1 and the characteristic value U3, thus the tumor shadow as shown in FIG. 2 If 7 is circular, V1 ≒ 1.0 and deviates from the circle, ie, if pixel P0 is within a linear shadow such as a rib shadow, V1 will deviate from 1.0.
上記差分の最大値Δ１（ｉ＝１〜８）を代表する第一
の特性値Ｕとしては、U1〜U4の最大値 Ｕ＝MAX（U1,U2,U3,U4） …（11） が採用され、上記差分の最大値Δ１（ｉ＝１〜８）のば
らつきを表わす第二の特性値Ｖとしては、V1〜V4の最大
値 Ｖ＝MAX（V1,V2,V3,V4） …（12） が採用される。このようにして第一および第二の特性値
U,Vが求められると、所定の画素P0が腫瘍影内の画素で
あるか否かを判定するための特性値C2として、これら第
一および第二の特性値の比率が採用され、この特性値C2が所定のしきい値Th2と比較
され、C2≧Th2であるか、C2＜Th2であるかにより、画素
P0がそれぞれ腫瘍影内の画素であるか否かが判定され
る。As the first characteristic value U representing the maximum value Δ1 (i = 1 to 8) of the difference, the maximum value U1 to U4 U = MAX (U1 , U2 , U3 , U4 ) (11) is adopted, and as the second characteristic value V representing the variation of the maximum value Δ1 (i = 1 to 8) of the difference, the maximum value V1 to V4 V = MAX (V1 ,V 2, V 3, V 4 ) ... (12) is employed. Thus, the first and second characteristic values
U, the V is obtained, as a characteristic value C2 for a given pixel P0 to determine whether the pixel of tumor Kagenai, the ratio of these first and second characteristic values There is employed, the characteristic value C2 is compared with a predetermined threshold value Th2, or a C2 ≧ Th2, depending on whether it is C2 <Th2, the pixel
It is determined whether each of P0 is a pixel in the tumor shadow.
以上、例を示したように、第４図に示すコンピュータ
システム40の候補抽出手段１では実空間フィルタを用い
てＸ線画像上を走査することにより、腫瘍影と考えられ
る円形パターンの抽出が行なわれる。As described above, as shown in the example, the candidate extracting means 1 of the computer system 40 shown in FIG. 4 extracts a circular pattern considered as a tumor shadow by scanning the X-ray image using a real space filter. It is.
尚、上記各フィルタ例においては、第１図，第２図に
示すように８本の線分L1〜L8上の画素Pi（Pij）を含む
各周辺領域Qi（Qij）に対応する画像データの平均値Qi
（Qij）を用いたが、この線分は８本である必要はな
く、たとえば16本等であってもよいことはもちろんであ
る。また、第２図を用いて説明した実施例ではr1,r2,r3
の３つの距離について演算を行なったが、これについて
も３つの距離に限るものでもなく、種々の大きさの腫瘍
影をさらに精度よく抽出するために、距離をr1からr3ま
で連続的に変えて演算を行なってもよい。In each of the above filter examples, as shown in FIGS. 1 and 2, each peripheral area Qi (Qij ) including pixels Pi (Pij ) oneight line segments L1 to L8. Average value Qi of the image data corresponding to
Although (Qij ) is used, it is needless to say that this line segment does not need to be eight, and may be, for example, sixteen. In the embodiment described with reference to FIG.2 , r1 , r2 , r3
The calculation was performed for the three distances, but this is not limited to the three distances as well, and in order to more accurately extract tumor shadows of various sizes, the distances were continuously changed from r1 to r3. The operation may be changed.
以上の実施例は、蓄積性蛍光体を用いて得られた人体
の乳房Ｘ線画像に典型的には円形として現われる腫瘍影
を抽出する例であるが、本発明は腫瘍影の抽出に限られ
るものではなく、また乳房Ｘ線画像に限られるものでも
なく胸部Ｘ線画像等にも適用することができ、さらに蓄
積性蛍光体を用いるシステムに限られるものでもなく、
被写体の放射線画像を表わす画像データに基づいて該放
射線画像内の所定の画素P0が、該放射線画像を構成する
円形パターン内の画素であるか否かを判定する際に広く
用い得るものである。The above embodiment is an example of extracting a tumor shadow that typically appears as a circle in a human breast X-ray image obtained using a stimulable phosphor, but the present invention is limited to the extraction of a tumor shadow. It is not limited to a breast X-ray image and can be applied to a chest X-ray image and the like, and is not limited to a system using a stimulable phosphor,
It can be widely used in determining whether or not a predetermined pixel P0 in the radiation image is a pixel in a circular pattern forming the radiation image based on image data representing the radiation image of the subject. .
また、本発明は上記のようなフィルタを用いて放射線
画像上を走査することはかならずしも必要ではなく、た
とえばオペレータが放射線画像上の所定点をマニュアル
指定し、その所定点が円形パターン内の点であるか否か
を判定してもよい。In the present invention, it is not always necessary to scan on a radiation image using the above-described filter.For example, an operator manually designates a predetermined point on a radiation image, and the predetermined point is a point in a circular pattern. It may be determined whether or not there is.
（発明の効果） 以上詳細に説明したように、本発明の円形パターン判
定方法および装置は、中央領域の平均的な値Q0と各周辺
領域の平均的な値Qiとの各差分Δｉを求め、これらの差
分Δｉを代表する第一の特性値とこれらの差分Δｉのば
らつきを表わす第二の特性値とを求めて、これらの特性
値の比率の大小に応じて所定の画素P0が前記円形パター
ン内の画素であるか否かを判定するようにしたため、上
記構成の単一のフィルタのみで精度の良い判定が行なわ
れる。(Effects of the Invention) As described above in detail, the circular pattern determination method and apparatus of the present invention provide each difference Δi between the average value Q0 of the central region and the average value Qi of each peripheral region. asking, seeking a second characteristic value representing the variation of the first characteristic value and the difference between them deltai representing these differences deltai, predetermined pixels in accordance with the magnitude of the ratio of these characteristic values since P0 was made to determine whether the pixels in the circular pattern, may determine accurate with only a single filter having the above structure is performed.
また、各線分Li毎に複数の平均的な値Qijを求め、こ
れらの代表値Qiを求め、上記平均的な値Qiに代えてこの
代表値Qiを用いるようにすると、小さな円形パターンか
ら大きな円形パターンまでの種々の円形パターンが混在
していても、またたとえば楕円パターンのような変形し
た円形パターンが種々存在していても、これらに対して
十分な判定能力を備えることができる。Further, obtains a plurality of average values Qij for each line segment Li, obtains these representative values Qi, when so using the representative value Qi in place of the average value Qi, small Even if various circular patterns from a circular pattern to a large circular pattern coexist, and even if there are various deformed circular patterns such as an elliptical pattern, for example, it is necessary to provide sufficient judgment ability for these. it can.
第１図は、円形の腫瘍影を抽出するための、本発明の一
実施例としての空間フィルタを説明するために、Ｘ線画
像上の所定の画素P0を中心に該画像上に仮想的に描いた
図、 第２図は、円形の腫瘍影を抽出するための、本発明の他
の実施例としての空間フィルタを説明するために、Ｘ線
画像上の所定の画素P0を中心に該画像上に仮想的に描い
た図、 第３図は、Ｘ線画像撮影装置の一例の概略図、 第４図は、Ｘ線画像読取装置の一例と、本発明の円形パ
ターン判定装置の一実施例を内包するコンピューターシ
ステムとを表わした斜視図である。 10……Ｘ線撮影装置、14……蓄積性蛍光体シート 20……Ｘ線画像読取装置 23……レーザ光源、26……回転多面鏡 29……輝尽発光光、30……光ガイド 31……フォトマルチプライヤ 40……コンピュータシステムThe first figure is for extracting circular tumor shadow to illustrate the spatial filter as an embodiment of the present invention, virtually on the image around a predetermined pixel P0 on the X-ray image FIG. 2 is a diagram centered on a predetermined pixel P0 on an X-ray image for explaining a spatial filter as another embodiment of the present invention for extracting a circular tumor shadow. FIG. 3 is a diagram schematically drawn on the image, FIG. 3 is a schematic diagram of an example of an X-ray image photographing device, FIG. 4 is an example of an X-ray image reading device, and FIG. 1 is a perspective view showing a computer system including an embodiment. 10 X-ray imaging device, 14 stimulable phosphor sheet 20 X-ray image reading device 23 Laser light source 26 Rotating polygon mirror 29 Photostimulated emission light 30, Light guide 31 …… Photomultiplier 40 …… Computer system
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1162909AJP2631748B2 (en) | 1989-06-26 | 1989-06-26 | Circular pattern determination method and apparatus |
| US07/542,487US5033100A (en) | 1989-06-26 | 1990-06-22 | Method and apparatus for classifying picture elements in radiation images |
| DE69033380TDE69033380T2 (en) | 1989-06-26 | 1990-06-26 | Method and device for classifying picture elements in radiation images |
| EP90112147AEP0405457B1 (en) | 1989-06-26 | 1990-06-26 | Method and apparatus for classifying picture elements in radiation images |
| EP95116196AEP0702321B1 (en) | 1989-06-26 | 1990-06-26 | Method and apparatus for classifying picture elements in radiation images |
| DE69033403TDE69033403T2 (en) | 1989-06-26 | 1990-06-26 | Method and device for classifying picture elements in radiation images |
| DE69033666TDE69033666T2 (en) | 1989-06-26 | 1990-06-26 | Method and device for classifying picture elements in radiation images |
| EP95116193AEP0698855B1 (en) | 1989-06-26 | 1990-06-26 | Method and apparatus for classifying picture elements in radiation images |
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1162909AJP2631748B2 (en) | 1989-06-26 | 1989-06-26 | Circular pattern determination method and apparatus |
| Publication Number | Publication Date |
|---|---|
| JPH0328982A JPH0328982A (en) | 1991-02-07 |
| JP2631748B2true JP2631748B2 (en) | 1997-07-16 |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1162909AExpired - Fee RelatedJP2631748B2 (en) | 1989-06-26 | 1989-06-26 | Circular pattern determination method and apparatus |
| Country | Link |
|---|---|
| JP (1) | JP2631748B2 (en) |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62125481A (en)* | 1985-11-26 | 1987-06-06 | インターナショナル・ビジネス・マシーンズ・コーポレーション | Pattern recognition equipment |
| JP2522277B2 (en)* | 1987-01-09 | 1996-08-07 | 三菱電機株式会社 | Image recognition device |
| Publication number | Publication date |
|---|---|
| JPH0328982A (en) | 1991-02-07 |
| Publication | Publication Date | Title |
|---|---|---|
| JP2574181B2 (en) | Abnormal shadow detector | |
| JPH06304159A (en) | Method for recognizing picture photographing position | |
| JP2867055B2 (en) | Edge determination method and apparatus | |
| EP0405457B1 (en) | Method and apparatus for classifying picture elements in radiation images | |
| JP2631737B2 (en) | Method for determining image points in subject image | |
| JP2582669B2 (en) | Device for determining benign / malignant breast mass shadow | |
| JP2961437B2 (en) | Image output method | |
| JPH0398174A (en) | Method and device for process of radiation image | |
| JP2582667B2 (en) | Linear pattern width calculator | |
| JP2631748B2 (en) | Circular pattern determination method and apparatus | |
| JP2582664B2 (en) | Pattern recognition device | |
| JP2582666B2 (en) | Abnormal shadow detector | |
| JP2582665B2 (en) | Abnormal shadow detector | |
| JP2867067B2 (en) | Pattern determination method and apparatus | |
| JP2631747B2 (en) | Circular pattern determination method and apparatus | |
| JP2707355B2 (en) | Image recording position recognition device | |
| JP2867036B2 (en) | Abnormal shadow detector | |
| JP3165538B2 (en) | Radiation image alignment method and apparatus | |
| JP2987634B2 (en) | Abnormal shadow judgment device | |
| JP2582660B2 (en) | Abnormal shadow detector | |
| JP3165530B2 (en) | Radiation image reading condition determination method and apparatus | |
| JP2574185B2 (en) | Method for determining reading conditions and / or image processing conditions in breast radiation image with chest wall | |
| JPH04156829A (en) | Pattern recognition method and device | |
| JP2561158B2 (en) | Pattern recognizer | |
| JPH0271247A (en) | Irradiation field recognizing method |
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees | Free format text:JAPANESE INTERMEDIATE CODE: R250 | |
| LAPS | Cancellation because of no payment of annual fees |