(57)【特許請求の範囲】
【請求項１】
光源によって触発されるトリガー、光学イメージングシステム、DMDマイクロミラー配列、光学集束収集システム、光減衰器、単一光子計数器、ドライブ制御モジュール及び最適化アルゴリズムモジュールを備え、圧縮センシング原理に基づいて前記トリガーを触発する極微弱光源を採用して物体を照射して得られたイメージングスペクトルを秒レベル分解精度にし、更に時間に従って動的に変化する物体をイメージングし、時系列に従って配列された連続グレースケールビデオ画像フレームを出力する時間分解単一光子計数イメージングスペクトルシステムであって、
前記トリガーはその前端に位置する極微弱光触発源で触発されて、該トリガーの出力端がドライブ制御モジュールの入力端と接続されて、トリガーが触発された際ドライブ制御モジュールはドライブ制御信号を出力してその出力端と接続されたDMDマイクロミラー配列と単一光子計数器を作動し始めるように触発し、該DMDマイクロミラー配列が反転し始め、単一光子計数器が同時に計数し始め、光減衰器は光強度を減衰することに用いられ、光学集束収集システムは光線を集束収集することに用いられ、単一光子計数器は光を単一光子計数するためのものであり、
前記単一光子計数器の出力端が前記最適化アルゴリズムモジュールの１つの入力端と接続され、単一光子計数器の計数値を最適化アルゴリズムモジュールの１つの入力パラメータとし、該ドライブ制御モジュールが記憶する選定領域のランダム測定マトリックスを受信してそのもう１つの入力パラメータとするために、前記最適化アルゴリズムモジュールのもう１つの入力端が前記ドライブ制御モジュールの一出力端と接続されて、前記最適化アルゴリズムモジュールは入力した単一光子計数器の測定値により前記ランダム測定マトリックスとともにスパース信号を構成し、且つ前記スパース信号を変換して反転し、二次元光子密度により構成した画像を取得し、即ち画像の各画素の値はDMDマイクロミラー配列の対応する画素に入射した光子密度に等しく、その内、光子密度とは単位時間の単位面積に前記DMDマイクロミラー配列の反射波長範囲内で入射される光子数であり、Ｍ個のｔ時間間隔を経て、時系列に従って配列された時間分解可能な一連の二次元グレースケール画像ビデオフレームを構成することができ、
前記DMDマイクロミラー配列はデジタルマイクロミラーデバイスであり、
前記DMDマイクロミラー配列の反転をドライブ制御するために、前記ドライブ制御モジュールのもう１つの出力端が前記DMDマイクロミラー配列の入力端と接続されて、該ドライブ制御モジュールはDLPデジタル光処理技術に基づいたものであり、該ドライブ制御モジュールがDMDマイクロミラー配列の作動領域を選定した後、ランダム測定マトリックスをダウンロードして前記DMDマイクロミラー配列の反転をドライブ制御し、前記DMDマイクロミラー配列が反転すると同時に前記単一光子計数器へ同期信号を送信して、DMDマイクロミラー配列と単一光子計数器の間の同期を保証し、即ちDMDマイクロミラー配列が一回反転すると、単一光子計数器が該回の反転の時間間隔内の単一光子数を累積計数し、DMDマイクロミラー配列が反転を完成した後、単一光子計数器がクリアして累積計数を再び開始し、すべての計数と該選定領域のランダム測定マトリックスがいずれも最適化アルゴリズムモジュールに伝送される時間分解単一光子計数イメージングスペクトルシステム。
【請求項２】
カラーイメージングに用いられる際、単一光子計数器に代えて単一光子検出器ラインアレイを用い、且つ前記単一光子検出器ラインアレイと光学集束収集システムの間の光路上に分光光度計を設置することを特徴とする請求項１に記載の時間分解単一光子計数イメージングスペクトルシステム。
【請求項３】
前記光学イメージングシステムと光学集束収集システムはいずれも光学レンズ群を採用し、それぞれ光学イメージングと光学集束を担当し、極微弱光は光学イメージングシステムを通過した後、DMDマイクロミラー配列上に同じ大きさ又は縮小又は拡大の画像を形成することができ、実際の要求に従ってイメージング調整することを特徴とする請求項１又は２に記載の時間分解単一光子計数イメージングスペクトルシステム。
【請求項４】
前記光学集束収集システムはファイバー結合技術を用いて、即ち前記分光光度計で分光された光束をファイバーに結合し、ファイバー結合技術を用いて分光がそれぞれ対応の単一光子計数器上に収集され、
前記分光光度計は、光視準部、分光部、角度測定部、光度観察及び測定部を備え、スペクトルの分析と測定に用いられ、
前記分光光度計は、プリズム分光光度計又は格子分光光度計から成ることを特徴とする請求項２に記載の時間分解単一光子計数イメージングスペクトルシステム。
【請求項５】
光源、イメージングスペクトル測定ユニット、電気検出ユニット、システム制御ユニット、処理ユニットを備え、圧縮センシング原理に基づいて、レーザーパルスにより物体蛍光を励起し、光強度分布とイメージングスペクトルを取得し、且つ時間分解方策を用いてピコ秒レベルまでの時間分解精度を実現し、更に蛍光寿命の過度的な周期性物質を検出することに用いられる時間分解単一光子計数イメージングスペクトルシステムであって、
前記光源は、システム制御ユニットで送信した触発パルスの触発で、レーザーパルスを発生し、測定しようとする物体に当てて、該物体は該物体に関する情報が載せられた蛍光を発生し、
前記イメージングスペクトル測定ユニットは、光学イメージングシステム、DMDマイクロミラー配列、光集束収集部材を更に備え、測定しようとする物体情報が載せられた蛍光を前記DMDマイクロミラー配列にイメージングし、該DMDマイクロミラー配列はランダム光変調を採用して且つその上のイメージングを反射し、且つ反射光は光路上に配布された光集束収集部材により視準した後電気検出ユニットに入射することに用いられ、
前記電気検出ユニットは、時間分解方策に従って入射光を検出し且つ入射光の光子数を記録し、且つ各光子が到達する時間次元情報を記録することに用いられ、
前記処理ユニットは、入力したすべての各時間帯の光子数、光子の到達時間の次元情報及びランダム光変調マトリックス情報により、圧縮センシングとイメージングスペクトルアルゴリズムを用いてマルチパラメーター情報のイメージングスペクトル構成を実現し、且つ8種類のイメージングを出力することに用いられ、
前記システム制御ユニットは、各部材のイネーブルを備え、即ちドライブイネーブル信号を上記各部材に送信してそれを正常な動作を開始するようにし、且つ光源触発パルスと光子計数に対する同期を実現することに用いられ、且つマイクロミラー配列の反転とその上のランダムマトリックスの交換を制御し且つ過度的な周期の開始前又は終了後に反転をすることに用いられ、
前記DMDマイクロミラー配列はランダム数発生器によりベルヌーイバイナリランダムマトリックスをロードして、光に対するランダム光変調を実現し、前記時間分解方策は非周期変化の長い時系列過程に対してフレームごとに測定の方法を適用し、周期変化特性を有する過度的な過程に対して、最適化の時間分解方策を採用し、該最適化の時間分解方策は、具体的に、過度的な過程の周期を予め定められた数の検出サブ時間帯に分けて、点及び配列単一光子検出器は該検出サブ時間帯で測定しようとする対象を検出し、その光子数と光子が到達すべき時間の次元情報を記録し、最後に各最小時間単位を構成の対象とすることである時間分解単一光子計数イメージングスペクトルシステム。
【請求項６】
前記電気検出ユニットがランダム数発生ユニット、高精度時間間隔測定器、点及び配列単一光子検出器、多チャネルカウンタ、及び遅延器を更に備え、
ランダム数発生ユニットは、DMDマイクロミラー配列を変調するランダム数を発生することに用いられ、且つ該ランダム数発生ユニットは自然界のランダム源を収集してランダム数源として、真のランダム数を処理してランダム数を取得しDMDマイクロミラー配列に出力し、
高精度時間間隔測定器は、検出時間を時間長さ区分することに用いられ、時間座標系における2つの時刻の間の持続時間を記録し、時間次元情報を取得し、精度をピコ秒レベルに制御することができることに用いられ、
点及び配列単一光子検出器は、該検出器が異なる波長に対応するアバランシェダイオードからなり、前記アバランシェダイオードはガイガーモードで作動し、必用に応じてその一部のアバランシェダイオードをイネーブル設置することができ、非周期変化の長い時系列過程において各フレーム内に到達した光子を検出し、及び過度的な周期性過程において各周期に予定された時間帯内に到達した光子を検出し、且つパルス波形で出力することに用いられ、該予定された時間帯は周期持続時間長さのサブセットであり、
多チャネルカウンタは、パルスピーク数を同定及び累計し、点及び配列単一光子検出器が各チャネルで検出した光子数と光子到達時間を記録することに用いられ、
遅延器は、ある過度的な周期内にレーザーパルスに対する制御信号を送信すると同時に、点及び配列単一光子検出器又は高精度時間測定器に対するイネーブルゲート制御信号を送信し又は送信のために一定期間ゲート幅上昇端を遅延させ、該遅延の時間帯が過度的な周期内に対応する点及び単一光子検出器の検出サブ時間帯とすることができ、該時間分解の精度が20psであり、
前記遅延器は時間長さ区分モジュールで代替することができ、前記時間長さ区分モジュールが過度的な周期を複数のサブ時間帯に等分し、且つ各サブ時間帯を点及び配列単一光子検出器及び多チャネルカウンタが検出と計数を行う時間長さ単位とすることに用いられ、該時間分解精度が50psであり、
前記時間長さ区分モジュール又は遅延器は更に多チャネルカウンタ内に設置された時間幅変換器で代替することができ、該時間幅変換器が取得した光子の到達時間を電圧形態に変換して且つ対応のチャネルに記録し、且つ光子到達時間に従って光子数を区分し、統計して1つの周期において各時間間隔内の複数の累積計数を取得することに用いられ、時間分解精度が5psであることを特徴とする請求項５に記載の時間分解単一光子計数イメージングスペクトルシステム。
【請求項７】
光集束収集部材と、前記点及び配列単一光子検出器の光路上に設置された格子分光部材とを更に備え、該格子分光部材が分光する前に凹面鏡を用いて入射光線を視準してそれを平行光とし、分光した後レンズを用いて各波長光をそれぞれ視準して相応チャネルのアバランシェダイオード上に集束し、スペクトル次元情報、即ち各波長での光強度情報を取得することに用いられ、前記スペクトル次元情報はスペクトル波長パラメーターが含まれた画像を構成することに用いられ、
前記点及び配列単一光子検出器は前記格子分光部材の焦平面上に設置されることを特徴とする請求項５に記載の時間分解単一光子計数イメージングスペクトルシステム。
【請求項８】
前記処理ユニットは、更にデータ読み込み記憶器、及びアルゴリズム処理ユニットを備え、
データ読み込み記憶器は、入力した各サブ時間帯、各フレーム又は各イネーブルゲート制御信号時間帯内の光子数及び各光子到達時間の次元情報、相応のランダムマトリックスと検出器の各チャネルの対応の波長情報を記憶することに用いられ、それにより読み込みとキャッシュの役割を果たし、
アルゴリズム処理ユニットは、該ユニットが入力した多チャネルカウンタの計数値、空間光変調復調器を制御するランダムマトリックス、遅延器又は時間長さ区分モジュール又は時間幅変換器で記録した時間次元情報及び検出器の各チャネルの対応の波長情報に基づいて画像を構成することに用いられ、更に以下の8種のイメージングを出力する：
第1種、単一光子二次元イメージングI（x，y）を出力することに用いられる場合、該アルゴリズム処理ユニットは入力計数値、ランダムマトリックスにより、圧縮センシングアルゴリズムを用いて構成し、且つマトリックス充填理論を合わせて画像を後処理し、前記圧縮センシングアルゴリズムは、貪欲構成アルゴリズム、マッチング追跡アルゴリズムMP、直交マッチング追跡アルゴリズムOMP、ベース追跡アルゴリズムBP、LASSO、LARS、GPSR、ベイズ推定アルゴリズム、magic、IST、TV、StOMP、CoSaMP、LBI、SP、l1_ls、smpアルゴリズム、SpaRSAアルゴリズム、TwISTアルゴリズム、l0構成アルゴリズム、l1構成アルゴリズム、l2構成アルゴリズムを含み、
第2種、単一光子二次元イメージングスペクトルI（x，y，l）を出力することに用いられる場合、前記第1種のアルゴリズムのうえに検出器が取得した各チャネルの対応の波長情報を追加し、各波長での光強度空間分布画像を取得し、イメージングスペクトル装置に用いられることができ、また赤、緑、青の三原色での光強度空間分布を取って、カラー画像になることができ、
第3種、単一光子時間分解二次元イメージングI（x，y，t）を出力することに用いられる場合、前記第1種のうえに光子到達時間の次元情報を追加し、各サブ時間帯内の対象に対して画像の構成を行い、時間長さ全体内の画像が動的に変化する過程を構成し、
第4種、単一光子時間分解二次元イメージングスペクトルI（x，y，l，t）を出力することに用いられる場合、前記第2種のアルゴリズムのうえに光子到達時間の次元情報を追加し、各サブ時間帯内の各波長光での対象に対して画像の構成を行い、
第5種、単一光子三次元イメージングI（x，y，z）を出力することに用いられる場合、大規模時間及び蛍光を励起しない測定対象に対して、前記第3種のアルゴリズムのうえに光子が到達する時間間隔を用いて空間位置の光路差を換算し、前記光路差は空間第三次元の距離情報であって、時間分解次元の表示形式及びサブセットであり、予め定められた数の層の異なる空間距離に対応する画像フレームを取得し、
第6種、単一光子三次元イメージングスペクトルI（x，y，z，l）を出力することに用いられる場合、前記第5種のアルゴリズムのうえに点及び単一光子検出器が取得した各チャネルの対応の波長情報を追加し、各波長での光強度空間三次元分布を取得し、
第7種、単一光子時間分解三次元イメージングI（x，y，z，t）を出力することに用いられる場合、前記第5種のアルゴリズムのうえに大規模化した前記空間第三次元の距離情報と過度的な周期のサブ時間帯次元情報を同定且つ区別し、前記大規模化した空間第三次元の距離情報の時間長さが前記過度的な周期のサブ時間帯次元情報よりも大きくし、更に時間分解三次元イメージングを構成し、
第8種、単一光子時間分解三次元イメージングスペクトルI（x，y，z，l，t）を出力することに用いられる場合、前記第6種のアルゴリズムのうえに大規模空間第三次元距離情報と過度的な周期のサブ時間帯次元情報を同定且つ区別し、最終に時間分解三次元イメージングスペクトルを取得する、
ことを特徴とする請求項６に記載の時間分解単一光子計数イメージングスペクトルシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、時間分解単一光子計数イメージングスペクトルシステムに関する。
【背景技術】
【０００２】
本発明は、以前の研究に基づいて改善と革新されたものであり、この分野において、本研究所は既に『単一光子計数イメージングシステム及びその方法』（出願番号又は特許番号：201110103559.3，出願者又は特許権者：中国科学院空間科学と応用研究センター）と『極めて弱い光マルチスペクトルイメージング方法及びそのシステム』（出願番号又は特許番号：201110166471.6，出願者又は特許権者：中国科学院空間科学と応用研究センター）という二つの特許があり、この二つの特許は本所の以前の研究であり、特徴は、二つがいずれも圧縮センシング理論、空間光変調技術及び単一光子検出技術に基づいて、単一光子計数イメージングに用いられることであり、区別は、前のものが単一光子検出器のみを検出素子とし、単一光子レベルの極微弱光対象の二次元イメージングを実現し、後のものが単一光子カウンターラインアレイを検出素子とし、分光光度計を合わせて、単一光子レベルの極微弱光マルチスペクトルカラーイメージングを実現することである。二つの共通の欠点は、静的物体のみをイメージングすることができ、光強度分析、時間分解及び空間三次元分解なしで、アルゴリズムのシミュレーションのみ行い、システムノイズの影響を考慮しておらず、実際の再構成精度が悪く、計数システムとDMDの同期の問題を解決しなくて、まだいくつかの技術的な欠陥が存在する。時間分解問題と多次元パラメータの検出並びにイメージングの問題を解決するために、本発明は時間分解単一光子又は極微弱光多次元イメージングスペクトルシステム及び方法を提出して、以上の一連の欠点を解決する。
【０００３】
いわゆる時間分解とは時間次元での間隔を分解することであり、超高時間分解とは物理と化学の一時的な過程を観察し且つその時間を分解することができることを指し、液相において、多くの物理と化学過程、例えば分子のシス-トランス異性化と方向性のリラクゼーション、電荷とプロトンの移転、励起分子衝突の前期解離、エネルギー伝達と蛍光寿命及び水中の電子の溶媒和などは、ただ10-8秒で完了することができ、ピコ秒の時間分解精度の分析装置のみによりこれらの高速過程をタイムリーに観察することができる。本発明において、大規模又は非常に短い時間間隔内に単一光子を検出して計数サンプリングする。
【０００４】
なお、現在、国際で一時的な極微弱光（例えば蛍光寿命）の測定に用いる技術は主に単一分子の検出技術、時間分解技術及び超解像度測定技術を有する。（１）単一分子検出技術は主に広視野共焦点蛍光顕微鏡技術、走査型近接場光顕微鏡(SNOM)技術、全反射蛍光顕微鏡(TIRF)技術、光学原子間力顕微鏡(AFOM)及びラマン散乱技術を有し、（２）時間分解技術は主に蛍光寿命イメージング(FLIM)、二光子蛍光寿命顕微イメージング、蛍光寿命関連スペクトル(FCS)技術及び多次元蛍光寿命顕微鏡技術を有し、（３）超解像度測定技術は主に誘導放射抑制顕微鏡(STED)技術、位置敏感型顕微鏡(PALM、STORM、dSTORM、GSDIM)技術、光学変動顕微鏡(SOFI)技術及び蛍光共鳴エネルギー転移顕微鏡技術(FRET)を有する。生体高分子に対する蛍光寿命イメージング及び関連スペクトル定量測定方法は、先にFLIM又はFCSシステムを用いて単一点蛍光寿命及び関連スペクトルの測定を行い、そして、レーザビーム走査又はサンプルスキャンシステムを採用して生体高分子蛍光寿命及び関連スペクトルイメージングの測定を行うものである。ナノ変位走査プラットフォームの安定性が悪く、走査過程が複雑で、製造コストを増加しただけでなく、極めてナノ材料と生体高分子の試験時間を延長し、成功率も著しく影響される。ナノ材料高解像度マイクロ構造に対するイメージング測定方法は、常に電子走査顕微鏡を採用してグラフィック特性解析し、高エネルギー電子イオン化が測定されたサンプルを損傷するので、生物活性分子及びナノ材料の非破壊イメージング測定を行うことができない。以上のこれらの技術の共通の問題は同時に観測対象に対して空間蛍光寿命測定と関連スペクトルの分析作業を行うことができない。科技研究の需要が次第に高空間分解、高時間分解、マルチバンド、検出急速、光子励起等の方向へ発展するに従って、これらの機能はますます増長する実際の需要を満たすことができないようになる。
【０００５】
2011年8月Future Markets, Inc.と2011年2月Global Industry Analysts, Inc.の光学顕微の世界市場の分析報告により、2017年まで光学顕微の世界市場の年収益が41億ドルを超えることを予測する。現在と未来の新市場が先進的、高解像度、高精度顕微イメージングに対する需求は不断に増長し、且つ半導体と電子、工業、微小電気機械システム、生物医学と製薬、ナノ技術及びナノ材料の研究開発分野の技術の進歩が大幅に関連計器の市場需要の成長を促進した。一言で言えば、時間分解イメージングスペクトルは魅力的な市場発展の展望を有し、大幅に関連産業の発展を促進することができる。物理と化学研究におけるナノ材料、結晶材料、材料の非破壊検出、量子ドットシステム、フォトニック結晶、高速現象検出、高分解スペクトル測定、量子化学の基礎科学研究、及び生物物理学研究における医療診断、単一分子生物物理、ナノ生物効果、分子バイオニクス、脳機能と認知、プロテオミクス等の生物物理科学の最先端の問題もこのような新規イメージングシステムと方法が必要とされ、多く分野、マルチパラメータ、マルチスケール定量研究の突破を実現する。生命科学、材料科学、化学、エネルギー科学の多くの研究が研究対象に含まれた成分及び成分の変化を了解することが不可欠に必要であり、イメージングスペクトル分析は成分分析を取得する最も効果的な非破壊分析手段の一つであり、研究された成分がどこに分布しているか、どのように成分の変化を発生するか、即ち成分変化の空間分布情報はイメージングスペクトル分析によりこのような研究需求を実現することができ、現在の科学研究発展に従って、成分の空間分布情報を知っている必要するだけではなく、かつ成分が時間に従って変化する過程、特に成分変化過程において中間体があるかどうか、中間体の形成過程、消滅及び作用機構を知っている必要があり、イメージングスペクトルに対して時間分解の要求は増大している。定成分、定量、定時、位置決めの４「定」の目標に達するためには、いずれの技術案から次元問題と感度問題という二つの主な技術的障害がある。従来のイメージングスペクトル技術は必ず一次元スペクトル、二次元平面画像におけるいずれの一次元が、補助走査の方式で実現する。このような作業方式による不利点は、サンプリング時間が必ず走査のために十分な空間を残す必要があり、それだけでなく、平面的センサは大量の感光性検出ユニットからなるので、超高の時間分解目標は電子学でも機器自身でも実現し難い。現在、国際で時間分解イメージングスペクトル機能を有すると公称する計器の実施案は、以下の二つがある：実施案１、時間分解スペクトル装置のうえに二次元空間走査を増加して実現し、該実施案の不利点は、時間分解と空間走査の間の時間の競合により効率低下を引き起こし、全体的性能がより高い科研需求を満たすことができないことであり、実施案２、画像スペクトル装置のうえに時間分解機能を増加し又は単一光子顕微イメージングのうえに光学チャネルを増加して実現し、該実施案の不利点は、時間分解能力を向上し難く、時間分解スペクトルの時間要求を達することができず、応用範囲が非常に有限であることである。
【０００６】
研究対象の光学信号強度の観点から分析すると、スペクトル装置機能に空間分解を導入することは、光センサによりキャプチャ可能なソース信号の強度が低下することを引き起こし、時間分解測定の要求がソース光信号の一層の低下を引き起こし、研究対象のソース信号の強度が低いほど、検出器がキャプチャする必要がある蛍光信号もそれに従って低下し、検出器の感度が時間分解イメージングスペクトル装置の主な性能であり、検出器の感度が高いほど、装置の性能も高くなり、単一光子検出器は光学検出のリミットである光最小エネルギーを検出することができ、単一光子検出技術はずっと時間分解イメージングスペクトル装置の究極の追求である。
【０００７】
原理と従来の検出技術の条件制限により、新原理、新方法のみにより本当の意味での時間分解単一光子多次元イメージングシステムを実現することができ、本発明は数学的研究の最新結果‐圧縮センシング理論を採用して、現代の洗練された技術的条件を合わせて、一つの単一光子点検出器で二次元空間の分解測定を実現し、勿論、このような方法は二次元空間分布での一次元を節約し、検出感度でも、波長範囲でも表面要素検出器と比較するとより多くの明らかな利点があり、感度が非常に高く、スループット測定が高く、信号対雑音比が高く、波長の相応範囲が広く且つコストが低く、点検出器を用いて単一光子計数イメージングを実現し、必ず未来の単一光子レベルイメージングの重要な発展傾向になる。単一光子計数イメージングは極微弱光検出技術であり、光子計数及び光子が検出された確率特徴付け光強度分布を記録することにより画像を反転する。
【０００８】
ラインアレイ又は配列ガイガーモードの単一光子検出器で、イメージングスペクトルを実現することができ、時間分解でピコ秒分解精度まで実現することができ、感度で単一光子検出レベルに達し、空間分解でナノスケールに達し、本当の意味での単一光子時間分解イメージングスペクトルを実現し、従来のイメージングスペクトルと本質的に異なる。前記イメージングスペクトルは精確的なスペクトル情報を取得と表示する重要な技術であり、スペクトル画像はより多いスペクトル情報を含み、且つマルチスペクトルイメージング技術は良好的にメタメリズム現象を克服する。前記単一光子が極微弱光であり、光の不可分のエネルギーの最小単位と考えられる。該離散パルス信号が検出可能なリミットであり、通常単一光子検出器を用いて検出する。計数型の単一光子検出器は飽和状態で作動すると、感度が単一光子レベルに達することができ、統計的方法を採用して光子密度画像を取得し、未飽和状態で作動すると、出力した電気信号の振幅が検出した光子数の変化に従って変化し、この電気信号により極微弱光画像を取得する。該単一光子計数方法はパルススクリーニング技術とデジタル計数技術を用いて極めて微弱な自然離散信号を識別して抽出し、不安定要素による影響が小さく、検出器の熱ノイズの大部分の影響を除去することができ、それにより大幅に信号対雑音比を向上し、デジタル信号を出力する。
【０００９】
本発明の理論基礎は圧縮センシング理論と空間光変調器技術であり、従来の技術的圧縮センシング理論と空間光変調器技術は具体的に以下のように説明する：
【００１０】
圧縮センシング原理はDonoho、Tao及びCandesなどで提出した１つの新たな数学理論であり、ランダムサンプリングの方式、一層少ないデータサンプリング数（ナイキスト／シャノンのサンプリング定理のリミットよりもはるかに低い）で完全にオリジナル信号を回復することを実現し、且つより高いロバスト性を有する。該原理は主に圧縮サンプリングステップ、スパース変換ステップ及びアルゴリズム再構成ステップに分けられることができる。圧縮サンプリングは、測定された信号が高次元から低次元へマッピングと収集の過程であり、前記のスパース変換は適切なΨを選択して、ｘがΨ変化を経て取得した値ｘ’がスパースであり、即ちｘがΨフレームでスパース表現することができ、前記のアルゴリズム再構成は観測データｙ、測定マトリックスＡ及びフレームΨが既知の条件でｙ＝ＡΨｘ’＋ｅを解く過程であり、最後に更に、
【数１】
0006002232.tif 000002
でｘを反転取得する。
【００１１】
前記空間光変調器（Spatial Light Modulator，SLMと略称する）は情報を一次元又は二次元の光学データフィールド上に載せて且つリアルタイムに光学情報処理する機器であり、時間に従って変化する電気駆動信号又は他の信号の制御で、空間での光分布の振幅又は強度、相位、偏光状態及び波長を変化し、又はインコヒーレント光をコヒーレント光に変換する。その最も典型的な代表はデジタルマイクロミラーデバイス（Digital Micro-mirror Device，DMDと略称する）であり、これは世界で最も精密的な光スイッチである。DMDの核心は数千のヒンジに設けられたマイクロミラーからなるマイクロミラー配列（主なDMDが1024×768の配列で構成され、最大で2048×1152に達することができる）であり、各ミラーのサイズが14μm×14μm（又は16μm×16μm）であり且つ１つの画素の光をオンオフすることができ、これらのマイクロミラーはすべて浮遊して、各ミラー下の記憶ユニットに対してすべてバイナリ平面信号で電子化アドレッシングを行うことにより、各ミラーを静電的に両側に10〜12°ほど（ここでは、+12°と-12°を取る）傾斜させ、この２種の状態を1及び0とそれぞれ標記し、それぞれ「オン」と「オフ」に対応し、ミラーが作動しない際、それらは0°の「停止」状態になる。
【発明の概要】
【発明が解決しようとする課題】
【００１２】
本発明は、ナノスケールとピコ秒スケールで、定成分、定量、定時、位置決めのように研究対象に含まれた成分及び成分の運動変化規律を「見る」且つ分析し、次元の増加による次元問題と感度問題を解決し、可視光と近赤外光の点及び配列単一光子検出技術、時間分解測定技術、空間光変調技術、圧縮センシング理論、単一光子計数イメージング技術、イメージングスペクトル技術、再構成アルゴリズム技術を合わせて、極微弱光単一光子レベルで物体を反射、散乱又は透過し、空間二次元、光強度、時間分解、スペクトル分解の５次元パラメータ情報を有するイメージングスペクトルを測定することを実現し、更に時間分解単一光子計数イメージングスペクトルシステムを提供し、国際関連分野の発展を促進するようになることを目的とする。
【００１３】
本発明は、時間分解単一光子計数イメージングスペクトルシステムを提出し、光源によって触発されるトリガー、光学イメージングシステム、DMDマイクロミラー配列、光学集束収集システム、光減衰器、単一光子計数器、ドライブ制御モジュール及び最適化アルゴリズムモジュールを備え、圧縮センシング原理に基づいて前記トリガーを触発する極微弱光源を採用して物体を照射して得られたイメージングスペクトルを秒レベル分解精度にし、更に時間に従って動的に変化する物体をイメージングし、時系列に従って配列された連続グレースケールビデオ画像フレームを出力する時間分解単一光子計数イメージングスペクトルシステムであって、
前記トリガーはその前端に位置する極微弱光触発源で触発されて、該トリガーの出力端がドライブ制御モジュールの入力端と接続されて、トリガーが触発された際ドライブ制御モジュールはドライブ制御信号を出力してその出力端と接続されたDMDマイクロミラー配列と単一光子計数器を作動し始めるように触発し、該DMDマイクロミラー配列が反転し始め、単一光子計数器が同時に計数し始め、光減衰器は光強度を減衰することに用いられ、光学集束収集システムは光線を集束収集することに用いられ、単一光子計数器は光を単一光子計数するためのものであり、
前記単一光子計数器の出力端が前記最適化アルゴリズムモジュールの１つの入力端と接続され、単一光子計数器の計数値を最適化アルゴリズムモジュールの１つの入力パラメータとし、該ドライブ制御モジュールが記憶する選定領域のランダム測定マトリックスを受信してそのもう１つの入力パラメータとするために、前記最適化アルゴリズムモジュールのもう１つの入力端が前記ドライブ制御モジュールの一出力端と接続されて、前記最適化アルゴリズムモジュールは入力した単一光子計数器の測定値により前記ランダム測定マトリックスとともにスパース信号を構成し、且つ前記スパース信号を変換して反転し、二次元光子密度により構成した画像を取得し、即ち画像の各画素の値はDMDマイクロミラー配列の対応する画素に入射した光子密度に等しく、その内、光子密度とは単位時間の単位面積に前記DMDマイクロミラー配列の反射波長範囲内で入射される光子数であり、Ｍ個のｔ時間間隔を経て、時系列に従って配列された時間分解可能な一連の二次元グレースケール画像ビデオフレームを構成することができ、
前記DMDマイクロミラー配列はデジタルマイクロミラーデバイスであり、
前記DMDマイクロミラー配列の反転をドライブ制御するために、前記ドライブ制御モジュールのもう１つの出力端が前記DMDマイクロミラー配列の入力端と接続されて、該ドライブ制御モジュールはDLPデジタル光処理技術に基づいたものであり、該ドライブ制御モジュールがDMDマイクロミラー配列の作動領域を選定した後、ランダム測定マトリックスをダウンロードして前記DMDマイクロミラー配列の反転をドライブ制御し、前記DMDマイクロミラー配列が反転すると同時に前記単一光子計数器へ同期信号を送信して、DMDマイクロミラー配列と単一光子計数器の間の同期を保証し、即ちDMDマイクロミラー配列が一回反転すると、単一光子計数器が該回の反転の時間間隔内の単一光子数を累積計数し、DMDマイクロミラー配列が反転を完成した後、単一光子計数器がクリアして累積計数を再び開始し、すべての計数と該選定領域のランダム測定マトリックスがいずれも最適化アルゴリズムモジュールに伝送される。
【００１４】
上記技術案において、カラーイメージングに用いられる際、単一光子計数器に代えて単一光子検出器ラインアレイを用い、且つ前記単一光子検出器ラインアレイと光学集束収集システムの間の光路上に分光光度計を設置する。
【００１５】
上記技術案において、前記光学イメージングシステムと光学集束収集システムはいずれも光学レンズ群を採用し、それぞれ光学イメージングと光学集束を担当し、極微弱光は光学イメージングシステムを通過した後、DMDマイクロミラー配列上に同じ大きさ又は縮小又は拡大の画像を形成することができ、実際の要求に従ってイメージング調整する。
【００１６】
上記技術において、前記光学集束収集システムはファイバー結合技術を用いて、即ち前記分光光度計で分光された光束をファイバーに結合し、ファイバー結合技術を用いて分光がそれぞれ対応の単一光子計数器上に収集され、
前記分光光度計は、光視準部、分光部、角度測定部、光度観察及び測定部を備え、スペクトルの分析と測定に用いられ、
前記分光光度計は、プリズム分光光度計又は格子分光光度計から成る。
【００１７】
本発明は、これに基づいて時間分解単一光子計数イメージングスペクトルシステムに発展し、光源、イメージングスペクトル測定ユニット、電気検出ユニット、システム制御ユニット、処理ユニットを備え、圧縮センシング原理に基づいて、レーザーパルスにより物体蛍光を励起し、光強度分布とイメージングスペクトルを取得し、且つ時間分解方策を用いてピコ秒レベルまでの時間分解精度を実現し、更に蛍光寿命の過度的な周期性物質を検出することに用いられる時間分解単一光子計数イメージングスペクトルシステムであって、
前記光源は、システム制御ユニットで送信した触発パルスの触発で、レーザーパルスを発生し、測定しようとする物体に当てて、該物体は該物体に関する情報（以下「物体情報」ということがある。）が載せられた蛍光を発生し、
前記イメージングスペクトル測定ユニットは、光学イメージングシステム、DMDマイクロミラー配列、光集束収集部材を更に備え、測定しようとする物体情報が載せられた蛍光を前記DMDマイクロミラー配列にイメージングし、該DMDマイクロミラー配列はランダム光変調を採用して且つその上のイメージングを反射し、且つ反射光は光路上に配布された光集束収集部材により視準した後電気検出ユニットに入射することに用いられ、
前記電気検出ユニットは、時間分解方策に従って入射光を検出し且つ入射光の光子数を記録し、且つ各光子が到達する時間次元情報を記録することに用いられ、
前記処理ユニットは、入力したすべての各時間帯の光子数、光子の到達時間の次元情報及びランダム光変調マトリックス情報により、圧縮センシングとイメージングスペクトルアルゴリズムを用いてマルチパラメーター情報のイメージングスペクトル構成を実現し、且つ8種類のイメージングを出力することに用いられ、
前記システム制御ユニットは、各部材のイネーブルを備え、即ちドライブイネーブル信号を上記各部材に送信してそれを正常な動作を開始するようにし、且つ光源触発パルスと光子計数に対する同期を実現することに用いられ、且つマイクロミラー配列の反転とその上のランダムマトリックスの交換を制御し且つ過度的な周期の開始前又は終了後に反転をすることに用いられ、
前記DMDマイクロミラー配列はランダム数発生器によりベルヌーイバイナリランダムマトリックスをロードして、光に対するランダム光変調を実現し、前記時間分解方策は非周期変化の長い時系列過程に対してフレームごとに測定の方法を採用し、周期変化特性を有する過度的な過程に対して、最適化の時間分解方策を採用し、該最適化の時間分解方策は、具体的に、過度的な過程の周期を予め定められた数の検出サブ時間帯に分けて、点及び配列単一光子検出器は該検出サブ時間帯で測定しようとする対象を検出し、その光子数と光子が到達すべき時間の次元情報を記録し、最後に各最小時間単位を構成の対象とすることである。
【００１８】
前記の時間分解単一光子計数イメージングスペクトルシステムは、前記電気検出ユニットがランダム数発生ユニット、高精度時間間隔測定器、点及び配列単一光子検出器、多チャネルカウンタ、及び遅延器を更に備え、
ランダム数発生ユニットは、DMDマイクロミラー配列を変調するランダム数を発生することに用いられ、且つ該ランダム数発生ユニットは自然界のランダム源を収集してランダム数源として、真のランダム数を処理してランダム数を取得しDMDマイクロミラー配列に出力し、
高精度時間間隔測定器は、検出時間を時間長さ区分することに用いられ、時間座標系における2つの時刻の間の持続時間を記録し、時間次元情報を取得し、精度をピコ秒レベルに制御することができることに用いられ、
点及び配列単一光子検出器は、該検出器が異なる波長に対応するアバランシェダイオードからなり、前記アバランシェダイオードはガイガーモードで作動し、必用に応じてその一部のアバランシェダイオードをイネーブル設置することができ、非周期変化の長い時系列過程において各フレーム内に到達した光子を検出し、及び過度的な周期性過程において各周期に予定された時間帯内に到達した光子を検出し、且つパルス波形で出力することに用いられ、該予定された時間帯は周期持続時間長さのサブセットであり、
多チャネルカウンタは、パルスピーク数を同定及び累計し、点及び配列単一光子検出器が各チャネルで検出した光子数と光子到達時間を記録することに用いられ、
遅延器は、ある過度的な周期内にレーザーパルスに対する制御信号を送信すると同時に、点及び配列単一光子検出器又は高精度時間測定器に対するイネーブルゲート制御信号を送信し又は送信のために一定期間ゲート幅上昇端を遅延させ、該遅延の時間帯が過度的な周期内に対応する点及び単一光子検出器の検出サブ時間帯とすることができ、該時間分解の精度が20psであり、
前記遅延器は時間長さ区分モジュールで代替することができ、前記時間長さ区分モジュールが過度的な周期を複数のサブ時間帯に等分し、且つ各サブ時間帯を点及び配列単一光子検出器及び多チャネルカウンタが検出と計数を行う時間長さ単位とすることに用いられ、該時間分解精度が50psであり、
前記時間長さ区分モジュール又は遅延器は更に多チャネルカウンタ内に設置された時間幅変換器で代替することができ、該時間幅変換器が取得した光子の到達時間を電圧形態に変換して且つ対応のチャネルに記録し、且つ光子到達時間に従って光子数を区分し、統計して1つの周期において各時間間隔内の複数の累積計数を取得することに用いられ、時間分解精度が5psであることを特徴とする。
【００１９】
前記時間分解単一光子計数イメージングスペクトルシステムは、光集束収集部材と、前記点及び配列単一光子検出器の光路上に設置された格子分光部材とを更に備え、該格子分光部材が分光する前に凹面鏡を用いて入射光線を視準してそれを平行光とし、分光した後レンズを用いて各波長光をそれぞれ視準して相応チャネルのアバランシェダイオード上に集束し、スペクトル次元情報、即ち各波長での光強度情報を取得することに用いられ、前記スペクトル次元情報はスペクトル波長パラメーターが含まれた画像を構成することに用いられ、且つ格子分光にゼロ次回折、１次回折、２次回折……があり、前記ゼロ次回折が原像であり、残りの次の回折がすべてスペクトルであり、１次回折スペクトルが最も明るく、１次回折スペクトルを検出、又は高次回折スペクトルをそれぞれ検出して更に格子回折公式を用いて全体の光強度波長分布を取得し、
前記点及び配列単一光子検出器は前記格子分光部材の焦平面上に設置される。
【００２０】
上記技術案において、前記処理ユニットは、更にデータ読み込み記憶器、及びアルゴリズム処理ユニットを備え、
データ読み込み記憶器は、入力した各サブ時間帯、各フレーム又は各イネーブルゲート制御信号時間帯内の光子数及び各光子到達時間の次元情報、相応のランダムマトリックスと検出器の各チャネルの対応の波長情報を記憶することに用いられ、それにより読み込みとキャッシュの役割を果たし、
アルゴリズム処理ユニットは、該ユニットが入力した多チャネルカウンタの計数値、空間光変調復調器を制御するランダムマトリックス、遅延器又は時間長さ区分モジュール又は時間幅変換器で記録した時間次元情報及び検出器の各チャネルの対応の波長情報に基づいて画像を構成することに用いられ、更に以下の8種のイメージングを出力する：
（１）単一光子二次元イメージングI（x，y）を出力することに用いられる場合、該アルゴリズム処理ユニットは入力計数値、ランダムマトリックスにより、圧縮センシングアルゴリズムを用いて構成し、且つマトリックス充填理論を合わせて画像を後処理し、前記圧縮センシングアルゴリズムは、貪欲構成アルゴリズム、マッチング追跡アルゴリズムMP、直交マッチング追跡アルゴリズムOMP、ベース追跡アルゴリズムBP、LASSO、LARS、GPSR、ベイズ推定アルゴリズム、magic、IST、TV、StOMP、CoSaMP、LBI、SP、l1_ls、smpアルゴリズム、SpaRSAアルゴリズム、TwISTアルゴリズム、l0構成アルゴリズム、l1構成アルゴリズム、l2構成アルゴリズムを含み、
（２）単一光子二次元イメージングスペクトルI（x，y，l）を出力することに用いられる場合、前記第1種のアルゴリズムのうえに検出器が取得した各チャネルの対応の波長情報を追加し、各波長での光強度空間分布画像を取得し、イメージングスペクトル装置に用いられることができ、また赤、緑、青の三原色での光強度空間分布を取って、カラー画像になることができ、
（３）単一光子時間分解二次元イメージングI（x，y，t）を出力することに用いられる場合、前記第1種のうえに光子到達時間の次元情報を追加し、各サブ時間帯内の対象に対して画像の構成を行い、時間長さ全体内の画像が動的に変化する過程を構成し、
（４）単一光子時間分解二次元イメージングスペクトルI（x，y，l，t）を出力することに用いられる場合、前記第2種のアルゴリズムのうえに光子到達時間の次元情報を追加し、各サブ時間帯内の各波長光での対象に対して画像の構成を行い、
（５）単一光子三次元イメージングI（x，y，z）を出力することに用いられる場合、大規模時間及び蛍光を励起しない測定対象に対して、前記第3種のアルゴリズムのうえに光子が到達する時間間隔を用いて空間位置の光路差を換算し、前記光路差は空間第三次元の距離情報であって、時間分解次元の表示形式及びサブセットであり、予め定められた数の層の異なる空間距離に対応する画像フレームを取得し、
（６）単一光子三次元イメージングスペクトルI（x，y，z，l）を出力することに用いられる場合、前記第5種のアルゴリズムのうえに点及び単一光子検出器が取得した各チャネルの対応の波長情報を追加し、各波長での光強度空間三次元分布を取得し、
（７）単一光子時間分解三次元イメージングI（x，y，z，t）を出力することに用いられる場合、前記第5種のアルゴリズムのうえに大規模化した前記空間第三次元の距離情報と過度的な周期のサブ時間帯次元情報を同定且つ区別し、前記大規模化した空間第三次元の距離情報の時間長さが前記過度的な周期のサブ時間帯次元情報よりも大きくし、更に時間分解三次元イメージングを構成し、
（８）単一光子時間分解三次元イメージングスペクトルI（x，y，z，l，t）を出力することに用いられる場合、前記第6種のアルゴリズムのうえに大規模空間第三次元距離情報と過度的な周期のサブ時間帯次元情報を同定且つ区別し、最終に時間分解三次元イメージングスペクトルを取得する。
【００２１】
前記圧縮センシングの該原理は主に圧縮サンプリングステップ、スパース変換ステップ及びアルゴリズム再構成ステップに分けられる。圧縮サンプリングは、測定された信号が高次元から低次元へマッピングと収集する過程であり、前記のスパース変換は適切なスパースフレームを選択する過程であり、前記のアルゴリズム再構成は観測データｙ、測定マトリックスＡ及びフレームΨが既知の条件で信号を回復する過程である。
【００２２】
圧縮センシング理論により、マイクロミラー配列が圧縮可能な観測対象情報が載せられた光をランダム反射する必要があり、マイクロミラー配列における単一のマイクロミラーが+12°で反転する際反射した光が点及び配列単一光子検出器で受信し、マイクロミラー配列における単一のマイクロミラーが-12°で反転する際反射光が点及び配列単一光子検出器で受信することができなくなり、それにより、測定された信号に対する圧縮サンプリングを完成し、同時にマイクロミラー配列のシェーディングマトリックスの真ランダム性を確保し、更に光集束収集部材に反射される確率がランダムとなるように極微弱光を制御する。
数学モデルは以下のようである：
ｘ∈ＲNが測定されたデータで、ｙ∈ＲKが観測データで、Ａ∈ＲRxNがランダム投影マトリックス（Ｋ＝Ｎ）で、ｅ∈ＲKがシステムノイズで、Ｋがｘにおける非ゼロ要素の数、いわゆるスパース度であると、圧縮サンプリングの過程は以下の（１）式で表現することができる：
【数２】
0006002232.tif 000003
ｘが圧縮可能な又はスパース表現可能なものであると、
【数３】
0006002232.tif 000004
であり、Ψ＝｛ψ1，ψ2，…，ψN｝がスパース変換マトリックスであり、（１）式は（２）式に変形される：
【数４】
0006002232.tif 000005
ＡΨがRIP基準を満たし、ＡがΨと関連しないほど、サンプリングが必要する測定数Ｋは小さくなり、計算量は小さくなり、このため、本発明において、（２）式においてΨがフレームであり、Ａがランダム測定マトリックスである。
【００２３】
二次元画像の画素総数をＮの場合、（１）式における測定マトリックスがＡ＝［ａ1，ａ2，…，ａN］であり、ａiがＡの第ｉ行であり、ｐ×ｑ画素の二次元画像の列をエンドツーエンド接続し、Ｎ×１（Ｎ＝ｐ×ｑ）の一次元列ベクトルになり、（１）式におけるｘに対応し、その各要素が相応位置の光子強度を表す。毎回の反転の際、マイクロミラー配列の行をエンドツーエンド接続し、１×Ｎの一次元行ベクトルになり、測定マトリックスＡにおける一行に対応し、その各要素は相応位置のマイクロミラーが光集束収集部材の主軸方向へ反転するかどうかを表し、測定マトリックスＡが合計Ｍ次元であり、即ちＭ行Ｎ列のマトリックスである。
【００２４】
システム制御ユニットはマイクロミラー配列をランダム反転し始め、毎回の点及び配列単一光子検出器が検出した光子数をｎと記し、それは光子強度画像とマイクロミラー配列のランダム測定配列の間の内積値に相当し、（１）式における観察ベクトルｙの１つの要素に対応し、
【数５】
0006002232.tif 000006
式において、Ｋ回の測定を繰り返し、全組の観測データｙ（ｙがＫ×１の一次元列ベクトルである）を取得することができる。
【００２５】
前記のスパース再構成は、観測データｙと測定マトリックスＡが既知の条件で（１）式におけるｘを解くものであり、これは１つのill-posed問題であり、一般的に圧縮センシングアルゴリズムを用いて解き、該アルゴリズムは以上のように多種のタイプに分けられ、例示として、そのうちの１種の一般的な表現方式を取り、（３）式と表現する：
【数６】
0006002232.tif 000007
||…||pがノルム演算子を表し、
【数７】
0006002232.tif 000008
であり、一項が制約付き最小二乗であり、ｆ（ｘ）と記し、第二項がｘスパース度に対する制約であり、２項の和が目標関数である。Ｍ≦Ｏ（Ｋ・log（Ｎ／Ｋ））回の測定のみで、完全に原信号を再構成することができる。
【００２６】
上記は光強度次元情報により単色光の画像を再構成するものであり、スペクトル次元情報を加えると、即ち各波長の信号をそれぞれ再構成し、スペクトル分析することができ、時間次元情報を加えると、即ち上記のうえに前記時間分解測定方法に従って時間間隔を分けて信号の再構成を行うことができ、時間分解を実現することができ、時間分解において取得した光子到達時間は、一定の条件の下で空間距離次元情報に変換することができ、即ち信号が距離に従って分類区分され、アルゴリズムが各距離での平面画像を再構成し、最後に各距離での再構成画像を一体に重ね合わせて、空間三次元立体画像を形成する。
【００２７】
本発明は、時間分解単一光子計数イメージングスペクトル方法を更に提出し、該方法が上記の時間分解単一光子計数イメージングスペクトルシステムに基づいたものであり、前記方法は、
極微弱光単一光子の条件の下で、物体が励起されて又は自己励起されることにより放出、反射、散乱、透過又は屈折した光がレンズ群を通してDMDマイクロミラー配列にイメージングし、
圧縮センシング理論により前記DMDマイクロミラー配列でのランダム変調ベースを設置し、且つその反射光を点又は配列単一光子検出器に収集し、時間分解の方策を採用して、圧縮センシングアルゴリズムとイメージングスペクトルアルゴリズムにより5次元パラメーター情報を有するイメージングスペクトルの構成を実現し、最終にナノスケールとピコ秒スケールでの単一光子時間分解イメージングスペクトルを実現し、
前記時間分解の測定方策は以下の3種の方法の中のいずれを採用することができる：
方策一、過度的な過程の周期が1.5ms〜5msである場合、光源をオンし、過度的な周期がTであると、該時間周期をd個の時間間隔に等分し、それぞれｔ1，ｔ2，ｔ3，…，ｔdと記し、該周期T内にDMDマイクロミラー配列でのランダム変調ベースマトリックスが不変することを保持し、点又は配列単一光子検出器はそれぞれｔi時間間隔内にある単一光子を検出し、ｉ＝１，２，…，ｄであり、多チャネルカウンタはそれぞれ各時間間隔内の単一光子数を記録し、且つ高精度時間測定器が記録したタイムコードと一体に合わせて1つのデータパケットとし、このように各計数の対応の時間間隔を取得することができ、同時に次回のレーザーパルスを発射する前に、即ち前の周期Tの対応のd個の時間間隔が丁度完全に測定し終わる時刻に、DMDマイクロミラー配列は瞬間に次のフレームに反転し、全部の上記操作をし、P回の反転を繰り返して、最終に各ｔi時間間隔が相応的にP個の計数を有し、それぞれP個のランダムマトリックスに対応し、この一対一の関係により、それぞれこのd個の時間間隔に対してアルゴリズムの構成を行い、1つの過度的な周期内のスペクトル強度変化過程を反転して取得することができ、光強度が極めて弱くすると、多回で測定して累計して相応の計数を大きくし、そのように操作し、
方策二、過度的な過程の周期が80ns〜1.5msである場合、光源に触発パルスを伝送する同時に、イネーブル制御信号を各部材に送信し、該信号が該点又は配列単一光子検出器と高精度時間測定器に到達する前に遅延器を経過する必要があり、具体的に、1）まず、DMDマイクロミラー配列が1つのフレームを固定して不変することを保持し、点又は配列単一光子検出器のゲート幅の開始端が過度的な周期の開始時刻と重なり合うことを保持し、ゲート幅が過度的な周期Tよりも小さくし、周期が開始する際、点又は配列単一光子検出器と高精度時間測定器は同時に測定し始め、該過度的な周期内に一回のみを検出し、測定した計数値がゲート幅と該過度的な周期の重なり時間帯内の単一光子数であり、順次にQ回を繰り返して、そして、統計原理に基づいて毎回の計数を合計して、その後、デジタル遅延器を用いてゲート幅に20psを加え、上記ステップに従って同様に1つの計数合計を取得することができ、第1個の合計を参考値とし、第2個の合計と第1個の合計の差を延長したゲート幅時間内の統計計数とし、この方法に従って参考点時刻から蛍光寿命終了時刻までの間のd個の部分の統計計数を取得することができ、2）また、ゲート幅が不変することを保持すると、ゲート幅の到達時刻が早めになり、以上の方法に従って取得したシリーズの計数差値が蛍光寿命開始時刻から参考点時刻までの間のd個の部分の統計計数であり、上記の2つの点は、すべて蛍光寿命全体内の各部分の統計計数を取得するものであり、DMDマイクロミラー配列が一回反転すると、全部の上記操作をし、P回の反転を繰り返して、この際の各ｔi時間間隔は相応的にP個の計数を有し、それぞれP個のランダムマトリックスに対応し、この一対一の関係により、それぞれこのd個の時間間隔に対してアルゴリズムの構成を行い、1つの周期内の一時的な過程を反転して取得することができ、光強度が極めて弱くすると、多回で測定して累計して相応の計数を大きくし、そのように操作し、
方策三、まず、DMDマイクロミラー配列が1つのフレームを固定して不変し、且つ光源に送信する触発パルスを時間幅変換器の参考パルスとし、そして時間幅変換器は光子を取得する時間を電圧形態で記録し、対応の時間チャネルに記録し、且つ光子が到達する時間に従って光子数を区分し、1つの周期内に各時間間隔内のd個の部分の累積計数を統計し、DMDマイクロミラー配列が一回反転すると、全部の上記操作をし、P回の反転を繰り返して、この際各ｔi時間間隔が相応的にP個の計数を有し、それぞれP個のランダムマトリックスに対応し、この一対一の関係により、それぞれこのd個の時間間隔に対してアルゴリズムの構成を行い、1つの周期内の過度的な過程を反転して取得することができ、光強度が極めて弱くすると、多回で測定して累計して相応の計数を大きくし、そのように操作する。
【００２８】
上記技術案において、前記単一光子検出器は、光電子増倍管、ICCD、EMCCD、線形モード又はガイガーモード単一光子検出器であってよく、
前記の光源は赤外又は可視光を選択使用し、且つ前記レーザーが適切な波長のフェムト秒又はナノ秒パルスレーザーを選択使用する。
【００２９】
上記技術案において、空間光変調器の正負級光路上に同時に収集検出部材が設置されて且つ同時測定をし、正負のランダム測定マトリックスの相互関連性を利用して、関連回復アルゴリズムを採用して、再構成画像の品質を向上させ及び回復アルゴリズムの実行時間を短縮することができる。
【００３０】
一言で言えば、本発明に記載の時間分解単一光子計数イメージングスペクトル方法は、蛍光寿命イメージング(FLIM)システムに単一光子計数圧縮センシングイメージング技術を導入し、現在、国際での蛍光寿命イメージングと蛍光イメージングスペクトルシステムは一般的に光電子増倍管、ICCD、EMCCDなどの高精度光検出技術を採用し、ナノ変位プラットフォーム又はレーザビーム走査技術を合わせて、本発明は、初めて圧縮センシング理論、単一光子検出技術及び蛍光寿命測定技術を組み合わせて、走査方式が蛍光寿命イメージングの安定性に対する影響を考えて、極微弱光映像をランダムマイクロミラー配列にマッピングし、更にその反射光を点及び配列単一光子検出器に集束し、複数の部分に区分し統計して計数し、更に該区分計数シーケンスを圧縮センシングすることにより再構成画像を形成し、同時に５次元パラメータを有する時間分解単一光子計数イメージングスペクトルを取得することができる。
【００３１】
従来の技術と比べて、本発明の利点は、本発明に記載のシステムが可視光と近赤外光の点及び配列単一光子検出技術、時間分解測定技術、空間光変調技術、圧縮センシング理論、単一光子計数イメージング技術、イメージングスペクトル技術、再構成アルゴリズム技術を組み合わせて、ナノ精度の空間二次元、光強度、ピコ秒精度の時間分解、スペクトル分解の５次元パラメータ情報を有するイメージングスペクトル測定を実現することである。
【００３２】
本発明の主な技術的な問題と革新は、以下を含む：
（１）可視光と近赤外光の点及び配列単一光子検出技術を用いて、取得した単一光子計数データを時間に従ってグループ区分し、圧縮センシング理論を用いて、グループを再構成して単一光子時間分解イメージングを実現する。
（２）可視光と近赤外光の点及び配列単一光子検出技術とイメージングスペクトル技術を用いて、それぞれ時間、スペクトル部分に従って単一光子計数データのグループ区分及び計算をし、単一光子時間分解イメージングスペクトルを実現し、イメージングとスペクトルは同時に測定することができる。
（３）空間距離次元は時間分解次元のサブセットであるので、時間分解で取得した光子到達時間に従って空間距離での光路差を推定し、空間第三次元の次元情報を取得し、三次元立体イメージングを実現することができ、縦方向の解像度がミリメートルレベルに達し、平面解像度がナノレベルに達することができる。
（４）一時的な周期対象の測定に対して３種の新しい圧縮センシング原理に基づいた時間分解測定方法を提出する：時間間隔測定に基づく時間分解方法、遅延測定に基づく時間分解方法及び光子到達時間に基づく時間分解方法であり、それにより時間分解精度がピコ秒レベルに達する。
（５）圧縮センシングで、適切なフレームを選択して、マトリックス填充の考えを導入してイメージング品質を向上させる。
（６）８種のイメージングに対して、圧縮センシング理論を組み合わせて、相応のアルゴリズムを設計する。
【００３３】
以上の利点と革新に基づいて、本発明が研究したシステムは単一分子生物物理学、材料欠陥の検出、ナノ材料、マイクロエレクトロニクス、量子ドット、生命科学及び新エネルギー光電変換材料などの多くの新たなハイテク産業分野に広く応用することができる。
【図面の簡単な説明】
【００３４】
【図１】図１は、本発明による時間分解単一光子計数イメージングスペクトルシステムの構造模式図である。
【図２】図２は、本発明による観測対象が自発光物体である際の時間分解単一光子計数イメージングスペクトルシステムの構造模式図である。
【図３】図３は、本発明による圧縮センシング理論に基づいた時間分解単一光子計数イメージングスペクトルシステムの構造模式図である。
【図４】図４は、本発明による観測対象が自発光物体である際の時間分解単一光子計数イメージングスペクトルシステムの構造模式図である。
【図５】図５は、本発明が提出した時間分解単一光子計数イメージングスペクトルシステムの全体構造図である。
【図６】図６は、本発明による空間光変調器がDMDを採用する際その上の単一のマイクロミラーの反射機構原理図である。
【図７】図７は、本発明実施例の正負画像構造の簡素化原理図である。
【図８】図８は、本発明実施例が採用した時間分解測定方法のブロック図である。
【図９】図９は、本発明による時間間隔を分けて単一光子数を記録する模式図である。
【図１０】図１０は本発明実施例が採用したもう１つの時間分解測定方法のブロック図である。
【図１１】図１１は本発明が統計して差を取得した後時間間隔内の単一光子の計数の模式図である。
【発明を実施するための形態】
【００３５】
以下、図面と実施例を参照して本発明を更に説明する。
【００３６】
本発明技術案が提出した時間分解単一光子計数二次元イメージングシステムは、該システムがトリガー触発技術、圧縮センシング理論、空間光変調技術、ファイバー結合技術及び単一光子検出技術に基づいて、大規模時間と極微弱光単一光子のレベルの下で、革新的にトリガー１８を採用して触発し、且つ毎回の触発後、検出して計数し始め、空間光変調技術を用いて画像信号をランダムに投影し、ランダムの光強度信号に変換し、単一光子検出器１６を検出素子とし、ｔ時間間隔（ｔが秒レベルである）毎に１組データ（１組データが１つの画像に対応する）をサンプリングし、累計の光子数を測定値として、最後に圧縮センシング及び関連アルゴリズム１５を用いてそれに対して再構成を行い、先に集中にサンプリングし後にバッチで再構成して計算する方式を採用して、単一光子レベルの極微弱光二次元イメージングの時間分解を実現し、大規模時間で点検出器を用いて測定しようとする物体の動的な変化を観測することができることを特徴とする。
【００３７】
図１と図２に示すように、本発明に記載の時間分解単一光子計数二次元イメージングシステムは主に異なる波長のレーザー又はハロゲンランプ１、トリガー１８、フィルタとアッテネータ２、光学イメージング部材３、DMD４、光集束収集部材７、単一光子検出器１６、計数器１７、システム制御プラットフォーム１３と圧縮センシング及び関連アルゴリズム１５からなる。
【００３８】
なお、先に単一光子検出器１６をイネーブルし、それを持続に検出し、新しく設置されたトリガー１８がシステム制御プラットフォーム１３と接続されて、トリガー１８が極微弱光で照射されると触発し始め、触発パルス形態でシステム制御プラットフォーム１３に通知し、システム制御プラットフォーム１３はドライブ制御信号を送信して空間光変調器（SLM）４を作動させ、計数器１７が計数し始める。
【００３９】
本発明のまた１つの革新は、フィルタの作用は極微弱光の迷光をフィルタ除去し、後続のシステムに進入した極微弱光を検出に必要な波長範囲内にさせ、波長に対して要求がない場合は省略することができる。アッテネータは光を単一光子検出器の作動範囲に減衰し、測定された光強度が大きすぎ又は単一光子検出器１６のゲート制御時間が長すぎることによる検出器の飽和ことを防止することに用いられる。光源の光強度が極めて微弱であると、更にアッテネータを設置する必要がない。光アッテネータを光源の後に設置し、これにより光源は極微弱光光源になり、システムノイズがそれに従って増大し、従って、設置の位置を考慮すべきである。
【００４０】
空間光変調器（SLM）４は+12°と-12°で（ある型番の空間光変調器（SLM）４が+10°と-10°で反転することができる）反転することができ、本システムにおいて、+12°を受信可能な反射角度とし、-12°で反転すると最後の単一光子検出器１６に進入可能な反射光が非常に少なく、無視することができ、このため、光集束収集部材７の主光軸と光学イメージング部材の主光軸の角度が24°である。空間光変調器（SLM）４の作動領域を選択した後、システム制御プラットフォーム１３に該領域のランダム数ファイルをダウンロードして且つ生成し、空間光変調器（SLM）４のランダム反転を制御し、空間光変調器（SLM）４は反転する同時に計数器へ同期信号を送信し、これは空間光変調器（SLM）４と単一光子検出器１６の間の同期を保証しており、即ち空間光変調器（SLM）４が一回反転すると、単一光子検出器１６は該回の反転の時間間隔内の単一光子数を累積計数し、空間光変調器（SLM）４の反転を完成した後、計数器がクリアして累積計数を再び開始し、すべての計数はすべてデータ回線により圧縮センシング及び関連アルゴリズム１５に伝送されて、ランダムマトリックスとともに計算する。
【００４１】
本発明の更の革新は、前記の光学イメージング部材４と光集束収集部材７が、すべて光学レンズ群であり、それぞれ光学イメージングと光学集束を担当し、従来の技術の単一のレンズを用いてイメージングと集束することに限定されたものではないので、本発明の適用範囲が更に広くなり、空間光変調器（SLM）４においてイメージングの面積は最大に2048×1152画素までに達することができ、従来の技術と異なり、前記の光学イメージング部材４において、極微弱光が光学イメージング部材４を通過した後、空間光変調器（SLM）４に同じ大きさ又は縮小又は拡大の画像を形成し、且つ実際の需求に従ってイメージングを行うことができることである。後続の光路収集がファイバー結合技術に関連し、即ち集束した後の光束をファイバーに結合する。その革新は、ファイバー結合技術を用いて集束光を単一光子検出器１６上に収集することであり、結合の良否が直接にイメージング品質を影響し、このため、結合の調整も更に困難になる。
【００４２】
前記の時間分解単一光子計数二次元イメージングシステムは、物体に対してグレースケール画像しか形成しなく、単一光子検出器ラインアレイ２０で前記単一光子検出器１６を切り替える場合、説明する必要があるのは、赤、緑、青の三原色の下での光強度空間分布を取得すればカラー画像を形成することができ、且つ前記単一光子検出器ラインアレイ２０と光集束収集部材７の間の光路上に分光光度計を設置することにより、上記システムは時間分解極微弱光マルチスペクトルイメージングシステムに発展することができることを特徴とする。
【００４３】
図３に示すように、本発明は、時間分解極微弱光マルチスペクトルイメージングシステムを更に提供し、該システムは、異なる波長のレーザー又はハロゲンランプ１、トリガー１８、フィルタとアッテネータ２、光学イメージング部材３、DMD４、光集束収集部材７、分光光度計１９、単一光子検出器ラインアレイ２０、多チャネルカウンタ１０、システム制御プラットフォーム１３と圧縮センシング及び関連アルゴリズム１５からなる。単一光子検出器ラインアレイ２０は異なる波長に対する複数の単一光子計数器からなる。
【００４４】
前記の時間分解極微弱光マルチスペクトルイメージングシステムは、前記光集束収集部材７がスペクトルにおける異なる波長の光をそれぞれ相応の単一光子検出器上に収集集束すること用いられ、自由空間結合方式又はファイバー視準方式を採用して、単一光子点検出器に対して優先にファイバー視準方式を選択することができ、アバランシェダイオード配列に対して優先に自由空間結合方式を選択することができ、
前記分光光度計１９は、光視準部、分光部、角度測定部、光度観察及び測定部を備え、スペクトル分析と測定に用いられ、前記分光光度計はプリズム分光光度計又は格子分光光度計であり、フィルタホイールの方式を採用してスペクトル分光してもよいことを特徴とする。
【００４５】
上記２種のシステムにおいて、極微弱光が物体に当てた後の光発出方式は、主に物体が励起されて又は自己励起されて放出、反射、散乱、透過又は屈折することである。補足説明として、観測対象は自己発光の物体である場合、時間分解単一光子計数二次元イメージングシステムにおける極微弱光光源即ち異なる波長のレーザー又はハロゲンランプ１を除去してよく、図２、４に示すように、トリガー１８、フィルタ及びアッテネータ２を自発光物体の右側の光路上に平行移動し、直接に観測する。
【００４６】
上記2つのシステムは、すべて従来の技術の上に革新的に時間分解方法を追加したものであり、それぞれは二次元平面画像フレームの時間分解とスペクトル次元を追加した平面画像フレームの時間分解であり、いずれもフレームずつ測定する方式を採用するが、秒レベルの大規模時間の動的な物体に対してのみイメージングすることができ、時間分解精度が比較的に低く、ただし、多くの物理と化学過程、例えば分子のシス-トランス異性化と方向性のリラクゼーション、電荷とプロトンの移転、励起分子の衝突前期解離、エネルギー伝達と蛍光寿命及び水中の電子の溶媒和などは、ただ10-8秒で完成することができ、ピコ秒の時間分解精度の分析装置のみによりこれらの高速過程をタイムリーに観察することができる。
【００４７】
図５に示すように、本発明は上記高精度時間分解問題を解決するために、時間分解単一光子多次元イメージングシステムを更に提出し、該システムは同様に圧縮センシング原理に基づいて、レーザーパルスにより物体蛍光を励起し、光強度分布とイメージングスペクトルを取得し、且つ時間分解方策を用いてピコ秒レベルまでの時間分解精度を実現し、更に蛍光寿命などの一時的な周期性物質を検出することに用いられることができ、前記多次元は光強度次元、平面二次元、スペクトル次元及び時間次元という５つの次元を含む。前記システムは、
光源、イメージングスペクトル測定ユニット、電気検出ユニット、システム制御ユニット、処理ユニットを備え、
本発明の技術ロードマップは主に可視光と近赤外光の点及び配列単一光子検出技術、時間分解測定技術、空間光変調技術、圧縮センシング理論、単一光子計数イメージング技術、イメージングスペクトル技術、再構成アルゴリズム技術によりシステムの研究を完成する。前記システムは、光源I、イメージングスペクトル測定ユニットII、電気検出ユニットIII、システム制御ユニットIV、アルゴリズムユニットVという５つのユニットを含む。
【００４８】
時間分解単一光子多次元イメージングシステムの全体構造図は図５に示すように、異なる波長のレーザー又はハロゲンランプ１は光を発出し、時間シーケンスを有する観測対象に照射し、それにフィルタとアッテネータ２を経て、光学イメージング部材３に進入する。特別なフィルタフィルタを用いて迷光を除去し、且つ単一色フィルタを採用して選択的にある波長の光を通過させ、光源が比較的に強い際、多組のアッテネータを組み合わせて光減衰し、単一光子レベルに達し、それにより後続の点及び配列単一光子検出器９が飽和することを防止する。前記光学イメージング部材３は観測対象の画像を空間光変調器（SLM）４にイメージングし、後のものはランダム数発生器８が発生したランダム数を載せることにより光をランダム変調し、その射出光を一定のランダム的な確率で後続の光集束収集部材７の方向へ偏向し、圧縮センシング理論により、そのランダム性が高いほど、イメージングの品質が良好になる。集束收光の方式は多種があり、目的はすべて迷光を集束視準し、最後に点及び配列単一光子検出器９で受信されることができ、なお、図５に記載の集束收光光路はただ概略図であり、唯一の評価基準とするものではない。スペクトル検出を行うために、更に格子分光部材６を増設する必要があり、光信号に対してスペクトル分解測定を行い、点及び配列単一光子検出器９における各単一光子点検出器はそれぞれ各特定波長の光をキャプチャする。
【００４９】
上記技術案の補足として、スペクトル分光機構の異なりにより、一般的にフィルタリング型と分散型の分光に分けられる。フィルタリング型の分光に対して、予定された制限な部分のスペクトルのみに検出を実現することができ、且つフィルタホイールを回転するなどの方式によりスペクトル次元で「走査」を行う。分散型分光に対して、分散素子（プリズム又は格子）は空間光変調器（SLM）４が反射した異なる波長の光照射野を短い波長から長い波長まで順次に点及び配列単一光子検出器９の異なる位置に照射し、走査する必要がなく、各スペクトル部分の画像を同時に取得する。該分光方式において、スペクトル解像度の高さが分散素子（プリズム又は格子）に到達する入射光の視準度と比例にし、視準性が良いほどスペクトル解像度が高くなり、且つスペクトル解像度と空間解像度は互いに制限される。ここで、説明のための例示として、ブレーズ格子を選択してスペクトル分散分光する。凹面鏡５の視準と光束拡張により空間光変調器（SLM）４の反射光照射野を視準してブレーズ格子面に投影し、光束をできるだけ全体の格子面に完全に照射する。更に光集束収集部材７における集束システムにより、準平行光は格子に入射し反射した後、集束され、異なる波長の光をその焦平面で空間的分離を実現する。焦平面に点及び配列単一光子検出器９を設置すれば、各スペクトル情報を記録することができ、圧縮センシング及び関連アルゴリズムを用いて全体のスペクトル範囲内の画像情報を取得することができる。
【００５０】
電気検出ユニットIIIにおいて、点及び配列単一光子検出器９は入力した光信号を有効的なパルス信号に処理して出力し、多チャネルカウンタ１０により特定の時間間隔内の単一光子数を記録し、同時に高精度時間測定器１１を並列接続することにより相応のタイムスタンプ（タイムコード）を記録する。最後に光子計数、タイムスタンプ（タイムコード）及びランダム数発生器８が出力したランダムマトリックスはともにパッケージ化されてアルゴリズムユニットVのデータパケット記憶器１４に進入する。
【００５１】
電気検出ユニットIVにおいて、システム制御プラットフォーム１３は各ユニット部材のイネーブルを制御し、且つ圧縮センシング原理に基づく時間分解測定において重大な役割を果たす。具体的に、それは一定の時間帯に空間光変調器（SLM）４の相応のランダムマトリックスの不変を維持し、且つ精確的にそれを制御してある特定な時刻に瞬間に次のフレームランダムマトリックスに反転することができる。各種の時間分解測定方法により、それは光源I、点及び配列単一光子検出器９、多チャネルカウンタ１０及び高精度時間測定器１１の間のペース協調を制御し、且つ非同期の時間差を除去する。なお、遅延測定に基づく時間分解方法においてデジタル遅延器１２を増設することにより、点及び配列単一光子検出器９に対するピコ秒レベルのゲート制御を完成する。光子到達時間に基づく時間分解方法において、多チャネルカウンタ１０における時間幅変換器は光子を取得する時間（start to stop）を電圧形態で記録し、対応の時間チャネルに記録し、且つ光子が到達する時間に従って光子数を区分することにより、１つの周期内での各時間間隔内の複数の累積計数を統計する。
【００５２】
システムは圧縮センシング原理により多回の測定を繰り返すと、シリーズの計数値、タイムスタンプ（タイムコード）及びランダムマトリックスを取得することができ、アルゴリズムユニットVにおける圧縮センシング及び関連アルゴリズム１５は最後にこれらの入力により必要するイメージングを再構成し、主に、単一光子二次元イメージングI（x，y）、単一光子三次元イメージングI（x，y，z）、単一光子二次元イメージングスペクトルI（x，y，λ）、単一光子時間分解二次元イメージング（例えば蛍光寿命）I（x，y，t）、単一光子三次元イメージングスペクトルI（x，y，z，λ）、単一光子時間分解三次元イメージングI（x，y，z，t）、単一光子時間分解二次元イメージングスペクトルI（x，y，λ，t）及び単一光子時間分解三次元イメージングスペクトルI（x，y，z，λ，t）がある。
【００５３】
なお、デジタルマイクロミラーデバイス（DMD）を１つの空間光変調器（SLM）４の特例として、その単一のマイクロミラーの反射機構原理を説明する。図６に示すように、入射光線がマイクロミラー配列の単一のマイクロミラー法線と24°になる際、反射光線も法線と24°になり、マイクロミラーが+12°で反転する際、マイクロミラーの法線も時計回りで+12°反転し、反射定律により、反射光線は時計回りで+24°反転する必要があり、即ち初期位置際の法線と同一の直線になり、該初始位置際の法線方向を後続の光集束収集部材７の受信方向に設置することができ、即ち射出光が最後に点及び配列単一光子検出器９に進入する。同様に、マイクロミラーが-12°反転する際、この際の反射光線は初始位置際の法線と-48°になり、後続の光集束収集部材にほぼ進入しなく、このため、マイクロミラーが-12°反転する際の反射光を無視してよい。ここには、時計回りの反転を正とし、反時計回りの反転を負とする。受信の方向が逆転して設置されてもよい。
【００５４】
上記技術案の補足として、図７に示す正負画像構造の簡素化原理図は、空間光変調器（SLM）４の正負級光路上にいずれも点及び配列単一光子検出器９を設置してよく、同期に符合測定をし、正負級ランダム測定マトリックスの相互関連性を利用して、関連回復アルゴリズムを採用し、再構成画像の品質を向上させるとともに回復アルゴリズムの実行時間を短縮することができることを特徴とする。
【００５５】
前記の方策一の時間分解測定方法は、時間間隔測定に基づいた時間分解方法であり、その原理ブロック図は図８に示すようなものであり、図９は時間間隔を分けて単一光子数を記録する模式図であり、前記の方策二の時間分解測定方法は遅延測定に基づいた時間分解方法であり、その原理ブロック図は図１０に示すようなものであり、図１１は統計して差を取得した後時間間隔内の単一光子計数の模式図であり、最後に時間間隔測定に基づく時間分解問題に変換し、同様に、前記の方策三の時間分解測定方法は光子到達時間に基づいた時間分解方法であり、最後に図９に示すような１つの周期内に各時間間隔内の複数の累積計数を統計するようになる。具体的な方法は発明内容において詳しく説明したので、ここで再び説明しない。
【００５６】
上記技術案の改善として、システムの信号対雑音比（signal to noise ratio，SNRと略称する）を分析するために、SNRを信号とシステムノイズの分散の比とする。システムノイズは、点及び配列単一光子検出器９の暗計数、環境ノイズ、電気ノイズ、光学ノイズなどを含む。分散は信号の変動状況と解することができ、システムノイズの変動が信号の変動をカバーすると、圧縮センシングアルゴリズムは働かなくなり、システムノイズの変動が信号の変動よりも小さくし又ははるかに小さいと、ほぼ完全に画像を再構成することができる。システムの信号対雑音比を向上させるために、システムに対してシールパッケージしてよく、且つ点及び配列単一光子検出器９の相応のパラメータを向上させ、システムの安定性を向上させ、これはイメージング品質の向上に大幅に寄与する。
【００５７】
以上のように、本発明が提出した時間分解単一光子又は極微弱光多次元イメージングスペクトルシステム及び方法は、新世代高性能時間分解単一光子多次元イメージングシステムに属し、それは圧縮感知（Compressive Sensing，CSと略称する）理論と単一光子検出技術に基づいて実現するものであり、次元の増加による次元問題と感度問題を効果的に解決することができる。圧縮センシング原理はシャノンサンプリング定理のリミットを超えて、ランダムサンプリングの方式で、より少ないデータサンプリング数（ナイキストサンプリング定理のリミットよりもはるかに少ない）で、完全にオリジナル信号を回復し、且つ高ロバスト性を有する。該原理は主に、圧縮サンプリングステップ、スパース変換ステップ及びアルゴリズム再構成ステップを含む。圧縮サンプリングは、測定された信号が高次元から低次元へマッピングして収集された過程であり、前記のスパース変換は適切なΨを選択して、ｘがΨ変化を経て取得した値ｘ’はスパースであり、即ちｘがΨフレームでスパース表現することができ、前記のアルゴリズム再構成は観測データｙ、測定マトリックスＡ及びフレームΨが既知の条件下でｙ＝ＡΨｘ’＋ｅを解く過程であり、最後に
【数８】
0006002232.tif 000009
でｘを反転して取得する。
【００５８】
本発明において、時間分解単一光子多次元イメージングシステムは強大な機能を実現することができ、国際での唯一の単一光子検出感度、ピコ秒時間解像度、ナノ空間解像度、広いスペクトルバンドを有し、点検出器でイメージングを実現する最先進の分析システムであり、速度が速く、解像度が高く、安定性が高く、効率が高く、試験成功の確率が大きいなどの利点を有し、材料学、物理、化学、生命学などの関連基礎先進分野の研究を行うことに必要する定量分析工具であり、且つ有効的に配列式ナノ材料構造と生体高分子及び細胞の微細構造を試験することができ、工業においてナノ半導体機器、生物製薬、臨床腫瘍組織及び産品のバッチ検出と品質監視などの科研分野にも適用する。
【００５９】
本発明は、同時に以前の研究の改善と革新に基づいたものであり、主に該分野の『単一光子計数イメージングシステム及びその方法』（中国出願番号又は特許番号：201110103559.3，出願者又は特許権者：中国科学院空間科学と応用研究センター）、『極微弱光マルチスペクトルイメージング方法及びそのシステム』という二つの特許に基づいたものである。本発明は、そのうえに、上記特定の場合に時間分解次元情報を用いて空間第三次元距離パラメータを導出し、８種類のイメージングに拡張する：単一光子二次元イメージングI（x，y）、単一光子三次元イメージングI（x，y，z）、単一光子二次元イメージングスペクトルI（x，y，λ）、単一光子時間分解二次元イメージング（例えば蛍光寿命）I（x，y，t）、単一光子三次元イメージングスペクトルI（x，y，z，λ）、単一光子時間分解三次元イメージングI（x，y，z，t）、単一光子時間分解二次元イメージングスペクトルI（x，y，λ，t）及び単一光子時間分解三次元イメージングスペクトルI（x，y，z，λ，t）である。本発明は初めて圧縮センシング理論、単一光子検出技術及び蛍光寿命測定技術を一体に組み合わせて、一言で言えば、本発明は時間分解単一光子又は極微弱光多次元イメージングスペクトルシステム及び方法を提供し、単一光子時間分解イメージングスペクトルの技術分野に属し、前記システムは可視光と近赤外光の点及び配列単一光子検出技術、時間分解測定技術、空間光変調技術、圧縮センシング理論、単一光子計数イメージング技術、イメージングスペクトル技術、再構成アルゴリズム技術を合わせて、極微弱光単一光子レベルで物体を反射、散乱又は透過する空間二次元、光強度、時間分解、スペクトル分解の５次元パラメータ情報を有するイメージングスペクトル測定を実現する。該システムは主に光源、イメージングスペクトル測定ユニット、電気検出ユニット、システム制御ユニット、アルゴリズムユニットからなり、物体情報を載せる光は空間光変調器にイメージングされて、圧縮センシング理論に従ってランダム変調されて、且つ点及び配列単一光子検出器を用いてその反射光を収集し、同時に光子数と光子到達時間を記録し、圧縮センシングアルゴリズムとイメージングスペクトル関連アルゴリズムを用いて再構成し、単一光子検出感度、高時間解像度、広いスペクトル範囲を有し、単一分子生物物理学、材料欠陥の検出、ナノ材料、マイクロエレクトロニクス、量子ドット、生命科学及び新エネルギー光電変換材料などの多くの新たなハイテク産業分野に広く応用することができる。
【００６０】
本発明は、単一光子検出器を用いて光子検出し、計数器を用いて該検出器が出力する離散パルスの数を記録する。時間分解技術を用いる際、タイムスケール変換（Time to Amplitude Converter，TACと略称する）技術を用いて時間分解測定方法を発展し、更に圧縮センシング理論との組み合わせに適しており、圧縮センシング理論を用いる際、核心アルゴリズムは共用であるが、該発明に適するようにアルゴリズムを発展させ、８種類のイメージングに対してそれぞれ８種類の相応の改善措置を提出し、具体的に発明内容を参考することができ、空間光変調器を用いる際は、主にランダム数発生器によりベルヌーイバイナリランダムマトリックスをロードして、光に対するランダム光変調を実現し、スペクトル分光を用いる際は公知の分光技術を用いる。
【００６１】
最後に説明する必要があるのは、具体的な実施形態において、前記アルゴリズムは大量の試験データの検証を経て、真実で信頼可能であり、ハードウェアを合わせて本発明の技術案を実現することができる。前記のすべての実施例は更に本発明を詳しく説明するものであり、絶対なものではなくて、相応的に発展することができる。当業者は、本発明の技術案に対する修正、添加、削除又は同等交換は、すべて本発明の技術案の主旨と範囲を逸脱しないものであり、いずれも本発明の請求の範囲に属すべきであることを理解すべきである。
【符号の説明】
【００６２】
１  異なる波長のレーザー又はハロゲンランプ
２  フィルタとアッテネータ
３  光学イメージング部材
４  空間光変調器（SLM）又はDMD
５  凹面鏡（選択可能）
６  格子分光部材
７  光集束収集部材
８  ランダム数発生器
９  点及び配列単一光子検出器
１０  多チャネルカウンタ
１１  高精度時間測定器
１２  デジタル遅延器（選択可能）
１３  システム制御プラットフォーム
１４  データパケット記憶器
１５  圧縮センシング及び関連アルゴリズム
１６  単一光子検出器
１７  計数器
１８  トリガー
１９  分光光度計
２０  単一光子検出器ラインアレイ
【図１】
0006002232.tif 000010
【図２】
0006002232.tif 000011
【図３】
0006002232.tif 000012
【図４】
0006002232.tif 000013
【図５】
0006002232.tif 000014
【図６】
0006002232.tif 000015
【図７】
0006002232.tif 000016
【図８】
0006002232.tif 000017
【図９】
0006002232.tif 000018
【図１０】
0006002232.tif 000019
【図１１】
0006002232.tif 000020