(57)【要約】【課題】  微分位相値のそれぞれが複数の画素の強度情報を用いて取得された微分位相像を積分して位相像を取得する演算装置であって、正確性が高い位相像を取得することができる演算装置を提供すること。【解決手段】  演算装置１７０は、微分位相像を積分して位相像を取得する。微分位相像は、微分干渉計により形成された１つの強度分布に含まれる複数の強度情報を用いて取得された１つの微分位相値を複数有する微分位相像である。演算手段１７０は、微分位相像のうち第１の領域の微分位相値を中央部の微分位相値よりも軽く重み付けすることで微分位相値像の重み付けを行う重み付け手段２１０と、重み付けされた微分位相像を積分して位相像を取得する積分手段２２０とを有する。  第１の領域は、微分位相像の端部のうち少なくとも１つ以上の画素を有する。【選択図】  図１
2015213574.tif 000002
【特許請求の範囲】
【請求項１】
微分位相像を積分して位相像を取得する演算装置であって、
微分干渉計により形成された１つの強度分布に含まれる複数の画素の強度情報を用いて取得された１つの微分位相値を複数有する微分位相像の重み付けを行うことで重み付けされた微分位相像を取得する重み付け手段と、
前記重み付けされた前記微分位相像を積分して位相像を取得する積分手段とを有し、
前記重み付け手段は、
前記微分位相像の端部のうち少なくとも１つ以上の画素を有する、第１の領域の微分位相値を、前記微分位相像の中央部の微分位相値よりも軽く重み付けすることを特徴とする演算装置。
【請求項２】
前記重み付け手段は、
入力された前記微分位相像から、前記第１の領域の微分位相値を削除することで、前記重み付けされた微分位相像を取得し、
前記重み付けされた微分位相像を前記積分手段へ出力することを特徴とする請求項１に記載の演算装置。
【請求項３】
前記重み付け手段は、
入力された前記微分位相像のうち、前記第１の領域の微分位相値を前記微分位相値より絶対値が小さな値に置き換えることで、前記重み付けされた微分位相像を取得し、
前記重み付けされた微分位相像を前記積分手段へ出力することを特徴とする請求項１に記載の演算装置。
【請求項４】
前記重み付け手段は、入力された前記微分位相像にフィルターをかけることで、前記重み付けされた微分位相像を取得し、
前記重み付けされた微分位相像を前記積分手段へ出力することを特徴とする請求項１に記載の演算装置。
【請求項５】
前記重み付け手段は、
前記重み付けされた微分位相像のうち前記第１の領域の微分位相値の絶対値の平均値が、前記入力された微分位相像のうち前記第１の領域の微分位相値の絶対値の平均値の１／２以下になるように、前記微分位相像の重み付けを行うことで前記重み付けされた微分位相像を取得することを特徴とする請求項３に記載の演算装置。
【請求項６】
前記重み付け手段は、前記中央部の微分位相値に掛けられるフィルターの値を１としたときに、前記第１の領域の微分位相値に掛けられるフィルターの値の平均値が０．５以下になるように、前記微分位相値の重み付けを行うことで前記重み付けされた微分位相像を取得することを特徴とする請求項４に記載の演算装置。
【請求項７】
前記微分位相像はｘ方向の微分位相像であり、
前記積分手段は、前記重み付けされた微分位相像をｘ方向に積分し、
前記重み付け手段は、前記微分位相像の端部のうち、前記ｘ方向と垂直に交わるｙ方向に配列した画素列の画素を１つ以上有する領域を前記第１の領域とすることを特徴とする請求項１乃至６のいずれか１項に記載の演算装置。
【請求項８】
前記微分位相像は、ｘ方向の微分位相像とｙ方向の微分位相像であり、
前記積分手段は、前記重み付けされたｘ方向の微分位相像をｘ方向に積分し、前記重み付けされたｙ方向の微分位相像をｙ方向に積分し、
前記重み付け手段は、前記ｘ方向の微分位相像と前記ｙ方向の微分位相像のそれぞれの端部のうち、前記ｘ方向と垂直に交わるｙ方向に配列した画素列の画素を１つ以上と、前記ｘ方向に配列した画素行の画素を１つ以上とを有する領域を、前記第１の領域とすることを特徴とする請求項１乃至６のいずれか１項に記載の演算装置。
【請求項９】
微分位相像を積分して位相像を取得する演算装置であって、
微分干渉計により形成された１つの強度分布に含まれる複数の画素の強度情報を用いて取得された１つの微分位相値を複数有する微分位相像の重み付けを行うことで重み付けされた微分位相像を取得する重み付け手段と、
前記重み付けされた微分位相像を積分して位相像を取得する積分手段とを有し、
前記重み付け手段は、
前記微分位相像うち、前記１つの強度分布に含まれる不連続な画素の強度情報を用いて取得された微分位相値を、前記微分位相像のうち、前記１つの強度分布に含まれる連続する画素の強度情報のみを用いて取得された微分位相値よりも軽く重み付けすることを特徴とする演算装置。
【請求項１０】
前記重み付け手段は、
入力された前記微分位相像から、前記１つの強度分布に含まれる不連続な画素の強度情報を用いて取得された微分位相値を削除することで、前記重み付けされた微分位相像を取得し、
前記重み付けされた微分位相像を前記積分手段へ出力することを特徴とする請求項９に記載の演算装置。
【請求項１１】
前記重み付け手段は、前記重み付けされた微分位相像のうち前記不連続な画素の強度情報を用いて取得された微分位相値の絶対値の平均値が、前記入力された微分位相像のうち前記不連続な画素の強度情報を用いて取得された微分位相値の絶対値の平均値の１／２以下になるように、前記微分位相像の重み付けを行うことで前記重み付けされた微分位相像を取得することを特徴とする請求項９に記載の演算装置。
【請求項１２】
前記重み付け手段は、前記連続する画素の強度情報のみを用いて取得された微分位相値に掛けられるフィルターの値の平均値を１としたときに、前記不連続な画素の強度情報を用いて取得された微分位相値に掛けられるフィルターの値の平均値が０．５以下になるように、入力された前記微分位相像にフィルターをかけることで前記微分位相像の重み付けを行い、前記重み付けされた微分位相像を取得することを特徴とする請求項９に記載の演算装置。
【請求項１３】
少なくとも前記１つの強度分布を用い、前記微分位相像を取得する手段を備えることを特徴とする請求項１乃至１２のいずれか１項に記載の演算装置。
【請求項１４】
微分位相像を積分して位相像を取得する演算装置であって、
強度分布を用いて微分位相像を取得する手段と、
前記微分位相像を積分して位相像を取得する積分手段とを有し、
前記微分位相像を取得する手段は、
前記微分位相像のそれぞれの微分位相値は、微分干渉計により形成された１つの強度分布のうち、連続する複数の画素の強度情報を用いて取得され、
前記１つの強度分布のうち、端部に対応する領域の微分位相値を取得しないことを特徴とする演算装置。
【請求項１５】
前記微分位相像を取得する手段は、
第１の位置の微分位相値を取得するために前記１つの強度分布から取得される、連続する複数の強度情報の数が、１つの強度分布から取得する必要がある強度情報の数よりも小さいとき、
前記第１の位置の微分位相値を取得しないことを特徴とする請求項１４に記載の演算装置。
【請求項１６】
前記１つの強度分布の少なくとも一部をフーリエ変換することにより前記微分位相像を取得することを特徴とする請求項１３乃至１５のいずれか１項に記載の演算装置。
【請求項１７】
トールボット干渉計又はトールボット・ロー干渉計と、
請求項１３乃至１６のいずれか１項に記載の演算装置とを備え、
前記演算装置は、前記トールボット干渉計又は前記トールボット・ロー干渉計により形成された強度分布から前記微分位相像を取得することを特徴とする位相値取得システム。
【請求項１８】
微分干渉計により形成された１つの強度分布に含まれる複数の強度情報を用いて取得された１つの微分位相値を複数有する微分位相像のうち、第１の領域の微分位相値を前記微分位相像の中央部の微分位相値よりも軽く重み付けする重み付け工程と、
重み付けされた前記微分位相像を積分して位相像を取得する積分工程とを有し、
前記第１の領域は、前記微分位相像の端部のうち少なくとも１つ以上の画素を有することを特徴とする位相像の取得方法。
【請求項１９】
微分位相像を取得する工程と、
前記微分位相像を積分して位相像を取得する積分工程とを有し、
前記微分位相像を取得する工程において、
前記微分位相像のそれぞれの微分位相値は、微分干渉計により形成された１つの強度分布のうち、連続する複数の画素の強度情報を用いて取得され、
前記１つの強度分布のうち、端部に対応する領域の微分位相値は取得されないことを特徴とする位相像の取得方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、位相像を取得する演算装置及び位相像の取得方法に関する。
【背景技術】
【０００２】
微分干渉計（シアリング干渉計とも呼ばれる）は、光源から照射されたコヒーレントな光を分割し、一方に被検体による波面の歪曲を生成し、もう一方の波面をわずかにずらすことで周期的な強度分布（いわゆる干渉像）を形成し、検出する干渉計である。この強度分布の変化から被検体による光の位相変化を取得することができる。尚、光の代わりにＸ線のような光以外の電磁波を用いたり、電子線を用いたりすることも可能である。尚、周期的な強度分布とは、強度分布中の周期が一定である強度分布に限らない。例えば、強度分布の中央に近いほど周期が大きかったり、小さかったりするように、周期が変化する強度分布であっても、明部と暗部とが配列した強度分布であれば周期的な強度分布であるとみなす。
【０００３】
強度分布の変化から被検体による光の位相変化に関する情報（つまり、被検体の位相情報）を取得する方法、すなわち位相回復法の１つとして、非特許文献１に記載されているようなフーリエ変換法が知られている。フーリエ変換法は、強度分布をフーリエ変換し、そのキャリア周波数と一致したスペクトルの周辺の情報から強度分布の位相変化量を取得する方法である。
【０００４】
微分干渉計の場合、強度分布の位相の変化量は、被検体による光の位相変化量の微分した値である。よって、位相回復により得られた強度分布の位相の変化量を積分することで、被検体による光の位相変化量を取得することができる。尚、本明細書では、被検体による光の位相変化量を微分した値の一つ一つを微分位相値、微分位相値の空間分布を微分位相像、被検体による光の位相変化量の一つ一つを位相値、位相値の空間分布を位相像と呼ぶ。
【０００５】
微分位相像を積分して位相像を取得する方法として、いくつかの方法が提案されている。例えば、単純に微分位相像の微分位相値を順次積算することで位相像を取得することが可能である。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】Ｊ．Ｏｐｔ．Ｓｏｃ．Ａｍ．，７２，１５６（１９８２）
【発明の概要】
【発明が解決しようとする課題】
【０００７】
フーリエ変換法は、各画素における強度分布の位相の変化量を取得するために、その画素の周辺の画素の強度情報も用いる。このように、各画素における強度分布の位相の変化量を取得するために、その画素の周辺の画素の強度情報も用いると、強度分布の端部においては各画素における微分位相値が正しく取得できないことがある。なぜならば、強度分布の端部においては、用いる周辺画素のデータの一部が不足するためである。フーリエ変換法では解析するデータに周期性を仮定しているため、例えば、強度分布の右端の画素における微分位相値を取得するために、強度分布の左端の画素の強度情報も用いる。このように、不連続な画素の強度情報を用いて取得された微分位相値は、連続する画素の強度情報のみを用いて取得された微分位相値に比べて、微分位相値の正確さが低くなることが理解できる。これは、１画素分の微分位相値を取得するために、複数の画素の強度情報を用いる方法であれば、フーリエ変換法以外の方法であっても同様である。
【０００８】
一方、位相像は、微分位相像を積分して取得されるため、微分位相像の各画素の微分位相値が正しい値からずれると、積分像においてはそのずれた値が積算されていく。よって、強度分布の端部における微分位相値の正確さの低下の影響が、位相像においては広範囲に広がる可能性がある。
【０００９】
そこで本発明は、微分位相値のそれぞれが複数の画素の強度情報を用いて取得された微分位相像を積分して位相像を取得する演算装置であって、位相像における微分位相値の正確性の低下の影響を低減することができる演算装置を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明の一側面としての演算装置は、微分位相像を積分して位相像を取得する演算装置であって、微分干渉計により形成された１つの強度分布に含まれる複数の画素の強度情報を用いて取得された１つの微分位相値を複数有する微分位相像の重み付けを行うことで重み付けされた微分位相像を取得する重み付け手段と、前記重み付けされた前記微分位相像を積分して位相像を取得する積分手段とを有し、前記重み付け手段は、前記微分位相像の端部のうち少なくとも１つ以上の画素を有する、第１の領域の微分位相値を、前記微分位相像の中央部の微分位相値よりも軽く重み付けすることを特徴とする。本発明のその他の側面については、以下で説明する実施の形態で明らかにする。
【発明の効果】
【００１１】
本発明により、微分位相値のそれぞれが複数の画素の強度情報を用いて取得された微分位相像を積分して位相像を取得する演算装置であって、位相像における微分位相値の正確性の低下の影響を低減することができる演算装置を提供することができる。
【図面の簡単な説明】
【００１２】
【図１】本発明に係る演算装置の機能ブロック図及び演算装置が行う取得方法の工程図。
【図２】第１実施形態の概要を説明する図。
【図３】第２実施形態の概要を説明する図。
【図４】比較例１〜２及び実施例１〜４のシミュレーションに使用した被検体を示す図。
【図５】比較例１〜２及び実施例１〜４のシミュレーションに使用した周期パターンを示す図。
【図６】比較例１により得られた微分位相像を示す図。
【図７】実施例１により得られた重み付けされた微分位相像及び位相像を示す図。
【図８】実施例２により得られた重み付けされた微分位相像及び位相像を示す図。
【図９】実施例３により得られた重み付けされた微分位相像及び位相像を示す図。
【図１０】実施例４により得られた重み付けされた微分位相像及び位相像を示す図。
【図１１】比較例２及び実施例３、４で使用したフィルターを示す図。
【図１２】比較例１により得られた位相像を示す図。
【図１３】比較例２により得られた微分位相像及び位相像を示す図。
【図１４】第１〜第４実施形態に係るトールボット・ロー干渉計の概略図。
【発明を実施するための形態】
【００１３】
以下、図面を用いて本発明の好ましい実施形態（１〜４）について説明する。図１（ａ）は、実施形態１〜３に係る演算装置の機能ブロック図である。実施形態１〜３に係る演算装置１７０は、微分位相像（微分位相値の空間分布）を取得する手段２５０と重み付け手段２１０と積分手段２２０から構成される。
【００１４】
微分位相像を取得する手段２５０は、強度分布の入力を受け、強度分布を用いて位相回復を行うことで、微分位相像を取得し、重み付け手段へ微分位相像を出力する。
【００１５】
重み付け手段２１０は、微分位相像の入力を受けて微分位相像を取得し、取得した微分位相像の重み付けを行うことで重み付けされた微分位相像を取得する。そして、重み付けされた微分位相像を積分手段へ出力する。取得する微分位相像が有する微分位相値のそれぞれは、１つの強度分布に含まれる複数の画素の強度情報を用いて取得される。重み付け手段２１０による重み付けは、入力された微分位相像の端部のうち、１つ以上の画素を有する領域の微分位相値を、中央部の微分位相値よりも軽く重み付けする。但し、本発明及び本明細書において端部とは、微分位相像の右端、左端、上端、下端の画素ことを指し、例えば、微分位相像が５１２×５１２画素であれば、端部は５１２×４−４（重複部分）＝２０４４画素存在することになる。また、本発明及び本明細書において中央部とは、微分位相値の中央の画素のことを指し、例えば、微分位相像が５１１×５１１画素であれば、上下左右の端から２５６個目の画素のことを指す。また、微分位相像の画素数が、５１２×５１２画素のように、縦横で偶数の場合、中央部とは、微分位相像の中央の点を囲む４つの画素のことを指す。また、画素数が偶数×奇数の場合、中央部とは、微分位相像の中央点を挟む２つの画素のことを指す。中央部の微分位相値よりも軽く重み付けする微分位相値は、端部の微分位相値を１つ以上含めば良く、例えば、端部に隣接する画素の微分位相値を中央部よりも軽く重み付けしても良い。本発明及び本明細書では、微分位相値を、中央部の微分位相値よりも軽く重み付けする領域のことを、単に第１の領域と呼んだり、軽く重み付けする領域と呼んだりすることがある。軽く重み付けする領域の微分位相値を中央部の微分位相値よりも軽く重み付けするとは、軽く重み付けする領域の微分位相値の絶対値を小さくすることに限定されない。例えば、中央部の微分位相値の絶対値を大きくすることで、軽く重み付けする領域の微分位相値を軽く重み付けしても良い。また、軽く重み付けする領域の微分位相値を中央部の微分位相値よりも軽く重み付けするとは、軽く重み付けする領域の微分位相値を０にすることと、軽く重み付けする領域の微分位相値自体を微分位相像から削除することも含む。
【００１６】
軽く重み付けする範囲は、微分位相像の端部だけでなく、その周辺を含んでいても良い。上述のように、不連続な画素の強度情報を用いて取得された微分位相値の正確性は、中央部の微分位相値のように、連続する画素の強度情報のみを用いて取得された微分位相値の正確性よりも低い。よって、不連続な画素の強度情報を用いて微分位相値を取得する領域の少なくとも一部を、軽く重み付けする領域とする。また、中央部でも軽く重み付けする領域でもない領域は中央部と同等に重み付けすることが好ましい。尚、中央部でも軽く重み付けする領域でもない領域は、連続する画素の強度情報のみを用いて微分位相値を取得する領域のうち、中央部以外の全ての領域を有することが好ましい。また、中央部でも軽く重み付けする領域でもない領域は、不連続な画素の強度情報を用いて微分位相値を取得する領域の一部を有していても良い。つまり、不連続な画素の強度情報を用いて取得された微分位相値の少なくとも一部が、連続する画素の強度情報のみを用いて取得された微分位相値よりも軽く重み付けされることが好ましい。
【００１７】
尚、本発明及び本明細書において連続する画素の強度情報のみを用いて取得された微分位相値とは、連続する画素の強度情報のみを用いて取得された微分位相値のことである。つまり、位相回復時に、空間的に連続する画素の強度情報のみを強度情報として用いて取得された微分位相値である。
【００１８】
空間的に連続する画素とは、ｘまたはｙ座標が、ｎ、ｎ−１、ｎ−２の様に、連続する整数である画素群のことを指す。例えば、１つの強度分布のみを用いて微分位相値を取得する場合は、１つの強度分布において連続する画素の強度情報のみを用いて取得された微分位相値のことを指す。一方、複数の強度分布を用いて微分位相値を取得する場合は、用いる複数の強度分布のうち、１つ１つの強度分布からは連続する画素の強度情報のみを抽出し、抽出したそれらの強度情報を用いて取得された微分位相値のことを指す。
【００１９】
積分手段２２０は、重み付けされた微分位相像の入力を受けて重み付けされた微分位相像を取得し、この重み付けされた微分位相像を積分して位相像を取得する。このように、正確性の低い微分位相値を軽く重み付けしてから積分を行うことで、これらの微分位相値が位相像に与える影響を軽減することができる。
【００２０】
演算装置１７０が行う位相像の取得方法について図１（ｂ）に示す。演算装置１７０が行う位相像の取得方法は、微分位相像を取得する工程Ｓ１０と重み付け工程Ｓ１１と積分工程Ｓ１２とを有する。微分位相像を取得する工程Ｓ１０は、微分位相像を取得する手段２５０により行われ、重み付け工程Ｓ１１は、重み付け手段２１０により行われ、積分工程Ｓ１２は、積分手段２２０により行われる。
【００２１】
尚、演算装置１７０は、強度分布の入力を受ける代わりに、微分位相像の入力を受けても良い。この場合、演算装置は微分位相像を取得する手段を備えず、微分位相像を取得する工程は、外部から微分位相像の入力を受けることで行われる。
【００２２】
実施形態４は、少なくとも強度分布の端部のうち、少なくとも１つ以上の画素の微分位相値を算出しない。言い換えると、実施形態４は、強度分布全体の微分位相像を経ることなく、直接重み付けされた微分位相像を取得し、これを積分して位相像を取得する。図１（ｃ）は、実施形態４に係る演算装置の機能ブロック図である。実施形態４に係る演算装置２３０は、微分位相像を取得する手段２４０と積分手段２２０から構成される。微分位相像を取得する手段２４０は、強度分布を用いて微分位相像を取得する。このとき、端部の微分位相値を取得しないことで、重み付けされた微分位相像を取得する点が、実施形態１〜３の微分位相像を取得する手段２５０とは異なる。
【００２３】
積分手段２２０は、実施形態１〜３における積分手段と同様であり、微分位相像を取得する手段２４０で取得された微分位相像を積分して位相像を取得する。実施形態４の演算装置が行う位相像の取得方法（不図示）は、微分位相像を取得する工程と、積分工程をと有する。微分位相像を取得する工程は微分位相像を取得する手段２４０により行われ、積分工程は、積分手段２２０により行われる。
【００２４】
各実施形態では、微分干渉計としてＸ線トールボット・ロー干渉計を用い、Ｘ線トールボット・ロー干渉計により得られる強度分布を用いて位相像を取得する形態について説明をするが、他の微分干渉計を用いることもできる。また、Ｘ線の代わりに光などの電磁波を用いても良いし、電子線を用いても良い。
【００２５】
（第１実施形態）
本実施形態では、軽く重み付けする領域の微分位相値を削除することで、微分位相像の重み付けを行い、これを積分することで位相像を取得する演算装置について説明をする。
【００２６】
本実施形態の演算装置は、図１４に示すトールボット・ロー干渉計により得られたモアレを強度分布として用いた位相回復を行うことにより、微分位相像を取得する手段を有する。この干渉計は、電磁波としてＸ線を用いる。トールボット・ロー干渉計は、トールボット干渉計の一種であり、ロー効果を用いることで、トールボット干渉計よりもコントラストが高い強度分布を形成することができる。トールボット・ロー干渉計については、例えば、Ｐｒｏｃ．ＳＰＩＥ  ６３１８，６３１８０Ｓ（２００６）に詳細が記載されているため、ここでは簡単に説明をする。
【００２７】
図１４は、本実施形態に係る、トールボット・ロー干渉計１の模式図である。トールボット・ロー干渉計は、Ｘ線源１１０からのＸ線を分割する線源格子１２０、線源格子からのＸ線で強度分布を形成する回折格子１４０、回折格子からのＸ線の一部を遮蔽する遮蔽格子１５０、遮蔽格子からのＸ線を検出する検出器１６０を備える。このトールボット・ロー干渉計２００は、線源格子にＸ線を照射するＸ線源１１０と、検出器１６０の検出結果に基づいて被検体の情報を取得する演算装置１７０とともに、位相値取得システム１００を構成することができる。
【００２８】
線源格子１２０はＸ線透過部とＸ線遮蔽部を有しており、Ｘ線源１１０からのＸ線をビーム状に分割する。Ｘ線源から回折格子に照射されるＸ線の可干渉性が、回折格子により干渉縞を形成できる程度に高ければ、線源格子は不要である。線源格子１２０により分割されたＸ線は被検体１３０に照射され、被検体１３０を透過したＸ線が回折格子により回折される。
【００２９】
回折格子１４０はＸ線を回折して自己像と呼ばれる干渉縞を形成する光学素子であり、例えば、Ｘ線の位相を周期的に変化させる位相型の回折格子又は、Ｘ線の振幅を周期的に変化させる振幅型の回折格子を用いることができる。回折格子１４０が被検体１３０を透過したＸ線を回折すると、トールボット距離と呼ばれる所定の距離にＸ線の干渉縞が形成される。
【００３０】
遮蔽格子１５０は、遮蔽格子１５０上に干渉縞が形成されるように、回折格子１４０とトールボット距離はなれた位置に配置される。遮蔽格子１５０はＸ線の遮蔽部とＸ線の透過部とが周期的に配列されている。遮蔽部と透過部とが配列された周期は、遮蔽格子１５０上に形成される干渉縞の周期とわずかに異なるため、遮蔽格子１５０を透過したＸ線は、モアレを形成する。本実施形態では、このモアレを強度分布として用い、微分位相像と位相像を取得する。検出器１６０はＸ線の強度を検出する画素を複数有しており、遮蔽格子１５０からのＸ線を検出することで、モアレの情報を取得する。尚、遮蔽格子の周期と遮蔽格子上に形成される干渉縞の周期とを同一にし、回転させることで回転モアレを形成しても良い。また、遮蔽格子を用いてモアレを形成する代わりに、干渉縞を直接検出器により検出しても良い。その場合、干渉縞の強度分布を取得することが可能な空間分解能を有する検出器を用いればよい。
【００３１】
演算装置１７０は、微分位相像を取得する手段２５０と、重み付け手段２１０と、積分手段２２０を有し、検出器１６０による検出結果を用いて位相像を取得する。図２を用いて、演算装置１７０が行う位相像の取得方法についてより詳細に説明する。図２は、演算装置による位相像の取得方法の概要を示す模式図である。図２（ａ）は検出器１６０により検出されたモアレを示す。ここでは、検出器として５１２×５１２画素のエリアセンサを用い、演算装置が検出器から５１２画素×５１２画素の強度分布１を取得したものとして説明をする。また、ここでは強度分布、微分位相像と位相像は白地で示す。微分位相像を取得する手段２５０は、フーリエ変換法等、複数の画素の強度情報を用いて１画素分の微分位相値を取得する位相回復法を用いて、微分位相像３を取得する工程Ｓ１０を行う（図２（ｂ））。重み付け手段２１０は、取得された微分位相像から、軽く重み付けする領域の微分位相値を削除することで、重み付け工程Ｓ１１を行う。削除方法は特に問わないが、例えば、微分位相像よりも小さな矩形のフィルターを用いても良い。また、微分位相像の座標を指定し、指定した範囲内の微分位相像を、積分工程（Ｓ１２）における重み付けされた微分位相像として指定するプログラムを組んでも良い。
【００３２】
軽く重み付けする領域は、端部のみでも良いし、その周辺の画素を含んでいても良い。強度分布から微分位相像を取得する方法によっては、端部の周辺の画素の微分位相値も、不連続な画素の強度情報を用いて取得されることがある。
【００３３】
例えば、強度分布をフーリエ変換し、取得したフーリエ空間からキャリア周波数のピークを中心とした領域をフィルター関数を用いて切り出し、切り出した領域を空白のフーリエ空間の原点へ移動して逆フーリエ変換を行うことで微分位相像を取得するとする。この場合、フーリエ空間からキャリア周波数のピークを中心とした領域を切り出す際に用いるフィルター関数の大きさによって、不連続な画素の強度情報を用いて取得される画素数が決まる。但し、本発明及び本明細書において、フィルター関数の大きさとは、実空間上における大きさであり、フーリエ空間上でフィルター関数を用いる場合は、そのフィルターを実空間へ変換したときの大きさのことを指す。
【００３４】
ｘ方向（横軸方向）の周期に由来するキャリア周波数のピークを中心として切り出す際、フィルター関数が３画素分の大きさを持つとき、（３−１）／２＝１で、端部から１列分（つまり端部のみ）の微分位相値が不連続な画素の強度情報を用いて取得される。同様に、フィルター関数の大きさが６画素分の大きさであるとき（６−１）／２＝２．５で、端部から３列分の微分位相値が不連続な画素の強度情報を用いて取得される。但し、フィルター関数の大きさとは、フィルターの最大値の１％以上の値となる範囲のことをさし、例えば、ガウス関数であれば、分散をσ２として、６σをフィルター関数の大きさとみなすことができる。一般的に、上述の方法で微分位相像を取得する場合、フィルター関数の大きさは３画素以上である。一方、特開２０１１−１５３９６９に記載されているように、２つの強度分布を用いて０次ピークをキャンセルしてからキャリア周波数のピークを中心とした領域を切り出す場合、フィルター関数の大きさは２画素とすることもできる。この場合、端部から１列分の微分位相値が、不連続な画素の強度情報を用いて取得される。また、“Ｗｉｎｄｏｗｅｄ  Ｆｏｕｒｉｅｒ  ｔｒａｎｓｆｏｒｍ  ｍｅｔｈｏｄ  ｆｏｒ  ｄｅｍｏｄｕｌａｔｉｏｎ  ｏｆ  ｃａｒｒｉｅｒ  ｆｒｉｎｇｅｓ，”  Ｏｐｔ．Ｅｎｇ．４３（７）  １４７２−１４７３（Ｊｕｌｙ  ２００４）に記載されているように、周期的な強度分布（周期パターンと呼ぶことがある）を複数の領域毎にフーリエ変換しても良い。この方法の場合、領域毎にフーリエ変換して取得された領域毎のフーリエスペクトルから、キャリア周波数のピークを中心として切り出す際に用いるフィルター関数の大きさによって、不連続な画素の強度情報を用いて取得される微分位相値の画素数が決まる。
【００３５】
また、特開２０１３−０５０４４１に記載されているように、複数の強度分布のそれぞれから、ある画素とその周辺の画素との強度情報を取り出し、１つの微分位相値を取得するとする。この場合、１つの強度分布から取り出す強度情報の数に応じて不連続な画素の強度情報を用いて取得される画素数が決まる。但し、強度分布毎に取り出す画素数が異なる場合、取り出す画素数の最大値に応じて不連続な画素の強度情報を用いて取得される画素数が決まる。例えば、第１の強度分布から２画素分の強度情報、第２の強度分布から１画素分の強度情報を取り出し、１つの微分位相値を取得する場合、端部から１列分の微分位相値が、不連続な画素の強度情報を用いて取得される。
【００３６】
しかしながら、不連続な画素の強度情報を用いて取得された微分位相値の全てを削除する必要はない。削除する微分位相値の数が多いほど、得られる位相像の画素数が減る。よって、正確性が低い微分位相値による影響をどの程度まで許容できるかと、得たい位相像の画素数とを考慮して削除する画素数を決定することができる。複数の列または行の微分位相値が、不連続な画素の強度情報を用いて取得された場合であっても、原理的には端部に近い画素の微分位相値ほど、取得される微分位相値の正確性が低い。よって、端部のみの微分位相値を削除することは、他の画素の微分位相値を削除することよりも、位相像の画素数減少というデメリットに対する位相像の正確性向上というメリットが大きい。よって、端部のうち１つ以上の画素は、軽く重み付けする領域とする。また、軽く重み付けする領域は、画素列または画素行単位で設定しても良いし、画素単位ごとに設定しても良い。
【００３７】
得たい位相像の画素数は、撮像範囲の大きさと、被検体の観察領域の大きさとの関係だけでなく、被検体の観察領域と撮像範囲における位置との関係を考慮して決めることが好ましい。例えば、被検体の観察領域が撮像範囲に対して小さくても、その観察領域が端部から３画素の領域にも存在するとすれば、軽く重み付けする領域を端部から３画素以内の画素列と画素行とすると、観察領域の情報が失われるため好ましくない。また、例えば、観察領域が撮像範囲中の右側に偏っている場合、微分位相像の右端から削除する画素の行数を、左端から削除する行数よりも小さくしても良い。また、削除する画素の列数と行数は適宜決定することができ、例えば左右５行、上下３列を削除の対象としても良い。つまり、削除する画素の行数及び列数は適宜決定することができ、行数同士が異なっていても良いし、列数同士が異なっていても良い。
【００３８】
また、１方向の微分位相像（位相像をｘ方向またはｙ方向に微分した微分位相像）を積分する場合、積分する方向に垂直な行または列を削除するだけでも良い。尚、本発明及び本明細書では、位相像をｘ方向に微分した微分位相像のことをｘ方向の微分位相像、位相像をｙ方向に微分した微分位相像のことをｙ方向の微分位相像と呼ぶ。たとえば、横方向をｘ方向とし、ｘ方向の微分位相像ｄΦ（ｘ，ｙ）／ｄｘを用いる場合は、積分手段は重み付けされた微分位相像をｘ方向に積分する。この場合、縦方向であるｙ方向には積分を行わないため、上下の端部の正確性の低下の影響は、位相像においても広がらない。よって、軽く重み付けする領域は、微分位相像の端部のうち、ｘ方向に垂直なｙ方向に配列した左右の画素列の画素を１つ以上有する領域とすることが好ましい。一方、２方向の微分位相像（ｘ方向の微分位相像とｙ方向の微分位相像）を積分して１つの位相像を取得する場合は、ｘ方向の微分位相像はｘ方向に、ｙ方向の微分位相像はｙ方向に積分する。よって、軽く重み付けする領域は、微分位相像の端部のうち、ｙ方向に配列した左端と右端の画素列の画素を１つ以上と、ｘ方向に配列した上端と下端の画素行の画素を１つ以上を有する領域とすることが好ましい。
【００３９】
軽く重み付けする領域の設定単位を、画素単位とする場合、例えば、端部のうちで周辺の画素と比較して微分位相値が急激に変化する画素と、そうではない画素とに分け、微分位相値が急激に変化する画素を軽く重み付けする領域とすればよい。この場合、隣の画素との微分位相値の差が任意の値以上の画素を軽く重み付けする領域としても良いし、周辺の画素の微分位相値の平均値との差が任意の値以上の画素を軽く重み付けする領域としても良い。
【００４０】
ここでは、具体例として端部から５画素ずつの画素行及び画素列を軽く重み付けする領域とする場合を説明する。演算装置の重み付け手段は、端部の微分位相値を削除することにより、５０２×５０２画素の重み付けされた微分位相像５（図２（ｃ））を取得する。そして、積分手段が、重み付けされた微分位相像５を積分することで、５０２×５０２画素の位相像７を取得することができる（図２（ｄ））。重み付けされた微分位相像の積分方法は特に問わない。例えば、単純に、右端又は左端から足し合わせていくことで、横軸方向（ｘ軸方向）への積分を行っても良いし、特開２０１３−１０２９５１に記載されているような、微分位相像を再度微分してからポワソン方程式を解き、積分を行っても良い。また、“Ｎｏｎｉｔｅｒａｔｉｖｅ  ｂｏｕｎｄａｒｙ−ａｒｔｉｆａｃｔ−ｆｒｅｅ  ｗａｖｅｆｒｏｎｔ  ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ  ｆｒｏｍ  ｉｔｓ  ｄｅｒｉｖａｒｉｖｅｓ，”Ａｐｐｌｉｅｄ  Ｏｐｔｉｃｓ，Ｖｏｌ．５１，Ｎｏ．２３，５６９８−５７０４（２０１２）に記載されているような積分を行っても良い。
【００４１】
得られた位相像は、演算装置と接続された画像表示装置（不図示）により表示されても良い。画像表示装置は、各種ディスプレイ又はプリンター等を用いることができる。また、位相像を画像としてではなく、位置情報が表示された位相値の一覧として表示しても良いし、位相像に対して更に演算を加え、その結果を表示しても良い。また、微分位相像と位相像の両方を表示する場合、重み付け前の微分位相像の方が画素数が多いため、微分位相像としては重み付けされる前の微分位相像を表示しても良い。
【００４２】
また、演算装置１７０は、複数の演算装置で構成されていても良い。たとえば、微分位相像の取得までを行う演算装置と、微分位相像の重み付けを行う演算装置と、重み付けされた微分位相像を積分して位相像を取得す演算装置の計３つの演算装置が演算装置１７０として機能しても良い。このように、演算装置１７０が複数の演算装置として機能するとき、演算装置のそれぞれは取得したデータの送受信ができれば、物理的に近接して配置する必要はない。
【００４３】
（第２実施形態）
本実施形態では、軽く重み付けする領域の微分位相値を別の値に置き換えることで、微分位相像の重み付けを行う演算装置について説明をする。  本実施形態における演算装置は、重み付け手段による微分位相像の重み付け方法以外は実施形態１と同様のため、説明は省略する。
【００４４】
演算装置による位相像の取得方法の概要を示す模式図を図３に示す。実施形態１と同様に、微分位相像を取得する手段は、強度分布１（図３（ａ）、５１２×５１２画素）を用いて、微分位相像３を取得する（図３（ｂ）、５１２×５１２画素）。
【００４５】
重み付け手段は、取得した微分位相像のうち、軽く重み付けする領域の微分位相値を別の値に置き換え、重み付けされた微分位相像を取得する。置き換える値は、０に近い値であることが好ましいが、置き換える前の微分位相値よりも０に近い（つまり、絶対値が小さい）値であれば本実施形態の効果を得ることができる。例えば、重み付け前の微分位相像から微分位相像全体の微分位相値の平均値を取得し、軽く重み付けする領域の微分位相値をその平均値に置き換えても良い。微分位相値の平均値は、軽く重み付けする領域の置き換え前の値よりも０に近いと推定されるためである。このように、０以外の値に置き換える場合、局所的にみると、置き換えにより微分位相値の絶対値が大きくなる領域が存在することがある。しかしながら、軽く重み付けする領域全体の、置き換え後の微分位相値の絶対値の平均値が、置き換え前の微分位相値の絶対値の平均値より小さければ、本実施形態の効果を得ることができる。言い換えると、本実施形態の重み付け手段は、重み付けされた微分位相像のうち軽く重み付けする領域の微分位相値の絶対値の平均値を、重み付け手段に入力された微分位相像のうち軽く重み付けする領域の微分位相値の絶対値の平均値より小さくする。重み付けされた微分位相像のうち軽く重み付けする領域の微分位相値の絶対値の平均値が、重み付け手段に入力された微分位相像のうち軽く重み付けする領域の微分位相値の絶対値の平均値の１／２以下になるように重み付けを行うことがより好ましい。重み付けされた微分位相像のうち軽く重み付けする領域の微分位相値の絶対値の平均値が、重み付け手段に入力された微分位相像のうち軽く重み付けする領域の微分位相値の絶対値の平均値の１／１０以下になるように重み付けを行うことが更に好ましい。
【００４６】
微分位相値を置き換えると、置き換えた範囲は被検体の正確な情報を持たない。よって、軽く重み付けする領域は、実施形態１の軽く重み付けする領域と同様に決定することができる。また、実施形態１と同様に、軽く重み付けする領域の設定単位は、列または行毎ででも良いし、画素毎でも良い。画素毎に軽く重み付けする領域を設定する場合も、実施形態１と同様に軽く重み付けする領域を決定することができる。また、実施形態１と同様に、軽く重み付けする領域は、左右対称でなくても良いし、上下対称でなくても良い。また、微分位相像と位相像の両方を表示する場合、微分位相像としては重み付けされる前の微分位相像を表示しても良い。
【００４７】
ここでは、端部から５画素ずつの画素行及び画素列を軽く重み付けする領域とし、軽く重み付けする領域の微分位相値を０に置き換える場合を例に説明をする。重み付け手段は、端部をいれて上下左右５画素ずつの行または列の微分位相値を０に置き換える。これにより、重み付けされた微分位相像１０（図３（ｃ））を取得する。重み付けされた微分位相像１０は、微分位相値が０の、５１２×５１２画素の領域９の中心と、５０２×５０２画素の微分位相像５の中心とが一致するように、領域９上に微分位相像５が配置された微分位相像とみなすこともできる。そして、重み付けされた微分位相像を積分することで、５１２×５１２画素の位相像１１を取得することができる（図３（ｄ））。
【００４８】
（第３実施形態）
本実施形態では、軽く重み付けする領域の微分位相値の絶対値が小さくなるようなフィルターをかけることで、微分位相像の重み付けを行う演算装置について説明をする。本実施形態における演算装置は、重み付け手段による微分位相像の重み付け方法以外は実施形態１と同様のため、説明は省略する。
【００４９】
微分位相像を取得する手段は、実施形態１と同様に、強度分布を用いて、微分位相像を取得する。
【００５０】
重み付け手段は、取得した微分位相像にフィルターをかけることで、中央部の微分位相値に対する軽く重み付けする領域の微分位相値を小さくする。言い換えると、軽く重み付けする領域の微分位相値の平均値／中央部の微分位相値（但し、複数ある場合はその平均値とする）を小さくする。フィルターは、中央部の微分位相値に対する軽く重み付けする領域の微分位相値が、フィルターをかける前よりも小さくなるフィルターであれば特に形状は問わない。例えば、軽く重み付けする領域の微分位相値の絶対値を均一に小さくする（例えば、軽く重み付けする領域の微分位相値の絶対値に対して１以下の所定の値をかけるような）フィルターを用いることができる。軽く重み付けする領域が端部の周辺の画素を含む場合、軽く重み付けする領域のなかでも端部に近いほど絶対値が小さくなるフィルターであることが好ましい。また、中央部の微分位相値に対する軽く重み付けする領域の微分位相値が小さくなればよいため、例えば、軽く重み付けする領域の微分位相値をそのままの値とし、中央部の微分位相値の絶対値が大きくなるフィルターでも良い。中央部の微分位相値に掛けられるフィルターの値（但し、中央部が複数ある場合はその平均値とする）を１としたときに、軽く重み付けする領域の微分位相値に掛けられるフィルターの値の平均値が０．５以下であることが好ましい。更に、連続する画素の強度情報のみを用いて取得された微分位相値に掛けられるフィルターの値の平均値を１としたときに、軽く重み付けする領域の微分位相値に掛けられるフィルターの値の平均値が０．５以下であることがより好ましい。尚、中央部でも軽く重み付けする領域でもない画素の微分位相値には、中央部の微分位相値にかけられるフィルターの値と同じ値が掛けられることが好ましい。また、フィルターは、微分位相像全体にかけても良いし、軽く重み付けする領域だけにかけても良い。
【００５１】
軽く重み付けする領域は、実施形態１と同様に決定することができる。また、実施形態１と同様に、軽く重み付けする領域の設定は、列または行単位でも良いし、画素単位でも良い。画素単位で軽く重み付けする領域とする場合も、実施形態１と同様に軽く重み付けする領域を決定することができる。また、実施形態１と同様に、軽く重み付けする領域は、左右対称でなくても良いし、上下対称でなくても良い。また、微分位相像と位相像の両方を表示する場合、微分位相像としては重み付けされる前の微分位相像を表示しても良い。
【００５２】
端部から５画素ずつの画素行及び画素列を軽く重み付けする領域とする場合、本実施形態は、図３（ｃ）における領域９の微分位相値が０ではなく、フィルターにより絶対値が小さくなった微分位相値である点が実施形態２（図３）と異なる。しかしながら、５１２×５１２画素の微分位相像を用いて５１２×５１２画素の位相像を取得する点は実施形態２と同様である。
【００５３】
（第４実施形態）
本実施形態では、微分位相像を取得する際に、端部のうち１つ以上の画素を有する領域の微分位相値を取得しない演算装置について説明をする。
【００５４】
本実施形態における演算装置は、実施形態１の演算装置と同様に、微分位相像を取得する手段２４０を有する。但し、本実施形態の微分位相像を取得する手段は、微分位相像を取得する際に、端部のうち１つ以上の画素を有する領域の微分位相値を取得しない点が実施形態１の微分位相像を取得する手段と異なる。言い換えると、強度分布（図２（ａ））から直接（撮像範囲全体の微分位相像（図２（２（ｂ））を経由することなく）、実施形態１で取得される重み付けされた微分位相像（図２（ｃ））を取得する。つまり、本実施形態の位相像を取得する手段２４０は、実施形態１の微分位相像を取得する手段２５０と重み付け手段２１０の両方の機能を有する。微分位相像を取得する方法としては、“Ｗｉｎｄｏｗｅｄ  Ｆｏｕｒｉｅｒ  ｔｒａｎｓｆｏｒｍ  ｍｅｔｈｏｄ  ｆｏｒ  ｄｅｍｏｄｕｌａｔｉｏｎ  ｏｆ  ｃａｒｒｉｅｒ  ｆｒｉｎｇｅｓ，”  Ｏｐｔ．Ｅｎｇ．４３（７）  １４７２−１４７３（Ｊｕｌｙ  ２００４）に記載されている方法や特開２０１３−０５０４４１に記載されている方法のように、周期パターンの領域毎に微分位相像を取得する方法であれば用いることができる。
【００５５】
例として、端部から５画素ずつの画素行及び画素列を微分位相値を取得しない領域とし、この領域の微分位相値を取得しない場合の説明をする。演算装置の微分位相像を取得する手段２４０は、検出器から５１２×５１２画素の強度分布（図２（ａ）と同じ）を受け取り、位相回復を行うことで微分位相像を取得する。この際、端部から５画素ずつの画素行及び画素列の微分位相値は取得しないため、取得される微分位相像は５０２×５０２画素であり、実施形態１で取得した重み付けされた微分位相像（図２（ｃ））と等しい。よって、本発明及び本明細書では、このように、検出器から取得した強度分布の端部に対応する微分位相値を取得せず、強度分布の画素数よりも小さな画素数の微分位相像を取得することも、端部の微分位相値を中央部の微分位相値よりも軽く重み付けすると呼ぶ。
【００５６】
積分手段２２０は、実施形態１の積分手段と同様に、重み付けされた微分位相像を積分することで、位相像を取得する。本実施形態により、実施形態１と同様の位相像を取得することができる。また、同じ位相回復方法を用いる場合は、本実施形態の方が実施形態１〜３よりも重み付けされた微分位相像の取得にかかる時間が短い。しかしながら、実施形態１〜３で取得される微分位相像よりも本実施形態で取得される重み付けされた微分位相像のほうが、画素数が小さい。よって、例えば、微分位相像を観察したい場合など、微分位相像自体の情報が必要であり、正確性が低くても本実施形態の微分位相値を取得しない領域に相当する領域の微分位相値を取得したい場合は実施形態１〜３を行うことが好ましい。
【００５７】
本実施形態も、微分位相値を取得しない領域は、実施形態１の軽く重み付けする領域と同様に決定することができる。また、実施形態１と同様に、微分位相値を取得しない領域の設定は、画素行または画素列単位でも良いし、画素単位でも良い。画素単位で設定する場合、例えば強度分布における被検体の位置に応じて微分位相値を取得しない画素を決めることができる。また、実施形態１と同様に、微分位相値を取得しない領域は、左右対称でなくても良いし、上下対称でなくても良い。
【００５８】
尚、演算装置は、実施形態１〜４のうち、複数の位相像取得方法を行うことが可能であっても良い。その場合、演算装置が行うことができる位相像取得方法をユーザーが適宜切り替えることができるように構成されていても良いし、微分位相像に応じて、演算装置が適宜位相像取得方法を選択しても良い。また、複数の位相像取得方法を行い、位相像取得方法が異なる複数の位相像を取得しても良い。この場合、複数の位相像を表示しても良い。
【実施例】
【００５９】
以下、実施例及び比較例について説明をする。
【００６０】
（比較例１）
本比較例では、“Ｎｏｎｉｔｅｒａｔｉｖｅ  ｂｏｕｎｄａｒｙ−ａｒｔｉｆａｃｔ−ｆｒｅｅ  ｗａｖｅｆｒｏｎｔ  ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ  ｆｒｏｍ  ｉｔｓ  ｄｅｒｉｖａｒｉｖｅｓ，”Ａｐｐｌｉｅｄ  Ｏｐｔｉｃｓ，Ｖｏｌ．５１，Ｎｏ．２３，５６９８−５７０４（２０１２）に記載されている積分法により位相像を取得する例を、シミュレーションを用いて説明する。
【００６１】
シミュレーションで設定した、被検体の影響がない強度分布ｇｒ（ｘ，ｙ）と被検体の影響を含む強度分布ｇｓ（ｘ，ｙ）を下記式に示す。ｇｒ（ｘ，ｙ）＝ａ（ｘ，ｙ）×（１＋ｂ（ｘ，ｙ）ｃｏｓ（２πｆ０ｘ））×（１＋ｂ（ｘ，ｙ）ｃｏｓ（２πｆ０ｙ））
ｇｓ（ｘ，ｙ）＝ａ（ｘ，ｙ）×（１＋ｂ（ｘ，ｙ）ｃｏｓ（２π（ｆ０ｘ＋ｐｘ（ｘ，ｙ））））×（１＋ｂ（ｘ，ｙ）ｃｏｓ（２π（ｆ０ｙ＋ｐｙ（ｘ，ｙ））））
ｆ０＝１／ｐｄ／ｐｍ
ｐｘ（ｘ，ｙ）＝α×（ｄΦ（ｘ，ｙ）／ｄｘ）
ｐｙ（ｘ，ｙ）＝α×（ｄΦ（ｘ，ｙ）／ｄｙ）
【００６２】
ただし、ｘとｙは画素単位で位置を表わす座標であり、ｐｄは画素の大きさを示し、本比較例では１０μｍとした。また、ｐｍは周期パターンの周期を画素単位で示したもので、本比較例では８画素とし、αは定数であり、本比較例では７２５００とした。また、計算領域は５００×５００画素、すなわち５×５ｍｍとした。また、ａ（ｘ，ｙ）、ｂ（ｘ，ｙ）は一定値とし、位相像Φ（ｘ，ｙ）は直径２．５ｍｍの球状の被検体が計算領域の端に位置するよう設定した。
【００６３】
シミュレーションに用いた球の位相像Φ（ｘ，ｙ）を図４（ｇ）に、図４（ｇ）に示した点線１２のラインプロファイルを図４（ｈ）に示した。また、ｄΦ（ｘ，ｙ）／ｄｘを図４（ａ）に、図４（ａ）に示した点線の枠ａ１内の拡大図を図４（ｃ）に、点線の枠ａ２内の拡大図を図４（ｄ）に示した。また、ｄΦ（ｘ，ｙ）／ｄｙを図４（ｂ）に、図４（ｂ）に示した点線の枠ｂ１内の拡大図を図４（ｅ）に、点線の枠ｂ２内の拡大図を図４（ｆ）に示した。また図５（ａ）にｇｒ（ｘ，ｙ）、図５（ｂ）にｇｓ（ｘ，ｙ）を示した。つまり、図４（（ａ）〜（ｈ））と図５（（ａ）、（ｂ））が、本比較例における入力値である。
【００６４】
図５（ａ）、（ｂ）で表わされた周期パターンを用いて、フーリエ変換法により取得した微分位相像を取得した。つまり、図５（ａ）をフーリエ変換し、フィルター関数を用いて、フーリエ空間上でｘ方向におけるキャリア周波数のスペクトルを切り出した。フィルター関数は、実空間で設定したガウス関数のフーリエスペクトルとし、ガウス関数の分散の平方根（σ）を２．５画素とした。切り出したスペクトルを別のフーリエ空間の原点へ移動し、これを逆フーリエ変換して、ｄΦｒ（ｘ，ｙ）／ｄｘを取得した。図５（ｂ）も同様にフーリエ変換から逆フーリエ変換までの工程を行い、ｄΦｓ（ｘ，ｙ）／ｄｘを取得した。ｄΦｓ（ｘ，ｙ）／ｄｘとｄΦｒ（ｘ，ｙ）／ｄｘの差分をとることで、ｘ方向の微分位相像ｄΦ（ｘ，ｙ）／ｄｘを取得した。同様に、フーリエ空間上でｙ方向におけるキャリア周波数のスペクトルを切り出し、ｄΦｓ（ｘ，ｙ）／ｄｙとｄΦｒ（ｘ，ｙ）／ｄｙの差分をとることで、ｙ方向の微分位相像ｄΦ（ｘ，ｙ）／ｄｙを取得した。ｄΦ（ｘ，ｙ）／ｄｘおよびｄΦ（ｘ，ｙ）／ｄｙを図６（ａ）、（ｂ）にそれぞれ示した。また、図４同様に、図６（ｃ）、（ｄ）に、図６（ａ）の点線の枠ａ１、ａ２内の拡大図をそれぞれ示した。また、図６（ｅ）、（ｆ）に、図６（ｂ）の点線の枠ｂ１、ｂ２内の拡大図をそれぞれ示した。図６（ａ）、（ｂ）に示した微分位相像の微分位相値のそれぞれは、６σ＝１５画素分の強度情報を用いて取得され、端部から、（１５−１）／２＝７行及び列の画素の微分位相値が不連続な画素の強度情報を用いて取得された。
【００６５】
ｄΦ（ｘ，ｙ）／ｄｘおよびｄΦ（ｘ，ｙ）／ｄｙを用い、“Ｎｏｎｉｔｅｒａｔｉｖｅ  ｂｏｕｎｄａｒｙ−ａｒｔｉｆａｃｔ−ｆｒｅｅ  ｗａｖｅｆｒｏｎｔ  ｒｅｃｏｎｓｔｒｕｃｔｉｏｎ  ｆｒｏｍ  ｉｔｓ  ｄｅｒｉｖａｒｉｖｅｓ，”Ａｐｐｌｉｅｄ  Ｏｐｔｉｃｓ，Ｖｏｌ．５１，Ｎｏ．２３，５６９８−５７０４（２０１２）に記載されている積分法を用いてΦ（ｘ，ｙ）を取得した。Φ（ｘ，ｙ）を図１２（ａ）に示し、図１２（ａ）に示した点線１２のラインプロファイルを図１２（ｂ）に示した。
【００６６】
図４（ｃ）と図６（ｃ）に示された、被検体の球の左端部の微分位相値ｄΦ（ｘ，ｙ）／ｄｘを比べると、本比較例で取得された微分位相像は入力値に基づく像（図４（ｃ））に対してボケていることが分かる。一方、図４（ｄ）と図６（ｄ）に示された球の右端部の微分位相値を比べると、本比較例で取得された微分位相像は入力値に基づく像（図４（ｄ））に対してボケているだけでなく、値が大きく異なっていることが分かる。これは、球の右端部が周期パターンの端部にあり、周期パターンの左端部の画素を用いて微分位相値が取得されたため、球の右端部の微分位相像の正確性が下がったものと理解できる。同様に、図４（ｅ）、（ｆ）と図６（ｅ）、（ｆ）に示したｄΦ（ｘ，ｙ）／ｄｙも、ｄΦ（ｘ，ｙ）／ｄｘのときと同様に、球の下端部では本比較例で取得された微分位相像（図６（ｆ））が入力値に基づく像（図４（ｆ））と大きく異なっていることが分かる。また、図１２（ａ）、（ｂ）を見ると、端部とその周辺だけでなく、点線１２のライン全体に渡って、入力値と取得された位相値とがずれており、端部とその周辺だけでなく、位相像全体における位相値が入力値とずれていることが分かる。
【００６７】
（比較例２）
本比較例では、周期パターンに、周期パターンの端部に向かうにつれ強度が滑らかに０になるようなフィルターを掛けてから、フーリエ変換法で解析したシミュレーション結果について説明をする。
【００６８】
図１１（ａ）に、本比較例で用いたフィルターとその拡大図を示した。また、図１１（ａ）の点線で囲まれた範囲の拡大図を図１１（ｂ）に示した。図１１（ａ）、（ｂ）に示したフィルターは、端部が０となり、端部を含めて６画素目で１となるｓｉｎ曲線を有する。つまり、端部を含めて５画素が重み付けされる。また、図１１（ｃ）に、被検体ありの周期パターンｇｓ（ｘ，ｙ）に、図１１（ａ）に示したフィルターを掛けて取得される周期パターンを示した。また、図１１（ｃ）の点線で囲まれた範囲の拡大図を図１１（ｄ）に示した。フィルターを掛けた周期パターンは、その端部（図１１中の右端と下端）で強度が滑らかにゼロになった。図１１（ｃ）に示した周期パターンを用いて、比較例１と同様に微分位相像をした。取得したｄΦ（ｘ，ｙ）／ｄｘおよびｄΦ（ｘ，ｙ）／ｄｙを、図１３（ａ）、（ｂ）にそれぞれ示した。また、図６同様に、図１３（ｃ）、（ｄ）に、図１３（ａ）の点線の枠ａ１、ａ２内の拡大図をそれぞれ示した。また、図１３（ｅ）、（ｆ）に、図１３（ｂ）の点線の枠ｂ１、ｂ２内の拡大図をそれぞれ示した。また、本比較例で取得したｄΦ（ｘ，ｙ）／ｄｘおよびｄΦ（ｘ，ｙ）／ｄｙを用い、比較例１と同じ方法を用いてΦ（ｘ，ｙ）を取得した。Φ（ｘ，ｙ）を図１３（ｇ）に示し、図１３（ｇ）に示した点線１２のラインプロファイルを図１３（ｈ）に示した。図１３（ｈ）において、実線が入力値（図４（ｈ）と同一）、太線が本比較例で取得したΦ（ｘ，ｙ）のラインプロファイル（取得値と呼ぶ）である。図１２（ａ）、（ｂ）と図１３（ｄ）、（ｆ）を比べると、周期パターンにフィルターを掛けることで、微分位相値の入力値とのズレが比較例１よりも緩和していることが分かり、周期パターンにフィルターを掛けることの効果が確認できた。
【００６９】
（実施例１）
実施例１では、実施形態１のより具体的な例を、シミュレーション結果を用いて説明をする。本実施例では、微分位相像のうち、軽く重み付けする領域の微分位相値を削除する例を説明する。
【００７０】
本実施例では、図６（ａ）、（ｂ）に示した微分位相像のうち、上下左右それぞれ５行もしくは５列を軽く重み付けする領域とし、この領域の微分位相値を削除することで、重み付けされた微分位相像を取得した。不連続な画素の強度情報を用いて取得された微分位相値の一部は中央部の微分位相値よりも軽く重み付けされずに残るが、残った微分位相値は強度情報を用いる１５画素のうち、１３画素以上が連続している。よって、微分位相値の正確性に与える影響は小さいと考えられる。
【００７１】
図６（ａ）、（ｂ）に示した５００×５００画素の微分位相像の上下左右を削除し、４９０×４９０画素の微分位相像を取得した。本実施例により取得される、重み付けされた微分位相像を図７の（ａ）、（ｂ）に示した。また、図６同様に、図７（ｃ）、（ｄ）に、図７（ａ）の点線の枠ａ１、ａ２内の拡大図をそれぞれ示した。また、図７（ｅ）、（ｆ）に、図７（ｂ）の点線の枠ｂ１、ｂ２内の拡大図をそれぞれ示した。本実施例で取得した重み付けされた微分位相像を用い、比較例１、２と同じ方法を用いて位相像Φ（ｘ，ｙ）を取得した。位相像を図７（ｇ）に示し、図７（ｇ）に示した点線１２のラインプロファイルを図７（ｈ）に示した。図７（ｈ）においても、実線が入力値、太線が本実施例の取得値である。
【００７２】
比較例１、２で取得した位相像とそのラインプロファイルと、本実施例で取得した位相像とそのラインプロファイルを比較すると、本実施例により、比較例１、２よりも入力値に近い位相像を取得できたことが分かった。
【００７３】
本実施例では、微分位相像の端部を行毎もしは列毎に強度情報を削除したが、実施形態１で説明をしたように、画素単位で微分位相値を削除しても良い。その場合、例えば、図７（ｄ）で拡大表示された範囲のみの右端を削除して重み付けされたｄΦ（ｘ，ｙ）／ｄｘを取得し、図７（ｆ）で拡大表示された範囲のみの下端を削除して重み付けされたｄΦ（ｘ，ｙ）／ｄｙを取得してもよい。
【００７４】
（実施例２）
実施例２では、実施形態２のより具体的な例を、シミュレーション結果を用いて説明をする。本実施例では、軽く重み付けする領域の微分位相値を０に置き換える例を説明する。
【００７５】
本実施例では、図６（ａ）、（ｂ）に示した微分位相像のうち、上下左右それぞれ５行もしくは５列を軽く重み付けする領域とし、この領域の微分位相値を０に置き換えることで、重み付けされた微分位相像を取得した。本実施例により取得される、重み付けされた微分位相像を図８の（ａ）、（ｂ）に示した。また、図６同様に、図８（ｃ）、（ｄ）に、図８（ａ）の点線の枠ａ１、ａ２内の拡大図をそれぞれ示した。また、図８（ｅ）、（ｆ）に、図８（ｂ）の点線の枠ｂ１、ｂ２内の拡大図をそれぞれ示した。本実施例で取得した重み付けされた微分位相像を用い、実施例１と同じ方法を用いて位相像Φ（ｘ，ｙ）を取得した。位相像を図８（ｇ）に示し、図８（ｇ）に示した点線１２のラインプロファイルを図８（ｈ）に示した。図８（ｈ）においても、実線が入力値、太線が本実施例の取得値である。
【００７６】
比較例１、２で取得した位相像とそのラインプロファイルと、本実施例で取得した位相像とそのラインプロファイルを比較すると、本実施例により、比較例１、２よりも入力値に近い位相像を取得できたことが分かった。
【００７７】
本実施例も、実施例１と同様に、画素単位で微分位相値を置き換えても良い。
【００７８】
（実施例３）
実施例３では、実施形態３のより具体的な例を、シミュレーション結果を用いて説明をする。本実施例では、微分位相像全体にフィルターをかける例を説明する。
【００７９】
本実施例では、図６（ａ）、（ｂ）に示した微分位相像全体に、図１１（ａ）のフィルターをかけることで、重み付けされた微分位相像を取得した。軽く重み付けする領域は、上下左右それぞれ５行もしくは５列とした。本実施例により取得される、重み付けされた微分位相像を図９の（ａ）、（ｂ）に示した。また、図６同様に、図９（ｃ）、（ｄ）に、図９（ａ）の点線の枠ａ１、ａ２内の拡大図をそれぞれ示した。また、図９（ｅ）、（ｆ）に、図９（ｂ）の点線の枠ｂ１、ｂ２内の拡大図をそれぞれ示した。本実施例で取得した重み付けされた微分位相像を用い、実施例１と同じ方法を用いて位相像Φ（ｘ，ｙ）を取得した。位相像を図９（ｇ）に示し、図９（ｇ）に示した点線１２のラインプロファイルを図９（ｈ）に示した。図９（ｈ）においても、実線が入力値、太線が本実施例の取得値である。
【００８０】
比較例１、２で取得した位相像とそのラインプロファイルと、本実施例で取得した位相像とそのラインプロファイルを比較すると、本実施例により、比較例１、２よりも入力値に近い位相像を取得できたことが分かった。
【００８１】
本実施例も、実施例１と同様に、画素単位で微分位相値を置き換えても良い。
【００８２】
（実施例４）
実施例４では、実施形態３のより具体的な例を、シミュレーション結果を用いて説明をする。本実施例は、周期パターンにフィルターをかけてから微分位相像を取得する点が実施例３と異なる。言い換えると、比較例１に実施形態３を適用した例が実施例３であり、比較例２に実施形態３を適用した例が実施例４である。
【００８３】
本実施例では、図１３（ａ）、（ｂ）に示した微分位相像全体に、図１１（ａ）のフィルターをかけることで、重み付けされた微分位相像を取得した。軽く重み付けする領域は、上下左右それぞれ５行もしくは５列とした。本実施例により取得される、重み付けされた微分位相像を図１０の（ａ）、（ｂ）に示した。また、図６同様に、図１０（ｃ）、（ｄ）に、図１０（ａ）の点線の枠ａ１、ａ２内の拡大図をそれぞれ示した。また、図１０（ｅ）、（ｆ）に、図１０（ｂ）の点線の枠ｂ１、ｂ２内の拡大図をそれぞれ示した。本実施例で取得した重み付けされた微分位相像を用い、実施例１と同じ方法を用いて位相像Φ（ｘ，ｙ）を取得した。位相像を図１０（ｇ）に示し、図９（ｇ）に示した点線１２のラインプロファイルを図１０（ｈ）に示した。図１０（ｈ）においても、実線が入力値、太線が本実施例の取得値である。
【００８４】
比較例１、２で取得した位相像とそのラインプロファイルと、本実施例で取得した位相像とそのラインプロファイルを比較すると、本実施例により、比較例１、２よりも入力値に近い位相像を取得できたことが分かった。このように、実施形態３は周期パターンにフィルターをかけてから微分位相像を取得することで、入力値（本来の値）に近い微分位相値を取得する技術と組み合わせることができる。尚、実施形態１、２も同様に、この技術と組み合わせることができる。
【００８５】
（その他の実施例）
本発明は、次に説明する処理でも実行可能である。その処理は、上述の実施形態の１以上の機能を実現するプログラムを、ネットワークまたは記憶媒体を介してシステム又は演算装置に供給し、そのシステム又は演算装置のコンピュータにおける１つ以上のプロセッサーがプログラムを読み出し実行する処理である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【００８６】
以上、本発明の好ましい実施形態について説明したが、本発明はこれらの実施形態に限定されず、その要旨の範囲内で種々の変形および変更が可能である。
【符号の説明】
【００８７】
１７０  演算装置
２１０  重み付け手段
２２０  積分手段
２３０  演算装置
２４０  微分位相像を取得する手段
【図１】
2015213574.tif 000003
【図２】
2015213574.tif 000004
【図３】
2015213574.tif 000005
【図１４】
2015213574.tif 000006
【図４】
2015213574.tif 000007
【図５】
2015213574.tif 000008
【図６】
2015213574.tif 000009
【図７】
2015213574.tif 000010
【図８】
2015213574.tif 000011
【図９】
2015213574.tif 000012
【図１０】
2015213574.tif 000013
【図１１】
2015213574.tif 000014
【図１２】
2015213574.tif 000015
【図１３】
2015213574.tif 000016