(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】
(43)【公表日】2022-10-13
(54)【発明の名称】光学イメージングパフォーマンス試験システムおよび方法
(51)【国際特許分類】
G02B 21/36 20060101AFI20221005BHJP
G02B 21/06 20060101ALI20221005BHJP
【FI】
G02B21/36
G02B21/06
【審査請求】有
【予備審査請求】未請求
(21)【出願番号】P 2022506600
(86)(22)【出願日】2020-08-07
(85)【翻訳文提出日】2022-03-28
(86)【国際出願番号】 US2020045525
(87)【国際公開番号】W WO2021026515
(87)【国際公開日】2021-02-11
(32)【優先日】2019-08-07
(33)【優先権主張国・地域又は機関】US
(81)【指定国・地域】
(71)【出願人】
【識別番号】399117121
【氏名又は名称】アジレント・テクノロジーズ・インク
【氏名又は名称原語表記】AGILENT TECHNOLOGIES, INC.
(74)【代理人】
【識別番号】110002310
【氏名又は名称】特許業務法人あい特許事務所
(72)【発明者】
【氏名】リッジウェイ,ウィリアム ケイ.
(72)【発明者】
【氏名】ダイク,アラン
(72)【発明者】
【氏名】チェン,ペンユエン
(72)【発明者】
【氏名】クルザヴァ,レベッカ ディ リッコ
【テーマコード(参考)】
2H052
【Fターム(参考)】
2H052AC05
2H052AD06
2H052AF25
(57)【要約】
【課題】光学イメージングシステムのイメージングパフォーマンス試験のための改善されたアプローチを提供する。
【解決手段】試験ターゲットを、光学システムの物体平面に置き、画像ビームを生成するために照射する。試験ターゲットの1つまたは複数の画像を画像ビームから取得する。取得した画像データから、試験ターゲット内の複数の場所でのエッジ広がり関数を演算し、点広がり関数のモデルを、エッジ広がり関数から構築する。点広がり関数に基づいて、複数の場所に対応する複数のイメージングパフォーマンス値を演算する。イメージングパフォーマンス値は、エンスクエアドエネルギー、エンサークルドエネルギー、またはストレールレシオに基づいている。
【選択図】
図1
【特許請求の範囲】
【請求項1】
光学システムのイメージングパフォーマンスを試験する方法であって、
試験ターゲットを前記光学システムの物体平面に位置付けることと、
前記光学システムを動作させて前記試験ターゲットを照明し、画像ビームを生成することと、
前記光学システムの焦点合わせステージを動作させて、複数の焦点ぼけの値に対応する前記画像ビームから前記試験ターゲットの複数の画像を取得することと、
各画像から、前記試験ターゲット内の複数の場所での複数のエッジ広がり関数を演算することと、
それぞれのエッジ広がり関数から複数の点広がり関数を構築することと、
前記点広がり関数に基づいて、前記複数の場所に対応する複数のイメージングパフォーマンス値を演算することであって、前記イメージングパフォーマンス値が、エンスクエアドエネルギー、エンサークルドエネルギー、およびストレールレシオからなる群から選択される計測基準に基づいている、演算することと、を含む、方法。
【請求項2】
(a)前記複数の場所は、前記試験ターゲット内の前記物体平面内の複数のフィールド座標(x,y)を含むことと、
(b)前記複数の場所は、前記試験ターゲットを通過する光軸に沿った複数の焦点位置(z)を含み、前記光学システムを動作させることは、異なる焦点位置(z)での前記試験ターゲットの複数の画像を取得することを含むことと、
のうち少なくとも1つを含む、請求項1に記載の方法。
【請求項3】
前記イメージングパフォーマンス値の組み合わせに基づいてイメージングパフォーマンスの1つまたは複数のマップを生成すること、を含む、請求項1に記載の方法。
【請求項4】
前記マップのうちの2つ以上を比較して、前記物体平面に対する各イメージングデバイスまたはチャネルの焦点面の相対的な位置合わせの尺度を提供すること、を含む、請求項3に記載の方法。
【請求項5】
前記1つまたは複数のマップは、異なるイメージングチャネルに対応し、前記異なるイメージングチャネルは、前記画像を取得するように動作する前記光学システムの異なるイメージングデバイス、または取得される前記画像の異なる波長、または異なるイメージングデバイスと異なる色との両方に対応する、請求項3に記載の方法。
【請求項6】
前記マップのうちの2つ以上を比較して、互いに対する各イメージングデバイスまたはチャネルの相対的な位置合わせの尺度を提供することと、前述の両方と、を含む、請求項5に記載の方法。
【請求項7】
前記1つまたは複数のマップを生成した後、前記光学システムの1つまたは複数の光学部品の位置を調整すること、または前記1つまたは複数のマップによって提供される情報に基づいて前記1つまたは複数の光学部品を交換すること、を含む、請求項3に記載の方法。
【請求項8】
前記1つまたは複数の光学部品は、前記イメージングデバイスのうちの1つまたは複数と、前記光学システムの対物レンズと、前記光学システムの1つまたは複数のチューブレンズと、1つまたは複数のミラーまたはダイクロイックミラーと、前述の2つ以上の組み合わせと、からなる群から選択される、請求項7に記載の方法。
【請求項9】
前記1つまたは複数のマップは、1つまたは複数の初期マップであり、前記1つまたは複数の光学部品を調整または交換した後、前記試験ターゲットの新しい画像を取得し、複数の新しいイメージングパフォーマンス値を演算し、イメージングパフォーマンスの1つまたは複数の新しいマップを生成すること、をさらに含む、請求項7に記載の方法。
【請求項10】
前記1つまたは複数の新しいマップを前記1つまたは複数の初期マップと比較して、イメージングパフォーマンスを最適化するために前記1つまたは複数の光学部品に対して行われる位置調整を決定すること、を含む、請求項9に記載の方法。
【請求項11】
(a)前記調整することまたは交換することで、前記光学システムにおいて共役像面の最適な一対が見つかることと、
(b)前記調整することまたは交換することで、前記イメージングデバイスの焦点マッチングと、イメージングデバイスの傾斜と、像面湾曲の平坦化と、非点収差の低減と、波長に依存する焦点シフトの低減と、前述の2つ以上の組み合わせと、からなる群から選択される属性が改善することと、
のうち少なくとも1つを含む、請求項7に記載の方法。
【請求項12】
前記イメージングパフォーマンス値の組み合わせに基づいてイメージングパフォーマンスの1つまたは複数のグローバルスコアを演算すること、を含む、請求項1に記載の方法。
【請求項13】
(a)前記1つまたは複数のグローバルスコアを変更して、前記試験ターゲット内の前記物体平面内のフィールド座標(x,y)の範囲にわたって、または前記物体平面の焦点位置(z)の範囲を通じて、または前述の両方において、前記イメージングパフォーマンス値の不均一性に対し、ペナルティを課すか、または報酬を与えることと、
(b)前記1つまたは複数のグローバルスコアを変更して、前記試験ターゲット内の前記物体平面内のフィールド座標(x,y)の関数、または前記物体平面の焦点位置(z)の関数、または前述の両方として、異なるイメージングチャネルの類似性に対し、ペナルティを課すか、または報酬を与えること、を含み、前記異なるイメージングチャネルは、前記画像を取得するように動作する前記光学システムの異なるイメージングデバイスに、または取得される前記画像の異なる波長に、または異なるイメージングデバイスと異なる色との両方に、対応することと、
のうち少なくとも1つを含む、請求項12に記載の方法。
【請求項14】
前記イメージングパフォーマンス値を演算した後、前記イメージングパフォーマンス値によって提供される情報に基づいて前記光学システムの1つまたは複数の光学部品の位置を調整すること、を含む、請求項1に記載の方法。
【請求項15】
前記イメージングパフォーマンス値は、初期イメージングパフォーマンス値であり、1つまたは複数の光学部品を調整した後、前記試験ターゲットの新しい画像を取得し、複数の新しいイメージングパフォーマンス値を演算すること、をさらに含む、請求項14に記載の方法。
【請求項16】
前記新しいイメージングパフォーマンス値を前記初期イメージングパフォーマンス値と比較して、イメージングパフォーマンスを最適化するために前記1つまたは複数の光学部品に対して行われる位置調整を決定すること、を含む、請求項15に記載の方法。
【請求項17】
(a)前記試験ターゲットを位置決めすることは、前記光学システムの1つまたは複数の光学部品と共有されるデータムに対して前記ターゲットを位置合わせすること、を含むことと、
(b)前記光学システムを動作させることは、前記画像ビーム内の対物レンズを利用すること、を含み、前記画像ビームの軸に沿って前記対物レンズの位置を調整して、異なる焦点位置(z)での前記試験ターゲットの複数の画像を取得すること、をさらに含むことと、
(c)前記光学システムを動作させることは、前記画像ビーム内の対物レンズを利用すること、を含み、前記画像ビームの軸に沿って前記対物レンズの位置を調整して、異なる焦点位置(z)での前記試験ターゲットの複数の画像を取得すること、をさらに含み、前記対物レンズは、前記対物レンズが無限共役顕微鏡用に構成されていることと、前記対物レンズが有限共役顕微鏡用に構成されていることと、からなる群から選択される構成を有することと、
(d)前記光学システムを動作させることは、前記試験ターゲットのそれぞれの画像を取得するために2つ以上のイメージングデバイスを動作させること、を含むことと、
(e)前記光学システムを動作させることは、前記試験ターゲットのそれぞれの画像を取得するために2つ以上のイメージングデバイスを動作させること、を含み、前記2つ以上のイメージングデバイスは、2つ以上の異なる波長で前記それぞれの画像を取得することと、
(f)前記光学システムを動作させることは、前記試験ターゲットのそれぞれの画像を取得するために2つ以上のイメージングデバイスを動作させること、を含み、前記2つ以上のイメージングデバイスは、2つ以上の異なる波長で前記それぞれの画像を取得し、前記試験ターゲットから伝搬する画像ビームを2つ以上の画像ビーム部分に分割することと、前記2つ以上の画像ビーム部分を前記2つ以上のイメージングデバイスにそれぞれ送信することと、をさらに含むことと、
(g)前記光学システムを動作させることは、選択された波長で前記画像ビームをフィルタリングするためにフィルタアセンブリを動作させること、を含むことと、
(h)前記光学システムを動作させることは、前記画像ビームにおいてチューブレンズを利用すること、を含み、複数の画像を取得するために前記チューブレンズ内の1つもしくは複数のレンズまたはレンズ群の相対位置を調整して、前記チューブレンズの異なる位置での前記試験ターゲットの複数の画像を取得すること、をさらに含むことと、
(i)前記試験ターゲットは、暗い材料と、前記暗い材料上に配置された明るい特性のアレイと、を含むことと、
(j)前記試験ターゲットは、暗い材料と、前記暗い材料上に配置された明るい特性のアレイと、を含み、前記明るい特性が、多角形であることと、
(k)前記試験ターゲットは、暗い材料と、前記暗い材料上に配置された明るい特性のアレイと、を含み、前記明るい特性は、多角形であり、前記明るい特性は、前記明るい特性のエッジが、前記画像を取得する前記光学イメージングシステムのピクセルアレイに対してある角度で配向されるように、傾斜していることと、
のうち少なくとも1つを含む、請求項1に記載の方法。
【請求項18】
光学イメージングパフォーマンス試験システムであって、
試験ターゲットを保持するように構成されたターゲットホルダと、
前記試験ターゲットを照明するように構成された光源と、
前記試験ターゲットの画像を取得するように構成されたイメージングデバイスと、
前記試験ターゲットと前記イメージングデバイスとの間のイメージング光路に位置付けられた対物レンズであって、前記対物レンズまたは前記ターゲットホルダのうち少なくとも1つの位置が、前記イメージング光路に沿って調整可能である、対物レンズと、
電子プロセッサとメモリとを備えたコントローラであって、請求項1に記載の方法の、前記複数のエッジ広がり関数を演算するステップと、前記複数の点広がり関数モデルを構築するステップと、前記複数のイメージングパフォーマンス値を演算するステップと、を制御するように構成された、コントローラと、
を備える、光学イメージングパフォーマンス試験システム。
【請求項19】
(a)前記対物レンズは、前記対物レンズが無限共役顕微鏡用に構成されていることと、前記対物レンズが有限共役顕微鏡用に構成されていることと、からなる群から選択される構成を有することと、
(b)前記イメージングデバイスは、複数のイメージングデバイスを備え、前記イメージングデバイスは、前記イメージング光路をそれぞれ前記イメージングデバイスに向けられた複数のイメージング光路に分割するように構成された画像分離ミラーをさらに備えることと、
(c)前記イメージングデバイスに伝播するために前記イメージング光路内の画像ビームの波長を選択するように構成されたフィルタアセンブリと、
(d)前記イメージング光路に位置付けられたチューブレンズであって、前記チューブレンズ内の1つもしくは複数のレンズまたはレンズ群の相対位置が調整可能である、チューブレンズと、
(e)前記試験ターゲットは、暗い材料と、前記暗い材料上に配置された明るい特性のアレイと、を含む、前記試験ターゲットと、
(f)前記試験ターゲットは、暗い材料と、前記暗い材料上に配置された明るい特性のアレイと、を含み、前記明るい特性は、多角形である、前記試験ターゲットと、
(g)前記試験ターゲットは、暗い材料と、前記暗い材料上に配置された明るい特性のアレイと、を含み、前記明るい特性は、多角形であり、前記明るい特性は、前記明るい特性のエッジが、前記画像を取得する前記光学イメージングシステムのピクセルアレイに対してある角度で配向されるように、傾斜している、前記試験ターゲットと、
のうち少なくとも1つを含む、請求項18に記載のシステム。
【請求項20】
プロセッサ上で実行されると、請求項1に記載の前記方法の、前記複数のエッジ広がり関数を演算するステップと、前記複数の点広がり関数モデルを構築するステップと、前記複数のイメージングパフォーマンス値を演算するステップと、を実行する、非一時的なコンピュータ可読媒体。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、全体的に、分析対象のサンプルからイメージングデータを取得するために顕微鏡検査で利用できるような、光学イメージングシステムのイメージングパフォーマンスの試験に関する。
【背景技術】
【0002】
顕微鏡で利用される自動化されたハイスループット光学イメージングシステムを含む光学イメージングシステムにおいて、得られるイメージングデータの品質の重要な測定基準は、得られる画像の鮮明さである。
【0003】
顕微鏡イメージングの多くのアプリケーション(例えば、DNAシーケンスにおける)は、標本の焦点と傾斜とが不完全であることに関係なく、1つまたは複数のチャネルまたはカメラの視野全体で均一に鮮明な画像を生成する能力に依存している。
【0004】
光学イメージングシステムのカメラおよび他の光学部品相互の、かつイメージングされる標本に対する、光学的位置合わせは、画像の鮮明さ、したがって生成されるイメージングデータの品質に大きく貢献する。
【0005】
したがって、光学イメージングシステムのイメージングパフォーマンス(つまり、品質)を試験および評価できることが望ましい。例えば、光学イメージングシステムが、光学イメージングシステムの所与の用途に許容可能であるとみなされる最低限の品質のある所定のレベルを満たすかまたは超えるかどうかを決定することが有用である。適切な製造および製造パフォーマンス測定基準の演算では、単一の正確な数値を使用して全体のイメージングパフォーマンスを記述することがさらに有用である。
【0006】
光学イメージングシステムによって取得された画像は、実際の標本(または対象物)[S(x,y)]と、顕微鏡において点広がり関数(Point Spread Function)[PSF(x,y)]と呼ばれる機器応答関数と、の混ざり合いとして特徴付けることができる。数学的には、画像I(x,y)は、次のように、標本とPSFとの畳み込みと、各ピクセルのノイズN(x,y)と、の合計である。
【0007】
【0008】
イメージングパフォーマンスの評価には、標本(試験ターゲット、つまり、試験目的で提供される対象物であってもよい)の画像を取得し、そのイメージングデータを用いてPSFを測定および定量することを含んでおり、所定の最小合格基準に対する結果と比較することができる。
【0009】
単一のスカラー数を使用してPSFの品質を記述することが望ましい。実空間PSFの半値全幅(FWHM)もしくはストレールレシオ(Strehl ratio)、単一周波数での波面誤差もしくは変調伝達関数(MTF)コントラスト、および実空間量エンサークルドエネルギー(Encircled Energy)もしくはエンスクエアドエネルギー(Ensquared Energy)などの、いくつかの測定基準が一般的に使用されている。各量の測定の推定の精度および正確さは、それらを演算するために使用される実験方法と強く関連している。それぞれの長所および短所は、それらが推定される精度、解釈の容易さ、ならびに不十分な画像解像度対ノイズ比(SNR)と不十分な画像信号対ノイズ比(SNR)との両方を捕捉する能力に関係する。
【0010】
上記の古典的な式(1)は、PSFが画像のフィールド全体で変化するという事実を反映しておらず、これは、画像が中央で鮮明(高品質)であり、エッジでぼやけている(品質が低い)可能性があるということと同じである。回折限界に近く、弱い光源の検出に最適化された光学システムは、4096×3000ピクセルの画像など、典型的にははるかに広い画像とは対照的に、約5×5ピクセルの領域に密度を含むコンパクトなPSFを備えている。
【0011】
小さなPSFは、画像全体で大幅に変化し得るため、視野全体でのPSFの複数回の独立した測定を行う必要がある。フィールド座標(x,y)の切り替えに加えて、PSFは、z方向の画像の焦点の影響を強く受け、高パフォーマンス、高開口数(NA)の対物レンズでは、焦点がぼけるとPSFのパフォーマンスが急速に低下し、これは、対物レンズの被写界深度が浅いということである。回折効果は、被写界深度に上限を設定するが、補正が不十分であったり、位置合わせが不十分であったりするシステムでは、被写界深度はさらに浅くなる。光学品質の低いシステムは、フィールドのすべての点で鮮明な画像を取得するために数回焦点を合わせる必要があり、システムのスループットが低下し、サンプルの光劣化が発生する可能性があるため、この実際的な制限は重要である。
【0012】
画像解析ソフトウェアによる画像のダウンストリーム処理は、画像の鮮明さに強く影響されうる結果を生み出す。高コントラストの点状物体のぼやけた画像では、信号対雑音比が根本的に低下し、解像度が低下する。例えば、DNAクラスタのぼやけた画像は、ヌクレオチド量の正確な測光を実行する能力を低下させ、観察できる小さな別個のクラスタの数を減らし、隣接するクラスタ間にクロストークを導入する。
【0013】
画像処理アルゴリズムは、周知のデコンボリューション手法を使用したいくつかの明示的なステップ、またはデータの選択的フィルタリングなどの暗黙的なステップのいずれかを実行することにより、ぼやけた画像を部分的に軽減できる。このようなアプローチでは、各画像または画像のバッチのPSF劣化の程度を推定する必要がある場合に演算コストが発生し、PSF劣化の推定自体がノイズを追加する可能性がある。
【0014】
より高速なアプローチでは、PSFの鮮明さが一定であることを想定しているが、実際のPSFと想定したPSFとの差異により、画像処理結果の品質が低下する可能性がある。したがって、入力画像全体のPSFの均一性は、画像処理アルゴリズムの有効性に影響を与える。
【発明の概要】
【発明が解決しようとする課題】
【0015】
既存の光学的位置合わせ方法は、労働集約的であり、光学システムのパフォーマンス全体を記述する迅速なフィードバックを得るのに苦労する。例えば、波面測定の場合、視野内の単一の点を正確かつ迅速に記述し、それら測定の作用がフィールド内の離れた場所にどのように影響するかをオペレータに知らせることを怠るかもしれない。
【0016】
直接オペレータへの高速予測フィードバックを使用することで、光学的位置合わせに関連する製造人件費を削減できる。したがって、システムパフォーマンスの高速なグローバル分析は、0~5秒程度の遅延を伴うフィードバックを提供するのに役立つことになる。
【0017】
さらに、製造歩留まりの分析(例えば、CpKプロセス能力統計を使用)では、各ユニットの品質を連続変数として推定する必要があり、そのような単一の変数を見て、ユニット全体のビルドを記述するのが理想的である。
【0018】
したがって、光学イメージングシステムのイメージングパフォーマンス試験のための改善されたアプローチに対する継続的な必要性が存在する。
【0019】
前述の問題の全部または一部、および/または当業者が観察しうる他の問題に対処するために、本発明は、以下に示す実施態様の例として説明されているように、方法、プロセス、システム、装置、機器、および/またはデバイスを提供する。
【課題を解決するための手段】
【0020】
本発明の一つの実施の形態によれば、光学システムのイメージングパフォーマンスを試験する方法は、試験ターゲットを光学システムの物体平面に位置付けることと、光学システムを動作させて試験ターゲットを照明し、画像ビームを生成することと、光学システムの焦点合わせステージを動作させて、複数の焦点ぼけ値に対応する試験ターゲットの複数の画像を取得することと、各画像から、試験ターゲット内の複数の場所での複数のエッジ広がり関数を演算することと、それぞれのエッジ広がり関数から複数の点広がり関数モデルを構築することと、点広がり関数に基づいて、焦点ぼけ位置の範囲を通して押し出された画像平面によって定義されるイメージングボリューム内の複数の場所に対応する複数のイメージングパフォーマンス値を演算することであって、イメージングパフォーマンス値が、エンスクエアドエネルギー、エンサークルドエネルギー、またはストレールレシオのうち少なくとも1つの計測基準に基づいている、演算することと、を含む、方法。
【0021】
一つの実施の形態では、試験ターゲットを照明することは、試験対象の光学システム以上の開口数(NA)で試験ターゲットを逆照明することを伴う。
【0022】
一つの実施の形態では、この方法は、イメージングボリューム全体にわたるイメージングパフォーマンス値の分布およびマッピングを評価および記述することを含み、さらに、そのような分布およびマップに基づいて要約スコアを演算することを含む。
【0023】
別の実施の形態によれば、光学イメージングパフォーマンス試験システムは、試験ターゲットを保持するように構成されたターゲットホルダと、試験ターゲットを照明するように構成された光源と、試験ターゲットの画像を取得するように構成されたイメージングデバイスと、試験ターゲットとイメージングデバイスとの間のイメージング光路に位置付けられた対物レンズであって、対物レンズの位置が、イメージング光路に沿って調整可能である、対物レンズと、電子プロセッサとメモリとを備えたコントローラであって、イメージングパフォーマンスを試験するための方法の、複数のエッジ広がり関数を演算することと、点広がり関数モデルを構築することと、複数のイメージングパフォーマンス値を演算することと、のステップまたは機能を少なくとも制御または実行するように構成された、コントローラと、を備えている。
【0024】
一つの実施の形態では、光学イメージングパフォーマンス試験システムは、試験ターゲットとイメージングデバイスとの間のイメージング光路に位置付けられたチューブレンズを備えている。
【0025】
一つの実施の形態では、光学イメージングパフォーマンス試験システムは、イメージングデバイスの角度方向および位置を互いに対しておよび対物レンズに対して変更することを可能にするように構成された他のまたは追加の機械的デバイスまたは構成要素を備えている。
【0026】
一つの実施の形態では、コントローラは、イメージングパフォーマンスを試験するための方法の、イメージングボリュームにわたるイメージングパフォーマンス値の分布およびマッピングを評価および記述することと、そのような分布とマップについて要約スコアを演算することと、のステップまたは機能を制御または実行するように構成されている。
【0027】
別の実施の形態によれば、非一時的なコンピュータ可読媒体は、プロセッサ上で実行されると、イメージングパフォーマンスを試験するための方法の、複数のエッジ広がり関数を演算することと、点広がり関数のモデルを構築することと、複数のイメージングパフォーマンス値を演算することと、のステップまたは機能を少なくとも制御または実行する、格納された命令を含む。
【0028】
別の実施の形態によれば、光学システムのイメージングパフォーマンスを試験するためのシステムは、コンピュータ可読記憶媒体を含む。
【0029】
本発明の他のデバイス、装置、システム、方法、特徴および利点は、以下の図面および詳細な説明の検討により当業者に明らかであるか、または明らかになるであろう。すべてのそのような追加のシステム、方法、特徴および利点は、この説明の中に含まれ、実施の形態の範囲内にあり、添付の特許請求の範囲によって保護されることが意図されている。
【0030】
本発明は、以下の図を参照することにより、よりよく理解できる。図面内の構成要素は必ずしも縮尺通りではなく、むしろ、本発明の原理を示すにあたり強調されている。図中において、同様の符号は、異なる図面を通して対応する部分を示す。
【図面の簡単な説明】
【0031】
【
図1】本開示の一つの実施の形態による光学イメージングパフォーマンス試験システムの一例の模式図である。
【
図2】一つの実施の形態による光学イメージングパフォーマンス試験システムのターゲット固定具の一例の模式斜視断面図である。
【
図3A】一つの実施の形態による光学イメージングパフォーマンス試験システムの試験ターゲットの一例の模式平面図である。
【
図3B】
図3Aに示した試験ターゲットのあるセクションの詳細図である。
【
図4A】本明細書で開示されたシステムと方法を実装することによって生成された、真の適度に鮮明な点広がり関数またはPSF(真のエンスクエアドエネルギー値が0.594のPSF)の3次元(3D)モデルの例であり、強度レベルは色分けされており、異なる色は異なるレベルの陰影で表されている。
【
図4B】線広がり関数(LSF×LSF)近似による線広がり関数から導出された、
図4Aに関連する同じPSFの3Dモデルの例である。
【
図4C】PSFモデルがエッジ広がり関数(ESF)データにフィットした、
図4Aに関連する同じPSFの3Dモデルの例である。
【
図5A】真の比較的鮮明なPSF(真のEE値が0.693)の、
図4Aと同様な3Dモデルの別の例である。
【
図5B】LSF×LSF近似から導出された、
図5Aに関連する同じPSFの3Dモデルの例である。
【
図5C】PSFモデルがエッジ広がり関数(ESF)データにフィットした、
図5Aに関連する同じPSFの3Dモデルの例である。
【
図6A】真の比較的不鮮明なPSF(真のEE値が0.349)の、それ以外は
図4Aと同様な3Dモデルの例である。
【
図6B】LSF×LSF近似から導出された、
図6Aに関連する同じPSFの3Dモデルの例である。
【
図6C】PSFモデルがエッジ広がり関数(ESF)データにフィットした、
図6Aに関連する同じPSFの3Dモデルの例である。
【
図7A】焦点合わせによる様々なPSFのLSF×LSF近似の精度のプロットである。
【
図7B】焦点合わせによる様々なPSFのESFデータへのEEフィットの精度のプロットである。
【
図7C】焦点合わせによるあるPSFのEE推定の精度のプロットである。
【
図8A】
図3Aに示した試験ターゲットのターゲット特性のうちの1つから取得した画像の例である。
【
図8B】
図8Aの画像から演算したESFの例である。ESFは、外れ値の除去を使用して整えられている。
【
図9】
図8BのESFを誤差関数の線形結合からなるモデルにフィットさせ、その結果のモデルの係数からEEを演算した結果の例である。
【
図10】ターゲット特性313の4つのエッジについてのESFを演算し、ESFを誤差関数の線形結合にフィットさせ、対応するフィットから4つの独立したEE値を演算する例である。
【
図11】本開示の一つの実施の形態による光学システムのイメージングパフォーマンスを試験する方法の一例を示す流れ図である。
【
図12】本開示の一つの実施の形態による画像のZスタックの取得を示す流れ図である。
【
図13】本開示の一つの実施の形態による個々の画像の処理を示す流れ図である。
【
図14】本開示の一つの実施の形態による、方法に従って生成できる一組のEEパフォーマンスマップの例を示す。
【
図15】本開示の一つの実施の形態による、方法に従って生成できる組み合わせEEのスコアおよびパフォーマンスマップの例を示す。
【
図16】
図1に示されたシステムなどの本開示の実施の形態による光学イメージングパフォーマンス試験システムの一部であるか、またはそれと通信することができるシステムコントローラ(またはコントローラ)の非限定的な例の模式図である。
【発明を実施するための形態】
【0032】
本開示は、光学イメージングシステムのイメージングパフォーマンスを試験するためのシステムおよび方法を提供する。
【0033】
図1は、本開示の一つの実施の形態による光学イメージングパフォーマンス試験システム100の一例を示す模式図である。全体的に、光学イメージングパフォーマンス試験システム100の様々な構成要素の構造および動作は、当業者に知られており、したがって、開示されている主題を理解するために必要な場合に、本明細書ではごく簡単に説明される。説明に役立つために、
図1にはデカルト(XYZ)座標系が含まれており、その原点は描かれている図に対して任意に位置している。
【0034】
全体的に、試験システム100は、試験対象の光学システム104、試験ターゲットステージアセンブリ106、および光学システム104および試験ターゲットステージアセンブリ106の適切な構成要素との(必要に応じて、例えば、電力、データ、測定信号、制御信号等を送信するための)電気的な通信におけるシステムコントローラ108を含む。
【0035】
光学システム104は、商業的入手可能性を目的とした製品など、試験される実際の光学システムであってもよい。あるいは、光学システム104は、試験システム100によって支持された、または試験システム100(例えば、試験システム100の光学ベンチ上)に装着された、実際の光学システムをエミュレートするように配置された、構成要素の集合体であってもよい(例えば、それらの構成要素が、互いに相対的な距離および方向で位置付けられ決定されて、目的の動作の構成に従って互いに位置合わせされた状態で)。
【0036】
全体的に、光学システム104は、分析対象の物体またはサンプルの光学画像を取得するように構成された任意のシステム(または機器、デバイス等)であってもよい。そのようなサンプルは、生物学的なもの(例えば、胞子、真菌、カビ、細菌、ウイルス、生物細胞または例えば核酸(DNA、RNA等)などの細胞内成分、皮膚細胞、有機堆積物などの生物学的に誘導された粒子等)または非生物学的なものであってもよい。
【0037】
典型的には、光学システム104は、ある範囲の倍率で光学画像を取得するように構成されており、したがって、一種の顕微鏡であってもよい(またはその一部であってもよい)。
【0038】
光学システム104は、(実際にイメージングされるサンプルの代わりに)試験ターゲット116を保持するように構成されたターゲットホルダ112(例えば、サンプルステージ、またはサンプルステージのエミュレーション)を含んでもよい、または含まなくてもよい。ターゲットホルダ112が分離している場合、全体的に、正確に位置する方法でターゲットホルダ112の取り付けを可能にする機械的データム(datum)160が設けられてもよい。
【0039】
光源120は、試験ターゲット116を照明するように構成されており、光源120は、試験対象の光学システム104の一部であってもよいし、一部でなくてもよい。カメラなどの1つまたは複数のイメージングデバイス124および128(例えば、図示の実施の形態では第1のイメージングデバイス124および第2のイメージングデバイス128)は、試験ターゲット116を撮像するように構成されている。光学システム104は、当業者によって理解されるように、光源120からターゲットホルダ112および試験ターゲット116まで(および図示の実施の形態では試験ターゲット116を通って)照明光路およびターゲットホルダ112および試験ターゲット116からイメージングデバイス124および128へのイメージング光路を定義するために必要な場合の他の中間光学系(光学部品)をさらに含む。
【0040】
特に、光学システム104は、典型的には、ターゲットホルダ112および試験ターゲット116とイメージングデバイス124および128との間のイメージング光路に位置付けられた対物レンズ132を含む。対物レンズ132は、イメージングデバイス124および128によって撮像される画像を拡大および焦点を合わせるように構成されている。無限遠補正顕微鏡光学系という任意選択の場合では、対物レンズ132は、チューブレンズ152および156と組み合わせてこのタスクを実行する。
【0041】
図示の実施の形態では、対物レンズ132は、対物レンズステージ(または位置決めデバイス)に装着された対物レンズシステムを含む対物レンズアセンブリを模式的に示している。対物レンズステージは、Z軸に沿って対物レンズシステムを移動(平行移動)するように構成されている。したがって、対物レンズ132は、試験目的の必要に応じて、試験ターゲット116の厚さの任意の所望の深さ(高さ)で焦点を合わせるまたは焦点をぼかすように調整可能である。対物レンズ132は異なる高さに焦点を合わせることができるので、調整可能性は、当業者によって理解されるように、光学システム104が画像の「zスタック」を生成することを可能にする。あるいは、焦点を変化させるために、対物レンズ132を空間に固定し、標本をz方向に沿って移動させることができる。あるいは、対物レンズ132および標本を所定の位置に保持することができ、チューブレンズ152および156を、イメージングデバイス124および128に対して移動させることができる。
【0042】
本実施例では、Z軸は、照明光路およびイメージング光路が全体的に存在する試験ターゲット116であって、それに沿って対物レンズ132が焦点合わせおよび焦点ぼけを調整することができる、試験ターゲット116を通る光軸(
図1の視点から垂直)となるようにとられる。したがって、X軸およびY軸は、Z軸に直交する横断面にある。したがって、試験ターゲット116がターゲットホルダ112によって固定されている物体平面上の場所は、(X,Y)座標を使用してマッピングすることができ、深さ(焦点位置)は、Z座標にマッピングすることができる。
【0043】
図1に示す実施の形態では、光学システム104は、対物レンズ132がターゲットホルダ112および試験ターゲット116の上方に位置付けられているトップリード構成を有する。代替の実施の形態では、光学システム104は、対物レンズ132がターゲットホルダ112および試験ターゲット116の下方に位置付けられる、反転またはボトムリード構成を有してもよい。
【0044】
光学システム104は、試験ターゲット116からの透過光または反射光を捕捉するように構成できる。図示の実施の形態では、光学システム104は、照明光路およびイメージング光路がターゲットホルダ112および試験ターゲット116の反対側にあるように、透過照明用に構成されている。代替の実施の形態では、光学システム104は、照明光路およびイメージング光路が、ターゲットホルダ112および試験ターゲット116の同じ側にあるように構成できる。例えば、対物レンズ132は、当業者によって理解されるように、落射照明または落射蛍光構成などにおける、照明光路およびイメージング光路の両方に位置付けることができる。
【0045】
光学システム104は、ターゲットホルダ112がその一部であるターゲット固定具またはアセンブリ136を含んでもよい。ターゲットホルダ112は、ターゲットステージ(または位置決めデバイス)140に装着されるか、取り付けられるか、さもなければ機械的に参照される。正確なデータム160を利用して、光学システム104を試験ターゲットステージアセンブリ106と嵌合させてもよい。ターゲットステージ140は、電動式または非電動式(および自動および/または手動)の作動を介して、公知の方法で正確に動くように構成できる。ターゲットステージ140は、1つまたは複数の軸に沿って平行移動および/またはそれら軸を中心に回転するように構成できる。例えば、ターゲットステージ140は、5自由度(DOF)での移動のために構成されてもよく、その場合、ターゲットステージ140は、Z軸に沿って移動する自由を欠く(焦点を合わせる動きを欠く)。
【0046】
照明光が(図示のように)ターゲットホルダ112の下方から向けられる実施の形態では、ターゲットホルダ112は、照明光の通過を可能にするための窓または開口188を含んでもよく、試験ターゲット116は、そのような窓または開口に装着できる。
【0047】
光源120は、光学顕微鏡法に適した任意のインコヒーレント光源であってもよい。対照的に、コヒーレント光源は、現在の分析方法と互換性のないESFを生成する。インコヒーレント光源120の例には、広帯域光源(例えば、ハロゲンランプ、白熱灯等)、発光ダイオード(LED)、またはLED、レーザ、白熱灯等によって光学的に励起されるリン型材料が含まれるが、これらに限定されない。いくつかの実施の形態では、光源120は、異なる波長で光を出射するように構成された複数の発光ユニット(例えば、LED)、および所与の時間に発光ユニットの選択をアクティブにすることを可能にするように構成された機構を含んでもよい。
【0048】
図示の実施の形態では、光源120は、ターゲット固定具136と統合され(したがって、それと共に移動可能になる)、Z軸に沿ってターゲットホルダ112および試験ターゲット116から一定の距離に位置付けられている。他の実施の形態では、光源120は、ターゲット固定具136から分離されていてもよい。
【0049】
いくつかの実施の形態では、当業者によって理解されるように、コンデンサ(図示せず)を、光源120とターゲットホルダ112および試験ターゲット116との間に配置して、光源120からの照明光を集中させて、試験ターゲット116の照明を増強させることができる。照明光線の角度、特に照明光線の開口数(NA)と放射強度の均一性は、測定されたPSF(Point Spread Function)の精度に影響を与える。理想的な光源のNAは、光学システムの対物NAをわずかに上回り、上記の角度全体で均一な放射強度を示す。
【0050】
イメージングデバイス124および128は、センサがピクセルアレイとなっているマルチピクセル(またはピクセル化)イメージングデバイスである(例えば、電荷結合デバイス(CCD)またはアクティブピクセルセンサ(APS)技術に基づく)。したがって、イメージングデバイス124および128は、顕微鏡検査および他のイメージング用途で一般的に利用されるタイプのカメラであってもよい。
【0051】
図示の実施の形態では、イメージングデバイス124および128のそれぞれは、イメージングデバイスステージ(または位置決めデバイス)に装着されたイメージングデバイスユニット(センサ、つまりピクセルアレイを含む)を含むイメージングデバイスアセンブリを模式的に示す。各イメージングデバイスステージは、6つのDOF移動、つまり3つすべての(X,Y,およびZ)軸に沿った平行移動とそれら軸を中心とした回転の両方に対して構成できる。
【0052】
イメージングデバイス124および128は、互いに対して任意の公称角度および距離で位置付けることができる。図示の実施の形態では、第1のイメージングデバイス124は、対物レンズ132と共に軸上に(Z軸に沿って)位置付けられており、第2のイメージングデバイス128は、Z軸(図示の実施の形態ではX軸上)に対して角度を付けられた(図示の実施の形態では90度)軸上に位置付けられている。
【0053】
光学システム104は、対物レンズ132と第1のイメージングデバイス124との間の軸上に位置付けられた画像分離ミラー144(例えば、ビームスプリッタまたはダイクロイックミラー)を含む。画像分離ミラー144は、試験ターゲット116から伝播する画像ビームを、画像分離ミラー144から(画像分離ミラー144を介した透過を介して)第1のイメージングデバイス124への第1の分割イメージング光路に沿って伝播する第1の画像ビーム部分と、画像分離ミラー144から(画像分離ミラー144からの反射を介して)第2のイメージングデバイス128への第2の分割イメージング光路に沿って伝播する第2の画像ビーム部分と、に分割するように構成されている。画像ビームの第1の画像ビーム部分と第2の画像ビーム部分とへの分離または分割は、動作の電磁スペクトル範囲内の波長(つまり、色)に基づいてもよい。例えば、第1の画像ビーム部分は、画像ビームの青色部分であってもよく、第2の画像ビーム部分は、赤色部分であってもよい。
【0054】
光学システム104は、対物レンズ132とイメージングデバイス124および128との間、または具体的には図示の実施の形態では対物レンズ132と画像分離ミラー144との間のイメージング光路内に位置付けられたフィルタアセンブリ148(つまり、波長セレクタ)を含むことができる。
【0055】
図示の実施の形態では、フィルタアセンブリ148は、
図1の両方向の水平矢印によって示されるように、線形平行移動(例えば、フィルタスライド)または回転(例えば、フィルタホイール)によって、イメージング光路へのおよびイメージング光路からの異なるフィルタを選択的に切り替えるように構成されたフィルタステージに装着された複数のフィルタを示すことができ、それによって、イメージングデバイス124および128によってそれぞれ受光される画像ビーム部分において通過される波長の選択を可能にしてもよい。
【0056】
フィルタアセンブリ148は、各チャネルが異なる波長(つまり、色)に関連付けられ、各イメージングデバイス124および128が1つまたは複数のチャネルに利用される、マルチチャネルイメージングを実施するために利用できる。
【0057】
ある非排他的な例として、光学システム104は、第1のイメージングデバイス124が、2つのチャネル(例えば、青および緑)について指定され、第2のイメージングデバイス128が他の2つのチャネル(例えば、赤および琥珀色)について指定された、4つのチャネル(例えば、青、緑、赤、および琥珀色)で画像を取得するように構成できる。フィルタアセンブリ148の切り替え(位置決め)は、4つのチャネルで画像を順次取得するために必要に応じて、イメージングデバイス124および128ならびに光学システム104の他の構成要素の動作と連係できる。ある特定の例では、4チャネル構成は、異なる蛍光標識が分析対象のサンプル中の異なる核酸塩基(アデニン、シトシン、グアニン、チミン)に関連付けられているDNA配列決定において有用である。
【0058】
試験対象の実施の形態に応じて、光学システム104は、照明光路および/またはイメージング光路に1つまたは複数の他のタイプの光学部品をさらに含むことができる。例には、フィールドレンズ、リレーレンズ、ビームエキスパンダ、ビームコリメータ、開口、スリット、ピンホール、共焦点ディスク等が含まれるが、これらに限定されない。
【0059】
図示の実施の形態では、光学システム104は、画像分離ミラー144と第1のイメージングデバイス124との間の第1の分割イメージング光路に位置付けられた第1のチューブレンズアセンブリ152と、画像分離ミラー144と第2のイメージングデバイス128との間の第2の分割イメージング光路に位置付けられた第2のチューブレンズアセンブリ156と、を含む。
【0060】
図示の実施の形態では、第1のチューブレンズアセンブリ152は、
図1の両方向の垂直矢印で示されているように、第1の分割イメージング光路の軸に沿った第1のチューブレンズ内の要素の相対位置を調整(平行移動)するように構成されたチューブレンズステージに装着された第1のチューブレンズを示してもよい。第2のチューブレンズアセンブリ156は、
図1の両方向の水平矢印で示されているように、第2の分割イメージング光路の軸に沿った第2のチューブレンズ内の要素の相対位置を調整(平行移動)するように構成されたチューブレンズステージに装着された第2のチューブレンズを示してもよい。当業者によって理解されるように、調整可能なチューブレンズを利用して、無限遠補正光学系を実現できる。
【0061】
光学システム104は、ターゲットホルダ112/試験ターゲット116、対物レンズ132等などの、光学システム104の様々な調整可能な構成要素の位置を定義するための基準データムとして機能するように構成された適切な基準データムを含んでもよい。
【0062】
図示の実施の形態では、光学システム104の様々な構成要素が収容されているエンクロージャ160の表面などの光学システム104の固定位置構造が、基準データムとして機能してもよい。
【0063】
全体的に、システムコントローラ(またはコントローラ)108は、光学システム104の様々な構成要素の動作を制御するように構成される。この制御は、光学システム104の調整可能な構成要素の位置を監視および制御(調整)することと、調整可能な構成要素の動作を光源120およびイメージングデバイス124および128などの他の構成要素と調整または同期させることと、を含む。
【0064】
システムコントローラ108は、イメージングデバイス124および128によって出力されたイメージングデータを処理する(例えば、画像をデジタル化および記録/保存すること、コンピュータスクリーンなどのディスプレイデバイス上に表示するための画像をフォーマットすること等を含むデータ取得および信号分析)ようにさらに構成されている。
【0065】
システムコントローラ108は、本明細書に開示される方法であって、その方法に関連付けられた任意のアルゴリズムを実行することを含む、方法に従って光学システム104のパフォーマンス試験を実行するようにさらに構成されている。
【0066】
システムコントローラ108は、光学システム104の動作およびパフォーマンス試験の実施態様を容易にするために必要に応じてユーザインターフェースを提供するようにさらに構成されている。そのようなすべての機能について、システムコントローラ108は、光学システム104と通信する電子モジュール164と、電子モジュール164と通信するコンピューティングデバイス168と、を含むものとして
図1に模式的に示されている。電子モジュール164およびコンピューティングデバイス168は、当業者に理解されるように、そのようなすべての機能を実行するのに適切なすべてのハードウェア(マイクロプロセッサ、メモリ、電気回路、周辺機器等)、ファームウェア、およびソフトウェアコンポーネントを示すと理解される。
【0067】
一つの実施の形態では、電子機器モジュール164は、光学システム104を制御することとパフォーマンス試験を実行することとに専用の構成要素を示してもよく、一方、コンピューティングデバイス168は、光学システム104を制御することとパフォーマンス試験を実行することとを容易にするために電子機器モジュール164とインターフェースするように構成可能またはプログラム可能なより汎用のコンピューティングデバイスを意味してもよい。
【0068】
システムコントローラ108(電子モジュール164および/またはコンピューティングデバイス168)は、本明細書に開示される方法を実行するための非一時的な命令を含む非一時的なコンピュータ可読媒体を含んでもよい。システムコントローラ108(電子モジュール164および/またはコンピューティングデバイス168)は、全体的な制御を提供する主電子プロセッサと、専用の制御動作または特定の信号処理タスクのために構成された1つまたは複数の電子プロセッサと、を含んでもよい。
【0069】
システムコントローラ108(典型的には、少なくともコンピューティングデバイス168)は、1つまたは複数のユーザ入力デバイス(例えば、キーパッド、タッチスクリーン、マウスなど)、ユーザ出力デバイス(例えば、ディスプレイスクリーン、プリンタ、視認インジケータもしくはアラート、可聴インジケータもしくはアラートなど)、ソフトウェアによって制御されるグラフィカルユーザインターフェース(GUI)、および電子プロセッサによって読み取り可能なメディアをロードするためのデバイス(例えば、ソフトウェア、データなどで具体化される非一時的な論理命令)を含んでもよい。システムコントローラ108(典型的には、少なくともコンピューティングデバイス168)は、システムコントローラ108の様々な機能を制御および管理するためのオペレーティングシステム(例えば、マイクロソフトウィンドウズ(登録商標)ソフトウェア)を含んでもよい。
【0070】
図1が、本明細書に開示される光学イメージングパフォーマンス試験システム100の高レベルの模式図であることが理解されよう。当業者によって理解されるように、追加の構造、デバイス、および電子機器などの他の構成要素が、所与の用途のために試験システム100がどのように構成されるかに応じて、実際の実施態様に必要に応じて含まれてもよい。
【0071】
図2は、一つの実施の形態によるターゲット固定具136の例を示す模式斜視断面図である。ターゲットホルダ112は、ターゲット支持ブロック172、クランププレート176、および1つまたは複数のファスナ180を含んでいる。光源120(図示の例ではLED)がターゲット固定具136と統合される図示の実施の形態では、ターゲット固定具136は、光源120を支持および包含する(および任意選択で光源120のヒートシンクとして機能する)ための装着部材184を含んでもよい。装着部材184は、ターゲット支持ブロック172の一部であるか、またはそれに取り付けられていてもよい。ターゲット支持ブロック172は、試験ターゲット116と、試験ターゲット116と光源120との間に提供される任意の適切な光学系(例えば、窓188)とを支持するための環状ショルダを含む。
【0072】
試験ターゲット116をターゲット支持ブロック172に装着した後、カバーガラス(またはカバースリップ)192を試験ターゲット116上に配置してもよい。一例として、カバーガラス192は、500±5μmの厚さを有してもよい。次に、クランププレート176(ばねとして構成できる)は、ファスナ180(例えば、ターゲット支持ブロック172のねじ穴に係合するねじであってもよい)を使用してターゲット支持ブロック172に固定されて、試験ターゲット116およびカバーガラス192をターゲット固定具136内の固定位置に固定する。
【0073】
図3Aは、一つの実施の形態による試験ターゲット116の模式平面図である。試験ターゲット116は、暗い(光吸収)材料309と、暗い材料309上に配置された(またはそれと統合された)明るい(光反射)ターゲット特性313の2次元アレイ(パターン)と、で構成される(または少なくともそれらで構成される外面を含む)平面ジオメトリ(例えば、板状)の基板であってもよいか、その基板を含んでもよい。
【0074】
本明細書では、「暗い」と「明るい」という用語は互いに関連している。すなわち、反射光照明を使用する実施態様の場合、暗い材料は、ターゲット特性313よりも暗く(より光を吸収し)、ターゲット特性313は、暗い材料309よりも明るい(より光を反射する)。透過光照明を特徴とする実施態様の場合、暗い材料は、より光を吸収し、明るい材料は、光をよりよく透過する(より多くの光を透過する)。
【0075】
ターゲット特性313は、好ましくは多角形であり、したがって、明るい特性313と周囲の暗い材料309との間のコントラストを示すエッジを提供している。特に、ターゲット特性313は、直線的(正方形または長方形)であってもよい。
【0076】
一つの実施の形態では、図示のように、ターゲット特性313は、水平方向および垂直方向に対して、特に、イメージングデバイス124および128のピクセルアレイに対して、試験ターゲット116の平面(光軸に鉛直な平面)において傾斜(回転)している。これは、
図8Aにさらに示されており、これは、イメージングデバイス124および128のいずれかのピクセルアレイ上で実現された、試験ターゲット116のあるセクションの詳細図である。
【0077】
非排他的な一例では、暗い材料309は、溶融シリカ基板であり(またはその一部であり)、ターゲット特性313は、基板上に配置されたパターン化されたクロム層である。ターゲット特性313のアレイは、42μmのエッジを有する14×10の正方形のアレイである。ターゲット特性313は、画像の長軸上で92μm、画像の短軸上で94μm離れて(中心から中心まで)離間している。ターゲット特性313は、光軸に対して100mRad(約5.73°)傾斜している。画像は、暗視野に対して19%明るい。クロムは、対物レンズの方を指している。
【0078】
試験システム100は、光学システム104のイメージングパフォーマンスを試験するための方法を実行するために利用できる。この方法によれば、試験ターゲット116は、試験ターゲット116をターゲットホルダ112に固定することによって、光学システム104の物体平面に位置付けられている。様々な調整可能な(段階的な)光学部品(例えば、ターゲットホルダ112、対物レンズ132、チューブレンズ152および156、イメージングデバイス124および128等)の位置決めに対する初期調整がなされてもよい。さらに、イメージングデバイス124および128に渡される波長は、フィルタアセンブリ148を調整することによって選択できる。
【0079】
次に、光源120がアクティブにされて、試験ターゲット116を照明するための照明光を生成し、その結果、試験ターゲット116から出射されるイメージング光がもたらされる。イメージング光は、対物レンズ132ならびにチューブレンズ152および156によって、(画像分離ミラー144によって分割される)画像ビーム部分として、それぞれのイメージングデバイス124および128のピクセルアレイ上に焦点が合わされて、試験ターゲット116の画像(対物レンズ132の視野内)をそれぞれの画像ビーム部分から取得している。次に、イメージングデータは、処理のためにイメージングデバイス124および128からシステムコントローラ108に転送される。
【0080】
この方法によれば、イメージングデータを利用して、試験ターゲット116内の複数の場所での複数のエッジ広がり関数を演算する。次に、演算されたそれぞれのエッジ広がり関数から点広がり関数のモデルが構築される。次に、これらの点広がり関数から、エンスクエアドエネルギーEEの複数の(推定)値が演算される。あるいはまたはさらに、エンサークルドエネルギーおよび/またはストレールレシオが演算できる。
【0081】
傾斜エッジターゲットからの鮮明さ(エンスクエアドエネルギーを含む)の推定
一つの実施の形態では、この方法は、点広がり関数のモデルを生成するためにフィットパラメータを利用できるように選択された基底系を使用して実験的エッジ広がり関数(Edge Spread Function; ESF)をフィットさせることを伴う。そのモデルから、鮮明さの様々な尺度を導き出すことができる。
【0082】
鮮明さ測定基準の定義
エンスクエアドエネルギー
点(x=0,y=0)を中心とする所与のPSF(x,y)について、エンスクエアドエネルギーは、次のように、同じ点を中心とする面積aの正方形内に収まるPSFの割合である。
【0083】
【0084】
エンサークルドエネルギー
点(x=0,y=0)を中心とする特定のPSF(x,y)について、エンサークルドエネルギーは、次のように、同じ点を中心とする半径rの円内に収まるPSFの割合である。
【0085】
【0086】
ストレールレシオ
所与のPSFについて、ストレールレシオは、正規化されたPSFの最大値(ここでは(x=0;y=0)で発生すると仮定)と完全な回折限界PSFの最大値との比である。円形瞳孔の場合、回折限界PSFは、エアリーディスクであり、そのストレールレシオは次のとおりである。
【0087】
【0088】
【0089】
線広がり関数(Line Spread Function; LSF)およびエッジ広がり関数(Edge Spread Function; ESF)
PSFを測定する方法の1つは、バックライトされた穴のアレイを点光源としてイメージングすることであり、穴の直径は、エアリーディスクよりも大幅に小さく、穴間の間隔は隣接する点光源からのPSFが実体のない程度に重なるほどに十分に大きい。そのようなアプローチの問題は、高ODターゲット背景に対するそのような小さな特性の正確な製造(特に高NA顕微鏡を試験する際に)と、PSFが1/r2としての2次元において減衰し、正規信号がすぐにノイズに陥るという事実と、を伴う。
【0090】
ダイナミックレンジが70dBを超えかつウェルの深さが10,000e-1を超える現在利用可能なセンサでも、PSFの中心から数ピクセル離れたPSF密度を特定できない可能性があり、その領域の密度はシステムの信号/ノイズ比(SNR)を急速に低下させることができる。代替の実験的なプロトコルには、量子ドットまたは単一フルオロフォアの希薄な単分散コーティングの蛍光照明が伴い、これにより、信頼性の高い点光源を生成できるが、依然としてホールターゲットと同じSNR制約の犠牲に陥る。
【0091】
一般的な解決策は、ピンホールの代わりに拡張エッジターゲットを使用することであり、これは、変調伝達関数(MTF)測定値を生成するために一般的に使用される方法のいわゆる傾斜エッジファミリである。そのようなターゲットは、線対称性を利用して、次のように線広がり関数(LSF)として知られる点広がり関数の1次元投影を調べる。
【0092】
【0093】
【0094】
LSFは1/rで減衰するため、LSFの解析では、中心から数ピクセル離れたPSFの拡張密度をより適切に捉えることができる。ピンホールを使用してPSFを測定するのではなく、エッジを使用してLSFを測定し、結果として得られる実験的なプロファイルは、エッジ広がり関数(ESF)として知られる。エッジは、イメージングシステムの軸に揃えるか、傾斜させることができる。
【0095】
まず、y軸に位置合わせされ、次の伝達関数E(x,y)で記述されるエッジについて考える。
【0096】
【0097】
このイメージは、PSFとE(x,y)の2次元畳み込み、または同等にLSF(x)とE(x)の1次元畳み込みである。
【0098】
【0099】
【0100】
【0101】
【0102】
概念的には、次の関係に従って、観測されたESF(x)の導関数としてLSF(x)を見つけるのは簡単である。
【0103】
【0104】
しかしながら、実際には、ノイズの多いデータを数値微分すると、望ましくないレベルのノイズが発生する。代わりに、本開示の一態様によれば、数式12は、以下に説明するように直接フィットされる。
【0105】
分離可能な関数を用いたPSFの近似
PSFではなくLSFを測定する方が安定しているが、最終的な目標はPSFについて知ることである。直交する方向に沿って記録された2つの1次元LSFから完全な2次元PSFを再構築することは形式上不可能である。代わりに、本開示の一態様によれば、この方法(例えば、制御デバイスによって実行されるソフトウェアで具体化できるアルゴリズム)は、次のようにPSFの近似を演算する。
【0106】
【0107】
【0108】
この仮定は、次の2つの場合に受け入れられない。(1)例えば、エアリーディスク内のヌルの場所を調べるなど、適切に修正されたPSFの非常に正確なモデルが必要である場合、および(2)評価対象のPSFが、完全にはほど遠く(>1/4波RMS波面誤差)、|y|≒|x|が存在するフィールドの対角コーナで発生するコマ収差などのy軸またはx軸に対しても実質的な非対称性が含まれている場合。
【0109】
LSF×LSF近似の有限精度を
図7Aで評価する(焦点合わせによる様々なPSFに対するLSF×LSF近似の精度)。分析のために、いくつかの公知のPSFのまとまりを水平および垂直LSFに投影し、数式14を使用してそれらの外積から近似PSFを再演算した。
【0110】
エンスクエアドエネルギー(EE4)を、2×2ピクセル領域にわたる元のPSFと近似PSFについて演算した。元のPSFは、高NA顕微鏡システムのいくつかのフィールド点の実験的な波面誤差測定から得られた実世界のPSFであるとした。9つの波面誤差測定のそれぞれについて、瞳孔全体で積分して画像空間PSFを取得する前に、焦点ぼけ項を追加することによって31個のPSFを生成した。より鮮明なPSFについて、LSF×LSF近似は、実際のPSFのEEを一貫してほぼ0:05で過小評価するが、多項式フィットまたは他の内挿法を使用した経験的補正を実行して精度を回復できるのに正確で十分な方法でそのように過小評価する。より不鮮明なPSFについて、LSF×LSF近似と真のPSFとの間の平均対応は良好であるが、精度はEE<0.1の場合の平均よりも低い。不鮮明なPSFの相対精度は低くなるが、不鮮明なPSFは、そのように明確に識別され、鮮明なPSFと混同されない。
【0111】
全体として、近似は、試験対象のシステムとPSF鮮明さ仕様との有意義な比較を可能にするのに十分正確である。
【0112】
ESFフィットプロセスの詳細を説明する前に、ピクセル化された入力データの形式を説明し、画像レジストレーションに関連する実際的な懸念事項に触れることは役立つ。
【0113】
傾斜エッジターゲット
傾斜エッジターゲットの画像は、この方法のPSFフィッティングの唯一のデータソースであり、エッジはそれぞれ4つのエッジを持つ正方形のアレイに配置されている。傾斜エッジESFターゲットの使用は十分に確立されている。単一のエッジについてのESFを考えると、エッジを中心とし、エッジの両側の任意の数のピクセルに進むデータ領域を抽出することができ、その長さは、適度に不鮮明なPSFのESF全体を撮像するのに十分な長さであり、また、画像の完全に暗い部分および完全に明るい部分のベースライン強度レベルの良好な推定値を提供する。傾斜エッジターゲットでは、複数のESFを組み合わせることができ、例えば、エッジがほとんど垂直であるが、1/aラジアンの傾斜がある場合(aは約10の値)、画像の隣接する行のまとまりを、各エッジが、隣接する行に対して1/aピクセルだけシフトされたESFの一意の複製を含有するように、収集してプールすることができる。次に、結果として得られるESFを一緒に表示し、xで(行の高さ)/aだけシフトして、各ESF遷移の中点を揃えることができる。
【0114】
このデータセットは、1/aピクセルの間隔でESFをサンプリングするが、特定のピクセルでの各観測は依然としてピクセルのアクティブ領域全体の平均であるため、それ自体は解像度を上げない。使用されている特定のセンサでは、ピッチ以外のピクセルのアクティブ領域およびジオメトリがわかっていなくてもよい。ピクセルのアクティブ領域は、100%のフィルファクタであることを想定している。したがって、垂直に傾斜したESFの周囲の観測画像は、次のようにモデル化される。
【0115】
【0116】
現在の慣例において、a=10の傾斜エッジターゲットの32ピクセル(エッジの両側16ピクセル)にまたがる64行の画像を、分析のために一緒にプールした。外れ値の除去は、l2ノルムとして演算された距離と、一緒にすべての行の中央値からの距離に基づいて行の4分の1を取り除いた。
【0117】
画像レジストレーション
粗い画像レジストレーション
エッジを、最初に、垂直方向と水平方向の両方での画像の中心の投影の単一周波数フーリエ変換を演算することによって識別した。そのフーリエ変換の位相を、公称倍率を想定して、ターゲットの既知の寸法を使用することにより、画像の全体的なシフトを推定するために使用した。
【0118】
その画像を、隣接するタイルに分割し、各タイルの平均値とともに、各タイルのxとyの各タイルの重心を演算した。FFTの粗いシフトから、画像の各正方形のおおよその場所が十分にわかっており、次いで事前に演算された重心値から各正方形の重心を演算し、各正方形の場所の早期の推定を提供した。
【0119】
正方形(x;y)座標のグローバルセットを、回転、平行移動、拡大の線形モデルにフィットさせ、グローバルモデルに基づいて正方形の場所の修正された推定値を生成した。
【0120】
細かい画像レジストレーション
各正方形を、四つのエッジの中心を特定し、場所情報を平均化することで、正確に位置付けた。垂直エッジへのフィットを、各行に沿って数値微分を演算し、その結果を3乗してエッジを識別し、次いでその積を重心演算の重みとして使用して、エッジが行座標に沿ってどこにあるかを推定することによって演算した。エッジ座標と行番号のまとまりを、中央値のロバスト統計中央値アプローチを使用して、勾配とオフセットの線形モデルにフィットさせた。この手法を、行と列を水平エッジについて交換して、4つのエッジすべてに対して繰り返した。エッジオフセットのまとまりを使用して正方形の中心を推定し、勾配の中央値を使用して、傾斜エッジの演算に使用する正方形の回転角度を演算した。
【0121】
次に、正方形の細かくレジストレーションされたフィット値を、画像全体の平行移動、回転、および拡大のグローバルモデルに再度フィットさせた。この結果を適用して、マルチカメラシステム内の2台のカメラの相対的な向き、具体的には相対的な画像の平行移動Δx;Δyおよび回転θを演算した。拡大パラメータを使用して、各画像の全体的な拡大を演算した。歪みは、特定の歪みモデルを含むようにフィット関数を変更するか、または特定の歪みモデル(例えば、歪みなし)を指定するモデルへのフィットの残差を集合的に分析することによって演算することができる。
【0122】
ESFデータへのPSFモデルのフィッティング
本開示の実施の形態によれば、いくつかの方法を利用して、ESFデータを取得し、次のように、LSF×LSF近似(数式14)と一致するPSFを作成することができる。
【0123】
外積(直積)
LSFx(x)およびLSFy(y)は、垂直および水平ESFデータの数値微分(数式13)から直接演算することができ、その外積を取ることができる。その結果として得られる行列は、エッジが1/aラジアンだけ傾斜している場合、PSFを係数aでオーバーサンプリングし、積分によって画像をダウンサンプリングして、PSFのピクセル化されたバージョンを生成することができる。
【0124】
しかしながら、別の実施の形態では、このアプローチは、いくつかの理由で利用されない。第一に、それはデータの多くの平均化を提供せず、全体的に高い数値微分のノイズに敏感である。第二に、これはピクセル化されたPSFの推定値であり、設計仕様および試験では、ピクセル化されていないPSFを一般的な光学モデリングソフトウェアの出力として参照することがよくある。
【0125】
ESFの補間
数値微分のノイズは、ESFを公知の導関数を使用してN個の滑らかな連続関数にフィットさせることで減らすことができる。ここで、Nは、実験的なESFの合計データ点数よりも小さい。
【0126】
【0127】
【0128】
【0129】
これは、ノイズに役立つが、前の方法のように、ピクセル化されていないPSFには直接フィットしない。
【0130】
PSFモデルへのESFのフィッティング
ESFを滑らかな関数にフィットさせる上での改善点は、最初にPSFaで表されるPSFの特定のモデルを想定し、次いでそのモデルを使用してESFの分析式を導出することである。これの鍵となるのは、PSFモデルがxとyで分離可能な関数であるということである。ガウス基底系は、対数間隔の逆幅biを使用して選択されるが、原則として他の基底系を使用することもできる。
【0131】
【0132】
【0133】
【0134】
【0135】
数式中、(x,y)は、試験ターゲット116内の物体平面上の点のフィールド座標であり、Ngは、PSFをモデル化するために利用される基底系を含むように用いられるガウス分布の数であり、aiは、数式20のモデルにおける各ガウス分布iの寄与を重み付けする係数の組であり、biは、各ガウス分布の幅の逆二乗に比例する。
【0136】
ピクセル化されていないESFは、係数の観点から直接演算することができる。
【0137】
【0138】
【0139】
【0140】
数式中、erf(x)は、ガウス誤差関数、つまり各ガウス分布の定積分であり、数値ライブラリで一般的に利用可能である。
【0141】
ピクセル化されたESFは、次のように、ピクセル化されていないESFの数値積分によって演算することができる。
【0142】
【0143】
【0144】
なお、LSFをフィットするために使用されるフィッティング関数は、LSFx(x-x0)=LSFx(x0-x)となるように、x0について対称である。この対称性は、厳密には必要ではないが、フィット空間を減らし、上昇ESFと下降ESFとの対比、つまり正方形の右側と正方形の左側との対比からのESFの同等性を高めるために使用した。高ゲインセンサからの画像はショットノイズが制限されたカウントプロセスであり、絶対的には、ESFデータの下部の暗い部分の分散は、頂部でのESFデータの分散よりも小さいため、ESFの上部フィット領域と下部フィット領域は、統計的にかなり異なる。ESFの重み付き線形最小二乗法が各ESFデータ点の信号レベル(光電子の単位)を、その分散に相当する重みとして使用して実行される場合、これにより、原点について非対称であるLSFを伴うPSFについての左右のESFへのフィットが異なる。
【0145】
代替のアルゴリズムは、ここで使用されている明白な対称性を無視する数式20における異なる基底系を利用することができる。このような基底系は、Bスプライン、多項式、または、非対称LSFのフィットを可能にする、様々な周波数の奇数乗または正弦関数の多項式と一緒に数式20で使用されるガウス分布などの組み合わせなどの、補間に一般的に使用される関数とすることができる。上昇および下降ESF(例えば、正方形の左エッジおよび右エッジ)の同時線形最小二乗法では、ショットノイズから導出された重みを使用して、非対称LSFを非対称基底関数にフィッティングすることができる。
【0146】
数値フィッティング
実際には、画像データは、線形最小二乗法を使用してフィッティングされる。QR分解を、インテルMKLマスカーネルライブラリ(INTEL MKL Math kernel library)に実装されている公知のDGELS LAPACK(線形代数パッケージ)ルーチンの一部として使用する。画像データImgを、計画行列Dを使用してフィッティングし、フィット係数Aiを得る。9つの異なるオフセットx0を使用して、Dのいくつかの事前演算されたバージョンが存在し、データのすべての行は、データのx0値に最も近いx0値に対応するDの行を使用する。
【0147】
【0148】
【0149】
【0150】
【0151】
【0152】
フィット係数X[0]は、バックグラウンドノイズと試験ターゲットのブリードスルーによる画像のベースライン強度を示す。フィットであるデータがESFx(x)となる垂直エッジである場合、係数X[1];X[2];:::;X[Ng]は、数式20において係数ax;1,ax;1,:::,ax;Ngを示す。 フィットであるデータがESFy(y)となる水平エッジである場合、係数X[1];X[2];:::;X[Ng]は、数式20において係数ay;1,ay;1,:::,ay;Ngを示す。
【0153】
PSFモデルの再構築と鮮明さの演算
PSFのモデルは、前セクションで得たフィット係数aiから再構築することができる。いくつかの異なるモデルは、異なるレベルの知識を示す。最も完全なモデルは、数式20を使用して垂直LSFと水平LSFとを組み合わせたものであり、エンスクエアドエネルギーは、係数の合計の単純な積から直接演算することができる。
【0154】
【0155】
【0156】
この方法の精度については、以下に説明する。
【0157】
エンサークルドエネルギーも同様に数値的に数式3および20を評価することで演算することができる。ストレールレシオは、次のように、単純な形を取る。
【0158】
【0159】
代替シンプルPSFモデル
PSFは基本的に2次元の概念であるが、垂直LSFと水平LSFの両方が同一であると仮定して決定されたPSFのEEを表現することは有用であった。そのような用途には、水平ESFから推定されたEEと垂直ESFから推定されたEEとの対比を最大化する焦点位置を比較することによって画像の非点収差を探すことが含まれる。この目的のため、簡単なPSFモデルは、フィット係数aiの非常に限られた数で採用されている。
【0160】
【0161】
【0162】
【0163】
【0164】
この形式のデータフィードバックは、システムの位置合わせが実行された主要な方法であって、動作において非常に安定していることが証明されている主要な方法である。
【0165】
数式37のモデルは、1次元LSFに関する情報のみを含むESFデータからPSFの2次元再構成を取得できるように、xおよびyに関して分離可能であるように選択された。そのため、このモデルでは、高度に収差または焦点がぼけたシステムから生じる複雑なPSFを示せない場合があるが、十分に補償された焦点の合ったシステムのPSFを記述するには十分である。ESFは、検出器MTFの寄与によるPSFの劣化を反映するようにさらに変更でき、これは、最悪の場合、各ピクセルの全長にわたって数式38の正確なESFを平均化することによって達成される。
【0166】
例えば線形最小二乗回帰に基づくデータフィッティングは、数式38を実験データに対してフィッティングすることによって実行されて係数aiを得る。PSFモデルは、フィットから得られたaiの値を使用して数式37を評価することによって構築される。次に、EE値は、正方形の領域にわたって数式37を積分することによって分析的に演算されて、数式40を得る。
【0167】
アルゴリズムの精度
アルゴリズム全体の精度を、LSF×LSF近似の精度を評価するために使用されたのと同じ実世界のPSFを使用して、既知のPSFを使用した合成画像を使用して評価した。次に、PSFを、無限に鮮明な試験ターゲットの画像で畳み込んだ。畳み込みは、10×14の傾斜した正方形のグリッドを持つ理想的なノイズのない画像を示す。正方形は、カメラのピクセルグリッドに対する各正方形の配置が、正方形の中心がどこにあるかによってサブピクセルレベルで異なるため、互いに重ね合わせることができないという点で、様々な試験を提供する。データに対するショットノイズの影響を理解するために、最初に画像をスケーリングしてターゲットの実際の画像の強度を一致させ、ベースラインを追加してシステムのバックグラウンド強度と試験ターゲットのブリードスルーの影響を加えた。スケーリングされた画像は、光電子の単位で表した。各ピクセルに対して特別に演算された分布から各ピクセルの強度をサンプリングすることにより、このスケーリングされた画像からショットノイズサンプリング画像を作成した。ショットノイズのポアソン分布は、スケーリングされた画像の平均に対応する平均と、画像の平均にも等しい分散と、を持った正規分布として近似した。画像は、最終的に216-1の画像カウントあたり10,550個の光電子のスケールファクタを使用することによって、16ビットTIFFの上位ビットとして12ビットの深さにスケーリングした。
【0168】
各画像を作成した後、完全な分析パイプラインを実行して、数式35を使用して各正方形のエンスクエアドエネルギーを推定した。さらに、PSFのモデルを、フィット係数a
x;iおよびa
y;iから作成し、直接の比較を可能にした。
図4A~
図4Fは、真のEEが0.594の適度に鮮明なPSFの結果を示している。PSFの3Dプロットは、真のPSF(
図4A)、LSF×LSF近似から導出されたPSF(
図4B)、および最後に合成画像の分析とESFデータのフィッティングとから復元されたPSFモデル(
図4C)を示している。3つのPSFの差異を、同じ垂直スケールを使用して
図4D~
図4Fに示す。予想どおり、PSFモデルは、LSF×LSF近似と非常によく似ている。なお、ベースでの真のPSFの密度は、プロットの左前よりもプロットの右前で速く低下する。そのプロパティは、LSF×LSFモデルによって保持され、ESFデータからのPSFモデルも忠実に回復する。LSF×LSF近似は、ピークの過小評価と、ピークから発生する水平と垂直の行と列の過大評価を引き起こす。これにより、真のEEが過小評価され、避けられない近似となる。PSFを上から見ると、特徴的な「+」字型のPSFモデルとなる。それにもかかわらず、LSF×LSFモデルの忠実な再構築と、フィット係数a
x;iおよびa
y;iのみからのEEの正確な演算とは、ESFフィッティング手順が期待どおりに動作していることを示唆している。
【0169】
2つの追加のPSFを同じ方法でフィットし、それぞれ鮮明なPSFと不鮮明なPSFを示す
図5A~
図5Fと
図6A~
図6Fとに示す。なお、予測どおり、不鮮明なPSFの非対称性はLSF×LSFモデルによっては適切に撮像されない。
【0170】
図7Aで使用したものと同じPSFを使用して、31個所の異なる焦点位置の範囲にわたって演算された9つの異なる実験的なPSFに対して合成画像分析を繰り返した。数式35を使用して推定されたEEと既知の開始PSFのEEとの対比の散布プロットを
図7Bに示す。散布プロットの各点は、10×14の正方形で構成される画像を示し、各正方形は、ショットノイズとサブピクセルの配置に関して異なる複数の正方形が表示されたときに露出するアルゴリズムエラー(存在する場合)との両方から生じるノイズとともに、それら正方形がEEの同じ推定値を生成することを期待して、個別にサンプリングされる。エラーバーは、140回の測定すべての標準偏差を示す。点は、各PSFの最適な焦点の周りに一連の同様のPSFを作成するPSFの生成から生じる様々な場所での自然なクラスタリングを示す。
【0171】
EEは、アルゴリズムによって過小評価されており、これは、コアLSF×LSF近似になると予想される。EEの推定は、LSF×LSFの近似と正確に一致すると予想されるが、この場合、わずかに高くなっている。これにより、全体的な推定がより正確になるが、ESFフィッティングおよびEE推定の数値手法のパフォーマンスを向上させる機会も得られる。
【0172】
図7BのEE推定の精度は、機器試験の目的にとっては非常に許容できるものであり、本明細書に開示された方法によって採用された基本的なアプローチが、使用可能なデータと有用な試験プラットフォームを生成できることを示している。
【0173】
アルゴリズムの正確さ
図7BのEE推定の標準偏差が比較的厳密であることから判断すると、アルゴリズムの正確さは許容できるものである。ただし、光学的位置合わせのためにフィードバックを使用すると、ちょうど厳密な標準偏差を超えた追加レベルの要件が提示される。典型的には、システムに小さな摂動を加えて最適化を実行する場合、フィードバック信号が調整変数の関数として飛び交わないことが重要である。
図7Cは、EEと焦点位置との対比のプロットであり、アルゴリズムが焦点ぼけの摂動に対して非常に滑らかで予測可能な応答曲線を提供することを示している。赤線と黒線(またはそれぞれ円形と三角形のデータ点)の間の距離は、EEが過小評価されていることを示しているが、それは、体系的かつ一貫して過小評価されており、焦点ぼけ位置が点ごとに変化しても、その方法の不正確さは大きく飛び交うことはない。したがって、この曲線は、EEが最大化される焦点ぼけ位置を決定するためなど、より高いレベルの解釈とフィッティングに適している。
【0174】
非対称LSFのフィッティングを可能にする基底系を含め、LSFを近似するために代替の基底系を利用できることは上記したとおりである。そのような方法を利用した場合、
図7Bに示すよりも高い精度が得られる可能性がある。ただし、
図7Cの形式の応答曲線が、最適化の使用に十分滑らかであることを検証する必要がある。証明されていない懸念は、より多くのフィッティングパラメータを使用すると滑らかさが損なわれることである。ただし、この仮定は試験されないままであり、より高精度のモデルは期待に反し、実際には同じように滑らかまたはより滑らかになる可能性がある。
【0175】
EEスコア
エンスクエアドエネルギー(EE)値は、試験対象の光学システム104のイメージングパフォーマンスの測定基準として利用できる。ユーザは、システムの光学部品に対して初期設定または調整を行ってもよい(例えば、部品間、および部品と試験ターゲット116との間の位置合わせ)。次いで、ユーザは、本明細書に記載の方法に従って、光学システム104の試験を実行してもよい。試験の結果に基づいて、ユーザは、その結果を改善するためにさらに調整を行う必要があるかどうかを判定してもよい(例えば、位置合わせおよび焦点合わせをさらに最適化する)。さらに調整を行った後、ユーザは、試験を再度実行してもよい。試験と調整の繰り返しを、結果が許容できると見なされるまで続けてもよい。
【0176】
光学システム104は、画像内のすべてのフィールド点が所定の最小閾値EE値を超えた場合(例えば、EE>0.51)、イメージングパフォーマンス試験に合格したと判定できる。ただし、すべてのEE値が合格したかどうかのみに基づいてイメージングパフォーマンスを評価することは、状況によっては不利と見なされる場合がある。このアプローチは、光学システム104がいくつかのぼんやりしたスポットを示すが、そうでなければ鮮明な(高画質)システムである場合には許されないかもしれない。さらに、位置合わせ手順の進行中に光学パフォーマンスが継続的に評価されるため、閾値処理の使用から生じる不連続性により、閾値処理方法を、連続測定基準が好ましい最適化プロセスとして使用するのに問題が生じ(例えば、レーザ出力に対する位置合わせ、または最小二乗平均平方根(RMS)波面誤差に対する位置合わせ)、閾値処理効果が製造歩留まりの演算に役立たなくなる場合がある。
【0177】
前述の潜在的な不利な点を克服するために、ここで開示されている方法のさらなる実施の形態では、EEスコアであるESを、EE値に基づいて演算する。EEスコアは、イメージングパフォーマンス試験に合格するために所定の最小閾値を超える(または満たす、または超える)必要がある単一の連続変数を使用してフィールド全体を記述する。例えば、EEスコアはES≧100が合格スコアとなるように重み付けしてもよい。EEスコアには次の2つの要素がある。<EE>は、指定された被写界深度にわたってすべてのフィールド点で観測されたすべてのEE値の空間平均であり、重なり積分EOは、異なるイメージングチャネル間(例えば、画像がフィルタリングされる4つの色に対応する4つのチャネル、ここで、図示の実施の形態では、2つのイメージングデバイス124および128が調整可能なフィルタアセンブリ148と組み合わせて利用されて、それぞれ4つのチャネルのうちの2つで画像を取得する)のEEの一貫性を測定するものである。フィールド位置(x,y)と焦点位置(z)との両方で鮮明な領域および不鮮明な領域が一致する場合、EOに報酬が与えられる。EOには、4つのチャネルが矛盾している場合、例えば、2つのイメージングデバイス124および128が互いに対して十分に焦点が合っていない場合、または青色チャネルが左上隅で鮮明であり、右下隅で不鮮明である一方で赤のチャネルが反対側である場合等に、ペナルティを課す。EEスコアとその構成要素は、以下の数式41~44に従って演算される。
【0178】
【0179】
【0180】
【0181】
【0182】
数式中、xおよびyは、サイズNx×Nyのグリッドに配列されたフィールド点を示し、cは、色であり(例えば、4チャネルの実施の形態では、色1、2、3、および4は、青、緑、赤、琥珀色に対応してもよい)、zは、焦点位置であり、nmは、ナノメートルである。
【0183】
EEスコアは、すべてのチャネルからの寄与を考慮するが、一度に単一チャネルを考慮し、フィールド位置および焦点合わせ位置に対する均一性のみを考慮した同様のスコアを持つことも有用である。チャネルごとのEEスコアECcは、単一チャネルの不均一性が、平均二乗の分散と同様な、数式46に捕捉されているような例である。
【0184】
【0185】
【0186】
【0187】
【0188】
数式46の積分がゼロの場合、不均一性EHによって表されるペナルティはゼロになる。
【0189】
EE対Zのフィッティング
数式42、43、および46の積分は、焦点に対してイメージングパフォーマンスを特徴づけることを目的としており、実験データの単純な合計を使用して演算することができる。あるいは、最初にEE値の実験的なペアとそれに対応するz場所をガウスプロファイルにフィッティングすることによって積分を演算することもできる(数式49)。次に、積分は、以下の数式49~55ごとのガウスフィット関数の定積分の分析解として演算される。
【0190】
【0191】
【0192】
【0193】
【0194】
【0195】
【0196】
【0197】
数式中、aは、焦点に対して得られたピークEE値、wは、ガウス幅であり、z0は、ピークEE値に対応するz値である。これは、z0をzステップサイズの数分の1以内に正確に補足しながら、比較的粗いzステップをとる(したがってデータ取得が少なくなる)ことを可能にするという点で特に有用である。数式55は、上記の数式42の積分を演算するために利用してもよい。
【0198】
EO演算の場合、数式44の異なるチャネル間の重なり積分は、数式56~63に示すように、実験的なデータのガウスフィットの分析積分を使用して演算することもできる。
【0199】
【0200】
【0201】
【0202】
【0203】
【0204】
【0205】
【0206】
【0207】
数式中、W1およびW2は、チャネル1および2に対応するガウス分布の幅である。2つのガウス分布の積は、幅WPを持つ、Z0Pを中心とした第三のガウス分布である。数式56~63を利用して、上記の数式47を演算してもよい。
【0208】
ペナルティは、様々な光学的欠陥によってどれだけのパフォーマンスが失われたかを示す。
【0209】
光学システム104に対する現在のEEスコアを推定する以外に、この方法は、1つもしくは複数のアーチファクトまたは光学的欠陥(例えば、収差)に起因してどれだけの潜在的スコアが失われるかをユーザに伝えることを伴ってもよい。このようにして、この方法は、どの特定の欠陥がイメージングパフォーマンスに最も悪影響を及ぼしているかについてユーザに通知してもよい。この目的のために、この方法は、1つまたは複数のそれぞれのアーチファクトに関連付けられた1つまたは複数のペナルティを演算してもよい。ペナルティは、現在のEEスコアと、欠陥が説明された場合に得られるスコアとの差分である。この方法は、ユーザに表示される画像のフィールド全体にわたるペナルティに対応する情報を表示することを含んでもよい。EEスコアへの<EE>の寄与のみを考慮する場合(重なりスコアEOなし)、ペナルティ演算は常により大きな<EE>量になるはずである。ペナルティは、実際のスコアと人為的に増加したスコアとの差分として伝達され、これは、常に負の数である必要がある。重なりスコアEOを考慮する場合、ペナルティ演算で正の数になることがあってもよい。
【0210】
ペナルティの一例は、イメージングデバイス(例えば、カメラ)のZシフトによるペナルティである。このペナルティを演算する際に、この方法は、イメージングデバイス124および128の両方が互いに対して最適に焦点を合わせることができると想定してもよいが、イメージングデバイス124および128内の軸上色収差がゼロであることは想定しない。この方法は、第1のイメージングデバイス124が取得する2つの波長フィルタリングされたイメージングデータの平均(例えば、青および緑のデータの平均)として、第1のイメージングデバイス124にとっての理想的な焦点位置を演算する。この方法はまた、第2のイメージングデバイス128が取得する2つの波長フィルタリングされたイメージングデータの平均(例えば、赤および琥珀色のデータの平均)として、第2のイメージングデバイス128にとっての理想的な焦点位置を演算する。
【0211】
ペナルティの別の例は、イメージングデバイス(例えば、カメラ)の傾斜によるペナルティである。各チャネルについて、この方法は、両方の色(例えば、チャネル1の青と緑)の平均として焦点面を構築し、線形最小二乗法を使用して平面にフィットさせ、各イメージングデバイス124および128の焦点差を維持しながら焦点から除去される先端部分および傾斜部分を見つける。
【0212】
ペナルティの別の例は、非点収差によるペナルティである。水平LSFおよび垂直LSFのEEスコアは、個別に演算され、焦点面の傾斜および高さの各方向への最適化を可能にする。次に、マップを平均化して、最終的な推定を作成する。
【0213】
ペナルティの別の例は、軸上色収差によるペナルティである。このペナルティは、チャネルごとに1回ずつ、2回演算される。チャネル1の場合、2つの色(例えば、青と緑)の最適な焦点がそれぞれについて見つかり、その差分が分割され、各チャネルはその量だけオフセットされる。結果として得られる焦点位置は、カメラ間の焦点の差異が不完全なため、依然として理想的ではない可能性がある。
【0214】
ペナルティの別の例は、フィールドの平坦性によるペナルティである。このペナルティは、チャネルごとに1回ずつ、2回演算される。チャネル1の場合、2つの色の露出(例えば、青と緑の露出)の最適な焦点面が平均化され、最適なフィット平面が(例えば、青と緑の)データの焦点位置から差し引かれる一方、チャネル2は、正常に処理される。このペナルティは、解釈が矛盾する可能性がある。例えば、両方のチャネルが等しく非フラットである場合、両方とも<EE>の観点から問題が発生するが、重なりEOは生産的な方法で強化される。チャネルをフラット化すると、両方の影響がなくなる。通常、<EE>の影響が支配的であるが、これは保証されない。
【0215】
図8Aは、試験ターゲット116のターゲット特性313のうちの1つから取得された画像の例である。
図8Bは、
図8Aの画像から演算されたESFの例である。
図8Bでは、画像座標ゼロは、イメージングされたターゲット特性313のエッジのうちの1つに対応する。
図8Bのデータ点の一部は外れ値である。図示のように、これらの外れ値のうちの1つまたは複数は、ターゲット特性313上の暗いアーチファクト(例えば、製造からの汚れまたは破片)に対応してもよい。この方法は、ESFデータの外れ値を除去する(例えば、無視または削除する)ことを含んでもよい。
【0216】
図9は、
図8BのESFを誤差関数にフィットさせ、ESFからEEを演算した結果の例である。
【0217】
図10は、ターゲット特性313の4つのエッジについてのESFを演算し、ESFを誤差関数にフィットさせ、対応するESFから4つの独立したEE値を演算する例である。名目上、左エッジと右エッジは、同じ測定値になり(EE
left=EE
right)、上エッジと下エッジは、同じ測定値になる(EE
up=EE
down)。ただし、左エッジと右エッジの測定値は、必ずしも上エッジと下エッジの測定値と同じであるとは限らない。
【0218】
代替の実施の形態では、エンスクエアドエネルギーの代わりに、エンサークルドエネルギーをイメージングパフォーマンス測定基準として利用してもよい。さらなる代替の実施の形態では、エンスクエアドエネルギーの代わりに、ストレールレシオをイメージングパフォーマンス測定基準として利用してもよい。
【0219】
図11は、本開示の一つの実施の形態による、本明細書に記載の光学システム104などの光学システムのイメージングパフォーマンスを試験する方法の一例を示す流れ
図1100である。図示のように、この方法は、対物レンズのZステージをおおよその焦点位置に移動することによって開始される。次に、ユーザは、Zスタック取得のプロセスを開始し、これにより、試験ターゲットの複数の画像が異なる焦点位置で取得される。次に、光学システムのハードウェアがZスタックの取得を実行し、その実施の形態を
図12に示す。次に、Zスタック範囲を、Zスタック取得中の焦点パフォーマンスと比較する。次に、Zスタックの結果に基づいてZ焦点位置を更新する。次に、焦点がZスタック範囲の中央の3分の1にあるかどうかを判定する。「はい」である場合、この方法は、現在取得されているデータが使用可能であると判定し、グローバル演算をZスタックに対し実行し、そのように生成された新しいデータがユーザに通知される。
【0220】
一方、Zスタック範囲の中央3分の1に焦点が合わないと判定された場合は、現在取得されているデータからイメージングパフォーマンス値(例えば、本明細書で説明するエンスクエアドエネルギーの推定値)の演算が行われるが、ユーザに対し、範囲の中央に焦点が見つからなかったためにデータが不完全または信頼できないという警告がなされる。次に、このZスタックから得られた焦点位置のより多くの情報に基づいた推測を用いて、Zスタック取得のプロセスを繰り返してもよい。
【0221】
図12は、本開示の一つの実施の形態による画像のZスタックの取得を示す流れ
図1200である。ユーザがZスタック取得を開始した後(
図11)、初期の焦点推測の周囲でZ位置の範囲を演算する。FOVの組内の各視野FOV
iについて、Zステージは、FOV
iに対応する焦点ぼけ位置に移動するように命令される。次に、Zステージの実際の位置を照会し、記録する。次に、すべてのイメージングデバイスからの画像を同時に取得する。次に、各画像自体の高速処理を実行し、その実施の形態を
図13にさらに示す。次に、画像データをプールする。次に、ESFごとに1つずつの複数のPSFモデルを、本明細書で説明するように、ガウス誤差関数の基底系にESFをフィットすることによって構築する。
【0222】
図13は、本開示の一つの実施の形態による、画像自体の高速処理を示す流れ
図1300である。その画像を、強度、水平重心、および垂直重心の表現にダウンサンプリングする。次に、ダウンサンプリングされた画像を使用して、粗いグローバルレジストレーションを実行する。各正方形(または試験ターゲットの他のタイプのターゲット特性)は、重心データを使用して位置付けされる。次に、複数の正方形の位置を、平行移動、回転、および拡大のグローバルモデルにフィットする。次に、各正方形の各エッジからESFを抽出し、エッジの角度と中心を演算する。次に、ESFデータから外れ値を削除する。次に、複数のPSFモデルを、本明細書で説明するように、ガウス誤差関数の基底系に独立して各ESFをフィットすることによって構築する。次に、各ESFについて、イメージングパフォーマンス値(例えば、本明細書で説明するエンスクエアドエネルギーの推定値)を演算する。データは、平均化できる(EE
horizontal=(EE
left=EE
right)/2、EE
vertical=(EE
up=EE
down)/2)。
図13に示されるプロセスは、複数の画像に対して繰り返されてもよく、結果として得られた画像データは、プールされ(
図12)、さらに処理できる。
【0223】
一つの実施の形態では、流れ
図1100、1200、および1300のうちの1つまたは複数は、上記で説明され、
図11、
図12、および
図13のうちの1つまたは複数に示されているステップまたは機能を実施(例えば、制御または実行)するように構成された光学イメージングパフォーマンス試験システムまたは光学イメージングパフォーマンス試験システムの一部を示してもよい。この目的のために、プロセッサ、メモリ、および当業者に理解される他の構成要素を含むコントローラ(例えば、
図1に示されるコントローラ108)を設けて、ステップまたは機能の実施に関与する試験システムのコンポーネントを制御することなどによって、ステップまたは機能のパフォーマンスを制御してもよい。これは、当業者によって理解されるように、ユーザ入力および/または機械実行可能命令を受信および処理することに関与してもよい。
【0224】
図14は、方法に従って生成されてもよく、また、コンピューティングデバイス(例えば、コンピューティングデバイス168)のディスプレイスクリーン上でユーザに表示できるEEパフォーマンスマップのセットの例を示している。マップの列は異なる色に対応しており、本実施例では、左から右に、青、緑、琥珀色、および赤である。最も上の行(行1)では、各マップに視野(FOV)の点ごとに平均EEが表示されている。次の行(行2)では、各マップに、FOVの点ごとに焦点ぼけ値で見つかった最良のEEが表示されている。次の行(行3)では、各マップに、EEによって判定されたものとして、最適な焦点の表面が表示されている。次の行(行4)では、各マップに、EEによって判定されたものとして、水平方向のイメージングデータと垂直方向のイメージングデータとの対比の非点収差が表示されている。最も下の行(行5)では、各マップに、EEが最小値を超えるz範囲として被写界深度が表示されている。
【0225】
図15は、方法に従って生成されてもよく、また、コンピューティングデバイス(例えば、コンピューティングデバイス168)のディスプレイスクリーン上でユーザに表示できる組み合わせEEのスコアおよびパフォーマンスマップの例を示している。最も上のマップ(1)は、すべてのチャネルおよび焦点ぼけ位置で使用可能なEEを示すグローバルEEマップである。最も上のマップ(1)の下において、マップのグループ(2)は、上記のような様々な光学的欠陥の寄与を示している。したがって、これらのマップ(2)は、ユーザが問題の欠陥に対処した場合に、イメージングパフォーマンスがどれだけ向上するかをユーザに通知する。本実施例では、マップ(2)のグループの最初の行は、異なる色、本実施例では、左から右に、青、緑、琥珀、赤のグローバルスコアを示している。マップ(2)のグループの2行目は、左から右に、重なり、イメージングデバイスZシフトによるEEへのペナルティ、チャネル1の軸上色収差によるEEへのペナルティ、およびチャネル2の軸上色収差によるEEへのペナルティを示している。マップ(2)のグループの3行目は、左から右に、非点収差によるEEへのペナルティ、イメージングデバイスの傾斜によるEEへのペナルティ、チャネル1のフィールドの平坦性によるEEへのペナルティ、およびチャネル2のフィールドの平坦性によるEEへのペナルティを示している。
【0226】
図16は、本開示の一つの実施の形態による、光学イメージングパフォーマンス試験システムの一部であるか、または光学イメージングパフォーマンス試験システムと通信できるシステムコントローラ(またはコントローラ)1600の非限定的な例の模式図である。例えば、システムコントローラ1600は、上記で説明され、
図1に示されている試験システム100のシステムコントローラ108に対応してもよい。
【0227】
図示の実施の形態では、システムコントローラ1600は、全体的な制御を提供する主要な電子プロセッサを代表するプロセッサ1602(典型的は電子機器ベース)と、専用制御動作または特定の信号処理タスクのために構成された1つまたは複数の電子プロセッサ(例えば、グラフィックス処理ユニットまたはGPU、デジタルシグナルプロセッサまたはDSP、特定用途向け集積回路またはASIC、フィールドプログラマブルゲートアレイまたはFPGA等)と、を含む。システムコントローラ1600はまた、データおよび/またはソフトウェアを格納するための1つまたは複数のメモリ1604(揮発性および/または不揮発性)を含む。システムコントローラ1600はまた、1つまたは複数のタイプのユーザインターフェースデバイスを制御し、ユーザインターフェースデバイスと、ユーザインターフェースデバイスと通信するシステムコントローラ1600のコンポーネントと、の間のインターフェースを提供するための1つまたは複数のデバイスドライバ1606を含んでもよい。そのようなユーザインターフェースデバイスは、ユーザ入力デバイス1608(例えば、キーボード、キーパッド、タッチスクリーン、マウス、ジョイスティック、トラックボールなど)およびユーザ出力デバイス1610(例えば、ディスプレイスクリーン、プリンタ、視覚的インジケータまたはアラート、可聴インジケータまたはアラートなど)を含んでもよい。
【0228】
様々な実施の形態では、システムコントローラ1600は、1つまたは複数のユーザ入力デバイス1608および/またはユーザ出力デバイス1610を含むものとして、または少なくともそれらと通信するものとして見なされてもよい。システムコントローラ1600はまた、メモリおよび/または1つまたは複数のタイプのコンピュータ可読媒体1614に含まれる1つまたは複数のタイプのコンピュータプログラムまたはソフトウェア1612を含んでもよい。コンピュータプログラムまたはソフトウェアは、試験システム100の様々な動作を制御または実行するための非一時的な命令(例えば、論理命令)を含んでもよい。コンピュータプログラムまたはソフトウェアは、アプリケーションソフトウェアおよびシステムソフトウェアを含んでもよい。システムソフトウェアは、ハードウェアとアプリケーションソフトウェアとの間の相互作用を含む、システムコントローラ1600の様々な機能を制御および管理するためのオペレーティングシステム(例えば、Microsoft Windows(登録商標)オペレーティングシステム)を含んでもよい。特に、オペレーティングシステムは、ユーザ出力デバイス1610を介して表示可能であり、かつユーザがユーザ入力デバイス1608を使用して対話できるグラフィカルユーザインターフェース(GUI)を提供してもよい。
【0229】
システムコントローラ1600はまた、GUIによるグラフ形式での提示のためのデータのフォーマットを含む、試験対象の光学システム(例えば、上記で説明され、
図1に示されている光学システム104によって出力された信号(例えば、イメージングデータ)を受信および処理するための1つまたは複数のデータ取得/信号調整コンポーネント(DAQ)1616(ハードウェア、ファームウェアおよび/またはソフトウェアで具体化できるものとして)を含んでもよい。DAQ1616はまた、制御信号を光学システムに送信して、本明細書に記載されているような調整可能な光学部品の動き/位置決めを制御するように構成できる。DAQ1616は、上記で説明され、
図1に示されている電子機器モジュール164の全部または一部に対応してもよい。
【0230】
システムコントローラ1700は、本開示の全体にわたって説明するように、イメージングパフォーマンス測定基準、スコア、マップなどを含む、光学システムから出力された信号を処理し、そこからデータを生成するように構成されたデータアナライザ(またはモジュール)1618をさらに含んでもよい。したがって、データアナライザ1618は、本明細書に開示される方法のいずれかの全部または一部を実装(制御または実行)するように構成できる。これらの目的のために、データアナライザ1618は、当業者によって理解されるように、ソフトウェアおよび/または電子機器(ハードウェアおよび/またはファームウェア)で具体化できる。データアナライザ1618は、上記で説明され、
図1に示されている電子機器モジュール164および/またはコンピューティングデバイス168の全部または一部に対応してもよい。
【0231】
図16は、本開示と一致するシステムコントローラ1600の例の高レベルの模式図であることが理解されよう。追加の構造、デバイス、電子機器、およびコンピュータ関連または電子プロセッサ関連の構成要素などの他の構成要素が、実際の実施態様の必要に応じて含まれてもよい。システムコントローラ1600は、備えることができる構造(例えば、回路、メカニズム、ハードウェア、ファームウェア、ソフトウェア等)を示すことを意図した機能ブロックとして
図16に模式的に表されていることも理解されよう。様々な機能ブロックおよびそれらの間の信号リンクは、図示のみを目的として任意に位置付けられており、いかなる方法においても制限するものではない。当業者であれば、実際には、システムコントローラ1600の機能が、様々な方法で実施できるが、必ずしも
図16に示され、本明細書で例によって説明されているままの方法であるとは限らないことを理解するであろう。
【0232】
例示的な実施の形態
ここで開示されている主題に従って提供される例示的な実施の形態には、以下が含まれるが、これらに限定されない。
【0233】
1.光学システムのイメージングパフォーマンスを試験する方法であって、試験ターゲットを光学システムの物体平面に位置付けることと、光学システムを動作させて試験ターゲットを照明し、画像ビームを生成することと、光学システムの焦点合わせステージを動作させて、複数の焦点ぼけの値に対応する画像ビームから試験ターゲットの複数の画像を取得することと、各画像から、試験ターゲット内の複数の場所での複数のエッジ広がり関数を演算することと、それぞれのエッジ広がり関数から複数の点広がり関数を構築することと、点広がり関数に基づいて、複数の場所に対応する複数のイメージングパフォーマンス値を演算することであって、イメージングパフォーマンス値が、エンスクエアドエネルギー、エンサークルドエネルギー、およびストレールレシオからなる群から選択される計測基準に基づいている、演算することと、を含む、方法。
【0234】
2.複数の場所は、試験ターゲット内の物体平面内の複数のフィールド座標(x,y)を含む、実施の形態1に記載の方法。
【0235】
3.複数の場所は、試験ターゲットを通過する光軸に沿った複数の焦点位置(z)を含み、光学システムを動作させることは、異なる焦点位置(z)での試験ターゲットの複数の画像を取得すること、を含む、実施の形態1または2に記載の方法。
【0236】
4.イメージングパフォーマンス値の組み合わせに基づいてイメージングパフォーマンスの1つまたは複数のマップを生成することを含む、実施の形態1~3のいずれか1つに記載の方法。
【0237】
5.マップのうちの2つ以上を比較して、物体平面に対する各イメージングデバイスまたはチャネルの焦点面の相対的な位置合わせの尺度を提供すること、を含む、実施の形態4に記載の方法。
【0238】
6.1つまたは複数のマップは、異なるイメージングチャネルに対応し、異なるイメージングチャネルは、画像を取得するように動作する光学システムの異なるイメージングデバイス、または取得される画像の異なる波長、または異なるイメージングデバイスと異なる色との両方に対応する、実施の形態4または5に記載の方法。
【0239】
7.マップのうちの2つ以上を比較して、互いに対する各イメージングデバイスまたはチャネルの相対的な位置合わせの尺度を提供することと、前述の両方と、を含む、実施の形態6に記載の方法。
【0240】
8.1つまたは複数のマップを生成した後、光学システムの1つまたは複数の光学部品の位置を調整すること、または1つまたは複数のマップによって提供される情報に基づいて1つまたは複数の光学部品を交換すること、を含む、実施の形態4~7のうちのいずれか1つに記載の方法。
【0241】
9.1つまたは複数の光学部品は、イメージングデバイスのうちの1つまたは複数と、光学システムの対物レンズと、光学システムの1つまたは複数のチューブレンズと、1つまたは複数のミラーまたはダイクロイックミラーと、前述の2つ以上の組み合わせと、からなる群から選択される、実施の形態8に記載の方法。
【0242】
10.1つまたは複数のマップは、1つまたは複数の初期マップであり、1つまたは複数の光学部品を調整または交換した後、試験ターゲットの新しい画像を取得し、複数の新しいイメージングパフォーマンス値を演算し、イメージングパフォーマンスの1つまたは複数の新しいマップを生成すること、をさらに含む、実施の形態8または9に記載の方法。
【0243】
11.1つまたは複数の新しいマップを1つまたは複数の初期マップと比較して、イメージングパフォーマンスを最適化するために1つまたは複数の光学部品に対して行われる位置調整を決定すること、を含む、実施の形態10に記載の方法。
【0244】
12.調整することまたは交換することで、光学システムにおいて共役像面の最適な一対が見つかる、実施の形態8~11のいずれか1つに記載の方法。
【0245】
13.調整することまたは交換することで、イメージングデバイスの焦点マッチングと、イメージングデバイスの傾斜と、像面湾曲の平坦化と、非点収差の低減と、波長に依存する焦点シフトの低減と、前述の2つ以上の組み合わせと、からなる群から選択される属性が改善する、実施の形態8~12のいずれか1つに記載の方法。
【0246】
14.イメージングパフォーマンス値の組み合わせに基づいてイメージングパフォーマンスの1つまたは複数のグローバルスコアを演算すること、を含む、実施の形態の1~13のいずれか1つに記載の方法。
【0247】
15.1つまたは複数のグローバルスコアを変更して、試験ターゲット内の物体平面内のフィールド座標(x,y)の範囲にわたって、または物体平面の焦点位置(z)の範囲を通じて、または前述の両方において、イメージングパフォーマンス値の不均一性に対し、ペナルティを課すか、または報酬を与えること、を含む、実施の形態14に記載の方法。
【0248】
16.1つまたは複数のグローバルスコアを変更して、試験ターゲット内の物体平面内のフィールド座標(x,y)の関数、または物体平面の焦点位置(z)の関数、または前述の両方として、異なるイメージングチャネルの類似性に対し、ペナルティを課すか、または報酬を与えること、を含み、異なるイメージングチャネルが、画像を取得するように動作する光学システムの異なるイメージングデバイスに、または取得される画像の異なる波長に、または異なるイメージングデバイスと異なる色との両方に、対応する、実施の形態14または15に記載の方法。
【0249】
17.イメージングパフォーマンス値を演算した後、イメージングパフォーマンス値によって提供される情報に基づいて光学システムの1つまたは複数の光学部品の位置を調整すること、を含む、実施の形態の1~16のいずれか1つに記載の方法。
【0250】
18.イメージングパフォーマンス値は、初期イメージングパフォーマンス値であり、1つまたは複数の光学部品を調整した後、試験ターゲットの新しい画像を取得し、複数の新しいイメージングパフォーマンス値を演算すること、をさらに含む、実施の形態17に記載の方法。
【0251】
19.新しいイメージングパフォーマンス値を初期イメージングパフォーマンス値と比較して、イメージングパフォーマンスを最適化するために1つまたは複数の光学部品に対して行われる位置調整を決定すること、を含む、実施の形態18に記載の方法。
【0252】
20.試験ターゲットを位置決めすることは、光学システムの1つまたは複数の光学部品と共有されるデータムに対してターゲットを位置合わせすること、を含む、実施の形態1~20のいずれか1つに記載の方法。
【0253】
21.光学システムを動作させることは、画像ビーム内の対物レンズを利用すること、を含み、画像ビームの軸に沿って対物レンズの位置を調整して、異なる焦点位置(z)での試験ターゲットの複数の画像を取得すること、をさらに含む、実施の形態1~20のいずれか1つに記載の方法。
【0254】
22.対物レンズは、対物レンズが無限共役顕微鏡用に構成されていることと、対物レンズが有限共役顕微鏡用に構成されていることと、からなる群から選択される構成を有する、実施の形態21に記載の方法。
【0255】
23.光学システムを動作させることは、試験ターゲットのそれぞれの画像を取得するために2つ以上のイメージングデバイスを動作させること、を含む、実施の形態1~22のいずれか1つに記載の方法。
【0256】
24.2つ以上のイメージングデバイスは、2つ以上の異なる波長でそれぞれの画像を取得する、実施の形態23に記載の方法。
【0257】
25.試験ターゲットから伝搬する画像ビームを2つ以上の画像ビーム部分に分割することと、2つ以上の画像ビーム部分を2つ以上のイメージングデバイスにそれぞれ送信することと、を含む、実施の形態24に記載の方法。
【0258】
26.光学システムを動作させることは、選択された波長で画像ビームをフィルタリングするためにフィルタアセンブリを動作させること、を含む、実施の形態1~25のいずれか1つに記載の方法。
【0259】
27.光学システムを動作させることは、画像ビームにおいてチューブレンズを利用すること、を含み、複数の画像を取得するためにチューブレンズ内の1つもしくは複数のレンズまたはレンズ群の相対位置を調整して、チューブレンズの異なる位置での試験ターゲットの複数の画像を取得すること、さらに含む、実施の形態1~26のいずれか1つに記載の方法。
【0260】
28.試験ターゲットは、暗い材料と、暗い材料上に配置された明るい特性のアレイと、を含む、実施の形態1~27のいずれか1つに記載の方法。
【0261】
29.明るい特性は、多角形である、実施の形態28に記載の方法。
【0262】
30.明るい特性は、明るい特性のエッジが、画像を取得する光学イメージングシステムのピクセルアレイに対してある角度で配向されるように、傾斜している、実施の形態29に記載の方法。
【0263】
31.試験ターゲットを保持するように構成されたターゲットホルダと、試験ターゲットを照明するように構成された光源と、試験ターゲットの画像を取得するように構成されたイメージングデバイスと、試験ターゲットとイメージングデバイスとの間のイメージング光路に位置付けられた対物レンズであって、対物レンズまたはターゲットホルダのうち少なくとも1つの位置が、イメージング光路に沿って調整可能である、対物レンズと、電子プロセッサとメモリとを備えたコントローラであって、先行する実施の形態のいずれか1つに記載の方法の、複数のエッジ広がり関数を演算するステップと、複数の点広がり関数モデルを構築するステップと、複数のイメージングパフォーマンス値を演算するステップと、を制御するように構成された、コントローラと、を備える、光学イメージングパフォーマンス試験システム。
【0264】
32.対物レンズは、対物レンズが無限共役顕微鏡用に構成されていることと、対物レンズが有限共役顕微鏡用に構成されていることと、からなる群から選択される構成を有する、実施の形態31に記載のシステム。
【0265】
33.イメージングデバイスは、複数のイメージングデバイスを備え、システムは、イメージング光路をそれぞれイメージングデバイスに向けられた複数のイメージング光路に分割するように構成された画像分離ミラーをさらに備える、実施の形態31または32に記載のシステム。
【0266】
34.イメージングデバイスに伝播するためにイメージング光路内の画像ビームの波長を選択するように構成されたフィルタアセンブリを備える、実施の形態31~33のいずれか1つに記載のシステム。
【0267】
35.イメージング光路に位置付けられたチューブレンズを備え、チューブレンズ内の1つもしくは複数のレンズまたはレンズ群の相対位置が調整可能である、実施の形態31~34のいずれか1つに記載のシステム。
【0268】
36.試験ターゲットを備え、試験ターゲットは、暗い材料と、暗い材料上に配置された明るい特性のアレイと、を備える、実施の形態31~35のいずれか1つに記載のシステム。
【0269】
37.明るい特性は、多角形である、実施の形態36に記載のシステム。
【0270】
38.明るい特性は、明るい特性のエッジが、画像を取得する光学イメージングシステムのピクセルアレイに対してある角度で配向されるように、傾斜している、実施の形態37に記載のシステム。
【0271】
39.プロセッサ上で実行されると、先行する実施の形態のいずれか1つに記載の方法の、複数のエッジ広がり関数を演算するステップと、複数の点広がり関数モデルを構築するステップと、複数のイメージングパフォーマンス値を演算するステップと、を実行する、格納された命令を含む、非一時的なコンピュータ可読媒体。
【0272】
40.実施の形態39に記載のコンピュータ可読記憶媒体を備える、光学システムのイメージングパフォーマンスを試験するためのシステム。
【0273】
本明細書に記載のプロセス、サブプロセス、およびプロセスステップのうちの1つまたは複数は、ハードウェア、ファームウェア、ソフトウェア、または前述の2つ以上の組み合わせによって、1つもしくは複数の電子機器またはデジタル制御デバイス上で実行できることが理解される。ソフトウェアは、例えば、
図1または16に模式的に示されているコントローラ108または1600などの適切な電子処理コンポーネントまたはシステムにおけるソフトウェアメモリ(図示せず)に常駐してもよい。ソフトウェアメモリは、論理機能(すなわち、デジタル回路もしくはソースコードなどのデジタル形式で、またはアナログの電気、音、もしくはビデオ信号などのアナログソースなどのアナログ形式で、実装できる「論理」)を実装するための実行可能命令の順序付けられたリストを含んでもよい。
【0274】
命令は、例えば、1つもしくは複数のマイクロプロセッサ、汎用プロセッサ、プロセッサの組み合わせ、デジタル信号プロセッサ(DSP)、フィールドプログラマブルゲートアレイ(FPGA)、または特定用途向け集積回路(ASIC)を含む処理モジュール内で実行できる。さらに、模式図とは、アーキテクチャまたは機能の物理的レイアウトによって制限されない物理的(ハードウェアおよび/またはソフトウェア)実施態様を有する機能の論理的分割を説明するものである。本明細書で説明されるシステムの例は、様々な構成で実装されてもよく、単一のハードウェア/ソフトウェアユニットにおいて、または別々のハードウェア/ソフトウェアユニットにおいて、ハードウェア/ソフトウェアコンポーネントとして動作してもよい。
【0275】
実行可能命令は、電子システムの処理モジュール(例えば、
図1または16のコントローラ108または1600)によって実行されるときに、電子システムに対して命令の実行を指示する、格納された命令を有するコンピュータプログラム製品として実装できる。コンピュータプログラム製品は、例えば、電子コンピュータベースシステム、プロセッサ搭載システムなどの命令実行システム、装置、またはデバイス、もしくは命令実行システム、装置、またはデバイスから命令を選択的にフェッチして命令を実行できる他のシステムによって、またはそれと関連して使用するための任意の非一時的なコンピュータ可読記憶媒体において選択的に具体化できる。
【0276】
本明細書において、コンピュータ可読記憶媒体は、命令実行システム、装置、またはデバイスによって、またはそれに関連して使用するためのプログラムを格納できる任意の非一時的な手段である。非一時的なコンピュータ可読記憶媒体は、選択的に、例えば、電子、磁気、光学、電磁、赤外線、もしくは半導体のシステム、装置、またはデバイスであってもよい。非一時的なコンピュータ可読媒体のより具体的な例の非網羅的なリストには、1つまたは複数のワイヤ(電子)を有する電気接続、ポータブルコンピュータディスケット(磁気)、ランダムアクセスメモリ(電子)、読み取り専用メモリ(電子)、例えば、フラッシュメモリ(電子)などの消去可能なプログラム可能な読み取り専用メモリ、例えば、CD-ROM、CD-R、CD-RW(光)などのコンパクトディスクメモリ、およびデジタル多用途ディスクメモリ、つまりDVD(光)が含まれる。なお、プログラムは、例えば、紙もしくは他の媒体の光学的走査を介して電子的にキャプチャされ、次にコンパイルされ、読み取られ、必要に応じて適切な方法で処理されて、コンピュータメモリもしくはマシンメモリに保存できるので、非一時的なコンピュータ可読記憶媒体は、プログラムが印刷される紙または別の適切な媒体であってもよい。
【0277】
本明細書で使用される「信号通信において」という用語は、2つ以上のシステム、デバイス、コンポーネント、モジュール、またはサブモジュールが、何らかのタイプの信号路を通って伝搬する信号を介して互いに通信できることを意味することも理解されよう。信号は、通信、電力、データ、またはエネルギー信号であり得、これらの信号は、情報、電力、またはエネルギーを第1のシステム、デバイス、コンポーネント、モジュール、またはサブモジュールから第2のシステム、デバイス、コンポーネント、モジュール、またはサブモジュールを、または第1のシステム、デバイス、コンポーネント、モジュール、またはサブモジュールと第2のシステム、デバイス、コンポーネント、モジュール、またはサブモジュールとの間の信号路に沿って、通信してもよい。信号路には、物理的、電気的、磁気的、電磁気的、電気化学的、光学的、有線、または無線接続が含まれてもよい。信号路にはまた、第1のシステム、デバイス、コンポーネント、モジュール、またはサブモジュールと第2のシステム、デバイス、コンポーネント、モジュール、またはサブモジュールとの間の追加のシステム、デバイス、コンポーネント、モジュール、またはサブモジュールが含まれてもよい。
【0278】
より一般的には、「通信する」および「~通信中」(例えば、第1のコンポーネントが第2のコンポーネントと「通信する」または「通信中である」)などの用語は、構造的、機能的、機械的、電気的、信号的、光学的、磁気的、電磁気的、イオン的、または2つ以上のコンポーネントもしくは要素間の流体的関係を示すために本明細書で使用される。したがって、あるコンポーネントが第2のコンポーネントと通信すると言うことは、追加のコンポーネントが第1のコンポーネントと第2のコンポーネントとの間に存在する、および/または互いに動作可能に関連付けられるか、もしくは従事する可能性を排除することを意図するものではない。
【0279】
本発明の範囲を逸脱することなく、本発明の様々な態様または詳細を変更することができることを理解されたい。さらに、上述の説明は、例示のみを目的としており、限定を目的とするものではなく、本発明は、特許請求の範囲によって定義される。
関連出願
本出願は、2019年8月7日に出願され、「光学イメージングパフォーマンス試験システムおよび方法」と題された仮特許出願第62/884,116号の米国特許法第119条(e)に基づく利益を主張し、その内容は、参照によりその全体が本明細書に組み込まれる。
【手続補正書】
【提出日】2022-03-28
【手続補正1】
【補正対象書類名】明細書
【補正対象項目名】0152
【補正方法】変更
【補正の内容】
【0152】
フィット係数X[0]は、バックグラウンドノイズと試験ターゲットのブリードスルーによる画像のベースライン強度を示す。フィットであるデータがESFx(x)となる垂直エッジである場合、係数X[1];X[2];:::;X[Ng]は、数式20において係数ax;1,ax;2,:::,ax;Ngを示す。 フィットであるデータがE
SFy(y)となる水平エッジである場合、係数X[1];X[2];:::;X[Ng]は、数式20において係数ay;1,ay;2,:::,ay;Ngを示す。
【手続補正2】
【補正対象書類名】明細書
【補正対象項目名】0163
【補正方法】削除
【補正の内容】
【手続補正3】
【補正対象書類名】明細書
【補正対象項目名】0166
【補正方法】変更
【補正の内容】
【0166】
例えば線形最小二乗回帰に基づくデータフィッティングは、数式38を実験データに対してフィッティングすることによって実行されて係数aiを得る。PSFモデルは、フィットから得られたaiの値を使用して数式37を評価することによって構築される。次に、EE値は、正方形の領域にわたって数式37を積分することによって分析的に演算されて、数式39を得る。
【手続補正4】
【補正対象書類名】明細書
【補正対象項目名】0230
【補正方法】変更
【補正の内容】
【0230】
システムコントローラ
1600は、本開示の全体にわたって説明するように、イメージングパフォーマンス測定基準、スコア、マップなどを含む、光学システムから出力された信号を処理し、そこからデータを生成するように構成されたデータアナライザ(またはモジュール)1618をさらに含んでもよい。したがって、データアナライザ1618は、本明細書に開示される方法のいずれかの全部または一部を実装(制御または実行)するように構成できる。これらの目的のために、データアナライザ1618は、当業者によって理解されるように、ソフトウェアおよび/または電子機器(ハードウェアおよび/またはファームウェア)で具体化できる。データアナライザ1618は、上記で説明され、
図1に示されている電子機器モジュール164および/またはコンピューティングデバイス168の全部または一部に対応してもよい。
【国際調査報告】