(57)【要約】【課題】より計算時間を短縮することのできるホログラムデータ作成方法及び、ホログラムデータ作成プログラムを提供する。【解決手段】本発明の一観点に係るホログラムデータ作成プログラムは、コンピュータに、ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、ウェーブレット領域データのうち一部を選択するステップ、選択された一部の前記ウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を実行させるためのものである。【選択図】図６
2019015935.tif 000002
【特許請求の範囲】
【請求項１】
コンピュータに、
ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、
前記ウェーブレット領域データのうち一部を選択するステップ、
選択された一部の前記ウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を実行させるためのホログラムデータ作成プログラム。
【請求項２】
前記ウェーブレット領域データを得るステップ、及び、前記ウェーブレット領域データのうち一部を選択するステップを複数繰り返し、
それぞれにおいて、選択された一部の前記ウェーブレット領域データの足し合わせを行うステップ、を備える請求項１記載のホログラムデータ作成プログラム。
【請求項３】
ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、
前記ウェーブレット領域データのうち一部を選択するステップ、
選択された一部の前記ウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を備えるホログラムデータ作成方法。
【請求項４】
前記ウェーブレット領域データを得るステップ、及び、前記ウェーブレット領域データのうち一部を選択するステップを複数繰り返し、
それぞれにおいて、選択された一部の前記ウェーブレット領域データの足し合わせを行うステップ、を備える請求項３記載のホログラムデータ作成方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ホログラムデータ作成プログラム及びホログラムデータ作成方法に関する。
【背景技術】
【０００２】
ホログラフィ技術によって作製されるホログラムは三次元ディスプレイ、特殊な光学素子、大容量メモリ、セキュリティといった様々な分野で利用される極めて有用なものである。
【０００３】
一方で、ホログラムを作製するホログラフィ技術では、三次元物体や画像からの光波を計算するが、この光波の計算にはゾンマーフェルト回折積分、フレネル回折積分等が用いられる。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、上記積分計算では、高速フーリエ変換を用いる方法が考えられる。
【０００５】
例えば、伝搬元、伝搬先ともにＮ×Ｎ画素の面間の回折計算を積分で計算する場合を考えると、その計算量はＮ４に比例して増加する。例えばＮ＝２０４８の場合、一般的なコンピュータを用いた場合でも１時間程度必要となる。これに対し、高速フーリエ変換を用いた場合、その計算量はＮ２ｌｏｇＮまで減らすことができる。すなわち、上記Ｎ＝２０４８の場合、１秒程度にまで短くすることができる。
【０００６】
しかしながら、高速フーリエ変換による光波の計算は平面間の光波伝搬しか計算できず、平面間に限定される結果、三次元的な光波の分布を求めるには高速フーリエ変換を複数回使用する必要があり、上記の時間でも十分短いとは言えない。
【０００７】
そこで、本発明は、上記課題に鑑み、より計算時間を短縮することのできるホログラムデータ作成方法及び、ホログラムデータ作成プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決する本発明の一観点に係るホログラムデータ作成プログラムは、コンピュータに、ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、ウェーブレット領域データのうち一部を選択するステップ、選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を実行させるためのものである。
【０００９】
また、本発明の他の一観点に係るホログラムデータ作成方法は、ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、ウェーブレット領域データのうち一部を選択するステップ、選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を備えるものである。
【発明の効果】
【００１０】
以上、本発明によって、より計算時間を短縮することのできるホログラムデータ作成方法及び、ホログラムデータ作成プログラムを提供することができる。
【図面の簡単な説明】
【００１１】
【図１】ウェーブレット係数の分布の例を示す図である。
【図２】ＰＳＦの例を示す図である。
【図３】ＰＳＦのウェーブレット分解の計算結果の例を示す図である。
【図４】ＰＳＦの例を示す図である。
【図５】図４の例における強度プロファイルを示す図である。
【図６】実施例にて作成したホログラムの例を示す図である。
【図７】実施例にて作成したホログラムの例を示す図である。
【図８】実施例にて作成したホログラムの例を示す図である。
【発明を実施するための最良の形態】
【００１２】
以下、本発明の実施形態について図面を用いて詳細に説明する。ただし、本発明は多くの異なる形態による実施が可能であり、以下に示す実施形態、実施例における具体的な例示にのみ限定されるわけではない。
【００１３】
本実施形態に係るホログラムデータ作成プログラム（以下「本プログラム」という。）は、
コンピュータに、（Ｓ１）ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、（Ｓ２）ウェーブレット領域データのうち一部を選択するステップ、（Ｓ３）上記（Ｓ１）及び（Ｓ２）を複数回繰り返し、それぞれにおいて、選択された一部のウェーブレット領域データの足し合わせを行うステップ、（Ｓ４）選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を実行させるためのものである。
【００１４】
限定されるわけではないが、本プログラムは、情報処理装置、いわゆるコンピュータのハードディスクやＣＤ−ＲＯＭ等の記録媒体に記録され、実行されることにより実現できる。より具体的に説明すると、中央演算装置（ＣＰＵ）、ハードディスクやフラッシュメモリ等の不揮発性の記憶装置、ＲＡＭ等の揮発性の記憶装置、これらを電気的に接続するバスライン、キーボードやマウス等の入力装置、液晶ディスプレイ等の表示装置等を備えたコンピュータを用い、上記不揮発性の記憶装置に本プログラムを格納し、揮発性の記憶装置に読み込み実行させることで所望の処理を実行することができる。
【００１５】
すなわち、上記の記載から明らかなとおり、上記プログラムを実行することで、（Ｓ１）ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、（Ｓ２）ウェーブレット領域データのうち一部を選択するステップ、（Ｓ３）上記（Ｓ１）及び（Ｓ２）を複数回繰り返し、それぞれにおいて、選択された一部のウェーブレット領域データの足し合わせを行うステップ、（Ｓ４）選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を実行するホログラムデータ作成方法（以下「本方法」という。）を提供することができる。以下、本方法の詳細について説明する。
【００１６】
光波の計算は、原理的には点光源から発する光波の足し合わせと考えることができ、本方法はこの光波の足し合わせをウェーブレット変換領域において行うことで、光波計算の計算量を大幅に削減することができる。
【００１７】
まず、本方法は、（Ｓ１）ウェーブレット変換を行い、ウェーブレット領域データを得るステップを備える。本ステップは、いわゆる事前計算である。
【００１８】
本ステップでは、点広がり関数（ＰＳＦ）ｕｚｊ（ｘｈ、ｙｈ）を物体点（０，０，ｚｊ）から計算し、その後、空間領域のＰＳＦをウェーブレット領域のＰＳＦに変換を行う。この場合においてウェーブレット変換は高速ウェーブレット変換であることが好ましい。
【００１９】
なおここで、ＰＳＦｕｚｊ（ｘｈ、ｙｈ）は、スケーリング係数ｓ（ｌ）ｍ，ｎを使用して以下のように表される。なおｌはウェーブレット分解のレベルを意味する。
【数１】
2019015935.tif 000003
【００２０】
また上記式中、ｍ、ｎは整数であり、φはスケーリング関数をそれぞれ意味する。
【００２１】
なお、スケーリング係数のレベルがゼロである場合、元の信号であると考えることができる。すなわち、下記を満たす。
【数２】
2019015935.tif 000004
【００２２】
また、ウェーブレット変換は、下記式で示されるように、ｕｚｊ（ｘｈ、ｙｈ）を、スケーリング係数ｓ（ｌ）ｍ，ｎとウェーブレット係数ｗ（ｌ）ＬＨ，ｍ、ｎ’、ｗ（ｌ）ＨＬ，ｍ、ｎ’、ｗ（ｌ）ＨＨ，ｍ、ｎ’、に分解する。
【数３】
2019015935.tif 000005
【００２３】
本方法において、ｌは、限定されるわけではないが、２以下であることが経験的に質の高い画像を得ることができることを確認しており、好ましいが、ｍ、ｎの値に応じて異なる。
【００２４】
ここで、ウェーブレット変換後のスケーリング係数及びウェーブレット係数の分布について図１に、図２にＰＳＦを、図３に、このＰＳＦのウェーブレット分解の計算結果の例について示しておく。なお、これらの図では、波長６３３ｎｍにおいて、画素のピッチ１０μｍ、ｚｊ＝１００ｍｍとした。また、元々のスケーリング係数及びウェーブレット係数の数は、ホログラムの画素数と同じであり、Ｎｈ２とした。
【００２５】
ところで、後述の足し合わせを行うステップにおいて、そのままの数を足し合わせるとその計算量は膨大となる。そこで、本方法は、（Ｓ２）ウェーブレット領域データのうち一部を選択するステップを備える。具体的には、上記（Ｓ１）において行ったウェーブレット変換において得られるデータのうち、ウェーブレット係数の大きなものだけを選択する。なお本ステップも事前計算の一環として行われる。
【００２６】
本ステップにおいて、一部を選択するステップは、計算量を削減できるものであれば適宜採用可能であるが、小さなウェーブレット係数のものを削除する（ウェーブレット係数の大きなものを残す）ことにより行う。なおこの係数の大きさは振幅により評価・選択することが好ましい。
【００２７】
また、本ステップにおいては、計算量を削減するために下記Ｎｒ係数を定義し、一部を選択する。なお下記式中、Ｗｊ２はＰＳＦｕｚｊの半径であり、ｒは選択率を意味し、振幅の大きなものからどの程度選択するのかを意味する。
【００２８】
この効果については後述の記載から明らかであるが、例えば選択する数をｒ＝１％とする場合、計算量は１／１００となる。
【数４】
2019015935.tif 000006
【００２９】
また、選択された係数は、例えば下記で定義されるベクトルｖｚにより記録しておくことが好ましい。
【数５】
2019015935.tif 000007
【００３０】
なお、上記式中、αｚ，ｋは下記で表現される。
【数６】
2019015935.tif 000008
【００３１】
また、本方法は、上記の通り（Ｓ３）上記（Ｓ１）及び（Ｓ２）を複数回繰り返し、それぞれにおいて、選択された一部のウェーブレット領域データの足し合わせを行うステップを備える。言い換えると、上記（Ｓ１）、（Ｓ２）を複数回繰り返すことで一部のウェーブレット領域データを複数求め、これらを足し合わせることで、足し合わせ後のウェーブレット領域データを得る。
【００３２】
本ステップにおいて、具体的に（Ｓ１）及び（Ｓ２）は、それぞれ物体点の高さ（ｚ）を異ならせて平面的なＰＳＦを複数求め、そのそれぞれに対してウェーブレット変換を行い、更に、それぞれに一部のウェーブレット領域データを選択し、そのそれぞれを足し合わせることにより行う。
【００３３】
ここで、足し合わせは、限定されるわけではないが、下記式に基づき行うことが好ましい。
【数７】
2019015935.tif 000009
【００３４】
ここでδ（ｍ，ｎ）はディラックδ関数である。また、上記式において、位置（ｍ，ｎ）における上記係数ｃｚ，ｋは、重ね合わせにおいて下記の式に従い、水平方向及び垂直方向にシフトし、また、ウェーブレット空間ψ（ｍ，ｎ）に加えられていく。なおｘ１、ｙ１はそれぞれ当初（シフト以前に）計算した物体点の座標位置を意味する。
【数８】
2019015935.tif 000010
【００３５】
なお、上記足し合わせの式において、足し合わせは、同一の領域におけるｃｚｊ，ｋを意味する。例えば、上記図１において、点ａのｃｚｊ，ｋが点ｂに移動した場合は、他幸せを行う一方、上記点ａがｃ又はｄに移動した場合は、足し合わせを行わない。これは、ａとｂは同じ領域であり、ａとｃ、ａとｄはそれぞれ違う領域だからである。
【００３６】
また、（Ｓ４）選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を備える。
【００３７】
本ステップの「逆ウェーブレット変換」は、上記（Ｓ１）において用いたウェーブレット変換と同様の変換の逆変換を用いることができる。
【００３８】
なお図４（ａ）にオリジナルのＰＳＦを、（ｂ）に、本方法をｒ＝２０％の条件で用いて求めたＰＳＦを、（ｃ）に、本方法をｒ＝１０％の条件で用いて求めたＰＳＦを、（ｄ）に、本方法をｒ＝５％の条件で用いて求めたＰＳＦを、それぞれ示す。
【００３９】
また、図４（ｅ）は、二つの物体点から求めた回折パターンを、（ｆ）に、（ｅ）に対し本方法をｒ＝２０％の条件で用いて求めた回折パターンを、（ｇ）に、（ｅ）に対し本方法をｒ＝１０％の条件で用いて求めた回折パターンを、（ｈ）に、（ｅ）に対し本方法をｒ＝５％の条件で用いて求めた回折パターンを、それぞれ示す。
【００４０】
また図５（ａ）は、上記図４（ａ）及び（ｄ）の強度プロファイルを、図５（ｂ）は、上記図４（ｅ）及び（ｈ）の強度プロファイルを、それぞれ示しておく。
【００４１】
以上の各ステップを実行することによって、本方法では、例えばｒ＝１％とした場合、計算量を、通常の高速フーリエ変換を行った場合に比べ、更に１／１００程度にすることが可能となる。すなわち、画質を維持したまま、光波の計算の高速化が可能となる。つまり、より計算時間を短縮することのできるホログラムデータ作成方法及び、ホログラムデータ作成プログラムを提供することができる。
【実施例】
【００４２】
ここで、実際に上記方法について、具体的な例を用いて本発明の効果について確認を行った。具体的には、従来手法（Ｓ．Ｃ．Ｋｉｍら，“Ｅｆｆｅｃｔｉｖｅ  ｇｅｎｅｒａｔｉｏｎ  ｏｆ  ｄｉｇｉｔａｌ  ｈｏｌｏｇｒａｍｓ  Ｏｏｆ  ｔｈｒｅｅ−ｄｉｍｅｎｓｉｏｎａｌ  ｏｂｊｅｃｔｓ  ｕｓｉｎｇ  ａ  ｎｏｖｅｌ  ｕｓｉｎｇ  ａ  ｌｏｏｋ−ｕｐ  ｔａｂｌｅ”，Ａｐｐｌ．Ｏｐｔ．４７，Ｄ５５−Ｄ６２，２００８）と、上記方法（ｒ＝５％、ｒ＝１％）との差異について確認を行った。
【００４３】
本実施例では、３種類の画像（物体点１万点、１０万点、１００万点）を用い、しかもそれぞれｚ＝０．５ｍ、０．２ｍ、０．１ｍ、０．０５ｍと異ならせ、そのそれぞれに対し上記方法による効果を確認した。この結果のうち、ｚ＝０．５ｍ、０．０５ｍの結果を図６乃至８に示す。また、これらの計算の結果について下記表１乃至３に示しておく。なお、この場合においてホログラムのサイズは２０４８×２０４８画素とした。
【表１】
2019015935.tif 000011
【表２】
2019015935.tif 000012
【表３】
2019015935.tif 000013
【００４４】
この結果によると、本手法ではいずれもＰＳＮＲが概ね３０ｄＢであり、画像処理において問題ないといわれている範囲に収まっていることが確認できた。また、計算時間においても、三次元物体間の距離や形状などによって異なるが、従来技術に比べ大幅に削減できていることを確認した。
【産業上の利用可能性】
【００４５】
本発明は、及びホログラムデータ作成プログラム及びホログラムデータ作成方法として産業上の利用可能性がある。
【図１】
2019015935.tif 000014
【図２】
2019015935.tif 000015
【図３】
2019015935.tif 000016
【図４】
2019015935.tif 000017
【図５】
2019015935.tif 000018
【図６】
2019015935.tif 000019
【図７】
2019015935.tif 000020
【図８】
2019015935.tif 000021