JPS594745B2 - 細胞診装置 - Google Patents

細胞診装置

Info

Publication number
JPS594745B2
JPS594745B2 JP51121658A JP12165876A JPS594745B2 JP S594745 B2 JPS594745 B2 JP S594745B2 JP 51121658 A JP51121658 A JP 51121658A JP 12165876 A JP12165876 A JP 12165876A JP S594745 B2 JPS594745 B2 JP S594745B2
Authority
JP
Japan
Prior art keywords
cells
cell
concentration
histogram
density
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired
Application number
JP51121658A
Other languages
English (en)
Other versions
JPS5346797A (en
Inventor
尚 恒川
英雄 沼上
隆志 近藤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Tokyo Shibaura Electric Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tokyo Shibaura Electric Co Ltd filed Critical Tokyo Shibaura Electric Co Ltd
Priority to JP51121658A priority Critical patent/JPS594745B2/ja
Publication of JPS5346797A publication Critical patent/JPS5346797A/ja
Publication of JPS594745B2 publication Critical patent/JPS594745B2/ja
Expired legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Optical investigation techniques, e.g. flow cytometry
    • G01N15/1468Optical investigation techniques, e.g. flow cytometry with spatial resolution of the texture or inner structure of the particle

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

【発明の詳細な説明】 本発明は細胞診における例えばガン細胞等の疑惑性の細
胞を選択抽出して細胞診を能率良く行うことのできる細
胞診装置に関する。
近年、予防医学上、ガンの集団検診が行われるようにな
つてきた。
この集団検診は被検者1人当り約10万個の細胞を検査
してその細胞中から例えばガン細胞等の疑惑性の細胞を
抽出するものである。そして上記検診に於て疑惑性の細
胞が抽出された被検者に対してガン検査を行うものであ
る。し力走、上記検診、即ち細胞診には検査技師の大き
な負担と共に多くの検査時間を要し、また診断の信頼性
も十分なものではなかつた。そこで本発明者らは先に第
1図に示す如き系統の細胞診装置を提唱した。即ち、採
取された細胞をスライドガラス上に染色、塗抹し、この
塗抹された細胞を例えば撮像機等からなる細胞画像入力
部1で撮像して細胞画像信号を得る。そしてこの細胞画
像信号を入力する細胞画像検出部2は例えばマスク等を
用いて細胞らしき画像を検出する。この検出された細胞
らしき画像信号について特徴抽出部3は細胞の特徴とな
る核面積や、核内のクロマチン構造あるいは細胞核面積
対細胞質面積の比、即ちN/C比を抽出する。そしてこ
の抽出された特徴に基づいて診断部4は疑惑性の細胞で
あるか否かを診断している。しかしこのような従来の細
胞診装置では検出される細胞数が多い為に個々の細胞に
ついて十分に信頼性の高い診断を行うことができないと
いう問題があつた。また特徴抽出部3において、N/C
比を算出する為の細胞質濃度や細胞核濃度等の特徴パラ
メータの閾値を適確に設定することが難しかつた。即ち
、細胞核および細胞質をそれぞれ正確に抽出するべく濃
度閾値を適正に設定することが困難であつた為、こ0結
果前記N/Cを求めてもその信頼性が低く、結局疑悪性
の細胞だけを抽出することが甚だ困難であつた。この為
に診断の信頼性を十分高めることができなかつた。また
上記装置には大きなハードウェアを必要とし、現実的な
装置ではなかつた。本発明はこのような事情を考慮して
なされたもので、その目的とするところは、採取された
多くの細胞中から疑悪性の細胞を抽出して細胞診を行う
ようにし、能率良くしかも信頼性高く細胞診を行うこと
のできる細胞診装置を提供することにある。
本発明の概要は、採取された細胞の画像信号を得、この
画像信号から上記細胞の濃度ヒストグラムを作成する。
そしてこの濃度ヒストグラムから背景と細胞質との変化
点を見出し、この変化点によつて区分される前記濃度ヒ
ストグラムの背景部分を除く部分の濃度分布偏度、つま
り濃度分布の偏り度合を求め、この濃度分布偏度と所定
の閾値とを比較することによつて、前記細胞が疑悪性で
あるか否かを判定するようにしたものである。換言すれ
ば、上記濃度分布偏度がN/Cと対応関係にあることを
利用して、N/Cを求めることなしに上記濃度分布偏度
から細胞の疑悪性判定を行うようにしたものである。そ
して疑悪性と判定された細胞についてだけ細胞の特徴抽
出を行い、細胞診の能率化をはかると共に信頼性を高め
たものである。以下、本発明の一実施例を示す図面を参
照して詳細を説明する。
第2図は装置の概略系統図であり、前記第1図に示す系
統図の細胞画像検出部2と特徴抽出部3との間に細胞選
択部5を挿入したものである。即ち、採取された細胞は
細胞画像入力部1で細胞画像信号に変換されて細胞画像
検出部2へ供給される。そしてこの検出部2で血球や粘
液等の画像信号と分離され1つの細胞を示す細胞画像信
号だけが細胞選択部5へ転送される。この細胞選択部5
は前記検出部2から転送されてきた細胞画像信号につい
て、細胞濃度ヒストグラムを作成する。そしてこのヒス
トグラムから濃度分布偏度を作成して診断を要する細胞
か否かを判定する。そして診断を要すると判定された細
胞の画像信号は特徴抽出部3へ転送され特徴抽出される
。この特徴抽出された結果に基づいて診断部4は診断を
行い疑悪性細胞の抽出を行うようになつている。ところ
で前記細胞選択部5は第3図に示す如く構成になつてい
る。
前記細胞画像検出部2から転送された細胞画像信号は例
えばランダムアクセスICメモリ等の画像メモリ11に
一旦書き込まれる。この画像メモl月1に書き込まれた
画像信号は一画素づつ順次読み出されて濃度ヒストグラ
ム計算回路12に入力される。このヒストグラム計算回
路12は例えば複数のカウンタ等で構成され、上記読み
出された画素の濃度に対応したカウンタを逐次歩進する
。従つて前記画像メモリ11から全ての画素情報が読み
出されたとき、前記ヒストグラム計算回路12の各カウ
ンタにはそれぞれ分類された濃度毎に画素情報の存在個
数が累積計数されている。これらの各累積計数された値
の例を第4図の濃度ヒストグラムに示す。そしてヒスト
グラム計算回路12で計算された各カウンタの計数値は
順次読み出されて細胞質レベル計算回路13に供給され
る。このレベル計算回路13は例えば比較器等によつて
構成されるもので上記計算されたヒストグラムから例え
ば細胞画像信号の背景と細胞質との境界レベルや細胞質
と細胞核との境界レベル等を抽出する。そしてこれらの
レベルから細胞質濃度閾値を算出する。この算出された
濃度閾値に基づいて例えば乗算器、加算器、除算器等か
らなる分布偏度計算回路14は前記ヒストグラム計算回
路12から各濃度の計数値を読み出して分布偏度を計算
する。この分布偏度計算回路14で計算された分布偏度
は比較器15に出力される。この比較器15にはあらか
じめ設定された分布偏度判定用の閾値θMが入力されて
おり、この閾値θMと前記記算された分布偏度とを比較
している。そして前記分布偏度計算回路14で計算され
た分布偏度が上記閾値θMを超えたとき、前記画像メモ
リ11に書き込まれた細胞画像信号は診断を必要とする
と判定する。この判定結果は細胞選択信号として前記特
徴抽出部3へ送出されるようになつている。また上記細
胞選択信号と共に前記画像メモリ11に書き込まれた細
胞画像信号も特徴抽出部3へ転送される。また第3図中
、16は計算回路制御装置で、前記画像メモリ11、各
計算回路12,13,14及び比較器15等の各動作タ
イミングを制御するものである。そしてこの計算回路制
御装置16は本細胞選択部5における細胞選択が終了し
たことを示す終了信号を前記特徴抽出部3へ送出するよ
うになつている。以下、第4図に示す細胞の濃度ヒスト
グラムを例に本装置の細胞選択の原理を説明する。一般
に細胞画像の濃度ヒストグラムを作成すると3つのピー
クができる。これらの各ピークのうち濃度レベルの低い
方の第1のピークは、細胞の存在しない領域、即ち背景
に対応したピークであり、濃度レベルの中央部に位置す
る第2のピークは細胞質に対応したピークである。そし
て濃度レベルの高い方の第3のピークは細胞核に対応し
たピークである。このような濃度ヒストグラムから第1
及び第2のピーク間に存在する谷を調べる。この谷に位
置する濃度レベルDcは細胞質濃度の閾値を示し、濃度
レベルが背景から細胞質へ変つたことを示している。こ
の谷の点によつて、ヒストグラム中の背景部分と細胞質
部分とを容易に区別することが可能となる。また濃度レ
ベルDeは細胞画像信号のうちの最大濃度レベルを示し
、細胞核の一番濃度が高い部分に相当したものである。
これらの濃度レベルDc,deを検出することによつて
、細胞の濃度を示す領域は濃度レベルDcから濃度レベ
ルDeに相当することを見出すことができる。そこで上
記濃度レベルDc,deを基準として、細胞の濃度分布
偏度Mを次のようにして算出する。次式は、分布偏度算
出の為の計算式で、式中Diは濃度レベルDcからDe
間の濃度レベルでiはcからeまで変化する変数である
。但し、Niは各濃度レベルの画素数である。
この式において、分母は算出された細胞を示す濃度が全
て一様であり、細胞全てが核であると判定される場合の
基準値である。そして分子は実際に算出された各濃度の
画素数から求められる例えばモーメント的な値である。
従つて細胞が全て細胞核に相当する場合には分子で計算
される値は分母と等しくなり、分布偏度MはIと算出さ
れることになる。また細胞が全て細胞質に相当する場合
には濃度レベルDe,dc間に存在する画素の数は0(
零)となり分子は0(零)となる。更に上式は単調増加
関数であり、通常細胞診で用いられる細胞核面積と細胞
質面積との比、所謂N/C比によく近似したものである
。なお第5図は上式によつて求められた分布偏度M.l
!:.N/C比との関係を示したもので、計算機により
擬似細胞画像を作り、この擬似細胞画像から求めたもの
である。従つて算出された分布偏度MはN/C比に対応
しており、細胞の特徴的形態を示す一つのパラメータと
見なすことができる。そして上記算出された分布偏度ノ
Mを細胞の異常性を示す目安としての閾値θMと比較し
、その比較結果に基づいて細胞の選択を行なう。
第6図は前記第3図に示した細胞選択部5において上記
計算式に基づいて分布偏度Mを算出して細胞の選択を行
うようにした回路例である。
この図に示されるように濃度ヒストグラム計算回路12
はデコーダ21とカウンタレジスタ22によつて構成さ
れている。また細胞レベル計算回路13は比較器23、
レジスタ24、フリツプフロツプ(FF)25,26等
によつて構成されている。また分布偏度計算回路14は
乗算器28、加算器29,30及び除算器31,32等
によつて構成されている。また図中33,34はフリツ
プフロツプであり、35,36はアドレスカウンタであ
り、37はデコーダである。このような構成の装置によ
れば画像メモリ11に細胞画像信号が書き込まれるとス
タート信号が供給される。
このスタート信号によつてアドレスカウンタ35、カウ
ンタレジスタ22、FF25,26等をりセツトして初
期状態に設定し、FF33をセツトする。このFF33
のセツトによりゲート回路38,39が開かれて、クロ
ツク1はアドレスカウンタ35に入力され、クロツク2
はカウンタレジスタ22へ供給される。そしてアドレス
カウンタ38による画像メモリ11の番地指定が行われ
、上記メモリ11に書き込まれた細胞画像信号が順次読
み出される。この読み出された画像信号はデコーダ21
により濃度レベルを示す信号に変換されてカウンタレジ
スタ22に供給される。このカウンタレジスタ22は濃
度レベルに対応した複数のカウンタ22a,22b,・
・・?22nから構成されたものである。そして前述の
如くカウンタレジスタ22は濃度に応じたカウンタ22
a,22b,・・・,22nを歩進して上記細胞画像の
濃度ヒストグラムを作成する。この濃度ヒストグラムの
作成が終了したとき前記アドレスカウンタ35からキヤ
リ一信号が出力されて前記FF33をりセツトする。こ
のりセツトによつて前記ゲート回路38,39が閉成さ
れ、ゲート回路40,41が開成される。このゲート回
路40,41の開成によつて前記カウンタレジスタ22
で作成された濃度ヒストグラムの各画素数はアドレスカ
ウンタ36により濃度レベルの低い方から順次読み出さ
れ比較器23と共にレジスタ24aに入力される。そし
てこのレジスタ24aに収納されたデータは次のクロツ
ク信号によつてレジスタ24bに転送される。また前記
比較器23は先にレジスタ24bに収納されていた画素
数のデータとカウンタレジスタ22から入力したデータ
とを逐次比較している。そしてカウンタ24bのデータ
よりカウンタレジスタ22から読み出されたデータの値
が小さくなつたとき、即ち前述した第1のピークを検出
したときに信号を出力してFF25をセツトする。この
FF25のセツト状態において前記比較器23が先のデ
ータよりレジスタ22から読み出されたデータの値が大
きくなつたとき、即ち前述した谷を検出したとき、再度
信号を出力して前記FF25から出力される信号との論
理積信号によつてFF26をセツトする。このFF25
のセツトにより前記濃度レベルDcが設定される。そし
て前記FF25のセツトによつてゲート回路42が開成
されカウンタ27が歩進される。つょり、カウンタ27
は濃度レベルDcから歩進され、その計数値は上式にお
ける(Di−Dc)を示している。そしてこのカウンタ
27から出力される信号と前記カウンタレジスタ22か
ら出力される信号とを乗算器28は乗算し、Niホ(D
i−De)を算出している。この乗算器28から出力さ
れる信号は加C算器29によつて累積加算され、.ΣN
iネ(Di−De)1:eを算出している。
そしてこの加算器29から算出される信号は除算器31
において前記カウンタ27からの信号によつて除算され
ている。従つて除算”C器31(力股後?(ΣNiネ(
Di−De)/(Dc−De)を算出丁:eする。
一方、加算器30は前記カウンタレジスタC22から出
力される信号を累積加算して.?Niを算出している。
そして除算器32は前記除算器31から出力された信号
を加算器30から出力された信号で除算してを算出する
しかして除算器32からは前述した分布偏度Mを示す信
号が出力されることになる。そしてこの除算器32から
出力される信号は比較器15に出力され、閾値θMと比
較されるようになつている。このように本装置によれば
従来の細胞診の特徴パラメータの一つであるN/C比を
用いることなく細胞を選択することができる。そしてN
/C比はその基準レベルの設定に多くの不碓定要素を含
み信頼性の高い特徴抽出を行うことができなかつたこと
に比して非常に簡単な原理によつて信頼性の高い特徴抽
出を行うことが可能となつた。従つて従来の装置に比し
て細胞診の信頼度を高くすることができると共に上述し
た如く簡単な回路構成において特徴抽出ができ、装置の
簡略化をはかることができる。また本装置は極めて高速
度に分布偏度Mの算出を行うことができ、更に細胞画像
信号に混在する雑音や画像信号の所謂ぼけ等の悪影響を
受けにくい。そして従来のN/C比算出手段のように濃
度閾値の変動による大きな誤差を生じることもない。更
に細胞選択部5において細胞診を要する細胞を選択して
、この選択された細胞についてだけ細胞診を行う為に極
めて能率良く細胞診を行うことができる。なお本発明は
上記実施例に限定されるものではない。
例えば濃度分布偏度の計算として、より明確な分布差を
もたせるのに濃度レベルDc−Deの重心をDmとして
として算出するようにしてもよい。
また簡単な細胞診手段として上記分布偏度M,質から直
接細胞の疑悪性を判断することもできる。そして上記濃
度分布偏度を算出する回路構成も用途、仕様等に応じ適
宜変形して用いてもよい。また濃度ヒストグラム作成に
おいて細胞の背景に汚れが多量に含まれる場合には、例
えばあらかじめ細胞の孤立性を検出し、この検出によつ
て背景に相当する領域の濃度レベル信号を除去する。し
かるのちに濃度ヒストグラムを作成して、上記汚れ等に
よる誤差を除去するようにしてもよい。また画像メモリ
11等も用途に応じて他の記憶素子他えばCCD等を用
い、濃度レベルをアナログ的に記憶しておいてもよい。
以上詳述したように本発明に係る細胞診装置は採取され
た細胞の細胞画像信号を得る手段と、この手段で得られ
た細胞画像信号から上記細胞の濃度ヒストグラムを作成
する手段と、この手段で作成された濃度ヒストグラムの
濃度分布偏度を算出する分布偏度算出装置と、この算出
装置で算出された濃度分布偏度に基づいて診断に供する
細胞を選択抽出する手段とを具備するものである。
従つて本発明によれば細胞の形態的特徴を信頼性高く、
しかも高速度で検出することができ、細胞診の能率を高
めることができると共に信頼性の向上をはかることがで
きる等の種々格別の利点を有する細胞診装置を提供する
ことができる。
【図面の簡単な説明】
第1図は従来の細胞診装置の一例を示す概略系統図、第
2図は本発明の一実施例を示す概略系統図、第3図は同
実施例における細胞選択部の概略構成図、第4図及び第
5図はそれぞれ本発明を説明する為の図、第6図は上記
細胞選択部の構成図である。 11・・・画像メモリ、12・・・濃度ヒストグラム計
算回路、13・・・細胞レベル計算回路、14・・・分
布偏度計算回路、15・・・比較器。

Claims (1)

    【特許請求の範囲】
  1. 1 採取された細胞から細胞画像信号を得る手段と、こ
    の細胞画像信号の濃度ヒストグラムを作成する手段と、
    この濃度ヒストグラムから背景と細胞質との変化点を判
    定し、この変化点で区別される前記濃度ヒストグラムの
    背景部分を除く部分の濃度分布偏度を求める手段と、こ
    の濃度分布偏度と所定の閾値とを比較して前記細胞が疑
    悪性か否かを判定する手段とを具備したことを特徴とす
    る細胞診装置。
JP51121658A 1976-10-09 1976-10-09 細胞診装置 Expired JPS594745B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP51121658A JPS594745B2 (ja) 1976-10-09 1976-10-09 細胞診装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP51121658A JPS594745B2 (ja) 1976-10-09 1976-10-09 細胞診装置

Publications (2)

Publication Number Publication Date
JPS5346797A JPS5346797A (en) 1978-04-26
JPS594745B2 true JPS594745B2 (ja) 1984-01-31

Family

ID=14816689

Family Applications (1)

Application Number Title Priority Date Filing Date
JP51121658A Expired JPS594745B2 (ja) 1976-10-09 1976-10-09 細胞診装置

Country Status (1)

Country Link
JP (1) JPS594745B2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020166469A1 (ja) * 2019-02-15 2020-08-20 コニカミノルタ株式会社 情報提供方法、情報提供装置及びプログラム

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5856588A (ja) * 1981-09-29 1983-04-04 Toa Medical Electronics Co Ltd 画像分析分法およびその装置
JPWO2022059508A1 (ja) * 2020-09-18 2022-03-24

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020166469A1 (ja) * 2019-02-15 2020-08-20 コニカミノルタ株式会社 情報提供方法、情報提供装置及びプログラム

Also Published As

Publication number Publication date
JPS5346797A (en) 1978-04-26

Similar Documents

Publication Publication Date Title
US4523278A (en) Method of automatic detection of cells and determination of cell features from cytological smear preparations
WO2021003821A1 (zh) 一种肾小球病理切片图像的细胞检测方法、装置及设备
US8737714B2 (en) System and method for assessing image interpretability in anatomic pathology
CN110288599B (zh) 一种坏点检测方法、装置、电子设备及存储介质
CN109086734A (zh) 一种对人眼图像中瞳孔图像进行定位的方法及装置
CN109510946B (zh) Hdr场景检测方法及系统
WO2023130661A1 (zh) 一种二维频谱多普勒超声心动图像的处理方法和装置
CN106529545A (zh) 一种基于图像特征描述的散斑图像质量识别方法及系统
CN104396231A (zh) 图像处理装置和图像处理方法
US7760939B2 (en) Method for automatically generating a background mask for images with noisy background areas, applications to that end and a nuclear magnetic resonance tomography apparatus for carrying out the method
CN111192241B (zh) 一种人脸图像的质量评估方法、装置及计算机存储介质
CN111354001A (zh) 脑肿瘤图像分割方法、装置和存储介质
US4254779A (en) Brain wave analysis apparatus for use in diagnosing living body conditions
KR20190067477A (ko) 아밀로이드 페트(pet) 영상을 이용한 알츠하이머 치매의 병기 측정 방법 및 장치
Hadian-Jazi et al. Data reduction for serial crystallography using a robust peak finder
CN108426521A (zh) 一种构件的质量检测方法及装置
JPS594745B2 (ja) 細胞診装置
CN113288214A (zh) 超声多普勒频谱数据的处理方法、设备及可读存储介质
WO2023202662A1 (zh) 基于排列随机森林的cest磁共振分析方法及装置
CN106709516B (zh) 一种基于朴素贝叶斯方法的模糊图像检测方法
US11928184B2 (en) Generation of a second object model based on a first object model for use in object matching
KR102254970B1 (ko) 전립선암 강조영상 생성 방법 및 장치
CN110991370B (zh) 多通道信息融合atm面板遗留物检测方法
JPH0225134B2 (ja)
JPS5830049B2 (ja) 網赤血球自動計測装置