6951738 ホログラムデータ作成プログラム及びホログラムデータ作成方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6951738

(24)【登録日】2021年9月29日

(45)【発行日】2021年10月20日

(54)【発明の名称】ホログラムデータ作成プログラム及びホログラムデータ作成方法

(51)【国際特許分類】

   G03H 1/08 20060101AFI20211011BHJP

【ＦＩ】

   G03H1/08

【請求項の数】4

【全頁数】10

(21)【出願番号】特願2017-135103(P2017-135103)

(22)【出願日】2017年7月10日

(65)【公開番号】特開2019-15935(P2019-15935A)

(43)【公開日】2019年1月31日

【審査請求日】2020年7月2日

【新規性喪失の例外の表示】特許法第３０条第２項適用 平成２９年１月９日に発行されたＯｐｔｉｃｓ Ｅｘｐｒｅｓｓ，Ｖｏｌｕｍｅ ２５，Ｉｓｓｕｅ １に発表

【新規性喪失の例外の表示】特許法第３０条第２項適用 平成２９年２月２７日の日刊工業新聞に発表

【新規性喪失の例外の表示】特許法第３０条第２項適用 平成２９年２月２０日に発行されたＯｐｔｉｃｓ Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ，Ｖｏｌｕｍｅ ３９３（２０１７）に発表

【新規性喪失の例外の表示】特許法第３０条第２項適用 平成２９年５月２９日韓国チェジュハンラ大学において開催されたＯＳＡ Ｄｉｇｉｔａｌ Ｈｏｌｏｇｒａｐｈｙ ＆ ３−Ｄ Ｉｍａｇｉｎｇで発表

【新規性喪失の例外の表示】特許法第３０条第２項適用 平成２９年７月６日千葉大学において開催された３次元画像コンファレンス ２０１７で発表

(73)【特許権者】

【識別番号】304021831

【氏名又は名称】国立大学法人千葉大学

(74)【代理人】

【識別番号】100098729

【弁理士】

【氏名又は名称】重信 和男

(74)【代理人】

【識別番号】100204467

【弁理士】

【氏名又は名称】石川 好文

(74)【代理人】

【識別番号】100148161

【弁理士】

【氏名又は名称】秋庭 英樹

(74)【代理人】

【識別番号】100195833

【弁理士】

【氏名又は名称】林 道広

(72)【発明者】

【氏名】下馬場 朋禄

(72)【発明者】

【氏名】伊藤 智義

【審査官】 中村 説志

(56)【参考文献】

【文献】 特表２００２−５３２７７１（ＪＰ，Ａ）

【文献】 米国特許出願公開第２０１６／０３２７９０５（ＵＳ，Ａ１）

【文献】 米国特許第０６０７０１３３（ＵＳ，Ａ）

【文献】 特開平１０−０５５４３３（ＪＰ，Ａ）

【文献】 特表２００６−５０８７７４（ＪＰ，Ａ）

【文献】 特表２０１４−５１２５１４（ＪＰ，Ａ）

【文献】 特開２０１４−２０３４１２（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０３Ｈ １／０８

(57)【特許請求の範囲】

【請求項1】

コンピュータに、
空間領域の点広がり関数（ＰＳＦ）をウェーブレット空間領域のＰＳＦに変換するウェーブレット変換を行い、ウェーブレット領域データを得るステップ、
前記ウェーブレット領域データのうちウェーブレット係数の大きなものだけを選択するステップ、
選択された前記ウェーブレット領域データのうち同一の領域に移動したものだけを足し合わせるステップ、
足し合わせた後の前記ウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を実行させるためのホログラムデータ作成プログラム。

【請求項2】

物体点の高さを異ならせて求めたＰＳＦに対して前記ウェーブレット領域データを得るステップ、及び、前記ウェーブレット領域データのうちウェーブレット係数の大きなものだけを選択するステップを複数繰り返し行う請求項１に記載のホログラムデータ作成プログラム。

【請求項3】

空間領域の点広がり関数（ＰＳＦ）をウェーブレット空間領域のＰＳＦに変換するウェーブレット変換を行い、ウェーブレット領域データを得るステップ、
前記ウェーブレット領域データのうちウェーブレット係数の大きなものだけを選択するステップ、
選択された前記ウェーブレット領域データのうち同一の領域に移動したものだけを足し合わせるステップ、
足し合わせた後の前記ウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を備えるホログラムデータ作成方法。

【請求項4】

物体点の高さを異ならせて求めたＰＳＦに対して前記ウェーブレット領域データを得るステップ、及び、前記ウェーブレット領域データのうちウェーブレット係数の大きなものだけを選択するステップを複数繰り返し行う請求項３に記載のホログラムデータ作成方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、ホログラムデータ作成プログラム及びホログラムデータ作成方法に関する。

【背景技術】

【0002】

ホログラフィ技術によって作製されるホログラムは三次元ディスプレイ、特殊な光学素子、大容量メモリ、セキュリティといった様々な分野で利用される極めて有用なものである。

【0003】

一方で、ホログラムを作製するホログラフィ技術では、三次元物体や画像からの光波を計算するが、この光波の計算にはゾンマーフェルト回折積分、フレネル回折積分等が用いられる。

【発明の概要】

【発明が解決しようとする課題】

【0004】

ところで、上記積分計算では、高速フーリエ変換を用いる方法が考えられる。

【0005】

例えば、伝搬元、伝搬先ともにＮ×Ｎ画素の面間の回折計算を積分で計算する場合を考えると、その計算量はＮ^４に比例して増加する。例えばＮ＝２０４８の場合、一般的なコンピュータを用いた場合でも１時間程度必要となる。これに対し、高速フーリエ変換を用いた場合、その計算量はＮ^２ｌｏｇＮまで減らすことができる。すなわち、上記Ｎ＝２０４８の場合、１秒程度にまで短くすることができる。

【0006】

しかしながら、高速フーリエ変換による光波の計算は平面間の光波伝搬しか計算できず、平面間に限定される結果、三次元的な光波の分布を求めるには高速フーリエ変換を複数回使用する必要があり、上記の時間でも十分短いとは言えない。

【0007】

そこで、本発明は、上記課題に鑑み、より計算時間を短縮することのできるホログラムデータ作成方法及び、ホログラムデータ作成プログラムを提供することを目的とする。

【課題を解決するための手段】

【0008】

上記課題を解決する本発明の一観点に係るホログラムデータ作成プログラムは、コンピュータに、ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、ウェーブレット領域データのうち一部を選択するステップ、選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を実行させるためのものである。

【0009】

また、本発明の他の一観点に係るホログラムデータ作成方法は、ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、ウェーブレット領域データのうち一部を選択するステップ、選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を備えるものである。

【発明の効果】

【0010】

以上、本発明によって、より計算時間を短縮することのできるホログラムデータ作成方法及び、ホログラムデータ作成プログラムを提供することができる。

【図面の簡単な説明】

【0011】

【図1】ウェーブレット係数の分布の例を示す図である。

【図2】ＰＳＦの例を示す図である。

【図3】ＰＳＦのウェーブレット分解の計算結果の例を示す図である。

【図4】ＰＳＦの例を示す図である。

【図5】図４の例における強度プロファイルを示す図である。

【図6】実施例にて作成したホログラムの例を示す図である。

【図7】実施例にて作成したホログラムの例を示す図である。

【図8】実施例にて作成したホログラムの例を示す図である。

【発明を実施するための最良の形態】

【0012】

以下、本発明の実施形態について図面を用いて詳細に説明する。ただし、本発明は多くの異なる形態による実施が可能であり、以下に示す実施形態、実施例における具体的な例示にのみ限定されるわけではない。

【0013】

本実施形態に係るホログラムデータ作成プログラム（以下「本プログラム」という。）は、
コンピュータに、（Ｓ１）ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、（Ｓ２）ウェーブレット領域データのうち一部を選択するステップ、（Ｓ３）上記（Ｓ１）及び（Ｓ２）を複数回繰り返し、それぞれにおいて、選択された一部のウェーブレット領域データの足し合わせを行うステップ、（Ｓ４）選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を実行させるためのものである。

【0014】

限定されるわけではないが、本プログラムは、情報処理装置、いわゆるコンピュータのハードディスクやＣＤ−ＲＯＭ等の記録媒体に記録され、実行されることにより実現できる。より具体的に説明すると、中央演算装置（ＣＰＵ）、ハードディスクやフラッシュメモリ等の不揮発性の記憶装置、ＲＡＭ等の揮発性の記憶装置、これらを電気的に接続するバスライン、キーボードやマウス等の入力装置、液晶ディスプレイ等の表示装置等を備えたコンピュータを用い、上記不揮発性の記憶装置に本プログラムを格納し、揮発性の記憶装置に読み込み実行させることで所望の処理を実行することができる。

【0015】

すなわち、上記の記載から明らかなとおり、上記プログラムを実行することで、（Ｓ１）ウェーブレット変換を行い、ウェーブレット領域データを得るステップ、（Ｓ２）ウェーブレット領域データのうち一部を選択するステップ、（Ｓ３）上記（Ｓ１）及び（Ｓ２）を複数回繰り返し、それぞれにおいて、選択された一部のウェーブレット領域データの足し合わせを行うステップ、（Ｓ４）選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を実行するホログラムデータ作成方法（以下「本方法」という。）を提供することができる。以下、本方法の詳細について説明する。

【0016】

光波の計算は、原理的には点光源から発する光波の足し合わせと考えることができ、本方法はこの光波の足し合わせをウェーブレット変換領域において行うことで、光波計算の計算量を大幅に削減することができる。

【0017】

まず、本方法は、（Ｓ１）ウェーブレット変換を行い、ウェーブレット領域データを得るステップを備える。本ステップは、いわゆる事前計算である。

【0018】

本ステップでは、点広がり関数（ＰＳＦ）ｕ_ｚｊ（ｘ_ｈ、ｙ_ｈ）を物体点（０，０，ｚ_ｊ）から計算し、その後、空間領域のＰＳＦをウェーブレット領域のＰＳＦに変換を行う。この場合においてウェーブレット変換は高速ウェーブレット変換であることが好ましい。

【0019】

なおここで、ＰＳＦｕ_ｚｊ（ｘ_ｈ、ｙ_ｈ）は、スケーリング係数ｓ^（ｌ）_ｍ，ｎを使用して以下のように表される。なおｌはウェーブレット分解のレベルを意味する。

【数1】

【0020】

また上記式中、ｍ、ｎは整数であり、φはスケーリング関数をそれぞれ意味する。

【0021】

なお、スケーリング係数のレベルがゼロである場合、元の信号であると考えることができる。すなわち、下記を満たす。

【数2】

【0022】

また、ウェーブレット変換は、下記式で示されるように、ｕ_ｚｊ（ｘ_ｈ、ｙ_ｈ）を、スケーリング係数ｓ^（ｌ）_ｍ，ｎとウェーブレット係数ｗ^（ｌ）_{ＬＨ，ｍ、ｎ’}、ｗ^（ｌ）_{ＨＬ，ｍ、ｎ’}、ｗ^（ｌ）_{ＨＨ，ｍ、ｎ’}、に分解する。

【数3】

【0023】

本方法において、ｌは、限定されるわけではないが、２以下であることが経験的に質の高い画像を得ることができることを確認しており、好ましいが、ｍ、ｎの値に応じて異なる。

【0024】

ここで、ウェーブレット変換後のスケーリング係数及びウェーブレット係数の分布について図１に、図２にＰＳＦを、図３に、このＰＳＦのウェーブレット分解の計算結果の例について示しておく。なお、これらの図では、波長６３３ｎｍにおいて、画素のピッチ１０μｍ、ｚ_ｊ＝１００ｍｍとした。また、元々のスケーリング係数及びウェーブレット係数の数は、ホログラムの画素数と同じであり、Ｎ_ｈ^２とした。

【0025】

ところで、後述の足し合わせを行うステップにおいて、そのままの数を足し合わせるとその計算量は膨大となる。そこで、本方法は、（Ｓ２）ウェーブレット領域データのうち一部を選択するステップを備える。具体的には、上記（Ｓ１）において行ったウェーブレット変換において得られるデータのうち、ウェーブレット係数の大きなものだけを選択する。なお本ステップも事前計算の一環として行われる。

【0026】

本ステップにおいて、一部を選択するステップは、計算量を削減できるものであれば適宜採用可能であるが、小さなウェーブレット係数のものを削除する（ウェーブレット係数の大きなものを残す）ことにより行う。なおこの係数の大きさは振幅により評価・選択することが好ましい。

【0027】

また、本ステップにおいては、計算量を削減するために下記Ｎ_ｒ係数を定義し、一部を選択する。なお下記式中、Ｗ_ｊ^２はＰＳＦｕ_ｚｊの半径であり、ｒは選択率を意味し、振幅の大きなものからどの程度選択するのかを意味する。

【0028】

この効果については後述の記載から明らかであるが、例えば選択する数をｒ＝１％とする場合、計算量は１／１００となる。

【数4】

【0029】

また、選択された係数は、例えば下記で定義されるベクトルｖ_ｚにより記録しておくことが好ましい。

【数5】

【0030】

なお、上記式中、α_ｚ，ｋは下記で表現される。

【数6】

【0031】

また、本方法は、上記の通り（Ｓ３）上記（Ｓ１）及び（Ｓ２）を複数回繰り返し、それぞれにおいて、選択された一部のウェーブレット領域データの足し合わせを行うステップを備える。言い換えると、上記（Ｓ１）、（Ｓ２）を複数回繰り返すことで一部のウェーブレット領域データを複数求め、これらを足し合わせることで、足し合わせ後のウェーブレット領域データを得る。

【0032】

本ステップにおいて、具体的に（Ｓ１）及び（Ｓ２）は、それぞれ物体点の高さ（ｚ）を異ならせて平面的なＰＳＦを複数求め、そのそれぞれに対してウェーブレット変換を行い、更に、それぞれに一部のウェーブレット領域データを選択し、そのそれぞれを足し合わせることにより行う。

【0033】

ここで、足し合わせは、限定されるわけではないが、下記式に基づき行うことが好ましい。

【数7】

【0034】

ここでδ（ｍ，ｎ）はディラックδ関数である。また、上記式において、位置（ｍ，ｎ）における上記係数ｃ_ｚ，ｋは、重ね合わせにおいて下記の式に従い、水平方向及び垂直方向にシフトし、また、ウェーブレット空間ψ（ｍ，ｎ）に加えられていく。なおｘ_１、ｙ_１はそれぞれ当初（シフト以前に）計算した物体点の座標位置を意味する。

【数8】

【0035】

なお、上記足し合わせの式において、足し合わせは、同一の領域におけるｃ_ｚｊ，ｋを意味する。例えば、上記図１において、点ａのｃ_ｚｊ，ｋが点ｂに移動した場合は、足し合わせを行う一方、上記点ａがｃ又はｄに移動した場合は、足し合わせを行わない。これは、ａとｂは同じ領域であり、ａとｃ、ａとｄはそれぞれ違う領域だからである。

【0036】

また、（Ｓ４）選択された一部のウェーブレット領域データに対して逆ウェーブレット変換を行うステップ、を備える。

【0037】

本ステップの「逆ウェーブレット変換」は、上記（Ｓ１）において用いたウェーブレット変換と同様の変換の逆変換を用いることができる。

【0038】

なお図４（ａ）にオリジナルのＰＳＦを、（ｂ）に、本方法をｒ＝２０％の条件で用いて求めたＰＳＦを、（ｃ）に、本方法をｒ＝１０％の条件で用いて求めたＰＳＦを、（ｄ）に、本方法をｒ＝５％の条件で用いて求めたＰＳＦを、それぞれ示す。

【0039】

また、図４（ｅ）は、二つの物体点から求めた回折パターンを、（ｆ）に、（ｅ）に対し本方法をｒ＝２０％の条件で用いて求めた回折パターンを、（ｇ）に、（ｅ）に対し本方法をｒ＝１０％の条件で用いて求めた回折パターンを、（ｈ）に、（ｅ）に対し本方法をｒ＝５％の条件で用いて求めた回折パターンを、それぞれ示す。

【0040】

また図５（ａ）は、上記図４（ａ）及び（ｄ）の強度プロファイルを、図５（ｂ）は、上記図４（ｅ）及び（ｈ）の強度プロファイルを、それぞれ示しておく。

【0041】

以上の各ステップを実行することによって、本方法では、例えばｒ＝１％とした場合、計算量を、通常の高速フーリエ変換を行った場合に比べ、更に１／１００程度にすることが可能となる。すなわち、画質を維持したまま、光波の計算の高速化が可能となる。つまり、より計算時間を短縮することのできるホログラムデータ作成方法及び、ホログラムデータ作成プログラムを提供することができる。

【実施例】

【0042】

ここで、実際に上記方法について、具体的な例を用いて本発明の効果について確認を行った。具体的には、従来手法（Ｓ．Ｃ．Ｋｉｍら，“Ｅｆｆｅｃｔｉｖｅ ｇｅｎｅｒａｔｉｏｎ ｏｆ ｄｉｇｉｔａｌ ｈｏｌｏｇｒａｍｓ Ｏｏｆ ｔｈｒｅｅ−ｄｉｍｅｎｓｉｏｎａｌ ｏｂｊｅｃｔｓ ｕｓｉｎｇ ａ ｎｏｖｅｌ ｕｓｉｎｇ ａ ｌｏｏｋ−ｕｐ ｔａｂｌｅ”，Ａｐｐｌ．Ｏｐｔ．４７，Ｄ５５−Ｄ６２，２００８）と、上記方法（ｒ＝５％、ｒ＝１％）との差異について確認を行った。

【0043】

本実施例では、３種類の画像（物体点１万点、１０万点、１００万点）を用い、しかもそれぞれｚ＝０．５ｍ、０．２ｍ、０．１ｍ、０．０５ｍと異ならせ、そのそれぞれに対し上記方法による効果を確認した。この結果のうち、ｚ＝０．５ｍ、０．０５ｍの結果を図６乃至８に示す。また、これらの計算の結果について下記表１乃至３に示しておく。なお、この場合においてホログラムのサイズは２０４８×２０４８画素とした。

【表1】

【表2】

【表3】

【0044】

この結果によると、本手法ではいずれもＰＳＮＲが概ね３０ｄＢであり、画像処理において問題ないといわれている範囲に収まっていることが確認できた。また、計算時間においても、三次元物体間の距離や形状などによって異なるが、従来技術に比べ大幅に削減できていることを確認した。

【産業上の利用可能性】

【0045】

本発明は、及びホログラムデータ作成プログラム及びホログラムデータ作成方法として産業上の利用可能性がある。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】