遺伝子配列解析装置および遺伝子配列解析方法
【課題】剛性が不明の作業台に設置された場合でも、簡易な構成で、性能を維持することができる遺伝子配列解析装置および遺伝子配列解析方法を提供する。
【解決手段】筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力のうち少なくとも一方の上限値を計算する演算部と、前記分注ロボットに要する駆動力が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または前記ステージ移動に要する駆動力が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、を備える。
【解決手段】筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力のうち少なくとも一方の上限値を計算する演算部と、前記分注ロボットに要する駆動力が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または前記ステージ移動に要する駆動力が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、を備える。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、遺伝子配列解析装置および遺伝子配列解析方法に関し、例えば、筐体内に、試薬を分注する分注ロボットと、前記試薬を載置し、二次元に移動するステージと、前記試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置に関する。
【背景技術】
【０００２】
遺伝子配列解析装置（以下、「ＤＮＡシーケンサー」という）は、ＤＮＡの塩基配列を決定する装置である。ＤＮＡシーケンサーでは、１ミクロンのビーズ上にＤＮＡを固定し、エマルジョンＰＣＲによりＤＮＡを増幅させ、これらのビーズをフローセルのセルに固定し、セル内に試薬を送液し、ＤＮＡの伸長あるいはライゲーション反応によって蛍光基質をＤＮＡ鎖に取り込ませる。この操作により、フローセルに多数のビーズが固定されるため、励起光でフローセルの撮影領域を照射すると、ビーズから蛍光が発せられる。全ビーズの蛍光情報は、ＸＹステージを逐次移動させ、ＣＣＤカメラでフローセルの撮影領域を撮影することによって得られる。このようにして得られた蛍光情報に基づき、ＤＮＡの塩基配列を決定する。
【０００３】
ＤＮＡシーケンサーで行う塩基配列の決定を低コスト化するには、塩基配列を決定するまでの時間を短縮し、スループットを増加させることが求められる。また、スループットを増加させる一方、装置の性能は維持しなければならない。装置の性能を低下させる要因として、装置の振動が挙げられる。装置が振動すると、ＣＣＤカメラの位置とフローセルの位置との間に相対的な位置ずれが発生する。この状態のままフローセルを撮影すれば、ＣＣＤカメラの視野に散在する蛍光の位置がずれて撮影されるので、正確な蛍光情報が得られなくなる。また、試料をフローセルに注入する作業に支障が生じる。したがって、装置の性能を維持するには、装置の振動を抑制する必要がある。
【０００４】
特許文献１には、装置の構造解析と特性評価を実施することで、所望の性能に達しない場合は設計データを修正し、短期間に所望の性能を実現するリソグラフィ装置が開示されている。
【０００５】
特許文献２には、装置が設置される床を加振させることにより床構造の応答情報を測り、床の剛性が不足している場合には床構造を補強することで、所望の性能を実現するリソグラフィの製造方法が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第３９２８２００号公報
【特許文献２】特開２００５−２４３８８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
蛍光の撮影中、ＤＮＡシーケンサーには、様々な要因で振動が発生する。例えば、装置に人がぶつかって発生する振動、オペレータが情報を入力することで発生する振動、装置内のステージ移動や分注動作によって発生する振動である。
【０００８】
また、ＤＮＡシーケンサーは、一般的な研究室の作業台等に設置されることが多い。ＤＮＡシーケンサーが一般的な作業台に設置される場合、設置される作業台の剛性は不明である。例えば、作業台の剛性が大きいほど振動は起きにくいが、作業台の剛性が小さいほど、装置の振動は大きくなり、性能が低下する傾向がある。
【０００９】
したがって、設置される作業台の剛性を問わず、装置の性能を維持しなければならない。
【００１０】
上記特許文献１および２は、いずれも作業台の剛性を考慮していない。よって、剛性の小さい作業台に設置された場合は、振動により性能が低下するという問題がある。
【００１１】
本発明は上記課題に鑑みてなされたものであり、剛性が不明の作業台に設置された場合でも、簡易な構成で、性能を維持することができる遺伝子配列解析装置および遺伝子配列解析方法を提供する。
【課題を解決するための手段】
【００１２】
上記課題を解決するために本発明は、画像に生じる蛍光の位置ずれの許容値から機構の駆動力の上限値を求め、この上限値を上回らないように振動源となる機構を制御する。
【００１３】
また、蛍光の位置ずれの許容値から筐体の加速度の上限値を求め、筐体に設置した加速度センサがこの上限値を上回るときは、所定の整定時間を待つか、振動源となる機構を制御する。
【００１４】
すなわち、本発明は、筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力のうち少なくとも一方の上限値を計算する演算部と、前記分注ロボットに要する駆動力が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または前記ステージ移動に要する駆動力が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、を備えることを特徴とする。
【００１５】
また、本発明は、筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算する演算部と、前記筐体の加速度を検出し、前記検出した加速度が前記上限値を上回る場合、所定の整定時間の経過を待つ加速度検出部と、を備えることを特徴とする。
【００１６】
また、本発明は、筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算する演算部と、前記筐体の加速度を検出する加速度検出部と、前記加速度が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または前記加速度が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、を備えることを特徴とする。
【発明の効果】
【００１７】
剛性が不明の作業台に設置された場合でも、簡易な構成で、性能を維持することができる遺伝子配列解析装置および遺伝子配列解析方法を提供する。
【図面の簡単な説明】
【００１８】
【図１】図１は、一般的なＤＮＡシーケンサーの概略構成図である。
【図２】図２は、一般的な作業台２０１に設置された本発明の実施形態に係るＤＮＡシーケンサーである。
【図３】図３は、本発明の実施形態に係るＤＮＡシーケンサーのシミュレーションモデルである。
【図４】図４は、本発明の実施形態に係る固有角速度ωを求める概念図である。
【図５】図５は、分注ロボットの駆動力を制御することで画像に生じる蛍光の位置ずれを抑制する本発明の実施形態に係るフローチャートである。
【図６】図６は、筐体の加速度を検出することで画像に生じる蛍光の位置ずれを抑制する本発明の実施形態に係るフローチャートである。
【図７】図７は、筐体の加速度を検出することで画像に生じる蛍光の位置ずれを抑制する本発明の実施形態に係るフローチャートである。
【発明を実施するための形態】
【００１９】
本発明は、筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置および遺伝子配列解析方法に関するものである。以下、添付図面を参照して本発明の実施形態について説明する。ただし、本実施形態は、本発明を実施するための一例に過ぎず、本発明の技術的範囲を限定するものではないことに注意すべきである。また、各図において共通の構成には同一の参照番号が付されている。
【００２０】
＜ＤＮＡシーケンサー＞
図１は、一般的なＤＮＡシーケンサーの概略構成図である。ＤＮＡシーケンサーは、底部に防振ゴム１０１と側面に排熱ファン１０２を有する筐体１０３内に、除振マウント１０４上に設置された光学検出部１０５、試薬ラック１０６に並べられた試薬１０７を注入し、フローセル１０８に分注する分注ロボット１０９、使用済みの試薬を回収する廃液タンク１１０、底部に防振ゴム１０１を有し、試薬を送液するポンプ１１１を備える。光学検出部１０５は、フローセル１０８を載置するステージ１１２、ステージ１１２上の試薬を照射する光源１１３、試薬から発する蛍光を撮影するＣＣＤカメラ等の画像取得装置１１４を備える。分注ロボット１０９およびステージ１１２は、分注ロボット移動装置およびステージ移動装置によって少なくとも水平方向に移動可能である。防振ゴム１０１は、常時駆動しているポンプ１１１、排熱ファン１０２からの振動が筐体１０３に接する他の機構に伝達されないように、振動を減衰させるものである。また、マウント１０４は、筐体からの振動が光学検出部に伝達されないように、振動を減衰させるものである。
【００２１】
このＤＮＡシーケンサーは、一般的な作業台に設置され、そこでＤＮＡ検出が行われる。このような使用環境で、特に問題になるのは、分注ロボット１０９の分注作業をしながら、または高速でステージ移動をしながら、フローセル上のビーズから発せられる蛍光を撮影するときである。作業台の剛性が低ければ、分注作業やステージ移動によって筐体に発生する振動が大きくなり、この振動によってＣＣＤカメラの視野で蛍光の位置がずれ、正確な蛍光情報が得られなくなるためである。蛍光の位置ずれは、画像取得装置１１４の分解能に対して、充分に小さな値にする必要がある。例えば、画像取得装置１１４の分解能が０．３８μｍであれば、蛍光の位置ずれの許容値は０．１μｍ以下となる。
【００２２】
図２は、一般的な作業台２０１に設置された本発明の実施形態に係るＤＮＡシーケンサーである。作業台２０１は、ゴム製のキャスタ２０２を有しており、キャスタ２０２を介して床に設置される。本発明のＤＮＡシーケンサーは、ＤＮＡシーケンサを制御する制御手段２０３、制御手段２０３に接続された表示装置２０４、ポインティングデバイス２０５、および記憶装置２０６、筐体に設置された加速度センサ２０７を設けている。その他の構成は、特記しない限り図１と同じである。
【００２３】
制御手段２０３は、分注ロボット移動装置の駆動力を制御する分注ロボット制御部２０８、ステージ移動装置を制御するステージ制御部２０９、画像取得装置１１４を制御する画像制御部２１０、加速度センサ２０７からの加速度を検出する加速度検出部２１１、所定の関係式を算出し、関係式に基づき計算を行う演算部２１２を備える。
【００２４】
記憶装置２０６には、ＤＮＡシーケンサーによる検査中の操作履歴およびエラー履歴を記録したログファイル２１３、演算部２１２が算出した関係式２１４、および関係式に代入する測定値２１５が格納されている。
【００２５】
＜分注ロボットによる振動の影響＞
ここでは、分注ロボットによって発生した振動が、蛍光の位置ずれをもたらす原理について説明する。
【００２６】
この原理は、振動する機構とそれに接する機構との間には、所定のバネ定数および減衰定数があるとみなせるという概念に基づく。分注ロボット１０９は筐体１０３に設置されているため、分注ロボット１０９が駆動力を発生すれば、筐体１０３の振動は増幅し、所定の固有角周波数ωで振動する。また、分注ロボット１０９が駆動力を停止すれば筐体１０３の振動は減衰する。したがって、分注ロボット１０９と筐体１０３の間には、所定のバネ定数および減衰定数があるとみなせる。同様の理由で、筐体１０３と作業台２０１の間、筐体１０３と光学検出部１０５の間、光学検出部１０５とステージ１１２の間には、それぞれ所定のバネ定数および減衰定数があるとみなせる。
【００２７】
図３は、本発明の実施形態に係るＤＮＡシーケンサーのシミュレーションモデルである。図示するように、各機構と機構のバネ定数および減衰定数については、床と作業台２０１の間をｋ０、ｃ０、作業台２０１と筐体１０３の間をｋ１およびｃ１、筐体１０３と光学検出器１０５の間をｋ２およびｃ２、光学検出器１０５とステージ１１２の間をｋ３およびｃ３、筐体１０３と分注ロボット１０９の間をｋ４およびｃ４とする。各機構の質量については、作業台２０１をｍ０、筐体１０３をｍ１、光学検出器１０５をｍ２、ステージ１１２をｍ３、分注ロボット１０９をｍ４とする。各機構の中心座標からの変位については、作業台２０１をｘ０、筐体１０３をｘ１、光学検出器１０５をｘ２、ステージ１１２をｘ３、分注ロボット１０９をｘ４とする。各機構が発生する駆動力については、作業台２０１をｈ０、筐体１０３をｈ１、光学検出器１０５をｈ２、ステージ１１２をｈ３、分注ロボット１０９をｈ４とする。
【００２８】
以上の設定に基づき、各機構ごとに運動方程式ｆ０〜ｆ４を立てると、次式（１）が成立する。
【００２９】
【数１】
【００３０】
この運動方程式を微分方程式によって解くため、駆動力ｈ（ｔ）を次式（２）のように置く。
【数２】
【００３１】
式（２）より、式（３）が導ける。
【数３】
【００３２】
ここで、変位Ｘおよび変位の位相差φは次式から求められる。
【数４】
【００３３】
式（２）および（３）を式（１）に代入し、運動方程式を次式（５）のように書き換える。
【００３４】
【数５】
【００３５】
式（５）を行列に書き換えると次式（６）になる。
【数６】
【００３６】
ここで、Ａ〜Ｍは、次式（７）のように置いた。
【数７】
【００３７】
式（６）の逆行列をＮ−１とおくと、次式（８）になる。
【数８】
【００３８】
したがって、各機構の変位Ｘ０〜Ｘ４は、Ｎ−１の因子となるバネ定数ｋ０〜ｋ４、減衰定数ｃ０〜ｃ４、および各機構で共通の固有角速度ω、ならびに各機構の駆動力Ｈ０〜Ｈ４を求めることで、計算可能である。なお、Ｘ０〜Ｘ４は複素数で表記される。
【００３９】
次に、バネ定数ｋ、減衰定数ｃ、および固有角速度ωの求め方について説明する。
図４は、本発明の実施形態に係る固有角速度ωを求める概念図である。固有角速度ωは、筐体１０３に設置された加速度センサ２０７によって、周波数ごとに測定する。分注ロボット１０９を一定の駆動力で動作させ、周波数を変えながら（たとえば１〜２０Ｈｚ）、分注ロボット１０９を水平方向に正弦波駆動（スイープ加振）をさせる。周波数ごとの加速度測定結果から、加速度振幅が最大となる固有周波数ｆを求め、固有角速度ωを計算する。図４の場合は、固有周波数ｆは１０Ｈｚ、固有角速度ωは６２．８ｒａｄ／ｓである。
【００４０】
バネ定数ｋは、次式（９）により計算可能である。
【数９】
【００４１】
次に、減衰係数ｃを求める。まず、加速度の測定結果から、固有周波数における共振倍率Ｑを求める。
【００４２】
【数１０】
【００４３】
減衰比ζは、次式（１１）から求められる。
【数１１】
【００４４】
また、作業台のキャスター（またはアジャスター）の設置面はゴムであることが多いので、減衰比ζをゴムの損失係数（η＝０．１）から求めても良い。この場合、減衰比は次式（１１）のようになる。
【００４５】
【数１２】
【００４６】
減衰係数ｃは、減衰比ζを求めた後、次式により計算可能である。
【数１３】
【００４７】
以上の式から求めた各パラメータの値の一例を挙げておく。
作業台の固有周波数 ｆ＝１０Ｈｚ
作業台の固有角速度 ω＝ ２πｆ＝ ６２．８ （ｒａｄ／ｓ）
作業台の質量 ｍ１＝１００（ｋｇ）
作業台のばね定数 ｋ１＝４００（ｋＮ／ｍ）
作業台の共振倍率 Ｑ＝０．１／０．０１＝１０
作業台の減衰比 ζ＝０．０５
作業台の減衰係数 ｃ１＝１９００（Ｎ／（ｍ／ｓ））
【００４８】
＜第１の実施形態＞
画像の撮影中でも、分注ロボットは分注作業を行っている。分注ロボットからの振動は、筐体１０３、光学検出部１０５（画像取得装置１１４）、ステージ１１２の順に伝わる。したがって、画像取得装置１１４とステージ１１２との間に相対的な位置ずれが発生し、これによって画像取得装置１１４が撮影する時に蛍光の位置がずれる。
【００４９】
ここでは、式（７）を用いて、分注ロボット１０９の駆動力ｈ４（ｔ）と蛍光の位置ずれの関係式を求める。この関係式に位置ずれの許容値を代入し、分注ロボット１０９の駆動力の上限値ｈ４ｍａｘを計算する。よって、撮影中、分注ロボット１０９の駆動力が上限値ｈ４ｍａｘ以上とならないように、振動源となる分注ロボット１０９の駆動力を制御すれば、位置ずれは許容値を満足する。以下、分注ロボットの駆動力ｈ４（ｔ）と位置ずれの関係式（１５）を導出する。
【００５０】
蛍光の位置ずれは、ステージの変位と光学検出器の変位との差（ｘ３−ｘ２）に等しいとみなせる。また、分注ロボット以外の機構は駆動力を発生しないとし、Ｈ０〜Ｈ３＝０とする。これと、式（２）および（８）から、次式（１４）が求まる。ここで、Ｂ（ω）は、ｋおよびｃの関数とする。
【００５１】
【数１４】
【００５２】
また、式（３）より、次式（１５）が求まる。
【数１５】
【００５３】
式（１４）および（１５）より、次式（１６）が導ける。
【数１６】
【００５４】
ここで、Ｘ２〜Ｘ４およびＢ（ω）は、式（８）またはｋ、ｃ、ωにより計算可能である。したがって、式（１６）に、ステージの変位と光学検出器の変位との差（ｘ３−ｘ２）の許容値（例えば０．１μｍ）を代入すれば、分注ロボットの駆動力ｈ４（ｔ）の上限値が計算可能である。
【００５５】
図５は、分注ロボット１０９の駆動力を制御することで画像ずれを抑制するフローチャートである。
【００５６】
Ｓ５００では、演算部２１２が、式（１６）等を用いて、分注ロボット１０９の駆動力ｈ４（ｔ）の上限値ｈ４ｍａｘを計算する。ステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）は、ユーザが入力するか、演算部２１２が画像の分解能から計算する。また、固有角速度ω、バネ定数ｋ、減衰定数ｃは、予め値が分かっている場合はその値を式（１６）に代入するが、値が分からなければ次のように求める。まず、分注ロボット制御部２０８が、分注ロボット１０９を一定の駆動力で動作させ、周波数を変えながら、分注ロボット１０９を水平方向に正弦波駆動（スイープ加振）させる。次に、加速度検出部２１１が、筐体に設けた加速度センサ２０７により、周波数ごとの加速度測定結果を求める。そして、演算部２１２が、加速度振幅が最大となる固有周波数ｆを求め、固有角速度ω、バネ定数ｋ、減衰定数ｃを計算する。
【００５７】
Ｓ５０１およびＳ５０２では、分注ロボット制御部２０８が、分注ロボット１０９の駆動を開始し、駆動力を制御する。分注ロボット１０９は、例えばタイミングベルトを介してパルスモータにより制御される。
【００５８】
Ｓ５０３では、分注ロボット制御部２０８が、予め計算した上限値ｈ４ｍａｘと駆動力とを比較する。駆動力が上限値ｈ４ｍａｘより小さければ、Ｓ５０４に進む。駆動力が上限値ｈ４ｍａｘ以上であれば、Ｓ５０５に進む。
【００５９】
Ｓ５０４では、画像制御部２１０が、画像取得装置１１４による撮影をする。
Ｓ５０５では、分注ロボット制御部２０８が、駆動力超過のエラーをログファイル２１３に記録する。その後、Ｓ５０２に戻り、分注ロボット制御部２０８が、駆動力が小さくなるように分注ロボット１０９を制御する。
【００６０】
以上の通り、本実施形態では、式（１６）にステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）、固有角速度ω、バネ定数ｋ、減衰定数ｃを代入し、分注ロボット１０９の駆動力の上限値ｈ４ｍａｘを計算する。そして、分注ロボット１０９の駆動力が上限値ｈ４ｍａｘ以上とならないように、振動源となる分注ロボット１０９の駆動力を制御する。また、撮影中に分注ロボット１０９の駆動力が上限値を超える場合には、駆動力超過のエラーをログファイル２１３に記録する。
【００６１】
これにより、分注ロボットによる分注作業を行いながら画像を撮影しても、画像ずれを許容値以内に抑えることが可能となる。また、ログファイルを参照すれば、撮影した画像と、駆動力超過との因果関係を確認することができる。
【００６２】
したがって、剛性が不明の作業台に設置されたＤＮＡシーケンサーの性能を維持することができる。また、画像の取得に失敗する回数が低減されるので、高価な試薬と検査時間を費やさなくてすむ。
【００６３】
なお、ここでは、分注ロボットの駆動力を制御する方法について述べたが、式（７）に基づき、その他の振動源となる機構の駆動力を制御しても良い。例えば、ステージ移動装置の駆動力ｈ３（ｔ）も制御可能である。
【００６４】
＜第２の実施形態＞
画像の撮影中でも、装置に人がぶつかったり、オペレータが情報を入力したりすることによって装置に振動が発生している。装置に発生した振動は、筐体１０３を介して、光学検出部１０５（画像取得装置１１４）およびステージ１１２に伝わる。したがって、画像取得装置１１４とステージ１１２との間に相対的な位置ずれが発生し、これによって画像取得装置１１４が撮影する画像にずれが発生する。
【００６５】
ここでは、式（７）を用いて、筐体の加速度ａ１（ｔ）と画像ずれの関係式を求める。この関係式に画像ずれの許容値を代入し、筐体１０３の加速度の上限値ａ１ｍａｘを計算する。よって、撮影中、筐体１０３に設置した加速度センサ２０７が上限値ａ１ｍａｘ以上を検出した場合は、加速度が上限値ａ１ｍａｘを下回るまで画像取得装置１１４の撮影を停止する。これにより、画像ずれは許容値を満足する。以下、筐体１０３の加速度ａ１（ｔ）と画像ずれの関係式（１９）を導出する。
【００６６】
画像ずれは、ステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）に等しいとみなせる。また、式（３）から、次式（１７）および（１８）が求まる。
【００６７】
【数１７】
【００６８】
【数１８】
【００６９】
式（１７）および（１８）から、次式（１９）が求まる。
【数１９】
【００７０】
ここで、ωは測定から求め、Ｘ１〜Ｘ３は式（７）により計算可能である。したがって、式（１９）に、ステージの変位と光学検出器１０５の変位との差（ｘ３−ｘ２）の許容値（例えば０．１μｍ）を代入すれば、筐体１０３の加速度ａ１（ｔ）の上限値が計算可能である。
【００７１】
図６は、筐体１０３の加速度を検出することで画像ずれを抑制するフローチャートである。
【００７２】
Ｓ６００では、演算部２１２が、式（１９）等を用いて、筐体１０３の加速度ａ１（ｔ）の上限値をａ１ｍａｘ計算する。なお、ステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）は、ユーザが入力してもよいし、演算部２１２が画像の分解能から計算しても良い。固有角速度ω、バネ定数ｋ、減衰定数ｃは、前述のＳ５００で説明したとおりである。
【００７３】
Ｓ６０１および６０２では、加速度検出部２１１が、筐体１０３の加速度を検出し、予め計算した上限値ａ１ｍａｘと検出した加速度を比較する。検出した加速度が上限値より小さければ、Ｓ６０３に進む。検出した加速度が上限値以上になれば、Ｓ６０４に進む。
【００７４】
Ｓ６０３では、画像制御部２１０が、画像取得装置１１４による撮影をする。
Ｓ６０４では、加速度検出部２１１が、加速度超過のエラーをログファイル２１３に記録する。
【００７５】
Ｓ６０５では、加速度検出部２１１が、所定の整定時間の経過を待った後、Ｓ６０１に戻り、再度筐体の加速度を検出する。ここで整定時間とは、例えば、筐体１０３の固有周波数ｆから求めた周期の定数倍に設定してもよいし、タイマー等で任意の時間に設定してもよい。なお、所定の整定時間の経過を待たずに、検査を中止してもよい。
【００７６】
以上の通り、本実施形態では、式（１９）にステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）、固有角速度ω、バネ定数ｋ、減衰定数ｃを代入し、筐体１０３の加速度の上限値ａ１ｍａｘを計算する。そして、筐体に設置した加速度センサが上限値ａ１ｍａｘ以上を検出した場合は、加速度が上限値ａ１ｍａｘを下回るまで画像取得装置１１４の撮影を停止する。また、撮影中に筐体１０３の加速度が上限値を超える場合には、加速度超過のエラーをログファイル２１３に記録する。
【００７７】
これにより、装置に人がぶつかったり、オペレータが情報を入力したりすることによって装置に振動が発生しても、画像ずれを許容値以内に抑えることが可能となる。また、ログファイルを参照すれば、撮影した画像と、加速度超過との因果関係を確認することができる。
【００７８】
したがって、剛性が不明の作業台に設置されたＤＮＡシーケンサーの性能を維持することができる。また、画像の取得に失敗する回数が低減されるので、高価な試薬と検査時間を費やさなくてすむ。
【００７９】
＜第３の実施形態＞
画像の撮影中でも、分注ロボット１０９の分注作業、ステージ移動等によって振動が発生している。各機構からの振動は、筐体１０３を介して、光学検出部１０５（画像取得装置１１４）およびステージ１１２に伝わる。したがって、画像取得装置１１４とステージ１１２の間に相対的な位置ずれが発生し、これによって画像取得装置１１４が撮影する画像にずれが発生する。
【００８０】
ここでは、式（８）を用いて、筐体１０３の加速度ａ１（ｔ）と画像ずれの関係式を求める。この関係式に画像ずれの許容値を代入し、筐体１０３の加速度の上限値ａ１ｍａｘを計算する。よって、撮影中、筐体１０３に設置した加速度センサ２０７が上限値ａ１ｍａｘ以上を検出した場合は、加速度が上限値ａ１ｍａｘを下回るまで画像取得装置１１４の撮影を停止する。これにより、画像ずれは許容値を満足する。以下、筐体１０３の加速度ａ１（ｔ）と画像ずれの関係式（１９）を導出する。
【００８１】
図７は、筐体の加速度を検出することで画像ずれを抑制するフローチャートである。ここでは、Ｓ７０５の処理のみが、図６のＳ６０５と異なる。
【００８２】
Ｓ７０５では、ステージ制御部２０９または分注ロボット制御部２０８が、振動源となるステージ移動装置または分注ロボット移動装置等の駆動力が小さくなるように制御する。ステージ移動装置の制御は、例えば、（パルス信号で制御された）パルスモータの回転角変位をボールネジにより直線変位に変換して駆動位置決めを行う方法が挙げられる。分注ロボットの制御は、例えば、パルスモータの回転角変位をタイミングベルトにより直線変位に変換して駆動位置決めする方法が挙げられる。
【００８３】
以上の通り、本実施形態では、式（１９）にステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）、固有角速度ω、バネ定数ｋ、減衰定数ｃを代入し、筐体１０３の加速度の上限値ａ１ｍａｘを計算する。そして、筐体１０３の加速度が上限値ａ１ｍａｘ以上とならないように、振動源となるステージ移動装置または分注ロボット移動装置の駆動力を制御する。また、撮影中に筐体１０３の加速度が上限値を超える場合には、加速度超過のエラーをログファイルに記録する。
【００８４】
これにより、分注ロボット１０９による分注作業を行いながら、画像を撮影しても、蛍光の位置ずれを許容値以内に抑えることが可能となる。また、ステージ移動を行った後、装置が完全に停止するまで待たなくても、迅速に画像を撮影することができ、蛍光の位置ずれを許容値以内に抑えることが可能となる。また、ログファイル２１３を参照すれば、撮影した画像と、加速度超過との因果関係を確認することができる。
【００８５】
したがって、剛性が不明の作業台に設置されたＤＮＡシーケンサーの性能を維持することができる。また、画像の取得に失敗する回数が低減されるので、高価な試薬と検査時間を費やさなくてすむ。
【符号の説明】
【００８６】
１０３・・・筐体
１０５・・・光学検出部
１０６・・・試薬ラック
１０９・・・分注ロボット
１１０・・・廃液タンク
１１２・・・ステージ
１１４・・・画像取得部
２０１・・・作業台
２０２・・・キャスター
２０３・・・制御手段
２０４・・・表示装置
２０５・・・ポインティングデバイス
２０６・・・記憶装置
２０７・・・加速度センサ
２０８・・・分注ロボット制御部
２０９・・・ステージ制御部
２１０・・・画像制御部
２１１・・・加速度検出部
２１２・・・演算部
２１３・・・ログファイル
２１４・・・関係式
【技術分野】
【０００１】
本発明は、遺伝子配列解析装置および遺伝子配列解析方法に関し、例えば、筐体内に、試薬を分注する分注ロボットと、前記試薬を載置し、二次元に移動するステージと、前記試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置に関する。
【背景技術】
【０００２】
遺伝子配列解析装置（以下、「ＤＮＡシーケンサー」という）は、ＤＮＡの塩基配列を決定する装置である。ＤＮＡシーケンサーでは、１ミクロンのビーズ上にＤＮＡを固定し、エマルジョンＰＣＲによりＤＮＡを増幅させ、これらのビーズをフローセルのセルに固定し、セル内に試薬を送液し、ＤＮＡの伸長あるいはライゲーション反応によって蛍光基質をＤＮＡ鎖に取り込ませる。この操作により、フローセルに多数のビーズが固定されるため、励起光でフローセルの撮影領域を照射すると、ビーズから蛍光が発せられる。全ビーズの蛍光情報は、ＸＹステージを逐次移動させ、ＣＣＤカメラでフローセルの撮影領域を撮影することによって得られる。このようにして得られた蛍光情報に基づき、ＤＮＡの塩基配列を決定する。
【０００３】
ＤＮＡシーケンサーで行う塩基配列の決定を低コスト化するには、塩基配列を決定するまでの時間を短縮し、スループットを増加させることが求められる。また、スループットを増加させる一方、装置の性能は維持しなければならない。装置の性能を低下させる要因として、装置の振動が挙げられる。装置が振動すると、ＣＣＤカメラの位置とフローセルの位置との間に相対的な位置ずれが発生する。この状態のままフローセルを撮影すれば、ＣＣＤカメラの視野に散在する蛍光の位置がずれて撮影されるので、正確な蛍光情報が得られなくなる。また、試料をフローセルに注入する作業に支障が生じる。したがって、装置の性能を維持するには、装置の振動を抑制する必要がある。
【０００４】
特許文献１には、装置の構造解析と特性評価を実施することで、所望の性能に達しない場合は設計データを修正し、短期間に所望の性能を実現するリソグラフィ装置が開示されている。
【０００５】
特許文献２には、装置が設置される床を加振させることにより床構造の応答情報を測り、床の剛性が不足している場合には床構造を補強することで、所望の性能を実現するリソグラフィの製造方法が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第３９２８２００号公報
【特許文献２】特開２００５−２４３８８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
蛍光の撮影中、ＤＮＡシーケンサーには、様々な要因で振動が発生する。例えば、装置に人がぶつかって発生する振動、オペレータが情報を入力することで発生する振動、装置内のステージ移動や分注動作によって発生する振動である。
【０００８】
また、ＤＮＡシーケンサーは、一般的な研究室の作業台等に設置されることが多い。ＤＮＡシーケンサーが一般的な作業台に設置される場合、設置される作業台の剛性は不明である。例えば、作業台の剛性が大きいほど振動は起きにくいが、作業台の剛性が小さいほど、装置の振動は大きくなり、性能が低下する傾向がある。
【０００９】
したがって、設置される作業台の剛性を問わず、装置の性能を維持しなければならない。
【００１０】
上記特許文献１および２は、いずれも作業台の剛性を考慮していない。よって、剛性の小さい作業台に設置された場合は、振動により性能が低下するという問題がある。
【００１１】
本発明は上記課題に鑑みてなされたものであり、剛性が不明の作業台に設置された場合でも、簡易な構成で、性能を維持することができる遺伝子配列解析装置および遺伝子配列解析方法を提供する。
【課題を解決するための手段】
【００１２】
上記課題を解決するために本発明は、画像に生じる蛍光の位置ずれの許容値から機構の駆動力の上限値を求め、この上限値を上回らないように振動源となる機構を制御する。
【００１３】
また、蛍光の位置ずれの許容値から筐体の加速度の上限値を求め、筐体に設置した加速度センサがこの上限値を上回るときは、所定の整定時間を待つか、振動源となる機構を制御する。
【００１４】
すなわち、本発明は、筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力のうち少なくとも一方の上限値を計算する演算部と、前記分注ロボットに要する駆動力が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または前記ステージ移動に要する駆動力が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、を備えることを特徴とする。
【００１５】
また、本発明は、筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算する演算部と、前記筐体の加速度を検出し、前記検出した加速度が前記上限値を上回る場合、所定の整定時間の経過を待つ加速度検出部と、を備えることを特徴とする。
【００１６】
また、本発明は、筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算する演算部と、前記筐体の加速度を検出する加速度検出部と、前記加速度が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または前記加速度が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、を備えることを特徴とする。
【発明の効果】
【００１７】
剛性が不明の作業台に設置された場合でも、簡易な構成で、性能を維持することができる遺伝子配列解析装置および遺伝子配列解析方法を提供する。
【図面の簡単な説明】
【００１８】
【図１】図１は、一般的なＤＮＡシーケンサーの概略構成図である。
【図２】図２は、一般的な作業台２０１に設置された本発明の実施形態に係るＤＮＡシーケンサーである。
【図３】図３は、本発明の実施形態に係るＤＮＡシーケンサーのシミュレーションモデルである。
【図４】図４は、本発明の実施形態に係る固有角速度ωを求める概念図である。
【図５】図５は、分注ロボットの駆動力を制御することで画像に生じる蛍光の位置ずれを抑制する本発明の実施形態に係るフローチャートである。
【図６】図６は、筐体の加速度を検出することで画像に生じる蛍光の位置ずれを抑制する本発明の実施形態に係るフローチャートである。
【図７】図７は、筐体の加速度を検出することで画像に生じる蛍光の位置ずれを抑制する本発明の実施形態に係るフローチャートである。
【発明を実施するための形態】
【００１９】
本発明は、筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置および遺伝子配列解析方法に関するものである。以下、添付図面を参照して本発明の実施形態について説明する。ただし、本実施形態は、本発明を実施するための一例に過ぎず、本発明の技術的範囲を限定するものではないことに注意すべきである。また、各図において共通の構成には同一の参照番号が付されている。
【００２０】
＜ＤＮＡシーケンサー＞
図１は、一般的なＤＮＡシーケンサーの概略構成図である。ＤＮＡシーケンサーは、底部に防振ゴム１０１と側面に排熱ファン１０２を有する筐体１０３内に、除振マウント１０４上に設置された光学検出部１０５、試薬ラック１０６に並べられた試薬１０７を注入し、フローセル１０８に分注する分注ロボット１０９、使用済みの試薬を回収する廃液タンク１１０、底部に防振ゴム１０１を有し、試薬を送液するポンプ１１１を備える。光学検出部１０５は、フローセル１０８を載置するステージ１１２、ステージ１１２上の試薬を照射する光源１１３、試薬から発する蛍光を撮影するＣＣＤカメラ等の画像取得装置１１４を備える。分注ロボット１０９およびステージ１１２は、分注ロボット移動装置およびステージ移動装置によって少なくとも水平方向に移動可能である。防振ゴム１０１は、常時駆動しているポンプ１１１、排熱ファン１０２からの振動が筐体１０３に接する他の機構に伝達されないように、振動を減衰させるものである。また、マウント１０４は、筐体からの振動が光学検出部に伝達されないように、振動を減衰させるものである。
【００２１】
このＤＮＡシーケンサーは、一般的な作業台に設置され、そこでＤＮＡ検出が行われる。このような使用環境で、特に問題になるのは、分注ロボット１０９の分注作業をしながら、または高速でステージ移動をしながら、フローセル上のビーズから発せられる蛍光を撮影するときである。作業台の剛性が低ければ、分注作業やステージ移動によって筐体に発生する振動が大きくなり、この振動によってＣＣＤカメラの視野で蛍光の位置がずれ、正確な蛍光情報が得られなくなるためである。蛍光の位置ずれは、画像取得装置１１４の分解能に対して、充分に小さな値にする必要がある。例えば、画像取得装置１１４の分解能が０．３８μｍであれば、蛍光の位置ずれの許容値は０．１μｍ以下となる。
【００２２】
図２は、一般的な作業台２０１に設置された本発明の実施形態に係るＤＮＡシーケンサーである。作業台２０１は、ゴム製のキャスタ２０２を有しており、キャスタ２０２を介して床に設置される。本発明のＤＮＡシーケンサーは、ＤＮＡシーケンサを制御する制御手段２０３、制御手段２０３に接続された表示装置２０４、ポインティングデバイス２０５、および記憶装置２０６、筐体に設置された加速度センサ２０７を設けている。その他の構成は、特記しない限り図１と同じである。
【００２３】
制御手段２０３は、分注ロボット移動装置の駆動力を制御する分注ロボット制御部２０８、ステージ移動装置を制御するステージ制御部２０９、画像取得装置１１４を制御する画像制御部２１０、加速度センサ２０７からの加速度を検出する加速度検出部２１１、所定の関係式を算出し、関係式に基づき計算を行う演算部２１２を備える。
【００２４】
記憶装置２０６には、ＤＮＡシーケンサーによる検査中の操作履歴およびエラー履歴を記録したログファイル２１３、演算部２１２が算出した関係式２１４、および関係式に代入する測定値２１５が格納されている。
【００２５】
＜分注ロボットによる振動の影響＞
ここでは、分注ロボットによって発生した振動が、蛍光の位置ずれをもたらす原理について説明する。
【００２６】
この原理は、振動する機構とそれに接する機構との間には、所定のバネ定数および減衰定数があるとみなせるという概念に基づく。分注ロボット１０９は筐体１０３に設置されているため、分注ロボット１０９が駆動力を発生すれば、筐体１０３の振動は増幅し、所定の固有角周波数ωで振動する。また、分注ロボット１０９が駆動力を停止すれば筐体１０３の振動は減衰する。したがって、分注ロボット１０９と筐体１０３の間には、所定のバネ定数および減衰定数があるとみなせる。同様の理由で、筐体１０３と作業台２０１の間、筐体１０３と光学検出部１０５の間、光学検出部１０５とステージ１１２の間には、それぞれ所定のバネ定数および減衰定数があるとみなせる。
【００２７】
図３は、本発明の実施形態に係るＤＮＡシーケンサーのシミュレーションモデルである。図示するように、各機構と機構のバネ定数および減衰定数については、床と作業台２０１の間をｋ０、ｃ０、作業台２０１と筐体１０３の間をｋ１およびｃ１、筐体１０３と光学検出器１０５の間をｋ２およびｃ２、光学検出器１０５とステージ１１２の間をｋ３およびｃ３、筐体１０３と分注ロボット１０９の間をｋ４およびｃ４とする。各機構の質量については、作業台２０１をｍ０、筐体１０３をｍ１、光学検出器１０５をｍ２、ステージ１１２をｍ３、分注ロボット１０９をｍ４とする。各機構の中心座標からの変位については、作業台２０１をｘ０、筐体１０３をｘ１、光学検出器１０５をｘ２、ステージ１１２をｘ３、分注ロボット１０９をｘ４とする。各機構が発生する駆動力については、作業台２０１をｈ０、筐体１０３をｈ１、光学検出器１０５をｈ２、ステージ１１２をｈ３、分注ロボット１０９をｈ４とする。
【００２８】
以上の設定に基づき、各機構ごとに運動方程式ｆ０〜ｆ４を立てると、次式（１）が成立する。
【００２９】
【数１】
【００３０】
この運動方程式を微分方程式によって解くため、駆動力ｈ（ｔ）を次式（２）のように置く。
【数２】
【００３１】
式（２）より、式（３）が導ける。
【数３】
【００３２】
ここで、変位Ｘおよび変位の位相差φは次式から求められる。
【数４】
【００３３】
式（２）および（３）を式（１）に代入し、運動方程式を次式（５）のように書き換える。
【００３４】
【数５】
【００３５】
式（５）を行列に書き換えると次式（６）になる。
【数６】
【００３６】
ここで、Ａ〜Ｍは、次式（７）のように置いた。
【数７】
【００３７】
式（６）の逆行列をＮ−１とおくと、次式（８）になる。
【数８】
【００３８】
したがって、各機構の変位Ｘ０〜Ｘ４は、Ｎ−１の因子となるバネ定数ｋ０〜ｋ４、減衰定数ｃ０〜ｃ４、および各機構で共通の固有角速度ω、ならびに各機構の駆動力Ｈ０〜Ｈ４を求めることで、計算可能である。なお、Ｘ０〜Ｘ４は複素数で表記される。
【００３９】
次に、バネ定数ｋ、減衰定数ｃ、および固有角速度ωの求め方について説明する。
図４は、本発明の実施形態に係る固有角速度ωを求める概念図である。固有角速度ωは、筐体１０３に設置された加速度センサ２０７によって、周波数ごとに測定する。分注ロボット１０９を一定の駆動力で動作させ、周波数を変えながら（たとえば１〜２０Ｈｚ）、分注ロボット１０９を水平方向に正弦波駆動（スイープ加振）をさせる。周波数ごとの加速度測定結果から、加速度振幅が最大となる固有周波数ｆを求め、固有角速度ωを計算する。図４の場合は、固有周波数ｆは１０Ｈｚ、固有角速度ωは６２．８ｒａｄ／ｓである。
【００４０】
バネ定数ｋは、次式（９）により計算可能である。
【数９】
【００４１】
次に、減衰係数ｃを求める。まず、加速度の測定結果から、固有周波数における共振倍率Ｑを求める。
【００４２】
【数１０】
【００４３】
減衰比ζは、次式（１１）から求められる。
【数１１】
【００４４】
また、作業台のキャスター（またはアジャスター）の設置面はゴムであることが多いので、減衰比ζをゴムの損失係数（η＝０．１）から求めても良い。この場合、減衰比は次式（１１）のようになる。
【００４５】
【数１２】
【００４６】
減衰係数ｃは、減衰比ζを求めた後、次式により計算可能である。
【数１３】
【００４７】
以上の式から求めた各パラメータの値の一例を挙げておく。
作業台の固有周波数 ｆ＝１０Ｈｚ
作業台の固有角速度 ω＝ ２πｆ＝ ６２．８ （ｒａｄ／ｓ）
作業台の質量 ｍ１＝１００（ｋｇ）
作業台のばね定数 ｋ１＝４００（ｋＮ／ｍ）
作業台の共振倍率 Ｑ＝０．１／０．０１＝１０
作業台の減衰比 ζ＝０．０５
作業台の減衰係数 ｃ１＝１９００（Ｎ／（ｍ／ｓ））
【００４８】
＜第１の実施形態＞
画像の撮影中でも、分注ロボットは分注作業を行っている。分注ロボットからの振動は、筐体１０３、光学検出部１０５（画像取得装置１１４）、ステージ１１２の順に伝わる。したがって、画像取得装置１１４とステージ１１２との間に相対的な位置ずれが発生し、これによって画像取得装置１１４が撮影する時に蛍光の位置がずれる。
【００４９】
ここでは、式（７）を用いて、分注ロボット１０９の駆動力ｈ４（ｔ）と蛍光の位置ずれの関係式を求める。この関係式に位置ずれの許容値を代入し、分注ロボット１０９の駆動力の上限値ｈ４ｍａｘを計算する。よって、撮影中、分注ロボット１０９の駆動力が上限値ｈ４ｍａｘ以上とならないように、振動源となる分注ロボット１０９の駆動力を制御すれば、位置ずれは許容値を満足する。以下、分注ロボットの駆動力ｈ４（ｔ）と位置ずれの関係式（１５）を導出する。
【００５０】
蛍光の位置ずれは、ステージの変位と光学検出器の変位との差（ｘ３−ｘ２）に等しいとみなせる。また、分注ロボット以外の機構は駆動力を発生しないとし、Ｈ０〜Ｈ３＝０とする。これと、式（２）および（８）から、次式（１４）が求まる。ここで、Ｂ（ω）は、ｋおよびｃの関数とする。
【００５１】
【数１４】
【００５２】
また、式（３）より、次式（１５）が求まる。
【数１５】
【００５３】
式（１４）および（１５）より、次式（１６）が導ける。
【数１６】
【００５４】
ここで、Ｘ２〜Ｘ４およびＢ（ω）は、式（８）またはｋ、ｃ、ωにより計算可能である。したがって、式（１６）に、ステージの変位と光学検出器の変位との差（ｘ３−ｘ２）の許容値（例えば０．１μｍ）を代入すれば、分注ロボットの駆動力ｈ４（ｔ）の上限値が計算可能である。
【００５５】
図５は、分注ロボット１０９の駆動力を制御することで画像ずれを抑制するフローチャートである。
【００５６】
Ｓ５００では、演算部２１２が、式（１６）等を用いて、分注ロボット１０９の駆動力ｈ４（ｔ）の上限値ｈ４ｍａｘを計算する。ステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）は、ユーザが入力するか、演算部２１２が画像の分解能から計算する。また、固有角速度ω、バネ定数ｋ、減衰定数ｃは、予め値が分かっている場合はその値を式（１６）に代入するが、値が分からなければ次のように求める。まず、分注ロボット制御部２０８が、分注ロボット１０９を一定の駆動力で動作させ、周波数を変えながら、分注ロボット１０９を水平方向に正弦波駆動（スイープ加振）させる。次に、加速度検出部２１１が、筐体に設けた加速度センサ２０７により、周波数ごとの加速度測定結果を求める。そして、演算部２１２が、加速度振幅が最大となる固有周波数ｆを求め、固有角速度ω、バネ定数ｋ、減衰定数ｃを計算する。
【００５７】
Ｓ５０１およびＳ５０２では、分注ロボット制御部２０８が、分注ロボット１０９の駆動を開始し、駆動力を制御する。分注ロボット１０９は、例えばタイミングベルトを介してパルスモータにより制御される。
【００５８】
Ｓ５０３では、分注ロボット制御部２０８が、予め計算した上限値ｈ４ｍａｘと駆動力とを比較する。駆動力が上限値ｈ４ｍａｘより小さければ、Ｓ５０４に進む。駆動力が上限値ｈ４ｍａｘ以上であれば、Ｓ５０５に進む。
【００５９】
Ｓ５０４では、画像制御部２１０が、画像取得装置１１４による撮影をする。
Ｓ５０５では、分注ロボット制御部２０８が、駆動力超過のエラーをログファイル２１３に記録する。その後、Ｓ５０２に戻り、分注ロボット制御部２０８が、駆動力が小さくなるように分注ロボット１０９を制御する。
【００６０】
以上の通り、本実施形態では、式（１６）にステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）、固有角速度ω、バネ定数ｋ、減衰定数ｃを代入し、分注ロボット１０９の駆動力の上限値ｈ４ｍａｘを計算する。そして、分注ロボット１０９の駆動力が上限値ｈ４ｍａｘ以上とならないように、振動源となる分注ロボット１０９の駆動力を制御する。また、撮影中に分注ロボット１０９の駆動力が上限値を超える場合には、駆動力超過のエラーをログファイル２１３に記録する。
【００６１】
これにより、分注ロボットによる分注作業を行いながら画像を撮影しても、画像ずれを許容値以内に抑えることが可能となる。また、ログファイルを参照すれば、撮影した画像と、駆動力超過との因果関係を確認することができる。
【００６２】
したがって、剛性が不明の作業台に設置されたＤＮＡシーケンサーの性能を維持することができる。また、画像の取得に失敗する回数が低減されるので、高価な試薬と検査時間を費やさなくてすむ。
【００６３】
なお、ここでは、分注ロボットの駆動力を制御する方法について述べたが、式（７）に基づき、その他の振動源となる機構の駆動力を制御しても良い。例えば、ステージ移動装置の駆動力ｈ３（ｔ）も制御可能である。
【００６４】
＜第２の実施形態＞
画像の撮影中でも、装置に人がぶつかったり、オペレータが情報を入力したりすることによって装置に振動が発生している。装置に発生した振動は、筐体１０３を介して、光学検出部１０５（画像取得装置１１４）およびステージ１１２に伝わる。したがって、画像取得装置１１４とステージ１１２との間に相対的な位置ずれが発生し、これによって画像取得装置１１４が撮影する画像にずれが発生する。
【００６５】
ここでは、式（７）を用いて、筐体の加速度ａ１（ｔ）と画像ずれの関係式を求める。この関係式に画像ずれの許容値を代入し、筐体１０３の加速度の上限値ａ１ｍａｘを計算する。よって、撮影中、筐体１０３に設置した加速度センサ２０７が上限値ａ１ｍａｘ以上を検出した場合は、加速度が上限値ａ１ｍａｘを下回るまで画像取得装置１１４の撮影を停止する。これにより、画像ずれは許容値を満足する。以下、筐体１０３の加速度ａ１（ｔ）と画像ずれの関係式（１９）を導出する。
【００６６】
画像ずれは、ステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）に等しいとみなせる。また、式（３）から、次式（１７）および（１８）が求まる。
【００６７】
【数１７】
【００６８】
【数１８】
【００６９】
式（１７）および（１８）から、次式（１９）が求まる。
【数１９】
【００７０】
ここで、ωは測定から求め、Ｘ１〜Ｘ３は式（７）により計算可能である。したがって、式（１９）に、ステージの変位と光学検出器１０５の変位との差（ｘ３−ｘ２）の許容値（例えば０．１μｍ）を代入すれば、筐体１０３の加速度ａ１（ｔ）の上限値が計算可能である。
【００７１】
図６は、筐体１０３の加速度を検出することで画像ずれを抑制するフローチャートである。
【００７２】
Ｓ６００では、演算部２１２が、式（１９）等を用いて、筐体１０３の加速度ａ１（ｔ）の上限値をａ１ｍａｘ計算する。なお、ステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）は、ユーザが入力してもよいし、演算部２１２が画像の分解能から計算しても良い。固有角速度ω、バネ定数ｋ、減衰定数ｃは、前述のＳ５００で説明したとおりである。
【００７３】
Ｓ６０１および６０２では、加速度検出部２１１が、筐体１０３の加速度を検出し、予め計算した上限値ａ１ｍａｘと検出した加速度を比較する。検出した加速度が上限値より小さければ、Ｓ６０３に進む。検出した加速度が上限値以上になれば、Ｓ６０４に進む。
【００７４】
Ｓ６０３では、画像制御部２１０が、画像取得装置１１４による撮影をする。
Ｓ６０４では、加速度検出部２１１が、加速度超過のエラーをログファイル２１３に記録する。
【００７５】
Ｓ６０５では、加速度検出部２１１が、所定の整定時間の経過を待った後、Ｓ６０１に戻り、再度筐体の加速度を検出する。ここで整定時間とは、例えば、筐体１０３の固有周波数ｆから求めた周期の定数倍に設定してもよいし、タイマー等で任意の時間に設定してもよい。なお、所定の整定時間の経過を待たずに、検査を中止してもよい。
【００７６】
以上の通り、本実施形態では、式（１９）にステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）、固有角速度ω、バネ定数ｋ、減衰定数ｃを代入し、筐体１０３の加速度の上限値ａ１ｍａｘを計算する。そして、筐体に設置した加速度センサが上限値ａ１ｍａｘ以上を検出した場合は、加速度が上限値ａ１ｍａｘを下回るまで画像取得装置１１４の撮影を停止する。また、撮影中に筐体１０３の加速度が上限値を超える場合には、加速度超過のエラーをログファイル２１３に記録する。
【００７７】
これにより、装置に人がぶつかったり、オペレータが情報を入力したりすることによって装置に振動が発生しても、画像ずれを許容値以内に抑えることが可能となる。また、ログファイルを参照すれば、撮影した画像と、加速度超過との因果関係を確認することができる。
【００７８】
したがって、剛性が不明の作業台に設置されたＤＮＡシーケンサーの性能を維持することができる。また、画像の取得に失敗する回数が低減されるので、高価な試薬と検査時間を費やさなくてすむ。
【００７９】
＜第３の実施形態＞
画像の撮影中でも、分注ロボット１０９の分注作業、ステージ移動等によって振動が発生している。各機構からの振動は、筐体１０３を介して、光学検出部１０５（画像取得装置１１４）およびステージ１１２に伝わる。したがって、画像取得装置１１４とステージ１１２の間に相対的な位置ずれが発生し、これによって画像取得装置１１４が撮影する画像にずれが発生する。
【００８０】
ここでは、式（８）を用いて、筐体１０３の加速度ａ１（ｔ）と画像ずれの関係式を求める。この関係式に画像ずれの許容値を代入し、筐体１０３の加速度の上限値ａ１ｍａｘを計算する。よって、撮影中、筐体１０３に設置した加速度センサ２０７が上限値ａ１ｍａｘ以上を検出した場合は、加速度が上限値ａ１ｍａｘを下回るまで画像取得装置１１４の撮影を停止する。これにより、画像ずれは許容値を満足する。以下、筐体１０３の加速度ａ１（ｔ）と画像ずれの関係式（１９）を導出する。
【００８１】
図７は、筐体の加速度を検出することで画像ずれを抑制するフローチャートである。ここでは、Ｓ７０５の処理のみが、図６のＳ６０５と異なる。
【００８２】
Ｓ７０５では、ステージ制御部２０９または分注ロボット制御部２０８が、振動源となるステージ移動装置または分注ロボット移動装置等の駆動力が小さくなるように制御する。ステージ移動装置の制御は、例えば、（パルス信号で制御された）パルスモータの回転角変位をボールネジにより直線変位に変換して駆動位置決めを行う方法が挙げられる。分注ロボットの制御は、例えば、パルスモータの回転角変位をタイミングベルトにより直線変位に変換して駆動位置決めする方法が挙げられる。
【００８３】
以上の通り、本実施形態では、式（１９）にステージ１１２の変位と光学検出器１０５の変位との差（ｘ３−ｘ２）、固有角速度ω、バネ定数ｋ、減衰定数ｃを代入し、筐体１０３の加速度の上限値ａ１ｍａｘを計算する。そして、筐体１０３の加速度が上限値ａ１ｍａｘ以上とならないように、振動源となるステージ移動装置または分注ロボット移動装置の駆動力を制御する。また、撮影中に筐体１０３の加速度が上限値を超える場合には、加速度超過のエラーをログファイルに記録する。
【００８４】
これにより、分注ロボット１０９による分注作業を行いながら、画像を撮影しても、蛍光の位置ずれを許容値以内に抑えることが可能となる。また、ステージ移動を行った後、装置が完全に停止するまで待たなくても、迅速に画像を撮影することができ、蛍光の位置ずれを許容値以内に抑えることが可能となる。また、ログファイル２１３を参照すれば、撮影した画像と、加速度超過との因果関係を確認することができる。
【００８５】
したがって、剛性が不明の作業台に設置されたＤＮＡシーケンサーの性能を維持することができる。また、画像の取得に失敗する回数が低減されるので、高価な試薬と検査時間を費やさなくてすむ。
【符号の説明】
【００８６】
１０３・・・筐体
１０５・・・光学検出部
１０６・・・試薬ラック
１０９・・・分注ロボット
１１０・・・廃液タンク
１１２・・・ステージ
１１４・・・画像取得部
２０１・・・作業台
２０２・・・キャスター
２０３・・・制御手段
２０４・・・表示装置
２０５・・・ポインティングデバイス
２０６・・・記憶装置
２０７・・・加速度センサ
２０８・・・分注ロボット制御部
２０９・・・ステージ制御部
２１０・・・画像制御部
２１１・・・加速度検出部
２１２・・・演算部
２１３・・・ログファイル
２１４・・・関係式
【特許請求の範囲】
【請求項１】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、
作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力のうち少なくとも一方の上限値を計算する演算部と、
前記分注ロボットに要する駆動力が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または
前記ステージ移動に要する駆動力が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、
を備えることを特徴とする遺伝子配列解析装置。
【請求項２】
さらに、前記筐体の加速度を検出する加速度検出部を備え、
前記関係式は、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力と、前記画像取得装置の変位と前記ステージの変位の差分値との関係式であり、
前記加速度検出部は、前記分注ロボットまたはステージ移動の少なくとも一方によって筐体に生じる加速度を検出し、
前記演算部は、前記検出した加速度から計算した固有角速度、ならびに前記作業台と床との間および前記筐体内の各機構間のバネ定数および減衰定数を前記関係式に代入することで、前記上限値を計算することを特徴とする請求項１記載の遺伝子配列解析装置。
【請求項３】
前記分注ロボット制御部またはステージ制御部の少なくとも一方は、前記駆動力が前記上限値を上回る場合は、駆動力超過のエラーをログファイルに記録することを特徴とする請求項１記載の遺伝子配列解析装置。
【請求項４】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、
作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算する演算部と、
前記筐体の加速度を検出し、前記検出した加速度が前記上限値を上回る場合、所定の整定時間の経過を待つ加速度検出部と、
を備えることを特徴とする遺伝子配列解析装置。
【請求項５】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、
作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算する演算部と、
前記筐体の加速度を検出する加速度検出部と、
前記加速度が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または
前記加速度が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、
を備えることを特徴とする遺伝子配列解析装置。
【請求項６】
前記関係式は、前記加速度と、前記画像取得装置の変位と前記ステージの変位の差分値との関係式であり、
前記加速度検出部は、前記筐体に生じる加速度を検出し、
前記演算部は、前記検出した加速度から計算した固有角速度、ならびに前記作業台と床との間および前記筐体内の各機構間のバネ定数および減衰定数を前記関係式に代入することで、前記上限値を計算することを特徴とする請求項４または５記載の遺伝子配列解析装置。
【請求項７】
前記加速度検出部は、前記加速度が前記上限値を上回る場合は、加速度超過のエラーをログファイルに記録することを特徴とする請求項４または５記載の遺伝子配列解析装置。
【請求項８】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析方法であって、
演算部が、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力のうち少なくとも一方の上限値を計算するステップと、
分注ロボット制御部が、前記分注ロボットに要する駆動力が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制するステップ、または
ステージ制御部が、前記ステージ移動に要する駆動力が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステップのうち少なくとも一方のステップと、
を備えることを特徴とする遺伝子配列解析方法。
【請求項９】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析方法であって、
演算部が、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算するステップと、
加速度検出部が、前記筐体の加速度を検出し、前記検出した加速度が前記上限値を上回る場合、所定の整定時間の経過を待つステップと、
を備えることを特徴とする遺伝子配列解析方法。
【請求項１０】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する際に画像取得装置とを備える遺伝子配列解析方法であって、
演算部が、作業台に載置された状態で、前記画像取得装置が撮影する画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算するステップと、
加速度検出部が、前記筐体の加速度を検出するステップと、
分注ロボット制御部が、前記加速度が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制するステップ、または
ステージ制御部が、前記加速度が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステップのうち少なくとも一方のステップと、
を備えることを特徴とする遺伝子配列解析方法。
【請求項１】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、
作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力のうち少なくとも一方の上限値を計算する演算部と、
前記分注ロボットに要する駆動力が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または
前記ステージ移動に要する駆動力が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、
を備えることを特徴とする遺伝子配列解析装置。
【請求項２】
さらに、前記筐体の加速度を検出する加速度検出部を備え、
前記関係式は、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力と、前記画像取得装置の変位と前記ステージの変位の差分値との関係式であり、
前記加速度検出部は、前記分注ロボットまたはステージ移動の少なくとも一方によって筐体に生じる加速度を検出し、
前記演算部は、前記検出した加速度から計算した固有角速度、ならびに前記作業台と床との間および前記筐体内の各機構間のバネ定数および減衰定数を前記関係式に代入することで、前記上限値を計算することを特徴とする請求項１記載の遺伝子配列解析装置。
【請求項３】
前記分注ロボット制御部またはステージ制御部の少なくとも一方は、前記駆動力が前記上限値を上回る場合は、駆動力超過のエラーをログファイルに記録することを特徴とする請求項１記載の遺伝子配列解析装置。
【請求項４】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、
作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算する演算部と、
前記筐体の加速度を検出し、前記検出した加速度が前記上限値を上回る場合、所定の整定時間の経過を待つ加速度検出部と、
を備えることを特徴とする遺伝子配列解析装置。
【請求項５】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析装置であって、
作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算する演算部と、
前記筐体の加速度を検出する加速度検出部と、
前記加速度が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制する分注ロボット制御部、または
前記加速度が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステージ制御部のうち少なくとも一方と、
を備えることを特徴とする遺伝子配列解析装置。
【請求項６】
前記関係式は、前記加速度と、前記画像取得装置の変位と前記ステージの変位の差分値との関係式であり、
前記加速度検出部は、前記筐体に生じる加速度を検出し、
前記演算部は、前記検出した加速度から計算した固有角速度、ならびに前記作業台と床との間および前記筐体内の各機構間のバネ定数および減衰定数を前記関係式に代入することで、前記上限値を計算することを特徴とする請求項４または５記載の遺伝子配列解析装置。
【請求項７】
前記加速度検出部は、前記加速度が前記上限値を上回る場合は、加速度超過のエラーをログファイルに記録することを特徴とする請求項４または５記載の遺伝子配列解析装置。
【請求項８】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析方法であって、
演算部が、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記分注ロボットに要する駆動力または前記ステージ移動に要する駆動力のうち少なくとも一方の上限値を計算するステップと、
分注ロボット制御部が、前記分注ロボットに要する駆動力が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制するステップ、または
ステージ制御部が、前記ステージ移動に要する駆動力が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステップのうち少なくとも一方のステップと、
を備えることを特徴とする遺伝子配列解析方法。
【請求項９】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する画像取得装置とを備える遺伝子配列解析方法であって、
演算部が、作業台に載置された状態で、前記画像取得装置が撮影する際に画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算するステップと、
加速度検出部が、前記筐体の加速度を検出し、前記検出した加速度が前記上限値を上回る場合、所定の整定時間の経過を待つステップと、
を備えることを特徴とする遺伝子配列解析方法。
【請求項１０】
筐体内に、試料または試薬を分注する分注ロボットと、前記試料または試薬を載置し、二次元に移動するステージと、前記試料または試薬からの蛍光を撮影する際に画像取得装置とを備える遺伝子配列解析方法であって、
演算部が、作業台に載置された状態で、前記画像取得装置が撮影する画像に生じる蛍光の位置ずれを計算するための所定の関係式に基づき、前記筐体の加速度の上限値を計算するステップと、
加速度検出部が、前記筐体の加速度を検出するステップと、
分注ロボット制御部が、前記加速度が前記上限値を上回る場合、前記分注ロボットに要する駆動力を抑制するステップ、または
ステージ制御部が、前記加速度が前記上限値を上回る場合、前記ステージの移動に要する駆動力を抑制するステップのうち少なくとも一方のステップと、
を備えることを特徴とする遺伝子配列解析方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【公開番号】特開２０１２−９５５８２（Ｐ２０１２−９５５８２Ａ）
【公開日】平成２４年５月２４日（２０１２．５．２４）
【国際特許分類】
【出願番号】特願２０１０−２４５１５２（Ｐ２０１０−２４５１５２）
【出願日】平成２２年１１月１日（２０１０．１１．１）
【出願人】（５０１３８７８３９）株式会社日立ハイテクノロジーズ (4,325)
【Ｆターム（参考）】
【公開日】平成２４年５月２４日（２０１２．５．２４）
【国際特許分類】
【出願日】平成２２年１１月１日（２０１０．１１．１）
【出願人】（５０１３８７８３９）株式会社日立ハイテクノロジーズ (4,325)
【Ｆターム（参考）】
[ Back to top ]