JP2004504659A - 基板における所望領域の位置決め法 - Google Patents
基板における所望領域の位置決め法 Download PDFInfo
- Publication number
- JP2004504659A JP2004504659A JP2002512708A JP2002512708A JP2004504659A JP 2004504659 A JP2004504659 A JP 2004504659A JP 2002512708 A JP2002512708 A JP 2002512708A JP 2002512708 A JP2002512708 A JP 2002512708A JP 2004504659 A JP2004504659 A JP 2004504659A
- Authority
- JP
- Japan
- Prior art keywords
- pixel
- intensity
- matrix
- value
- substrate
- 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.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10056—Microscopic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20092—Interactive image processing based on input by user
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30072—Microarray; Biochip, DNA array; Well plate
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
基板内の少なくとも一領域の信号強度であってバックグランド強度を有するものをコンピュータによって評価する方法において、基板の散乱パラメータを決定する。また、領域及びその周囲を囲む評価ウィンドウ内に位置する画素強度のマトリックス(IS)が決定される。マトリックス(IS)の画素強度は画素強度のヒストグラムを得るのに用いる。この評価ヒストグラムは2つの分布ピークを示す。低い強度の第1の分布ピークは周囲画素に対応し、高い強度の第2の分布ピークは領域画素に対応する。2つの分布ピークを有する曲線はこの評価ヒストグラムでフィッティングされる。散乱パラメータを用いて、散乱についてマトリックス(IS)又はフィッティングされた曲線を補正する。領域での画素信号強度は、画素強度ヒストグラムでフィッティングされた曲線から得られたデータによって決定される。
Description
【0001】
【発明の属する技術分野】
本発明は、基板に埋め込まれている基板の少なくとも一の領域の信号強度であってバックグランド強度を有するものをコンピュータによって評価する方法、及び、基板に埋め込まれている基板の関心有力領域を位置決めする方法に関するものである。
【0002】
【従来の技術題】
米国特許出願第5,795,716号明細書は、コンピュータ支援画像技術において使用されるこのタイプの方法及び塩基配列評価の分析装置を開示している。通常、基板は領域のアレイを備えている。ここで、各領域は、分析物と特別に結合することができる、公知の結合物質あるいはプローブを有する。アレイを用いることができる分析(検定)には、交雑(ハイブリダイゼーション)、免疫検査(測定)、レセプター/配位子検査等が含まれてもよい。例えば、アレイは、多数の分析物の存在のために、血液のような生物学的試料を遮蔽するために用いてもよい。基板を一又は二以上の分析物を含む液体に接触させるならば、反応パターンが生成されて、アレイの結合物質について分析物の特定の親和性を示すものでもよい。アレイは核酸プローブを含む領域から成ってもよい。アレイは、ウィルスあるいはバクテリアの核酸の検出及び/又は型別のため、又は突然変異の検出のために用いてもよい。同様に、アレイは免疫検査を実施するのに用いることもできる。その場合には、結合物質あるいはプローブは抗原(ペプチド)あるいは抗体であってもよい。分析物を結合物質に結合することによって、蛍光信号のような検出可能な信号は、例えば、第2の標識(ラベル付き)試薬に結合することによって生成される。例えば、スキャナーはイメージファイルを生成し、このイメージファイルを評価して各領域の信号強度を決定する。正確な情報を得るために、結合物質が配置する場所である当該領域の信号強度を正確に決定することが非常に重要である。公知の方法では、この信号強度は“ブランク”領域のバックグランド強度を用いることによって決定される。ここで、ブランク領域はバックグランド強度だけを与えると仮定されている。しかしながら、確実にバックグランド強度だけを提供する基板における領域を見つけることは非常に困難である。さらに、公知の方法の多くは、この方法で決定された一のバックグランド強度は、基板の表面に沿ったバックグランド強度の変動を全く考慮していない。さらに、基板における散乱に起因する測定信号強度における効果も全く考慮されていない。
【0003】
【発明が解決しようとする課題】
本発明の目的は、上述のタイプの改良型の方法を提供することである。
【0004】
【課題を解決するための手段】
本発明では、基板に埋め込まれている基板の少なくとも一の領域の信号強度であってバックグランド強度を有するものをコンピュータによって評価する方法において、基板の散乱パラメータを決定し、領域及びその周囲部を囲む評価ウィンドウ内に位置する画素の画素強度のマトリックス(IS)を決定し、ここで、マトリックス(IS)の画素強度は画素強度の評価ヒストグラムを得るのに用い、前記評価ヒストグラムは2つの分布ピークを示し、低い強度の第1の分布ピークは周囲画素に対応し、高い強度の第2の分布ピークは領域画素に対応し、2つの分布ピークを有する曲線は前記評価ヒストグラムでフィッティングされ、散乱パラメータを用いて、散乱についてマトリックス(IS)又はフィッティングされた曲線を補正し、領域での画素信号強度は前記評価ヒストグラムでフィッティングされた曲線から得られたデータによって決定されることを特徴とする。
【0005】
このように、基板の散乱特性について、領域とその周囲の両方の画素強度を補正することによって、向上した精度で領域における画素の平均信号強度を決定する方法が得られる。本発明の方法では、バックグランド強度の評価は当該領域の信号強度の評価の本質的な一部である。それによって、バックグランド強度を決定するために基板の一部を位置決めする独立の段階は必要とされず、さらに、基板の表面に沿ったバックグランド強度の変動は評価結果には影響しないように、各領域のバックグランド強度が対応する領域に近接する位置で決定される。
【0006】
本発明はさらに、基板周囲に埋め込まれている基板の所望の有力領域をコンピュータによって位置決めする方法において、基板のイメージファイルをローパスフィルタアルゴリズムで処理してイメージファイルの全画素について局所的平均値のマトリックスを決定し、局所的平均値のマトリックスはイメージファイルの実際の画素値のマトリックスと結合して、実際の画素値に依存する高又は低のいずれかが対応する局所的平均値より上又は低いのいずれかである高/低画素値の第1のマトリックスを得、第1のマトリックスの高/低画素値をさらに中央フィルタアルゴリズムで処理して中央画素値の第2のマトリックスを得、各中央画素値は第1のマトリックスの高/低画素値の多数と等しく、中央画素値を列毎及び行毎に処理して平均列値及び平均行値をそれぞれ決定し、最大平均列値及び最大平均行値は所望の有力領域の列中心線と行中心線の評価として選択し、ここで、それら列中心線と行中心線の交点が所望の有力領域の中心であることを特徴とする。
【0007】
【発明の実施の形態】
基板周囲に埋め込まれている基板の領域の信号強度を評価する本発明の方法は、コンピュータで使用するコンピュータプログラムとして使用される。この方法は、例えば、交雑された核酸プローブの信号強度を分析することによる塩基配列評価のために用いることができる。しかしながら、上述のように、基板は、免疫検査(測定)、レセプター/配位子検査等の他の分析について用いることができる。本発明の方法は信号強度の評価だけに関するので、塩基配列評価の方法は詳述しない。例えば、国際公開第9902266号パンフレットを参照されたい。
【0008】
このような塩基配列評価においては、国際公開第9902266号パンフレットに記載されているような、例えば、図3で示しような場合において、各領域あるいはドットが公知の塩基配列を有するプローブを有する6×4領域1のアレイを有する基板が使用できる。蛍光ラベル(標識)を有する材料が基板の表面に接触するとき、標識材料の異なる濃度が異なるドット領域で見つけられる。基板の照明は、各ドット領域1に起因する蛍光を与え、図3はCCDカメラによって得られたイメージを示す。この場合には蛍光が得られるが、本発明の方法は蛍光強度を評価することには限定されない。本発明の方法は、周囲に埋め込まれている領域からの信号強度の種類を評価することを含む。図3で示したイメージをコンピュータ、例えば、ハードディスクに、イメージファイルとして格納する。図3からわかるように、周囲からドット領域1への鋭い遷移(トランジッション)はなく、そのため、ドット領域1の近傍におけるバックグランド強度を決定するために位置を見つけることは困難である。ドット領域から周囲への漸進的なあるいはぼやけた遷移は標識材料の蛍光の散乱に起因する。この散乱効果を補正するために、基板の散乱パラメータを本発明の方法で決定する。
【0009】
図1は例えば、基板の散乱パラメータを決定する方法を示している。図1は、この例では、60μmの層厚を有する基板を側方から照射したイメージである。図1では、ライン(線)2は画素強度が決定される線に沿って描かれている。“画素強度”の語は、CCDカメラあるいはアレイのイメージを提供するのに用いられる他の種類のイメージング装置の一の画素から得られた信号強度を示すために用いられることに留意されたい。このライン2に沿って測定された画素強度から、図2に示したような、この場合は33.8μの特徴的な減衰を有する指数散乱関数を作ることができる。この特徴的な散乱減衰は基板の散乱効果について測定された画素強度を補正するために本発明の方法で用いる。この方法で決定された散乱減衰を基板の側方照射によって測定するときは、通常の測定環境において使用されるイメージング装置に直接向かう散乱されない照射は考慮されない。散乱なしで直接基板に起因するこのような信号強度を補正するために、例えば、全画素について受ける信号強度の予め決めた割合を差し引くことができる。
【0010】
図1では乾燥(ドライ)基板を測定される点に留意されたい。ウェット状態についての同じ基板の散乱パラメータは対応する方法で決定することができる。異なる種類の基板についてさらなる散乱パラメータを決定してコンピュータに保存することができる。評価法では、格納された散乱パラメータを、使用される基板の種類及び状態に応じてコンピュータのメモリから選択することができる。適当な入力装置を介して使用される基板の公知の散乱パラメータを入れることも可能である。この点では、基板の散乱パラメータを決定する段階が、記載された方法での使用のために前に決定した散乱パラメータの入力の方法を包含することに留意されたい。
【0011】
ドットの領域1はドット毎に評価される。この評価のために、一のドット領域1は、図5で示した評価(推定)ウィンドウ3によって基板イメージの残りの部分から、そのすぐ周囲部と共に孤立されている。評価ウィンドウ3のサイズは、全ドット領域1とそのすぐの周囲部とが評価ウィンドウ3内に配置するようになっている。図5の例では、窓3は55×55画素を備える。ドット領域1を孤立させる際に、まずコンピュータプログラムは図3の全イメージを評価し、各ドット領域1の中心に配置し、使用される評価ウィンドウ3のサイズを決定する。適当なユーザーインターフェースによって、各ドット領域の中心位置及び/又はウィンドウのサイズを変更することができる。ドットの中心を配置するために、図3の全イメージを評価する好適な方法を以下に記載する。
【0012】
評価ウィンドウ3内の全画素の強度は画素強度のマトリックスISにおいて挿入する。すなわち、この強度は基板の散乱効果を含む信号強度を表す。これらの画素強度はこの記載では散乱画素強度と称する。基板の散乱パラメータが公知であるとき、散乱効果を数学的に補正することができる。記載した方法の第1の実施形態では、散乱強度をデコンボルーション(逆重畳)によって考慮される。このようなデコンボルーション法は数学においては公知である。基板の散乱パラメータで散乱画素強度のマトリックスISを逆重畳した結果は、非散乱画素強度マトリックスIDである。すなわち、これらの非散乱画素強度マトリックスは散乱効果について補正された信号強度を表す。図4は、デコンボルーション後、図3のイメージファイルのコンピュータモニターにおける表示を示している。この処理されたイメージファイルの表示においては、ドット領域及びその直ぐ周囲での変わり目での散乱のぼやけ効果が除去されることがわかる。
【0013】
この実施形態の次の段階では、図6に示した画素強度評価ヒストグラムを決定するためにマトリックスIDの非散乱画素強度を用いる。明らかに、画素強度のこのヒストグラムは2つの分布ピーク4及び5を示し、ここでは、低い強度を有する第1の分布ピーク4は周囲部の画素を表し、高い強度を有する第2の分布ピーク5はドット領域1の画素を表している。
【0014】
これら両方の分布4及び5は通常標準分布と似ているので、標準のフィッティング法を用いて2つの通常の分布曲線6及び7を分布ピーク4及び5にフィットさせることができる。このようなフィッティング法はこれ自体公知であるので、詳細な説明を必要がないだろう。例えば、最小二乗平均法を用いて、最適なフィッティング曲線を見つけることができる。通常の分布曲線6及び7の決定に際して、マトリックス(IS)の画素強度にあるいは図3の基板の全イメージファイルの全画素強度に存在するノイズのノイズパラメータ、特に、このノイズの標準偏差を用いる。この標準偏差σnoiseは例えば、平均値0でかつ標準偏差σ=σnoise*√2の分布となる隣接画素間の差をとることによって決定される。
【0015】
最終的に、第1のフィットされた分布曲線6の平均値は、純粋なバックグランド強度の平均画素強度についての最良の評価としてとる。第2のフィットされた分布曲線7の平均値は、ドット領域1の平均画素強度、すなわち、標識材料の蛍光強度と合わせたバックグランド強度についての最良の評価としてとる。従って、両方の平均値の差が、ドット領域1すなわち標識材料における画素の平均蛍光強度を与える。
【0016】
本発明の方法の第2の実施形態では、基板の散乱パラメータが異なる方法で考慮される。デコンボルーションはCCDカメラから得られたデータに存在するノイズを増幅することが公知である。低強度信号の場合、デコンボルーションによる評価が低質の結果を与える。
【0017】
本発明の第2の実施形態では、散乱パラメータは、図7Aから図7Cで模式的に示したような予め決められた理論による標識領域の散乱応答を決定することによって考慮される。図7Aでは、異なる散乱パラメータ、すなわち、異なる基板及び異なる基板条件についての散乱パラメータについては、異なる散乱応答が示されている。X=1での直線は非散乱応答曲線に、すなわち、理論的な標識領域内の強度Iでその領域の周囲部での強度0であることに対応する。この非散乱応答を異なる散乱パラメータに結合すると、図7で示した異なる散乱応答になる。図7Aでは、規格化された標識領域及び規格化された強度について、すなわち、領域の半径Rと強度Iとを1として、散乱応答を示している。図7Bは、対応する散乱応答曲線8と共に理論的標識領域の散乱応答を示す図5において使用される評価ウィンドウと同じサイズの評価ウィンドウを示している。この散乱応答曲線8は、評価ウィンドウの端まで理論的領域の中心から一の半径に沿って強度を示すので、図7Cで示したような画素強度の規格化されたヒストグラムは、散乱応答曲線8と標識領域の中心を囲繞する画素の数とによって与えられた画素強度値を用いることによって作成することができる。
【0018】
図8は、本発明の方法の第2の実施形態における次の2段階を示す。図8Aによれば、図7Cの規格化されたヒストグラムは、特定のバックグランド強度IBと領域強度IAとを選択することによって倍率が変えられる。この倍率が変えられたヒストグラムは、上述の方法で得られたノイズ分布パラメータでたたみ込む。この方法では、理論散乱ヒストグラムを得て、この理論散乱ヒストグラムは散乱で生じたバックグランド強度と領域強度とを表す2つの分布ピーク9及び10を示す。理論散乱ヒストグラムは、マトリックスISの散乱された画素強度から得られた評価ヒストグラム上でフィッティングされる。標準的なフィッティング法を用いて、少なくともバックグランド強度IB及び領域強度IAの値を変えることによって最良のフィッティング理論散乱ヒストグラムを見つけることができる。他の変形として、理論領域の半径、ノイズの標準偏差及び使用される散乱パラメータ値を変えることも可能である。
【0019】
図9は、最良のフィッティング散乱ヒストグラムを有するマトリックスISの散乱画素強度から得られた評価ヒストグラムを示す。最良のフィッティング散乱ヒストグラム曲線を見つけるとき、対応するバックグランド強度と領域画素強度の値を、バックグランド強度及び領域画素強度についての最良の評価として用いる。領域画素強度とバックグランド強度との差はドット領域1における画素の蛍光強度に対応する。
【0020】
記載された本方法の上述の例では、フィッティングされた曲線のピークの平均値だけを評価に用いられることに留意されたい。フィッティングされた曲線はさらに、検査用に用いることができる他の情報も含む。例えば、各ピークの表面領域を、ドット領域及び周囲バックグランドにおける画素数を決定するのに用いることができる。
【0021】
図10Aから図10Dと図11Aから図11Fとは、ドット領域1の中心を配置し、上述の方法において使用される評価ウィンドウ3のサイズを決定するために、図3に示されたような基板のイメージファイルを評価する方法の好適な実施形態における段階を示している。全イメージファイルを評価するこの方法は、一ドット領域を評価する方法と組み合わせて使用することができ、従って、この方法は上述の方法や他の特定のドット領域評価方法には限定されないことに留意されたい。
【0022】
図10Aは、図3の基板に対応する基板のイメージファイルの表示を示している。このイメージファイルは明らかにノイズを有するが、当該有力(可能な、予定の)領域のいくつかはイメージファイルの画素の未処理画素値の表示におけるドット領域1として認識することができる。当該有力領域のすべてを配置するために、イメージファイルの画素値をローパスフィルタアルゴリズムにおいて処理してイメージファイルの全画素について局所的な平均値のマトリックスを決定する。イメージファイルの全画素について局所的な平均値の得られたマトリックスを図10Bに示した。これらの局所的平均値は、ドット領域の信号強度を評価するために用いられる評価ウィンドウの好適には25%から50%のサイズを有する図示しない第2のウィンドウを用いることによって得られる。この第2のウィンドウは、画素毎にイメージファイルの画素に沿ってずらす。第2のウィンドウの各位置では、第2のウィンドウ内のイメージファイルの全画素の平均値を、第2のウィンドウの中心での画素の局所的平均値としてとる。第2のウィンドウをイメージファイルの全面に沿ってずらすことによって、イメージファイルの各画素について局所的平均値を決定する。このようにして、イメージファイルの実際の全高周波情報を除去する。
【0023】
イメージファイルの全画素について高/低画素値の第1のマトリックスを得るために、イメージファイルの実際の画素値のマトリックスを上述の方法で得られた局所的平均値と組み合わせる。これらの高/低画素値は、実際の画素値が画素の対応する局所的平均値より上か下かに依存して高いか低い。これは、実際の画素値が対応する局所的平均値より上であれば、第1のマトリックスの高/低画素値は高く、他方、実際の画素値が対応する局所的平均値より下であれば、第1のマトリックスの高/低画素値は低いことを意味する。得られた第1のマトリックスを図10Cに示す。
【0024】
第1のマトリックスの高/低画素値を中央(メジアン)域フィルタアルゴリズムで処理して、イメージファイルの全画素について中央画素値の第2のマトリックスを得る。中央画素値の得られたマトリックスを図10Dに示す。この中央域フィルタアルゴリズムによって、第1のマトリックスの各画素の高/低画素値をこの画素をすぐ囲む画素の高/低画素値と比較する。高/低画素値は、これらの周囲画素の高/低画素値の多数の値と等しく作られ、この新しい値を中央画素値と称する。
【0025】
図10Dで示すように、得られたマトリックスは当該可能な領域すべてを明瞭に示している。これらの可能な当該領域の中心の配置についての第1の粗い評価は、列毎及び行毎に平均列値及び平均行値を決定することによって得られる。これらの平均列値及び平均行値は、列毎及び行毎に、第2のマトリックスの全中央画素値を加えること、及び、得られた値を画素数で割ることによって得られる。これらの平均値は、図10Dのイメージファイルの表示の下部及ぶ右側で曲線10及び11として示されている。6×4個の有力領域を有する基板を用いることは公知であるので、4個の最高平均列値と6個の最高平均行値とが、列中心線12と行中心線13の評価として選択される。これらの列及び行中心線12,13の交差点は、全有力領域の中心14の位置として用いられる。
【0026】
図11Aから図11Fはさらに、当該有力領域の中心を配置する方法の段階を示すものである。図11Aは、列中心線12と行中心線13とを有する当該領域の評価ウィンドウを示す。さらなる段階では、線11及び12の交差点14が真に当該領域の画素の重心15にあるか否かを決定する。この結果、重心は、図11Aで示したような評価ウィンドウ内の第2のマトリックスの画素の中央画素値から決定される。この重心が交差点14に一致しないならば、中心線12,13は、交差点が図11Bで示したような重心15に一致するようにずらされる。
【0027】
高い値を有する評価ウィンドウ内の画素が真に当該領域を表すか又はただのランダムノイズであるか、すなわち、基板のこの位置での結合実体は分析物に結合しなかったか否かを決定するために、以下の段階を実施する。図11C及び図11Eは、評価ウィンドウ内の2つの有力領域を示す。第1に、第2のマトリックスの評価ウィンドウ内の高い値を有する画素の表面を決定し、この表面から当該領域の半径を決定する。さらに、この表面の半径をこの表面の周縁から決定する。図11Cに示したように、高い値を有する画素は正確な円を表さず、そのため、この領域の周縁から決定された半径は図11Cで示した正確な円16の半径より大きい。しかしながら、2つの半径の比は、所定の参照値例えば、2より小さいべきである。図11Dは、高い値を有する画素の表面の周縁での高い値を有する全画素を示す。周縁から得られた半径と表面から得られた半径との比はこの場合、1.47である。これは、当該領域が存在することを意味する。図11Fでは、図11Eで示した表面の周縁での高い値を有する全画素を示し、周縁から得られた“半径”と図11Eで示した表面から得られた半径との比はこの場合、4.47である。これは、図11Eの場合における評価ウィンドウが当該領域を囲繞しない、すなわち、当該領域が存在しないことを意味する。
【0028】
当該領域の存在又は欠落を決定するために、代替の段階あるいはさらなる段階として、当該可能な領域の重心を、図11で示したような評価ウィンドウ内の高い値を有する全画素の表面に対応する表面を有する虚の(想像上の)円ウィンドウの中心としてとることも可能である。その後、この虚の円ウィンドウは、ある位置を見つけるために重心に対して移動でき、ここで、虚の円ウィンドウは、例えば、全画素の少なくとも90%で高い値を有する所定数の画素をカバーする。このような位置を見つけることができないならば、当該領域は存在しない。このような位置を見つけることができるならば、当該領域が存在する。
【0029】
当該可能な領域を調べるため、かつ、当該領域の信号強度を評価するためにのに用いられるサイズは、高い値を有する画素の数が低い数の画素の数に一致するように決定される。これは、ドット領域1の存在する場合では、ドット画素の数とバックグランド画素の数とが等しい、すなわち、全画素数の50%であることを意味する。評価ウィンドウ内にノイズだけが存在するならば、このような分布も存在し、すなわち、評価ウィンドウ内の画素の50%が平均値より上であり、50%が平均値より下である。このような評価ウィンドウは、ドット領域の信号強度を評価する上述の方法において高い信頼性でフィッティングを行うことができるという利点を示す。また、ユーザーがドット領域の中心を移動したいあるいは評価ウィンドウのサイズを変更したいならば、ユーザーは適当なユーザーインターフェースによってそうすることができる。
【0030】
本発明は、上述の実施形態に限定されず、特許請求の範囲内で多くの変更が可能である。
【図面の簡単な説明】
【図1】基板の散乱パラメータを決定するために用いる画像を示す説明図である。
【図2】図1の画像から得られた散乱減衰関数を示す図である。
【図3】各領域が公知のプローブを有する6×4領域のアレイを有する基板であって、基板上に標識材料の堆積後のイメージを示す図である。
【図4】散乱効果を機械的に除去した後の図3のイメージを示す図である。
【図5】図4の領域のうちの一つ上の評価ウィンドウを示す図である。
【図6】2つのフィッティングされた曲線を有する
【図7】本発明の方法の第2の実施形態においてフィッティングのために用いられる理論的散乱ヒストグラム決定を示す図である。
【図8】本発明の方法の第2の実施形態においてフィッティングのために用いられる理論的散乱ヒストグラム決定を示す図である。
【図9】図7及び図8で示された方法で得られるフィッティングされる理論的ヒストグラムと共に、図3の領域のマトリックス画素強度ISから得られた評価ヒストグラムを示す図である。
【図10】(A)図3で示した基板のイメージファイルを評価するために本発明に係る方法の好適な一実施形態の一段階の説明図である。(B)他の段階の説明図である。(C)他の段階の説明図である。(D)他の段階の説明図である。
【図11】(A)図3で示した基板のイメージファイルを評価するために本発明に係る方法の好適な一実施形態の一段階の説明図である。(B)他の段階の説明図である。(C)他の段階の説明図である。(D)他の段階の説明図である。(E)他の段階の説明図である。(F)他の段階の説明図である。
【符号の説明】
1 ドット領域
2 ライン
3 評価ウィンドウ
4 第1の分布ピーク
5 第2の分布ピーク
6,7 通常の分布曲線
8 散乱応答曲線
12 列中心線
13 行中心線
14 領域の中心
15 重心
【発明の属する技術分野】
本発明は、基板に埋め込まれている基板の少なくとも一の領域の信号強度であってバックグランド強度を有するものをコンピュータによって評価する方法、及び、基板に埋め込まれている基板の関心有力領域を位置決めする方法に関するものである。
【0002】
【従来の技術題】
米国特許出願第5,795,716号明細書は、コンピュータ支援画像技術において使用されるこのタイプの方法及び塩基配列評価の分析装置を開示している。通常、基板は領域のアレイを備えている。ここで、各領域は、分析物と特別に結合することができる、公知の結合物質あるいはプローブを有する。アレイを用いることができる分析(検定)には、交雑(ハイブリダイゼーション)、免疫検査(測定)、レセプター/配位子検査等が含まれてもよい。例えば、アレイは、多数の分析物の存在のために、血液のような生物学的試料を遮蔽するために用いてもよい。基板を一又は二以上の分析物を含む液体に接触させるならば、反応パターンが生成されて、アレイの結合物質について分析物の特定の親和性を示すものでもよい。アレイは核酸プローブを含む領域から成ってもよい。アレイは、ウィルスあるいはバクテリアの核酸の検出及び/又は型別のため、又は突然変異の検出のために用いてもよい。同様に、アレイは免疫検査を実施するのに用いることもできる。その場合には、結合物質あるいはプローブは抗原(ペプチド)あるいは抗体であってもよい。分析物を結合物質に結合することによって、蛍光信号のような検出可能な信号は、例えば、第2の標識(ラベル付き)試薬に結合することによって生成される。例えば、スキャナーはイメージファイルを生成し、このイメージファイルを評価して各領域の信号強度を決定する。正確な情報を得るために、結合物質が配置する場所である当該領域の信号強度を正確に決定することが非常に重要である。公知の方法では、この信号強度は“ブランク”領域のバックグランド強度を用いることによって決定される。ここで、ブランク領域はバックグランド強度だけを与えると仮定されている。しかしながら、確実にバックグランド強度だけを提供する基板における領域を見つけることは非常に困難である。さらに、公知の方法の多くは、この方法で決定された一のバックグランド強度は、基板の表面に沿ったバックグランド強度の変動を全く考慮していない。さらに、基板における散乱に起因する測定信号強度における効果も全く考慮されていない。
【0003】
【発明が解決しようとする課題】
本発明の目的は、上述のタイプの改良型の方法を提供することである。
【0004】
【課題を解決するための手段】
本発明では、基板に埋め込まれている基板の少なくとも一の領域の信号強度であってバックグランド強度を有するものをコンピュータによって評価する方法において、基板の散乱パラメータを決定し、領域及びその周囲部を囲む評価ウィンドウ内に位置する画素の画素強度のマトリックス(IS)を決定し、ここで、マトリックス(IS)の画素強度は画素強度の評価ヒストグラムを得るのに用い、前記評価ヒストグラムは2つの分布ピークを示し、低い強度の第1の分布ピークは周囲画素に対応し、高い強度の第2の分布ピークは領域画素に対応し、2つの分布ピークを有する曲線は前記評価ヒストグラムでフィッティングされ、散乱パラメータを用いて、散乱についてマトリックス(IS)又はフィッティングされた曲線を補正し、領域での画素信号強度は前記評価ヒストグラムでフィッティングされた曲線から得られたデータによって決定されることを特徴とする。
【0005】
このように、基板の散乱特性について、領域とその周囲の両方の画素強度を補正することによって、向上した精度で領域における画素の平均信号強度を決定する方法が得られる。本発明の方法では、バックグランド強度の評価は当該領域の信号強度の評価の本質的な一部である。それによって、バックグランド強度を決定するために基板の一部を位置決めする独立の段階は必要とされず、さらに、基板の表面に沿ったバックグランド強度の変動は評価結果には影響しないように、各領域のバックグランド強度が対応する領域に近接する位置で決定される。
【0006】
本発明はさらに、基板周囲に埋め込まれている基板の所望の有力領域をコンピュータによって位置決めする方法において、基板のイメージファイルをローパスフィルタアルゴリズムで処理してイメージファイルの全画素について局所的平均値のマトリックスを決定し、局所的平均値のマトリックスはイメージファイルの実際の画素値のマトリックスと結合して、実際の画素値に依存する高又は低のいずれかが対応する局所的平均値より上又は低いのいずれかである高/低画素値の第1のマトリックスを得、第1のマトリックスの高/低画素値をさらに中央フィルタアルゴリズムで処理して中央画素値の第2のマトリックスを得、各中央画素値は第1のマトリックスの高/低画素値の多数と等しく、中央画素値を列毎及び行毎に処理して平均列値及び平均行値をそれぞれ決定し、最大平均列値及び最大平均行値は所望の有力領域の列中心線と行中心線の評価として選択し、ここで、それら列中心線と行中心線の交点が所望の有力領域の中心であることを特徴とする。
【0007】
【発明の実施の形態】
基板周囲に埋め込まれている基板の領域の信号強度を評価する本発明の方法は、コンピュータで使用するコンピュータプログラムとして使用される。この方法は、例えば、交雑された核酸プローブの信号強度を分析することによる塩基配列評価のために用いることができる。しかしながら、上述のように、基板は、免疫検査(測定)、レセプター/配位子検査等の他の分析について用いることができる。本発明の方法は信号強度の評価だけに関するので、塩基配列評価の方法は詳述しない。例えば、国際公開第9902266号パンフレットを参照されたい。
【0008】
このような塩基配列評価においては、国際公開第9902266号パンフレットに記載されているような、例えば、図3で示しような場合において、各領域あるいはドットが公知の塩基配列を有するプローブを有する6×4領域1のアレイを有する基板が使用できる。蛍光ラベル(標識)を有する材料が基板の表面に接触するとき、標識材料の異なる濃度が異なるドット領域で見つけられる。基板の照明は、各ドット領域1に起因する蛍光を与え、図3はCCDカメラによって得られたイメージを示す。この場合には蛍光が得られるが、本発明の方法は蛍光強度を評価することには限定されない。本発明の方法は、周囲に埋め込まれている領域からの信号強度の種類を評価することを含む。図3で示したイメージをコンピュータ、例えば、ハードディスクに、イメージファイルとして格納する。図3からわかるように、周囲からドット領域1への鋭い遷移(トランジッション)はなく、そのため、ドット領域1の近傍におけるバックグランド強度を決定するために位置を見つけることは困難である。ドット領域から周囲への漸進的なあるいはぼやけた遷移は標識材料の蛍光の散乱に起因する。この散乱効果を補正するために、基板の散乱パラメータを本発明の方法で決定する。
【0009】
図1は例えば、基板の散乱パラメータを決定する方法を示している。図1は、この例では、60μmの層厚を有する基板を側方から照射したイメージである。図1では、ライン(線)2は画素強度が決定される線に沿って描かれている。“画素強度”の語は、CCDカメラあるいはアレイのイメージを提供するのに用いられる他の種類のイメージング装置の一の画素から得られた信号強度を示すために用いられることに留意されたい。このライン2に沿って測定された画素強度から、図2に示したような、この場合は33.8μの特徴的な減衰を有する指数散乱関数を作ることができる。この特徴的な散乱減衰は基板の散乱効果について測定された画素強度を補正するために本発明の方法で用いる。この方法で決定された散乱減衰を基板の側方照射によって測定するときは、通常の測定環境において使用されるイメージング装置に直接向かう散乱されない照射は考慮されない。散乱なしで直接基板に起因するこのような信号強度を補正するために、例えば、全画素について受ける信号強度の予め決めた割合を差し引くことができる。
【0010】
図1では乾燥(ドライ)基板を測定される点に留意されたい。ウェット状態についての同じ基板の散乱パラメータは対応する方法で決定することができる。異なる種類の基板についてさらなる散乱パラメータを決定してコンピュータに保存することができる。評価法では、格納された散乱パラメータを、使用される基板の種類及び状態に応じてコンピュータのメモリから選択することができる。適当な入力装置を介して使用される基板の公知の散乱パラメータを入れることも可能である。この点では、基板の散乱パラメータを決定する段階が、記載された方法での使用のために前に決定した散乱パラメータの入力の方法を包含することに留意されたい。
【0011】
ドットの領域1はドット毎に評価される。この評価のために、一のドット領域1は、図5で示した評価(推定)ウィンドウ3によって基板イメージの残りの部分から、そのすぐ周囲部と共に孤立されている。評価ウィンドウ3のサイズは、全ドット領域1とそのすぐの周囲部とが評価ウィンドウ3内に配置するようになっている。図5の例では、窓3は55×55画素を備える。ドット領域1を孤立させる際に、まずコンピュータプログラムは図3の全イメージを評価し、各ドット領域1の中心に配置し、使用される評価ウィンドウ3のサイズを決定する。適当なユーザーインターフェースによって、各ドット領域の中心位置及び/又はウィンドウのサイズを変更することができる。ドットの中心を配置するために、図3の全イメージを評価する好適な方法を以下に記載する。
【0012】
評価ウィンドウ3内の全画素の強度は画素強度のマトリックスISにおいて挿入する。すなわち、この強度は基板の散乱効果を含む信号強度を表す。これらの画素強度はこの記載では散乱画素強度と称する。基板の散乱パラメータが公知であるとき、散乱効果を数学的に補正することができる。記載した方法の第1の実施形態では、散乱強度をデコンボルーション(逆重畳)によって考慮される。このようなデコンボルーション法は数学においては公知である。基板の散乱パラメータで散乱画素強度のマトリックスISを逆重畳した結果は、非散乱画素強度マトリックスIDである。すなわち、これらの非散乱画素強度マトリックスは散乱効果について補正された信号強度を表す。図4は、デコンボルーション後、図3のイメージファイルのコンピュータモニターにおける表示を示している。この処理されたイメージファイルの表示においては、ドット領域及びその直ぐ周囲での変わり目での散乱のぼやけ効果が除去されることがわかる。
【0013】
この実施形態の次の段階では、図6に示した画素強度評価ヒストグラムを決定するためにマトリックスIDの非散乱画素強度を用いる。明らかに、画素強度のこのヒストグラムは2つの分布ピーク4及び5を示し、ここでは、低い強度を有する第1の分布ピーク4は周囲部の画素を表し、高い強度を有する第2の分布ピーク5はドット領域1の画素を表している。
【0014】
これら両方の分布4及び5は通常標準分布と似ているので、標準のフィッティング法を用いて2つの通常の分布曲線6及び7を分布ピーク4及び5にフィットさせることができる。このようなフィッティング法はこれ自体公知であるので、詳細な説明を必要がないだろう。例えば、最小二乗平均法を用いて、最適なフィッティング曲線を見つけることができる。通常の分布曲線6及び7の決定に際して、マトリックス(IS)の画素強度にあるいは図3の基板の全イメージファイルの全画素強度に存在するノイズのノイズパラメータ、特に、このノイズの標準偏差を用いる。この標準偏差σnoiseは例えば、平均値0でかつ標準偏差σ=σnoise*√2の分布となる隣接画素間の差をとることによって決定される。
【0015】
最終的に、第1のフィットされた分布曲線6の平均値は、純粋なバックグランド強度の平均画素強度についての最良の評価としてとる。第2のフィットされた分布曲線7の平均値は、ドット領域1の平均画素強度、すなわち、標識材料の蛍光強度と合わせたバックグランド強度についての最良の評価としてとる。従って、両方の平均値の差が、ドット領域1すなわち標識材料における画素の平均蛍光強度を与える。
【0016】
本発明の方法の第2の実施形態では、基板の散乱パラメータが異なる方法で考慮される。デコンボルーションはCCDカメラから得られたデータに存在するノイズを増幅することが公知である。低強度信号の場合、デコンボルーションによる評価が低質の結果を与える。
【0017】
本発明の第2の実施形態では、散乱パラメータは、図7Aから図7Cで模式的に示したような予め決められた理論による標識領域の散乱応答を決定することによって考慮される。図7Aでは、異なる散乱パラメータ、すなわち、異なる基板及び異なる基板条件についての散乱パラメータについては、異なる散乱応答が示されている。X=1での直線は非散乱応答曲線に、すなわち、理論的な標識領域内の強度Iでその領域の周囲部での強度0であることに対応する。この非散乱応答を異なる散乱パラメータに結合すると、図7で示した異なる散乱応答になる。図7Aでは、規格化された標識領域及び規格化された強度について、すなわち、領域の半径Rと強度Iとを1として、散乱応答を示している。図7Bは、対応する散乱応答曲線8と共に理論的標識領域の散乱応答を示す図5において使用される評価ウィンドウと同じサイズの評価ウィンドウを示している。この散乱応答曲線8は、評価ウィンドウの端まで理論的領域の中心から一の半径に沿って強度を示すので、図7Cで示したような画素強度の規格化されたヒストグラムは、散乱応答曲線8と標識領域の中心を囲繞する画素の数とによって与えられた画素強度値を用いることによって作成することができる。
【0018】
図8は、本発明の方法の第2の実施形態における次の2段階を示す。図8Aによれば、図7Cの規格化されたヒストグラムは、特定のバックグランド強度IBと領域強度IAとを選択することによって倍率が変えられる。この倍率が変えられたヒストグラムは、上述の方法で得られたノイズ分布パラメータでたたみ込む。この方法では、理論散乱ヒストグラムを得て、この理論散乱ヒストグラムは散乱で生じたバックグランド強度と領域強度とを表す2つの分布ピーク9及び10を示す。理論散乱ヒストグラムは、マトリックスISの散乱された画素強度から得られた評価ヒストグラム上でフィッティングされる。標準的なフィッティング法を用いて、少なくともバックグランド強度IB及び領域強度IAの値を変えることによって最良のフィッティング理論散乱ヒストグラムを見つけることができる。他の変形として、理論領域の半径、ノイズの標準偏差及び使用される散乱パラメータ値を変えることも可能である。
【0019】
図9は、最良のフィッティング散乱ヒストグラムを有するマトリックスISの散乱画素強度から得られた評価ヒストグラムを示す。最良のフィッティング散乱ヒストグラム曲線を見つけるとき、対応するバックグランド強度と領域画素強度の値を、バックグランド強度及び領域画素強度についての最良の評価として用いる。領域画素強度とバックグランド強度との差はドット領域1における画素の蛍光強度に対応する。
【0020】
記載された本方法の上述の例では、フィッティングされた曲線のピークの平均値だけを評価に用いられることに留意されたい。フィッティングされた曲線はさらに、検査用に用いることができる他の情報も含む。例えば、各ピークの表面領域を、ドット領域及び周囲バックグランドにおける画素数を決定するのに用いることができる。
【0021】
図10Aから図10Dと図11Aから図11Fとは、ドット領域1の中心を配置し、上述の方法において使用される評価ウィンドウ3のサイズを決定するために、図3に示されたような基板のイメージファイルを評価する方法の好適な実施形態における段階を示している。全イメージファイルを評価するこの方法は、一ドット領域を評価する方法と組み合わせて使用することができ、従って、この方法は上述の方法や他の特定のドット領域評価方法には限定されないことに留意されたい。
【0022】
図10Aは、図3の基板に対応する基板のイメージファイルの表示を示している。このイメージファイルは明らかにノイズを有するが、当該有力(可能な、予定の)領域のいくつかはイメージファイルの画素の未処理画素値の表示におけるドット領域1として認識することができる。当該有力領域のすべてを配置するために、イメージファイルの画素値をローパスフィルタアルゴリズムにおいて処理してイメージファイルの全画素について局所的な平均値のマトリックスを決定する。イメージファイルの全画素について局所的な平均値の得られたマトリックスを図10Bに示した。これらの局所的平均値は、ドット領域の信号強度を評価するために用いられる評価ウィンドウの好適には25%から50%のサイズを有する図示しない第2のウィンドウを用いることによって得られる。この第2のウィンドウは、画素毎にイメージファイルの画素に沿ってずらす。第2のウィンドウの各位置では、第2のウィンドウ内のイメージファイルの全画素の平均値を、第2のウィンドウの中心での画素の局所的平均値としてとる。第2のウィンドウをイメージファイルの全面に沿ってずらすことによって、イメージファイルの各画素について局所的平均値を決定する。このようにして、イメージファイルの実際の全高周波情報を除去する。
【0023】
イメージファイルの全画素について高/低画素値の第1のマトリックスを得るために、イメージファイルの実際の画素値のマトリックスを上述の方法で得られた局所的平均値と組み合わせる。これらの高/低画素値は、実際の画素値が画素の対応する局所的平均値より上か下かに依存して高いか低い。これは、実際の画素値が対応する局所的平均値より上であれば、第1のマトリックスの高/低画素値は高く、他方、実際の画素値が対応する局所的平均値より下であれば、第1のマトリックスの高/低画素値は低いことを意味する。得られた第1のマトリックスを図10Cに示す。
【0024】
第1のマトリックスの高/低画素値を中央(メジアン)域フィルタアルゴリズムで処理して、イメージファイルの全画素について中央画素値の第2のマトリックスを得る。中央画素値の得られたマトリックスを図10Dに示す。この中央域フィルタアルゴリズムによって、第1のマトリックスの各画素の高/低画素値をこの画素をすぐ囲む画素の高/低画素値と比較する。高/低画素値は、これらの周囲画素の高/低画素値の多数の値と等しく作られ、この新しい値を中央画素値と称する。
【0025】
図10Dで示すように、得られたマトリックスは当該可能な領域すべてを明瞭に示している。これらの可能な当該領域の中心の配置についての第1の粗い評価は、列毎及び行毎に平均列値及び平均行値を決定することによって得られる。これらの平均列値及び平均行値は、列毎及び行毎に、第2のマトリックスの全中央画素値を加えること、及び、得られた値を画素数で割ることによって得られる。これらの平均値は、図10Dのイメージファイルの表示の下部及ぶ右側で曲線10及び11として示されている。6×4個の有力領域を有する基板を用いることは公知であるので、4個の最高平均列値と6個の最高平均行値とが、列中心線12と行中心線13の評価として選択される。これらの列及び行中心線12,13の交差点は、全有力領域の中心14の位置として用いられる。
【0026】
図11Aから図11Fはさらに、当該有力領域の中心を配置する方法の段階を示すものである。図11Aは、列中心線12と行中心線13とを有する当該領域の評価ウィンドウを示す。さらなる段階では、線11及び12の交差点14が真に当該領域の画素の重心15にあるか否かを決定する。この結果、重心は、図11Aで示したような評価ウィンドウ内の第2のマトリックスの画素の中央画素値から決定される。この重心が交差点14に一致しないならば、中心線12,13は、交差点が図11Bで示したような重心15に一致するようにずらされる。
【0027】
高い値を有する評価ウィンドウ内の画素が真に当該領域を表すか又はただのランダムノイズであるか、すなわち、基板のこの位置での結合実体は分析物に結合しなかったか否かを決定するために、以下の段階を実施する。図11C及び図11Eは、評価ウィンドウ内の2つの有力領域を示す。第1に、第2のマトリックスの評価ウィンドウ内の高い値を有する画素の表面を決定し、この表面から当該領域の半径を決定する。さらに、この表面の半径をこの表面の周縁から決定する。図11Cに示したように、高い値を有する画素は正確な円を表さず、そのため、この領域の周縁から決定された半径は図11Cで示した正確な円16の半径より大きい。しかしながら、2つの半径の比は、所定の参照値例えば、2より小さいべきである。図11Dは、高い値を有する画素の表面の周縁での高い値を有する全画素を示す。周縁から得られた半径と表面から得られた半径との比はこの場合、1.47である。これは、当該領域が存在することを意味する。図11Fでは、図11Eで示した表面の周縁での高い値を有する全画素を示し、周縁から得られた“半径”と図11Eで示した表面から得られた半径との比はこの場合、4.47である。これは、図11Eの場合における評価ウィンドウが当該領域を囲繞しない、すなわち、当該領域が存在しないことを意味する。
【0028】
当該領域の存在又は欠落を決定するために、代替の段階あるいはさらなる段階として、当該可能な領域の重心を、図11で示したような評価ウィンドウ内の高い値を有する全画素の表面に対応する表面を有する虚の(想像上の)円ウィンドウの中心としてとることも可能である。その後、この虚の円ウィンドウは、ある位置を見つけるために重心に対して移動でき、ここで、虚の円ウィンドウは、例えば、全画素の少なくとも90%で高い値を有する所定数の画素をカバーする。このような位置を見つけることができないならば、当該領域は存在しない。このような位置を見つけることができるならば、当該領域が存在する。
【0029】
当該可能な領域を調べるため、かつ、当該領域の信号強度を評価するためにのに用いられるサイズは、高い値を有する画素の数が低い数の画素の数に一致するように決定される。これは、ドット領域1の存在する場合では、ドット画素の数とバックグランド画素の数とが等しい、すなわち、全画素数の50%であることを意味する。評価ウィンドウ内にノイズだけが存在するならば、このような分布も存在し、すなわち、評価ウィンドウ内の画素の50%が平均値より上であり、50%が平均値より下である。このような評価ウィンドウは、ドット領域の信号強度を評価する上述の方法において高い信頼性でフィッティングを行うことができるという利点を示す。また、ユーザーがドット領域の中心を移動したいあるいは評価ウィンドウのサイズを変更したいならば、ユーザーは適当なユーザーインターフェースによってそうすることができる。
【0030】
本発明は、上述の実施形態に限定されず、特許請求の範囲内で多くの変更が可能である。
【図面の簡単な説明】
【図1】基板の散乱パラメータを決定するために用いる画像を示す説明図である。
【図2】図1の画像から得られた散乱減衰関数を示す図である。
【図3】各領域が公知のプローブを有する6×4領域のアレイを有する基板であって、基板上に標識材料の堆積後のイメージを示す図である。
【図4】散乱効果を機械的に除去した後の図3のイメージを示す図である。
【図5】図4の領域のうちの一つ上の評価ウィンドウを示す図である。
【図6】2つのフィッティングされた曲線を有する
【図7】本発明の方法の第2の実施形態においてフィッティングのために用いられる理論的散乱ヒストグラム決定を示す図である。
【図8】本発明の方法の第2の実施形態においてフィッティングのために用いられる理論的散乱ヒストグラム決定を示す図である。
【図9】図7及び図8で示された方法で得られるフィッティングされる理論的ヒストグラムと共に、図3の領域のマトリックス画素強度ISから得られた評価ヒストグラムを示す図である。
【図10】(A)図3で示した基板のイメージファイルを評価するために本発明に係る方法の好適な一実施形態の一段階の説明図である。(B)他の段階の説明図である。(C)他の段階の説明図である。(D)他の段階の説明図である。
【図11】(A)図3で示した基板のイメージファイルを評価するために本発明に係る方法の好適な一実施形態の一段階の説明図である。(B)他の段階の説明図である。(C)他の段階の説明図である。(D)他の段階の説明図である。(E)他の段階の説明図である。(F)他の段階の説明図である。
【符号の説明】
1 ドット領域
2 ライン
3 評価ウィンドウ
4 第1の分布ピーク
5 第2の分布ピーク
6,7 通常の分布曲線
8 散乱応答曲線
12 列中心線
13 行中心線
14 領域の中心
15 重心
Claims (19)
- 基板に埋め込まれている基板の少なくとも一の領域の信号強度であってバックグランド強度を有するものをコンピュータによって評価する方法において、
基板の散乱パラメータを決定し、前記領域及びその周囲部を囲む評価ウィンドウ内に位置する画素の画素強度のマトリックス(IS)を決定するものであり、ここで、マトリックス(IS)の画素強度は画素強度の評価ヒストグラムを得るのに用いられ、前記評価ヒストグラムは2つの分布ピークを示し、低い強度の第1の分布ピークは前記周囲の画素に対応し、高い強度の第2の分布ピークは前記領域の画素に対応し、2つの分布ピークを有する曲線は前記評価ヒストグラムでフィッティングされ、ここで、散乱パラメータを用いて散乱に対してマトリックス(IS)又はフィッティングされた曲線のいずれかが補正され、前記領域での画素の信号強度は前記評価ヒストグラムでフィッティングされた曲線から得られたデータによって決定されることを特徴とする方法。 - 画素強度のマトリックス(IS)を散乱パラメータを用いて処理してマトリックス(IS)のデコンボルーションを決定して非散乱画素強度のマトリックス(ID)を得するものであり、ここで、前記非散乱画素強度を用いて評価ヒストグラムが得られ、マトリックス(IS)の画素強度から得られたノイズパラメータによって決定された通常の分布曲線は前記の2つの分布のいずれかでフィッティングされ、ここで、第1のフィッティングされた分布曲線の平均値はバックグランド強度の平均画素強度とされ、かつ、第2のフィッティングされた分布曲線の平均値は領域強度の平均画素強度とされ、ここで、平均値の差は領域における画素の信号強度に対応する請求項1に記載の方法。
- 画素強度の理論上の散乱ヒストグラムが、予め決められた信号強度(IA)を有する理論上の領域と、マトリックス(IS)の値と散乱パラメータとマトリックス(IS)の画素強度から得られるノイズパラメータとを決定するために用いられる評価ウィンドウに対応する評価ウィンドウ内に位置するバックグランド強度(IB)を有する周囲とを用いて決定されるものであり、ここで、前記評価ヒストグラムがマトリックス(IS)の画素強度から成り、理論上の散乱ヒストグラムが少なくとも前記バックグランドと前記の理論上の散乱ヒストグラムとを変えることによって前記評価ヒストグラムでフィッティングされ、領域強度とフィッティングされた理論上のヒストグラムとの差が評価される領域における画素の信号強度に対応する請求項1に記載の方法。
- 規格化された散乱応答曲線の表を、規格化された理論上の領域と規格化された予め決められた強度とバックグランド強度なしの周囲とを用いて作成するものであり、ここで、各散乱応答曲線は、異なる基板と基板状態とに対して多数の散乱パラメータのそれぞれについて、前記の理論上の領域の中心から前記評価ウィンドウの端までの半径に沿って決定され、ここで、散乱応答曲線は、基板と使用される基板の基板状態とに対応する表から選択され、ここで、選択された散乱応答曲線を用いて使用される画素強度を決定して散乱ヒストグラムを作成し、ノイズパラメータをこの散乱ヒストグラムを用いてたたみ込んで理論上の散乱ヒストグラムを得る請求項3に記載の方法。
- 前記ノイズパラメータはマトリックス(IS)の画素強度に存在するノイズの標準偏差であり、ここで、ノイズの標準偏差は例えば、2つの隣接画素の間の差の分布を決定することによって決定するのが好ましい請求項2から4のいずれか一項に記載の方法。
- 散乱パラメータは、前記表面に平行な方向で基板の側部から照明することによって決定するものであり、ここで、前記照明方向において前記側部から表面上の一の線に沿った画素強度を決定して指数散乱減衰関数を得る請求項1から5のいずれか一項に記載の方法。
- 散乱パラメータは散乱なしの場合の信号強度に対して補正される請求項6に記載の方法。
- 基板の異なる状態について異なる基板の散乱パラメータが決定されかつ格納されるものであり、ここで、散乱パラメータは基板の状態に対応する評価における使用に対して選択される請求項1から7のいずれか一項に記載の方法。
- 基板は領域のアレイを有し、画素強度のマトリックスを各領域及び対応する周囲について決定する請求項1から8のいずれか一項に記載の方法。
- 基板の全表面を評価して、各領域の中心を位置決めしかつ評価ウィンドウのサイズを決定するものであり、各領域の位置決めはユーザーインターフェースを通して補正することができ、その後各領域の信号強度を評価する請求項9に記載の方法。
- 基板周囲に埋め込まれている基板の所望の有力領域をコンピュータを用いて位置決めする方法において、基板のイメージファイルをローパスフィルタアルゴリズムで処理してイメージファイルの全画素に対して局所的平均値のマトリックスを決定するものであり、ここで、局所的平均値のマトリックスはイメージファイルの実際の画素値のマトリックスと結合して、実際の画素値に依存する高又は低のいずれかが局所的平均値より上又は低いのいずれかに対応するところの高/低画素値の第1のマトリックスを得、ここで、第1のマトリックスの高/低画素値をさらに中央フィルタアルゴリズムで処理して中央画素値の第2のマトリックスを得、ここで、各中央画素値は第1のマトリックスの高/低画素値の多数と等しく、ここで、中央画素値を列毎及び行毎に処理して平均列値及び平均行値をそれぞれ決定し、ここで、最大平均列値及び最大平均行値は所望の有力領域の列中心線と行中心線の評価値として選択され、ここで、所望の有力領域の中心がそれら列中心線と行中心線の交点がであることを特徴とする方法。
- 局所的平均値のマトリックスを決定するローパスフィルタアルゴリズムは、所望の予想有力領域及びその周囲を囲繞する一のサイズを有する第2のウィンドウを用いること、イメージマトリックスに沿って画素毎に前記第2のウィンドウをずらすこと、及び、第2のウィンドウの各位置で第2のウィンドウ内の全画素の平均値を第2のウィンドウの中心での各画素の局所的平均値とすることを含む請求項11に記載の方法。
- 中央フィルタアルゴリズムは、各画素の高/低画素値を周囲画素の高/低画素値と比較することを含むものであり、ここで、高/低画素値は周囲画素の高/低画素の多数と等しくされている請求項11又は12のいずれかに記載の方法。
- 所望の有力領域の各中心について、第2のウィンドウ内の高い値を有する中央画素値から重心を決定するものであり、ここで、重心を所望の有力領域の中心とする請求項11から13のいずれか一項に記載の方法。
- 所望の有力領域の各中心について、高い値を有する中央画素値の表面とこの表面の周縁とから所望の有力領域の半径を決定するものであり、ここで、所望の領域の存在あるいは欠如を決めるために2つの半径の比を決定する請求項14に記載の方法。
- 所望の有力領域の重心を、高い値を有する中央画素値の表面に対応する表面を有する虚の円ウィンドウの中心とするものであり、ここで、虚の円ウィンドウは重心に対して移動して高い値を有する所定数の画素をカバーする位置を見つけ、ここで、所望の領域の存在あるいは欠如についての決定はそのような位置が見つかるか否かに依存して行われる請求項14又は15のいずれかに記載の方法。
- 所望の領域を評価する評価ウィンドウのサイズは、高い値を有する画素数が低い値を有する画素数に一致するように決定される請求項11から16のいずれか一項に記載の方法。
- 請求項1から17のいずれか一項に記載の方法を実施するように適合されたコンピュータプログラムコードを備え、コンピュータ上で作動するコンピュータプログラム。
- コンピュータ読込可能媒体あるいはコンピュータにダウンロード可能なファイルに格納された請求項17に記載のコンピュータプログラム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP00202551 | 2000-07-18 | ||
PCT/EP2001/008012 WO2002006854A1 (en) | 2000-07-18 | 2001-07-11 | Method for locating areas of interest of a substrate |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2004504659A true JP2004504659A (ja) | 2004-02-12 |
Family
ID=8171817
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002512708A Withdrawn JP2004504659A (ja) | 2000-07-18 | 2001-07-11 | 基板における所望領域の位置決め法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20040019433A1 (ja) |
EP (1) | EP1301806A1 (ja) |
JP (1) | JP2004504659A (ja) |
AU (1) | AU2001279739A1 (ja) |
WO (1) | WO2002006854A1 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8009889B2 (en) | 2006-06-27 | 2011-08-30 | Affymetrix, Inc. | Feature intensity reconstruction of biological probe array |
WO2011046614A2 (en) * | 2009-10-16 | 2011-04-21 | The Regents Of The University Of California | Methods and systems for phylogenetic analysis |
WO2012143909A1 (en) | 2011-04-20 | 2012-10-26 | Dynex Technologies, Inc. | Method of analysing reagent beads |
US20150103181A1 (en) * | 2013-10-16 | 2015-04-16 | Checkpoint Technologies Llc | Auto-flat field for image acquisition |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5790692A (en) * | 1994-09-07 | 1998-08-04 | Jeffrey H. Price | Method and means of least squares designed filters for image segmentation in scanning cytometry |
US5845007A (en) * | 1996-01-02 | 1998-12-01 | Cognex Corporation | Machine vision method and apparatus for edge-based image histogram analysis |
AU2899599A (en) * | 1998-03-05 | 1999-09-20 | Universal Health-Watch, Inc. | Optical imaging system for diagnostics |
-
2001
- 2001-07-11 US US10/333,365 patent/US20040019433A1/en not_active Abandoned
- 2001-07-11 AU AU2001279739A patent/AU2001279739A1/en not_active Abandoned
- 2001-07-11 JP JP2002512708A patent/JP2004504659A/ja not_active Withdrawn
- 2001-07-11 WO PCT/EP2001/008012 patent/WO2002006854A1/en not_active Application Discontinuation
- 2001-07-11 EP EP01957952A patent/EP1301806A1/en not_active Withdrawn
Also Published As
Publication number | Publication date |
---|---|
US20040019433A1 (en) | 2004-01-29 |
EP1301806A1 (en) | 2003-04-16 |
WO2002006854A1 (en) | 2002-01-24 |
AU2001279739A1 (en) | 2002-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3093270B2 (ja) | 光学的試料分析装置及び方法 | |
Wang et al. | Quantitative quality control in microarray image processing and data acquisition | |
JP4532261B2 (ja) | 生物学的試料のための光線画像分析 | |
US6341182B1 (en) | Method and apparatus for analyzing an image | |
US8712118B2 (en) | Automated measurement of concentration and/or amount in a biological sample | |
AU758339B2 (en) | Agglutination assays | |
US4741043A (en) | Method of and an apparatus for image analyses of biological specimens | |
CN110462057B (zh) | 数字分子测定 | |
EP2302363A2 (en) | Method for normalization of assay data | |
US20040246252A1 (en) | Method and apparatus for visualizing data | |
EP2227774A1 (en) | Method and system for image analysis | |
BE1025903B1 (nl) | Overspraakcorrectie bij multiplexing-analyse van biologisch monster | |
JP2005502369A (ja) | 分画細胞計数方法ならびにそれを実行するための関連する装置およびソフトウェア | |
US20140354986A1 (en) | Evaluation of focusing performance in an assay analysis system | |
CN115768559A (zh) | 用于提高灵敏度的诱导聚集测定 | |
CA2398795C (en) | Image metrics in the statistical analysis of dna microarray data | |
CN117388227B (zh) | 使用检测微球测定样品中目标分子浓度的方法、计算机可读介质和分析设备 | |
JP2014530594A (ja) | インサイツ分析用に調製された癌細胞の核病理学的特徴を測定するためのスペクトルイメージング | |
JP2014530594A5 (ja) | ||
JP2004504659A (ja) | 基板における所望領域の位置決め法 | |
CN117110270A (zh) | 使用检测微球测定样品中目标分子浓度的方法、计算机可读介质和分析设备 | |
WO2020005768A1 (en) | Plasmonic swarm biosensing system and methods of use | |
KR20220156564A (ko) | 자동화된 형광 이미징 및 단일 셀 분할 | |
US6832163B2 (en) | Methods of identifying heterogeneous features in an image of an array | |
Aung et al. | High-Plex Assessment of Biomarkers in Tumors |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20081007 |