ニュースリリース
令和 3 年 11 ⽉ 25 ⽇
国⽴⼤学法⼈千葉⼤学
国⽴⼤学法⼈京都⼯芸繊維⼤学

光伝播のスローモーション動画記録技術の 55 万倍⾼速シミュレーションに成功〜散乱媒質深部における光の振る舞いの可視化への第⼀歩〜
千葉⼤学⼤学院⼯学研究院 ⾓江崇助教、下⾺場朋禄教授、伊藤智義教授、京都⼯芸繊維⼤学電気電⼦⼯学系 粟辻安浩教授、同⼤学院⼯芸科学研究科博⼠後期課程 井上智好⽒らの研究グループは、光の伝播をスローモーション動画で記録できる技術を、従来と⽐較して 55 万倍⾼速にコンピューターシミュレーション可能な計算アルゴリズムの開発に成功しました。
この成果により、当該技術において⼤きな課題となっていたスローモーション動画に含まれる時空間歪み注1）を除去し、光の振る舞いを忠実に再現可能にする⾒通しが⽴ちました。光の伝播をスローモーション動画で記録できる技術は、⽣体細胞などの散乱注 2)媒質内を伝播する光の振る舞いの可視化を可能にするため、本成果は、散乱により乱された光の情報を復元し、媒質の内部状態や媒質と光との相互作⽤を明らかにすることを⽬指す「散乱透視学」の実現へ向けた第⼀歩であると期待できます。
この研究成果は、2021 年 11 ⽉ 23 ⽇に、⽶国光学会（Optica）が刊⾏する科学論⽂雑誌の「Journal ofthe Optical Society of America A」に掲載されました。

図 1 計算時間の計測結果。いずれも 10 回計測した際の平均値を表す。グラフの横軸は光の情報を画像として表現した場合の⼀辺あたりの画素数、縦軸は計算時間。緑のひし形の点が従来の計算アルゴリズム、⻘の丸の点が開発した計算アルゴリズムによる計算時間の測定結果を表す。⾚の⽮印上に書かれている数値は、従来の計算アルゴリズムに対して開発した計算アルゴリズムが計算時間を何倍⾼速化できているかを表す。

■研究の背景超短パルスレーザー注 3)を利⽤する最先端技術のさらなる発展のためには、レーザーから発せられた光が伝播する様⼦を精緻に観察、評価できる技術の確⽴が求められています。⼀⽅、アインシュタインの相対性理論
では光はこの世の中で最も速く、真空中を毎秒 30 万 km の速さで伝播します。光の⾼速性は世界最⾼速級の⾼速度カメラの性能をはるかに凌駕しており、光の振る舞いを捉えることは極めて困難です。
研究グループは、3 次元画像技術であるホログラフィー注 4)と超短パルスレーザーを組み合わせた、光の伝播をスローモーション動画で記録できる技術に関する研究を⾏ってきました。どのような時空間歪みが⽣じているかについては、従来の計算アルゴリズムによるコンピューターシミュレーションですでに明らかにされています。
-1-しかし、光の伝播を精緻に表現するには、光の時間的な振る舞いと空間的な振る舞いの双⽅を逐次的かつ精密に計算する必要があります。これまでの技術では、スローモーション動画から再現される光の振る舞いに時空間歪みが含まれており、光の伝播の様⼦を忠実に再現できていませんでした。また、⾼性能なコンピューターを⽤いたとしても膨⼤な計算時間を要するため、この技術のコンピューターシミュレーションはほとんどなされておらず、時空間歪みなく光の振る舞いを再現することには⾄っていませんでした。
さらに、従来の計算アルゴリズムでは、光の回折注 5)を数値計算する際に必要な⾜し合わせの計算を直接⾏っており、その計算量は、光の情報を画像として表現する場合、画像を構成する画素の数に⽐例していたため、莫⼤な時間を要していました。通常の光の回折を数値計算する際に必要な⾜し合わせの計算は、⾼速フーリエ変換注 6)を利⽤して⾼速で処理可能であることが知られていましたが、この⼿法は、光の伝播をスローモーシ
ョン動画で記録できる技術における光の回折の数値計算には利⽤できないと考えられていました。
■研究の内容と成果
本研究で開発した計算アルゴリズムに基づいて、光の伝播をスローモーション動画で記録できる技術のコンピューターシミュレーションを⾏いました。計算効率の観点から開発した計算アルゴリズムの有効性を⽰すため、シミュレーションに要する時間を計測し、従来の計算アルゴリズムにおける結果と⽐較しました（図 1）。
その結果、光の情報を 4,096×4,096 画素の画像で表現した場合において、従来の計算アルゴリズムでは約19,000 秒（約 5.3 時間）を要していたのに対し、開発した計算アルゴリズムではわずか 0.034 秒しか必要としませんでした。つまり、開発した計算アルゴリズムにより、従来の約 55 万倍の⾼速化を達成できました。
グラフからは、⼀辺あたりの画素数が増えるほど、⾼速化の割合も⾼くなることが分かりました。従来の計算アルゴリズムでは⼀辺当たりの画素数をさらに多くすると⼀⽇単位で時間を要するため、⼀辺あたり 4,096画素よりも多い条件で⽐較はできていませんが、グラフの傾向からは⾼速化の割合を 55 万倍以上に⾼められることが予想できます。
さらに本研究では、開発した計算アルゴリズムに基づくコン
ピューターシミュレーションで得られたスローモーション動画（図 2-a）を、光学実験により得られたもの（図 2-b）と⽐較しました。拡散板注 7)上を伝播する超短光パルスを記録対象とし、光が伝播する様⼦を分かりやすくするために、拡散板には解像⼒テストチャート（図 2-c）のパターンを貼り付けました。

図 2 拡散板上を伝播する超短光パルスのスローモーション動画。
(a)の⼀部。超短光パルスの形状は直線状であるが、動画においては時空間歪みの影響

光学実験、コンピューターシミュレーションのいずれの結果においても、⼸状の形をした超短光パルスの像が再⽣されていま

コンピューターシミュレーションで⽣成したホログラム画像から再⽣される動画により⼸状に歪んで再⽣される。

(b)

光学実験で得られたホログラム画像から再⽣される動画の⼀部。

(c)

超短光パルスが伝播する様⼦の模式図。超短光パルスは向かって左から右へ伝播する。

す。これは、再⽣される動画に時空間歪みが発⽣している影響によるものです。この時空間歪みも含めて、コンピューターシミュレーションによる結果は光学実験による結果とよく⼀致しており、今回開発した計算アルゴリズムが有効であることが⽰されました。
-2-

本研究では、光が通る経路の⻑さ（光路⻑）に基づいて計算条件を決定することで、光の回折を⾼速フーリエ変換に基づいて数値計算する⼿法を、光の伝播をスローモーション動画で記録できる技術へ適⽤できることを新たに⾒いだしました。光路⻑を考慮することは光の時間的な振る舞いを取り扱うこととも⾔えます。つ
まり、今回⾒いだした計算アルゴリズムは、光の時間的な情報を考慮する点と⾼速な数値計算を可能にする点の両⽅を同時に満たすことを初めて可能にしました。
■今後の展望本研究で提案した計算アルゴリズムにより、光の伝播をスローモーション動画で記録できる技術を⾼速でコンピューターシミュレーションできるようになりました。この成果は、実験で得られたスローモーション動画の解析、評価を従来よりも格段に効率化できるだけでなく、スローモーション動画に⽣じている時空間歪みの影響を除去するアルゴリズムの開発へもつながります。
本成果を発展させ、より複雑な状況における光の伝播のコンピューターシミュレーションを可能にして、将来的には⽣体細胞に代表される散乱媒質内の深部を伝播する光の時空間的な振る舞いを精緻に捉え、媒質の内部状態や媒質と光との相互作⽤を明らかにする未踏技術の確⽴を⽬指します。
■研究プロジェクトについて
本研究は、独⽴⾏政法⼈⽇本学術振興会（JSPS） 学術変⾰領域研究（A）20H05887、特別研究員奨励費20J23542、公益財団法⼈ ⾼柳健次郎財団、公益財団法⼈ 島津科学技術振興財団の⽀援により遂⾏されました。
■⽤語解説
注 1）時空間歪み︓スローモーション動画に含まれる時間的な誤差と、その誤差によって超短光パルスの像に⽣じる空間的な歪みを表します。これらの誤差や歪みにより、スローモーション動画では実際の光の振る舞いとは異なる様⼦が再現されてしまいます。通常の動画撮影に例えるなら、実際に⽬で⾒た様⼦に対して、撮影した動画では対象がブレたり形状が歪んだりして、元の様⼦を正しく把握できなくなってしまうことに相当します。
注 2）散乱︓光がなんらかの物質中を伝播する際、その物質によって光の進⾏⽅向が様々な⽅向へと変化させられる現象です。
注 3) 超短パルスレーザー︓発光時間の極めて短い光を放つレーザーです。⼀般的なレーザーは光を持続的に放ち続けるのに対して、パルスレーザーは細かい時間間隔で途切れ途切れに光を放ちます。パルスレーザーの中でも 10 兆分の１秒程度以下の時間だけ光を放つことができるものは超短パルスレーザーと呼ばれ、フェムト秒パルスレーザーなどがあります。超⾼速現象の誘起や材料の微細加⼯に⽤いられます。フェムトという接頭辞は 1000 兆分の 1 を表します。2018 年にノーベル物理学賞を受賞し、最先端の情報通信、微細加⼯、医療など幅広い分野で利⽤されています。
注 4) ホログラフィー︓光の⼲渉と回折を利⽤して、物体からやってくる光のすべての情報を記録、再⽣できる 3 次元画像技術です。私たちが物体を⾒るときに認識している、物体を透過または物体で反射した光を物体光と呼びます。物体光と基準となる光（参照光）とを⼲渉させ、⼲渉で⽣じる光の明るさ分布（⼲渉縞）を感光材料やイメージセンサーで画像として記録します。⼲渉縞の画像が記録されている媒体がホログラムです。記録に⽤いた参照光と同じ光をホログラムに照射すると、ホログラムに記録されている⼲渉縞により光が回折して物体光が再⽣され、奥⾏きの情報を含む物体の 3 次元像を観察できます。

-3-注 5) 光の回折︓光が伝播する際、光の直進性を考慮すると本来は到達不可能である領域、例えば障害物の背後などへも回り込む現象のことです。
注 6) フーリエ変換︓数学的変換の⼀つであり、データ解析⼿法として⼯学をはじめとした様々な分野において広く⽤いられています。⾼速フーリエ変換は、フーリエ変換をコンピューター上で⾼速に処理するためのアルゴリズムです。
注７）拡散板︓光を様々な⽅向へと散乱させる機能を持つ板状の材料です。拡散板の⾝近な例としてすりガラスがあげられます。
■論⽂情報
・論⽂タイトル︓
FFT-based simulation of hologram-recording process for light-in-flight recording by holography
・著者︓
Takashi Kakue, Tomoyoshi Inoue, Tomoyoshi Shimobaba, Tomoyoshi Ito, and Yasuhiro Awatsuji
・掲載誌︓Journal of the Optical Society of America A
・DOI︓https://doi.org/10.1364/JOSAA.441585-4-