JP2002312761A - 細胞画像の画像処理方法 - Google Patents
細胞画像の画像処理方法Info
- Publication number
- JP2002312761A JP2002312761A JP2001113692A JP2001113692A JP2002312761A JP 2002312761 A JP2002312761 A JP 2002312761A JP 2001113692 A JP2001113692 A JP 2001113692A JP 2001113692 A JP2001113692 A JP 2001113692A JP 2002312761 A JP2002312761 A JP 2002312761A
- Authority
- JP
- Japan
- Prior art keywords
- image
- luminance
- cell
- pixel
- binarization
- 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
Landscapes
- Investigating Or Analysing Biological Materials (AREA)
- Apparatus Associated With Microorganisms And Enzymes (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】
【課題】 細胞画像から正しく形状抽出を行うことがで
きる細胞画像の画像処理方法を提供することを目的とす
る。 【解決手段】 透過照明にて得られた細胞画像を画像処
理して細胞本体２４ａから突出する突起２４ｂを検出す
る突起抽出処理において、細胞画像から形状抽出用の探
査用画像を取得するために行われる動的二値化処理に先
だって、原画像中の輝度分布から透過照明光の回折によ
って細胞本体２４ａの周囲に生じた高輝度部分を検出
し、検出した高輝度部分の輝度値を平均値Ｍに置換して
これらの高輝度部分を平均輝度部２３’として消去する
処理を行う。これにより動的二値化処理において高輝度
部分の影響による二値化処理の誤差を排除して、正しい
形状抽出を行うことができる。
きる細胞画像の画像処理方法を提供することを目的とす
る。 【解決手段】 透過照明にて得られた細胞画像を画像処
理して細胞本体２４ａから突出する突起２４ｂを検出す
る突起抽出処理において、細胞画像から形状抽出用の探
査用画像を取得するために行われる動的二値化処理に先
だって、原画像中の輝度分布から透過照明光の回折によ
って細胞本体２４ａの周囲に生じた高輝度部分を検出
し、検出した高輝度部分の輝度値を平均値Ｍに置換して
これらの高輝度部分を平均輝度部２３’として消去する
処理を行う。これにより動的二値化処理において高輝度
部分の影響による二値化処理の誤差を排除して、正しい
形状抽出を行うことができる。
Description
【０００１】
【発明の属する技術分野】本発明は、動植物の細胞画像
を画像処理することにより所定のデータ抽出を行う細胞
画像の画像処理方法に関するものである。
を画像処理することにより所定のデータ抽出を行う細胞
画像の画像処理方法に関するものである。
【０００２】
【従来の技術】製薬業界等で行われる創薬スクリーニン
グ作業においては、動植物の細胞や微生物などをスクリ
ーニング対象の薬剤が種々の条件で投与された状態で培
養し、経時的変化を観察することにより薬剤の効果を確
認する作業が行われる。例えば細胞の神経突起など逐次
伸長する部分を観察対象とする場合には、これらを経時
的に追跡して観察・記録し、データ化する作業が必要と
なる。
グ作業においては、動植物の細胞や微生物などをスクリ
ーニング対象の薬剤が種々の条件で投与された状態で培
養し、経時的変化を観察することにより薬剤の効果を確
認する作業が行われる。例えば細胞の神経突起など逐次
伸長する部分を観察対象とする場合には、これらを経時
的に追跡して観察・記録し、データ化する作業が必要と
なる。
【０００３】この観察作業において、従来一般的に行わ
れていた作業者の目視観察に依存する方法に替えて、近
年この観察作業やデータ化処理を画像処理によって自動
的に行う試みが進行している。この方法は、対象となる
細胞の顕微画像を取得し、この画像から必要部分のみを
抽出する処理を行い、抽出された画像から所定のデータ
を取得するものである。この必要部分の抽出において
は、取得画像中の画素を輝度によって白黒に二分する二
値化処理が行われ、これによりデータ取得のための形状
抽出対象となる探査用画像が作成される。
れていた作業者の目視観察に依存する方法に替えて、近
年この観察作業やデータ化処理を画像処理によって自動
的に行う試みが進行している。この方法は、対象となる
細胞の顕微画像を取得し、この画像から必要部分のみを
抽出する処理を行い、抽出された画像から所定のデータ
を取得するものである。この必要部分の抽出において
は、取得画像中の画素を輝度によって白黒に二分する二
値化処理が行われ、これによりデータ取得のための形状
抽出対象となる探査用画像が作成される。
【０００４】
【発明が解決しようとする課題】ところが従来の細胞画
像処理において用いられる二値化処理には、以下のよう
な不具合があった。画像処理対象の顕微画像取得時に
は、撮像対象を透過照明によって照明した状態で、顕微
鏡による撮像が行われる。すなわち撮像対象の細胞の背
後から投射された照明光がカメラに入射することによ
り、細胞が撮像される。
像処理において用いられる二値化処理には、以下のよう
な不具合があった。画像処理対象の顕微画像取得時に
は、撮像対象を透過照明によって照明した状態で、顕微
鏡による撮像が行われる。すなわち撮像対象の細胞の背
後から投射された照明光がカメラに入射することによ
り、細胞が撮像される。
【０００５】このとき、細胞本体の周囲を通過した光が
カメラに入射する際の光の回折により、取得された画像
上で細胞本体の周囲に実際には存在しない高輝度部分が
輪状に発生する場合がある。そして画像にこのような高
輝度部分が存在すると、前述の二値化処理に際して誤っ
た形状抽出が行われる。
カメラに入射する際の光の回折により、取得された画像
上で細胞本体の周囲に実際には存在しない高輝度部分が
輪状に発生する場合がある。そして画像にこのような高
輝度部分が存在すると、前述の二値化処理に際して誤っ
た形状抽出が行われる。
【０００６】すなわち、二値化処理においては二値化処
理対象範囲毎にこの範囲内の輝度の平均値に基づいてし
きい値を設定し、このしきい値を用いて範囲内画素の白
黒の切り分けが行われる。このとき、対象範囲内に本来
存在しない高輝度部分が存在する場合には画素の平均輝
度が上昇することから、設定されたしきい値は高輝度側
に偏ったものとなる。
理対象範囲毎にこの範囲内の輝度の平均値に基づいてし
きい値を設定し、このしきい値を用いて範囲内画素の白
黒の切り分けが行われる。このとき、対象範囲内に本来
存在しない高輝度部分が存在する場合には画素の平均輝
度が上昇することから、設定されたしきい値は高輝度側
に偏ったものとなる。
【０００７】そしてこのしきい値を用いて二値化を行う
と、本来は黒とされるべきでない画素が黒として設定さ
れる結果、二値化処理によって得られた探査画像から抽
出される形状は本来の正しい形状を示さず、誤った形状
抽出結果を与える。このように従来の細胞画像の画像処
理においては、撮像時の回折光の影響により形状抽出結
果に誤差を生じるという問題点があった。
と、本来は黒とされるべきでない画素が黒として設定さ
れる結果、二値化処理によって得られた探査画像から抽
出される形状は本来の正しい形状を示さず、誤った形状
抽出結果を与える。このように従来の細胞画像の画像処
理においては、撮像時の回折光の影響により形状抽出結
果に誤差を生じるという問題点があった。
【０００８】そこで本発明は、細胞画像から正しく形状
抽出を行うことができる細胞画像の画像処理方法を提供
することを目的とする。
抽出を行うことができる細胞画像の画像処理方法を提供
することを目的とする。
【０００９】
【課題を解決するための手段】請求項１記載の細胞画像
の画像処理方法は、透過照明にて得られた細胞画像にお
いて、光の回折によって細胞本体の周囲に生じた高輝度
部分を検出する工程と、検出した前記高輝度部分の輝度
の値を所定の値に置換してこれらの高輝度部分を消去す
る工程と、前記高輝度部分が消去された細胞画像に対し
て動的二値化処理を行う工程とを含む。
の画像処理方法は、透過照明にて得られた細胞画像にお
いて、光の回折によって細胞本体の周囲に生じた高輝度
部分を検出する工程と、検出した前記高輝度部分の輝度
の値を所定の値に置換してこれらの高輝度部分を消去す
る工程と、前記高輝度部分が消去された細胞画像に対し
て動的二値化処理を行う工程とを含む。
【００１０】請求項２記載の細胞画像の画像処理方法
は、前記高輝度部分の検出は、少なくとも動的二値化処
理の対象となる複数の画素の輝度の分布より高輝度検出
用のしきい値を設定し、輝度がこのしきい値を超えた画
素を高輝度部分とする。
は、前記高輝度部分の検出は、少なくとも動的二値化処
理の対象となる複数の画素の輝度の分布より高輝度検出
用のしきい値を設定し、輝度がこのしきい値を超えた画
素を高輝度部分とする。
【００１１】請求項３記載の細胞画像の画像処理方法
は、前記高輝度部分の輝度の値と置換する所定の値は、
少なくとも動的二値化処理の対象となる複数の画素の輝
度の値の平均値である。
は、前記高輝度部分の輝度の値と置換する所定の値は、
少なくとも動的二値化処理の対象となる複数の画素の輝
度の値の平均値である。
【００１２】請求項４記載の細胞画像の画像処理方法
は、前記動的二値化処理は、二値化処理の対象となる一
つの画素とその周囲の複数の画素の輝度の分布から二値
化しきい値を設定する処理と、前記二値化処理の対象と
なる一つの画素の輝度の値と前記二値化しきい値とを比
較する処理とを含む。
は、前記動的二値化処理は、二値化処理の対象となる一
つの画素とその周囲の複数の画素の輝度の分布から二値
化しきい値を設定する処理と、前記二値化処理の対象と
なる一つの画素の輝度の値と前記二値化しきい値とを比
較する処理とを含む。
【００１３】請求項５記載の細胞画像の画像処理方法
は、前記二値化しきい値が、二値化処理の対象となる一
つの画素とその周囲の複数の画素の輝度の値の平均値で
ある。
は、前記二値化しきい値が、二値化処理の対象となる一
つの画素とその周囲の複数の画素の輝度の値の平均値で
ある。
【００１４】上記構成の発明によれば、細胞画像から形
状抽出用の探査用画像を取得するために行われる二値化
処理に先だって、透過照明光の回折によって細胞本体の
周囲に生じた高輝度部分を検出し、検出した高輝度部分
の輝度の値を所定の値に置換してこれらの高輝度部分を
消去する処理を行うことにより、高輝度部分の影響を排
除して正しい形状抽出を行うことができる。
状抽出用の探査用画像を取得するために行われる二値化
処理に先だって、透過照明光の回折によって細胞本体の
周囲に生じた高輝度部分を検出し、検出した高輝度部分
の輝度の値を所定の値に置換してこれらの高輝度部分を
消去する処理を行うことにより、高輝度部分の影響を排
除して正しい形状抽出を行うことができる。
【００１５】
【発明の実施の形態】次に本発明の一実施の形態を図面
を参照して説明する。図１は本発明の一実施の形態の画
像処理装置の構成を示すブロック図、図２は本発明の一
実施の形態の細胞画像の画像処理方法のメインフローを
示すフロー図、図３は本発明の一実施の形態の細胞画像
の画像処理方法における二値化処理のフロー図、図４は
本発明の一実施の形態の細胞画像の画像処理方法におけ
るサーチライン設定処理のフロー図、図５は本発明の一
実施の形態の細胞画像の画像処理方法における突起検出
処理のフロー図、図６は本発明の一実施の形態の細胞画
像の画像処理方法における突起判定処理のフロー図、図
７は本発明の一実施の形態の細胞画像を示す画像図、図
８、図９は本発明の一実施の形態の細胞画像の画像処理
方法における二値化処理の説明図、図１０は本発明の一
実施の形態の細胞画像の画像処理方法におけるサーチラ
イン設定処理の説明図、図１１、図１２、図１３は本発
明の一実施の形態の細胞画像の画像処理方法における突
起検出処理の説明図、図１４は本発明の一実施の形態の
細胞画像の画像処理方法における突起判定処理の説明図
である。
を参照して説明する。図１は本発明の一実施の形態の画
像処理装置の構成を示すブロック図、図２は本発明の一
実施の形態の細胞画像の画像処理方法のメインフローを
示すフロー図、図３は本発明の一実施の形態の細胞画像
の画像処理方法における二値化処理のフロー図、図４は
本発明の一実施の形態の細胞画像の画像処理方法におけ
るサーチライン設定処理のフロー図、図５は本発明の一
実施の形態の細胞画像の画像処理方法における突起検出
処理のフロー図、図６は本発明の一実施の形態の細胞画
像の画像処理方法における突起判定処理のフロー図、図
７は本発明の一実施の形態の細胞画像を示す画像図、図
８、図９は本発明の一実施の形態の細胞画像の画像処理
方法における二値化処理の説明図、図１０は本発明の一
実施の形態の細胞画像の画像処理方法におけるサーチラ
イン設定処理の説明図、図１１、図１２、図１３は本発
明の一実施の形態の細胞画像の画像処理方法における突
起検出処理の説明図、図１４は本発明の一実施の形態の
細胞画像の画像処理方法における突起判定処理の説明図
である。
【００１６】まず図１を参照して細胞画像の画像処理装
置の構成を説明する。図１において位置決めステージ１
上には、観察対象の細胞を収容するマイクロタイタープ
レート２が載置されている。マイクロタイタープレート
２は多数のウェル２ａを備えており、ウェル２ａ内に
は、観察対象の動植物の細胞を含む試料３が収容されて
いる。位置決めステージ１の下方には顕微観察を行う光
学系４が配設されており、光学系４の下方にはカメラ５
が配置されている。位置決めステージ１の上方に配設さ
れた照明装置６を点灯してウェル２ａを照明した状態
で、光学系４を介してカメラ５によってウェル２ａ内の
試料３中の動植物の細胞を撮像することにより、カメラ
５は細胞画像（図７参照）を取得する。
置の構成を説明する。図１において位置決めステージ１
上には、観察対象の細胞を収容するマイクロタイタープ
レート２が載置されている。マイクロタイタープレート
２は多数のウェル２ａを備えており、ウェル２ａ内に
は、観察対象の動植物の細胞を含む試料３が収容されて
いる。位置決めステージ１の下方には顕微観察を行う光
学系４が配設されており、光学系４の下方にはカメラ５
が配置されている。位置決めステージ１の上方に配設さ
れた照明装置６を点灯してウェル２ａを照明した状態
で、光学系４を介してカメラ５によってウェル２ａ内の
試料３中の動植物の細胞を撮像することにより、カメラ
５は細胞画像（図７参照）を取得する。
【００１７】カメラ５は画像処理部１０と接続されてお
り、画像処理部１０は以下に説明する各部より構成され
る。画像記憶部１１は、カメラ５によって取得された細
胞画像（原画像）データを記憶する。二値化処理部１２
は、突起抽出データ取得のための形状抽出対象となる探
査用画像を原画像より生成する。二値化画像記憶部１３
は、生成された二値化画像を記憶する。
り、画像処理部１０は以下に説明する各部より構成され
る。画像記憶部１１は、カメラ５によって取得された細
胞画像（原画像）データを記憶する。二値化処理部１２
は、突起抽出データ取得のための形状抽出対象となる探
査用画像を原画像より生成する。二値化画像記憶部１３
は、生成された二値化画像を記憶する。
【００１８】細胞本体検出処理部１４は、原画像データ
を周波数変換することにより、細胞画像中の突起を消去
した細胞本体画像を検出する。サーチライン設定処理部
１５は、細胞本体画像に基づいて、細胞本体画像の周囲
に突起検出用のサーチラインを設定する。サーチライン
記憶部１６は、設定されたサーチラインの位置情報を記
憶する。突起情報記憶部１８は、突起検出処理において
サーチラインを横切るエッジとして検出された突出部分
についての突起情報を記憶する。突起判定処理部１９
は、記憶された突出部分が、検出目的に合致した突起に
該当するか否かを判定する。
を周波数変換することにより、細胞画像中の突起を消去
した細胞本体画像を検出する。サーチライン設定処理部
１５は、細胞本体画像に基づいて、細胞本体画像の周囲
に突起検出用のサーチラインを設定する。サーチライン
記憶部１６は、設定されたサーチラインの位置情報を記
憶する。突起情報記憶部１８は、突起検出処理において
サーチラインを横切るエッジとして検出された突出部分
についての突起情報を記憶する。突起判定処理部１９
は、記憶された突出部分が、検出目的に合致した突起に
該当するか否かを判定する。
【００１９】画像処理部１０を構成するハードウェアと
しては、パーソナルコンピュータのようにコンピュータ
プログラムによって作動する汎用的な処理装置が用いら
れる。本実施の形態ではこのような処理装置に、画像処
理部１０として機能させるためのコンピュータプログラ
ム及び二値化処理部１２、細胞本体検出処理部１４、サ
ーチライン設定処理部１５、突起検出処理部１７、突起
判定処理部１９の機能を実現させるためのプログラムを
インストールして細胞画像の処理装置を構成している。
しては、パーソナルコンピュータのようにコンピュータ
プログラムによって作動する汎用的な処理装置が用いら
れる。本実施の形態ではこのような処理装置に、画像処
理部１０として機能させるためのコンピュータプログラ
ム及び二値化処理部１２、細胞本体検出処理部１４、サ
ーチライン設定処理部１５、突起検出処理部１７、突起
判定処理部１９の機能を実現させるためのプログラムを
インストールして細胞画像の処理装置を構成している。
【００２０】次に細胞画像の画像処理方法について、各
図を参照して説明する。この画像処理方法は、神経細胞
などの細胞本体から外側に突出する突起を、画像処理手
法によって細胞画像から抽出するものである。ここで
は、観察対象の細胞を撮像して得られた細胞画像（図７
参照）中の細胞本体２４ａの周囲をサーチして突起２４
ｂの根本を検出し、この根本から先端へ向かって連続す
るエッジを検出する方法が用いられる。
図を参照して説明する。この画像処理方法は、神経細胞
などの細胞本体から外側に突出する突起を、画像処理手
法によって細胞画像から抽出するものである。ここで
は、観察対象の細胞を撮像して得られた細胞画像（図７
参照）中の細胞本体２４ａの周囲をサーチして突起２４
ｂの根本を検出し、この根本から先端へ向かって連続す
るエッジを検出する方法が用いられる。
【００２１】まず図２のフロー図に沿って画像処理方法
のメインフローについて説明する。図２のフロー図は、
画像処理部１０の機能を実現するための処理プログラム
である。ここでは、上記処理の概要について述べ、詳細
な説明を要するステップについては個別に後述する。図
２において、取得された細胞画像（原画像）の二値化処
理が行われる（ＳＴ１）。次に細胞本体２４ａの外周に
突起の根本を検出するためのサーチラインを設定するサ
ーチライン設定処理が行われる（ＳＴ２）。
のメインフローについて説明する。図２のフロー図は、
画像処理部１０の機能を実現するための処理プログラム
である。ここでは、上記処理の概要について述べ、詳細
な説明を要するステップについては個別に後述する。図
２において、取得された細胞画像（原画像）の二値化処
理が行われる（ＳＴ１）。次に細胞本体２４ａの外周に
突起の根本を検出するためのサーチラインを設定するサ
ーチライン設定処理が行われる（ＳＴ２）。
【００２２】そして二値化処理により得られた二値化画
像（探査用画像）においてサーチラインＳＬに該当する
位置（画素）に沿って突起抽出が行われる（図１１参
照）。すなわち、サーチラインに沿って輝度の変化を確
認することによりサーチライン上を横切る突起の根本の
エッジを検出し、このエッジを起点として突起の輪郭を
検出する。
像（探査用画像）においてサーチラインＳＬに該当する
位置（画素）に沿って突起抽出が行われる（図１１参
照）。すなわち、サーチラインに沿って輝度の変化を確
認することによりサーチライン上を横切る突起の根本の
エッジを検出し、このエッジを起点として突起の輪郭を
検出する。
【００２３】図１１に示すように探査用画像を構成する
格子配列の各画素Ａには、位置情報を示すラベルｉが付
されている。すなわち、ラベルｉを指定することによ
り、その画素を注目画素とし、またラベルｉの値で左上
端に位置する画素Ａ［０］から右下端に位置するＡ［Ｓ
ＩＺＥ−１］まで、各画素の位置を特定することができ
る。上述のサーチライン検出のための画像処理において
は、処理の対象となる注目画素Ａ［ｉ］を各画素位置に
移動させながら所定の処理ステップが反復実行される。
格子配列の各画素Ａには、位置情報を示すラベルｉが付
されている。すなわち、ラベルｉを指定することによ
り、その画素を注目画素とし、またラベルｉの値で左上
端に位置する画素Ａ［０］から右下端に位置するＡ［Ｓ
ＩＺＥ−１］まで、各画素の位置を特定することができ
る。上述のサーチライン検出のための画像処理において
は、処理の対象となる注目画素Ａ［ｉ］を各画素位置に
移動させながら所定の処理ステップが反復実行される。
【００２４】このため、（ＳＴ２）以後のメインフロー
においては、まず探査用画像の画像サイズ入力が行わ
れ、探査用画像の画素数（ＳＩＺＥ）を画像サイズとし
て入力し（ＳＴ３）、ラベルｉを０とする（ＳＴ４）こ
とにより、最初の画素Ａ［０］から順次各画素について
以下の処理ステップが行われる。
においては、まず探査用画像の画像サイズ入力が行わ
れ、探査用画像の画素数（ＳＩＺＥ）を画像サイズとし
て入力し（ＳＴ３）、ラベルｉを０とする（ＳＴ４）こ
とにより、最初の画素Ａ［０］から順次各画素について
以下の処理ステップが行われる。
【００２５】まず対象となる当該注目画素Ａ［ｉ］はサ
ーチラインに該当するか否かを判断し（ＳＴ５）、サー
チラインに該当すると判断されたならば、さらに当該サ
ーチラインは未だ以降の処理の対象となっていない未検
出のサーチラインであるか否かを判断し（ＳＴ６）、未
検出のサーチラインであれば、突起抽出を行う対象とな
るサーチラインが検出されたと判断して、突起検出処理
（ＳＴ７）を行う。
ーチラインに該当するか否かを判断し（ＳＴ５）、サー
チラインに該当すると判断されたならば、さらに当該サ
ーチラインは未だ以降の処理の対象となっていない未検
出のサーチラインであるか否かを判断し（ＳＴ６）、未
検出のサーチラインであれば、突起抽出を行う対象とな
るサーチラインが検出されたと判断して、突起検出処理
（ＳＴ７）を行う。
【００２６】この（ＳＴ５）〜（ＳＴ７）のサーチライ
ン検出処理（及び検出されたサーチラインについて行わ
れる突起検出処理）は、探査用画像の全ての画素が注目
画素とされ、所要処理ステップが完了するまで反復され
る。すなわち、ラベルｉをｉ＋１に置き換え（ＳＴ
８）、ここでｉが画像サイズとして入力された画素数
（ＳＩＺＥ）に到達したか否かを判断し、到達していな
ければ（ＳＴ５）に戻ってサーチライン検出処理を反復
する。
ン検出処理（及び検出されたサーチラインについて行わ
れる突起検出処理）は、探査用画像の全ての画素が注目
画素とされ、所要処理ステップが完了するまで反復され
る。すなわち、ラベルｉをｉ＋１に置き換え（ＳＴ
８）、ここでｉが画像サイズとして入力された画素数
（ＳＩＺＥ）に到達したか否かを判断し、到達していな
ければ（ＳＴ５）に戻ってサーチライン検出処理を反復
する。
【００２７】そして、ラベルｉが画素数（ＳＩＺＥ）に
到達している場合には、全画素を注目画素としたサーチ
ライン検出処理及び検出されたサーチラインについて行
われる突起検出処理が完了したと判断し、突起判定処理
（ＳＴ１０）に移行する。ここでは（ＳＴ７）にて検出
された突起が検出目的に該当する突起であるか否かを判
定する。そして突起判定の後、細胞画像の画像処理を終
了する。
到達している場合には、全画素を注目画素としたサーチ
ライン検出処理及び検出されたサーチラインについて行
われる突起検出処理が完了したと判断し、突起判定処理
（ＳＴ１０）に移行する。ここでは（ＳＴ７）にて検出
された突起が検出目的に該当する突起であるか否かを判
定する。そして突起判定の後、細胞画像の画像処理を終
了する。
【００２８】以上が細胞画像の画像処理方法のメインフ
ローであり、以下個別処理の詳細について説明する。ま
ず図８，図９を参照して、二値化処理について図３のフ
ロー図に沿って説明する。この二値化処理は、透過照明
によって細胞を撮像して得られた細胞画像（原画像）か
ら、形状抽出のための探査用画像を生成するものであ
り、後述するように同一対象画像内において、二値化対
象画素に応じてしきい値を変動させる動的二値化処理が
用いられる。
ローであり、以下個別処理の詳細について説明する。ま
ず図８，図９を参照して、二値化処理について図３のフ
ロー図に沿って説明する。この二値化処理は、透過照明
によって細胞を撮像して得られた細胞画像（原画像）か
ら、形状抽出のための探査用画像を生成するものであ
り、後述するように同一対象画像内において、二値化対
象画素に応じてしきい値を変動させる動的二値化処理が
用いられる。
【００２９】次に二値化処理の対象となる原画像につい
て説明する。図７に示すように原画像では、各種のノイ
ズ成分を含む背景部分２２の中に細胞２４の画像が現れ
ている。細胞２４は、細胞本体２４ａとこの細胞本体か
ら突出する突起２４ｂより成る。
て説明する。図７に示すように原画像では、各種のノイ
ズ成分を含む背景部分２２の中に細胞２４の画像が現れ
ている。細胞２４は、細胞本体２４ａとこの細胞本体か
ら突出する突起２４ｂより成る。
【００３０】この原画像においては、細胞本体２４ａの
周囲には輪状の高輝度部分２３が現出している。このよ
うな高輝度部分は、照明装置６からの透過照明光によっ
てカメラ５で細胞を撮像する場合に、照明光が細胞本体
２４ａの近傍を通過する際に回折することによって生じ
るものである。そしてこのような本来存在しない高輝度
部分を含む画像をそのまま動的二値化処理の対象とする
と、形状抽出用として適切な二値化画像を得ることがで
きない場合が生じるため、本実施の形態に示す二値化処
理においては、以下に説明する前処理を施した画像を動
的二値化処理の対象とするようにしている。
周囲には輪状の高輝度部分２３が現出している。このよ
うな高輝度部分は、照明装置６からの透過照明光によっ
てカメラ５で細胞を撮像する場合に、照明光が細胞本体
２４ａの近傍を通過する際に回折することによって生じ
るものである。そしてこのような本来存在しない高輝度
部分を含む画像をそのまま動的二値化処理の対象とする
と、形状抽出用として適切な二値化画像を得ることがで
きない場合が生じるため、本実施の形態に示す二値化処
理においては、以下に説明する前処理を施した画像を動
的二値化処理の対象とするようにしている。
【００３１】図３のフロー図は、二値化処理部１２の機
能を実現するための二値化プログラムを示している。ま
ず図８（ａ）に示すように二値化対象の全画素の輝度の
分布を求める（ＳＴ１１）。図８（ａ）は輝度の分布を
示すためのヒストグラムであり、横軸に輝度、縦軸には
度数（画素数）をとっている。次いで輝度の分布から輝
度の平均値Ｍを求めるとともに、高輝度検出用のしきい
値ＴＨを設定し（ＳＴ１２）、輝度値Ｂがこのしきい値
ＴＨを超える範囲ＨＢに対応する画素を、高輝度部分と
して抽出する（ＳＴ１３）。このしきい値ＴＨの設定方
法としては、しきい値を超える画素の全画素に対する割
合が所定の値を超えるような値に設定する方法等、様々
な手法が適用できる。透過照明にて得られた細胞画像に
おいて、光の回折によって細胞本体の周囲に生じた高輝
度部分を検出する。これにより、図７に示す原画像中の
高輝度部分２３が画素単位で特定される。
能を実現するための二値化プログラムを示している。ま
ず図８（ａ）に示すように二値化対象の全画素の輝度の
分布を求める（ＳＴ１１）。図８（ａ）は輝度の分布を
示すためのヒストグラムであり、横軸に輝度、縦軸には
度数（画素数）をとっている。次いで輝度の分布から輝
度の平均値Ｍを求めるとともに、高輝度検出用のしきい
値ＴＨを設定し（ＳＴ１２）、輝度値Ｂがこのしきい値
ＴＨを超える範囲ＨＢに対応する画素を、高輝度部分と
して抽出する（ＳＴ１３）。このしきい値ＴＨの設定方
法としては、しきい値を超える画素の全画素に対する割
合が所定の値を超えるような値に設定する方法等、様々
な手法が適用できる。透過照明にて得られた細胞画像に
おいて、光の回折によって細胞本体の周囲に生じた高輝
度部分を検出する。これにより、図７に示す原画像中の
高輝度部分２３が画素単位で特定される。
【００３２】そして特定された高輝度部分２３の輝度の
値を、（ＳＴ１２）で求めた平均値Ｍと置換する（ＳＴ
１４）。すなわち高輝度部分の輝度の値を所定の値に置
換してこれらの高輝度部分を消去する。これにより、図
８（ｂ）に示すように原画像のうち高輝度部分２３の輝
度を平均値Ｍで置き換えて平均輝度部２３’とした前処
理画像が生成される。この後、このようにして得られた
前処理画像を対象として、動的二値化処理によって全画
素を二値化処理し、得られた二値化画像を二値化画像記
憶部１３に記憶させる（ＳＴ１５）。
値を、（ＳＴ１２）で求めた平均値Ｍと置換する（ＳＴ
１４）。すなわち高輝度部分の輝度の値を所定の値に置
換してこれらの高輝度部分を消去する。これにより、図
８（ｂ）に示すように原画像のうち高輝度部分２３の輝
度を平均値Ｍで置き換えて平均輝度部２３’とした前処
理画像が生成される。この後、このようにして得られた
前処理画像を対象として、動的二値化処理によって全画
素を二値化処理し、得られた二値化画像を二値化画像記
憶部１３に記憶させる（ＳＴ１５）。
【００３３】ここで動的二値化処理について説明する。
この動的二値化処理においては、図８（ｃ）に示すよう
に、二値化対象画素Ａ［ｊ］の周囲に複数の画素Ａを含
む参照範囲Ｚを設定し、この参照範囲Ｚ内の画素Ａの輝
度の値に基づいて、二値化対象画素Ａ［ｊ］の二値化し
きい値を設定する。ここでは、参照範囲Ｚ内の画素Ａの
輝度の平均値を二値化しきい値として用いている。そし
て二値化対象画素Ａ［ｊ］の輝度を二値化しきい値と比
較することにより、二値化対象画素Ａ［ｊ］の輝度値
は、「１」または「０」のいずれかに置き換えられる。
この動的二値化処理においては、図８（ｃ）に示すよう
に、二値化対象画素Ａ［ｊ］の周囲に複数の画素Ａを含
む参照範囲Ｚを設定し、この参照範囲Ｚ内の画素Ａの輝
度の値に基づいて、二値化対象画素Ａ［ｊ］の二値化し
きい値を設定する。ここでは、参照範囲Ｚ内の画素Ａの
輝度の平均値を二値化しきい値として用いている。そし
て二値化対象画素Ａ［ｊ］の輝度を二値化しきい値と比
較することにより、二値化対象画素Ａ［ｊ］の輝度値
は、「１」または「０」のいずれかに置き換えられる。
【００３４】これにより、図９（ａ）に示すように、背
景画像を明部（輝度値「１」）とし、細胞本体２４ａ、
突起２４ｂを暗部（輝度値「０」）とする二値化画像を
得る。この原画像の二値化において前述のような動的二
値化処理を用いることにより、同一の画像内において領
域によって輝度値の分布に偏りがある場合にあっても、
画像内の全ての領域について正しい形状抽出結果を与え
る探査用画像を生成することができる。
景画像を明部（輝度値「１」）とし、細胞本体２４ａ、
突起２４ｂを暗部（輝度値「０」）とする二値化画像を
得る。この原画像の二値化において前述のような動的二
値化処理を用いることにより、同一の画像内において領
域によって輝度値の分布に偏りがある場合にあっても、
画像内の全ての領域について正しい形状抽出結果を与え
る探査用画像を生成することができる。
【００３５】ここで、この動的二値化処理において、原
画像中の高輝度部分を消去する前処理の意義について説
明する。図９（ｂ）は、同様の細胞画像を対象として、
このような前処理を行わない従来の動的二値化処理によ
って得られた二値化画像を参照のために示したものであ
る。この場合においても、細胞本体２４ａの周囲には透
過照明光の回折によって高輝度部分２３が輪状に現れて
いる。
画像中の高輝度部分を消去する前処理の意義について説
明する。図９（ｂ）は、同様の細胞画像を対象として、
このような前処理を行わない従来の動的二値化処理によ
って得られた二値化画像を参照のために示したものであ
る。この場合においても、細胞本体２４ａの周囲には透
過照明光の回折によって高輝度部分２３が輪状に現れて
いる。
【００３６】そしてこの高輝度部分２３の外側近傍の画
素を動的二値化処理の対象とする場合において、この範
囲にノイズなどにより輝度値が幾分周囲よりも低い部分
（グレーゾーン）が存在する場合には、以下に説明する
不具合が生じる。すなわち二値化しきい値を設定するた
めに参照範囲Ｚ内の複数画素の輝度値の平均値を求める
と、高輝度部分２３が参照範囲Ｚ内に含まれていること
から、平均値は高輝度部分２３が存在しない場合と比較
して高輝度側に偏る。この結果、設定された二値化しき
い値が上述のグレーゾーンの輝度値よりも高くなる場合
が発生する。
素を動的二値化処理の対象とする場合において、この範
囲にノイズなどにより輝度値が幾分周囲よりも低い部分
（グレーゾーン）が存在する場合には、以下に説明する
不具合が生じる。すなわち二値化しきい値を設定するた
めに参照範囲Ｚ内の複数画素の輝度値の平均値を求める
と、高輝度部分２３が参照範囲Ｚ内に含まれていること
から、平均値は高輝度部分２３が存在しない場合と比較
して高輝度側に偏る。この結果、設定された二値化しき
い値が上述のグレーゾーンの輝度値よりも高くなる場合
が発生する。
【００３７】そしてこのような場合には、グレーゾーン
を二値化対象画素とした動的二値化処理においては、本
来ならば輝度値「１」に区分されるべき場合にあって
も、二値化しきい値が高輝度側に偏って設定されている
ことから、当該画素は輝度値「０」に区分される。この
ため、高輝度部分２３の近傍にグレーゾーンが存在する
場合には、二値化画像では本来存在しない暗部２５とし
て現れ、形状抽出においてこのような暗部２５は突起２
４ｂと誤認されやすいことから、突起検出精度を低下さ
せることとなっていた。
を二値化対象画素とした動的二値化処理においては、本
来ならば輝度値「１」に区分されるべき場合にあって
も、二値化しきい値が高輝度側に偏って設定されている
ことから、当該画素は輝度値「０」に区分される。この
ため、高輝度部分２３の近傍にグレーゾーンが存在する
場合には、二値化画像では本来存在しない暗部２５とし
て現れ、形状抽出においてこのような暗部２５は突起２
４ｂと誤認されやすいことから、突起検出精度を低下さ
せることとなっていた。
【００３８】これに対し、本実施の形態に示す細胞画像
の画像処理においては、前述のように原画像中の高輝度
部分２３を予め消去した前処理画像を動的二値化処理の
対象としていることから、上述の高輝度部分２３の存在
に起因する二値化処理上の誤差を生じることなく、正し
い形状抽出を行うことができる。
の画像処理においては、前述のように原画像中の高輝度
部分２３を予め消去した前処理画像を動的二値化処理の
対象としていることから、上述の高輝度部分２３の存在
に起因する二値化処理上の誤差を生じることなく、正し
い形状抽出を行うことができる。
【００３９】次に図１０を参照して、サーチライン設定
処理について図４のフロー図に沿って説明する。図４の
フロー図は、サーチライン設定処理部１５の機能を実現
するためのサーチライン設定プログラムを示している。
まず図１０（ａ）に示すように、原画像より突起２４ｂ
を消去し細胞本体２４ａのみとした画像（細胞本体画
像）を細胞本体検出処理部１４より入手する（ＳＴ２
１）。この細胞本体画像は、細胞本体検出処理部１４に
て生成され、原画像中の細胞画像をウェーブレット変換
によって周波数変換して得られる低周波側の成分から細
胞本体画像が生成される。なお、周波数変換としては、
ウェーブレット変換以外にもアダマール変換や離散コサ
イン変換（ＤＣＴ）、フーリエ変換等を使用したものを
用いてもよい。
処理について図４のフロー図に沿って説明する。図４の
フロー図は、サーチライン設定処理部１５の機能を実現
するためのサーチライン設定プログラムを示している。
まず図１０（ａ）に示すように、原画像より突起２４ｂ
を消去し細胞本体２４ａのみとした画像（細胞本体画
像）を細胞本体検出処理部１４より入手する（ＳＴ２
１）。この細胞本体画像は、細胞本体検出処理部１４に
て生成され、原画像中の細胞画像をウェーブレット変換
によって周波数変換して得られる低周波側の成分から細
胞本体画像が生成される。なお、周波数変換としては、
ウェーブレット変換以外にもアダマール変換や離散コサ
イン変換（ＤＣＴ）、フーリエ変換等を使用したものを
用いてもよい。
【００４０】そしてこの細胞本体画像から突起の根本を
検出するためのサーチラインＳＬが設定される。すなわ
ち図１０（ｂ）に示すように、この細胞本体画像におい
て細胞本体２４ａを数画素太らせた画像（細胞本体２４
ａ’の画像）を生成して、サーチライン記憶部１６に格
納し（ＳＴ２２）、メインフローに戻る。そしてここで
求められた拡大された細胞本体２４ａ’の外側のエッジ
が、サーチラインＳＬとして設定される。すなわち細胞
本体画像の細胞本体２４ａのエッジを所定画素外側に拡
大して得られる細胞本体のエッジが突起の根本を検出す
るためのサーチラインＳＬとして設定される。
検出するためのサーチラインＳＬが設定される。すなわ
ち図１０（ｂ）に示すように、この細胞本体画像におい
て細胞本体２４ａを数画素太らせた画像（細胞本体２４
ａ’の画像）を生成して、サーチライン記憶部１６に格
納し（ＳＴ２２）、メインフローに戻る。そしてここで
求められた拡大された細胞本体２４ａ’の外側のエッジ
が、サーチラインＳＬとして設定される。すなわち細胞
本体画像の細胞本体２４ａのエッジを所定画素外側に拡
大して得られる細胞本体のエッジが突起の根本を検出す
るためのサーチラインＳＬとして設定される。
【００４１】なお、ここでサーチラインＳＬを求める処
理として、図１０（ｃ）に示すように細胞本体画像にお
いて細胞本体２４ａのエッジの数画素外側に直接サーチ
ラインＳＬを設定し、これをサーチライン記憶部１６に
格納するようにしてもよい（ＳＴ２２’）。すなわちこ
の場合にはサーチラインＳＬは、細胞本体画像の細胞本
体２４ａのエッジから、根本検出のサーチのために適切
な距離に相当する所定画素だけ外側に設定される。
理として、図１０（ｃ）に示すように細胞本体画像にお
いて細胞本体２４ａのエッジの数画素外側に直接サーチ
ラインＳＬを設定し、これをサーチライン記憶部１６に
格納するようにしてもよい（ＳＴ２２’）。すなわちこ
の場合にはサーチラインＳＬは、細胞本体画像の細胞本
体２４ａのエッジから、根本検出のサーチのために適切
な距離に相当する所定画素だけ外側に設定される。
【００４２】次に図１１、図１２、図１３を参照して、
突起検出処理について図５のフロー図に沿って説明す
る。図５のフロー図は、突起検出処理部１７の機能を実
現するための突起検出処理プログラムを示している。こ
の突起検出処理は、図１１に示すように二値化処理によ
って求められた二値化画像（探査用画像）において、サ
ーチラインＳＬに該当する位置（画素）に沿ってエッジ
をサーチして突起の根本を検出し、検出された根本から
突起の先端へ向かって連続するエッジを検出するもので
ある。これらの処理は、サーチラインＳＬに該当する位
置（画素）に沿って注目画素を移動させながら行われ
る。そして以下に示すステップが注目画素に対して適用
される。
突起検出処理について図５のフロー図に沿って説明す
る。図５のフロー図は、突起検出処理部１７の機能を実
現するための突起検出処理プログラムを示している。こ
の突起検出処理は、図１１に示すように二値化処理によ
って求められた二値化画像（探査用画像）において、サ
ーチラインＳＬに該当する位置（画素）に沿ってエッジ
をサーチして突起の根本を検出し、検出された根本から
突起の先端へ向かって連続するエッジを検出するもので
ある。これらの処理は、サーチラインＳＬに該当する位
置（画素）に沿って注目画素を移動させながら行われ
る。そして以下に示すステップが注目画素に対して適用
される。
【００４３】図５においてまず画素Ａ［ｉ］（図１２参
照）を注目画素とする（ＳＴ３１）。次いで注目画素Ａ
［ｉ］の輝度を検出する（ＳＴ３２）。そしてこの注目
画素Ａ［ｉ］が根本条件に合致しているか否かを判断す
る（ＳＴ３３）。根本検出条件とは、サーチラインＳＬ
上で注目画素を移動させて突起を検出する過程におい
て、当該注目画素が細胞本体から外側に突出する突起の
根本のエッジに到達したか否かを判定するための条件で
あり、具体的には図１２（ａ）に示す点Ｐ１から点Ｐ２
に注目画素が移動した場合のように、現在の注目画素
（点Ｐ２）の輝度値「１」であり、かつサーチラインＳ
Ｌ上の１つ手前の注目画素（点Ｐ１）の輝度値が「０」
となる場合、現在の注目画素（点Ｐ２）を根本検出条件
に合致するものとする。
照）を注目画素とする（ＳＴ３１）。次いで注目画素Ａ
［ｉ］の輝度を検出する（ＳＴ３２）。そしてこの注目
画素Ａ［ｉ］が根本条件に合致しているか否かを判断す
る（ＳＴ３３）。根本検出条件とは、サーチラインＳＬ
上で注目画素を移動させて突起を検出する過程におい
て、当該注目画素が細胞本体から外側に突出する突起の
根本のエッジに到達したか否かを判定するための条件で
あり、具体的には図１２（ａ）に示す点Ｐ１から点Ｐ２
に注目画素が移動した場合のように、現在の注目画素
（点Ｐ２）の輝度値「１」であり、かつサーチラインＳ
Ｌ上の１つ手前の注目画素（点Ｐ１）の輝度値が「０」
となる場合、現在の注目画素（点Ｐ２）を根本検出条件
に合致するものとする。
【００４４】そして根本検出条件に合致しているなら
ば、当該注目画素Ａ［ｉ］を突起のエッジとして突起情
報記憶部１８に記憶させる。また（ＳＴ３３）にて根本
検出条件に合致していない場合には、サーチラインＳＬ
に沿って注目画素を１つだけ移動する（ＳＴ３９）。そ
して突起の根本が検出されたならば、この注目画素に連
続する突起のエッジを検出する（ＳＴ３５）。すなわち
図１２（ｂ）に示すように突起２４ｂの根本として検出
された注目画素（矢印（イ））を起点として、突起２４
ｂの輪郭であるエッジを連続して検出するエッジ追跡が
行われる。ここでは、図１３に示すように注目画素Ａ
［ｉ］の廻りを時計回りにサーチし、輝度値が「０」→
「１」となる画素を見いだすことにより、エッジに該当
する画素を検出する。図１３においてＷが新たなエッジ
として検出される。
ば、当該注目画素Ａ［ｉ］を突起のエッジとして突起情
報記憶部１８に記憶させる。また（ＳＴ３３）にて根本
検出条件に合致していない場合には、サーチラインＳＬ
に沿って注目画素を１つだけ移動する（ＳＴ３９）。そ
して突起の根本が検出されたならば、この注目画素に連
続する突起のエッジを検出する（ＳＴ３５）。すなわち
図１２（ｂ）に示すように突起２４ｂの根本として検出
された注目画素（矢印（イ））を起点として、突起２４
ｂの輪郭であるエッジを連続して検出するエッジ追跡が
行われる。ここでは、図１３に示すように注目画素Ａ
［ｉ］の廻りを時計回りにサーチし、輝度値が「０」→
「１」となる画素を見いだすことにより、エッジに該当
する画素を検出する。図１３においてＷが新たなエッジ
として検出される。
【００４５】このエッジ検出の過程において、検出対象
とする突起２４ｂの輪郭が全て追跡されたか否かを確認
する目的で、以下のステップが実行される。すなわちエ
ッジが検出される度に、検出されたエッジはサーチライ
ンＳＬに該当するか否かを判断する（ＳＴ３６）。エッ
ジ追跡が未だ完了しておらず注目画素がエッジ上の矢印
（ロ）の位置にあって、サーチラインＳＬに該当しない
ならば、検出した当該エッジの画素を次の処理の対象と
なる注目画素とし（ＳＴ３７）、（ＳＴ３４）に戻る。
そして以降のステップを反復し、エッジ追跡を継続して
行う。このようにして（ＳＴ３４）〜（ＳＴ３７）のス
テップを反復することにより、突起２４ｂの根本から先
端に向かって連続するエッジが検出される。
とする突起２４ｂの輪郭が全て追跡されたか否かを確認
する目的で、以下のステップが実行される。すなわちエ
ッジが検出される度に、検出されたエッジはサーチライ
ンＳＬに該当するか否かを判断する（ＳＴ３６）。エッ
ジ追跡が未だ完了しておらず注目画素がエッジ上の矢印
（ロ）の位置にあって、サーチラインＳＬに該当しない
ならば、検出した当該エッジの画素を次の処理の対象と
なる注目画素とし（ＳＴ３７）、（ＳＴ３４）に戻る。
そして以降のステップを反復し、エッジ追跡を継続して
行う。このようにして（ＳＴ３４）〜（ＳＴ３７）のス
テップを反復することにより、突起２４ｂの根本から先
端に向かって連続するエッジが検出される。
【００４６】また（ＳＴ３６）において注目画素がサー
チラインＳＬに該当するならば、すなわちエッジとサー
チラインＳＬの交点を示す矢印（ハ）に注目画素が戻っ
たならば、注目画素を根本検出条件に合致した最新の画
素（矢印（イ））に戻す（ＳＴ３８）。そしてこの後サ
ーチラインＳＬに沿って注目画素を１つだけ移動させて
（ＳＴ３９）、新たな突起の根本検出を行う。すなわち
エッジ検出工程で検出したエッジの軌跡が、突起２４ｂ
の縁を周回して細胞本体２４ａ側に向かって戻り、サー
チラインＳＬに到達したならば、エッジ検出工程を終了
して根本検出工程を再開する。
チラインＳＬに該当するならば、すなわちエッジとサー
チラインＳＬの交点を示す矢印（ハ）に注目画素が戻っ
たならば、注目画素を根本検出条件に合致した最新の画
素（矢印（イ））に戻す（ＳＴ３８）。そしてこの後サ
ーチラインＳＬに沿って注目画素を１つだけ移動させて
（ＳＴ３９）、新たな突起の根本検出を行う。すなわち
エッジ検出工程で検出したエッジの軌跡が、突起２４ｂ
の縁を周回して細胞本体２４ａ側に向かって戻り、サー
チラインＳＬに到達したならば、エッジ検出工程を終了
して根本検出工程を再開する。
【００４７】そしてサーチラインＳＬに沿って順次移動
させた注目画素が、画素Ａ［ｉ］であるか否か、すなわ
ちサーチラインＳＬ上の起点に戻ったか否かを判定し、
戻っていなければ（ＳＴ３２）に戻って同一サーチライ
ンＳＬ上で新たな突起を検出するためのステップを反復
する。また起点に戻ったならば、当該サーチラインＳＬ
についての突起検出処理を終了しメインフローに戻る。
させた注目画素が、画素Ａ［ｉ］であるか否か、すなわ
ちサーチラインＳＬ上の起点に戻ったか否かを判定し、
戻っていなければ（ＳＴ３２）に戻って同一サーチライ
ンＳＬ上で新たな突起を検出するためのステップを反復
する。また起点に戻ったならば、当該サーチラインＳＬ
についての突起検出処理を終了しメインフローに戻る。
【００４８】上記説明したように、この細胞画像の画像
処理におけるエッジ検出工程は、細胞本体２４ａからス
タートして細胞本体２４ａから離れ、再び細胞本体２４
ａに戻る連続したエッジを検出するものである。そして
細胞本体２４ａからスタートして再び細胞本体２４ａに
戻る形態の１つとして、ここでは細胞本体の外側に設定
されたサーチラインＳＬからスタートしてサーチライン
ＳＬに戻る形態となっている。
処理におけるエッジ検出工程は、細胞本体２４ａからス
タートして細胞本体２４ａから離れ、再び細胞本体２４
ａに戻る連続したエッジを検出するものである。そして
細胞本体２４ａからスタートして再び細胞本体２４ａに
戻る形態の１つとして、ここでは細胞本体の外側に設定
されたサーチラインＳＬからスタートしてサーチライン
ＳＬに戻る形態となっている。
【００４９】次に図１４を参照して、突起判定処理につ
いて図６のフロー図に沿って説明する。ここでは、突起
検出処理において検出された突起が、検出目的に該当す
る突起であるか否かが判定される。まず突起２４ｂの連
続するエッジのデータを抽出し（ＳＴ４１）、連続する
エッジの軌跡は所定距離内で隣接しているか否かを判断
する（ＳＴ４２）。
いて図６のフロー図に沿って説明する。ここでは、突起
検出処理において検出された突起が、検出目的に該当す
る突起であるか否かが判定される。まず突起２４ｂの連
続するエッジのデータを抽出し（ＳＴ４１）、連続する
エッジの軌跡は所定距離内で隣接しているか否かを判断
する（ＳＴ４２）。
【００５０】すなわち図１４に示すように、突起２４ｂ
の縁に沿ったエッジが連続して検出された軌跡のうち、
細胞本体２４ａから遠ざかる方向の往路の軌跡と細胞本
体２４ａ側に向かって戻る復路の軌跡との間の距離Ｄが
所定距離以内であって隣接していると判断されたなら
ば、当該エッジは検出目的に該当する突起であると判定
する。そしてさらに未判定のエッジがあるか否かを判断
し（ＳＴ４４）、未判定のエッジがある場合には、（Ｓ
Ｔ４１）に戻って、当該エッジについて同様の突起判定
処理を行う。また未判定のエッジが存在しないならば、
メインフローに戻る。
の縁に沿ったエッジが連続して検出された軌跡のうち、
細胞本体２４ａから遠ざかる方向の往路の軌跡と細胞本
体２４ａ側に向かって戻る復路の軌跡との間の距離Ｄが
所定距離以内であって隣接していると判断されたなら
ば、当該エッジは検出目的に該当する突起であると判定
する。そしてさらに未判定のエッジがあるか否かを判断
し（ＳＴ４４）、未判定のエッジがある場合には、（Ｓ
Ｔ４１）に戻って、当該エッジについて同様の突起判定
処理を行う。また未判定のエッジが存在しないならば、
メインフローに戻る。
【００５１】これに対し（ＳＴ４２）において、図１４
に示す細胞本体２４ａから外側へ幾分突出した凸部２４
ｃのように、軌跡形状が緩やかな弧状であり往路の軌跡
と復路の軌跡が離れているような場合には、すなわち連
続するエッジの軌跡が隣接していないならば、当該エッ
ジは検出目的に該当せずと判定して当該エッジのデータ
を突起情報記憶部１８から削除し（ＳＴ４３）、この後
（ＳＴ４４）に進む。
に示す細胞本体２４ａから外側へ幾分突出した凸部２４
ｃのように、軌跡形状が緩やかな弧状であり往路の軌跡
と復路の軌跡が離れているような場合には、すなわち連
続するエッジの軌跡が隣接していないならば、当該エッ
ジは検出目的に該当せずと判定して当該エッジのデータ
を突起情報記憶部１８から削除し（ＳＴ４３）、この後
（ＳＴ４４）に進む。
【００５２】すなわち、ここでは細胞本体から外側へ突
出した部分のエッジのうち、エッジの往復の軌跡が互い
に所定距離内で隣接したものを突起のエッジとする。こ
れにより、突起抽出処理において本来の検出目的に合致
する形状の突起のみを正しく抽出することができる。
出した部分のエッジのうち、エッジの往復の軌跡が互い
に所定距離内で隣接したものを突起のエッジとする。こ
れにより、突起抽出処理において本来の検出目的に合致
する形状の突起のみを正しく抽出することができる。
【００５３】上記説明したように、細胞画像から神経突
起などの突起部分のみを検出する突起抽出処理におい
て、細胞本体の周囲のみを突起検出のための探査対象と
することにより、画像に含まれるノイズの影響を低減し
て、突起抽出を適切にかつ効率よく行うことができる。
また、探査用画像の生成において、動的二値化処理に先
立って撮像時に照明光の回折によって細胞本体の周囲に
生じた高輝度部分を検出し、検出した高輝度部分の輝度
の値を所定の値に置換してこれらの高輝度部分を消去す
る処理を行うようにしたので、高輝度部分の影響を排除
して正しい形状抽出を行うことができる。
起などの突起部分のみを検出する突起抽出処理におい
て、細胞本体の周囲のみを突起検出のための探査対象と
することにより、画像に含まれるノイズの影響を低減し
て、突起抽出を適切にかつ効率よく行うことができる。
また、探査用画像の生成において、動的二値化処理に先
立って撮像時に照明光の回折によって細胞本体の周囲に
生じた高輝度部分を検出し、検出した高輝度部分の輝度
の値を所定の値に置換してこれらの高輝度部分を消去す
る処理を行うようにしたので、高輝度部分の影響を排除
して正しい形状抽出を行うことができる。
【００５４】
【発明の効果】本発明によれば、細胞画像から形状抽出
用の探査画像を取得するために行われる二値化処理に先
だって、光の回折によって細胞本体の周囲に生じた高輝
度部分を検出し、検出した高輝度部分の輝度の値を所定
の値に置換してこれらの高輝度部分を消去する処理を行
うようにしたので、高輝度部分の影響を排除して正しい
形状抽出を行うことができる。
用の探査画像を取得するために行われる二値化処理に先
だって、光の回折によって細胞本体の周囲に生じた高輝
度部分を検出し、検出した高輝度部分の輝度の値を所定
の値に置換してこれらの高輝度部分を消去する処理を行
うようにしたので、高輝度部分の影響を排除して正しい
形状抽出を行うことができる。
【図面の簡単な説明】
【図１】本発明の一実施の形態の画像処理装置の構成を
示すブロック図
示すブロック図
【図２】本発明の一実施の形態の細胞画像の画像処理方
法のメインフローを示すフロー図
法のメインフローを示すフロー図
【図３】本発明の一実施の形態の細胞画像の画像処理方
法における二値化処理のフロー図
法における二値化処理のフロー図
【図４】本発明の一実施の形態の細胞画像の画像処理方
法におけるサーチライン設定処理のフロー図
法におけるサーチライン設定処理のフロー図
【図５】本発明の一実施の形態の細胞画像の画像処理方
法における突起検出処理のフロー図
法における突起検出処理のフロー図
【図６】本発明の一実施の形態の細胞画像の画像処理方
法における突起判定処理のフロー図
法における突起判定処理のフロー図
【図７】本発明の一実施の形態の細胞画像を示す画像図
【図８】本発明の一実施の形態の細胞画像の画像処理方
法における二値化処理の説明図
法における二値化処理の説明図
【図９】本発明の一実施の形態の細胞画像の画像処理方
法における二値化処理の説明図
法における二値化処理の説明図
【図１０】本発明の一実施の形態の細胞画像の画像処理
方法におけるサーチライン設定処理の説明図
方法におけるサーチライン設定処理の説明図
【図１１】本発明の一実施の形態の細胞画像の画像処理
方法における突起検出処理の説明図
方法における突起検出処理の説明図
【図１２】本発明の一実施の形態の細胞画像の画像処理
方法における突起検出処理の説明図
方法における突起検出処理の説明図
【図１３】本発明の一実施の形態の細胞画像の画像処理
方法における突起検出処理の説明図
方法における突起検出処理の説明図
【図１４】本発明の一実施の形態の細胞画像の画像処理
方法における突起判定処理の説明図
方法における突起判定処理の説明図
２ マイクロタイタープレート ３ 試料 ５ カメラ １０ 画像処理部 １２ 二値化処理部 １４ 細胞本体検出処理部 １５ サーチライン設定処理部 １７ 突起検出処理部 １９ 突起判定処理部 ２４ 細胞 ２４ａ 細胞本体 ２４ｂ 突起
フロントページの続き (51)Int.Cl.7 識別記号 ＦＩ テーマコート゛(参考） // Ｃ１２Ｍ 1/34 Ｃ１２Ｍ 1/34 Ａ Ｆターム(参考） 2G045 AA40 CB01 FA16 JA01 4B029 AA07 AA27 BB11 CC02 FA15 4B063 QA05 QQ08 QR69 QS24 QS39 QX01 5B057 AA10 BA02 CA08 CA12 CA16 CB08 CB12 CB16 CC03 CE12 CF10 DA08 DB02 DB09 DC16 DC22
Claims (10)
- 【請求項１】透過照明にて得られた細胞画像において光
の回折によって細胞本体の周囲に生じた高輝度部分を検
出する工程と、検出した前記高輝度部分の輝度の値を所
定の値に置換してこれらの高輝度部分を消去する工程
と、前記高輝度部分が消去された細胞画像に対して動的
二値化処理を行う工程とを含むことを特徴とする細胞画
像の画像処理方法。 - 【請求項２】前記高輝度部分の検出は、少なくとも動的
二値化処理の対象となる複数の画素の輝度の分布より高
輝度検出用のしきい値を設定し、輝度がこのしきい値を
超えた画素を高輝度部分とすることを特徴とする請求項
１記載の細胞画像の画像処理方法。 - 【請求項３】前記高輝度部分の輝度の値と置換する所定
の値は、少なくとも動的二値化処理の対象となる複数の
画素の輝度の値の平均値であることを特徴とする請求項
１記載の細胞画像の画像処理方法。 - 【請求項４】前記動的二値化処理は、二値化処理の対象
となる一つの画素とその周囲の複数の画素の輝度の分布
から二値化しきい値を設定する処理と、前記二値化処理
の対象となる一つの画素の輝度の値と前記二値化しきい
値とを比較する処理とを含むことを特徴とする請求項１
記載の細胞画像の画像処理方法。 - 【請求項５】前記二値化しきい値が、二値化処理の対象
となる一つの画素とその周囲の複数の画素の輝度の値の
平均値であることを特徴とする請求項４記載の細胞画像
の画像処理方法。 - 【請求項６】コンピュータプログラムによって作動する
処理装置に細胞画像の画像処理を実現させるためのプロ
グラムであって、透過照明にて得られた細胞画像におい
て光の回折によって細胞本体の周囲に生じた高輝度部分
を検出する工程と、検出した前記高輝度部分の輝度の値
を所定の値に置換してこれらの高輝度部分を消去する工
程と、前記高輝度部分が消去された細胞画像に対して動
的二値化処理を行う工程とを含むことを特徴とする細胞
画像の画像処理を実現させるためのプログラム。 - 【請求項７】前記高輝度部分の検出は、少なくとも動的
二値化処理の対象となる複数の画素の輝度の分布より高
輝度検出用のしきい値を設定し、輝度がこのしきい値を
超えた画素を高輝度部分とすることを特徴とする請求項
６記載の細胞画像の画像処理を実現させるためのプログ
ラム。 - 【請求項８】前記高輝度部分の輝度の値と置換する所定
の値は、少なくとも動的二値化処理の対象となる複数の
画素の輝度の値の平均値であることを特徴とする請求項
６記載の細胞画像の画像処理を実現させるためのプログ
ラム。 - 【請求項９】前記動的二値化処理は、二値化処理の対象
となる一つの画素とその周囲の複数の画素の輝度の分布
から二値化しきい値を設定する処理と、前記二値化処理
の対象となる一つの画素の輝度の値と前記二値化しきい
値とを比較する処理とを含むことを特徴とする請求項６
記載の細胞画像の画像処理を実現させるためのプログラ
ム。 - 【請求項１０】前記二値化しきい値が、二値化処理の対
象となる一つの画素とその周囲の複数の画素の輝度の値
の平均値であることを特徴とする請求項９記載の細胞画
像の画像処理を実現させるためのプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001113692A JP2002312761A (ja) | 2001-04-12 | 2001-04-12 | 細胞画像の画像処理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001113692A JP2002312761A (ja) | 2001-04-12 | 2001-04-12 | 細胞画像の画像処理方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2002312761A true JP2002312761A (ja) | 2002-10-25 |
Family
ID=18964884
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001113692A Pending JP2002312761A (ja) | 2001-04-12 | 2001-04-12 | 細胞画像の画像処理方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2002312761A (ja) |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005333889A (ja) * | 2004-05-27 | 2005-12-08 | Cybox Co Ltd | 細胞壊死方法及び細胞壊死装置 |
| WO2006051813A1 (ja) * | 2004-11-09 | 2006-05-18 | Hitachi Medical Corporation | 細胞培養装置、画像処理装置及び細胞検出システム |
| WO2006104201A1 (ja) * | 2005-03-29 | 2006-10-05 | Olympus Corporation | 細胞画像解析方法、細胞画像解析プログラム、細胞画像解析装置、スクリーニング方法およびスクリーニング装置 |
| JP2007020449A (ja) * | 2005-07-14 | 2007-02-01 | Olympus Corp | スクリーニング方法およびスクリーニング装置 |
| WO2007069414A1 (ja) * | 2005-12-14 | 2007-06-21 | The University Of Tokyo | 線状の形態を有する細胞等を解析する方法及び神経細胞解析方法並びにそれら方法を実行する装置及びプログラム |
| WO2010074179A1 (ja) * | 2008-12-24 | 2010-07-01 | Ｔａｋｕｍｉ Ｖｉｓｉｏｎ株式会社 | 画像処理方法及びコンピュータプログラム |
| JP2012039927A (ja) * | 2010-08-18 | 2012-03-01 | Univ Of Fukui | 細胞画像解析装置及び方法 |
| JP2013094127A (ja) * | 2011-11-01 | 2013-05-20 | Hioki Ee Corp | 光透過性粒子の撮像方法および光透過性粒子の撮像装置 |
| JP2013236564A (ja) * | 2012-05-11 | 2013-11-28 | Nikon Corp | 細胞評価装置、細胞評価方法およびプログラム |
| WO2016027542A1 (ja) * | 2014-08-18 | 2016-02-25 | 株式会社Ｓｃｒｅｅｎホールディングス | 閾値決定方法、画像処理方法および画像処理装置 |
| CN105463054A (zh) * | 2016-01-29 | 2016-04-06 | 广州恒迪生物科技有限公司 | 一种肿瘤药物敏感性试验抑制率计算方法 |
| JP2016214269A (ja) * | 2016-09-29 | 2016-12-22 | 株式会社ニコン | 細胞評価装置、細胞評価方法およびプログラム |
| WO2018084002A1 (ja) * | 2016-11-01 | 2018-05-11 | 倉敷紡績株式会社 | 細胞測定方法 |
| JP2019158352A (ja) * | 2018-03-07 | 2019-09-19 | 株式会社島津製作所 | 粒子画像解析装置及び粒子画像解析方法 |
| JP2021027836A (ja) * | 2020-11-16 | 2021-02-25 | 株式会社ニコン | 細胞評価装置、インキュベータおよびプログラム |
-
2001
- 2001-04-12 JP JP2001113692A patent/JP2002312761A/ja active Pending
Cited By (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005333889A (ja) * | 2004-05-27 | 2005-12-08 | Cybox Co Ltd | 細胞壊死方法及び細胞壊死装置 |
| WO2006051813A1 (ja) * | 2004-11-09 | 2006-05-18 | Hitachi Medical Corporation | 細胞培養装置、画像処理装置及び細胞検出システム |
| CN101048492B (zh) * | 2004-11-09 | 2013-01-09 | 株式会社钟化 | 细胞培养装置、图像处理装置和细胞检测系统 |
| US8064661B2 (en) | 2004-11-09 | 2011-11-22 | Kaneka Corporation | Cell culture device, image processing device and cell detecting system |
| WO2006104201A1 (ja) * | 2005-03-29 | 2006-10-05 | Olympus Corporation | 細胞画像解析方法、細胞画像解析プログラム、細胞画像解析装置、スクリーニング方法およびスクリーニング装置 |
| JP2007020449A (ja) * | 2005-07-14 | 2007-02-01 | Olympus Corp | スクリーニング方法およびスクリーニング装置 |
| US8077945B2 (en) | 2005-12-14 | 2011-12-13 | The University Of Tokyo | Method of analyzing cell or the like having linear shape, method of analyzing nerve cell and apparatus and program for performing these methods |
| WO2007069414A1 (ja) * | 2005-12-14 | 2007-06-21 | The University Of Tokyo | 線状の形態を有する細胞等を解析する方法及び神経細胞解析方法並びにそれら方法を実行する装置及びプログラム |
| JPWO2007069414A1 (ja) * | 2005-12-14 | 2009-05-21 | 国立大学法人 東京大学 | 線状の形態を有する細胞等を解析する方法及び神経細胞解析方法並びにそれら方法を実行する装置及びプログラム |
| US8660379B2 (en) | 2008-12-24 | 2014-02-25 | Rohm Co., Ltd. | Image processing method and computer program |
| JP2010152536A (ja) * | 2008-12-24 | 2010-07-08 | Rohm Co Ltd | 画像処理方法及びコンピュータプログラム |
| WO2010074179A1 (ja) * | 2008-12-24 | 2010-07-01 | Ｔａｋｕｍｉ Ｖｉｓｉｏｎ株式会社 | 画像処理方法及びコンピュータプログラム |
| JP2012039927A (ja) * | 2010-08-18 | 2012-03-01 | Univ Of Fukui | 細胞画像解析装置及び方法 |
| JP2013094127A (ja) * | 2011-11-01 | 2013-05-20 | Hioki Ee Corp | 光透過性粒子の撮像方法および光透過性粒子の撮像装置 |
| JP2013236564A (ja) * | 2012-05-11 | 2013-11-28 | Nikon Corp | 細胞評価装置、細胞評価方法およびプログラム |
| WO2016027542A1 (ja) * | 2014-08-18 | 2016-02-25 | 株式会社Ｓｃｒｅｅｎホールディングス | 閾値決定方法、画像処理方法および画像処理装置 |
| JP2016041032A (ja) * | 2014-08-18 | 2016-03-31 | 株式会社Ｓｃｒｅｅｎホールディングス | 閾値決定方法、画像処理方法および画像処理装置 |
| CN105463054A (zh) * | 2016-01-29 | 2016-04-06 | 广州恒迪生物科技有限公司 | 一种肿瘤药物敏感性试验抑制率计算方法 |
| JP2016214269A (ja) * | 2016-09-29 | 2016-12-22 | 株式会社ニコン | 細胞評価装置、細胞評価方法およびプログラム |
| WO2018084002A1 (ja) * | 2016-11-01 | 2018-05-11 | 倉敷紡績株式会社 | 細胞測定方法 |
| CN109890975A (zh) * | 2016-11-01 | 2019-06-14 | 仓敷纺绩株式会社 | 细胞测定方法 |
| JPWO2018084002A1 (ja) * | 2016-11-01 | 2019-09-19 | 倉敷紡績株式会社 | 細胞測定方法 |
| US10920262B2 (en) | 2016-11-01 | 2021-02-16 | Kurashiki Boseki Kabushiki Kaisha | Cell measurement method |
| JP7088839B2 (ja) | 2016-11-01 | 2022-06-21 | 倉敷紡績株式会社 | 細胞測定方法 |
| JP2019158352A (ja) * | 2018-03-07 | 2019-09-19 | 株式会社島津製作所 | 粒子画像解析装置及び粒子画像解析方法 |
| JP7017100B2 (ja) | 2018-03-07 | 2022-02-08 | 株式会社島津製作所 | 粒子画像解析装置及び粒子画像解析方法 |
| JP2021027836A (ja) * | 2020-11-16 | 2021-02-25 | 株式会社ニコン | 細胞評価装置、インキュベータおよびプログラム |
| JP7235027B2 (ja) | 2020-11-16 | 2023-03-08 | 株式会社ニコン | 細胞評価装置、インキュベータおよびプログラム |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2002312761A (ja) | 細胞画像の画像処理方法 | |
| US7133572B2 (en) | Fast two dimensional object localization based on oriented edges | |
| EP3477586A1 (en) | Image processing device, image processing method, and image processing program | |
| US8300937B2 (en) | Image segmentation method and system based on region features of pixels | |
| JP2007299144A (ja) | ロゴ判別装置および方法ならびにプログラム | |
| CN115100104A (zh) | 玻璃油墨区的缺陷检测方法、装置、设备及可读存储介质 | |
| JP2007025902A (ja) | 画像処理装置、画像処理方法 | |
| CN111476074A (zh) | 基于毫米波图像的人体异物检测方法 | |
| WO2001008098A1 (en) | Object extraction in images | |
| JP3743316B2 (ja) | 細胞の突起検出方法 | |
| US9305235B1 (en) | System and method for identifying and locating instances of a shape under large variations in linear degrees of freedom and/or stroke widths | |
| JP3743317B2 (ja) | 細胞の突起検出方法 | |
| CN116596921B (zh) | 一种焚烧炉渣分选方法及系统 | |
| CN1497494A (zh) | 用于低景深图像分割的方法和装置 | |
| CN112435256A (zh) | 基于图像的cnv活动性病灶检测方法及装置、电子设备 | |
| JPH06207909A (ja) | 表面欠陥検査装置 | |
| CN115546219A (zh) | 检测板式生成方法、板卡缺陷检测方法、装置及产品 | |
| JP2001311730A (ja) | 細胞系譜抽出方法 | |
| CN108875612B (zh) | 一种足迹图像的处理方法 | |
| CN112634266B (zh) | 喉镜图像的半自动标注方法、介质、设备及装置 | |
| US20040037475A1 (en) | Method and apparatus for processing annotated screen capture images by automated selection of image regions | |
| CN113793316A (zh) | 一种超声扫查区域提取方法、装置、设备和存储介质 | |
| JP2001258599A (ja) | 細胞系譜抽出方法 | |
| JP2000082146A (ja) | 動画像処理装置 | |
| AU754884B2 (en) | Automatic masking of objects in images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050412 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050419 |
|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20050701 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20050816 |