画像の選択的デコンボリューション
画像の特徴間のクロストークを低減するためにデコンボリューションを選択的に使用する方法およびシステムを提供する。デコンボリューションのために画像の領域を選択する方法は、ａ）複数の特徴を有する画像を提供するステップであって、各特徴が少なくとも１つの値（ｖ）に関連する、ステップと、ｂ）画像の既知の低値ゾーンに隣接する高値特徴であるテスト特徴を識別するステップであって、テスト特徴が、テスト特徴の値（ｖｔ）の画像の隣接する低値ゾーンの値（Ｖｏ）に対する比率であるテール率（ｒｔ）を有する、ステップと、ｃ）テスト特徴のテール率（ｒｔ）の関数である閾値ｔを計算するステップと、ｄ）画像の選択領域を識別するステップであって、選択領域が、隣接する特徴の間の値（ｖ）の比率が上記閾値（Ｔ（ｒｔ））を上回る領域である、ステップと、を含む。通常、本発明の方法は、さらに、画像の選択領域に対してデコンボリューションを行うステップを含む。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理に関し、特に、画像の特徴間のクロストークを低減するためにデコンボリューションを選択的に使用することに関する。本発明は、デコンボリューション、すなわち通常集中的な計算を必要とするプロセスのために、関連領域を選択することにより、優れた画質を提供するために必要な計算労力を大幅に低減することができる。
【背景技術】
【０００２】
特許文献１は、画像の重心積分（ｃｅｎｔｒｏｉｄ ｉｎｔｅｇｒａｔｉｏｎ）の方法を開示している。特許文献２は、画像のオートグリッド解析について開示している。特許文献３は、画像の自動閾値処理の方法について開示している。
【特許文献１】米国特許第６，４７７，２７３号明細書
【特許文献２】米国特許第６，６３３，６６９号明細書
【発明の開示】
【課題を解決するための手段】
【０００３】
簡単に、本発明は、デコンボリューションのために画像の領域を選択する方法であって、ａ）複数の特徴を有する画像を提供するステップであって、各特徴が少なくとも１つの値（ｖ）に関連する、ステップと、ｂ）画像の既知の低値ゾーンに隣接する高値特徴であるテスト特徴を識別するステップであって、テスト特徴が、テスト特徴の値（ｖｔ）の画像の隣接する低値ゾーンの値（Ｖｏ）に対する比率であるテール率（ｒｔ）を有する、ステップと、ｃ）テスト特徴のテール率（ｒｔ）の関数である閾値ｔを計算するステップと、ｄ）画像の選択領域を識別するステップであって、選択領域が、隣接する特徴の間の値（ｖ）の比率が上記閾値（Ｔ（ｒｔ））を上回る領域である、ステップと、を含む、方法を提供する。画像は、通常、グリッドで配置された特徴を含む。通常、オートグリッド解析によって擬似画像が形成される。通常、ステップｂ）は、テール率（ｒｔ）を計算する前に、テスト特徴の値（ｖｔ）と画像の隣接する低値ゾーンの値（Ｖｏ）との両方から背景定数を減算することをさらに含む。背景定数を、任意に、いかなるテール効果も回避するようにいかなる特徴からも十分に離れている画像の低値ゾーンの（ｖｂ）の値であるようにとってもよく、そのゾーンは、任意に、いかなる特徴からも特徴間の平均距離の少なくとも２倍離れている画像の低値ゾーンであってもよい。通常、閾値（Ｔ（ｒｔ））は、テスト特徴のテール率（ｒｔ）の倍数である。通常、本発明の方法は、さらに、画像の選択領域に対してデコンボリューションを行うステップを含む。
【０００４】
別の態様では、本発明は、デコンボリューションのために画像の領域を選択するシステムであって、ａ）デジタル化画像を提供する画像装置と、ｂ）データ記憶装置と、ｃ）画像装置からデジタル化画像を受け取り、データ記憶装置に対してデータを書き込みかつ読み出すことができる中央処理装置であって、
ｉ）画像装置からデジタル化画像を受け取り、
ｉｉ）複数の特徴を識別し各特徴を少なくとも１つの値（ｖ）に関連付け、
ｉｉｉ）画像の既知の低値ゾーンに隣接する高値特徴であるテスト特徴を識別し、テスト特徴が、テスト特徴の値（ｖｔ）の画像の隣接する低値ゾーンの値（ｖｏ）に対する比率であるテール率（ｒｔ）を有し、
ｉｖ）テスト特徴のテール率（ｒｔ）の関数である閾値ｔを計算し、
ｖ）画像の選択領域を識別し、選択領域が画像全体に満たない領域を含み、選択領域が、隣接する特徴の間の値（ｖ）の比率が該閾値（Ｔ（ｒｔ））を上回る領域である
ようにプログラムされた中央処理装置と、
を備える、システムを提供する。画像は、通常、グリッドで配置された特徴を含む。通常、中央処理装置は、オートグリッド解析によって擬似画像を形成するようにさらにプログラムされる。通常、ステップｉｉｉ）は、テール率（ｒｔ）を計算する前に、テスト特徴の値（ｖｔ）と画像の隣接する低値ゾーンの値（Ｖｏ）との両方から背景定数を減算することをさらに含む。背景定数を、任意に、いかなるテール効果も回避するようにいかなる特徴からも十分に離れている画像の低値ゾーンの（ｖｂ）の値であるようにとってもよく、そのゾーンは、任意に、いかなる特徴からも特徴間の平均距離の少なくとも２倍離れている画像の低値ゾーンであってもよい。通常、閾値（Ｔ（ｒｔ））は、テスト特徴のテール率（ｒｔ）の倍数である。通常、中央処理装置は、画像の選択領域に対してデコンボリューションを行うようにさらにプログラムされる。
【発明の効果】
【０００５】
本発明の効果は、画像から高品質データを導出するために必要な計算労力を低減する方法を提供することである。
【発明を実施するための最良の形態】
【０００６】
本発明は、デコンボリューションのために画像の領域を選択する方法を提供する。反復法およびブラインド法を含む、本技術分野において既知のデコンボリューションの任意の適当な方法を使用することができる。反復法には、リチャードソン・ルーシー（Ｒｉｃｈａｒｄｓｏｎ−Ｌｕｃｙ）および反復条件付きチホノフ・ミラー（Ｉｔｅｒａｔｉｖｅ Ｃｏｎｓｔｒａｉｎｅｄ Ｔｉｋｈｏｖａｎ−Ｍｉｌｌｅｒ）法がある。ブラインド法には、ウィナーフィルタ（Ｗｅｉｎｅｒ Ｆｉｌｔｅｒｉｎｇ）、シミュレーテッドアニーリング（Ｓｉｍｕｌａｔｅｄ Ａｎｎｅａｌｉｎｇ）および最大尤度推定器（Ｍａｘｉｍｕｍ Ｌｉｋｅｌｉｈｏｏｄ Ｅｓｔｉｍａｔｏｒ）法がある。デコンボリューションは、相対的に暗い特徴が明るい特徴に近接することによって間違って明るくなること等、画像の特徴の間のクロストークを低減することができる。
【０００７】
選択の方法は、ａ）複数の特徴を含む画像を提供するステップであって、各特徴が少なくとも１つの値（ｖ）に関連する、ステップと、ｂ）画像の既知の低値ゾーンに隣接する高値特徴であるテスト特徴を識別するステップであって、テスト特徴が、テスト特徴の値（ｖｔ）の画像の隣接する低値ゾーンの値（ｖｏ）に対する比率であるテール率（ｒｔ）を有する、ステップと、ｃ）テスト特徴のテール率（ｒｔ）の関数である閾値ｔを計算するステップと、ｄ）画像の選択領域を識別するステップであって、選択領域が、隣接する特徴の間の値（ｖ）の比率が上記閾値（Ｔ（ｒｔ））を上回る領域である、ステップと、を含む。通常、１つまたは複数のステップは自動化される。より一般的には、すべてのステップが自動化される。
【０００８】
画像を提供するステップは、任意の適当な方法で実現することができる。通常、このステップは自動化される。画像を、ビデオカメラ、デジタルカメラ、光化学カメラ、マイクロスコープ、テレスコープ、視覚的走査システム、プローブ走査システム、または２次元アレイでデータ点を生成する他の検知装置を使用して収集することができる。通常、目標画像は、明確な特徴を含む画像であることが期待されるが、それはさらに雑音を含む可能性もある。通常、特徴は、行および列を含むグリッドで配置される。本明細書では、「列」は、一方向における特徴の全般的な配列を示すために使用し、「行」は、列に対して略直交する方向における特徴の全般的な配列を示すために使用する。いずれの方向が列であるか、およびいずれの方向が行であるかは完全に任意であり、そのため、一方の用語を用いて他方の用語を用いないことには、いかなる重要性も置くべきではなく、かつ、行および列は完全に直線状でなくてもよい、ということが理解されよう。別法として、グリッドは、三角形配置または六角形配置等、特徴の他の何らの繰返しの幾何学的配置を含んでもよい。別法として、特徴を、天文学上の画像におけるような事前に確定されていないパターンで配置してもよい。画像は、最初に、画像取込機器または画像作成機器によりデジタル形式で作成されていない場合、通常、画素になるようにデジタル化される。通常、本明細書で説明する方法は、中央処理装置またはコンピュータを使用して実現される。
【０００９】
図１は、本発明をともに使用することができるスキャニングシステムを示す。図１のシステムでは、集束光線が物体を横切って移動し、システムは、結果として生じる反射光または蛍光を検出する。これを行うために、光源１０からの光は、光源光学系１２を通して集束され、ミラー１４によってここではサンプル３×４アッセイプレート１６として示す物体上に偏向される。モータ２４を使用してミラー１４の位置を変更することにより、光源１０からの光をサンプルの種々の場所に向けることができる。蛍光を発するかまたはサンプル１６から反射する光は、通常ハーフミラーであるミラー１５を介して検出光学系１８に戻る。別法として、米国特許第５，９００，９４９号明細書に示されているように、光源を中央に与えることができ、かつ放射光または蛍光をシステムの側方から検出することができ、もしくは、光源をシステムの側方から与えることができ、かつ放出光または蛍光を中央で検出することができ、もしくは他の同様の変形とすることも可能である。検出光学系１８を通過する光は、テレビカメラ、ＣＣＤ、レーザ反射システム、光電子増倍管、アバランシェフォトダイオード、フォトダイオードまたは単光子計数モジュール等の任意の適当な画像取込システム２０を使用して検出され、それら画像取込システムの出力は、解析のためかつシステム全体を制御するようにプログラムされたコンピュータ２２に提供される。コンピュータ２２は、通常、プログラムを実行する中央処理装置と、データ格納のためのＲＡＭ、ハードドライブ等のシステムと、を含む。この説明は例示のみを目的とするものであることが理解されよう。すなわち、本発明を、光ベースの画像だけではなく磁気センサまたは触覚センサから生成される「シミュレートされた」画像で、かつサンプル１６だけではなく検査される任意の物体で等しく使用することができる。
【００１０】
米国特許第６，４７７，２７３号明細書および同第６，６３３，６６９号明細書に記載されているように、画像に対して、さらなる解析の前に重心積分およびオートグリッド解析をほどこしてもよい。各特徴に対し、そこではその「値」として提供される積分された強度を割り当ててもよく、または、特徴値としての局所極大値等の選択を含む場合もある他の任意の適当な方法によって値を割り当ててもよい。オートグリッド解析によって形成される擬似画像を生成してもよい。
【００１１】
本明細書では、「高値」および「低値」を、写真画像における明るい特徴および暗い特徴に関連して使用する。「高値」、「低値」および「値」という用語を、画像および画像を収集するために使用される装置の性質に応じて、限定なしに色値、ｘ線透過値、電波放射値等を含む、画像に表されうる任意の特徴に適用してもよい。通常、「高値」は、画像収集装置の性質に応じて、隣接する「低値」特徴においてクロストークをもたらす傾向にある特徴を指す。
【００１２】
テスト特徴を識別するステップを、任意の適当な方法で実現してもよい。通常、このステップは自動化される。テスト特徴は、画像の既知の低値ゾーンに隣接する高値特徴である。低値ゾーンは、低値特徴、またはエッジ領域、または特徴が期待される領域の外側であることが分かっている他の領域等、低値であることが分かっている領域であってもよい。一実施形態では、特徴の期待されるグリッドのエッジを構成する特徴が検査され、明るいエッジ特徴がテスト特徴として選択される。テスト特徴として選択される特徴は、一組の候補の最高値であってもよく、または、事前に選択された閾値を超える、最初に検査された値であってもよい。別の実施形態では、画像形成される物体に、基準点としての役割を果たすように隣接する高値特徴および低値特徴が設けられる。
【００１３】
テスト特徴の値（ｖｔ）を画像の隣接する低値ゾーンの値（ｖｏ）で除算することによって、テール率（ｔａｉｌ ｒａｔｉｏ）（ｒｔ）が計算される。通常、テール率（ｒｔ）を計算する前に、テスト特徴の値（ｖｔ）と画像の隣接する低値ゾーンの値（ｖｏ）との両方から背景定数が減算される。背景定数を、任意の適当な方法によって確定してもよい。背景定数を、いかなるテール効果も回避するようにいかなる特徴からも十分に離れている画像の低値ゾーンの（ｖｂ）の値であるようにとってもよい。特徴がグリッドで配置される場合、その離れている低値ゾーンは、通常、いかなる特徴からも特徴間の平均距離の少なくとも２倍離れている。別法として、背景定数は、所与の装置に適当であるように事前に確定された一定の値であってもよい。
【００１４】
テスト特徴のテール率（ｒｔ）の関数である閾値ｔが計算される。算術、対数、指数、三角法等の関数を含む任意の適当な関数を使用してもよい。通常、閾値（Ｔ（ｒｔ））は、単にテール率（ｒｔ）の倍数であり、すなわち、Ａが任意の適当な数であるがもっとも一般的には２と２０との間である場合、Ｔ（ｒｔ）＝Ａ×ｒｔである。
【００１５】
そして、閾値ｔを使用して、任意の適当な方法により画像の選択された領域が識別される。通常、このステップは自動化される。最も一般的には、選択された領域は、隣接する特徴の間の値（ｖ）の比率が上記閾値（Ｔ（ｒｔ））を上回る領域である。
【００１６】
本発明は、特に、ＤＮＡ検出または分類のような自動化された解析プロセスに含まれる可能性のある、トレー、スライド等の上のサンプル点のマトリクスの自動読取りにおける、光学情報の自動読取りに有用である。別法として、本発明は、天文学、医用画像処理、リアルタイム画像解析等において有用であり得る。特に、本発明は、過度の計算なしに画像のデコンボリューションにより空間クロストークを低減するのに有用である。
【００１７】
本発明の目的および利点について、以下の例によってさらに説明するが、これらの例で列挙する方法ステップの特定の順序および詳細は、他の条件および詳細と同様に、本発明を不当に限定するものと解釈されるべきではない。
【００１８】
（例）
この例で使用する対象画像を図２に示す。画像は、７４×６２画素のサイズであり、１０列および９行で配置された特徴を示す。各画素の明度は、強度値によって表される。
【００１９】
まず、米国特許第６，４７７，２７３号明細書および同第６３３，６６９号明細書に示すように、米国特許第６，６３３，６６９号明細書に記載された「曲げ（ｆｌｅｘｉｎｇ）」を含むオートグリッド解析を画像に対して行うことにより、図３に示す解析グリッドを作成し、各特徴に対して積分強度を割り当てた。表Ｉは、各列および行位置に対する積分強度値を報告するものである。
【００２０】
表Ｉ
【００２１】
テスト特徴として、列１、行Ｅにおける明るいエッジ特徴を選択した。図４は、各画素から背景定数を減算した後のこの特徴および隣接する暗いゾーンの拡大図である。背景定数を、任意の明るい特徴からの略最大距離にある画像のエッジにおける画素の小さいグループの平均強度値として取得した。図５は、ｘ方向におけるテスト特徴のテールを示すグラフである。各ｘ位置に対し、グラフは、ｙ方向における４つの画素にわたって積分された強度値を報告している。このテスト特徴に対するテール率は、テスト特徴から１特徴幅（５画素）離れて中心を置く隣接する暗いゾーンのある領域に亙る積分強度（２５、図５の画素２〜５に亙って積分）と、テスト特徴に亙る積分強度（１４８９、図５の画素７〜１０に亙って積分）と、の比であり、すなわち０．０１６８である。
【００２２】
閾値を、テール率の１０倍、すなわち０．１６８とした。このため、目的は、隣接する明るい特徴からの予測される寄与の１０倍未満の明るさ、すなわち１０×隣接する特徴の明度（ａ）×テール率未満の強度（ｂ）を有する特徴を選択することである。この条件を、式Ｉ：ｂ＜ａ×１０×（テール率）またはｂ＜ａ×（閾値）に表すことができる。
【００２３】
連続した演算を簡略化するために、積分強度値および閾値を対数に変換した。表ＩＩは、各列および行位置に対して表Ｉで報告された積分強度値の自然対数を含む。ｌｎ（閾値）の値は−１．７８であった。式Ｉを、式ＩＩにおいて対数に関して表し、すなわちｌｎ（ｂ）＜ｌｎ（ａ）＋ｌｎ（閾値）となり、これは−ｌｎ（閾値）＜ｌｎ（ａ）−ｌｎ（ｂ）となる。明／暗遷移と暗／明遷移との両方を検出するために明度差の絶対値をとると、式ＩＩは式ＩＩＩ：−ｌｎ（閾値）＜｜ｌｎ（ａ）−ｌｎ（ｂ）｜となる。
【００２４】
表ＩＩ
【００２５】
表ＩＩＩは、ｘ方向における表ＩＩの隣接する値の間の差の絶対値、すなわち｜ｌｎ（ａ）−ｌｎ（ｂ）｜を報告している。したがって、表ＩＩＩは、９つの列および９つの行を含む。表ＩＩＩの値を、表の最大値、すなわち２．９１１で除算することにより、１．０００に正規化した。正規化された値は表ＩＶに報告されている。１．７８という−ｌｎ（閾値）の値を、１．７８／２．９１１＝０．６１に正規化した。正規化された閾値を表ＩＶに適用することにより表Ｖを生成した。この表Ｖは、−ｌｎ（閾値）すなわち０．６１を下回る値に対して０を報告し、−ｌｎ（閾値）すなわち０．６１を上回る値に対して１を報告している。
【００２６】
表ＩＩＩ
【００２７】
表ＩＶ
【００２８】
表Ｖ
【００２９】
表ＶＩは、ｙ方向における表ＩＩの隣接する値の差の絶対値、すなわち｜ｌｎ（ａ）−ｌｎ（ｂ）｜を報告している。したがって、表ＶＩは、１０の列および８つの行を含む。表ＶＩの値を、表の最大値、すなわち３．２７５１で除算することにより、１．０００に正規化した。正規化された値は表ＶＩＩに報告されている。１．７８の−ｌｎ（閾値）の値を、１．７８／３．２７５１＝０．５４に正規化した。正規化された閾値を表ＶＩＩに適用することにより表ＶＩＩＩを生成した。この表ＶＩＩＩは、−ｌｎ（閾値）すなわち０．５４を下回る値に対して０を報告し、−ｌｎ（閾値）すなわち０．５４を上回る値に対して１を報告している。
【００３０】
表ＶＩ
【００３１】
表ＶＩＩ
【００３２】
表ＶＩＩＩ
【００３３】
表Ｖを、カーネル
［１ １］
でコンボリューションすることにより、９×１０マトリクスを作成した（表ＩＸ）。ここでは、非ゼロエントリは、ｘ方向における明から暗への遷移または暗から明への遷移を示す。
【００３４】
表ＩＸ
【００３５】
表ＶＩＩＩをカーネル
でコンボリューションすることにより、９×１０マトリクスを作成した（表Ｘ）。ここでは、非ゼロエントリは、ｙ方向における明から暗への遷移または暗から明への遷移を示す。
【００３６】
表Ｘ
【００３７】
表ＩＸおよびＸによって表されるマトリクスを加算した結果、表ＸＩとして報告されるマトリクスとなった。
【００３８】
表ＸＩ
【００３９】
表ＸＩの非ゼロ値のすべてを包含する４つの矩形領域（Ａ１：Ｂ３、Ｄ５：Ｆ７、Ｈ１：Ｉ３、Ｉ９：Ｉ１０）をデコンボリューションのために選択した。選択された領域は、９０の特徴のうちの２３を含み、デコンボリューションの多くの方法によって、計算労力の程度が解析される領域のサイズに応じて指数関数的に上昇するようになるため、画像全体のデコンボリューションに必要となる計算労力の少なくとも約７４％、おそらくははるかにそれより多くが省かれる。
【００４０】
当業者には、本発明の範囲および原理から逸脱することなく本発明のさまざまな変更および改変が明らかとなろう。また、本発明は、上述した例示的な実施形態に不当に限定されるべきではない、ということが理解されるべきである。
【図面の簡単な説明】
【００４１】
【図１】本発明を使用してもよい典型的なスキャニングシステムの概略図である。
【図２】例で使用される対象画像である。
【図３】例で説明するような図２の画像の解析グリッドである。
【図４】図２の第１列第５行の特徴を含む図２の詳細である。
【図５】図４のセグメントのｘ位置に対してプロットされたｙ方向における４画素にわたって積分された画素強度のグラフである。
【技術分野】
【０００１】
本発明は、画像処理に関し、特に、画像の特徴間のクロストークを低減するためにデコンボリューションを選択的に使用することに関する。本発明は、デコンボリューション、すなわち通常集中的な計算を必要とするプロセスのために、関連領域を選択することにより、優れた画質を提供するために必要な計算労力を大幅に低減することができる。
【背景技術】
【０００２】
特許文献１は、画像の重心積分（ｃｅｎｔｒｏｉｄ ｉｎｔｅｇｒａｔｉｏｎ）の方法を開示している。特許文献２は、画像のオートグリッド解析について開示している。特許文献３は、画像の自動閾値処理の方法について開示している。
【特許文献１】米国特許第６，４７７，２７３号明細書
【特許文献２】米国特許第６，６３３，６６９号明細書
【発明の開示】
【課題を解決するための手段】
【０００３】
簡単に、本発明は、デコンボリューションのために画像の領域を選択する方法であって、ａ）複数の特徴を有する画像を提供するステップであって、各特徴が少なくとも１つの値（ｖ）に関連する、ステップと、ｂ）画像の既知の低値ゾーンに隣接する高値特徴であるテスト特徴を識別するステップであって、テスト特徴が、テスト特徴の値（ｖｔ）の画像の隣接する低値ゾーンの値（Ｖｏ）に対する比率であるテール率（ｒｔ）を有する、ステップと、ｃ）テスト特徴のテール率（ｒｔ）の関数である閾値ｔを計算するステップと、ｄ）画像の選択領域を識別するステップであって、選択領域が、隣接する特徴の間の値（ｖ）の比率が上記閾値（Ｔ（ｒｔ））を上回る領域である、ステップと、を含む、方法を提供する。画像は、通常、グリッドで配置された特徴を含む。通常、オートグリッド解析によって擬似画像が形成される。通常、ステップｂ）は、テール率（ｒｔ）を計算する前に、テスト特徴の値（ｖｔ）と画像の隣接する低値ゾーンの値（Ｖｏ）との両方から背景定数を減算することをさらに含む。背景定数を、任意に、いかなるテール効果も回避するようにいかなる特徴からも十分に離れている画像の低値ゾーンの（ｖｂ）の値であるようにとってもよく、そのゾーンは、任意に、いかなる特徴からも特徴間の平均距離の少なくとも２倍離れている画像の低値ゾーンであってもよい。通常、閾値（Ｔ（ｒｔ））は、テスト特徴のテール率（ｒｔ）の倍数である。通常、本発明の方法は、さらに、画像の選択領域に対してデコンボリューションを行うステップを含む。
【０００４】
別の態様では、本発明は、デコンボリューションのために画像の領域を選択するシステムであって、ａ）デジタル化画像を提供する画像装置と、ｂ）データ記憶装置と、ｃ）画像装置からデジタル化画像を受け取り、データ記憶装置に対してデータを書き込みかつ読み出すことができる中央処理装置であって、
ｉ）画像装置からデジタル化画像を受け取り、
ｉｉ）複数の特徴を識別し各特徴を少なくとも１つの値（ｖ）に関連付け、
ｉｉｉ）画像の既知の低値ゾーンに隣接する高値特徴であるテスト特徴を識別し、テスト特徴が、テスト特徴の値（ｖｔ）の画像の隣接する低値ゾーンの値（ｖｏ）に対する比率であるテール率（ｒｔ）を有し、
ｉｖ）テスト特徴のテール率（ｒｔ）の関数である閾値ｔを計算し、
ｖ）画像の選択領域を識別し、選択領域が画像全体に満たない領域を含み、選択領域が、隣接する特徴の間の値（ｖ）の比率が該閾値（Ｔ（ｒｔ））を上回る領域である
ようにプログラムされた中央処理装置と、
を備える、システムを提供する。画像は、通常、グリッドで配置された特徴を含む。通常、中央処理装置は、オートグリッド解析によって擬似画像を形成するようにさらにプログラムされる。通常、ステップｉｉｉ）は、テール率（ｒｔ）を計算する前に、テスト特徴の値（ｖｔ）と画像の隣接する低値ゾーンの値（Ｖｏ）との両方から背景定数を減算することをさらに含む。背景定数を、任意に、いかなるテール効果も回避するようにいかなる特徴からも十分に離れている画像の低値ゾーンの（ｖｂ）の値であるようにとってもよく、そのゾーンは、任意に、いかなる特徴からも特徴間の平均距離の少なくとも２倍離れている画像の低値ゾーンであってもよい。通常、閾値（Ｔ（ｒｔ））は、テスト特徴のテール率（ｒｔ）の倍数である。通常、中央処理装置は、画像の選択領域に対してデコンボリューションを行うようにさらにプログラムされる。
【発明の効果】
【０００５】
本発明の効果は、画像から高品質データを導出するために必要な計算労力を低減する方法を提供することである。
【発明を実施するための最良の形態】
【０００６】
本発明は、デコンボリューションのために画像の領域を選択する方法を提供する。反復法およびブラインド法を含む、本技術分野において既知のデコンボリューションの任意の適当な方法を使用することができる。反復法には、リチャードソン・ルーシー（Ｒｉｃｈａｒｄｓｏｎ−Ｌｕｃｙ）および反復条件付きチホノフ・ミラー（Ｉｔｅｒａｔｉｖｅ Ｃｏｎｓｔｒａｉｎｅｄ Ｔｉｋｈｏｖａｎ−Ｍｉｌｌｅｒ）法がある。ブラインド法には、ウィナーフィルタ（Ｗｅｉｎｅｒ Ｆｉｌｔｅｒｉｎｇ）、シミュレーテッドアニーリング（Ｓｉｍｕｌａｔｅｄ Ａｎｎｅａｌｉｎｇ）および最大尤度推定器（Ｍａｘｉｍｕｍ Ｌｉｋｅｌｉｈｏｏｄ Ｅｓｔｉｍａｔｏｒ）法がある。デコンボリューションは、相対的に暗い特徴が明るい特徴に近接することによって間違って明るくなること等、画像の特徴の間のクロストークを低減することができる。
【０００７】
選択の方法は、ａ）複数の特徴を含む画像を提供するステップであって、各特徴が少なくとも１つの値（ｖ）に関連する、ステップと、ｂ）画像の既知の低値ゾーンに隣接する高値特徴であるテスト特徴を識別するステップであって、テスト特徴が、テスト特徴の値（ｖｔ）の画像の隣接する低値ゾーンの値（ｖｏ）に対する比率であるテール率（ｒｔ）を有する、ステップと、ｃ）テスト特徴のテール率（ｒｔ）の関数である閾値ｔを計算するステップと、ｄ）画像の選択領域を識別するステップであって、選択領域が、隣接する特徴の間の値（ｖ）の比率が上記閾値（Ｔ（ｒｔ））を上回る領域である、ステップと、を含む。通常、１つまたは複数のステップは自動化される。より一般的には、すべてのステップが自動化される。
【０００８】
画像を提供するステップは、任意の適当な方法で実現することができる。通常、このステップは自動化される。画像を、ビデオカメラ、デジタルカメラ、光化学カメラ、マイクロスコープ、テレスコープ、視覚的走査システム、プローブ走査システム、または２次元アレイでデータ点を生成する他の検知装置を使用して収集することができる。通常、目標画像は、明確な特徴を含む画像であることが期待されるが、それはさらに雑音を含む可能性もある。通常、特徴は、行および列を含むグリッドで配置される。本明細書では、「列」は、一方向における特徴の全般的な配列を示すために使用し、「行」は、列に対して略直交する方向における特徴の全般的な配列を示すために使用する。いずれの方向が列であるか、およびいずれの方向が行であるかは完全に任意であり、そのため、一方の用語を用いて他方の用語を用いないことには、いかなる重要性も置くべきではなく、かつ、行および列は完全に直線状でなくてもよい、ということが理解されよう。別法として、グリッドは、三角形配置または六角形配置等、特徴の他の何らの繰返しの幾何学的配置を含んでもよい。別法として、特徴を、天文学上の画像におけるような事前に確定されていないパターンで配置してもよい。画像は、最初に、画像取込機器または画像作成機器によりデジタル形式で作成されていない場合、通常、画素になるようにデジタル化される。通常、本明細書で説明する方法は、中央処理装置またはコンピュータを使用して実現される。
【０００９】
図１は、本発明をともに使用することができるスキャニングシステムを示す。図１のシステムでは、集束光線が物体を横切って移動し、システムは、結果として生じる反射光または蛍光を検出する。これを行うために、光源１０からの光は、光源光学系１２を通して集束され、ミラー１４によってここではサンプル３×４アッセイプレート１６として示す物体上に偏向される。モータ２４を使用してミラー１４の位置を変更することにより、光源１０からの光をサンプルの種々の場所に向けることができる。蛍光を発するかまたはサンプル１６から反射する光は、通常ハーフミラーであるミラー１５を介して検出光学系１８に戻る。別法として、米国特許第５，９００，９４９号明細書に示されているように、光源を中央に与えることができ、かつ放射光または蛍光をシステムの側方から検出することができ、もしくは、光源をシステムの側方から与えることができ、かつ放出光または蛍光を中央で検出することができ、もしくは他の同様の変形とすることも可能である。検出光学系１８を通過する光は、テレビカメラ、ＣＣＤ、レーザ反射システム、光電子増倍管、アバランシェフォトダイオード、フォトダイオードまたは単光子計数モジュール等の任意の適当な画像取込システム２０を使用して検出され、それら画像取込システムの出力は、解析のためかつシステム全体を制御するようにプログラムされたコンピュータ２２に提供される。コンピュータ２２は、通常、プログラムを実行する中央処理装置と、データ格納のためのＲＡＭ、ハードドライブ等のシステムと、を含む。この説明は例示のみを目的とするものであることが理解されよう。すなわち、本発明を、光ベースの画像だけではなく磁気センサまたは触覚センサから生成される「シミュレートされた」画像で、かつサンプル１６だけではなく検査される任意の物体で等しく使用することができる。
【００１０】
米国特許第６，４７７，２７３号明細書および同第６，６３３，６６９号明細書に記載されているように、画像に対して、さらなる解析の前に重心積分およびオートグリッド解析をほどこしてもよい。各特徴に対し、そこではその「値」として提供される積分された強度を割り当ててもよく、または、特徴値としての局所極大値等の選択を含む場合もある他の任意の適当な方法によって値を割り当ててもよい。オートグリッド解析によって形成される擬似画像を生成してもよい。
【００１１】
本明細書では、「高値」および「低値」を、写真画像における明るい特徴および暗い特徴に関連して使用する。「高値」、「低値」および「値」という用語を、画像および画像を収集するために使用される装置の性質に応じて、限定なしに色値、ｘ線透過値、電波放射値等を含む、画像に表されうる任意の特徴に適用してもよい。通常、「高値」は、画像収集装置の性質に応じて、隣接する「低値」特徴においてクロストークをもたらす傾向にある特徴を指す。
【００１２】
テスト特徴を識別するステップを、任意の適当な方法で実現してもよい。通常、このステップは自動化される。テスト特徴は、画像の既知の低値ゾーンに隣接する高値特徴である。低値ゾーンは、低値特徴、またはエッジ領域、または特徴が期待される領域の外側であることが分かっている他の領域等、低値であることが分かっている領域であってもよい。一実施形態では、特徴の期待されるグリッドのエッジを構成する特徴が検査され、明るいエッジ特徴がテスト特徴として選択される。テスト特徴として選択される特徴は、一組の候補の最高値であってもよく、または、事前に選択された閾値を超える、最初に検査された値であってもよい。別の実施形態では、画像形成される物体に、基準点としての役割を果たすように隣接する高値特徴および低値特徴が設けられる。
【００１３】
テスト特徴の値（ｖｔ）を画像の隣接する低値ゾーンの値（ｖｏ）で除算することによって、テール率（ｔａｉｌ ｒａｔｉｏ）（ｒｔ）が計算される。通常、テール率（ｒｔ）を計算する前に、テスト特徴の値（ｖｔ）と画像の隣接する低値ゾーンの値（ｖｏ）との両方から背景定数が減算される。背景定数を、任意の適当な方法によって確定してもよい。背景定数を、いかなるテール効果も回避するようにいかなる特徴からも十分に離れている画像の低値ゾーンの（ｖｂ）の値であるようにとってもよい。特徴がグリッドで配置される場合、その離れている低値ゾーンは、通常、いかなる特徴からも特徴間の平均距離の少なくとも２倍離れている。別法として、背景定数は、所与の装置に適当であるように事前に確定された一定の値であってもよい。
【００１４】
テスト特徴のテール率（ｒｔ）の関数である閾値ｔが計算される。算術、対数、指数、三角法等の関数を含む任意の適当な関数を使用してもよい。通常、閾値（Ｔ（ｒｔ））は、単にテール率（ｒｔ）の倍数であり、すなわち、Ａが任意の適当な数であるがもっとも一般的には２と２０との間である場合、Ｔ（ｒｔ）＝Ａ×ｒｔである。
【００１５】
そして、閾値ｔを使用して、任意の適当な方法により画像の選択された領域が識別される。通常、このステップは自動化される。最も一般的には、選択された領域は、隣接する特徴の間の値（ｖ）の比率が上記閾値（Ｔ（ｒｔ））を上回る領域である。
【００１６】
本発明は、特に、ＤＮＡ検出または分類のような自動化された解析プロセスに含まれる可能性のある、トレー、スライド等の上のサンプル点のマトリクスの自動読取りにおける、光学情報の自動読取りに有用である。別法として、本発明は、天文学、医用画像処理、リアルタイム画像解析等において有用であり得る。特に、本発明は、過度の計算なしに画像のデコンボリューションにより空間クロストークを低減するのに有用である。
【００１７】
本発明の目的および利点について、以下の例によってさらに説明するが、これらの例で列挙する方法ステップの特定の順序および詳細は、他の条件および詳細と同様に、本発明を不当に限定するものと解釈されるべきではない。
【００１８】
（例）
この例で使用する対象画像を図２に示す。画像は、７４×６２画素のサイズであり、１０列および９行で配置された特徴を示す。各画素の明度は、強度値によって表される。
【００１９】
まず、米国特許第６，４７７，２７３号明細書および同第６３３，６６９号明細書に示すように、米国特許第６，６３３，６６９号明細書に記載された「曲げ（ｆｌｅｘｉｎｇ）」を含むオートグリッド解析を画像に対して行うことにより、図３に示す解析グリッドを作成し、各特徴に対して積分強度を割り当てた。表Ｉは、各列および行位置に対する積分強度値を報告するものである。
【００２０】
表Ｉ
【００２１】
テスト特徴として、列１、行Ｅにおける明るいエッジ特徴を選択した。図４は、各画素から背景定数を減算した後のこの特徴および隣接する暗いゾーンの拡大図である。背景定数を、任意の明るい特徴からの略最大距離にある画像のエッジにおける画素の小さいグループの平均強度値として取得した。図５は、ｘ方向におけるテスト特徴のテールを示すグラフである。各ｘ位置に対し、グラフは、ｙ方向における４つの画素にわたって積分された強度値を報告している。このテスト特徴に対するテール率は、テスト特徴から１特徴幅（５画素）離れて中心を置く隣接する暗いゾーンのある領域に亙る積分強度（２５、図５の画素２〜５に亙って積分）と、テスト特徴に亙る積分強度（１４８９、図５の画素７〜１０に亙って積分）と、の比であり、すなわち０．０１６８である。
【００２２】
閾値を、テール率の１０倍、すなわち０．１６８とした。このため、目的は、隣接する明るい特徴からの予測される寄与の１０倍未満の明るさ、すなわち１０×隣接する特徴の明度（ａ）×テール率未満の強度（ｂ）を有する特徴を選択することである。この条件を、式Ｉ：ｂ＜ａ×１０×（テール率）またはｂ＜ａ×（閾値）に表すことができる。
【００２３】
連続した演算を簡略化するために、積分強度値および閾値を対数に変換した。表ＩＩは、各列および行位置に対して表Ｉで報告された積分強度値の自然対数を含む。ｌｎ（閾値）の値は−１．７８であった。式Ｉを、式ＩＩにおいて対数に関して表し、すなわちｌｎ（ｂ）＜ｌｎ（ａ）＋ｌｎ（閾値）となり、これは−ｌｎ（閾値）＜ｌｎ（ａ）−ｌｎ（ｂ）となる。明／暗遷移と暗／明遷移との両方を検出するために明度差の絶対値をとると、式ＩＩは式ＩＩＩ：−ｌｎ（閾値）＜｜ｌｎ（ａ）−ｌｎ（ｂ）｜となる。
【００２４】
表ＩＩ
【００２５】
表ＩＩＩは、ｘ方向における表ＩＩの隣接する値の間の差の絶対値、すなわち｜ｌｎ（ａ）−ｌｎ（ｂ）｜を報告している。したがって、表ＩＩＩは、９つの列および９つの行を含む。表ＩＩＩの値を、表の最大値、すなわち２．９１１で除算することにより、１．０００に正規化した。正規化された値は表ＩＶに報告されている。１．７８という−ｌｎ（閾値）の値を、１．７８／２．９１１＝０．６１に正規化した。正規化された閾値を表ＩＶに適用することにより表Ｖを生成した。この表Ｖは、−ｌｎ（閾値）すなわち０．６１を下回る値に対して０を報告し、−ｌｎ（閾値）すなわち０．６１を上回る値に対して１を報告している。
【００２６】
表ＩＩＩ
【００２７】
表ＩＶ
【００２８】
表Ｖ
【００２９】
表ＶＩは、ｙ方向における表ＩＩの隣接する値の差の絶対値、すなわち｜ｌｎ（ａ）−ｌｎ（ｂ）｜を報告している。したがって、表ＶＩは、１０の列および８つの行を含む。表ＶＩの値を、表の最大値、すなわち３．２７５１で除算することにより、１．０００に正規化した。正規化された値は表ＶＩＩに報告されている。１．７８の−ｌｎ（閾値）の値を、１．７８／３．２７５１＝０．５４に正規化した。正規化された閾値を表ＶＩＩに適用することにより表ＶＩＩＩを生成した。この表ＶＩＩＩは、−ｌｎ（閾値）すなわち０．５４を下回る値に対して０を報告し、−ｌｎ（閾値）すなわち０．５４を上回る値に対して１を報告している。
【００３０】
表ＶＩ
【００３１】
表ＶＩＩ
【００３２】
表ＶＩＩＩ
【００３３】
表Ｖを、カーネル
［１ １］
でコンボリューションすることにより、９×１０マトリクスを作成した（表ＩＸ）。ここでは、非ゼロエントリは、ｘ方向における明から暗への遷移または暗から明への遷移を示す。
【００３４】
表ＩＸ
【００３５】
表ＶＩＩＩをカーネル
でコンボリューションすることにより、９×１０マトリクスを作成した（表Ｘ）。ここでは、非ゼロエントリは、ｙ方向における明から暗への遷移または暗から明への遷移を示す。
【００３６】
表Ｘ
【００３７】
表ＩＸおよびＸによって表されるマトリクスを加算した結果、表ＸＩとして報告されるマトリクスとなった。
【００３８】
表ＸＩ
【００３９】
表ＸＩの非ゼロ値のすべてを包含する４つの矩形領域（Ａ１：Ｂ３、Ｄ５：Ｆ７、Ｈ１：Ｉ３、Ｉ９：Ｉ１０）をデコンボリューションのために選択した。選択された領域は、９０の特徴のうちの２３を含み、デコンボリューションの多くの方法によって、計算労力の程度が解析される領域のサイズに応じて指数関数的に上昇するようになるため、画像全体のデコンボリューションに必要となる計算労力の少なくとも約７４％、おそらくははるかにそれより多くが省かれる。
【００４０】
当業者には、本発明の範囲および原理から逸脱することなく本発明のさまざまな変更および改変が明らかとなろう。また、本発明は、上述した例示的な実施形態に不当に限定されるべきではない、ということが理解されるべきである。
【図面の簡単な説明】
【００４１】
【図１】本発明を使用してもよい典型的なスキャニングシステムの概略図である。
【図２】例で使用される対象画像である。
【図３】例で説明するような図２の画像の解析グリッドである。
【図４】図２の第１列第５行の特徴を含む図２の詳細である。
【図５】図４のセグメントのｘ位置に対してプロットされたｙ方向における４画素にわたって積分された画素強度のグラフである。
【特許請求の範囲】
【請求項１】
デコンボリューションのために画像の領域を選択する方法であって、
ａ）複数の特徴を有する画像を提供するステップであって、各特徴が少なくとも１つの値（ｖ）に関連する、ステップと、
ｂ）テスト特徴を識別するステップであって、該テスト特徴が該画像の既知の低値ゾーンに隣接する高値特徴であり、該テスト特徴がテール率（ｒｔ）を有し、該テール率が、該テスト特徴の値（ｖｔ）の該画像の該隣接する低値ゾーンの値（Ｖｏ）に対する比率である、ステップと、
ｃ）閾値ｔを計算するステップであって、該閾値（Ｔ（ｒｔ））が、該テスト特徴のテール率（ｒｔ）の関数である、ステップと、
ｄ）該画像の選択領域を識別するステップであって、該選択領域が画像全体に満たない領域を含み、該選択領域が、隣接する特徴の間の値（ｖ）の比率が該閾値（Ｔ（ｒｔ））を上回る領域である、ステップと、
を含む、方法。
【請求項２】
ステップｂ）が、前記テール率（ｒｔ）を計算する前に、前記テスト特徴の値（ｖｔ）と前記画像の前記隣接する低値ゾーンの前記値（Ｖｏ）との両方から背景定数を減算することをさらに含む、請求項１に記載の方法。
【請求項３】
ｅ）オートグリッド解析により擬似画像を形成するステップ
をさらに含む、請求項１に記載の方法。
【請求項４】
前記閾値（Ｔ（ｒｔ））が、前記テスト特徴のテール率（ｒｔ）の倍数である、請求項１に記載の方法。
【請求項５】
ｆ）前記画像の前記選択領域に対してデコンボリューションを行うステップ
をさらに含む、請求項１に記載の方法。
【請求項６】
デコンボリューションのために画像の領域を選択するシステムであって、
ｂ）デジタル化画像を提供する画像装置と、
ｃ）データ記憶装置と、
ｄ）該画像装置から該デジタル化画像を受け取り、該データ記憶装置に対してデータを書き込みかつ読み出すことができる中央処理装置であって、
ｉ）該画像装置からデジタル化画像を受け取り、
ｉｉ）複数の特徴を識別し各特徴を少なくとも１つの値（ｖ）に関連付け、
ｉｉｉ）テスト特徴を識別し、該テスト特徴が該画像の既知の低値ゾーンに隣接する高値特徴であり、該テスト特徴がテール率（ｒｔ）を有し、該テール率が、該テスト特徴の値（ｖｔ）の該画像の該隣接する低値ゾーンの値（ｖｏ）に対する比率であり、
ｉｖ）閾値ｔを計算し、該閾値（Ｔ（ｒｔ））が該テスト特徴のテール率（ｒｔ）の関数であり、
ｖ）該画像の選択領域を識別し、該選択領域が該画像全体に満たない領域を含み、該選択領域が、隣接する特徴の間の値（ｖ）の比率が該閾値（Ｔ（ｒｔ））を上回る領域である
ようにプログラムされた中央処理装置と、
を備える、システム。
【請求項７】
前記中央処理装置が、前記テール率（ｒｔ）を計算する前に、前記テスト特徴の値（ｖｔ）と前記画像の前記隣接する低値ゾーンの前記値（Ｖｏ）との両方から背景定数を減算するようにさらにプログラムされる、請求項６に記載のシステム。
【請求項８】
前記中央処理装置が、オートグリッド解析により擬似画像を形成するようにさらにプログラムされる、請求項６に記載のシステム。
【請求項９】
前記閾値（Ｔ（ｒｔ））が、前記テスト特徴のテール率（ｒｔ）の倍数である、請求項６に記載のシステム。
【請求項１０】
前記中央処理装置が、前記画像の前記選択領域に対してデコンボリューションを行うようにさらにプログラムされる、請求項６に記載のシステム。
【請求項１】
デコンボリューションのために画像の領域を選択する方法であって、
ａ）複数の特徴を有する画像を提供するステップであって、各特徴が少なくとも１つの値（ｖ）に関連する、ステップと、
ｂ）テスト特徴を識別するステップであって、該テスト特徴が該画像の既知の低値ゾーンに隣接する高値特徴であり、該テスト特徴がテール率（ｒｔ）を有し、該テール率が、該テスト特徴の値（ｖｔ）の該画像の該隣接する低値ゾーンの値（Ｖｏ）に対する比率である、ステップと、
ｃ）閾値ｔを計算するステップであって、該閾値（Ｔ（ｒｔ））が、該テスト特徴のテール率（ｒｔ）の関数である、ステップと、
ｄ）該画像の選択領域を識別するステップであって、該選択領域が画像全体に満たない領域を含み、該選択領域が、隣接する特徴の間の値（ｖ）の比率が該閾値（Ｔ（ｒｔ））を上回る領域である、ステップと、
を含む、方法。
【請求項２】
ステップｂ）が、前記テール率（ｒｔ）を計算する前に、前記テスト特徴の値（ｖｔ）と前記画像の前記隣接する低値ゾーンの前記値（Ｖｏ）との両方から背景定数を減算することをさらに含む、請求項１に記載の方法。
【請求項３】
ｅ）オートグリッド解析により擬似画像を形成するステップ
をさらに含む、請求項１に記載の方法。
【請求項４】
前記閾値（Ｔ（ｒｔ））が、前記テスト特徴のテール率（ｒｔ）の倍数である、請求項１に記載の方法。
【請求項５】
ｆ）前記画像の前記選択領域に対してデコンボリューションを行うステップ
をさらに含む、請求項１に記載の方法。
【請求項６】
デコンボリューションのために画像の領域を選択するシステムであって、
ｂ）デジタル化画像を提供する画像装置と、
ｃ）データ記憶装置と、
ｄ）該画像装置から該デジタル化画像を受け取り、該データ記憶装置に対してデータを書き込みかつ読み出すことができる中央処理装置であって、
ｉ）該画像装置からデジタル化画像を受け取り、
ｉｉ）複数の特徴を識別し各特徴を少なくとも１つの値（ｖ）に関連付け、
ｉｉｉ）テスト特徴を識別し、該テスト特徴が該画像の既知の低値ゾーンに隣接する高値特徴であり、該テスト特徴がテール率（ｒｔ）を有し、該テール率が、該テスト特徴の値（ｖｔ）の該画像の該隣接する低値ゾーンの値（ｖｏ）に対する比率であり、
ｉｖ）閾値ｔを計算し、該閾値（Ｔ（ｒｔ））が該テスト特徴のテール率（ｒｔ）の関数であり、
ｖ）該画像の選択領域を識別し、該選択領域が該画像全体に満たない領域を含み、該選択領域が、隣接する特徴の間の値（ｖ）の比率が該閾値（Ｔ（ｒｔ））を上回る領域である
ようにプログラムされた中央処理装置と、
を備える、システム。
【請求項７】
前記中央処理装置が、前記テール率（ｒｔ）を計算する前に、前記テスト特徴の値（ｖｔ）と前記画像の前記隣接する低値ゾーンの前記値（Ｖｏ）との両方から背景定数を減算するようにさらにプログラムされる、請求項６に記載のシステム。
【請求項８】
前記中央処理装置が、オートグリッド解析により擬似画像を形成するようにさらにプログラムされる、請求項６に記載のシステム。
【請求項９】
前記閾値（Ｔ（ｒｔ））が、前記テスト特徴のテール率（ｒｔ）の倍数である、請求項６に記載のシステム。
【請求項１０】
前記中央処理装置が、前記画像の前記選択領域に対してデコンボリューションを行うようにさらにプログラムされる、請求項６に記載のシステム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図２】
【図３】
【図４】
【図５】
【公表番号】特表２００８−５０１１８７（Ｐ２００８−５０１１８７Ａ）
【公表日】平成２０年１月１７日（２００８．１．１７）
【国際特許分類】
【出願番号】特願２００７−５１５１０６（Ｐ２００７−５１５１０６）
【出願日】平成１７年４月２９日（２００５．４．２９）
【国際出願番号】ＰＣＴ／ＵＳ２００５／０１４８２３
【国際公開番号】ＷＯ２００５／１１９５９３
【国際公開日】平成１７年１２月１５日（２００５．１２．１５）
【出願人】（５９９０５６４３７）スリーエム イノベイティブ プロパティズ カンパニー (1,802)
【Ｆターム（参考）】
【公表日】平成２０年１月１７日（２００８．１．１７）
【国際特許分類】
【出願日】平成１７年４月２９日（２００５．４．２９）
【国際出願番号】ＰＣＴ／ＵＳ２００５／０１４８２３
【国際公開番号】ＷＯ２００５／１１９５９３
【国際公開日】平成１７年１２月１５日（２００５．１２．１５）
【出願人】（５９９０５６４３７）スリーエム イノベイティブ プロパティズ カンパニー (1,802)
【Ｆターム（参考）】
[ Back to top ]