３Ｄデータ解析装置および３Ｄデータ解析方法
【課題】多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、容易かつ直感的なゲーティング操作が可能で、高精度に解析を行うことができるデータ解析装置を提供。
【解決手段】微小粒子の測定データを保存するデータ格納部１３０と、前記測定データから独立した３種の変数を選択する入力部１４１と、前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成するデータ処理部１２０と、前記３Ｄ立体画像を表示する表示部１４２と、を有し、入力部１４１からの入力信号に基づいて前記座標空間内に分画領域を変更可能に設定し、前記３Ｄ立体画像中に表示する３Ｄデータ解析装置１を提供する。
【解決手段】微小粒子の測定データを保存するデータ格納部１３０と、前記測定データから独立した３種の変数を選択する入力部１４１と、前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成するデータ処理部１２０と、前記３Ｄ立体画像を表示する表示部１４２と、を有し、入力部１４１からの入力信号に基づいて前記座標空間内に分画領域を変更可能に設定し、前記３Ｄ立体画像中に表示する３Ｄデータ解析装置１を提供する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、３Ｄデータ解析装置および３Ｄデータ解析方法に関する。より詳しくは、解析対象とする微小粒子の指定を３Ｄ立体画像中で行うことができる３Ｄデータ解析装置等に関する。
【背景技術】
【０００２】
細胞、微生物およびリポソーム等の生体関連粒子や、ラテックス粒子、ゲル粒子および工業用粒子等の合成粒子などの微小粒子を解析するため、微小粒子の分散液を流路内に導入して、微小粒子を光学的、電気的あるいは磁気的に測定する微小粒子測定装置が用いられている。
【０００３】
一例として、合成粒子を大きさや形状に応じて判別するパーティクルアナライザーがある。パーティクルアナライザーで測定可能なパラメータ（変数）には、微小粒子の元素組成および粒径などがある。
【０００４】
また、生体関連粒子の解析では、フローサイトメータ（フローサイトメトリー）が用いられている。フローサイトメータで測定可能なパラメータには、微小粒子の前方散乱光（FS）、側方散乱（SS）、蛍光（FL）およびインピーダンスなどがある。前方散乱光（FS）、側方散乱（SS）および蛍光（FL）は細胞や微生物（以下、単に「細胞」という）の光学的特性を示すパラメータとして、インピーダンスは細胞の電気的特性を示すパラメータとして用いられる。
【０００５】
具体的には、まず、前方散乱光は、レーザー光の軸に対して前方向の小さい角度で散乱する光であり、細胞の表面で生じるレーザー光の散乱光や回折光、屈折光からなる。前方散乱光は、主として細胞の大きさを示すパラメータとして用いられている。次に、側方散乱は、レーザー光の軸に対して約９０度の角度で散乱する光であり、細胞の内部の顆粒や核などで生じるレーザー光の散乱光である。側方散乱は、主として細胞の内部構造を示すパラメータとして用いられている。また、蛍光は、細胞に標識された蛍光色素から発生する光であり、蛍光色素標識抗体が認識する細胞表面抗原の有無や蛍光色素が結合した核酸の量などを示すパラメータとして用いられている。また、インピーダンスは、電気抵抗法によって測定され、細胞の容積を示すパラメータとして用いられている。
【０００６】
フローサイトメータにおける測定データの解析のため、これらの測定パラメータを軸にとって各細胞の測定値をプロットすることにより、細胞集団中における細胞の特性分布を表す図を作成して表示するデータ解析装置が用いられている。測定パラメータを１つとした一次元分布図はヒストグラムとも称され、測定パラメータをＸ軸にとり、細胞数（カウント）をＹ軸にとったグラフとして作成される。また、測定パラメータを２つとした二次元分布図はサイトグラムとも称され、一方の測定パラメータをＸ軸とし、他方の測定パラメータをＹ軸とした座標面内に、各細胞をその測定値に基づいてプロットしていくことにより作成される。
【０００７】
サンプルとする細胞集団の中には、解析の対象としない不要な細胞も含まれているため、測定データの解析は、サンプルとする細胞集団の中から解析対象とする細胞小集団を選び出した後に行われている。解析対象とする細胞小集団の選択は、ヒストグラム上あるいはサイトグラム上において該細胞小集団が存在する領域を指定することによって行われている。この操作は、ヒストグラム上あるいはサイトグラム上で指定する領域内に目的とする細胞を囲い込むことから、「ゲーティング」と呼ばれている。
【０００８】
一つの測定パラメータを軸としたヒストグラム上あるいは一つの組み合わせの測定パラメータを軸としたサイトグラム上においては、解析対象とする細胞小集団と不要な細胞とが重なり合う領域に存在している場合がある。例えば、ヒト抹消血をサンプルとしてリンパ球の解析を行う際、前方散乱光（FS）と側方散乱（SS）を軸にとったサイトグラム上では、一部の単球がリンパ球と同一の領域に存在する場合がある。従って、ゲーティングを行う際、ユーザは、単球が囲い込まれないようにしてリンパ球のみが存在する領域を指定してやる必要がある。
【０００９】
不要な細胞を囲い込むことなく、解析対象とする細胞小集団のみを囲い込むことができるように領域を指定するため、従来、ユーザは複数のヒストグラムあるいはサイトグラムを参照しながらゲーティングを行う必要があった。フローサイトメータの性能向上に伴って測定可能なパラメータ数も増加しており、ユーザはより多くのヒストグラムあるいはサイトグラムを参照する必要が生じていた。また、この際、ユーザは、２つのサイトグラムを組み合わせた立体的な分布図（三次元分布図）を想起しながらゲーティング操作を行うことが求められていた。
【００１０】
ユーザのゲーティング操作を支援するため、特許文献１には、「分析物から第１、第２および第３の測定データを取得する測定データ取得手段と、前記第１、第２および第３の測定データを軸とし、分析物に含まれる有形成分の分布を示す３次元分布図を作成する３次元分布図作成手段と、前記３次元分布図上に分画領域を変更可能に設定する領域設定手段と、前記領域設定手段によって設定された前記分画領域に属する有形成分について、前記第１および第２の測定データを軸とする２次元分布図および前記第１の測定データを軸とする度数分布図の少なくとも一方を作成する参考用分布図作成手段と、を備える分析装置。」が提案されている（当該文献請求項９参照）。この分析装置によれば、３次元分布図とともに表示される２次元分布図（サイトグラム）および度数分布図（ヒストグラム）を参照しながら、３次元分布図上に分画領域を設定することができる。なお、この分析装置の３次元分布図は、ディスプレイ上に平面的に表示されるものであり、立体視されるものではない。
【００１１】
本発明に関連して、二眼式ステレオ立体画像技術（３Ｄ立体画像技術）について説明する。二眼式ステレオ立体映像では、まず、対象物を右眼および左眼で見たときの２つの映像を用意する。そして、これらの映像を同時に表示した上で、右眼用映像を右眼だけに、左眼用映像を左眼だけに提示する。これにより、三次元空間で対象物を見ている際に目に映る映像を再現し、ユーザに対象物を立体視させる。
【００１２】
立体視が可能な３Ｄディスプレイには、（ａ）メガネ方式、（ｂ）裸眼方式、（ｃ）ビューア方式が主に採用されている。（ａ）メガネ方式には、さらにアナグリフ方式、偏光フィルタ方式、時分割方式がある。また、（ｂ）裸眼方式にはパララックスバリア方式、レンティキュラ方式があり、（ｃ）ビューア方式にはステレオスコープ方式、ヘッドマウント方式がある。
【先行技術文献】
【特許文献】
【００１３】
【特許文献１】特開２００６−１７４９７号公報
【発明の概要】
【発明が解決しようとする課題】
【００１４】
上述のように、従来、フローサイトメータに用いられているデータ解析装置では、ゲーティングの際、不要な細胞を囲い込むことなく、解析対象とする細胞小集団のみを囲い込むため、ユーザは多数のヒストグラムあるいはサイトグラムを参照したり、２つのサイトグラムを組み合わせた立体的な分布図（三次元分布図）を想起したりしながら操作を行う必要があった。このため、従来のデータ解析装置では、ゲーティングの操作に時間がかかり、ユーザが操作に習熟していない場合には不要な細胞が混入してしまい、解析精度に問題が生じていた。
【００１５】
そこで、本発明は、多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、容易かつ直感的なゲーティング操作が可能で、高精度に解析を行うことができるデータ解析装置を提供することを主な目的とする。
【課題を解決するための手段】
【００１６】
上記課題解決のため、本発明は、微小粒子の測定データを保存するデータ格納部と、前記測定データから独立した３種の変数を選択する入力部と、前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成するデータ処理部と、前記３Ｄ立体画像を表示する表示部と、を有し、前記入力部からの入力信号に基づいて前記座標空間内に分画領域を変更可能に設定し、前記３Ｄ立体画像中に表示する３Ｄデータ解析装置を提供する。
この３Ｄデータ解析装置によれば、任意に選択した３種のパラメータを座標軸とした三次元分布図を立体視しながら、その座標空間において解析対象とする微小粒子ならびに微小粒子小集団のゲーティングを行うことができる。
【００１７】
また、本発明は、微小粒子の測定データから独立した３種の変数を選択する手順と、前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成する手順と、前記３Ｄ立体画像を表示させる手順と、前記３Ｄ立体画像中の前記座標空間内に分画領域を設定する手順と、を含む３Ｄデータ解析方法をも提供する。
【００１８】
本発明において、「微小粒子」には、細胞、微生物およびリポソーム等の生体関連粒子、ラテックス粒子、ゲル粒子、工業用粒子等の合成粒子などの微小粒子が広く含まれるものとする。
細胞には、動物細胞（血球系細胞など）および植物細胞が含まれる。微生物には、大腸菌等の細菌類、タバコモザイクウイルス等のウイルス類、イースト菌等の菌類などが含まれる。生体関連粒子には、各種細胞を構成する染色体、リポソーム、ミトコンドリア、オルガネラ(細胞小器官)なども含まれる。さらに、生体関連粒子には、核酸やタンパク質、これらの複合体などの生体関連高分子も包含され得るものとする。工業用粒子は、有機もしくは無機高分子材料、金属などであってもよい。有機高分子材料には、ポリスチレン、スチレン・ジビニルベンゼンおよびポリメチルメタクリレートなどが含まれる。無機高分子材料には、ガラス、シリカおよび磁性体材料などが含まれる。金属には、金コロイドおよびアルミなどが含まれる。これら微小粒子の形状は、一般には球形であるが、非球形であってもよい。また、微小粒子の大きさや質量なども特に限定されない。
【発明の効果】
【００１９】
本発明により、多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、容易かつ直感的なゲーティング操作が可能で、高精度に解析を行うことができるデータ解析装置が提供される。
【図面の簡単な説明】
【００２０】
【図１】フローサイトメータに連設された本発明に係る３Ｄデータ解析装置の構成を説明するブロック図である。
【図２】本発明に係る３Ｄデータ解析装置の機能的構成を説明するブロック図である。
【図３】本発明に係る３Ｄデータ解析装置により表示される三次元分布図を説明する模式図である。
【図４】本発明に係る３Ｄデータ解析装置により表示される二眼式ステレオ立体画像（３Ｄ立体画像）を説明する模式図である。
【図５】本発明に係る３Ｄデータ解析装置において設定され得る分画領域の具体例を説明する模式図である。
【図６】本発明に係る３Ｄデータ解析装置において設定され得る分画領域の具体例を説明する模式図である。
【図７】３Ｄ立体画像中の微小粒子に対応する図形の形状を説明する模式図である。
【図８】シェード処理された図形の立体観察像を説明する概念図である。
【図９】シェード処理の処理方法を説明する模式図である。
【図１０】座標軸の立体観察像を説明する概念図である。
【図１１】微小粒子に対応する図形を揺動させた動画像の立体観察像を説明する概念図である。
【発明を実施するための形態】
【００２１】
以下、本発明を実施するための好適な形態について図面を参照しながら説明する。なお、以下に説明する実施形態は、本発明の代表的な実施形態の一例を示したものであり、これにより本発明の範囲が狭く解釈されることはない。なお、説明は以下の順序で行う。
１．３Ｄデータ解析装置の構成
２．３Ｄ立体画像の表示
３．ゲーティング
４．３Ｄ立体画像の特徴
（４−１）図形の形状
（４−２）図形のシェード処理
（４−３）座標軸
（４−４）動画像
【００２２】
１．３Ｄデータ解析装置の構成
図１に、本発明に係る３Ｄデータ解析装置の装置構成を示す。ここでは、３Ｄデータ解析装置を微小粒子測定装置に連設して微小粒子解析システムとして構成した実施形態を示す。また、図２に、この微小粒子解析システムの機能的構成を示す。以下では、微小粒子測定装置としてフローサイトメータを用いる場合を例に説明する。
【００２３】
図中符号１で示す３Ｄデータ解析装置は、通信ケーブル４によりフローサイトメータ２と接続されて微小粒子解析システム３を構成している。３Ｄデータ解析装置１は、中央処理装置（ＣＰＵ）１０、メモリ２０、ハードディスク３０、ユーザインターフェイスなどを含んでいる。ハードディスク３０には、３Ｄデータ解析プログラム３１、微小粒子の測定データ３２およびオペレーティングシステム（ＯＳ）３３などが格納・保持されている。ユーザインターフェイスには、ユーザからの情報の入力を受け付けるマウス４１およびキーボード４２などの入力デバイスと、ユーザに対して情報を出力するディスプレイ４３およびプリンタ４４などの出力デバイスが含まれる。なお、マウス４１およびキーボード４２に替えて、あるいはこれらとともに、スティックコントローラおよびペンタブレットなどの入力デバイスを設けてもよい。
【００２４】
データ格納部１３０（ハードディスク３０）は、フローサイトメータ２から出力される微小粒子（細胞）の測定データ３２を保存する。フローサイトメータ２の入出力インターフェイス２５０から出力される測定データは、通信ケーブル４を介して３Ｄデータ解析装置１の入出力インターフェイス１５０に入力され、データ格納部３０（ハードディスク３０）に格納される。
【００２５】
測定データ３２は、データ処理部１２０において処理される。データ処理部１２０は、入力部１４１（マウス４１あるいはキーボード４２など）からのユーザの入力を受けて処理を開始する。すなわち、データ処理部１２０は、ユーザによって測定データ３２から独立した３種の変数（パラメータ）が選択され入力されると、選択されたパラメータを座標軸として微小粒子の特性分布を表す三次元分布図を作成する。三次元分布図は、選択されたパラメータを座標軸とした座標空間内に微小粒子をプロットすることにより作成される。微小粒子のプロットは、選択されたパラメータの測定値から各微小粒子の座標空間内の位置と図形を計算し、算出された位置に算出された図形を描画することによって行われる。
【００２６】
また、データ処理部１２０は、作成した三次元分布図の座標空間内に、入力部１４１からのユーザの入力信号に基づいて微小粒子をゲーティングするための分画領域を設定する。入力部１４１からのユーザの入力は、三次元分布図中に表示されるポインタをマウス４１あるいはキーボード４２などの入力デバイスを用いて移動、変形、拡大または縮小させることによって行われる。
【００２７】
さらに、データ処理部１２０は、設定された分画領域内に存在する微小粒子の測定データの処理（クラスター解析）を行う。クラスター解析は、例えば、全微小粒子に占める分画領域内に存在する微小粒子の割合、画領域内に存在する微小粒子の数、平均値、標準偏差および変動係数などの統計データが算出される。
【００２８】
ここで、「独立したパラメータ」とは、微小粒子の前方散乱光（FS）、側方散乱（SS）、蛍光（FL）およびインピーダンスなどから選択される互いに異なるパラメータを意味する。蛍光（FL）は、微小粒子に標識された蛍光色素の波長域毎に異なるパラメータとして扱うことができ、ＦＬ１，ＦＬ２〜ＦＬｎ（ｎは３以上の整数）などと表記される。独立した３種のパラメータとしては、例えば、前方散乱光（FS）、側方散乱（SS）および蛍光（FL1）の組み合わせや、前方散乱光（FS）、側方散乱（SS）およびインピーダンスの組み合わせが例示される。この他、独立した３種のパラメータは、測定データから任意に選択して組み合わせることができる。
【００２９】
データ処理部１２０により作成された三次元分布図と、その座標空間内に設定された分画領域とは、出力部１４２（ディスプレイ４３）に３Ｄ立体画像として表示される。この３Ｄ立体画像は、次に詳しく説明する二眼式ステレオ立体画像とされる。出力部１４２への３Ｄ立体画像の表示は、入力部１４１（マウス４１あるいはキーボード４２など）からのユーザの入力信号に基づき任意に回転、拡大あるいは縮小させて行ってもよい。また、データ処理部１２０によるクラスター解析の結果も、出力部１４２（ディスプレイ４３およびプリンタ４４など）に出力される。統計データは、出力部１４２に数値あるいは図表によって表示、印刷等されてよい。
【００３０】
フローサイトメータ２は、従来公知の装置と同様の構成あるいはこれを適宜改変した構成とでき、具体的には、制御部２１０、フロー系２２０、検出系２３０、入出力インターフェイス２５０等から構成される。
【００３１】
フロー系２２０は、フローセルやマイクロチップに形成された流路内において、微小粒子を含むサンプル液層流をシース液層流の中心に流し、層流中に微小粒子を一列に配列させる。検出系２３０は、流路を通流する微小粒子の特性を示すパラメータ値を取得する。具体的には、光学検出部２３１は、通流する微小粒子に光を照射し、微小粒子から生じる散乱光や蛍光などを検出し、その強度を取得する。光学検出部２３１は、レーザー光源、レンズ、ミラー、フィルター、ＣＣＤおよびＣＭＯＳ素子等のエリア撮像素子あるいはＰＭＴ（photo multiplier tube）などを含んでなる。また、電気検出部２３２は、通流する微小粒子に対向して配された電極を含んでなり、微小粒子のインピーダンス、容量値（キャパシタンス）およびインダクタンスなどを取得する。フローサイトメータ２は、解析の結果所望の特性を有すると判定された微小粒子を分取するための分取系２４０を備えていてもよい。分取系２４０には、例えば、微小粒子を含む液滴をフローセル外の空間に吐出し、液滴の移動方向を制御することにより、所望の微小粒子のみを容器に回収する方式を採用できる。
【００３２】
検出系２３０において検出された散乱光および蛍光などの強度の測定値、あるいはインピーダンス、容量値（キャパシタンス）およびインダクタンスなどの測定値は、電気信号に変換され、測定データとして入出力インターフェイス２５０から出力される。
【００３３】
２．３Ｄ立体画像の表示
図３に、本発明に係る３Ｄデータ解析装置により表示される三次元分布図を模式的に示す。この三次元分布図は、出力部１４２に３Ｄ立体画像により表示され、ユーザによって立体的に視認され得るものである。
【００３４】
三次元分布図５は、ユーザにより選択された３種のパラメータを座標軸とした座標空間６内に微小粒子の特性分布を示す。三次元分布図５では、選択されたパラメータの測定値から算出される位置に各微小粒子に対応する図形７が描画される。
【００３５】
図では、３種のパラメータを、前方散乱光（FS-Lin：Ｘ軸）、側方散乱（SS-Lin：Ｙ軸）および第一の蛍光（FL1-Lin：Ｚ軸）の組み合わせとした場合を例示した。各座標軸にとるパラメータは、任意に選択される組み合わせとできる。例えば、Ｘ軸に第一の蛍光（FL1）、Ｙ軸に第二の蛍光（FL2）、Ｚ軸にインピーダンスをとることもできる。
【００３６】
三次元分布図５の座標空間６内には、微小粒子をゲーティングする際に分画領域を設定するために使用されるポインタ５０が表示されている。ポインタ５０も３Ｄ立体画像においてユーザにより立体的に視認され得るものである。ポインタ５０は、入力部１４１からのユーザの入力信号により、座標空間６内で任意に移動、変形、拡大または縮小される。
【００３７】
ポインタ５０の形状は、３Ｄ立体画像の立体観察時に所定の立体形状として視認されれば特に限定されないが、例えば、図３（Ａ）に示すような球体や、（Ｂ）に示すような立方体を採用できる。
【００３８】
三次元分布図５の３Ｄ立体表示は、二眼式ステレオ立体画像により行われる。図４に、本発明に係る３Ｄデータ解析装置により表示される二眼式ステレオ立体画像を模式的に示す。
【００３９】
データ処理部１２０は、ユーザによってパラメータが選択されると三次元分布図５を作成して、該分布図を左眼で見たときの画像（左眼用画像５Ｌ）と、右眼で見たときの画像（右眼用画像５Ｒ）を作成する。出力部１４２（ディスプレイ４３）は、左眼用画像５Ｌおよび右眼用画像５Ｒを同時に表示し、左眼用画像５Ｌを左眼だけに、右眼用画像５Ｒを右眼だけに分離提示する。
【００４０】
分離提示は、例えばメガネ方式の一つである時分割方式では、左眼用画像５Ｌと右眼用画像５Ｒとをわずかな時間差で交互に表示させ、これにシャッタメガネ８を同期させることによって行うことができる。この他、分離提示は、アナグリフ方式および偏光フィルタ方式などの他のメガネ方式、パララックスバリア方式およびレンティキュラ方式などの裸眼方式、ステレオスコープ方式およびヘッドマウント方式などのビューア方式を採用してもよい。
【００４１】
ディスプレイ４３は、左眼用画像５Ｌおよび右眼用画像５Ｒを分離提示することによって、三次元空間で三次元分布図を見ている際に目に映る映像を再現し、ユーザに分布図を立体視させる。
【００４２】
３．ゲーティング
微小粒子をゲーティングするための分画領域の設定は、三次元分布図中に表示されるポインタ５０をマウス４１あるいはキーボード４２などの入力デバイスを用いて移動、変形、拡大または縮小させることによって行われる。分画領域は、移動等されたポインタ５０の内部空間に設定される。図５に設定され得る分画領域の具体例を示す。
【００４３】
図５（Ａ）は、楕円の球体とされたポインタ５０の内部空間に設定された分画領域５１を示している。この分画領域５１は、図３（Ａ）に示したポインタ５０を、解析対象とする微小粒子に対応する図形７が存在する領域に移動し、正円から楕円の球体に変形させて拡大し、当該領域とポインタ５０の内部空間とを全部一致あるいは一部一致させることによって設定できる。分画領域５１は複数設定することができ、解析対象とする２以上の微小粒子小集団を各分画領域内に囲い込むことができる。
【００４４】
また、図５（Ｂ）は、移動されるポインタ５０の内部空間の軌跡にそって連続的に設定された分画領域５１を示している。この分画領域５１は、図３（Ａ）に示したポインタ５０を、その内部空間が、解析対象とする微小粒子に対応する図形７が存在する領域を通過するように移動させ、当該領域とポインタ５０の内部空間の軌跡とが全部一致あるいは一部一致するようにすることで設定できる。分画領域５１は不連続な２以上の軌跡により複数設定することができ、解析対象とする２以上の微小粒子小集団を各分画領域内に囲い込むことができる。
【００４５】
ポインタ５０は任意の形状に変形され、任意の大きさ（体積）に拡大あるいは縮小されて三次元分布図中を移動され、解析対象とする微小粒子に対応する図形７をポインタ５０の内部空間に設定される分画領域５１内に囲い込む。この際、分画領域内に囲い込まれた図形７は、分画領域５１外に存在する図形７と異なる色（色相、彩度、明度）によって３Ｄ立体画像中に表示されるようにしてもよい。図５では、分画領域５１内にゲートされた図形７を黒色で、分画領域外の図形７を白色で示した。
【００４６】
分画領域内に囲い込まれた図形７と分画領域５１外に存在する図形７は、図６に示すように、異なる形状によって３Ｄ立体画像中に表示されるようにしてもよい。図６では、分画領域５１内にゲートされた図形７を８面体（図７（Ｂ）参照）で、分画領域外の図形７を６面体（図７（Ａ）参照）で示した。
【００４７】
分画領域内に囲い込まれた図形７と分画領域５１外に存在する図形７は、色および形状の双方が異なるように表示されてもよい。例えば、分画領域内に囲い込まれた図形７に対応する微小粒子が白血球である場合は、図形を白血球の形態を模した色および形とすることができる。これに対して、例えば、分画領域５１外に存在する図形７が赤血球である場合には、図形を赤血球の形態を模した色および形とすることができる。
【００４８】
さらに、ゲートされた以外の図形７（分画領域５１外に存在する図形７）は、３Ｄ立体画像中に表示しないようにすることもできる。あるいは、逆に、解析対象としない微小粒子に対応する図形７をゲーティングする場合には、ゲートされた図形７（分画領域内に存在する図形７）を３Ｄ立体画像中に表示しないようにすることもできる。
【００４９】
このように、３Ｄデータ解析装置１では、ユーザは、任意に選択した３種のパラメータを座標軸とした三次元分布図とその座標空間に表示されるポインタを立体視しながら、解析対象とする微小粒子ならびに微小粒子小集団のゲーティングを行うことができる。このため、３Ｄデータ解析装置１では、従来のように多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、分布図上において解析対象とする微小粒子ならびに微小粒子小集団を容易かつ直感的に特定してゲーティングすることができる。また、座標軸に用いるパラメータを任意に組み合わせて三次元分布図を表示させることで、一つのグラフで微小粒子の３つの特性に関する情報を得ることができる。このため、３Ｄデータ解析装置１では、解析対象とする微小粒子ならびに微小粒子小集団を高精度にゲートすることができ、かつ、従来のヒストグラムあるいはサイトグラムによる表示に比べて参照すべきグラフ数を減らしてゲーティング操作を効率化できる。
【００５０】
４．３Ｄ立体画像の特徴
以下、本発明に係る３Ｄデータ解析装置により表示される３Ｄ立体画像の特徴について順に説明する。
【００５１】
（４−１）図形の形状
図３中符号７で示した微小粒子に対応する図形は、所定形状のポリゴンの組み合わせからなる多面体として計算され、３Ｄ立体画像中に表示される。上記のように、データ処理部１２０は、ユーザが選択したパラメータの測定値に基づいて各微小粒子の座標空間内の位置と図形７を計算し、三次元分布図を作成する。この際、図形７を所定形状のポリゴンの組み合わせからなる多面体として計算することで、データ処理部１２０における計算負荷を軽減できる。また、図形７を所定形状のポリゴンの組み合わせからなる多面体として３Ｄ立体画像中に表示することで、立体視した際の画像の立体感を増強することができる。
【００５２】
所定形状のポリゴンの組み合わせからなる多面体としては、例えば、図７（Ａ）に示すような、三角形のポリゴンが６つ組み合わされてなる６面体や、（Ｂ）に示す該ポリゴンが８つ組み合わされてなる８面体を採用できる。図形７の形状は、分画領域内外で異なる形状とされ得るものであり、所定形状のポリゴンの組み合わせからなる多面体である限りにおいて特に限定されないが、計算負荷の軽減と立体感の観点から６面体あるいは８面体が好ましい。
【００５３】
（４−２）図形のシェード処理
３Ｄ立体画像において、図形７は、立体視した際に手前に観察される図形ほど濃く、奥に観察される図形ほど淡く表示される。このように図形７の濃淡を変化させる処理を、以下「シェード処理」と称するものとする。
【００５４】
シェード処理された図形７の立体観察像（以下単に「立体像」という）の概念図を図８に示す。図中矢印の方向に従って、手前に観察される図形７ほど濃く、奥に観察される図形７ほど淡くなっている。このように、図形７のシェード処理を行うことによって、３Ｄ立体画像の立体像に奥行きを与えて立体感を向上できる。
【００５５】
図９を参照して、シェード処理の処理方法を説明する。ディスプレイ４３には左眼用画像と右眼用画像が同時に表示され、立体視した際にディスプレイ４３の画面の位置に観察される図形７０の左眼用画像と右眼用画像は重なった状態で表示される（図９（Ｂ）参照）。
【００５６】
ディスプレイ４３に表示される左眼用画像が右眼用画像よりも右側にある場合（図９（Ａ）参照）、図形は、ディスプレイ４３の画面の位置よりも手前に立体視される。画面位置から飛び出して観察される図形の立体像を図中符号７１で示し、ディスプレイ４３に表示される図形７１の左眼用画像を符号７１Ｌで、右眼用画像を符号７１Ｒで示す。一方、ディスプレイ４３に表示される左眼用画像が右眼用画像よりも左側にある場合（図９（Ｃ）参照）、図形は、ディスプレイ４３の画面の位置よりも奥に立体視される。画面位置から飛び出して観察される図形の立体像を図中符号７２で示し、ディスプレイ４３に表示される図形７１の左眼用画像を符号７２Ｌで、右眼用画像を符号７２Ｒで示す。
【００５７】
シェード処理では、手前に観察される図形７１の左眼用画像７１Ｌおよび右眼用画像７１Ｒをより濃く、奥に観察される図形７２の左眼用画像７２Ｌおよび右眼用画像７２Ｒをより淡く表示するようにする。
【００５８】
（４−３）座標軸
３Ｄ立体画像において、座標軸は、立体視した際に手前に観察される部分ほど太く、奥に観察される分ほど細く表示される。太さを変化させた座標軸の立体像の概念図を図１０に示す。このように、座標軸の太さを変化させることによって、３Ｄ立体画像の立体像に奥行きを与えて立体感を向上できる。
【００５９】
また、図１０に示すように、座標軸の目盛間隔を立体視した際に手前に観察される部分ほど広く、奥に観察される分ほど狭く表示することによって、立体像にさらに奥行きを与えられる。また、座標軸名（図ではSS-Lin）および目盛数値（図では200, 400, 600, 800, 1000）の文字を手前ほど大きく、奥ほど小さく表示したりすることによっても同様の効果が得られる。なお、座標軸の太さ、目盛間隔および文字の大きさを変化させる処理は、上述のシェード処理を応用して行うことが可能である。
【００６０】
（４−４）動画像
既に説明したように、出力部１４２（ディスプレイ４３）への３Ｄ立体画像の表示は、入力部１４１（マウス４１あるいはキーボード４２など）からのユーザの入力信号に基づき任意に回転、拡大あるいは縮小させて行ってもよい。３Ｄ立体画像を回転する際には、図３に示したように、座標空間６を構成する立体形状（図では立方体）の各辺に座標軸が表示されていることが好ましい。これらの座標軸により座標空間６の立体形状が明確になるため、３Ｄ立体画像を回転させた際の三次元分布図の向きの変化を認識し易くなる。
【００６１】
ディスプレイ４３に表示される３Ｄ立体画像は、ユーザからの入力により随意に回転される他、常時一定方向にゆっくりと回転していてもよい。３Ｄ立体画像を常時回転する動画像により表示することで、静止画像による表示に比べて立体感を高められる。
【００６２】
さらに、ディスプレイ４３に表示される３Ｄ立体画像は、微小粒子に対応する図形が揺動する動画像により表示されてもよい。この際、立体視した際に奥に観察される図形に比して手前に観察される図形を大きく揺動させて表示させる。揺動動作が付与された図形の立体像の概念図を図１０に示す。図形７１，７２は、図中矢印に示すように左右に揺動し、左右への揺動幅は手前に観察される図形７１で大きく、奥に観察される図形７２で小さくなっている。このように、立体視した際に奥に観察される図形に比して手前に観察される図形をより大きく揺動させて表示することで、３Ｄ立体画像の立体像に奥行きを与えて立体感を向上できる。
【００６３】
以上のように、本発明に係る３Ｄデータ解析装置は、表示される３Ｄ立体画像の立体感を向上させるための工夫がなされている。従って、点（微小粒子に対応する図形）と線（座標軸）のみからなる三次元分布図であっても、ユーザはその立体像を良好に視認しながら測定データの解析を行うことができ、分布図上において解析対象とする微小粒子ならびに微小粒子小集団を容易かつ直感的に特定できる。
【産業上の利用可能性】
【００６４】
本発明に係る３Ｄデータ解析装置によれば、多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、容易かつ直感的なゲーティング操作が可能で、高精度に解析を行うことができる。従って、本発明に係る３Ｄデータ解析装置は、例えばフローサイトメータとともに用いられて、医療分野、公衆衛生分野あるいは創薬分野等において、細胞や微生物の特性を容易かつ高精度に解析するために用いられ得る。
【符号の説明】
【００６５】
１：３Ｄデータ解析装置、１０：中央処理装置、１１０：制御部、１２０：データ処理部、１３０：データ格納部、１４１：入力部、１４２：表示部、１５０：入出力インターフェイス、２：フローサイトメータ、２１０：制御部、２２０：フロー系、２３０：検出系、
２３１：光学検出部、２３２：電気検出部、２４０：分取系、２５０：入出力インターフェイス、３：微小粒子解析システム、３０：ハードディスク、３１：３Ｄデータ解析プログラム、３２：測定データ、３３：オペレーティングシステム、４：通信ケーブル、４１：マウス、４２：キーボード、４３：ディスプレイ、４４：プリンタ、５：三次元分布図、５０：ポインタ、５１：分画領域、６：座標空間、７：図形、８：シャッタメガネ
【技術分野】
【０００１】
本発明は、３Ｄデータ解析装置および３Ｄデータ解析方法に関する。より詳しくは、解析対象とする微小粒子の指定を３Ｄ立体画像中で行うことができる３Ｄデータ解析装置等に関する。
【背景技術】
【０００２】
細胞、微生物およびリポソーム等の生体関連粒子や、ラテックス粒子、ゲル粒子および工業用粒子等の合成粒子などの微小粒子を解析するため、微小粒子の分散液を流路内に導入して、微小粒子を光学的、電気的あるいは磁気的に測定する微小粒子測定装置が用いられている。
【０００３】
一例として、合成粒子を大きさや形状に応じて判別するパーティクルアナライザーがある。パーティクルアナライザーで測定可能なパラメータ（変数）には、微小粒子の元素組成および粒径などがある。
【０００４】
また、生体関連粒子の解析では、フローサイトメータ（フローサイトメトリー）が用いられている。フローサイトメータで測定可能なパラメータには、微小粒子の前方散乱光（FS）、側方散乱（SS）、蛍光（FL）およびインピーダンスなどがある。前方散乱光（FS）、側方散乱（SS）および蛍光（FL）は細胞や微生物（以下、単に「細胞」という）の光学的特性を示すパラメータとして、インピーダンスは細胞の電気的特性を示すパラメータとして用いられる。
【０００５】
具体的には、まず、前方散乱光は、レーザー光の軸に対して前方向の小さい角度で散乱する光であり、細胞の表面で生じるレーザー光の散乱光や回折光、屈折光からなる。前方散乱光は、主として細胞の大きさを示すパラメータとして用いられている。次に、側方散乱は、レーザー光の軸に対して約９０度の角度で散乱する光であり、細胞の内部の顆粒や核などで生じるレーザー光の散乱光である。側方散乱は、主として細胞の内部構造を示すパラメータとして用いられている。また、蛍光は、細胞に標識された蛍光色素から発生する光であり、蛍光色素標識抗体が認識する細胞表面抗原の有無や蛍光色素が結合した核酸の量などを示すパラメータとして用いられている。また、インピーダンスは、電気抵抗法によって測定され、細胞の容積を示すパラメータとして用いられている。
【０００６】
フローサイトメータにおける測定データの解析のため、これらの測定パラメータを軸にとって各細胞の測定値をプロットすることにより、細胞集団中における細胞の特性分布を表す図を作成して表示するデータ解析装置が用いられている。測定パラメータを１つとした一次元分布図はヒストグラムとも称され、測定パラメータをＸ軸にとり、細胞数（カウント）をＹ軸にとったグラフとして作成される。また、測定パラメータを２つとした二次元分布図はサイトグラムとも称され、一方の測定パラメータをＸ軸とし、他方の測定パラメータをＹ軸とした座標面内に、各細胞をその測定値に基づいてプロットしていくことにより作成される。
【０００７】
サンプルとする細胞集団の中には、解析の対象としない不要な細胞も含まれているため、測定データの解析は、サンプルとする細胞集団の中から解析対象とする細胞小集団を選び出した後に行われている。解析対象とする細胞小集団の選択は、ヒストグラム上あるいはサイトグラム上において該細胞小集団が存在する領域を指定することによって行われている。この操作は、ヒストグラム上あるいはサイトグラム上で指定する領域内に目的とする細胞を囲い込むことから、「ゲーティング」と呼ばれている。
【０００８】
一つの測定パラメータを軸としたヒストグラム上あるいは一つの組み合わせの測定パラメータを軸としたサイトグラム上においては、解析対象とする細胞小集団と不要な細胞とが重なり合う領域に存在している場合がある。例えば、ヒト抹消血をサンプルとしてリンパ球の解析を行う際、前方散乱光（FS）と側方散乱（SS）を軸にとったサイトグラム上では、一部の単球がリンパ球と同一の領域に存在する場合がある。従って、ゲーティングを行う際、ユーザは、単球が囲い込まれないようにしてリンパ球のみが存在する領域を指定してやる必要がある。
【０００９】
不要な細胞を囲い込むことなく、解析対象とする細胞小集団のみを囲い込むことができるように領域を指定するため、従来、ユーザは複数のヒストグラムあるいはサイトグラムを参照しながらゲーティングを行う必要があった。フローサイトメータの性能向上に伴って測定可能なパラメータ数も増加しており、ユーザはより多くのヒストグラムあるいはサイトグラムを参照する必要が生じていた。また、この際、ユーザは、２つのサイトグラムを組み合わせた立体的な分布図（三次元分布図）を想起しながらゲーティング操作を行うことが求められていた。
【００１０】
ユーザのゲーティング操作を支援するため、特許文献１には、「分析物から第１、第２および第３の測定データを取得する測定データ取得手段と、前記第１、第２および第３の測定データを軸とし、分析物に含まれる有形成分の分布を示す３次元分布図を作成する３次元分布図作成手段と、前記３次元分布図上に分画領域を変更可能に設定する領域設定手段と、前記領域設定手段によって設定された前記分画領域に属する有形成分について、前記第１および第２の測定データを軸とする２次元分布図および前記第１の測定データを軸とする度数分布図の少なくとも一方を作成する参考用分布図作成手段と、を備える分析装置。」が提案されている（当該文献請求項９参照）。この分析装置によれば、３次元分布図とともに表示される２次元分布図（サイトグラム）および度数分布図（ヒストグラム）を参照しながら、３次元分布図上に分画領域を設定することができる。なお、この分析装置の３次元分布図は、ディスプレイ上に平面的に表示されるものであり、立体視されるものではない。
【００１１】
本発明に関連して、二眼式ステレオ立体画像技術（３Ｄ立体画像技術）について説明する。二眼式ステレオ立体映像では、まず、対象物を右眼および左眼で見たときの２つの映像を用意する。そして、これらの映像を同時に表示した上で、右眼用映像を右眼だけに、左眼用映像を左眼だけに提示する。これにより、三次元空間で対象物を見ている際に目に映る映像を再現し、ユーザに対象物を立体視させる。
【００１２】
立体視が可能な３Ｄディスプレイには、（ａ）メガネ方式、（ｂ）裸眼方式、（ｃ）ビューア方式が主に採用されている。（ａ）メガネ方式には、さらにアナグリフ方式、偏光フィルタ方式、時分割方式がある。また、（ｂ）裸眼方式にはパララックスバリア方式、レンティキュラ方式があり、（ｃ）ビューア方式にはステレオスコープ方式、ヘッドマウント方式がある。
【先行技術文献】
【特許文献】
【００１３】
【特許文献１】特開２００６−１７４９７号公報
【発明の概要】
【発明が解決しようとする課題】
【００１４】
上述のように、従来、フローサイトメータに用いられているデータ解析装置では、ゲーティングの際、不要な細胞を囲い込むことなく、解析対象とする細胞小集団のみを囲い込むため、ユーザは多数のヒストグラムあるいはサイトグラムを参照したり、２つのサイトグラムを組み合わせた立体的な分布図（三次元分布図）を想起したりしながら操作を行う必要があった。このため、従来のデータ解析装置では、ゲーティングの操作に時間がかかり、ユーザが操作に習熟していない場合には不要な細胞が混入してしまい、解析精度に問題が生じていた。
【００１５】
そこで、本発明は、多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、容易かつ直感的なゲーティング操作が可能で、高精度に解析を行うことができるデータ解析装置を提供することを主な目的とする。
【課題を解決するための手段】
【００１６】
上記課題解決のため、本発明は、微小粒子の測定データを保存するデータ格納部と、前記測定データから独立した３種の変数を選択する入力部と、前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成するデータ処理部と、前記３Ｄ立体画像を表示する表示部と、を有し、前記入力部からの入力信号に基づいて前記座標空間内に分画領域を変更可能に設定し、前記３Ｄ立体画像中に表示する３Ｄデータ解析装置を提供する。
この３Ｄデータ解析装置によれば、任意に選択した３種のパラメータを座標軸とした三次元分布図を立体視しながら、その座標空間において解析対象とする微小粒子ならびに微小粒子小集団のゲーティングを行うことができる。
【００１７】
また、本発明は、微小粒子の測定データから独立した３種の変数を選択する手順と、前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成する手順と、前記３Ｄ立体画像を表示させる手順と、前記３Ｄ立体画像中の前記座標空間内に分画領域を設定する手順と、を含む３Ｄデータ解析方法をも提供する。
【００１８】
本発明において、「微小粒子」には、細胞、微生物およびリポソーム等の生体関連粒子、ラテックス粒子、ゲル粒子、工業用粒子等の合成粒子などの微小粒子が広く含まれるものとする。
細胞には、動物細胞（血球系細胞など）および植物細胞が含まれる。微生物には、大腸菌等の細菌類、タバコモザイクウイルス等のウイルス類、イースト菌等の菌類などが含まれる。生体関連粒子には、各種細胞を構成する染色体、リポソーム、ミトコンドリア、オルガネラ(細胞小器官)なども含まれる。さらに、生体関連粒子には、核酸やタンパク質、これらの複合体などの生体関連高分子も包含され得るものとする。工業用粒子は、有機もしくは無機高分子材料、金属などであってもよい。有機高分子材料には、ポリスチレン、スチレン・ジビニルベンゼンおよびポリメチルメタクリレートなどが含まれる。無機高分子材料には、ガラス、シリカおよび磁性体材料などが含まれる。金属には、金コロイドおよびアルミなどが含まれる。これら微小粒子の形状は、一般には球形であるが、非球形であってもよい。また、微小粒子の大きさや質量なども特に限定されない。
【発明の効果】
【００１９】
本発明により、多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、容易かつ直感的なゲーティング操作が可能で、高精度に解析を行うことができるデータ解析装置が提供される。
【図面の簡単な説明】
【００２０】
【図１】フローサイトメータに連設された本発明に係る３Ｄデータ解析装置の構成を説明するブロック図である。
【図２】本発明に係る３Ｄデータ解析装置の機能的構成を説明するブロック図である。
【図３】本発明に係る３Ｄデータ解析装置により表示される三次元分布図を説明する模式図である。
【図４】本発明に係る３Ｄデータ解析装置により表示される二眼式ステレオ立体画像（３Ｄ立体画像）を説明する模式図である。
【図５】本発明に係る３Ｄデータ解析装置において設定され得る分画領域の具体例を説明する模式図である。
【図６】本発明に係る３Ｄデータ解析装置において設定され得る分画領域の具体例を説明する模式図である。
【図７】３Ｄ立体画像中の微小粒子に対応する図形の形状を説明する模式図である。
【図８】シェード処理された図形の立体観察像を説明する概念図である。
【図９】シェード処理の処理方法を説明する模式図である。
【図１０】座標軸の立体観察像を説明する概念図である。
【図１１】微小粒子に対応する図形を揺動させた動画像の立体観察像を説明する概念図である。
【発明を実施するための形態】
【００２１】
以下、本発明を実施するための好適な形態について図面を参照しながら説明する。なお、以下に説明する実施形態は、本発明の代表的な実施形態の一例を示したものであり、これにより本発明の範囲が狭く解釈されることはない。なお、説明は以下の順序で行う。
１．３Ｄデータ解析装置の構成
２．３Ｄ立体画像の表示
３．ゲーティング
４．３Ｄ立体画像の特徴
（４−１）図形の形状
（４−２）図形のシェード処理
（４−３）座標軸
（４−４）動画像
【００２２】
１．３Ｄデータ解析装置の構成
図１に、本発明に係る３Ｄデータ解析装置の装置構成を示す。ここでは、３Ｄデータ解析装置を微小粒子測定装置に連設して微小粒子解析システムとして構成した実施形態を示す。また、図２に、この微小粒子解析システムの機能的構成を示す。以下では、微小粒子測定装置としてフローサイトメータを用いる場合を例に説明する。
【００２３】
図中符号１で示す３Ｄデータ解析装置は、通信ケーブル４によりフローサイトメータ２と接続されて微小粒子解析システム３を構成している。３Ｄデータ解析装置１は、中央処理装置（ＣＰＵ）１０、メモリ２０、ハードディスク３０、ユーザインターフェイスなどを含んでいる。ハードディスク３０には、３Ｄデータ解析プログラム３１、微小粒子の測定データ３２およびオペレーティングシステム（ＯＳ）３３などが格納・保持されている。ユーザインターフェイスには、ユーザからの情報の入力を受け付けるマウス４１およびキーボード４２などの入力デバイスと、ユーザに対して情報を出力するディスプレイ４３およびプリンタ４４などの出力デバイスが含まれる。なお、マウス４１およびキーボード４２に替えて、あるいはこれらとともに、スティックコントローラおよびペンタブレットなどの入力デバイスを設けてもよい。
【００２４】
データ格納部１３０（ハードディスク３０）は、フローサイトメータ２から出力される微小粒子（細胞）の測定データ３２を保存する。フローサイトメータ２の入出力インターフェイス２５０から出力される測定データは、通信ケーブル４を介して３Ｄデータ解析装置１の入出力インターフェイス１５０に入力され、データ格納部３０（ハードディスク３０）に格納される。
【００２５】
測定データ３２は、データ処理部１２０において処理される。データ処理部１２０は、入力部１４１（マウス４１あるいはキーボード４２など）からのユーザの入力を受けて処理を開始する。すなわち、データ処理部１２０は、ユーザによって測定データ３２から独立した３種の変数（パラメータ）が選択され入力されると、選択されたパラメータを座標軸として微小粒子の特性分布を表す三次元分布図を作成する。三次元分布図は、選択されたパラメータを座標軸とした座標空間内に微小粒子をプロットすることにより作成される。微小粒子のプロットは、選択されたパラメータの測定値から各微小粒子の座標空間内の位置と図形を計算し、算出された位置に算出された図形を描画することによって行われる。
【００２６】
また、データ処理部１２０は、作成した三次元分布図の座標空間内に、入力部１４１からのユーザの入力信号に基づいて微小粒子をゲーティングするための分画領域を設定する。入力部１４１からのユーザの入力は、三次元分布図中に表示されるポインタをマウス４１あるいはキーボード４２などの入力デバイスを用いて移動、変形、拡大または縮小させることによって行われる。
【００２７】
さらに、データ処理部１２０は、設定された分画領域内に存在する微小粒子の測定データの処理（クラスター解析）を行う。クラスター解析は、例えば、全微小粒子に占める分画領域内に存在する微小粒子の割合、画領域内に存在する微小粒子の数、平均値、標準偏差および変動係数などの統計データが算出される。
【００２８】
ここで、「独立したパラメータ」とは、微小粒子の前方散乱光（FS）、側方散乱（SS）、蛍光（FL）およびインピーダンスなどから選択される互いに異なるパラメータを意味する。蛍光（FL）は、微小粒子に標識された蛍光色素の波長域毎に異なるパラメータとして扱うことができ、ＦＬ１，ＦＬ２〜ＦＬｎ（ｎは３以上の整数）などと表記される。独立した３種のパラメータとしては、例えば、前方散乱光（FS）、側方散乱（SS）および蛍光（FL1）の組み合わせや、前方散乱光（FS）、側方散乱（SS）およびインピーダンスの組み合わせが例示される。この他、独立した３種のパラメータは、測定データから任意に選択して組み合わせることができる。
【００２９】
データ処理部１２０により作成された三次元分布図と、その座標空間内に設定された分画領域とは、出力部１４２（ディスプレイ４３）に３Ｄ立体画像として表示される。この３Ｄ立体画像は、次に詳しく説明する二眼式ステレオ立体画像とされる。出力部１４２への３Ｄ立体画像の表示は、入力部１４１（マウス４１あるいはキーボード４２など）からのユーザの入力信号に基づき任意に回転、拡大あるいは縮小させて行ってもよい。また、データ処理部１２０によるクラスター解析の結果も、出力部１４２（ディスプレイ４３およびプリンタ４４など）に出力される。統計データは、出力部１４２に数値あるいは図表によって表示、印刷等されてよい。
【００３０】
フローサイトメータ２は、従来公知の装置と同様の構成あるいはこれを適宜改変した構成とでき、具体的には、制御部２１０、フロー系２２０、検出系２３０、入出力インターフェイス２５０等から構成される。
【００３１】
フロー系２２０は、フローセルやマイクロチップに形成された流路内において、微小粒子を含むサンプル液層流をシース液層流の中心に流し、層流中に微小粒子を一列に配列させる。検出系２３０は、流路を通流する微小粒子の特性を示すパラメータ値を取得する。具体的には、光学検出部２３１は、通流する微小粒子に光を照射し、微小粒子から生じる散乱光や蛍光などを検出し、その強度を取得する。光学検出部２３１は、レーザー光源、レンズ、ミラー、フィルター、ＣＣＤおよびＣＭＯＳ素子等のエリア撮像素子あるいはＰＭＴ（photo multiplier tube）などを含んでなる。また、電気検出部２３２は、通流する微小粒子に対向して配された電極を含んでなり、微小粒子のインピーダンス、容量値（キャパシタンス）およびインダクタンスなどを取得する。フローサイトメータ２は、解析の結果所望の特性を有すると判定された微小粒子を分取するための分取系２４０を備えていてもよい。分取系２４０には、例えば、微小粒子を含む液滴をフローセル外の空間に吐出し、液滴の移動方向を制御することにより、所望の微小粒子のみを容器に回収する方式を採用できる。
【００３２】
検出系２３０において検出された散乱光および蛍光などの強度の測定値、あるいはインピーダンス、容量値（キャパシタンス）およびインダクタンスなどの測定値は、電気信号に変換され、測定データとして入出力インターフェイス２５０から出力される。
【００３３】
２．３Ｄ立体画像の表示
図３に、本発明に係る３Ｄデータ解析装置により表示される三次元分布図を模式的に示す。この三次元分布図は、出力部１４２に３Ｄ立体画像により表示され、ユーザによって立体的に視認され得るものである。
【００３４】
三次元分布図５は、ユーザにより選択された３種のパラメータを座標軸とした座標空間６内に微小粒子の特性分布を示す。三次元分布図５では、選択されたパラメータの測定値から算出される位置に各微小粒子に対応する図形７が描画される。
【００３５】
図では、３種のパラメータを、前方散乱光（FS-Lin：Ｘ軸）、側方散乱（SS-Lin：Ｙ軸）および第一の蛍光（FL1-Lin：Ｚ軸）の組み合わせとした場合を例示した。各座標軸にとるパラメータは、任意に選択される組み合わせとできる。例えば、Ｘ軸に第一の蛍光（FL1）、Ｙ軸に第二の蛍光（FL2）、Ｚ軸にインピーダンスをとることもできる。
【００３６】
三次元分布図５の座標空間６内には、微小粒子をゲーティングする際に分画領域を設定するために使用されるポインタ５０が表示されている。ポインタ５０も３Ｄ立体画像においてユーザにより立体的に視認され得るものである。ポインタ５０は、入力部１４１からのユーザの入力信号により、座標空間６内で任意に移動、変形、拡大または縮小される。
【００３７】
ポインタ５０の形状は、３Ｄ立体画像の立体観察時に所定の立体形状として視認されれば特に限定されないが、例えば、図３（Ａ）に示すような球体や、（Ｂ）に示すような立方体を採用できる。
【００３８】
三次元分布図５の３Ｄ立体表示は、二眼式ステレオ立体画像により行われる。図４に、本発明に係る３Ｄデータ解析装置により表示される二眼式ステレオ立体画像を模式的に示す。
【００３９】
データ処理部１２０は、ユーザによってパラメータが選択されると三次元分布図５を作成して、該分布図を左眼で見たときの画像（左眼用画像５Ｌ）と、右眼で見たときの画像（右眼用画像５Ｒ）を作成する。出力部１４２（ディスプレイ４３）は、左眼用画像５Ｌおよび右眼用画像５Ｒを同時に表示し、左眼用画像５Ｌを左眼だけに、右眼用画像５Ｒを右眼だけに分離提示する。
【００４０】
分離提示は、例えばメガネ方式の一つである時分割方式では、左眼用画像５Ｌと右眼用画像５Ｒとをわずかな時間差で交互に表示させ、これにシャッタメガネ８を同期させることによって行うことができる。この他、分離提示は、アナグリフ方式および偏光フィルタ方式などの他のメガネ方式、パララックスバリア方式およびレンティキュラ方式などの裸眼方式、ステレオスコープ方式およびヘッドマウント方式などのビューア方式を採用してもよい。
【００４１】
ディスプレイ４３は、左眼用画像５Ｌおよび右眼用画像５Ｒを分離提示することによって、三次元空間で三次元分布図を見ている際に目に映る映像を再現し、ユーザに分布図を立体視させる。
【００４２】
３．ゲーティング
微小粒子をゲーティングするための分画領域の設定は、三次元分布図中に表示されるポインタ５０をマウス４１あるいはキーボード４２などの入力デバイスを用いて移動、変形、拡大または縮小させることによって行われる。分画領域は、移動等されたポインタ５０の内部空間に設定される。図５に設定され得る分画領域の具体例を示す。
【００４３】
図５（Ａ）は、楕円の球体とされたポインタ５０の内部空間に設定された分画領域５１を示している。この分画領域５１は、図３（Ａ）に示したポインタ５０を、解析対象とする微小粒子に対応する図形７が存在する領域に移動し、正円から楕円の球体に変形させて拡大し、当該領域とポインタ５０の内部空間とを全部一致あるいは一部一致させることによって設定できる。分画領域５１は複数設定することができ、解析対象とする２以上の微小粒子小集団を各分画領域内に囲い込むことができる。
【００４４】
また、図５（Ｂ）は、移動されるポインタ５０の内部空間の軌跡にそって連続的に設定された分画領域５１を示している。この分画領域５１は、図３（Ａ）に示したポインタ５０を、その内部空間が、解析対象とする微小粒子に対応する図形７が存在する領域を通過するように移動させ、当該領域とポインタ５０の内部空間の軌跡とが全部一致あるいは一部一致するようにすることで設定できる。分画領域５１は不連続な２以上の軌跡により複数設定することができ、解析対象とする２以上の微小粒子小集団を各分画領域内に囲い込むことができる。
【００４５】
ポインタ５０は任意の形状に変形され、任意の大きさ（体積）に拡大あるいは縮小されて三次元分布図中を移動され、解析対象とする微小粒子に対応する図形７をポインタ５０の内部空間に設定される分画領域５１内に囲い込む。この際、分画領域内に囲い込まれた図形７は、分画領域５１外に存在する図形７と異なる色（色相、彩度、明度）によって３Ｄ立体画像中に表示されるようにしてもよい。図５では、分画領域５１内にゲートされた図形７を黒色で、分画領域外の図形７を白色で示した。
【００４６】
分画領域内に囲い込まれた図形７と分画領域５１外に存在する図形７は、図６に示すように、異なる形状によって３Ｄ立体画像中に表示されるようにしてもよい。図６では、分画領域５１内にゲートされた図形７を８面体（図７（Ｂ）参照）で、分画領域外の図形７を６面体（図７（Ａ）参照）で示した。
【００４７】
分画領域内に囲い込まれた図形７と分画領域５１外に存在する図形７は、色および形状の双方が異なるように表示されてもよい。例えば、分画領域内に囲い込まれた図形７に対応する微小粒子が白血球である場合は、図形を白血球の形態を模した色および形とすることができる。これに対して、例えば、分画領域５１外に存在する図形７が赤血球である場合には、図形を赤血球の形態を模した色および形とすることができる。
【００４８】
さらに、ゲートされた以外の図形７（分画領域５１外に存在する図形７）は、３Ｄ立体画像中に表示しないようにすることもできる。あるいは、逆に、解析対象としない微小粒子に対応する図形７をゲーティングする場合には、ゲートされた図形７（分画領域内に存在する図形７）を３Ｄ立体画像中に表示しないようにすることもできる。
【００４９】
このように、３Ｄデータ解析装置１では、ユーザは、任意に選択した３種のパラメータを座標軸とした三次元分布図とその座標空間に表示されるポインタを立体視しながら、解析対象とする微小粒子ならびに微小粒子小集団のゲーティングを行うことができる。このため、３Ｄデータ解析装置１では、従来のように多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、分布図上において解析対象とする微小粒子ならびに微小粒子小集団を容易かつ直感的に特定してゲーティングすることができる。また、座標軸に用いるパラメータを任意に組み合わせて三次元分布図を表示させることで、一つのグラフで微小粒子の３つの特性に関する情報を得ることができる。このため、３Ｄデータ解析装置１では、解析対象とする微小粒子ならびに微小粒子小集団を高精度にゲートすることができ、かつ、従来のヒストグラムあるいはサイトグラムによる表示に比べて参照すべきグラフ数を減らしてゲーティング操作を効率化できる。
【００５０】
４．３Ｄ立体画像の特徴
以下、本発明に係る３Ｄデータ解析装置により表示される３Ｄ立体画像の特徴について順に説明する。
【００５１】
（４−１）図形の形状
図３中符号７で示した微小粒子に対応する図形は、所定形状のポリゴンの組み合わせからなる多面体として計算され、３Ｄ立体画像中に表示される。上記のように、データ処理部１２０は、ユーザが選択したパラメータの測定値に基づいて各微小粒子の座標空間内の位置と図形７を計算し、三次元分布図を作成する。この際、図形７を所定形状のポリゴンの組み合わせからなる多面体として計算することで、データ処理部１２０における計算負荷を軽減できる。また、図形７を所定形状のポリゴンの組み合わせからなる多面体として３Ｄ立体画像中に表示することで、立体視した際の画像の立体感を増強することができる。
【００５２】
所定形状のポリゴンの組み合わせからなる多面体としては、例えば、図７（Ａ）に示すような、三角形のポリゴンが６つ組み合わされてなる６面体や、（Ｂ）に示す該ポリゴンが８つ組み合わされてなる８面体を採用できる。図形７の形状は、分画領域内外で異なる形状とされ得るものであり、所定形状のポリゴンの組み合わせからなる多面体である限りにおいて特に限定されないが、計算負荷の軽減と立体感の観点から６面体あるいは８面体が好ましい。
【００５３】
（４−２）図形のシェード処理
３Ｄ立体画像において、図形７は、立体視した際に手前に観察される図形ほど濃く、奥に観察される図形ほど淡く表示される。このように図形７の濃淡を変化させる処理を、以下「シェード処理」と称するものとする。
【００５４】
シェード処理された図形７の立体観察像（以下単に「立体像」という）の概念図を図８に示す。図中矢印の方向に従って、手前に観察される図形７ほど濃く、奥に観察される図形７ほど淡くなっている。このように、図形７のシェード処理を行うことによって、３Ｄ立体画像の立体像に奥行きを与えて立体感を向上できる。
【００５５】
図９を参照して、シェード処理の処理方法を説明する。ディスプレイ４３には左眼用画像と右眼用画像が同時に表示され、立体視した際にディスプレイ４３の画面の位置に観察される図形７０の左眼用画像と右眼用画像は重なった状態で表示される（図９（Ｂ）参照）。
【００５６】
ディスプレイ４３に表示される左眼用画像が右眼用画像よりも右側にある場合（図９（Ａ）参照）、図形は、ディスプレイ４３の画面の位置よりも手前に立体視される。画面位置から飛び出して観察される図形の立体像を図中符号７１で示し、ディスプレイ４３に表示される図形７１の左眼用画像を符号７１Ｌで、右眼用画像を符号７１Ｒで示す。一方、ディスプレイ４３に表示される左眼用画像が右眼用画像よりも左側にある場合（図９（Ｃ）参照）、図形は、ディスプレイ４３の画面の位置よりも奥に立体視される。画面位置から飛び出して観察される図形の立体像を図中符号７２で示し、ディスプレイ４３に表示される図形７１の左眼用画像を符号７２Ｌで、右眼用画像を符号７２Ｒで示す。
【００５７】
シェード処理では、手前に観察される図形７１の左眼用画像７１Ｌおよび右眼用画像７１Ｒをより濃く、奥に観察される図形７２の左眼用画像７２Ｌおよび右眼用画像７２Ｒをより淡く表示するようにする。
【００５８】
（４−３）座標軸
３Ｄ立体画像において、座標軸は、立体視した際に手前に観察される部分ほど太く、奥に観察される分ほど細く表示される。太さを変化させた座標軸の立体像の概念図を図１０に示す。このように、座標軸の太さを変化させることによって、３Ｄ立体画像の立体像に奥行きを与えて立体感を向上できる。
【００５９】
また、図１０に示すように、座標軸の目盛間隔を立体視した際に手前に観察される部分ほど広く、奥に観察される分ほど狭く表示することによって、立体像にさらに奥行きを与えられる。また、座標軸名（図ではSS-Lin）および目盛数値（図では200, 400, 600, 800, 1000）の文字を手前ほど大きく、奥ほど小さく表示したりすることによっても同様の効果が得られる。なお、座標軸の太さ、目盛間隔および文字の大きさを変化させる処理は、上述のシェード処理を応用して行うことが可能である。
【００６０】
（４−４）動画像
既に説明したように、出力部１４２（ディスプレイ４３）への３Ｄ立体画像の表示は、入力部１４１（マウス４１あるいはキーボード４２など）からのユーザの入力信号に基づき任意に回転、拡大あるいは縮小させて行ってもよい。３Ｄ立体画像を回転する際には、図３に示したように、座標空間６を構成する立体形状（図では立方体）の各辺に座標軸が表示されていることが好ましい。これらの座標軸により座標空間６の立体形状が明確になるため、３Ｄ立体画像を回転させた際の三次元分布図の向きの変化を認識し易くなる。
【００６１】
ディスプレイ４３に表示される３Ｄ立体画像は、ユーザからの入力により随意に回転される他、常時一定方向にゆっくりと回転していてもよい。３Ｄ立体画像を常時回転する動画像により表示することで、静止画像による表示に比べて立体感を高められる。
【００６２】
さらに、ディスプレイ４３に表示される３Ｄ立体画像は、微小粒子に対応する図形が揺動する動画像により表示されてもよい。この際、立体視した際に奥に観察される図形に比して手前に観察される図形を大きく揺動させて表示させる。揺動動作が付与された図形の立体像の概念図を図１０に示す。図形７１，７２は、図中矢印に示すように左右に揺動し、左右への揺動幅は手前に観察される図形７１で大きく、奥に観察される図形７２で小さくなっている。このように、立体視した際に奥に観察される図形に比して手前に観察される図形をより大きく揺動させて表示することで、３Ｄ立体画像の立体像に奥行きを与えて立体感を向上できる。
【００６３】
以上のように、本発明に係る３Ｄデータ解析装置は、表示される３Ｄ立体画像の立体感を向上させるための工夫がなされている。従って、点（微小粒子に対応する図形）と線（座標軸）のみからなる三次元分布図であっても、ユーザはその立体像を良好に視認しながら測定データの解析を行うことができ、分布図上において解析対象とする微小粒子ならびに微小粒子小集団を容易かつ直感的に特定できる。
【産業上の利用可能性】
【００６４】
本発明に係る３Ｄデータ解析装置によれば、多数のヒストグラムあるいはサイトグラムを参照したり、三次元分布図を想起したりすることなく、容易かつ直感的なゲーティング操作が可能で、高精度に解析を行うことができる。従って、本発明に係る３Ｄデータ解析装置は、例えばフローサイトメータとともに用いられて、医療分野、公衆衛生分野あるいは創薬分野等において、細胞や微生物の特性を容易かつ高精度に解析するために用いられ得る。
【符号の説明】
【００６５】
１：３Ｄデータ解析装置、１０：中央処理装置、１１０：制御部、１２０：データ処理部、１３０：データ格納部、１４１：入力部、１４２：表示部、１５０：入出力インターフェイス、２：フローサイトメータ、２１０：制御部、２２０：フロー系、２３０：検出系、
２３１：光学検出部、２３２：電気検出部、２４０：分取系、２５０：入出力インターフェイス、３：微小粒子解析システム、３０：ハードディスク、３１：３Ｄデータ解析プログラム、３２：測定データ、３３：オペレーティングシステム、４：通信ケーブル、４１：マウス、４２：キーボード、４３：ディスプレイ、４４：プリンタ、５：三次元分布図、５０：ポインタ、５１：分画領域、６：座標空間、７：図形、８：シャッタメガネ
【特許請求の範囲】
【請求項１】
微小粒子の測定データを保存するデータ格納部と、
前記測定データから独立した３種の変数を選択する入力部と、
前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成するデータ処理部と、
前記３Ｄ立体画像を表示する表示部と、を有し、
前記入力部からの入力信号に基づいて前記座標空間内に分画領域を変更可能に設定し、前記３Ｄ立体画像中に表示する３Ｄデータ解析装置。
【請求項２】
前記分画領域は、前記３Ｄ立体画像の立体像観察時に所定の立体形状として視認されるポインタの内部空間に設定され、
前記入力部は、前記座標空間内で前記ポインタを移動、変形、拡大または縮小させる入力デバイスを含む請求項１記載の３Ｄデータ解析装置。
【請求項３】
前記分画領域は、前記座標空間内を移動される前記ポインタの軌跡に沿って連続的に設定される請求項２記載の３Ｄデータ解析装置。
【請求項４】
前記３Ｄ立体画像において、前記分画領域内に存在する前記図形を、前記分画領域外に存在する前記図形と異なる色で表示する請求項３記載の３Ｄデータ解析装置。
【請求項５】
前記３Ｄ立体画像において、前記分画領域内に存在する前記図形の形状を、前記分画領域外に存在する前記図形と異なる形状で表示する請求項３記載の３Ｄデータ解析装置。
【請求項６】
前記３Ｄ立体画像において、前記分画領域内に存在する前記図形を非表示とする請求項３記載の３Ｄデータ解析装置。
【請求項７】
前記表示手段がディスプレイであり、
該ディスプレイに表示される前記３Ｄ立体画像を立体視するためのメガネを備える請求項１〜６のいずれか一項に記載の３Ｄデータ解析装置。
【請求項８】
請求項１〜６に記載のいずれか一項に記載の３Ｄデータ解析装置と微小粒子測定装置とを連設してなる微小粒子解析システム。
【請求項９】
微小粒子の測定データから独立した３種の変数を選択する手順と、
前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成する手順と、
前記３Ｄ立体画像を表示させる手順と、
前記３Ｄ立体画像中の前記座標空間内に分画領域を設定する手順と、を含む３Ｄデータ解析方法。
【請求項１０】
前記３Ｄ立体画像の立体像観察時に所定の立体形状として視認されるポインタを前記座標空間内で移動、変形、拡大または縮小させ、前記ポインタの内部空間に前記分画領域を設定する請求項９記載の３Ｄデータ解析方法。
【請求項１】
微小粒子の測定データを保存するデータ格納部と、
前記測定データから独立した３種の変数を選択する入力部と、
前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成するデータ処理部と、
前記３Ｄ立体画像を表示する表示部と、を有し、
前記入力部からの入力信号に基づいて前記座標空間内に分画領域を変更可能に設定し、前記３Ｄ立体画像中に表示する３Ｄデータ解析装置。
【請求項２】
前記分画領域は、前記３Ｄ立体画像の立体像観察時に所定の立体形状として視認されるポインタの内部空間に設定され、
前記入力部は、前記座標空間内で前記ポインタを移動、変形、拡大または縮小させる入力デバイスを含む請求項１記載の３Ｄデータ解析装置。
【請求項３】
前記分画領域は、前記座標空間内を移動される前記ポインタの軌跡に沿って連続的に設定される請求項２記載の３Ｄデータ解析装置。
【請求項４】
前記３Ｄ立体画像において、前記分画領域内に存在する前記図形を、前記分画領域外に存在する前記図形と異なる色で表示する請求項３記載の３Ｄデータ解析装置。
【請求項５】
前記３Ｄ立体画像において、前記分画領域内に存在する前記図形の形状を、前記分画領域外に存在する前記図形と異なる形状で表示する請求項３記載の３Ｄデータ解析装置。
【請求項６】
前記３Ｄ立体画像において、前記分画領域内に存在する前記図形を非表示とする請求項３記載の３Ｄデータ解析装置。
【請求項７】
前記表示手段がディスプレイであり、
該ディスプレイに表示される前記３Ｄ立体画像を立体視するためのメガネを備える請求項１〜６のいずれか一項に記載の３Ｄデータ解析装置。
【請求項８】
請求項１〜６に記載のいずれか一項に記載の３Ｄデータ解析装置と微小粒子測定装置とを連設してなる微小粒子解析システム。
【請求項９】
微小粒子の測定データから独立した３種の変数を選択する手順と、
前記３種の変数を座標軸とする座標空間内における位置と図形を計算し、前記微小粒子の特性分布を表す３Ｄ立体画像を作成する手順と、
前記３Ｄ立体画像を表示させる手順と、
前記３Ｄ立体画像中の前記座標空間内に分画領域を設定する手順と、を含む３Ｄデータ解析方法。
【請求項１０】
前記３Ｄ立体画像の立体像観察時に所定の立体形状として視認されるポインタを前記座標空間内で移動、変形、拡大または縮小させ、前記ポインタの内部空間に前記分画領域を設定する請求項９記載の３Ｄデータ解析方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２０１２−１１７９９３（Ｐ２０１２−１１７９９３Ａ）
【公開日】平成２４年６月２１日（２０１２．６．２１）
【国際特許分類】
【出願番号】特願２０１０−２６９９４５（Ｐ２０１０−２６９９４５）
【出願日】平成２２年１２月３日（２０１０．１２．３）
【出願人】（０００００２１８５）ソニー株式会社 (34,172)
【Ｆターム（参考）】
【公開日】平成２４年６月２１日（２０１２．６．２１）
【国際特許分類】
【出願日】平成２２年１２月３日（２０１０．１２．３）
【出願人】（０００００２１８５）ソニー株式会社 (34,172)
【Ｆターム（参考）】
[ Back to top ]