(57)【特許請求の範囲】
【請求項１】
スポットを含む画像を入力する画像入力手段と、
前記画像入力手段により入力された画像に対して、カーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施すフィルタリング手段と、
前記フィルタリング手段によりフィルタリングが複数回施された各画像における輝度分布に基づいて前記スポットの特徴を定量するスポット定量手段と、
を有し、
前記カーネルが、ガウシアンカーネル、イパネクニコフカーネル、又はメキシカンハットカーネルであり、
前記フィルタリング手段は、前記カーネルのバンド幅を所定範囲内で変更してフィルタリングを複数回施し、
前記スポット定量手段は、前記フィルタリング手段によりフィルタリングが複数回施された各画像における正規化された輝度分布のピーク値を求め、該ピーク値が最大となる画像で使用されたカーネルのバンド幅を前記スポットの直径として定量する、スポット定量装置。
【請求項２】
スポットを含む画像を入力する画像入力手段と、
前記画像入力手段により入力された画像に対して、カーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施すフィルタリング手段と、
前記フィルタリング手段によりフィルタリングが複数回施された各画像における輝度分布に基づいて前記スポットの特徴を定量するスポット定量手段と、
を有し、
前記カーネルが同心円カーネルであり、
前記フィルタリング手段は、前記カーネルの内円半径を所定範囲内で変更してフィルタリングを複数回施し、
前記スポット定量手段は、前記フィルタリング手段によりフィルタリングが複数回施された各画像における正規化された輝度分布のピーク値を求め、該ピーク値が最大となる画像で使用されたカーネルの内円半径を前記スポットの半径として定量する、スポット定量装置。
【請求項３】
前記スポット定量手段は、
前記フィルタリング手段によりフィルタリングが複数回施された各画像における正規化された輝度分布のピーク値を求め、該ピーク値が最大となる画像で使用されたカーネルのパラメータの値に基づいて前記スポットの特徴を定量する、請求項１又は２に記載のスポット定量装置。
【請求項４】
前記フィルタリング手段は、前記カーネルのパラメータを、前記スポットの特徴量に対応するパラメータの値を含む所定範囲内で変更してフィルタリングを施す、請求項１乃至３のいずれか１項に記載のスポット定量装置。
【請求項５】
前記フィルタリング手段は、
前記画像入力手段により入力された画像に対して所定のカーネルを使用してフィルタリングを施す第１のフィルタリング手段と、
前記第１のフィルタリング手段によりフィルタリングが施された画像をスポットの中心間の距離に基づいてスポットを含む画像領域と含まない画像領域とに分割する画像領域分割手段と、
前記画像領域分割手段により分割されたスポットを含む画像領域におけるスポットの中心をスポット位置として検出するスポット検出手段と、
前記スポット検出手段により検出されたスポット位置を中心とした所定の画像領域を前記画像入力手段により入力された画像から抽出し、抽出した画像領域に対してカーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施す第２のフィルタリング手段と、
を有する、請求項１乃至４のいずれか１項に記載のスポット定量装置。
【請求項６】
前記画像領域分割手段は、互いの中心点が所定距離内にあるスポットを同一の画像領域に帰属させ、互いの中心点が所定距離より離れたスポットを異なる画像領域に帰属させることで、前記第１のフィルタリング手段によりフィルタリングが施された画像をスポットを含む画像領域と含まない画像領域とに分割する、請求項５に記載のスポット定量装置。
【請求項７】
前記スポット検出手段は、前記画像領域分割手段により分割されたスポットを含む画像領域における各画素を輝度分布に基づいて複数のクラスタに分類し、輝度値の平均値が最も高いクラスタに帰属する画素のうち最高輝度の画素の位置をスポット位置として検出する、請求項５又は６に記載のスポット定量装置。
【請求項８】
スポットを含む画像を入力する画像入力ステップと、
入力した画像に対してカーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施すフィルタリングステップと、
複数のフィルタリングが施された各画像における輝度分布に基づいて前記スポットの特徴を定量するスポット定量ステップと、を有し、
前記カーネルが、ガウシアンカーネル、イパネクニコフカーネル、又はメキシカンハットカーネルであり、
前記フィルタリングステップでは、前記カーネルのバンド幅を所定範囲内で変更してフィルタリングを複数回施し、
前記スポット定量ステップでは、複数のフィルタリングが施された各画像における正規化された輝度分布のピーク値を求め、該ピーク値が最大となる画像で使用されたカーネルのバンド幅を前記スポットの直径として定量する、スポット定量方法。
【請求項９】
スポットを含む画像を入力する画像入力ステップと、
入力した画像に対してカーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施すフィルタリングステップと、
複数のフィルタリングが施された各画像における輝度分布に基づいて前記スポットの特徴を定量するスポット定量ステップと、を有し、
前記カーネルが同心円カーネルであり、
前記フィルタリングステップでは、前記カーネルの内円半径を所定範囲内で変更してフィルタリングを複数回施し、
前記スポット定量ステップでは、複数のフィルタリングが施された各画像における正規化された輝度分布のピーク値を求め、該ピーク値が最大となる画像で使用されたカーネルの内円半径を前記スポットの半径として定量する、スポット定量方法。
【請求項１０】
コンピュータを、
スポットを含む画像に対してカーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施すフィルタリング手段、及び
前記フィルタリング手段によりフィルタリングが複数回施された各画像領域における輝度分布に基づいて前記スポットの特徴を定量するスポット定量手段、
として機能させるためのプログラムであって、
前記カーネルが、ガウシアンカーネル、イパネクニコフカーネル、又はメキシカンハットカーネルであり、
前記フィルタリング手段は、前記カーネルのバンド幅を所定範囲内で変更してフィルタリングを複数回施し、
前記スポット定量手段は、前記フィルタリング手段によりフィルタリングが複数回施された各画像における正規化された輝度分布のピーク値を求め、該ピーク値が最大となる画像で使用されたカーネルのバンド幅を前記スポットの直径として定量する、プログラム。
【請求項１１】
コンピュータを、
スポットを含む画像に対してカーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施すフィルタリング手段、及び
前記フィルタリング手段によりフィルタリングが複数回施された各画像領域における輝度分布に基づいて前記スポットの特徴を定量するスポット定量手段、
として機能させるためのプログラムであって、
前記カーネルが同心円カーネルであり、
前記フィルタリング手段は、前記カーネルの内円半径を所定範囲内で変更してフィルタリングを複数回施し、
前記スポット定量手段は、前記フィルタリング手段によりフィルタリングが複数回施された各画像における正規化された輝度分布のピーク値を求め、該ピーク値が最大となる画像で使用されたカーネルの内円半径を前記スポットの半径として定量する、プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像に含まれるスポットの特徴を定量する技術に関する。
【０００２】
特に、蛍光抗体によって標識された遺伝子あるいは染色体に対応する複数のスポットを含む画像におけるスポットの特徴を定量する技術に関する。
【背景技術】
【０００３】
バイオテクノロジーに基づいた医療分野においては、遺伝子あるいは染色体を蛍光抗体によって標識し、その様態を観察することにより、癌をはじめとする様々な病気の診断が行われている。
【０００４】
具体的には、蛍光抗体により標識化された遺伝子又は染色体を顕微鏡あるいは顕微鏡で観察すると、これらの遺伝子等に対応する蛍光領域（以下、「スポット」といういう）を含む画像が得られる。この画像をデジタル化してコンピュータに取り込み、診断対象とする遺伝子又は染色体に対応するスポットのみを検出し、検出したスポットの大きさや形状などの特徴を定量する。そして、定量化されたスポットの特徴に基づき、病気の診断を行う。
【０００５】
例えば、慢性骨髄性白血病の診断においては、９番染色体上にある発ガン遺伝子ｃ−ａｂｌ と２２番染色体上にある ｂｃｒ遺伝子との間での転座の有無が診断の重要な指標となるので、これらの遺伝子に対応するスポットの特徴量に基づき転座の有無を判断し、診断を行う。
【０００６】
検査技師が画像を肉眼で観察することにより、診断対象の遺伝子に対応するスポットを検出し、その特徴の定量を行うこともできるが、この方法では検査技師の主観が入り易く、客観性や再現性、定量性に問題があるので、病気の診断においては、スポットの特徴の定量をコンピュータ上で、自動的、客観的に行うことが試みられている。
【０００７】
コンピュータを用いたスポットの特徴の定量装置または方法としては、特許文献１に開示されているように、スポットの強度（輝度）がピークとなる点を求め、ピーク点を中心とするｎ本の線分を設定し、各線分の輝度分布について、輝度が所与の値(th)となるｎ個の点を求め、これらの点を結んだ、ｎ角形を最小楕円近似して、スポット径を測定するする装置又は方法がある。
【０００８】
また、特許文献２に開示されているように、細胞核を表す領域を残りの領域から分割し、分割した領域の輝度分布に基づき、スポット形状を定量する装置又は方法もある。
【０００９】
ここで、ノイズが生じる場合には、定量を正確に行うため、特許文献３に開示されているように所定のカーネルを使用してフィルタリングを施すことでノイズを低減する。或いは、特許文献４に開示さているように、一つのスポットを認識するための可変的なスポット認識領域を設定し、このスポット認識領域のサイズを変化させたときのシグナル強度分布の変化に基づいてスポットを検出することで、ノイズの影響を抑える。
【特許文献１】特開平６−２５９７８４号公報
【特許文献２】特表２００４−５３５５６９号公報
【特許文献３】特表２００６−５０５７８２号公報
【特許文献４】特開２００６−０８４２８１号公報
【発明の開示】
【発明が解決しようとする課題】
【００１０】
しかし、特許文献１に開示された装置又は方法では、輝度分布のピークを求める手順に加え、ｎ個の点を求める手順やｎ角形を最小楕円近似する手順を行う必要があり、特徴の測定手順は煩雑なものとなる。そして、特許文献３に開示された装置又は方法では、細胞核を表す領域が備えるパラメータのパターンについて、予め多数の実験データを取得する必要があり、簡便な方法ではない。
【００１１】
また、特許文献３又は４に開示された装置又は方法では、ノイズの影響を充分に抑制できず、ノイズや輝度の変動によりスポットの特徴の定量が不正確になるという問題があった。
【００１２】
上述した問題点を鑑み、本発明は、検出したスポットの特徴を簡易かつ正確に定量できる技術を提供することを目的とする。
【課題を解決するための手段】
【００１３】
上記目的を達成するために、本発明にかかるスポット定量装置は、スポットを含む画像を入力する画像入力手段と、前記画像入力手段により入力された画像に対して、カーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施すフィルタリング手段と、前記フィルタリング手段によりフィルタリングが複数回施された各画像における輝度分布に基づいて前記スポットの特徴を定量するスポット定量手段と、
を有する。
【００１４】
本発明にかかるスポット定量方法は、スポットを含む画像を入力し、入力した画像に対してカーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施し、複数のフィルタリングが施された各画像における輝度分布に基づいて前記スポットの特徴を定量する。
【００１５】
本発明にかかるプログラムは、コンピュータを、スポットを含む画像に対してでカーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施すフィルタリング手順、及び前記フィルタリング手順によりフィルタリングが複数回施された各画像領域における輝度分布に基づいて前記スポットの特徴を定量するスポット定量手順、として機能させる。
【発明の効果】
【００１６】
本発明によれば、検出したスポットの特徴を簡易かつ正確に定量できる技術を提供することができる。
【発明を実施するための最良の形態】
【００１７】
（第１実施形態）
本発明を実施するための第１実施形態について図１〜図９を参照して詳細に説明する。
【００１８】
図１は、本実施形態のスポット定量装置１の構成を示すブロック図である。図１を参照すると、スポット定量装置１は、画像入力部１０と、フィルタリング部１５と、特徴量算出部５０と、から構成されている。
【００１９】
画像入力部１０により、蛍光染色された遺伝子又は染色体に対応するスポットを含む画像がスポット定量装置１に入力される。
【００２０】
この画像入力部１０は、例えば、図２に示すように、光源１００と、励起フィルタ１０１と、バリアフィルタ１０２と、ＣＣＤカメラ１０３と、試料１０４と、対眼レンズ１０５と、対物レンズ１０６と、ダイクロックミラー１０７と、から構成される。
【００２１】
ユーザは、対物レンズ１０６、対眼レンズ１０５を通して、試料１０４の画像を観察し、蛍光試薬の波長特性に合致した励起フィルタ１０１及びバリアフィルタ１０２を選択する。
【００２２】
そして、光源１００から選択された励起フィルタ１０１を通して、ダイクロックミラー１０７で反射された励起光が試料１０４に照射される。この試料１０４の画像を、所定の波長域の光を分離する、ダイクロックミラー１０７及びバリアフィルタ１０２を通してＣＣＤカメラ１０３が撮影する。
【００２３】
図１に戻り、画像フィルタリング部１５は、画像入力部１０により入力された画像データに対してフィルタリングを施す。このフィルタリング部１５は、画像フィルタリング部２０と、領域分割部３０と、スポット検出部４０と、から構成される。
【００２４】
画像フィルタリング部２０は、画像入力部１０により入力された画像データに対して、所与のパラメータにより規定されるカーネル（マスク）を使用して、フィルタリングを施す。フィルタリングでは、カーネルを定義するフィルタ関数と入力画像とを畳み込み積分することでフィルタ後の画像を得る。
【００２５】
カーネルとしては、例えば、下記の（１）式又は（２）式で定義されるガウシアン（Gaussian）カーネルを用いることができる。
【００２６】
【数１】
0005347272.tif 000002
ここで、Ｋ（ｘ）はフィルタ関数、ｘは輝度の入力値、ｈはバンド幅である。
【００２７】
【数２】
0005347272.tif 000003
ここで、Ｋ（ｘ，ｙ）はフィルタ関数、ｘは輝度の入力値、ｈ1はｘ軸方向のバンド幅、ｈ2はｙ軸方向のバンド幅、ρは軸の方向を定めるパラメータである。
【００２８】
また、カーネルとして、下記の（３）式で定義されるイパネクニコフ（Epanechnikov）カーネルを用いることもできる。
【００２９】
【数３】
0005347272.tif 000004
ここで、Ｋ（ｘ）はフィルタ関数、ｘは輝度の入力値、|x|は入力ベクトルの絶対値、ｈはバンド幅である。Ｉｃは条件Ｃが成立するときは１となり、それ以外のときは０となる指示関数である。
【００３０】
また、カーネルとして、下記の（４）式で定義される同心円カーネルを用いることもできる。
【００３１】
【数４】
0005347272.tif 000005
ここで、Ｋ（ｘ）はフィルタ関数、ｘは輝度の入力値、|x|は入力ベクトルの絶対値、ｒ1は内心円の半径、ｒ2は外心円の半径である。
【００３２】
また、カーネルとして、下記の（５）式で定義されるメキシカンハットカーネルを用いることもできる。
【００３３】
【数５】
0005347272.tif 000006
ここで、Ｋ（ｘ）はフィルタ関数、ｘは輝度の入力値、|x|は入力ベクトルの絶対値、ｈはバンド幅である。
【００３４】
次に、領域分割部３０は、スポットの中心間の距離に基づいて、画像フィルタリング部２０によりフィルタリングが施された画像を、スポットを含む画像領域と含まない画像領域とに分割する。画像領域を分割する手続きの詳細については後述する。
【００３５】
スポット検出部４０は、領域分割部３０により分割されたスポットを含む画像領域に対して、各画像領域内の輝度分布を求め、そのピーク点をスポットの中心位置（以下、「スポット位置」という）として検出する。また、スポット検出部４０は、画像フィルタリング部２０によりフィルタリングが施される前の原画像を特徴量算出部５０に送信する。
【００３６】
特徴量算出部５０は、フィルタリングが施される前の原画像から、スポット検出部４０により検出されたスポット位置を中心とする所定の大きさの矩形領域を抽出する。そして、特徴量算出部５０は、抽出した矩形領域に対して、異なるパラメータにより規定されるカーネルを使用してフィルタリングを施し、スポットの特徴量を算出する。使用するカーネルとして、例えば、上述したカーネルのいずれかを選択する。スポットの特徴量の定量方法については後述する。
【００３７】
特徴量算出部５０で定量するスポットの特徴は、スポットの大きさ、面積、境界長又は形状等である。定量する形状は、長短軸の比率、長短軸のなす角度、円形度などである。
【００３８】
スポットの円形度は、例えば下記の（６）式により算出する。
【００３９】
【数６】
0005347272.tif 000007
ここで、Ｓはスポットの面積、Ｌはスポットの境界の長さである。スポット形状が真円の場合に円形度は１をとり、そうでない場合は１より小さくなる。
【００４０】
上述した画像フィルタリング部２０、領域分割部３０、スポット検出部４０、及び特徴量算出部５０は、例えば、図２に示すように、パソコン等のコンピュータ１０８および液晶表示装置等のモニタ１０９により構成される。
【００４１】
次に、スポット定量装置１の動作について図３及び図４を参照して、詳細に説明する。
【００４２】
図３は、第１実施形態における、スポット定量方法の処理手順を示すフローチャートである。
【００４３】
まず、画像入力部１０により、遺伝子又は染色体に対応するスポットを含む画像がスポット定量装置１に入力される（ステップＳ１）。
【００４４】
画像フィルタリング部２０は、入力された画像に対して上記（１）〜（５）式等により定義される所定のカーネルを使用して、フィルタリングを施す（ステップＳ２）。
【００４５】
領域分割部３０は、フィルタリングが施された画像をスポットを含む画像領域と含まない画像領域とに分割する（ステップＳ３）。
【００４６】
ステップＳ３で実行される領域分割処理について、図４を参照して詳細に説明する。図４は、領域分割方法の処理手順を示すフローチャートである。
【００４７】
同図を参照すると、まず、領域分割部３０は、領域番号ｋの初期値を０に初期化する（ステップＡ１）。この領域番号ｋは、領域分割処理（ステップＳ３）でスポットを含む領域として分割される各画像領域に対して割り当てられる整数である。
【００４８】
領域分割部３０は、入力された画像の輝度分布を求める。そして、入力画像において所定の閾値ｔ以上の輝度値を持つ画素の画素番号を集合Ｉとして保持する（ステップＡ２）。ここでの閾値ｔは、入力された画像の平均輝度に比較して充分小さな値であればよく、場合によっては０に設定することにより全ての画素の画素番号を集合Ｉに含めてもよい。集合Ｉは、入力画像のうち、スポット探索の対象とする画像領域内における、画素の画素番号の集合を意味する。また、画素番号としては、例えば、画素の座標を示す値を使用する。
【００４９】
領域分割部３０は、集合Ｉに含まれる画素の中から、最高輝度を持つ画素ｐを選択する（ステップＡ３）。
【００５０】
領域分割部３０は、ステップＡ３で選択された画素ｐを中心とした所与の半径ｒの円内に含まれる全ての画素の画素番号を求め、これらを集合Ｄp（以下、「近傍集合」という）として保持する（ステップＡ４）。この半径ｒの値は、検出対象とするスポットの大きさよりも十分大きな値であればよい。
【００５１】
領域分割部３０は、集合Ｒi（ｉ＝０，・・・，ｋ−１）と近傍集合Ｄpとが重なりを持つか否かを判別する（ステップＡ５）。ここで、Ｒi（ｉ＝０，・・・，ｋ−１）は、領域分割処理で抽出した画像領域内における、画素の画素番号の集合を意味し、ステップＡ５が最初に実行される際は（ｋ＝０のとき）、Ｒiはまだ得られていない。
【００５２】
ＲiとＤpとの間に重なりがない、若しくはＲiがまだ得られていない場合（ステップＡ５：ＮＯ）、領域分割部３０は、近傍集合Ｄpを集合Ｒkとする（ステップＡ６）。そして、領域分割部３０は、領域番号ｋの値を１だけ増加させる（ステップＡ７）。
【００５３】
ＲiとＤpとが重なりを持つ場合（ステップＡ５：ＹＥＳ）、領域分割部３０は近傍集合ＤpとＲiとの間の共通集合の要素数を重なりの大きさとして求め、この重なりの大きさが最大となるＲiにＤpを統合する（ステップＡ８）。重なりの大きさが最大となるＲiが複数ある場合、例えば、領域番号が最小のＲiを求める、あるいは領域サイズが最大であるＲiを求めるなどの方法で、統合すべきＲiを選択する。
【００５４】
ステップＡ７又はステップＡ８の後、領域分割部３０は、集合Ｉから既に探索を終えた近傍集合Ｄpを削除する（ステップＡ９）。
【００５５】
領域分割部３０は、集合Ｉに未探索の画素が含まれているか否か、即ち集合Ｉが空集合であるか否か（Ｉ＝Φ）を判別する（ステップＡ１０）。
【００５６】
集合Ｉが空集合である場合（ステップＡ１０：ＹＥＳ）、全ての領域が探索されたので領域分割部３０は領域分割処理を終了し、スポット検出部４０がステップＳ４を実行する。集合Ｉが空集合でない場合（ステップＡ１０：ＮＯ）、領域分割部３０は、ステップＡ３に戻る。
【００５７】
このように、領域分割処理（ステップＡ１〜Ａ１０）では、領域分割部３０は複数のピークのうち、互いに所定距離以上近い位置にあるピークの近傍領域を統合して一つの画像領域としている。一般に、複数のピークが存在する場合、互いに近い位置にあるピークのうち強度の比較的弱いものは、強度の比較的高いピークに対する副次的なピークか、ノイズである場合が多い。このため、統合した画像領域ごとに最高輝度を持つピークを検出することで、その領域内の副次的なピークやノイズを取り除くことができる。
【００５８】
また、領域分割処理では、領域分割部３０は互いに所定距離以上遠い位置にあるピークを異なる画像領域に帰属させている。一般に、複数のピークが存在する場合、あるピークに着目すると、そのピークから所定距離以上離れた位置にあるピークは、ノイズや副次的なピークではなく、独立した別のスポットに対応するピークである場合が多い。このため、抽出された各領域についてスポット検出を行うことで、強度が比較的弱いシグナルをもつピークであってもスポットとして検出することができる。
【００５９】
図３に戻り、スポット検出部４０は、抽出された各画像領域の輝度分布のピーク点を検出し、この点をスポット位置としてスポット検出を行う（ステップＳ４）。
【００６０】
特徴量算出部５０は、検出されたスポット位置を中心とした所定の矩形領域をフィルタリング前の原画像から抽出し、この矩形領域に対してしてパラメータをさまざまに変更したカーネルを使用してフィルタリングを施す(ステップＳ５)。
【００６１】
ここで、使用するパラメータの値としては、検出対象とするスポットの特徴量と一致するパラメータを予測し、その予測した値の前後の所定範囲内でパラメータを変化させる。例えばガウシアンカーネルの場合、検出対象とするスポットの大きさ（例えば、半径）をＲと仮定して、カーネルのパラメータ（例えば、バンド幅）をｈ（バンド幅）＝０．５Ｒからｈ＝２Ｒまで０．１ずつ変化させてフィルタリングを施す。
【００６２】
そして、特徴量算出部５０は、矩形領域におけるスポットの特徴量を算出する（ステップＳ６）。ステップＳ６の後、スポット定量装置１は、スポット定量処理を終了する。
【００６３】
具体的な特徴量の算出方法について説明する。ステップ５において、原画像から切り出した矩形領域内の輝度分布を正規化して、例えば矩形領域内の輝度の和が1となるようにする。この正規化された矩形領域に対して上述のフィルタリングを施し、フィルタリング後の輝度分布のピーク値を求める。
【００６４】
このフィルタリング及びピーク値の探索は、異なるパラメータ（例えば、バンド幅）により規定されるカーネルを用いて繰り返し行う。
【００６５】
一般に、スポットの輝度分布は、輝度が中心位置から周辺へ向けてなだらかに変化する形状を有することが多い。この輝度分布を持つスポットを含む画像に対して、単純に所定の閾値以上の輝度値の画素のみを抽出する方法を用いても、スポットの形状は正確に検出されない。そこで、例えば、（１）〜（５）式に示したように、中心位置から周辺に向けて重み付けを徐々に変化させる、フィルタ関数を使用してフィルタリングを行うことで、スポットの形状を正確に検出できる。このフィルタ関数のパラメータ（例えば、バンド幅）を変更して複数回フィルタリングを行い、スポットが最も明瞭に検出されるパラメータを求める。この最適のパラメータに対応するフィルタ関数の波形は、実際の蛍光処理された遺伝子等（スポット）の形状とよく一致するので、最適なパラメータの値から検出対象のスポットの形状が正確に求められる。
【００６６】
例えば、ガウシアンカーネル（（１）式）を使用した場合、複数回検出したピーク値の系列において、最大のピーク値が検出された際にカーネルを規定したバンド幅ｈがスポットの直径にほぼ一致する。
【００６７】
（２）式で定義されるガウシアンカーネルを使用する場合は、最大のピーク値に対応するｘ軸、ｙ軸のバンド幅（ｈ1、ｈ2）、ρは軸の方向を定めるパラメータσから、長短軸の比率、スポットの面積、長短軸のなす角度等を求めることができる。
【００６８】
（３）式、（５）式で定義されるカーネルを使用する場合も、バンド幅ｈが、スポットの直径に相当する。（４）式で定義されるカーネルを使用する場合、内円の半径ｒ1がスポットの半径にほぼ一致する。
【００６９】
このようにして、ステップ６において、この最大のピーク値に対応するパラメータを求めることで、スポットの特徴を簡易な手順により定量することができる。
【００７０】
次に、第１実施形態にかかる、スポット定量処理を実行した結果の一例を図５〜図９に示す。
【００７１】
図５は、スポット定量装置１に入力されたフィルタリング前の原画像の一例である。同図に示すように、この画像は、蛍光染色された遺伝子又は染色体に対応するスポットを１個、又は複数個含んでいる。
【００７２】
図６は、図５に示した画像におけるシグナル強度（輝度）の分布を３次元表示した図である。図５及び図６を参照すると、入力画像は、比較的強いシグナルのスポットＰ１や比較的弱いシグナルのスポットＰ２を含んでいる。そして、スポットＰ１の周辺には、ノイズが混入している。このノイズは、スポットの取り逃しや疑スポットの誤検出の要因となり得るのでノイズの低減を試みる。
【００７３】
ノイズを低減するため、例えば画像の輝度分布に基づいてスポットを検出する場合、特定の閾値以上のピークをスポットとして検出する。しかし、このノイズの強度は検出対象のスポットのシグナル強度よりも大きい場合があるので、ノイズを避けるために閾値を上げると弱いシグナルの見逃しが発生しやすくなり、逆に弱いシグナルを拾うために閾値を下げると、ノイズを拾って擬スポットの誤検出が生じ得るという問題がある。この問題は、図３に示したフィルタリング（ステップＳ２）及び図４に示した領域分割処理（ステップＡ１〜Ａ１０）を行うことで克服できる。
【００７４】
図７に、領域分割処理後の画像を示す。同図に示す、白い実践で囲まれた領域Ｉ’は、所定の閾値ｔ以上の輝度を持つ画素の画素番号の集合Ｉに対応する領域である（ステップＡ２）。白い点線で囲まれた円内の領域Ｄi’（ｉ＝１，２，３，４）は、それぞれｉ番目に輝度の高い画素の画素番号を含む近傍集合Ｄiに対応する領域である。
【００７５】
上述の領域分割処理において、まず、集合Ｉの中で最高輝度を持つ画素を中心とした半径ｒの円内の画素の画素番号を含む、近傍集合Ｄ1が保持される（ステップＡ４）。最初はＲ0がないので（ステップＡ５：ＮＯ）、このＤ1がそのままＲ0となる（ステップＡ６）。次にＤ1がＩから除かれ（ステップＡ９）、残りのＩの中で最高輝度（Ｉ全体の中では２番目に高い輝度）の画素を中心とした近傍集合Ｄ2が保持される（ステップＡ６）。Ｄ2は、既に抽出された集合Ｒ0と重なりがあるので（ステップＡ５：ＹＥＳ）、Ｒ0に統合され、更新された集合Ｒ0は、Ｒ0＝Ｄ1∪Ｄ2となる（ステップＡ８）。Ｉから、さらにＤ2が除かれ（ステップＡ９）、残りのＩの中で最高輝度（Ｉ全体の中では３番目に高い輝度）の画素を中心した近傍集合Ｄ3が保持される（ステップＡ６）。Ｄ2もＲ0と重なりがあるので（ステップＡ５：ＹＥＳ）、Ｒ0に統合され、更新された集合Ｒ0は、Ｒ0＝Ｄ1∪Ｄ2∪Ｄ3となる（ステップＡ８）。
【００７６】
Ｉから、さらにＤ3が除かれ（ステップＡ９）、残りのＩの中で最高輝度（Ｉ全体の中では４番目に高い輝度）の画素を中心した近傍集合Ｄ4が保持される（ステップＡ６）。このＤ3はＲ0と重なりがないので（ステップＡ５：ＮＯ）、このＤ4はＲ1となる（ステップＡ６）。ＩからＤ4を除くと（ステップＡ９）、Ｉが空集合となるので（ステップＡ１０：ＹＥＳ）、スポット定量装置１は領域分割処理を終了する。
【００７７】
図７を参照すると、上述の領域分割処理により、集合Ｒ0、Ｒ1に対応する２つの画像領域Ｒ0’、Ｒ1’が抽出されている。
【００７８】
そして、スポット検出部４０は、抽出された画像領域Ｒ0’、Ｒ1’における、それぞれの輝度分布のピーク点をスポット位置として検出する（ステップＳ４）。
【００７９】
画像領域Ｒ0’、Ｒ1’のそれぞれについてスポット検出を行うことで、領域Ｒ0’内の疑スポットの誤検出を抑制することができ、画像領域Ｒ1’内の弱いシグナルのスポットを確実に検出することができる。
【００８０】
続いて、特徴量算出部５０は検出されたスポット位置を中心とした所定の矩形領域を原画像から抽出し、異なるパラメータにより規定されるカーネルを使用したフィルタリングを複数回施して、スポットの特徴量を算出する（ステップＳ６）。
【００８１】
図８に、スポットの特徴を定量した結果を示す。ここでは、２つのスポットＰ６、Ｐ７が検出されたので、これらのスポットの位置を中心とした矩形領域に対し、外半径ｒ1を内半径ｒ2の１．５倍とし、ｒ2を０．０４から０．２まで０．０４刻みで変えた同心円カーネルを使用してフィルタを施している。同図を参照すると、スポットＰ６に対応する領域では、内半径（ｒ2）が０．１６のときに最大のピーク値０．７１をとり、スポットＰ７に対応する領域では、内半径（ｒ2）が０．０６のときに最大ピーク値１．７２をとっている。スポットのピーク値が最大となるときに、内半径（ｒ2）とスポットの半径とがほぼ一致するので、これらの結果から、スポットＰ６の直径として約０．３２、スポットＰ７の直径として約０．１２が算出される。これらの直径の値は、矩形領域の大きさで正規化した値である。
【００８２】
定量した特徴量が正確な値か否かを、スポットを陽に抽出することで確認する。図９（ａ）、（ｂ）にステップ４で検出されたスポットの輝度分布を２次元表示した一例を示す。同図（ａ）を参照すると、比較的面積の大きなスポットＰ６はおよそ０．３の直径を持ち、同図（ｂ）を参照すると、比較的面積の小さなスポットＰ７はおよそ０．１２の直径を持っている。よって、図８（ａ）、（ｂ）及び図９の結果より、定量したスポットの直径は、実際のスポットの直径とよく一致していることがわかる。
【００８３】
このように、フィルタを複数回施し、各画像における輝度分布のピーク値を調べるだけで、特徴量算出部５０はスポットの特徴を簡易かつ正確に定量することができる（ステップＳ６）。
【００８４】
以上説明したように、本実施形態によれば、カーネルのパラメータを変更してフィルタリングを施した複数の画像における輝度分布に基づいてスポットの特徴を定量するので、スポットの特徴を簡易かつ正確に定量することができる。
【００８５】
また、本実施形態によれば、入力された画像に対して所定のカーネルを使用してフィルタリングを施し、フィルタリング後の画像を領域分割し、領域分割された画像領域におけるスポット位置を検出し、検出した位置を中心とした所定の画像領域に対してカーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施す。このため、フィルタリングにおいてノイズがさらに低減され、ノイズや輝度の変化に対して頑強な定量が可能となる。
【００８６】
（第２実施形態）
本発明を実施するための第２実施形態について図１０及び図１１を参照して詳細に説明する。
【００８７】
図１０は、本実施形態のスポット定量装置２の構成を示すブロック図である。図１０を参照すると、スポット定量装置２の構成は、フィルタリング部１５にクラスタリング部６０を備えるほかは、第１実施形態にかかる、スポット定量装置１と同様である。
【００８８】
クラスタリング部６０は、領域分割部３０により分割された画像領域の各画素について輝度分布に基づきクラスタリングを行い、この画像領域から輝度が比較的高い領域を抽出する。
【００８９】
領域分割処理により分割された画像領域に対して、第２実施形態に示したように、単に各領域内で、輝度値の最も高い画素を抽出しても（図３におけるステップＳ４）、領域内に残存する副次的なピークやノイズ（例えば、図７おける領域Ｄ2’、Ｄ3’内のピーク）を誤検出する可能性はある。そこで、分割された領域から、クラスタリングにより輝度値の平均値が比較的高い領域（例えば、図７おける領域Ｄ1’）のみを抽出し、抽出した領域内においてスポット検出部４０が輝度分布のピーク点を検出することで誤検出を避けることができる。
【００９０】
クラスタリングにおいては、例えば、画像データ空間を粗視化した後、粗視化経験確率分布の算出、各ピクセルのクラス帰属確率の算出、クラスの属性を規定するパラメータの更新及び評価関数の算出を評価関数の変化がなくなるまで繰り返し、閾値を自動的に決定する、特開２００３−３０３３４４号公報に開示された画像抽出方法を使用する。
【００９１】
次に、スポット定量装置１の動作について図１１を参照して、詳細に説明する。
【００９２】
図１１は、第２実施形態にかかるスポット定量方法の処理手順を示すフローチャートである。同図を参照すると、本実施形態にかかるスポット定量方法は、ステップＳ３の後にステップＴ１を実行し、ステップＴ１の後にステップＳ４を実行する他は、第１実施形態にかかるステップ定量方法と同様である。
【００９３】
領域分割処理により画像領域が抽出された後（ステップＳ３）、クラスタリング部６０は、抽出された画像領域の各画素についてクラスタリングを行い、輝度値の比較的高い画像領域を抽出する（ステップＴ１）。
【００９４】
そして、スポット検出部４０は、ステップＴ１で抽出された画像領域の輝度分布のピーク点をスポットの位置として検出する（ステップＳ４）。
【００９５】
上述したように、本実施形態によれば、領域分割された画像領域における各画素を輝度分布に基づいて複数のクラスタに分類し、輝度値の平均値が最も高いクラスタに帰属する画素のうち最高輝度の画素の位置をスポット位置として検出する。このため、領域分割処理で分割した画像領域内に残存するノイズをさらに低減し、スポットの特徴をより正確に定量することができる。
【００９６】
なお、上述の実施形態では、細胞の蛍光染色画像を入力とし、蛍光抗体によって標識された遺伝子領域のスポットを定量する場合を例にとって説明したが、それ以外の染色手法、あるいは撮影方法であっても、細胞核内の染色体などに対応するスポットが局在し、輝度分布によって表示されている画像であれば、本実施形態にかかる方法を使用して解析できる。
【００９７】
それ以外の場合にも、例えば、核、遺伝子(または染色体)がＲＧＢ（Ｒｅｄ、Ｇｒｅｅｎ、Ｂｌｕｅ）の異なるチャネルで撮影されたカラー画像を入力とする場合であっても、原画像をＲＧＢに対応した３枚のグレースケール画像に分解し、解析対象となるスポットを含んだグレースケール画像に対して同様な手順を行えばよい。
【００９８】
フィルタリングで使用するカーネルは、上記（１）〜（５）式で定義したものに限られない。ノイズを低減し、スポットの特徴を定量するのに適したカーネルであれば、上述したもの以外のカーネルを使用することもできる。
【００９９】
カーネルのパラメータを所定範囲内で変更してフィルタリングを複数回施す場合（図３におけるステップＳ５）、上記（１）〜（５）式で定義したカーネルのうち１つのみ使用するのでなく、例えば、ガウシアンカーネル（（１）式）及びメキシカンハットカーネル（（５）式）など、複数の種類のカーネルを使用することもできる。この場合、それぞれのカーネルにおいてパラメータを所定範囲内で変更してフィルタリングを施し、輝度値のピーク値を求め、得られたピーク値のうち、最も高いピーク値に対応するパラメータの値に基づいてスポットの特徴を定量する。
【０１００】
また、図３のステップＳ２〜Ｓ６の処理の全部又は一部は、コンピュータプログラムで実現することもできる。
【図面の簡単な説明】
【０１０１】
【図１】第１実施形態におけるスポット定量装置の構成を示すブロック図である。
【図２】第１実施形態における、スポット定量装置の構成の一例を示した図である。
【図３】第１実施形態における、スポット定量方法の処理手順を示すフローチャートである。
【図４】第１実施形態における、領域分割方法の処理手順を示すフローチャートである。
【図５】第１実施形態における、入力画像の一例を示す図である。
【図６】第１実施形態における、入力画像の輝度分布を３次元表示した一例を示す図である。
【図７】第１実施形態における、領域分割処理により分割された画像領域の一例を示す図である。
【図８】第１実施形態における、同心円カーネルでフィルタリングを施して得られた特徴量の一例である。
【図９】（ａ）第１実施形態における、直径の比較的大きなスポットの輝度分布を二次元表示した一例を示す図である。（ｂ）第１実施形態における、直径の比較的小さなスポットの輝度分布を二次元表示した一例を示す図である。
【図１０】第２実施形態における、スポット定量装置の構成を示すブロック図である。
【図１１】第２実施形態における、スポット定量方法の処理手順を示すフローチャートである。
【符号の説明】
【０１０２】
１，２  スポット定量装置
１０    画像入力部
１５    フィルタリング部
２０    画像フィルタリング部
３０    領域分割部
４０    スポット検出部
５０    特徴量算出部
６０    クラスタリング部
１００    光源
１０１    励起フィルタ
１０２    バリアフィルタ
１０３    ＣＣＤカメラ
１０４    試料
１０５    対眼レンズ
１０６    対物レンズ
１０７    ダイクロックミラー
１０８    コンピュータ
１０９    モニタ
Ｓ１〜Ｓ６，Ａ１〜Ａ１０，Ｔ１    ステップ
Ｐ１，Ｐ２，Ｐ６，Ｐ７    スポット
Ｉ’，Ｄ1’，Ｄ2’，Ｄ3’，Ｄ4’，Ｒ0’ Ｒ1’    画像領域
【図１】
0005347272.tif 000008
【図２】
0005347272.tif 000009
【図３】
0005347272.tif 000010
【図４】
0005347272.tif 000011
【図８】
0005347272.tif 000012
【図９】
0005347272.tif 000013
【図１０】
0005347272.tif 000014
【図１１】
0005347272.tif 000015
【図５】
0005347272.tif 000016
【図６】
0005347272.tif 000017
【図７】
0005347272.tif 000018