"PS-OCTで得られた計測データを処理するコンピュータに搭載されるプログラムであって、コンピュータを、PS-OCTで計測したリターデーションについての位相差の分布のデータを、非線形の変換関数によって、対称な位相差の分布のデータに変換し、平均することにより、真の位相量を求める手段として機能させることを特徴とするプログラム。PS-OCTで得られた計測データを処理するコンピュータに搭載されるプログラムであって、コンピュータを、予め、0~πまでの複数の位相量それぞれについての分布から成る位相差の分布のデータのセットを、異なるESNR毎に複数のセット作成し記憶する手段、PS-OCTで計測した位相差の分布のデータを対称な位相差の分布のデータに変換する非線形の変換関数の係数を決める手段、PS-OCTで計測した位相差の分布のデータからESNRを見積もり、該見積もられたESNRに基づき、前記異なるESNR毎に作成した位相差の分布の複数のセットのいずれかを選択する手段、PS-OCTで計測した位相差の分布のデータを、前記選択された位相差の分布のセットに応じて決められた係数の下で、前記変換関数によって対称な位相差の分布のデータに変換する手段、及び前記対称な位相差の分布のデータを平均することにより、真の位相量を算出する手段、として機能させることを特徴とするプログラム。前記非線形の変換関数は、次の式2であることを特徴とする請求項2記載のプログラム。但し、式2中、δEは対称な位相差の分布であり、δMはPS-OCTで計測した位相量であり、bii=0、1、2、3、.........nは変換関数の係数であり、γはESNRである。前記変換関数の係数を決める手段は、前記非線形の変換関数について、0~πまで、π/mきざみで変え、次のm1の連立方程式である式3を作成し、該連立方程式を解いて、変換関数の係数biのテーブルを作る構成であることを特徴とする請求項3記載のプログラム。但し、式3中のmは、πをきざむ回数である。前記連立方程式を解いて前記変換関数の係数を決める手段においては、式4の行列DMの疑似逆行列DMを数値的に求めて、式5で計算し、最適の変換関数の係数biを一義的に決める構成であることを特徴とする請求項4記載のプログラム。但し、式4中のmは、πをきざむ回数である。前記変換関数の係数を決める手段は、変分法を用い、次の式6で示される二乗誤差が最小となるように最適の変換関数の係数biを一義的に決める構成であることを特徴とする請求項3記載のプログラム。但し、式6中のmは、πをきざむ回数である。PS-OCTと、該PS-OCTで得られた計測データを処理するコンピュータとを備えたPS-OCTシステムであって、前記コンピュータは、入力装置、出力装置、CPU及び記憶装置を備え、PS-OCTで計測した位相差の分布のデータを、非線形の変換関数によって、対称な位相差の分布のデータに変換し、平均することにより、真の位相量を求める手段として機能することを特徴とするPS-OCTシステム。PS-OCTと、該PS-OCTで得られた計測データを処理するコンピュータとを備えたPS-OCTシステムであって、前記コンピュータは、入力装置、出力装置、CPU及び記憶装置を備えており、予め、0~πまでの複数の位相量それぞれについての分布から成る複数の位相差の分布のセットを、異なるESNR毎に複数のセット作成し記憶装置に記憶させる手段、PS-OCTで計測した位相差の分布のデータを対称な位相差の分布のデータに変換する非線形の変換関数の係数を決める手段、PS-OCTで計測した位相差の分布のデータからESNRを見積もり、該見積もられたESNRに基づき、前記異なるESNR毎に作成した位相差の分布の複数のセットのいずれかを選択する手段、PS-OCTで計測した位相差の分布のデータを対称な位相差の分布を、前記選択された位相差の分布のセットに応じて決められた係数の下で、前記変換関数によって対称な位相差の分布のデータに変換する手段、及び前記対称な位相差の分布のデータを平均することにより、真の位相量を算出する手段、として機能することを特徴とするPS-OCTシステム。前記非線形の変換関数は、次の式2であることを特徴とする請求項8記載のPS-OCTシステム。但し、式2中、δEは対称な位相差の分布であり、δMはPS-OCTで計測した位相量であり、bii=0、2、3、.........nは変換関数の係数であり、γはESNRである。前記変換関数の係数を決める手段は、前記非線形の変換関数について、0~πまで、π/mきざみで変え、次のm1の連立方程式である式3を作成し、該連立方程式を解いて、変換関数の係数biのテーブルを作る構成であることを特徴とする請求項9記載のPS-OCTシステム。但し、式3中のmは、πをきざむ回数である。前記連立方程式を解いて変換関数の係数を決める手段においては、式4の行列DMの疑似逆行列DMを数値的に求めて、式5で計算し、最適の変換関数の係数biを一義的に決める構成であることを特徴とする請求項10記載のPS-OCTシステム。但し、式4中のmは、πをきざむ回数である。前記変換関数の係数を決める手段は、変分法を用い、次の式6で示される二乗誤差が最小となるように最適の変換関数の係数biを一義的に決める構成であることを特徴とする請求項9記載のPS-OCTシステム。但し、式6中のmは、πをきざむ回数である。一般的に、OCTを生体計測へ応用することを考える場合、物体の屈折率分布の微分構造は、非破壊、高分解能で捉えることはできるものの、繊維状の構造繊維の伸長方向等や歯のエナメル質の相違に起因する複屈折による偏光依存性を持つ生物試料の測定においては、解像度の低下とともに、構造を捉えられない。この点、PS-OCTは、ある特定部分からの散乱光成分とある偏光状態の参照光とを干渉させて、その干渉成分には偏光特性が強く反映され、位相遅延情報リターデーション量も得ることができ、その結果、奥行き方向の断面のある特定部分の偏光情報、複屈折情報を捉えることが可能となる。従って、PS-OCTは、眼科、循環器、呼吸器、消化器など広く非侵襲な組織判別を必要とする臨床分野における検査装置として有用である。このように、PS-OCTは、組織の偏光特性を非侵襲に可視化する光トモグラフィー装置として、組織コントラストという定性観察には有用である。しかしながら、本発明者らは、PS-OCTの高精度化を目的に鋭意研究開発する過程において、PS-OCTは、病変部位を観察し、病期の正確な分類などの定量観察においては精度が必ずしも十分とは言えず、さらなる改善すべき点があるということが判った。本発明は、上記従来のPS-OCTについて得られた改善すべき点を解決することを目的とするものであり、PS-OCTにおいて計測されたデータを非線形に補正し、PS-OCTの定量解析能力を高め、病期の正確な定量診断等を可能とするプログラム及び該プログラムを搭載したPS-OCTシステムを実現することを課題とする。より具体的には、PS-OCTで計測されるリターデーション位相遅延量、位相差の量は誤差を含んでいて、真値の周りに分布これをリターデーションの分布、或いは位相差の分布と称すしている。この分布は一般的に対称ではなく、必ずしも真値が中心にあるとは限らない非対称である場合も含まれる。従って、標準的な平均値をとる位相解析では正しいリターデーションを得ることができない。本発明は、このような従来のPS-OCTの問題を解決しようとするものであり、リターデーションが誤差を含みノイズとなり、その分布が正規分布でもなく、真の値の周りに対象に分布していない場合であっても、モンテカルロシミュレーションでノイズの特性を解析することによって得られた分布変換関数を用いて、計測データを変換することにより、系統的な誤差を除去し、ノイズに埋もれた真の値を推測し、PS-OCTの画像をより鮮明に補正するプログラムを実現するとともに、このプログラムを搭載したPS-OCTシステムを実現することを課題とする。本発明は上記課題を解決するために、PS-OCTで得られた計測データを処理するコンピュータに搭載されるプログラムであって、コンピュータを、PS-OCTで計測したリターデーションについての位相差の分布のデータを、非線形の変換関数によって、対称な位相差の分布のデータに変換し、平均することにより、真の位相量を求める手段として機能させることを特徴とするプログラムを提供する。本発明は上記課題を解決するために、PS-OCTで得られた計測データを処理するコンピュータに搭載されるプログラムであって、コンピュータを、予め、0~πまでの複数の位相量それぞれについての分布から成る位相差の分布のデータのセットを、異なるESNR毎に複数のセット作成し記憶する手段、PS-OCTで計測した位相差の分布のデータを対称な位相差の分布のデータに変換する非線形の変換関数の係数を決める手段、PS-OCTで計測した位相差の分布のデータからESNRを見積もり、該見積もられたESNRに基づき、前記異なるESNR毎に作成した位相差の分布の複数のセットのいずれかを選択する手段、PS-OCTで計測した位相差の分布のデータを、前記選択された位相差の分布のセットに応じて決められた係数の下で、前記変換関数によって対称な位相差の分布のデータに変換する手段、及び前記対称な位相差の分布のデータを平均することにより、真の位相量を算出する手段、として機能させることを特徴とするプログラム。本発明は上記課題を解決するために、PS-OCTと、該PS-OCTで得られた計測データを処理するコンピュータとを備えたPS-OCTシステムであって、前記コンピュータは、入力装置、出力装置、CPU及び記憶装置を備え、PS-OCTで計測した位相差の分布のデータを、非線形の変換関数によって、対称な位相差の分布のデータに変換し、平均することにより、真の位相量を求める手段として機能することを特徴とするPS-OCTシステムを提供する。本発明は上記課題を解決するために、PS-OCTと、該PS-OCTで得られた計測データを処理するコンピュータとを備えたPS-OCTシステムであって、前記コンピュータは、入力装置、出力装置、CPU及び記憶装置を備えており、予め、0~πまでの複数の位相量それぞれについての分布から成る複数の位相差の分布のセットを、異なるESNR毎に複数のセット作成し記憶装置に記憶させる手段、PS-OCTで計測した位相差の分布のデータを対称な位相差の分布のデータに変換する非線形の変換関数の係数を決める手段、PS-OCTで計測した位相差の分布のデータからESNRを見積もり、該見積もられたESNRに基づき、前記異なるESNR毎に作成した位相差の分布の複数のセットのいずれかを選択する手段、PS-OCTで計測した位相差の分布のデータを対称な位相差の分布を、前記選択された位相差の分布のセットに応じて決められた係数の下で、前記変換関数によって対称な位相差の分布のデータに変換する手段、前記対称な位相差の分布のデータを平均することにより、真の位相量を算出する手段、として機能することを特徴とするPS-OCTシステムを提供する。前記非線形の変換関数は、後記する式2である。前記変換関数の係数を決める手段は、前記非線形の変換関数について、0~πまで、π/mきざみで変え、次のm1の連立方程式である後記する式3を作成し、該連立方程式を解いて、変換関数の係数biのテーブルを作る構成である。但し、mは、πをきざむ回数である。前記連立方程式を解いて前記変換関数の係数を決める手段においては、後記する式4の行列DMの疑似逆行列DMを数値的に求めて、後記する式5で計算し、最適の変換関数の係数biを一義的に決める構成としてもよい。前記変換関数の係数を決める手段は、変分法を用い、後記する式6で示される二乗誤差が最小となるように最適の変換関数の係数biを一義的に決める構成としてもよい。本発明によれば、PS-OCTにおいて計測されたデータを非線形に補正し、PS-OCTの定量解析能力を高めることができるので、従来のPS-OCTは、病変部位の形態観察等の用途に限られていたが、本発明のPS-OCTシステムは、病変部位の病期を含め、正確な定量診断を可能とし、コンピュータ診断のための有用な手段として利用することが可能となる。本発明に係るPS-OCTの計測データを補正するプログラム及び該プログラムを搭載したPS-OCTシステムを実施するための形態を図面を参照して、以下説明する。本発明に係るPS-OCTシステムは、PS-OCTと、PS-OCTで得られた画像データを処理する画像処理装置とを備えている。画像処理装置は、通常のコンピュータが使用され、本発明に係るプログラムは、PS-OCTで得られた画像データを補正する手段として機能させるものである。PS-OCTPS-OCTは、すでに特許第4344829号公報等で知られているが、本発明の前提となる技術であるから、その概要を説明する。本発明に係る偏光感受光画像計測装置は、Bスキャンと同時に同期して光源からの偏光ビーム偏光子により直線的に偏光されたビームをEO変調器偏光変調器、電気光学変調器によって連続的に変調し、この連続的に偏光を変調した偏光ビームを分けて、一方を入射ビームとして走査して試料に照射し、その反射光物体光を得ると共に、他方を参照光として、両者のスペクトル干渉によりOCT計測を行うものである。そして、このスペクトル干渉成分のうち、垂直偏光成分Hと水平偏光成分Vを同時に2つの光検出器で測定することにより、試料の偏光特性を表すジョーンズベクトルを得るH画像とV画像構成を特徴とするものである。図1は、本発明に係る偏光感受光画像計測装置の光学系の全体構成を示す図である。図1に示す偏光感受光画像計測装置1は、光源2、偏光子3、EO変調器4、ファイバーカプラー光カプラー5、参照アーム6、試料アーム7、分光器8等の光学要素を備えている。この偏光感受光画像計測装置1の光学系は、光学要素が互いにファイバー9で結合されているが、ファイバーで結合されていないタイプの構造フリースペース型であってもよい。光源2は、広帯域スペクトルを有するスーパールミネッセントダイオードSLD:SuperLuminessentDiodeを使用する。なお、光源2は、パルスレーザでもよい。光源2には、コリメートレンズ11、光源2からの光を直線偏光にする偏光子3、進相軸を45°の方向にセットされたEO変調器偏光変調器、電気光学変調器4、集光レンズ13及びファイバーカプラー5が、順次、接続されている。EO変調器4は、進相軸を45°の方向に固定して、該EO変調器4にかける電圧を正弦的に変調することで、進相軸とそれに直交する遅相軸との間の位相差リターデーションを連続的に変えるもので、これにより、光源2から出て偏光子3で縦直線偏光となった光がEO変調器4に入射すると、上記変調の周期で、直線偏光→楕円偏光→直線偏光.........などのように変調される。EO変調器4は、市販されているEO変調器を使用すればよい。ファイバーカプラー5には分岐するファイバー9を介して、参照アーム6と試料アーム7が接続されている。参照アーム6には、偏波コントローラpolarizationcontroller10、コリメートレンズ11、偏光子12、集光レンズ13及び参照鏡固定鏡14が、順次、設けられている。参照アーム6の偏光子12は、上記のとおり偏光状態を変調しても参照アーム6から戻ってくる光の強度が変化しないような方向を選択するために用いている。この偏光子12の方向直線偏光の偏光方向の調整は偏波コントローラ10とセットで行う。試料アーム7では、偏波コントローラ15、コリメートレンズ11、固定鏡24、ガルバノ鏡16、集光レンズ13が、順次、設けられ、ファイバーカプラー5からの入射ビームが2軸のガルバノ鏡16により走査されて試料17に照射される。試料17からの反射光は物体光として再びファイバーカプラー5に戻り、参照光と重畳されて干渉ビームとして分光器8に送られる。分光器8は、順次接続される偏波コントローラ18、コリメートレンズ11、偏光感受型体積位相ホログラフィック回折格子19、フーリエ変換レンズ20、偏光ビームスプリッター21及び2つの光検出器22、23を備えている。この実施の形態では、光検出器22、23として、ラインCCDカメラ1次元CCDカメラを利用する。ファイバーカプラー5から送られてくる干渉ビームは、コリメートレンズ11でコリメートされ、回折格子19によって干渉スペクトルに分光される。回折格子19で分光された干渉スペクトルビームは、フーリエ変換レンズ20でフーリエ変換され偏光ビームスプリッター21で水平及び垂直成分に分けられ、それぞれ2つラインCCDカメラ光検出器22、23で検出される。この2つラインCCDカメラ22、23は、水平および垂直偏光信号両方の位相情報を検知するために使われるので、2つのラインCCDカメラ22、23は同一の分光器の形成に寄与するものでなくてはならない。なお、光源2、参照アーム6、試料アーム7及び分光器8には、それぞれ偏波コントローラ10、15、18が設けられているが、これらは、光源2から参照アーム6、試料アーム7、分光器8に送られるそれぞれのビームの初期偏光状態を調整して、EO変調器4で連続的に変調された偏光状態が、参照光と物体光においても互いに一定の振幅と一定の相対偏光状態の関係が維持され、さらにファイバーカプラー5に接続された分光器8において一定の振幅と一定の相対偏光状態を保たれるようにコントロールする。また、2つラインCCDカメラ22、23を含む分光器8を校正するときはEO変調器4は止める。参照光をブロックし、スライドガラスと反射鏡を試料アーム7におく。この配置は水平および垂直偏光成分のピークの位置が同じであることを保証する。そして、スライドガラスの後ろの面と反射鏡からのOCT信号は2つの分光器8で検知される。OCT信号のピークの位相差はモニターされる。この位相差はすべての光軸方向の深さでゼロであるべきである。次に、信号は2つラインCCDカメラ22、23を含む分光器8で複素スペクトルを得るために、ウィンドウされ逆フーリエ変換される。この位相差はすべての周波数でゼロであるべきなので、これらの値をモニターすることによって2つラインCCDカメラ22、23の物理的な位置は位相差が最小になるようにアライメントされる。PS-OCTの原理:本発明の特徴は、次のとおりである。光源2からの光を直線偏光し、この直線偏光されたビームをEO変調器4により連続的に偏光状態の変調を行う。即ち、EO変調器4は、進相軸を45°の方向に固定して、該EO変調器4にかける電圧を正弦的に変調することで、進相軸とそれに直交する遅相軸との間の位相差偏光角:リターデーションを連続的に変えるもので、これにより、光源2から出て直線偏光子で縦直線偏光となった光がEO変調器4に入射すると、上記変調の周期で、直線偏光→楕円偏光→直線偏光.........などのように変調される。そして、直線偏光された偏光ビームをEO変調器4により連続的に偏光状態の変調を行うと同時に、Bスキャンを同期して行う。即ち、1回のBスキャンの間に、EO変調器4による偏光の連続的な変調を複数周期行う。ここで、1周期とは、偏光角リターデーションφが0~2πと変化する期間である。要するに、この1周期の間に、偏光子からの光の偏光が、直線偏光垂直偏光→楕円偏光→直線偏光水平偏光.........などのように連続的に変調する。このように偏光ビームの偏光を連続的に変調しながら、試料アーム7では、入射ビームをガルバノ鏡16により試料17に走査してBスキャンを行い、分光器8において、その反射光である物体光と参照光の干渉スペクトルについて、その水平偏光成分および垂直偏光成分を2つのラインCCDカメラ22、23で検出する。これにより、1回のBスキャンによって、それぞれ水平偏光成分及び垂直偏光成分に対応する2枚のA-Bスキャン画像が得られる。上記のとおり、1回のBスキャンの間に、偏光ビームの偏光の連続的な変調を複数周期行うが、各周期1周期の連続的な変調の間に2つのラインCCDカメラ22、23で検出した水平偏光成分および垂直偏光成分それぞれの偏光情報が1画素分の偏光情報となる。1周期の連続的な変調の間に2つのラインCCDカメラ22、23で偏光情報を検出タイミング信号に同期して行い、1周期に検出回数取込回数を、4回、8回等、適宜決めればよい。このようにして1回のBスキャンの間に得た2枚のA-Bスキャン画像のデータを、Bスキャン方向に1次元フーリエ変換を行う。すると、0次、1次、-1次のピークが出る。ここで、0次のピークをそれぞれ抽出し、そのデータのみを用いて逆フーリエ変換すると、H0、V0画像が得られる。同様に、1次のピークをそれぞれ抽出し、そのデータのみを用いて逆フーリエ変換すると、H1、V1画像が得られる。H0、H1画像から、試料17の偏光特性を示すジョーンズマトリックスの成分のうち、J1,1およびJ1,2を求める事ができる。そして、V0、V1画像から、試料17の偏光特性を示すジョーンズマトリックスの成分のうち、J2,1およびJ2,2を求める事ができる。なお、ここでは、ジョーンズマトリックスの成分J1,1、J1,2、J2,1、J2,2を特許文献1の記載に基づいて記載したが、それぞれ、本発明の後記する式1のジョーンズマトリックスの成分I0Hz、I1Hz、I0Vz、I0Vzに相当する。このようにして、1回のBスキャンにおいて4つの偏光特性を含む情報が得られる。そして、この4つの情報をそれぞれ、通常のFD-OCTと同様にAスキャン方向にフーリエ変換すると、1次のピークが試料17の深さ方向の情報を有し、しかもそれぞれ偏光特性に応じた4枚のA-B画像の計測データが得られる。本発明のプログラム及びPS-OCTシステム上記構成のPS-OCTで得られた計測データは、画像処理装置として使用するコンピュータ30に入力される。このコンピュータ30は通常のコンピュータであり、図2に示すように、入力部31、出力部32、CPU33、記憶装置34及びデータバス35を備えている。本発明に係るPS-OCTシステムに搭載されるプログラムは、コンピュータ30の記憶装置34に記憶されるプログラムであって、コンピュータ30に入力されたPS-OCTで得られた画像の計測データを補正し、PS-OCTの画像をより鮮明に補正する手段として、コンピュータ30を機能させるプログラムである。このプログラムを搭載することにより、本発明に係るPS-OCTシステムは、画像をより鮮明に補正する手段を備えたシステムとなる。即ち、従来はPS-OCTで計測された生データから得られた画像を、直接、診断その他の各種の用途に用いていた。本発明のプログラム及びPS-OCTシステムでは、PS-OCTで計測された生データを、非線形関数で補正して高精度な生体複屈折計測を可能とするものであり、この非線形関数は、モンテカルロ法MonteCarloanalysisで、予め本発明のプログラムを用い、コンピュータ30においてシミュレーションで設計される。以下、本発明のプログラム及び該プログラム搭載したPS-OCTシステムを詳細に説明する。図3は、リターデーション位相差、位相遅延量、の分布図位相差の分布図である。横軸は位相であり、縦軸はリターデーションの生じる頻度である。図3の右欄は、PS-OCTで計測したリターデーションに基づく位相差の分布であり、図3の左欄は、シミュレーションで作成した位相差の分布である。位相差の分布は、横軸をπ/mずつ刻んだm1個のデータm1個の位相差の分布を1セットとして、ESNRを変えて複数のセット作成する。例えば、m=60として、横軸をπ/60ずつ刻んだ61個のデータを1セットとして、ESNRを変えて複数のセット作成する。図3では、2つの異なるESNRのセットについて、それぞれπ/2及びπについての位相差の分布を示す図である。ところで、PS-OCTによるリターデーションのノイズモデルは、ジョーンズマトリックスのそれぞれの要素に、加算的なノイズを設定して、次の式1で示される。ここで、Sは真のリターデーションの値であり、Nはノイズであり、Iは測定値である。添え字の0、1は、計測チャンネル図1の光検出器22、23参照である。前記段落0052に説明したジョーンズマトリックスの成分は、真の値SにノイズNが加わった測定値Iに対応する。ここでzは深さ方向の座標で、ジョーンズマトリックスのそれぞれの要素は表面z=0から深さzまでのリターデーション等の積分量に相当する。計測データのノイズを表す評価値として、通常は、S/N比SN比、SNR等とも呼ばれているを使用するが、本発明では、より効果的にノイズを表す評価値としてESNRを使用する。本明細書では、ESNRはγとして表し、次の式で定義される。1/γ=1/41/SNRS,01/SNRS,11/SNRr,01/SNRr,1なお、この式中、sはPS-OCTのサンプルアームであり、rは参照アームである。0と1は、回数である。複数の独立したSNRをもつものを、単一の値で示すためにESNRを導入した。添え字の0、1は、計測チャンネル図1の光検出器22、23参照である。SNRS,0、SNRS,1、SNRr,0、SNRr,1は、それぞれ添え字の組み合わせによって、サンプルアームから0チャンネルへ寄与するSNR、サンプルアームから1チャンネルへ寄与するSNR、参照アームから0チャンネルへ寄与するSNR、参照アームから1チャンネルへ寄与するSNRを表す。図3の右欄の上から2番目の図に示すように、PS-OCTで計測されるリターデーション位相遅延量、位相差は誤差を含んでいて、真値π/2の周りに分布している。そして、ESNRが小さくなると、図3の右欄の上から1番目の図に示すように、誤差の幅を大きくなる。また、この位相差の分布は、図3の右欄の上から3、4番目の図に示すように、必ずしも、対称ではなく、真値が中心にあるとは限らず、非対称である場合も含まれる。要するに、リターデーションが誤差を含みノイズとなり、位相差の分布が正規分布でもなく、真の位相量の周りに対象に分布していない場合も含まれる。このような場合は、標準的な平均値をとる位相解析では正しいリターデーションを得ることができない。本発明に係るプログラムは、モンテカルロシミュレーションでノイズの特性を解析することによって得られた分布変換関数を用いて、PS-OCTで計測された画像の計測データを変換することにより、系統的な誤差を除去し、ノイズに埋もれた真の位相量を推測し、PS-OCTの画像をより鮮明に補正する手段としてコンピュータ30を機能させるプログラムであり、本発明のシステムは、このプログラムを搭載し上記補正を行うPS-OCTシステムである。原理:本発明に係るプログラムは、PS-OCTで計測された画像の計測データによる実測のリターデーションの分布位相差の分布を補正し、真の位相量δT真のリターデーションの値を求める手段とし、コンピュータ30を機能させるものであるが、まず、その原理アルゴリズムを、以下に説明する。本発明では、PS-OCTで計測された画像の計測データによる実測の位相量δMから、非線形関数fを用いて対称な位相差の分布δEに変換する。要するに、対称な位相差の分布δEは、実測の位相量δMの冪関数であるδE=fδMとして表される。この関数を本明細書では、変換関数と言い、次の式2で表される。これは、計測値の非対称な分布δMを、正しい値を重みの中心とする分布δEに変換するために行われる非線形な写像を意味している。この式2中、biは変換関数の係数である。最適な変換関数の係数biiは、0、1、2、............nは、ESNRγの関数であるが、次のように、連立方程式を解いて決めるか、或いは最小二乗法を用いて決め、biのテーブルを作っておく。そのために、位相値を、0からπまでπ/mきざみで変え、モンテカルロ法によってδMを求める。係数biは、式3に示す連立方程式を解いて決める。式3中、δM,π/mは、π/mきざみ、例えば、π/60きざみでシミュレートした、0からπにおけるリターデーションのシミュレーション値δMを示す。ここで、mは、例えば、1、2、......60である。なお、式3中のの意味は、多数のモンテカルロシミュレーションによって得られた値の平均値を示す。具体的には、式3を、式4に示すように行列で書き下し、係数biを求める。そして、係数を決めるには、式4のシミュレーションで得られたδM,π/mnのn×mの行列DMの疑似逆行列DMを数値的に求めて、次の式5で計算し係数を求める。なお、式4中のの意味は、多数のモンテカルロシミュレーションによって得られた値の平均値を示す。最小二乗法を用いて係数を決める手段は、変分法を用い、次の式6で示す、二乗誤差R2が最小になるように、モンテカルロ法によって、最適の変換関数の係数biを一義的に決めることができる。ここでも、位相値を、0からπまでπ/mきざみで、例えばm=60として、π/60きざみで変え、モンテカルロ法によってシミュレーション値δMを求める。kは、位相をこきざみに計測した計測回を示し、0、1、2......mであり、例えばm=60等である。最適の変換関数の係数biが決まれば、対称な位相差の分布δEを表す変換関数δE=fδMが決まり、PS-OCTによる実測の位相量δMを代入すれば、δEが求まる。そして、δEを次の式7に示すように平均する式7中のは平均の意味ことにより、真の位相量δTを求めることができる。平均は、一つの位相量を求めるために行った実測回数でおこなう。構成:本発明の原理は以上のとおりであるが、この原理に基づき、本発明のPS-OCTのシステムにおいては、本発明に係るプログラムは、画像処理装置として使用するコンピュータ30を次の手段として機能させる。1位相差の分布の作成手段位相差の分布の作成手段は、予め、図3の左欄に示すようなリターデーションについての位相差の分布をシミュレーションで複数作成する手段である。具体的には、作成する位相差の分布は、π/m、例えばm=60として、横軸をπ/60ずつ刻んだ、61個のデータ61個の位相差の分布を1セットとして、ESNRを変えて位相差の分布のセットを複数セット作成する。この複数セットのデータは、作成後、記憶装置34に記憶される。2変換関数の係数の決定手段変換関数の係数の決定手段は、複数の位相差の分布のセットについて、それぞれ変換関数の係数biを決める手段である。この変換関数の係数biは、前記したとおり、式3の連立方程式を解いて決めるか、或いは式6において最小二乗法を用いて決める。biが決められたら、そのテーブルを作成する。この変換関数の係数biのテーブルは、記憶装置34に記憶される。3位相差の分布の選択手段位相差の分布の選択手段は、PS-OCTで計測された画像の実測の位相差の分布、位相量の補正に先立ち、その計測データ及びシステム設計等からESNRの値を見積もり、上記1に示す複数の位相差の分布図3の左欄のような位相差の分布のセットについて、どのセットを用いるかを決定する。4変換関数による位相差の分布の変換手段変換関数の算出手段は、前記記憶装置34に記憶された変換関数の係数bi又はそのテーブルのうち、選択手段で選択された位相差の分布のセットに対応する係数bi又はそのテーブルを用いて、式2のbiに代入し、PS-OCTで計測した実測のリターデーションの位相量δMから、対称な位相差の分布を表す変換関数δE=fδMを求める手段である。5真の位相量算出手段真の位相量算出手段は、上記のとおり求めた対称な位相差の分布を表すδE=fδMから、前記式7に基づき平均することにより、真の位相量δTを算出する手段である。試験例1本発明者らは、本発明のプログラム及び該プログラム搭載したPS-OCTシステムを評価するために、シミュレーション値と実測値を比較する比較試験を行った。この比較試験では、試料計測対象物は、ガラス板、1/8波長板及び1/4波長板であり、それぞれに生じるリターデーションの真値は、0、π/2及びπである。なお、シミュレーション値は、前記したとおり、予めbiの係数をモンテカルロシミュレーションで決めておき、これによって得られた係数を用いてこのシミュレーションデータを変換したものである。実測データを上記係数で変換したものが実測値である。モンテカルロシミュレーションにはランダムノイズが入っていることと、biを有限の項数で打ち切っているので、シミュレーションといえども完全には一致はしない。この試験では、PS-OCTによる計測のために、1.3μmのプローブを使用した。このPS-OCTは、空気中での深度解析は8.3μmであり、30,000ライン/秒の測定速度を有する。計測は、試料の前方に光強度を変えることの可能な非偏光フィルタを置くことによって、ESNRを変えて行った。図4は、シミュレーションと試験を比較した比較試験の結果を示す。点線は、真のリターデーションの値0、π/2及びπを示す。実線b、c、eは、本発明による補正を行わずに、リターデーションの実測値を単に平均のみをとったものである。実線a、d、fは、本発明によるシミュレーションの結果であり、係数を作るのに用いたシミュレーション数は65536個であり、これは、係数を求めるためのデータでもあり、結果を得るために用いた位相データでもある。実線b、c、e及びa、d、fの近辺に位置する●■のマークは、それぞれ試験によって得られたリターデーションの実測値及び該実測値を本発明により補正した試験補正値あって、ここでは64個の測定データを用いている。シミュレーションと試験補正値は、絶対値及びESNR依存性等において、ほぼ一致している。試験例2本発明者らは、本発明の評価をするために、PS-OCTの評価に通常使用されている鶏の胸肉ex-vivoを試料として計測する比較試験を行った。この試験では、試料の10mmの部分について、512本のAラインから成るBスキャンを行った。要するに、A方向深さ方向の計測で得られるAラインを、BスキャンA方向に直交するB方向のスキャンによって、512本得る。しかも、この計測は、試料について同じ箇所を128回行った。8dbより小さいESNRは不安定であるので除去した。図5は、試験例2の試験結果を示す画像図である。図5aはOCT強度像を示し、図5bは段落0064のESNRの値を画像濃度で表した像を示し、図5cは単純平均位相像を示し、図5dは本発明により補正され最適化された画像を示す。試料の深い部分では、図5cの単純平均位相像ではコントラストが低下するが、図5dはそれほどコントラストは低下しないという結果が得られた。図6は、図5のc、dにおける点線で示す断面像であり、細線はESNRの値を示しており、太い点線はリターデーションを示している。図5c、dにおいて、ESNRの値細線で示すは、試料の奥深い部分にいくほど低くなっていく、要するに信号強度が小さくなりノイズが大きくなることが示されている。リターデーションは、単純平均による図6cに示すものより、本発明を用いた図6dに示すものの方が、より明確に示されている。特に、図6dに示すものは、0及びπ付近に近い値が再生されている事がわかる。これによりPS-OCTの画像の解析において、本発明が有用であるということが評価された。また、変換を行わない処理図5c参照では0およびπ付近の位相データが欠落し位相変化の連続性が失われているが、本発明によって処理した結果図5d参照では、0およびπ付近の位相も再現されていて深さ方向のリターデーションの絶対値が正しく計測されていることがわかる。試験例3本発明者らは、本発明の評価をするために、PS-OCTにより、人の眼底in-vivoについて計測を行った。この計測は、1.0μmのプローブビームを用いて、11.0μmの眼底深さについて、30,000ライン/秒の測定速度を有するPS-OCTによって行った。図7は、この試験によって得られた画像を示す。図7aはこれまでの単純平均処理画像であり、図7bは本発明による補正処理をおこなったものである。図7bに示すように、本発明によれば、位相量が0の表面付近が正確に再生されていて、より真実に近い位相量、位相差の分布が得られているものと考えられる。以上、本発明に係るプログラム及び該プログラムを搭載したPS-OCTシステムを実施するための形態を説明したが、本発明はこのような実施の形態に限定されるものではなく、特許請求の範囲に記載された技術的事項の範囲内でいろいろな実施の形態があることは言うまでもない。"