昭和 37 年３月８日 第３種郵便物認可／昭和 37 年２月 23 日国鉄東局特別扱承認雑誌 1243号／ 平成 23 年 5 月 10 日発行 （毎月１回 10 日発行）第 50 巻第 5 号 ISSN 0453-4662

計
測

と
制
御

VOL.50 2011

計測と制御 5
ミニ特集

回折イメージング 〜位相回復の新展開〜●総論
新たなナノスケールイメージング─ 古くて新しい位相回復問題 ─●解説コヒーレントX線イメージング
光 ─ 回折強度からのレンズレスイメージング─
電子回折イメージング─入射ビームの角度拡がり─
位相回復 ─ 計算アルゴリズム─

●リレー解説
脳機能計測と生体信号入出力●ミニ特集

第5回：ダイポールレイヤーを用いた脳機能イメージングとその応用

位相回復による解析イメージング●学界だより
第53回自動制御連合講演会●部門だよりシステム・情報部門学術講演会2010（SSI 2010）

回折イメージング 〜位相回復の新展開〜

計
測
自
動
制
御
学
会

●支部だより
第19回中国支部学術講演会

公益社団法人計測自動制御学会
URL : http://www. sice. or. jp/

総論：ミニ特集
回折イメージング〜位相回復の新展開〜

新たなナノスケールイメージング—古くて新しい位相回復問題—
郷

原 一

寿＊

＊北海道大学大学院 工学研究院 北海道札幌市北区北 13 条西 8 丁目
＊Research Faculty of Engineering, Hokkaido University, Kita 13, Nishi 8,
Kita-ku, Sapporo, Hokkaido, Japan＊E-mail: gohara@eng.hokudai.ac.jp

キーワード：回折イメージング (diffractive imaging)，位相回復 (phaseretrieval)，X 線 (X-ray)，光 (photon)，電子 (electron)，コヒーレンス(coherence)．

見るということは，知覚の中でも格段に大きな割合を占以上を背景に，本ミニ特集号では，回折イメージングがめており，人間の認識において重要な役割を果たしている．

展開されている応用領域として，X 線，光，電子の 3 つの百聞は一見に如かず，見ることは理解することの第一歩で分野を取り上げ，さらにこの手法の基礎となるフーリエ反ある．しかし，人間が目で直接見ることができるものは，可復位相回復法のアルゴリズムについて，情報理論の観点か

視光の波長帯に限られている．見ることの魅力に憑りつから解説を加えた．全体は，それぞれに対応した，以下の 4れた人間は，知恵を働かせて，可視光よりもより長く，さ編の解説から構成されている．

らに短い領域の波長でも見ることを追い駆け続け，多くのイメージングに関する手法を考え出した．中でも，原子のC 2011 SICE
JL 0005/11/5005–0313 • X 線はこの手法が最も活発に展開されている分野である．X 線での初期の検証実験から深く関わり，積極スケールでイメージングすることは，分野を超えて大きな的に研究を進めている西野先生に，この手法の基礎，影響を与える．たとえば，DNA やカーボンナノチューブの原子スケール構造解明とその後の展開は，原子が識別でSPring-8 で展開されている X 線回折顕微鏡による研究，X 線自由電子レーザーによるフェムト秒コヒーレ

きる空間分解能で物質の構造を知ることが，いかに大きなントイメージングの現状と最新の研究などについて解

影響を社会に与えるかを示す典型的な例である．これらの構造解明には，X 線または電子による回折パターンが重要な役割を果たしている．

説していただいた．

• 光の分野は位相回復問題を基礎と応用に渡って最も広く研究している領域である．この分野で新たな手法を

原子スケールの構造を前提にする分野では，回折という現活発に提案し続けられている中島先生に，光の分野か

象が，周期性をもつ結晶に対してのみ問題になるという一般

ら見た歴史，X 線回折イメージングの各種方法の分類的な常識がある．しかし，“回折イメージング (Diﬀractiveと現状，そしてご自身の最近の研究である開口アレイ

Imaging)” または “回折顕微法（Diﬀraction Microscopy)”

フィルターによる位相回復などについて解説していたと呼ばれる新たなイメージング手法の出現によりこの常識

だいた．

よって計測された回折パターンを元に，計算機による数値• 電子もド・ブロイ波長で特徴付けられる波の性質をもつことから，X 線と並行して活発に研究が展開されて計算によって実像を得るイメージング手法であり，物理的いる．原子分解能でイメージングできる電子線領域でなレンズの機能をデジタル計算によって実現する “デジタのこの手法の可能性を追求している郷原が，入射ビー

ルレンズ” と言える．

ムの角度拡がりの問題を中心に取り上げ，電子顕微鏡

は大きく変わりつつある．回折イメージングとは，実験にこの方法が X 線領域で可能であることが，実験的に示さでの具体的な実験例について解説した．

れたのは 1999 年のことである．波長オーダーで十分な性• 最後に，この手法の基礎となる位相回復問題の数学的能をもつレンズがない X 線の分野では，結晶であることを定式化，位相回復の基本アルゴリズム，情報量による必要とせず，レンズも要らないこの手法は，大きな驚きと考察について塩谷・郷原が解説した．

期待で迎えられた．そして，これまでの約 10 年で飛躍的な回折イメージングは，物質科学と情報科学などの異なる発展を遂げ，X 線の分野では現在もその勢いをさらに増し学問の融合領域である．主要なターゲットは，ライフイノているように見える．

ベーション・グリーンイノベーションを支える重要物質の一見，にわかに登場したように見えるこの手法は，位相ナノスケール構造解明である．古くて新しい位相回復問題

回復問題と呼ばれる基本的な問題を含んだパラダイムに立を基礎にした回折イメージングに関する研究の一端をお伝脚しており，寄って立つところの土台に対して，新たな問えすることで，本特集号が，多くの分野で新たな挑戦をす

題提起をしている．

る方々のヒントになることを期待したい．
（2011 年 2 月 14 日受付）

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

313

解説：ミニ特集
回折イメージング〜位相回復の新展開〜

コヒーレント X 線イメージング
西

野 吉

則＊

＊北海道大学 電子科学研究所 北海道札幌市北区北 21 条西 10 丁目＊Research Institute for Electronic Science, Hokkaido University, Kita 21
Nishi 10, Kita-ku, Sapporo, Hokkaido, Japan＊E-mail: yoshinori.nishino@es.hokudai.ac.jp

1.

はじめに

キーワード：コヒーレント X 線 (coherent X-rays)，コヒーレントイメー
ジング (coherent imaging)，回折顕微法 (diffraction microscopy)，ホロ
グラフィー (holography)，自由電子レーザー (free-electron laser)．
C 2011 SICE
JL 0005/11/5005–0314 

メージングの概念図を示す．光を完全に透過する透明な試料 (位相物体) に対しても，試料内部に非一様な屈折率の分放射光や自由電子レーザーなどの短波長コヒーレント光布があれば，試料背後で光の波面がゆがみ位相シフトが生の光源技術の近年の目覚ましい発展に伴い，コヒーレント

じる．この位相シフトを各種の方法で検出することにより，光を利用したイメージング研究が世界的に活発に行われて位相コントラストイメージングが実現する．コヒーレント

いる．コヒーレント光とは，図 1 に示す入射波のように，イメージングは，位相板などの光学素子を必要としない理

波面の揃った光を意味する．

想的な位相コントラストイメージング法である．

回折顕微法1)〜5) やホログラフィーなどのコヒーレントイ

本稿では，コヒーレント X 線が可能にする先端的イメー

メージングでは，通常の結像型や走査型の顕微鏡で必要とジングについて解説する．

される像形成のためのレンズを用いず，試料からのコヒー2.

レント散乱パターンを計測し，計算機で試料像を再構成す
る．レンズを必要としないため，開口数の大きなレンズのX 線回折顕微法

構造解析法を，結晶以外の試料にも適用できるように拡張

2.1 スペックルとオーバーサンプリング条件
X 線の波が試料を通過する際に起こる波面のゆがみは，
遠方界では Fraunhofer 回折パターンとして観測される．
Fraunhofer 回折パターンから試料構造を再構成する手法がX 線回折顕微法である．図 2 に X 線回折顕微法の概念図をした手法と言える4) ．

示す．

製作が困難な，X 線領域の光に対して特に有効である．オ
ングストロームの波長をもつ X 線を用いると原理的には原子分解能でのイメージングも可能である．従来の X 線結晶

X 線は透過性に優れ，厚みのある試料の 3 次元観察に適している．しかし，従来の X 線顕微鏡の多くは，X 線の吸レーザー光においてはなじみ深い，スペックルと呼ばれる

収コントラストを可視化するレントゲン写真の延長にすぎ

斑点模様が観測される．回折パターン中のスペックル構造なかった．このため，X 線をほとんどすべて透過する細胞は，試料によって散乱された X 線の干渉によって起こるたや細胞小器官などの試料に対しては，像のコントラストがめ，試料構造の僅かな違いがスペックルの顕著な違いを生悪く，内部構造を観察することは困難であった．X 線を用

む．X 線回折顕微法は，回折パターンの試料構造に対する

無秩序な構造をもつ試料からの Fraunhofer 回折では，いこれら生体試料を高いコントラストで観察するには，透明な試料を可視化できる位相コントラストイメージングが求められる1) ．
図 1 に，コヒーレント光を利用した位相コントラストイ

図 1 コヒーレント光を利用した位相コントラストイメー
ジングの概念図．試料を通過した光の波のゆがみ (位相シフト) を可視化する．

314

図 2 X 線回折顕微法の概念図．試料にコヒーレント X 線を照射し遠方界で回折パターンを計測する．レンズの代わりに計算機を用い，回折パターンから試料像を再構成する．
3 次元イメージングでは，試料を回転させ，さまざまな入射角での回折パターンを計測する．

計測と制御

第 50 巻

第5号

2011 年 5 月号

この敏感性を利用している．

まず，逆空間での既知情報として，回折強度の測定デー

Fraunhofer 回折パターン中の 1 つ 1 つのスペックルの大きさは，試料 (散乱領域) の大きさに反比例する．これは，タを用いる．この際，位相は実験的に得られないため，各

小さなスリットを通った波が，スリットの開口幅に逆比例

実空間では，試料が存在するであろうと推定される領域

反復における計算値をそのまま採用する．

した角度広がりで伝播することからも理解できる．量子力(サポート) の外では電子密度がゼロであることや，物理的学的には，不確定性原理として知られている現象である．

に電子密度が非負であることなどを拘束条件として課すこスペックルの大きさは，光の波長を λ，試料の大きさを a，とができる．ここで，試料が存在するであろう領域は，ス

試料と検出器の間の距離を D とすると，およそ λD/a でペックルのサイズが試料サイズに反比例するため，スペッ

与えられる．

クルサイズから推定することができる．具体的な実空間拘実験では，回折パターン中の 1 つ 1 つのスペックルを解束条件にはさまざまな流儀があるが，Fienup により提案さ像できるよう細かい間隔で回折強度をサンプリング (計測)れた HIO (hybrid input-output) アルゴリズムが一般的にする．回折パターンのサンプリング間隔に対する要請は，用いられている8) ．

オーバーサンプリング条件と呼ばれ，検出器のピクセル面

情報量が，未知の試料構造の情報量を上回る条件を与える．

2.3 X 線回折顕微鏡装置
X 線回折顕微法の基礎となるアイデアは 50 年以上も昔に提唱されたが3) ，世界初の実験報告は 1999 年に行われた4) ．
実証までに長い年月を要したのは，コヒーレント X 線の利このため，オーバーサンプリング条件は，回折パターンか用が放射光の発展とともに近年になりやっと可能となった

ら試料構造を再構成する必要条件となる．

ためである．その後，研究は世界的な広がりをみせ，実験

2.2 反復的位相回復法
試料の電子密度分布 (屈折率分布) と Fraunhofer 回折波の複素振幅とは，フーリエ変換の関係にある．X 線回折波の複素振幅を計測することができれば，逆 Fourier 変換に装置やデータ解析法の改良が進められた．特に，波長の短積が，スペックルの面積の半分以下であることが求められ
る6) ．オーバーサンプリング条件は，測定した回折強度のい硬 X 線の領域では，日本の大型放射光施設 SPring-8 において多くの先駆的な研究が行われた．
図 4 にわれわれが開発した X 線回折顕微鏡装置の写真を

よって電子密度分布を再構成できる．しかし，測定可能な示す．X 線回折顕微鏡では通常，マイクロメートルほどの物理量は光子数から得られる複素振幅の絶対値のみで，複サイズの試料の内部構造を観察するが，そのような小さな

素振幅の位相は直接的には測定できない．このため，電子試料からの X 線散乱強度は極めて弱い．このため，空気に密度分布を再構成するためには，なんらかの方法で回折波よるノイズ散乱も無視できないため，測定は通常真空中での位相を回復する必要がある．これがよく知られた X 線回行われる．

折における位相問題である．

回折パターンは試料の下流に設置された 2 次元検出器に

X 線回折顕微法では，X 線回折パターンからの位相回復

よって計測される．われわれの装置では，試料と検出器のに Gerchberg と Saxton に端を発する反復法を用いる7) ．図間の距離 D を連続的に変えられる工夫がされている．こ3 に反復的位相回復法の概念図を示す．反復法では，ランダ

れにより，さまざまなサイズの試料に対しても，固定のピ

ムな電子密度分布 (またはランダムな回折波の位相) から出

クセルサイズをもった 2 次元検出器で，オーバーサンプリ

発し，Fourier 変換と逆 Fourier 変換を繰り返す．そして，ング条件を満たした測定を行うことができる．

各反復において，実空間および逆空間 (Fourier 空間) で既

回折パターンを測定する 2 次元検出器には，スペックルを解像する高い空間分解能 (点拡がり関数) と，シングル

知の情報を補っていく．

フォトンレベルの微弱な回折強度を計測する高効率と低ノ

図 3 反 復 的 位 相 回 復 法 の 模 式 図 ．実 空 間 と 逆 空 間(Fourier 空間) で既知の情報を与えながら，Fourier 変
換と逆 Fourier 変換を繰り返すことで，位相を回復する．

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

図 4 X 線回折顕微鏡装置．SPring-8 BL29XU に設置された様子を示す．

315

図 5 X 線回折顕微法によって観察したヒト染色体の 2 次
元投影像．イメージコントラストを人工的に高める試料処
理なしに，初めて染色体の軸状の構造が観察された．
イズが求められる．これらの条件を満たす 2 次元検出器として，現在多くの実験では，X 線を直接照射するタイプの

X 線 CCD 検出器が現在用いられている．一方，一般的な
X 線実験では通常，X 線を蛍光板に照射し可視光に変換するタイプの CCD 検出器が用いられる．蛍光板を使った場
合，X 線から可視光への変換効率を上げるため蛍光板を厚くすると，空間分解能が悪化するという問題がある．

図 6 X 線回折顕微法によって観察したヒト染色体の 3 次
元像．染色体の表面のみならず，内部構造も観察できる．

2.4 X 線回折顕微法による細胞内構造イメージング
顕微鏡研究においては，生体試料の観察は大きなテーマ

パク質に重金属に結合させ染色したり，蛍光分子で標識すである．細胞小器官の高空間分解能での構造観察では，クることによって，特定のタンパク質のイメージコントラス

ライオ電顕が現在最も有力な手法であるが，厚い細胞小器

トを人工的に高めていた．われわれの研究では，これらイ

官や細胞全体の可視化には，透過能の優れた X 線顕微鏡がメージコントラストを高める人工的な試料処理をせず，染潜在的優位性を有している9) ．X 線に対してこれらマイク

色体の軸状構造を観測した初めての例となる．さらに，わロメートルサイズの生体試料は，吸収によるコントラスト

れわれが得たヒト染色体の 2 次元投影像では，軸状の構造が得られないほぼ理想的な位相物体であるが，X 線回折顕がほぼ規則的に曲がりくねって，波状になっていることが微法は回折現象を利用することによって，試料による X 線観察できた．この構造は，蛍光顕微鏡を用いてこれまで観波の位相シフトを高感度で検知することができる．

察されていた，らせん状の軸構造との類似性を示している．

われわれは，X 線回折顕微法を用いた世界初の生体試料

図 6 に，再構成されたヒト染色体の 3 次元電子密度マッの観察として，2003 年に大腸菌の 2 次元投影像のイメージ

プを示す．3 次元再構成は，試料を回転させさまざまな入ングを報告した10) ．その後，世界のいくつかのグループに射角で測定した回折強度データを 3 次元逆空間 (Fourier 空より X 線回折顕微法による生体試料の観察報告がなされた間) に当てはめ，3 次元 Fourier 変換を使った反復的位相回が，いずれも 2 次元投影像の観察であり，3 次元的に生体復法を適用することによって行った．3 次元再構成像にお試料を観察する研究が待たれていた．

いても，軸付近に電子密度の高い構造が観察された．

近年，われわれは X 線回折顕微法による細胞小器官の 3

われわれの研究により，コヒーレント X 線を使った新し．X 線回折顕微法といい顕微鏡が，イメージコントラストを人工的に高める染色

う特定の手法に限らず，エネルギーの高い X 線を使って細や標識などの試料処理の必要なく，細胞小器官の内部構造

胞の中を 3 次元的に観察したのは世界初である．試料にはを高いイメージコントラストで観察できることが実証され

ヒト染色体が用いられた．

た．これにより，細胞内構造を観察する研究に新たな可能次元観察に世界で初めて成功した

11)

図 5 にコヒーレント X 線回折パターンから再構成された，ヒト染色体の 2 次元投影像を示す．再構成像では，2 つの姉妹染色分体が，明瞭に可視化されている．さらに，特筆すべき構造として，染色体の軸付近に，幅 200 nm 程のイメージ強度の高い領域が観察された．染色体の軸状構造
は，これまで，免疫電子顕微鏡や蛍光顕微鏡を用いて観察されてきたが，これら従来法では，軸付近に存在するタン316

性を開いた．

3.

自由電子レーザーを用いたコヒーレ
ントイメージング

3.1 X 線自由電子レーザー
一昨年 2009 年にアメリカの LCLS (Linac Coherent
Light Source) で，史上初のオングストローム・オーダーの計測と制御

第 50 巻

第5号

2011 年 5 月号

波長をもつ X 線レーザーが実現し，X 線科学の新しい幕がFEL のパルス幅はサブピコ秒ときわめて短く，シングル

開いた12) ．日本の X 線レーザー計画が続いており，兵庫県

ショットデータから試料像が再構成できれば，サブピコ秒のの播磨科学公園都市にある SPring-8 サイトで建設が進め

露光時間をもつコヒーレントイメージングが実現する．こられている．今年 2011 年の夏には X 線レーザーの発振がれにより，超高速で起こる材料科学現象も，ぶれることな

期待される

13)

．これらアメリカ，日本の施設に続き，ヨーロッパで X 線レーザーの建設が進められているほか，世界各地で建設に向けた議論が始まっている．

く撮影できる．
実験では，HERALDO (holography with extended ref-erence by autocorrelation linear diﬀerential operation)

通常のレーザーが束縛電子を用いるのとは異なり，これと呼ばれる近年開発された新しいホログラフィーの手法を

ら X 線レーザーは加速器で加速された自由電子を用いるた

用いた16) ．HERALDO は，古くから知られるフーリエ変め，自由電子レーザー (FEL) と呼ばれる．

換ホログラフィーを改良した手法である．フーリエ変換ホ

X 線 FEL (XFEL) は，高い干渉性，超短パルス性，高ログラフィーでは，試料の隣に参照光源となるピンホール

いピーク輝度をもっており，これら優れた特徴を用いるこを置くが，このピンホールの大きさがフーリエ変換ホログとで，コヒーレントイメージングの飛躍的な発展が期待さラフィーでの再構成像の空間分解能を決める．高い空間分れる．たとえば，コヒーレントイメージングを時間軸方向解能を得るためにピンホール径を小さくすると，参照光のに押し広げた 4 次元イメージングへの新たな可能性が開か

強度が弱くなり，結果としてホログラムデータの統計精度れる．X 線の短波長性と XFEL の超短パルス性を活用するが悪化するという問題があった．HERALDO ではこの問と，究極的には，原子分解能での超高速の構造ダイナミク

題を克服し，広がった参照光源を用いても，測定データのス測定が可能となる．図 7 に，XFEL を用いたポンプ・プ

統計精度を高く保ったまま，高い空間分解能が実現可能でローブ法による構造ダイナミクス測定の概念図を示す．測ある．

定した回折パターンの時間列からは，究極的には原子・分

図 8 に SCSS 試験加速器を用いた HERALDO 測定の様子の動画が得られる．

子を示す．波長 61 nm，パルス幅〜100 fs，パルスエネル

3.2 フェムト秒コヒーレントイメージング
われわれは，XFEL を利用した材料科学現象の超高速

ギー >10 µJ の FEL を試料に照射し，試料の下流 776 mmに置かれた CCD 検出器でホログラムを計測した．CDD 検

コヒーレントイメージングに向けた研究を進めている．

出器は市販のもので，全画素の読み出しに数秒を要する．

SPring-8 のサイト内に XFEL の原理実証機として作られた SCSS 試験加速器からの極紫外線 FEL (EUV-FEL)14)

これは，実験で用いた EUV-FEL の繰り返し周波数であるを用いて基礎実験を行っている．
構造ダイナミクス測定に向けた第一歩として，シングル
ショットでのコヒーレントイメージング実験を行った15) ．

図 7 X 線自由電子レーザーを用いた原子分解能での超高速の構造ダイナミクス測定の概念図．究極的には，原子・分子の動画が得られる．

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

図 8 SCSS 試験加速器からの極紫外自由電子レーザーを用いたフェムト秒シングルショット・ホログラフィー測定の様子 (上) と概念図 (下)．シングルショット計測は，加速器からのタイミング信号と CCD 検出器の露光信号をトリガーとして動作するパルスセレクターを用いて実現した．

317

図 9 試料の走査イオン顕微鏡像．試料となる SPring-8のロゴマークの近傍に，参照光を作り出す縦および横のス
リットが配置されている．

図 11 フェムト秒シングルショットポログラムから再構成した試料像．4 つのコピー像を適切に足し合わせることにより (c) に示す最終的な再構成像が得られる．
データを示す．縦および横に伸びる干渉縞がのったストリー

20Hz に比べ格段に遅い．読み出しの遅い CCD 検出器でシ

クは，それぞれ，横および縦のスリットによる回折の寄与

ングルショット計測を実現するため，われわれはビームラである．

インのパルスセレクターを利用した17) ．パルスセレクター

ホログラフィーでは位相情報が，ホログラムの実験デー

は，FEL の繰り返し周期と同期して，任意数のパルスを切タに含まれているので，位相回復の必要がなく，単純な数り出せる．本計測では，パルスセレクターの切り出しパル

学的操作で試料像を再構成できる．図 11 に，図 10 のホロ

ス数を 1 に設定し，外部トリガー端子に CCD の露光開始

グラムから再構成した試料のフェムト秒スナップショット

信号を入力した．

を示す．図 11 (c) に示す最終的な再構成像では，試料であ図 9 に測定に用いた試料の走査イオン顕微鏡像を示す．

る SPring-8 のロゴマークが明瞭に観察できる．

試料は 100 nm 厚の Si3 N4 に金を 800 nm 蒸着したのち，この実験により，FEL を用いたサブピコ秒露光でのコ

集束イオンビーム (FIB) を用いて貫通パターンを加工した．

ヒーレントイメージングが実現した．シングルショット計測

パターンは，試料となる SPring-8 のロゴマークと，縦およ

したホログラムからの試料像の再構成の成功を受けて，わび横の参照スリットで構成されている．パターンの大きさ

れわれは，フェムト秒レーザーによって引き起こされる構や形状は，シングルショット計測に適したホログラムデー

造ダイナミクスの測定を進めている．

タが得られるよう設計した．
図 10 にシングルショット計測したフェムト秒 HERALDO

4.

おわりに

放射光利用研究は，近年その裾野を大きく広げ，比較的小型の放射光専用施設が世界各地で数を増やしている．これら新しい放射光施設では，コヒーレント X 線の利用研究が大きな特色となっており，X 線回折顕微法をはじめとするコヒーレントイメージングを主目的とするビームラインの運用も始まりつつある．
また，新世代の X 線である X 線自由電子レーザーの利用も始まり，実験結果の報告が徐々に現われている．本稿では，XFEL が切りひらく新たな可能性として，フェムト秒の構造ダイナミクスイメージングを示したが，他にも多くの新しい発展が期待される．たとえば，結晶化が困難または不可能な種類の生体超分子試料をつぎつぎと XFEL に照射することにより，3 次元構造解析を行うというシナリオが議論されている．XFEL の超短パルス性能により，試料構造が破壊される前に X 線回折測定を行う可能性が開け18) ，
生体試料観察における従来の分解能限界を超えた高空間分
図 10 実験で測定したフェムト秒シングルショットポログ
ラム．100 フェムト秒ほどのパルス幅をもつ極紫外自由電子レーザーを 1 パルス照射して計測した．

318

解能が達成できると期待されている．XFEL を利用した測定は，技術的チャレンジを伴うが，これまで見ることのできなかった世界が開ける大きな可能性を秘めている．

計測と制御

第 50 巻

第5号

2011 年 5 月号

最後に，以下の方々をはじめとする共同研究者に感謝する (敬称略)：田中義人，伊藤基巳紀，高橋幸生，三村秀和，
松 山智 至 ，前 島 一 博，志 村ま り，城地 保 昌 ，苙口友隆，
別所義隆，竹内昌治，岡田真，岡谷基弘，野崎公彦，大路祐 介，松下雄多，堤 良輔 ，矢 橋牧 名，永 園充，富 樫格 ，
大橋治彦，松井真二，松原英一郎，山内和人，石川哲也．
本研究は科研費および X 線自由電子レーザー利用推進課
題の助成を受けたものである．
（2011 年 1 月 23 日受付）
参

考

文 献

1） 西野吉則：コヒーレント X 線が明かす細胞の内部世界，パリ
ティ(丸善)，24, 14 (2009)2） 西野吉則，石川哲也：X 線回折顕微法の原理，放射光，19, 3(2006)
3） D. Sayre: Some implications of a theorem due to Shannon,
Acta Crystallogr., 5, 843 (1952)
4） J. Miao, P. Charalambous, J. Kirz and D. Sayre: Extending
the methodology of X-ray crystallography to allow imaging of micrometre-sized non-crystalline specimens, Nature,400, 342 (1999)
5） I.K. Robinson, I.A. Vartanyants, G.J. Williams, M.A.
Pfeifer and J.A. Pitney: Reconstruction of the Shapesof Gold Nanocrystals Using Coherent X-Ray Diﬀraction,
Phys. Rev. Lett., 87, 195505 (2001)
6） J. Miao, D. Sayre and H.N. Chapman: Phase retrieval fromthe magnitude of the Fourier transforms of nonperiodic objects, J. Opt. Soc. Am. A, 15, 1662 (1998)
7） R.W. Gerchberg and W.O. Saxton: Phase determination
for image and diﬀraction plane pictures in the electron microscope, Optik (Stuttgart), 34, 275 (1971)
8） J.R. Fienup: Phase retrieval algorithms: a comparison,
Appl. Opt., 21, 2758 (1982)
9） R. Henderson: The potential and limitations of neutrons,
electrons and X-rays for atomic resolution microscopy of
unstained biological molecules, Q. Rev. Biophys., 28, 171(1995)
10） J. Miao, K.O. Hodgson, T. Ishikawa, C.A. Larabell, M.A.
LeGros and Y. Nishino: Imaging whole Escherichia coli bacteria by using single-particle x-ray diﬀraction, Proc. Natl.
Acad. Sci. USA, 100, 110 (2003)

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

11） Y. Nishino, Y. Takahashi, T. Ishikawa, N. Imamoto and K.
Maeshima: Three-Dimensional Visualization of a Human
Chromosome Using Coherent X-Ray Diﬀraction, Phys. Rev.
Lett., 102, 018101 (2009)12） P. Emma, et al.: First lasing and operation of an angstromwavelength free-electron laser, Nat. Photonics, 4, 641(2010)
13） 北村英男，新竹積，石川哲也：SPring-8 におけるオングス–ト
ローム–FEL–開発，放射光，16, 65 (2003)
14） T. Shintake, et al.: A compact free-electron laser for generating coherent radiation in the extreme ultraviolet region,Nat. Photonics, 2, 555, (2008)
15） Y. Nishino, Y. Tanaka, M. Okada, M. Okaya, Y. Uozaki, K.
Nozaki, M. Yabashi, M. Nagasono, K. Tono, H. Kimura, H.
Ohashi, S. Matsui, T. Ishikawa and E. Matsubara: Femtosecond Snapshot Holography with Extended Reference
Using Extreme Ultraviolet Free-Electron Laser, Appl. Phys.
Express, 3, 102701 (2010)16） M. Guizar-Sicairos and J.R. Fienup: Holography with extended reference by autocorrelation linear diﬀerential operation, Opt. Express, 15, 17592 (2007)17） T. Kudo, T. Hirono, M. Nagasono and M. Yabashi:
Vacuum-compatible pulse selector for free-electron laser,Rev. Sci. Instrum., 80, 093301 (2009)18） R. Neutze, R. Wouts, D. van der Spoel, E. Weckert and J.
Hajdu: Potential for biomolecular imaging with femtosecond X-ray pulses, Nature, 406, 752 (2000)

［著
にし

の

よし

のり

西

野

吉

則 君

者 紹 介］

1996 年大阪大学大学院理学研究科物理学専攻
博士課程修了．博士 (理学)．96〜2001 年高輝度光科学研究センター，98〜2000 年ドイツ電子シンクロトロン (DESY) 客員研究員，01〜10 年理化学研究所，10 年〜北海道大学電子科学研究所
教授，現在に至る．

319

解説：ミニ特集
回折イメージング〜位相回復の新展開〜

光—回折強度からのレンズレスイメージング—
中

島 伸

治＊

＊静岡大学工学部 静岡県浜松市中区城北 3–5–1＊Faculty of Engineering, Shizuoka University, 3–5–1 Johoku, Naka-ku,
Hamamatsu, Shizuoka, Japan＊E-mail: tsnnaka@ipc.shizuoka.ac.jp

1.

はじめに

キーワード：位相回復 (phase retrieval)，コヒーレント回折イメージン
グ (coherent diffractive imaging)，波面再生 (wavefront reconstruction)，
逆問題 (inverse problem)，フーリエ変換 (Fourier transform)．
C 2011 SICE
JL 0005/11/5005–0320 

される．本解説では，上記で述べた X 線の波長程度の分解能を目指した回折イメージング法に焦点を絞って，特に，一般に光学顕微鏡は，物体からの散乱・回折光を凸レン

光学分野から見た歴史と最近の動向について述べる．また，ズで集めて結像させる原理を用いている．そのため，分解

著者が最近提案した開口アレイフィルターを用いたイメー

能は，結像レンズの有限な広がりや収差の影響を受け劣化

ジング法についても紹介する．

する．分解能を照明光の波長程度まで上げるために要求さ2.

て厳しくなる．特に，X 線領域 (波長：約 0.01〜10 nm) で回折イメージング法の進展：光から
X 線へ

は，波長程度まで分解能を上げることは非常に困難である．

レーザーの発明によりホログラフィーが注目され始めた

れる結像レンズの精度は，照明光の波長が短くなるにつれ

古くから，X 線結晶解析というレンズを用いず原子像を1960 年代，光学の分野では，ホログラフィーで不可欠な得る手法がある．この手法は，回折強度分布のみから回折

参照波を利用せず，回折強度分布のみから位相分布を直接，波の位相分布 (ある点の位相を基準とした位相差分布) を結回復することを目指した位相回復問題と呼ばれる研究が始

晶の周期性を利用して求め，その振幅と位相から得られる

まった7) ．当初は，理論的な研究が中心で目立った進展は

回折波面の複素振幅関数を逆フーリエ変換して原子像を得なかったが，コンピュータの発展に伴って 1970 年代にフーる回折イメージング法である．このような手法では，レン

リエ反復法2), 3) という強力な方法が提案されてから研究がズによる制約がないため，物体の分解能を照明 X 線の波長活発になった．フーリエ反復法は，物体からの回折波の強程度まで上げることができるが，周期性をもたない物体に度分布と物体の既知情報 (広がり範囲，非負条件など) を用は適用できない．

いて，回折波の位相分布を計算機アルゴリズムによって反光学の分野では，レンズを用いず周期性のない孤立物体

復的に求めていく手法である．そのため，常に解の収束性・像を得る手法が古くから研究されていて，たとえば，ホロ

唯一性の問題が伴う．しかし，1 枚の強度分布からの 2 次グラフィー1) やフーリエ反復法2), 3) などの方法が提案され

元位相回復が簡単に取り扱えるという特長があるため，瞬ている．これらはレーザー光源からの高可干渉性 (コヒーレ

く間に種々の問題に応用されるようになった．

ント) 光を用いた測定を前提としているため，従来の干渉性

一方，解析的な方法で位相回復する試みも 60 年代からいが悪い X 線源による測定に利用することが難しかった．最ろいろと行われてきたが7)〜10) ，理論的な困難さがあり 80

近，第三世代のシンクロトロン放射光源によって空間的に年代の中頃においても 1 次元位相回復すら満足に扱えない

干渉性が高い高強度な X 線が得られるようになったため，状態であった．80 年代初め日本でも位相回復問題の解説が光学分野で開発されたレンズレス回折イメージング法を利学会誌に掲載されている11), 12) ．ようやく，80 年代末に，解用した物体再生実験が世界的に行われようになった．

析関数の数学的性質を利用した 2 次元位相回復法によって，による像と異なり，物体を透過した波の振幅と位相の 2 つ

3 枚の回折強度分布からの物体再生が実験的に可能となった13) ．しかし，90 年代には入ってからは，回折イメージンの情報をもっている点である．特に，軽元素からなる物体

グに関連した位相回復問題は世界的にあまり活発に研究が(生体軟部組織や有機材料など) では，X 線の位相変化率が吸収率の約 1,000 倍大きいという性質があり4) ，従来の X行われなくなった．その理由は 2 つある．1 つは，1 枚の回回折イメージング法で得られる再生像の特徴は，レンズ

折強度分布からの 2 次元位相回復問題を解析的に解く試み

線吸収像ではわからない物体構造を，回折イメージングにが行き詰まったことによる．もう 1 つは，可視光の波長領域

よる再生物体位相から得られる利点がある．

ではレーザーというコヒーレントな光源があり，かつ高精

X 線の位相を測定する方法としては，シェアリング5) や度な干渉計測技術が確立されているため，ホログラフィー

タルボ・ロー6) などの干渉計を用いる手法があるが，通常等の参照波を用いる干渉計測法に対して，参照波を用いな

これらは物体の X 線投影像に対して用いる方法であるため，

い位相回復法の優位な点が見出せなかったことによる．

分解能が光源サイズおよび検出器のピクセルサイズで制限

320

90 年代の中頃から，高可干渉性・高強度な第三世代シン

計測と制御

第 50 巻

第5号

2011 年 5 月号

クロトロン放射 X 線が利用できるようになり，X 線イメー
ジング・計測分野の研究者が光学分野で提案されていた位相回復法を用いて実験を始めたことで再び注目を集めるようになった14), 15) ．これは，X 線の波長が可視光より遙かに短いために，精度の良い X 線結像系の作製が困難である
こと，さらに，高可干渉性 X 線とはいえレーザー光のように広い領域で干渉縞が得られないため，干渉計測やホログ
ラフィーの適用が制限される問題があり，それらに対して，物体の回折強度分布測定だけの簡単な測定系を用いる位相回復法の優位性が認識されるようになったからである．

2000 年代に入ってからは，種々の方法による X 線回折
イメージング実験が世界的に活発に行われるようになった．
おもに用いられている方法は，ホログラフィー法，フーリ
エ反復法，強度輸送方程式による方法であり，最近は反復法の欠点をホログラフィー法で補う融合法が提案されている．次章では，各方法について概説する．

3.

主要な X 線回折イメージング法

3.1 ホログラフィー法
放射光源からの高可干渉性 X 線であっても，レーザー光のような空間的に広い範囲の干渉性がないため，ホログラ
フィーで必要な参照波の与え方としては，一般に図 1 (a)，(b) の 2 つがある16) ．図 1 (a) は，オフアクシス型ホログラ
フィーというシステムで，X 線を透過させない平板に物体

図1

ホログラフィー法(a) オフアクシス型，(b) インライン型，(c) 角をもつ参照波による微分型．

用の穴と参照波用の小さな穴が開けてある．平板を背面からコヒーレント X 線波で照明し，物体面から十分離れた面

上記のように，従来のホログラフィー法は参照波またはで測定した場合，波の回折性によって物体面上の複素振幅

照明波用の穴の大きさが小さいほど再生分解能が良くなる．

関数 (波面の振幅と位相の 2 次元関数) のフーリエ変換が観

最近これを覆す方法として，図 1 (c) ような角 (かど) をも測面で得られる．実際に測定されるのは，物体と参照用のつ広がった参照波用開口を用いた微分ホログラフィー法が穴からの回折波の干渉強度分布であり，この強度データを提案された17), 18) ．このシステムでは，測定した干渉強度分

計算機で逆フーリエ変換して物体関数と参照波の相関関数

布に線形振幅フィルターをかけて逆フーリエ変換する (微が得られる．参照波用の穴が物体に対して十分小さく近似分) 処理によって，参照開口関数の微分関数と物体関数の的にデルタ関数と見なせる場合，相関関数は物体関数に相相関関数が得られる．図のような角をもつ開口関数を水平当する．しかし，実際には参照波用の穴は有限な広がりを

および垂直方向に微分すると角の部分が急峻なパルス関数

もつため，再生物体関数の分解能は X 線の波長ではなく，となるため，参照開口サイズより小さな穴を使用したオフ

その穴の物理的サイズで制限されてしまう．

アクシス型ホログラフィー法と同等な分解能が良い再生象

図 1 (b) は，インライン型ホログラフィーというシステとなる．この方法は，角のある種々の形状が利用できるのムで，物体面の後方にある小さな穴からの発散波を照明波

で，カーボンナノチューブなどのナノサイズ物体を参照波として用い，物体からの回折波と直接透過波との干渉強度源に用いて，さらに分解能を上げられる可能性がある．

を観測面で測定する．その干渉強度から直接透過波の強度成分を引いて，逆フーリエ変換することで物体関数が再生

3.2 フーリエ反復法
1972 年，Gerchberg と Saxton2) によって，参照波を利できる．この方式はガボール型1) であるため真の像と共役

用せず物体の回折強度分布とレンズによって結像された物

像が重なって再生されるが，光源の穴から物体までの距離

体強度分布から，計算機による反復アルゴリズムを用いてを波長に比べて十分長く設定すると，真の像を再生した際，回折面と像面の両方の位相分布を決定する方法が提案され共役像は広がった背景光となり影響はほとんどない．しかた．この方法では，回折面と像面の 2 枚の強度分布が必要

し，参照波である直接透過波の観測面での広がりが穴の大となるが，1978 年に Fienup3) によって 1 枚の回折強度分

きさで決まるため，再生像の分解能はオフアクシス型と同

布と物体の既知情報 (広がり範囲，非負条件など) から回折様に穴のサイズによって制限される．

面の位相を回復して物体を再生する方法が提案された．特計測と制御

第 50 巻 第 5 号

2011 年 5 月号

321

能であることが示されている21) ．

3.3 強度輸送方程式による方法
強度輸送方程式とは，自由空間での波の伝搬を表わすヘルムホルツ方程式を伝搬方向に対して近軸近似し，伝搬方向に垂直な面上の波の強度と位相の関係を表わした微分方程式である10) ．この方程式を解くと，強度の伝搬方向微分から波の位相を導く解析的な式が得られる22) ．まず，コヒー
レント波照明された物体からの回折波をある距離 z 離れた
観測面 1 とさらに微小距離 ∆z 離れた観測面 2 でそれぞれ独立に測定し，その強度分布を I1 ，I2 とおくと，波の強度の伝搬方向微分を ∂I/∂z ∼= (I2 − I1 )/∆z と近似的に求図2

Fienup のフーリエ反復アルゴリズム

めることができる．この関係を解析的な式に代入して観測面 1 の位相が求まる．観測面 1 の振幅 (I1 の平方根) と回復位相からなる複素振幅関数を数値的に逆伝搬させることに，1 次元より 2 次元の位相回復の方が安定した物体再生になることが初めて示されたため，多くの研究者の注目を集め，位相回復分野の発展の契機となった．

で，物体関数が再生できる．
この方法は，反復法と異なり解析的に位相が求まるが，ホログラフィー法や反復法などで再生可能な渦位相 (vortex

図 2 が Fienup 反復法の概略図である．計算機で発生さphase) 状態を再生できないことが理論的および計算機実験せた乱数を初期物体関数として用い，そのフーリエ変換をで明らかにされている23) ．渦位相とは，波の強度がゼロに計算する．求めたフーリエ関数の位相はそのままにして，なる点の周りで波面の等位相面が螺旋状に渦になっている

振幅分布を既知データ (観測強度の平方根である振幅分布)状態をいう．一般的に，物体から回折した波は伝搬していと入れ替えて逆フーリエ変換し，物体面に戻す．得られた

くうちに渦位相が自然に発生するため，あまり物体から距関数に対して，物体広がり範囲外のデータおよび負値デー離をとって観測面を設定できず，回折強度は物体の投影像

タをゼロの値と入れ替えるなどの処理をしてつぎのサイクに近い状態で測定する必要がある．そのため，再生物体のルの入力物体とする．このサイクルを物体面とフーリエ面

分解能は，光源サイズおよび検出器のピクセルサイズでほの条件を満足するまで繰り返す．観測データは，参照波をぼ決まってしまい，波長程度まで上げられずマイクロメー

利用しない回折強度分布のみでよいため測定が容易である

トルサイズの比較的大きな物体に対して用いられている22) ．

という利点があり，材料や生物試料などの多くの X 線回折

なお，渦位相を含む位相分布は，伝搬距離を変えて観測しイメージング実験で利用されている．しかし，反復アルゴ

た 3 枚以上の回折強度分布に Gerchberg-Saxton 法2) を基

リズムには，解の収束性，唯一性の問題が常に存在する．

礎にした反復法を適用して再生することが可能である24) ．

Fienup 法は，物体に非負 (すなわち位相変化のない) 条件

3.4 反復法とホログラフィー法との融合が使える場合，収束性が良くほとんど正しい解に収束する

最近，反復法を基本としながら，反復法で再生が困難なことが知られているが，位相分布をもつ複素振幅物体の場

複素振幅物体の位相回復を，ホログラフィーの干渉データ

合，アルゴリズムの収束性が悪く，真の解から外れた結果を補助的に利用しながら，ホログラフィー法のみの場合よに留まる停滞現象が生じやすくなる．なお，物体関数の有りも高分解能に物体再生を行う方法が提案されている．以限な広がりが別の測定などである程度正確にわかっている

下では，2 つの手法について述べる．

場合は，その広がり範囲と回折強度分布から Fienup 法で（ 1 ） オフアクシス型との融合25)

複素振幅物体関数を再生できることがある．しかし，たと

3.1 節で述べたように，ホログラフィーの再生分解能は，

えば円形のような中心対称広がりをもつ複素振幅物体は正参照波用穴のサイズで決定する．言い換えると，参照波の確な広がりを入れても物体再生できない場合がある19) ．

観測面上の広がり (干渉データの範囲) で再生物体の分解能

反復法による複素振幅物体の再生を改善する方法として，
タイコグラフィー (ptychography) 法

20)が決まる．観測面上で参照波の広がりを大きくするため参がある．この方法

照波用穴のサイズを小さくすると光量不足になり，観測面

では，物体広がりよりも狭い有限開口を透して物体を照明上で中心から離れた場所の干渉データの信号対雑音比が低し，物体面上での照明領域が重なるように有限開口を移動

下してしまう．そこで，図 1 (a) のシステムで，比較的大きさせながら開口の各位置で測定した複数枚の回折強度分布

な参照波用穴を用い，それを閉じた場合と開いた場合の観から，有限開口の広がりおよび位置の情報を用いて反復的測強度分布 (それぞれ IC , IO とおく) を観測する．その 2に物体関数を再生する．実際の X 線イメージング実験につの強度データの差 IO − IC の逆フーリエ変換から参照波

よって，物体が複素振幅関数であっても安定した再生が可

用穴の自己相関関数を分離し，そのフーリエ変換によって

322

計測と制御

第 50 巻

第5号

2011 年 5 月号

参照波の観測面上の強度分布 IR を得る．IC , IO , IR の 3 つ測定した回折強度分布から物体関数の各部分を再生し，そのデータと物体面上での物体および参照用穴の広がり条件

れらを合成して全体を再生することができる．

を用いて，フーリエ反復法により複素振幅物体を再生する．

4.

この場合，干渉強度 IO の情報を追加したことで，フーリ
エ反復法のみで再生した場合よりも速く安定な再生が可能となる．

開口アレイフィルターによる位相回復

X 線回折イメージングでは，物体の X 線によるダメージを抑えるため，1 回の回折強度分布測定から 2 次元物体をさらに，得られた物体関数を初期値として，物体関数の再生できることが望ましい．さらに，参照波を用いない非みの回折強度分布 IC に対して新たに単独でフーリエ反復

干渉型システムの方が簡便で，再生分解能が参照波によっ法を適用することで，大きな参照波用穴を用いた干渉強度て制限される問題を回避できる．著者は，最近，参照波に分布 IO からは得られない高分解能な物体が得られる．

よる干渉測定や反復法を利用せず，開口アレイフィルター（ 2 ） インライン型との融合26)を用いて測定した 1 枚の回折強度分布のみから，物体の厳フーリエ反復法は物体面と観測面の間のフーリエ変換関密な広がり条件なしで，2 次元複素振幅物体関数を再生す係を用いているが，観測面を物体面に近づけると，フレネる方法を提案した28), 29) ．この方法は，有限な広がりの物体

ル回折積分を用いなければならない．このとき，物体面と関数と観測面上の関数との間のフーリエ変換関係を基にし観測面の間で同様な反復アルゴリズムを適用できる．

た解析関数の数学的性質を利用している．以下では，この従来から，フレネル反復法の方が一般に収束性が良く解方法について概説し，光を用いた実証実験例について示す．

が安定していることは知られていたが27) ，物体が小さい場

図 4 が観測システムの概略図である．コヒーレント波に合，物体に近いフレネル回折も小さな分布となり，特に X

よって照明された物体からの回折波を等間隔に並んだ開口

線では回折強度分布の観測が困難となる．そこで，図 3 のアレイフィルターを通して観測面上で強度分布を測定する．

ように，X 線ビームをフレネルゾーンプレートを用いて細その際，各開口からの回折波の主要部分が干渉せず，それ

く絞り，ビームの焦点から少しずれた位置に物体を置くこ

ぞれ独立に強度分布を観測できるように距離 z, l を調整すとで，拡大されたフレネル回折分布が観測面で得られ測定る．また，開口間隔は強度変化を追従できるように，サンが容易になる．しかし，拡大させるためビームの焦点位置

プリング定理を基にして，物体面上の照明開口広がりの逆に近づけると照明波の波面の曲率が大きくなり，反復法で数から決定した．図 4 に示すように，観測強度分布から抽物体再生する際，照明波の物体面上の複素振幅分布が無視出した離散的な各開口の回折強度データ列 I1 (実線) は，物

できなくなる．そこで，事前に物体がない場合の照明波の複

体の回折複素振幅分布と開口関数の相関関数を，開口間隔

素振幅分布を，物体面と観測面間の距離を変えて複数枚測

でサンプリングした強度データとなる．さらに，強度デー

定した照明波の回折強度分布から Gerchberg-Saxton 法2)

タ列 I1 の位置から既知の微小距離 τ ずれた別の強度デーを基礎にした反復法で求めておく必要がある．その照明波

タ列 I2 = I1 + τ (破線) を観測強度から抽出する．強度デーの複素振幅を既知情報として，物体を置いた場合の 1 枚のタ列 I2 は，I1 と同様に物体の回折複素振幅分布と開口関

フレネル回折強度分布から物体の複素振幅分布を反復的に数の相関関数のサンプリング強度となるが，図 4 の幾何学分離して求める．ただし，物体広がりの情報が必要となる的関係からわかるように，物体関数が相対的に下にシフトので，フレネル回折強度分布の中心部分の照明波と物体波したときの強度データとなる．

の干渉強度の逆フレネル変換から物体形状をある程度求め，理論的に，強度データの比 I2 /I1 の対数は，相関関数のそれを初期条件として反復再生しながら徐々に物体広がり位相に対応した式 (近似的には位相の微分) で表わすことがを狭めて厳密な範囲を定める手法を用いている26) ．

できる29) ．観測面上の縦と横のそれぞれのデータ列から得この方法は，オフアクシス型照明と異なり，照明領域にられる式を数値的に解いて，相関関数の 2 次元位相分布がよって物体再生範囲を規定できるので，領域よりも広がった物体に対しても，照明領域を物体面上で走査させながら

図4
図 3 集光ビームによるフレネル回折強度からの反復再生

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

開口アレイフィルターによる位相回復

各実線 (または破線) の終点が観測強度データ列 I1 (または I2 = I1 + τ )の位置を表わし，始点が各データ列と物体との幾何学的関係を表わしている．

323

験系の場合，分解能は 50 µm となるので，実験結果から明らかなように，開口幅 0.1 mm よりも細かい物体のライン
周期 88.5 µm が解像されていることがわかる．図 5(c) は，
1 枚の凸レンズを使って物体を結像させたときの振幅分布
であり，比較のため図 5 (b) と同じ分解能で表示してある．
図 5 (b) と (c) を比較すると，本方法による (b) の方が全体的に雑音の影響が強く出ている．これは，(c) が CCD の一部分で受けた像であるのに対して，(b) は CCD カメラ
全面で測定した回折強度分布から再生しているため，信号対雑音比が相対的に悪くなっているからである．また，(c)の右側に縦に並んだ 2 点が (b) では再生されていないことがわかる．この結果は，各点の広がりが物体面上で分解能

50 µm よりも小さく，点からの波が物体のフーリエ変換面に相当する観測面上で大きく広がった分布となり，位相計
算の際に必要な回折強度の変化が得られなかったことが原因である．一方，(c) のレンズによる像は，分解能が低い場
合でもレンズを通過していく光波の位相は残るので，2 点の像は分解能に対応した広がった像として現れる．
本方法による回折イメージングでは分解能より小さな物
体は再生されないが，観測面を広げれば分解能を上げることができる．ただし，観測面を広げていくと物体面との間
図 5 実験結果の回折波のフーリエまたはフレネル変換関係が成立しなく(a) 測定強度分布，(b) 測定強度からの再生物体振幅，(c) レンズによる物なり，厳密な回折積分式を使用する必要がある．本方法は，

体振幅像．

そのような厳密な場合へ拡張することもできるので30) ，原
理的に波長程度まで分解能を上げることが可能である．

求まる．再生位相と観測振幅からなる相関関数を数値的に逆伝搬させ，物体面上で既知の開口関数の影響を取り除く
ことで物体関数が再生される．

5.

おわりに

近年，世界的に X 線自由電子レーザーの開発が進んでい図 5 が，光を用いた実証実験結果の一例である．半導体て，その高強度で高可干渉性 X 線パルスを異なる方向からレーザー (波長：0.635 µm) で，直径 0.4 mm の円形開口

単一分子に照射し得られる複数枚の 2 次元回折強度分布に，Fienup 反復法を適用して 3 次元分子像の再生を行う手法が検討されている31) ．分子の場合，X 線を散乱させる各原を通して物体 (1951 USAF テストネガパターンのグループ

3/ライン 4：11.30 本/mm) を照明した．市販の回折格子(5 本/mm) 2 枚を直交させて重ね合わせた正方形格子 (1
辺 0.1 mm：64 × 64 個) を開口アレイフィルターとして用いた．このアレイフィルターを取り付けた冷却 CCD カメ
ラ (1024 × 1024：1 ピクセルサイズ 13 × 13 µm) を物体から z = 530 mm 離れた位置に置いて，物体からの回折強度分布を測定した．このとき，フィルターから CCD 受光面までの距離は，各開口からの回折波の主要部分が干渉しないように選んだ距離 l = 6.7 mm を用いた．

子の電子密度は常に正値となるので，非負性という反復法の強力な拘束条件を利用できる．さらに最近，分子の 3 次元構造を，一方向からの X 線照射による 1 枚の遠方球面上
回折強度分布から，反復法によって再生する方法が提案されており32) ，回折イメージングの究極目標である単一分子
イメージングが実現可能な段階に入りつつある．
一方，テーブルトップで使用できる X 線装置による高分
解能な回折イメージング顕微鏡の開発33) も，生物や材料の

図 5 (a) が，1 回の測定で得た回折強度分布であり，この研究において重要である．特に，観察試料を透過した X 線

強度から 2 次元位相を計算し，逆フーリエ変換によって再の位相情報が再生できれば，通常の結像システムでは得ら生した物体関数の振幅分布が図 5 (b) である．ここでは示しれない物体情報を得ることができる．

ていないが，再生物体の位相は元の物体位相と同様なほぼ

今後，単一分子から生物・材料に至る幅広い領域で，レ一定の値となった．この方法は，従来の干渉計測法のようンズレス回折イメージング法が新しい X 線顕微鏡システムにフィルターの各開口の大きさによって分解能が制限されとして実用化されていくことを期待したい．

ることがなく，分解能は波長と開口数 (物体からフィルター（2011 年 1 月 20 日受付）

面への最大見込角度の sin 関数) の比で決定される．本実

324

計測と制御

第 50 巻

第5号

2011 年 5 月号

参

考

文 献

1） D. Gabor: A New Microscope Principle, Nature, 161,
777/778 (1948)
2） R.W. Gerchberg and W.O. Saxton: A practical algorithm
for the determination of phase from image and diﬀraction
plane pictures, Optik, 35, 237/246 (1972)
3） J.R. Fienup: Reconstruction of an object from the modulus
of its Fourier transform, Opt. Lett., 3, 27/29 (1978)
4） 百生 敦：位相コントラスト X 線イメージング，光学，29, 287/294(2000)
5） K. Iwata: Phase imaging and refractive index tompgraphy
for X-rays and visible rays, Progress in Optics (E. Wolf,ed.), 47, 393/432, North-Holland, Amsterdam (2005)
6） 百生 敦：タルボ・ロー干渉計による X 線位相イメージング，光学，38, 510/515 (2009)7） A. Walther: The question of phase retrieval in optics, Opt.
Acta, 10, 41/49 (1963)
8） R.E. Burge, M.A. Fiddy, A.H. Greenaway and G. Ross: The
Phase Problem, Proc. R. Soc. London Ser. A 350, 191/212(1976)
9） N. Nakajima and T. Asakura: Study of zero location bymeans of an exponential ﬁlter in the phase retrieval problem, Optik, 60, 289/305 (1982)
10） M.R. Teague: Deterministic phase retrieval: A Green’s
function solution, J. Opt. Soc. Am., 73, 1434/1441 (1983)
11） 前田純治，村田和美：光学における位相回復問題，光学，11,
230/240 (1982)
12） 中島伸治，朝倉利光：ヒルベルト変換による位相回復，光学，
11, 241/252 (1982)
13） N. Nakajima: Reconstruction of phase objects from experimental far-ﬁeld intensities by exponentail ﬁltering, Appl.
Opt. 29, 3369/3374 (1990)
14） K.A. Nugent, T.E. Gureyev, D.F. Cookson, D. Paganin and
Z. Barnea: Quantitative phase imaging using hard x rays,
Phys. Rev. Lett., 77, 2961/2964 (1996)
15） J. Miao, P. Charalambous, J. Kirz and D. Sayre: Extending
the methodology of X-ray crystallography to allow imaging of micrometre-sized non-crystalline specimens, Nature，
400, 342/344 (1999)16） A.P. Mancuso, Th. Gorniak, F. Staier, O.M. Yefanov, R.
Barth, C. Christophis, B. Reime, J. Gulden, A. Singer,
M.E. Pettit, Th. Nisius, Th. Wilhein, C. Gutt, G. Grübel,
N. Guerassimova, R. Treusch, J. Feldhaus, S. Eisebitt, E.
Weckert, M. Grunze, A. Rosenhahn and I.A. Vartanyants:
Coherent imaging of biological samples with femtosecond
pulses at the free-electron laser FLASH, New J. Phys. 12,035003-1/14 (2010)17） S.G. Podorov, K.M. Pavlov and D.M. Paganin: A noniterative reconstruction method for direct and unambiguous
coherent diﬀractive imaging, Opt. Express, 15, 9954/9962(2007)
18） D. Zhu, M. Guizar-Sicairos, B. Wu, A. Scherz, Y. Acremann, T. Tylisztzak, P. Fischer, N. Friedenberger, K.
Ollefs, M. Farle, J.R. Fienup and J. Stöhr: High-resolution
X-ray lensless imaging by diﬀerential holographic encoding,Phys. Rev. Lett., 105, 043901-1/4 (2010)19） M. Guizar-Sicairos and J.R. Fienup: Phase retrieval withtransverse translation diversity: a nonlinear optimization
approach, Opt. Express, 16, 7264/7278 (2008)20） H.M.L. Faulkner and J.M. Rodenburg: Movable aperturelensless transmission microscopy: a novel phase retrieval

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

algorithm, Phys. Rev. Lett., 93, 023903-1/4 (2004)
21） J.M. Rodenburg, A.C. Hurst, A.G. Cullis, B.R. Dobson,
F. Pfeiﬀer, O. Bunk, C. David, K. Jeﬁmovs and I. Johnson:
Hard x-ray lensless imaging of extended objects, Phys. Rev.
Lett., 98, 034801-1/4 (2007)
22） D. Paganin and K.A. Nugent: Noninterferometric phaseimaging with partially coherent light, Phys. Rev. Lett., 80,
2586/2589 (1998)
23） L.J. Allen and M.P. Oxley: Phase retrieval from series ofimages obtained by defocus variation, Opt. Commun., 199,
65/75 (2001)24） L.J. Allen, H.M.L. Faulkner, K.A. Nugent, M.P. Oxley and
D. Paganin: Phase retrieval from images in the presence ofﬁrst-order vortices, Phys. Rev. E, 63, 037602-1/4 (2001)
25） D. Zhu, B. Wu, R. Rick, J. Stöhr and A. Scherz: Phase retrieval in x-ray lensless holography by reference beam tuning, Opt. Lett., 34, 2604/2606 (2009)26） B. Abbey, K.A. Nugent, G.J. Williams, J.N. Clark, A.G.
Peele, M.A. Pfeifer, M. De Jonge and I. McNulty: Keyhole
coherent diﬀractive imaging, Nat. Phys., 4, 394/398 (2008)27） R. Rolleston and N. George: Stationary phase approximations in Fresnel-zone magnitude-only reconstructions, J.
Opt. Soc. Am. A, 4, 148/153 (1987)
28） N. Nakajima: Noniterative phase retrieval from a single
diﬀraction intensity pattern by use of an aperture array,
Phys. Rev. Lett., 98, 223901-1/4 (2007)29） N. Nakajima: Lensless coherent imaging by a deterministic phase retrieval method with an aperture-array ﬁlter, J.
Opt. Soc. Am. A, 25, 742/750 (2008)30） N. Nakajima: Phase retrieval from a high-numericalaperture intensity distribution by use of an aperture-arrayﬁlter, J. Opt. Soc. Am. A, 26, 2172/2180 (2009)31） K.J. Gaﬀney and H.N. Chapman: Imaging atomic structureand dynamics with ultrafast x-ray scattering, Science, 316,
1444/1448 (2007)32） K.S. Raines, S. Salha, R.L. Sandberg, H. Jiang, J.A. Rodriguez, B.P. Fahimian, H.C. Kapteyn, J. Du and J. Miao:
Three-dimensional structure determination from a single
view, Nature, 463, 214/217 (2010)33） A. Ravasio, D. Gauthier, F.R.N.C. Maia, M. Billon, J-P.
Caumes, D. Garzella, M. Géléoc, O. Gobert, J-F. Hergott,
A-M. Pena, H. Perez, B. Carré, E. Bourhis, J. Gierak, A.
Madouri, D. Mailly, B. Schiedt, M. Fajardo, J. Gautier, P.
Zeitoun, P.H. Bucksbaum, J. Hajdu and H. Merdji: Singleshot diﬀractive imaging with a table-top femtosecond
soft x-ray laser-harmonics source, Phys. Rev. Lett., 103,
028104-1/5 (2009)

［著
なか

じま

のぶ

はる

中

島

伸

治君

者 紹 介］

1983 年北海道大学院工学研究科博士課程修了(工学博士)．同年，静岡大学工業短期大学部講師，
85 年同助教授，94 年静岡大学工学部助教授，99年同教授．光の位相回復とその応用に関する研究に従事．応用物理学会，日本光学会，アメリカ光学会の会員．

325

解説：ミニ特集
回折イメージング〜位相回復の新展開〜

電子回折イメージング—入射ビームの角度拡がり—
郷

原 一

寿＊

＊北海道大学大学院 工学研究院 北海道札幌市北区北 13 条西 8 丁目
＊Research Faculty of Engineering, Hokkaido University, Kita 13, Nishi 8,
Kita-ku, Sapporo, Hokkaido, Japan＊E-mail: gohara@eng.hokudai.ac.jp

1.

はじめに

回折パターンの計測データをもとに，計算機を用いた数

キーワード：回折イメージング (diffractive imaging)，位相回復 (phaseretrieval)，電子 (electron)，入射ビーム角度拡がり (incident beam divergence)，横コヒーレンス (transversal coherence)．
C 2011 SICE
JL 0005/11/5005–0326 

2.

入射ビームの角度拡がり「入射ビームの角度拡がり」とは，図 1 に示すように，物値計算により位相回復を行うことで，レンズを使わずに物

体 (Object) に対して異なる角度で平面波が入射すること

質のイメージングが可能であることは，1999 年に初めて実を指す (収束角と表現することもある)．回折イメージング

験的に示された1) ．これを分岐点として，回折パターンをもにおいて，入射ビームの角度拡がりが重要な課題の 1 つでとに，結晶性を前提とせず，対物レンズも必要としない新たあることは，早くから指摘されている12), 27) ．

なイメージング手法の開拓が始まり，これまでの約 10 年間

平面波が物体 f (r) によって 1 回のみ散乱され (第 1 ボで飛躍的に進展した．特に，X 線領域においては，次世代高ルン近似)，かつ散乱された波を十分遠方で観測する場合(フラウンフォーファー回折)，検出面で計測される回折強度 I0 (k) は，物体 f (r) のフーリエ変換 F {f (r)} の振幅|F (k)| の 2 乗に比例する．すなわち，簡単のため比例係数
を 1 とすれば，輝度高干渉性光源プロジェクトの中心的手法として大きな期待が寄せられている1)〜11) ．この手法は “回折イメージン
グ (Diﬀractive Imaging)” 又は “回折顕微法 (Diﬀraction

Microscopy)” と呼ばれている．本稿では，回折イメージン
グで統一する．
電子はド・ブロイ波長で特徴付けられる波の性質をもつ
ことから，新たなイメージング手法の 1 つとして，電子顕I0 (k) = |F (k)|2 = |F {f (r)}|2 .

(1)

である．

微鏡の領域においても回折イメージングに関する研究が行

また，異なる入射角をもつ平面波の間に干渉性のない場合われている．電子顕微鏡での実験的検証は，2002 年，透過(incoherent) には，検出面で計測される回折強度 Iconv (k)は，入射ビームの角度拡がりを表わす関数 h によるコンボ
リューションとして，以下のように表現できる28) : ∞I0 (k − k  )h(k  )dk  .
(2)
Iconv (k) =電子顕微鏡 (TEM) を用いて対物レンズの後焦点面で検出した回折パターンをもとに，位相回復を行った研究が最初である12) ．回折イメージングの条件から外れているため，得られた結果は通常の TEM 像と比較して分解能は劣っているが，電子顕微鏡においても，回折パターンをもとにして，位相回復を行うことにより実像が得られることを実験的に初めて示した報告である．その後，複数のグループが活発に研究を展開している13)〜19) ．著者らのグループもその 1 つであり，回折イメージングの基礎と応用の両面から−∞(1) 式および (2) 式によって，入射ビームの角度拡がり hがある場合に，物体 f (r) の散乱によって検出面で観測される回折強度 Iconv (k) が計算できる．なお，電子による散乱の場合には，物体 f (r) は静電ポテンシャルである．

研究を進めている20)〜25) ．
回折イメージングはフーリエ変換を基本としていることから，フーリエ変換の関係から外れる要素を特定して，それらを解消することが重要である．本稿では，電子顕微鏡に回折イメージングを適用する際の 1 つの問題として，おもに入射ビームの角度拡がりを取り上げ，その影響と解消法，解消法の有効性の検証について紹介する21), 23), 24) ．電子顕微鏡分野での回折イメージングに関する全般的な研究の進展について，さらに興味のある方は，著者らの別解説26) を参考にしていただければ幸いである．

326

図1

入射ビームの角度拡がりと回折

計測と制御

第 50 巻

第5号

2011 年 5 月号

3.

コンボリューション

数)．右の回復物体を見るとスリットの位置がシフトしているが，モデルのスリット幅 50 ピクセルを忠実に再現して

図 2 に示すように，フーリエ反復位相回復法は，物体にいる．下段は，角度拡がりがある場合 (σ = 1，標準偏差 1

対するオブジェクト拘束条件 (Object Constraint) と，回ピクセルに対応) を示している．左の回折パターンの拡大

折パターンによって得られた回折波の振幅をフーリエ拘束

図を上段と比べると，角度拡がりによるコンボリューショ

条件 (Fourier Constraint) とし，フーリエ変換 (FT ) と逆ンの結果，谷底の位置がゼロに落ちていないことがわかる．

フーリエ変換 (FT−1) を逐次的に交互に繰り返すことによっ

右の回復物体を見るとスリットの位置がシフトしているが，て位相を得る手法である29)〜31) ．なお，オブジェクト拘束条

モデルのスリット幅 50 ピクセルは再現している．しかし，件は，具体的なアルゴリズムで異なる．

スリットの形がモデルの矩形とは異なっており，上側の領域が不規則に削られ，変形していることがわかる．
図 4 に，2 ピクセル離れた幅 2 ピクセルのスリット 2 本
を 1 組として，さらに 120 ピクセル離してもう 1 組を配置した合計 4 本のスリットからなる複スリットを物体 f のモ
デルとして，図 3 と同様な条件で位相回復を行ったシミュ
レーション結果を示す．全体は 1024 ピクセルであり，サポートは 200 ピクセルとした．

図2

フーリエ反復位相回復法

図 3 は，回折イメージングに対する入射ビームの角度拡がりの影響の例をシミュレーションによって示している．

1024 ピクセルの中心に幅 50 ピクセルの単スリットを物体
f のモデルとし，角度拡がりをガウス分布として，(1) 式
および (2) 式によって回折強度 |F |2 を計算し，これを元に，図 2 の位相回復法を適用した．フーリエ拘束条件とし上段は，角度拡がりがゼロの場合 (σ = 0) を示しており，
左の回折パターンには大きな山・谷の形状が見られる．1
つの山を拡大してみると (拡大図参照)，その中にさらに微細な山・谷の構造がある．これは，2 ピクセルと 120 ピク
セルに対応した干渉パターンが重なって生じている (正確にはそれぞれを特徴付ける関数の積)．右の回復物体を見ると全体がシフトしているが，元の形状を忠実に再現している．下段は，角度拡がりがある場合を示している．左の回て，左側の回折パターンを用い，オブジェクト拘束条件と折パターンの拡大図を上段と比べると，角度拡がりによるして，実関数・非負性および中心から左右に 100 ピクセル

コンボリューションの結果，図 3 と同様に谷底の位置が落が物体を取り囲む領域 (サポート) として，それ以外の値をち込み切れていない．右の回復物体を見ると左側の 1 組のゼロとし，500 回までは HIO (Hybrid Input-Output)，そスリットが右側の 1 組に比べて大幅に低い値となっている．

れ以降 1000 回までは ER (Error Reduction) と呼ばれる更さらに拡大図を見ると最小幅のスリットの値も僅かに小さ

新アルゴリズムを適用して，右側の物体が得られた．

くなっており，4 本とも元の矩形とは異なっている．

上段は，角度拡がりがゼロの場合 (σ = 0) を示している．

以上から，入射ビームの角度拡がりによって，オリジナ

左の回折パターンを見ると干渉に特有な山・谷の特徴的な

ルの物体は変形され，大きなサイズの構造が小さなサイズ

形状が明瞭に見られる (矩形関数のフーリエ変換は sinc 関の構造よりもより強く影響を受けることがわかる．

図 3 単スリットの位相回復

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

図 4 複スリットの位相回復

327

4.

て量子ノイズが重畳した場合の位相回復への影響を示すシ

デコンボリューション

ミュレーションである．ここでは，入射ビームの角度拡が(2) 式の I0 は，デコンボリューションによって解析的に，I0 = F −1 {F {Iconv }/F {h}} と求めることができる．

りをガウス分布とし，量子ノイズをポアソン分布とし，総カウント数を 107 個とした．

しかし，実際の実験においては，必ずノイズが含まれてお上段の角度拡がりがゼロ (σ = 0) の回折パターン (左) をり，その場合のデコンボリューション法が複数提案されて見ると，強いコントラストが見られる．強度の弱い領域をいる32) ．著者らは種々の方法について考察し，逐次的なデ見るとノイズが重畳していることがわかる (右上の拡大図

コンボリューション法を位相回復法に組み込む方法が有効

参照)．この回折パターンをもとに位相回復によって得られであることを検証した24) ．

た右の物体を見ると，ノイズが少し重畳しているが，元のモデル物体とほぼ同一のイメージが得られている．下段の角度拡がりがある場合 (σ = 1) の回折パターンを見ると，
角度拡がりがない場合と比べて，ぼやけたコントラストを示していることがよくわかる (右上の拡大図参照)．これを元に位相回復によって得られた右の物体も，全体的にぼやけており，部分的にしか再構成できていない．

図 5 位相回復法に組み込んだデコンボリューション
図 5 は，その模式図を示している．図 2 と比べると，フーリエ拘束条件を適用する直前にデコンボリューションを行う
ことが特徴である．位相回復法に組み込んでデコンボリュー
ションを実行できるものは，Richardson-Lucy (RL)，Max-

imum Entropy Method (MEM)，およびそれらの変形手法があるが，ここでは，RL を用いた方法とその効果について説明する．RL は以下の式で逐次的に更新する方法である．


InoiseI(t + 1) = I(t) · h ∗h ∗ I(t)図 7 は (3) 式を使って，図 6 の σ = 1 の回折パターン
を用い，図 5 による位相回復過程を示している．グラフはHIO を 5000 回後，さらに ER を 5000 回繰り返したときの R ファクター ((3) 式の収束の度合いを評価する値) の変化を示している．HIO では増減を激しく繰り返し，ER に切り替えた直後で急激に下がり，その後は単調に徐々に減
少していき，ほぼ飽和に達している．これは一般的によく知られている，HIO‐ER アルゴリズムの典型的な変化を示している．グラフの下に 3000，6000，10000，それぞれの繰り返し回数時点での回折パターンの拡大図を示した．これを見ると，ぼやけたノイズの多いパターンから，徐々に(3)

ここで，I(t) は t 回目の回折強度，Inoise は入射ビームの角度拡がり h でコンボリューションされ，さらにノイズ

コントラストのハッキリしたパターンに変化していることがよくわかる．グラフ中の 10000 回目の回復物体を見ると，
図 6 の σ = 1 からは大きく改善して，σ = 0 にほぼ近い
物体イメージが得られている．

を含んだ回折強度である．また記号 ∗ はコンボリューショ
ンを表わす．
図 6 は，回折パターンに入射ビームの角度拡がりに加え

図 6 角度拡がりと量子ノイズが重畳した場合の位相回復

328

図 7 デコンボリューションを考慮した位相回復過程

計測と制御

第 50 巻

第5号

2011 年 5 月号

5.

図 9 に，実験データをもとに飽和領域のない回折パター

実験データへの適用

この章では，前章のデコンボリューションの有効性を，実際の電子顕微鏡を使った実験データによって検証した例を示す．実験は，透過電子顕微鏡 (JEM2010F，JEOL 社製)
を，高分散回折モードと呼ばれる対物レンズを切ったモードで使用した．このモードでは入射ビームの並行性を上げることができ，角度拡がりを小さく抑えることができるが，散乱角が小さくなるために，原子スケールの高分解能観察には使用できない．しかし，試料と検出面の距離を長くとれるため，電子顕微鏡において対物レンズレスで，サブミ
クロン程度の試料サイズでも容易に回折パターンが取得でき，回折イメージングの基本的実験が可能であるというメ
リットがある．ここでは，このメリットを生かした実験を紹介する．
図 8 は，試料に Ta/SiC 回折格子を用いて，加速電圧200 kV，カメラ長 40 m，照射時間 1s で，1 ピクセル 25 µmのイメージングプレート (IP) によって検出された回折パ
ターンである．上が回折パターン，真中に中心の拡大図，下に拡大図の水平ラインのプロファイルを示している．試料の SEM 像を右上に挿入した．上の回折パターンには，水平方向に特徴的な干渉パターンが見られる．拡大図を見ると，一様な黒い領域があり，ラインプロファイルにより，この領域の強度が飽和していることがわかる．この飽和領域は同じ IP を 2 回読み取り，1 回目と 2 回目の読み取り値を合成した回折パターンから，1024 × 1024 ピクセルを切り出し，位相回復を行った．

ンをフーリエ拘束条件としてそのまま用い，図 2 に示した
位相回復を行った結果 (上段) と，図 5 で示した組み込み型のデコンボリューションを適用して位相回復を行った結果(下段) を示している．デコンボリューションのための角度
拡がりは，実験により得られたダイレクトビームの形状を規格化して使用した．両方の位相回復とも，オブジェクト
拘束条件として，実関数・非負性および円形のサポート領
域，HIO-ER それぞれ 1000 回行った．
上段の回折パターンをそのまま用いて位相回復を行った回復像を見ると，中心付近に 2 本の棒状の強い強度領域がある．2 本とも周辺が不規則にかすれており，輪郭が不明
瞭であるように見える．物理的な再構成サイズ ∆x は，横コヒーレンス長によって制限されており，波長を λ，拡がり角度を α として，∆x ∼ λ/(2α) で見積もることができる24) ．ここで，加速電圧 200 kV に対するド・ブロイ波長，λ = 0.0025 nm，ダイレクトビームの測定から求めた拡がり角度 α = 5×10−6 rad を代入すると，∆x ∼ 250 nm とな
り，再構成サイズについては実験結果とほぼ対応している．
下段の最終的に得られた回折パターンを，上段の元の回
折パターン (初期パターン) と比べると，細くシャープな回
折ピークであり，両者に明確な違いがあり，デコンボリュー
ションの効果がよくわかる．これに対応した右側の回復像を見ると，4 本のスリットが再構成されており，上段の回復像と比較すると，再構成領域が拡大して本数が増えていること，左上・右下にわずかなかすれが残っているが，1 本

1 本のスリットの輪郭がハッキリしていることがわかる．4
本のスリットの下に配置した SEM 像と比較すると，スリット幅およびスリット間隔が定量的に一致している．

図 8 実験に用いた Ta/SiC 回折格子からの回折パターン

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

図 9 実験データをフーリエ拘束条件とした位相回復の結果 (上段) と組み込み型のデコンボリューションを適用した
位相回復の結果 (下段)

329

6.

る場合の具体的な例を示している．すなわち，従来得られ電子回折顕微鏡

著者らは日立製作所との共同研究により，従来の電子顕微鏡の性能を活かしつつ，結像のための対物レンズを使用
せずに回折パターンを計測可能とする，低加速電圧で動作
する回折イメージング専用の電子回折顕微鏡の試作機を開発した23) ．図 10 は試作機の概観を示している．試作機は，ている手法で荒い画像 (コースイメージ) を観察し，さらに高い分解能を求めたい場合は，回折イメージングで詳細な画像 (ファインイメージ) を観察するという相補的な使用が想定できる．この場合，STEM，収差補正機付 TEM など
もコースイメージの候補となる．従来の既存技術と回折イ
メージングをシームレスに接続すれば，これまでに得られSEM(走査電子顕微鏡) をベースとしており，SEM モードていない高分解能化を実現できる可能性を示している．

と回折イメージングモードに電子光学系を切り換えること7.

が可能である．

討論

この試作機によって，図 5 のデコンボリューション法も計測技術としての電子顕微鏡の重要な特徴の 1 つは，局取り入れた回折イメージングの実験を，共同研究の一環と所的な構造を高い空間分解能でイメージングできる点にあして行った23), 24) ．図 11 に，実験によって得られたイメー

る．最新技術のレンズ収差補正器の適用によって，たとえば

ジの例を示す．最初に，SEM モードによって，所望の領加速電圧 300 kV で空間分解能 0.05 nm に達しており，サ域を順次拡大して観察領域を絞り込み (下図)，最大倍率にブオングストロームの分解能が得られる時代に突入してい達した段階で，さらに拡大したい領域 (白丸) に並行ビームを照射して，回折パターンを計測し，それをもとに回折イる33), 34) ．しかし，300 kV のド・ブロイ波長は 0.002 nm であり，現在到達している空間分解能は波長の約 25 倍であメージングにより実像を再構成した (上図)．再構成された

る．原理的な波長 (回折) 限界とはまだまだ大きなギャップ

イメージから，この試料は 4 層のマルチウォールカーボンがある．

ナノチューブであり，外径，内径，ウォール間隔はそれぞ

回折イメージングは原理的に対物レンズを必要としない

れ 4.10 nm，2.10 nm，0.34 nm であることがわかった．

イメージング方法である．レンズ性能を上げることで，高

この例は 30 kV であるが，さらに低加速の 10 kV でも同分解能を追求する従来の電子顕微鏡とは異なる技術を必要

様な分解能でイメージング可能であることが明らかとなっとしている．前章で示したような，この方法に特化した研てきており，ノイズ低減，原子 1 つ 1 つを識別できる高分究・技術開発は開始されたばかりであり，これまでは，す解能化など，さらなる改善を進めている．この試作機では，

でに見えているものを，この方法でも見えることを検証しSEM によって初期画像およびサポート領域を得ることがている段階であるといえる．今後は，今までに見えていな

でき，回折イメージングの重要なオブジェクト拘束条件といものをイメージングすることが始まる．

して，実験データをそのままアルゴリズムに組み込むこと

回折イメージングの分解能に対する原理的限界は回折限ができる．回折イメージングを電子顕微鏡として装置化す界である．長期的な視点に立てば，技術の発展は原理的限界に向かって収束する．今後の展開が大いに期待される．

図 10

330

電子回折顕微鏡のプロトタイプ

図 11 マルチウォールカーボンナノチューブ

計測と制御

第 50 巻

第5号

2011 年 5 月号

謝辞：本稿に記した一部は，JST 育成研究プロジェクトとして実施されたことをここに記し，関係者の皆様に対する感謝に代えさせていただきます．
（2011 年 2 月 7 日受付）
参

考

文 献

1） J. Miao, P. Charalambous, J. Kirz and D. Sayre: Nature,400, 342 (1999)
2） J. Miao, T. Ishikawa, B. Johnson, E.H. Anderson, B. Lai
and K.O. Hodgson: Phys. Rev. Lett., 89, 88303 (2002)
3） J. Miao, K.O. Hodgson, T. Ishikawa, C.A. Larabell, M.A.
LeGros and Y. Nishino: PNAS, 100–1, 110 (2003)
4） D. Shapiro, P. Thibault, T. Beetz, V. Elser, M. Howells, C.
Jacobsen, J. Kirz, E. Lima, H. Miao, A.M. Neiman and D.
Sayre: PNAS, 102–43, 15343 (2005)
5） 西野吉則，石川哲也：放射光，19–1, 3 (2006)
6） K.J. Gaﬀney and H.N. Chapman: Science, 316, 1444 (2007)
7） R.L. Sandberg, A. Paul, D.A. Raymondson, S. Hadrich,
D.M. Gaudiosi, J. Holtsnider, R.I. Tobey, O. Cohen, M.M.
Murnane, H.C. Kapteyn, H.C. Kapteyn, C. Song, J. Miao,
Y. Liu and F. Salmassi: Phys. Rev. Lett., 99, 098103 (2007)8） 西野吉則：顕微鏡，44–1, 24 (2009)9） Y. Nishino, Y. Takahashi, N. Imamoto, T. Ishikawa and K.
Maeshima: Phys. Rev. Lett., 102, 018101 (2009)10） K.S. Raines, S. Salha, R.L. Sandberg, H. Jiang, J.A. Rodriguez, B.P. Fahimian, H.C. Kapteyn, J. Du and J. Miao:
Nature, 463, 14 (2010)11） M.M. Seibert, et al.: Nature, 470, 78 (2011)12） U. Weierstall, Q. Chen, J.C.H. Spence, M.R. Howells, M.
Isaacson and R.R. Panepucci: Ultramicroscopy, 90, 171(2002)
13） J.M. Zuo, I. Vartanyants, M. Gao, R. Zhang and L.A. Nagahara: Science, 300, 1419 (2003)
14） S. Morishita, J. Yamasaki, K. Nakamura, T. Kato and N.
Tanaka: Appl. Phys. Lett., 93, 183103 (2008)
15） W.J. Huang, J.M. Zuo, B. Jiang, K.W. Kwon and M. Shim:
Nature Phys., 5, 129 (2009)16） R. Dronyak, K.S. Liang, Y.P. Stetsko, Ting-Kuo Lee, ChiKai Feng, Jin-Sheng Tsai and Fu-Rong Chen: Appl. Phys.
Lett., 95, 111908 (2009)17） L.D. Caro, E. Carlino, G. Caputo, P.D. Cozzoli and C. Giannini: Nature Nanotechnology, 5, 360 (2010)18） J. Rodenburg and Maiden: Microscopy and Analysis, 24–1,
17 (2010)

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

19） F. Hue, J.M. Rodenburg, A.M. Maiden, F. Sweeney and
P.A. Midgley: Phys. Rev. B, 82, 121415 (2010)20） H. Shioya and K. Gohara: Opt. Commun., 266, 88 (2006)21） O. Kamimura, K. Kawahara, T. Doi, T. Dobashi, T. Abe
and K. Gohara: Appl. Phys. Lett., 92, 024106 (2008)22） H. Shioya and K. Gohara: J. Opt. Soc. Am. A, 25–11, 2846
(2008)
23） O. Kamimura, T. Dobashi, K. Kawahara, T. Abe and K.
Gohara: Ultramicroscopy, 110, 130 (2010)24） K. Kawahara, K. Gohara, Y. Maehara, T. Dobashi and O.
Kamimura: Phys. Rev. B, 81, 081404(R) (2010)25） H. Shioya, Y. Maehara and K. Gohara: J. Opt. Soc. Am.
A, 27–5, 1214 (2010)26） 郷原一寿，上村 理：顕微鏡，44–1, 69 (2009)27） I.A. Vartanyants and I.K. Robinson and J. Phys.: Condens.
Matter, 13, 10593 (2001)28） W.E. McBride, D.J.H. Cockayne and C.M. Goringe: Ultramicroscopy, 76, 115 (1999)
29） R.W. Gerchberg and W.O. Saxton: Optik, 35–2, 237 (1972)30） J.R. Fienup: Opt. Lett., 3–1, 27 (1978)
31） J.R. Fienup: Appl. Opt., 21–15, 2758 (1982)
32） P.A. Jansson: Deconvolution of images and spectra, Academic Press, New York (1997)
33） R. Erni, M.D. Rossell, C. Kisielowski and U. Dahmen:
Phys. Rev. Lett., 102, 096101 (2009)34） H. Sawada, Y. Tanishiro, N. Ohashi, T. Tomita F.
Hosokawa, T. Kaneyama, Y. Kondo and K. Takayanagi:
J. Electron Microsc., 58, 357 (2009)［著者 紹 介］

ごう

はら

かず

とし

郷

原

一

寿 君（正会員）
1982 年名古屋大学大学院工学研究科修士課程修了，工学博士．同大学助手，講師，中部大学助教授，北海道大学助教授，現在，工学研究院応用物理学部門教授．応用物理学，生物物理学，非線形ダイナミクスを専門とする．日本物理学会，日本顕微鏡学会，日本生物物理学会，日本神経科学
学会，北米神経科学学会，日本神経回路学会など
の会員．

331

解説：ミニ特集
回折イメージング〜位相回復の新展開〜

位相回復—計算アルゴリズム—
塩

谷 浩

之＊・郷 原一 寿＊＊＊室蘭工業大学 工学部 北海道室蘭市水元町 27–1＊＊北海道大学大学院 工学研究院 北海道札幌市北区北 13 条西 8 丁目
＊Faculty of Engineering, Muroran Institute of Technology, 27–1 Mizumotocho, Muroran, Hokkaido, Japan＊＊Research Faculty of Engineering, Hokkaido University, Kita 13, Nishi 8,
Kita-ku, Sapporo, Hokkaido, Japan＊E-mail: shioya@csse.muroran-it.ac.jp

1.

キーワード：位相回復 (phase retrieval)，逆問題 (inverse problem)，GS
アルゴリズム (GS-algorithm)，情報量 (information measure)．
C 2011 SICE
JL 0005/11/5005–0332 換である．フーリエ変換が 1 対 1 であることから，背景条

まえがき

件は素直であるが，与えられているのは振幅の自乗 |g(y)|2

波の計測において，強度を観測できるが位相を観測できであるので，悪条件問題と言える．もし位相が再構築できない状況は位相問題 と呼ばれる．たとえば，X 線で物質ると，フーリエ強度との逆フーリエ変換によって f を求めの構造を見る場合においては，試料で散乱した波の強度をることができる．位相を回復することは，未知の実像を求測定することはできるが，位相は失われている．位相問題

めることと等価である．

1)は，結晶のような原子レベルの小さなものを対象にする各

本特集の題目「回折イメージング」とは，回折強度から種の顕微鏡から，電波望遠鏡など天体を扱う世の中のすべ実像を求め，その像を顕微鏡像として扱うことであり，回てのスケールに関わっている．これまで，回折強度からの折パターンのもつ最大の分解能の実像が望まれる．それに実像のイメージングには，結晶のような周期的な構造をもは，フーリエ強度からの位相回復ができることが，この手つ場合が常識であったが，1999 年の Miao らの実験2) 以来，法の根本的条件となる．位相回復に対して，最初にその可回折強度を用いた非周期構造のイメージングの研究が精力能性を示唆したのは Sayre3) である．それはシャノン4) の的に進められている．

サンプリング定理に基づいて指摘されており，50 年代に位相問題は，逆問題の 1 つとして表現される．逆問題のおいて，位相回復を不完全データに対する情報科学の問題

一般的形式から簡単に述べると，未知の関数 f (x) があり，
X 上で定義されているとする．得られているのは，X × Y
上の核関数 K(x, y) による以下の積分である．
g(y) =
K(x, y)f (x)dx(1)と捉えている先見性が感じられる．Sayre から 20 年後にGerchberg と Saxton5) は，フーリエ変換の繰り返しによ問題は，g(y) が得られている条件で f (x) を求めることでない．本解説では，位相回復の計算アルゴリズムとそれにXるアルゴリズムを考案し，それは位相回復の現実的な解法として現在も輝いている．これは単純な正逆のフーリエ変換の繰り返しなのであるが，その本当の姿を捉えきれていある．核関数による変換が 1 対 1 であるときには自明とな付随する理論研究を紹介し，回折イメージングのための計るが，現実の問題としては，たとえば g(y) に加法的ノイズ

算数理の今後の方向性を展望する．

が入ることなどの条件では，解 f があるかどうかは保証さ2.

れない．図 1 においては，X, Y の関数空間 LX , LY にお

基本設定

いて，観測で得られている gobs からの f の探索を表わす．

位相回復の対象となるのは，文字通り “位相” であるが，位相回復をこの形式でいうと，実像の空間および逆空間

それによってわかる対象物が真の目的である．それを定義が複素ベクトル空間，核関数 K による積分がフーリエ変域上で複素数値をとる関数とする．信号であれば時間軸上の関数であり，構造解析であれば原子の配置と結合に対応
する．2 次元領域上の画像とすると理解しやすい．それを一般に実像 (Object) と呼び，その定義空間を実空間 (Ob-ject domain) という．有限体積をもつ関数を実像とし，そのフーリエ変換は，その核関数の定義域の一部となる逆空間 (Fourier domain) 上の複素数値をとる関数となる．実空間 X 上の実像を関数 ρ(x) とし，そのフーリエ変換は，逆空間 K 上の関数となる．実用的には，実空間は離散である場合が多く，その場合は離散フーリエ変換がよく用い
図1332

観測された gobs からの f の再構築

られる．そこで，実空間 X を 2 次元の離散空間 N × N とし，逆空間も同様とする．離散フーリエ変換によって与え計測と制御

第 50 巻

第5号

2011 年 5 月号

られる ρ のフーリエ変換は，以下のように与えられる．

F (k) =



ρ(x) exp{−i2πx · k/N }.

(2)

x∈X

F は振幅と位相によって |F | exp(iΦ) と表現される．位相問題としては，フーリエ強度 |F |2 は得られているが，位相 Φ は失われている．以下，強度に対して観測を表わすた
めに |F |2 を |Fobs |2 と表現する．図 2 に示すように，なんらかの方法で位相を推定したとしよう．その推定位相をΦ̂ とすると，以下の逆フーリエ変換によって，再構築された実像 ρ̂ を得ることができる．
1  
ρ̂(x) = 2
F (k) exp{i2πx · k/N },
(3)
N
k∈K

ただし F  = Fobs exp[iΦ] とする．

図 3 繰り返しフーリエ変換と拘束条件によって構成される GS ダイヤグラム
像関数 ρ を用意する．そのフーリエ変換 F に対して，条件として与えられている振幅 Fobs と入れ替える．すなわ
ち，F = |F | exp(iΦ) とすると，F  = Fobs exp(iΦ) を，
逆空間拘束条件を適用した関数とする．つぎに，F  を逆フーリエ変換したものを ρ = |ρ| exp(iθ) とし，実空間拘束条件である ρnorm を同様に適用することで，更新された
実像関数 ρnorm exp(iθ) が得られる．更新数を n で表わすと，ρn → Fn → Fn → ρn → ρn+1 を繰り返し行い，ρnと ρn+1 が一致したときに，両空間の拘束条件を満たす実像が得られたことになる．図 3 は GS ダイアグラムと呼ばれ，この繰り返し操作は GS アルゴリズムと呼ばれている．
図 2 推定位相と再構築実像
しかしながら，このままでは逆空間上の位相関数を決める手立てがまったくない．実用においては，実像について
何か事前情報があることが多い．信号や画像であれば，画像空間や時間軸のサポート（台）であったりする．これは一般的には，値がない領域が “ある程度” わかっていることに相当する．そのような実像に関する条件を実空間拘束条件(Object-domain constraint) といい，位相回復ではフーリ
エ強度が条件として与えられるが，そのような逆空間において与えられる条件を，逆空間拘束条件 (Fourier-domain

constraint) という．逆空間拘束条件として与えられる未知関数のフーリエ変換の強度，および実像の事前情報に相当する実空間拘束条件を満たす関数を求めることと表現でき設定される実空間における拘束条件として，実数や正値性などがある．それらは物理的に対応がつく条件であることから，有用な拘束条件の設定と言える．位相回復のアルゴ
リズムを考えるときに，解探索の途中の仮説実像について，
実空間 X は，拘束条件を満たす領域と満たさない領域に分けられる．領域によって更新をそれぞれ行うことが，以下に述べるアルゴリズムの有用さを支える．

3.2 ER と HIO
GS から約 10 年後，Gerchberg と Saxton の方法に基づいて，Fienup は，GS ダイヤグラムにおける ρ から ρ への実像の更新が逆空間拘束条件を満たす最小化における最急降下法に相当することを，以下の導関数から示した6) ．

2∂  
|F (k)| − |F  (k)| = 2(ρ − ρ ),
∂ρ

(4)

k∈K

る．つまり，制約条件を満たす関数を推定する問題である．

3.

ただし，ここで用いられている実空間拘束条件は，求める

位相回復アルゴリズム

および ρ (x) が実数とする．GS における実空間拘束条件

3.1 GS アルゴリズム
位相回復とは，先に述べたように実空間拘束条件と逆空間拘束条件を満たす未知関数を探すことである．Gerchbergと Saxton5) は，実像とフーリエ像の強度が与えられている
ときに，フーリエ位相を求めるために，図 3 に表わされる
ように，拘束条件を適用しつつ，フーリエ変換と逆フーリ
エ変換を繰り返す手法を提案した．
実像とフーリエ像の振幅をそれぞれ ρnorm と Fobs として，実空間および逆空間拘束条件とする．まず，任意の実計測と制御

第 50 巻 第 5 号

実像が実数であり，上式における任意の x について ρ(x)2011 年 5 月号

について，誤差を下げるものとして，実空間拘束条件における n 番目の仮説実像から n + 1 番目の実像への以下の更新方法を Error reduction (ER) と名付けた．

ρn+1 (x) =ρn (x)
0

x∈
/D
,
x∈D

(5)

ただし D は実空間拘束条件を満たさない領域であり，逆空間から得られる ρ に起因して定まる領域である．

333

Fienup はさらに，D における更新方法として，以下のHybrid Input Output (HIO) という位相回復アルゴリズ
ムを提案した．これは位相回復においては，欠かせない中心的アルゴリズムとなっている6) ．

ρn+1 (x) =ρn (x)ρn (x) − βρn (x)ρ2 (x) =ρ(x)0

if ρ(x) < δ.
otherwise(9)

そして，δ を基準に ρ1 + ρ2 と符号を反転させた ρ1 − ρ2を用いる．実像 ρn から ρn+1 への更新として表現すると

x∈
/D
,
x∈D

(6)

ただし β は正の定数であり，効果のある設定値が知られてつぎのようになる．

ρn+1 (x) =ρn (x)−ρn (x)if ρn (x) ≥ δ.
if ρn (x) < δ(10)いる．HIO は，ρ に対し得られている ρ からの更新方法の

符号を反転させることから，このアルゴリズムは Charge

模索の中で示されたものであるが，ER における最急降下

ﬂipping と名付けられている8) ．このアルゴリズムにおけ法のような整然とした根底は，その後のいろいろな研究展

るしきい値の設定については，試行錯誤的な要素を排除す

開もあるにしても，その良さの真意は不明確のままである．

ることはできないが，その有効性から，結晶構造解析にお位相回復の良さは，仮説として得られる実像のフーリエ強度 |Fcal |2 が，観測値である |Fobs |2 とどのくらい一致す
るかにある．その指標として，以下の R-ファクターがよく用いられる．

k∈K

|Fobs (k)|3.4 画像復元との関連失ったフーリエ位相の回復には，与えられる限定的な拘束条件を利用する．実空間における実像関数の更新は，拘|Fcal (k)| − |Fobs (k)|R(Fcal , Fobs ) =いて多く使われている．

束条件を満たす集合へのある種の射影と言える．Fienup が.

(7)

HIO を世に送り出したと同じ時期に，Youla は凸集合上への射影（POCS: Projection onto Convex Sets）による画像修復法を展開した9) ．ヒルベルト空間上で，拘束条件となk∈K

GS アルゴリズムを元にした ER と HIO は，数値シミュ
レーションにおいてはもちろんのこと，1999 年の Miao から始まった実験で得られる回折像からの位相回復においても，その有用さが発揮されている．X 線や電子線による実験と位相回復による回折イメージングについては，本特集の記事で先端研究がご覧いただける．

3.3 HIO 以降の理論とアルゴリズム
HIO や ER による有限回の繰り返し写像計算によって得られる位相回復像について，収束や一意性に関する理論研究の中で，位相回復アルゴリズムを位置付けていくことは，
回折イメージングの進展に重要である．Elser は，差分写像

る複数の凸集合 Cj (j = 1, · · · , m) が与えられている．それらの積集合C0 = ∩m
j=1 Cj ,(11)

が空集合ではないとして，拘束条件をすべて満たす C0 に含まれる解を探索する．POCS とは，凸集合 Cj への射影

Pj を，j = 1, 2, · · · , m で繰り返していく方法である．図 4に凸集合 A;B に対する POCS を示す．射影を行うことは位相回復においても同様であるが，拘束条件を満たす集合の凸性が必ずしも保証されないために，POCS の理論を位相回復にそのまま適用はできない．

(Diﬀerential map) の枠組の中で，Fienup の ER や HIOなどを位置づけている7) ．
物体領域を完全に与えることは，対象の小ささから実験では難しいが，より良い実空間拘束条件を与えるために，物体サポート領域をタイトにすることは，位相回復の計算において望ましい．未知なサポートへの対応する 1 つの方法

図4

凸集合 A, B に対する POCS

が，Ozlányi と Süto によって示された．それは，GS アル
ゴリズムにおける実空間拘束条件の適用に特徴があり，仮トモグラフィは医療など広く応用されている．復元の根説の実像にしきい値を設定し，その値によって更新式を反本となる投影切断面定理とは，射影された像のフーリエ変転させる方式である．簡単に述べると，逆空間拘束条件か換は，実像のフーリエ変換の射影角度方向に相当することら得られた関数の逆フーリエ変換について，実数条件を適である．Gerchberg-Papoulis アルゴリズム10), 11) は，正逆用した ρ について，しきい値 δ によってつぎの 2 つの関のフーリエ変換によるアルゴリズムであり，GS アルゴリ

数 ρ1 と ρ2 に分ける．

ズムと構造的共通性がある．つまり，トモグラフィなどのρ1 (x) =

334

ρ(x)
0

if ρ(x) ≥ δ,otherwise

画像復元と位相回復とは，与えられる拘束条件の違いはあ(8)るが，GS ダイヤグラムを共通な土台として両者を見ることができる．

計測と制御

第 50 巻

第5号

2011 年 5 月号

4.

それらの間の情報量が共に小さいことに相当する．よって，情報量に基づく位相回復

GS ダイヤグラムに情報量を導入することで，図 5 として

4.1 実空間と逆空間の情報量表現できる．

位相回復は拘束条件の下での探索問題であり，関数の推定問題でもある．失ったフーリエ位相を探すためには，計算としてアルゴリズムが必要となるが，アルゴリズムとは，事前の仮説からつぎの事後の仮説を与える方式と考えると，2 項間の関係を導入することが求められる．たとえば確率分布の場合の 2 項間の関係として情報量があり，情報理論における主役である．そこで，つぎに情報量を出発点とすることで，位相回復のアルゴリズムに言及する12) ．
位相回復において，対象の構造を表わす物体関数 ρ とそのフーリエ変換 F が対に関係する．ρ と F は実空間 X

図 5 情報量と GS ダイヤグラム

および逆空間 K 上で複素数値をとる関数である（簡単のため，X, K は離散集合とする）．X 上の体積有限な関数全体の集合を P とする．
よく用いられる実空間拘束条件として，実数と正値性を仮定する．その条件を満たす P の要素の集合を P̄ とおく．
P̄ の 2 つの関数 ρ と τ 間の差異を測るために以下の量を導入する．

Dγ (ρ, τ ) =

1
x∈X

−

γ

ρ をどのように選ぶかが，位相回復アルゴリズムの定式化につながる．ρ のフーリエ変換 F と逆空間拘束 Fobs を適用した F  との間の情報量が小さいことが，望ましい ρ の条件となる．そこで，実空間と逆空間の情報量の和として以下を導入する．

ρ(x)(ρ(x)γ − τ (x)γL(λ) = Dγ (ρ, τ ) + λ E(F, F  ).

1
(ρ(x)1+γ − τ (x)1+γ ) ,
1+γ

(12)

ただし γ ∈ [0, 1] とする．この量は，Density-power di-vergence と呼ばれる情報量であり，ρ と τ が確率分布ときに，つまりx ρ(x) = 1 およびx τ (x) = 1 のときに，分布間距離の最小化の観点から，この量の統計的性質
が Basu によって示されている13) ．さらには，情報理論における情報量の解析的視点から，Csiszár の f -divergenceの拡張形式に基づいて位置付けられている14), 15) ．

τ と ρ をそれぞれ事前・事後の物体関数として，ρ のフーリエ変換 F と，F の位相をそのまま用いて逆空間拘束条
件 Fobs を振幅として入れ替えた F  との差異については，
以下の 2 乗ノルムがよく用いられる．
E(F, F  ) =|F (k) − F  (k)|2 .
(13)
k∈K

F と F  の位相は等しいことから，E(F, F  ) = E(|F |, |F  |)となる．つまり，比較する関数については，非負で体積有

限な実関数であるので，情報量 Dγ に F と F を導入し，
Dγ (F, F ) を定式化することができる．それと (13) 式とは以下の関係にある．

lim Dγ (F, F  ) = E(F, F  ).

γ→1

第 50 巻 第 5 号

2011 年 5 月号

(15)

ρ に関する上式の 2 項目の導関数を考慮するために，逆空間のほうは D1 = E としている．∂L(λ)/∂ρ = 0 から以下が得られる．

τ (x)γ − 1
ρ(x)γ − 1=
+ λ(ρ (x) − ρ(x)),
γ
γ

(16)

ただし，∂E(F, F  )/∂ρ については，Fienup の (4) 式を用いている．さらに，|ρ − τ | が十分小と仮定し，その近似とτ と ρ の関係に基づいて，cγ を正の定数とする ρn から
ρn+1 への更新を，以下のように定めることができる．
1

ρn+1 (x) = {ρn (x)γ + cγ (ρn (x) − ρn (x))} γ . (17)これは，逆空間拘束条件を満たさない方向としたときに得られる ρn から ρn+1 への更新を元に，時定数 n を逆転させることで解くことができる．

(17) 式において，γ = 1 のときには，
ρn+1 (x) = (1 − cγ )ρn (x) + cγ ρn (x),(18)

となり，さらに cγ = 1 とすることで，位相回復における

ER アルゴリズムが得られる．γ → 0 のときには，
ρn+1 (x) = ρn (x) exp[cγ (ρn (x) − ρn (x))],(19)が得られる．このとき，(15) 式における第 1 項は，(14)

4.2 情報量最小化と位相回復
GS アルゴリズム（図 3）においては，ρ と ρ および Fと F  とそれぞれの差異が小さいことが望ましい．つまり，計測と制御事前の物体関数 τ が与えられたときに，事後の物体関数

lim Dγ (ρ, τ ) =γ→0ρ(x) ln

x∈X

−




ρ(x)
+
τ (x)
τ (x)
x∈X

ρ(x),

(20)

x∈X

335

となる．この情報量は I-divergence16) と呼ばれ，τ と ρ をただし，F0 はユニットセルの総電子数とする．Fobs (k),確率分布としたときには，情報理論でよく登場するカルバッ

σ(k)2 , 係数 λ1 および λ2 が与えられると電子密度関数がクの情報量17) である．

得られる．

4.3 MEM 結晶構造解析と位相回復
最大エントロピー法 (Maximum Entropy Method :MEM) は，画像復元などでおなじみな方法であるが，その結晶構造解析と位相回復との関係について，1990 年に有用性の一端として，結晶の原子配置を見出す結晶構造解

Millane は，結晶構造解析と位相回復は互いに別々に研究が進展したことを指摘している21) ．その関係について，Collinsの MEM 結晶構造解析と情報量に基づいた位相回復との関析においても用いられている．Collins18), 19) による MEM

わりから，つぎに述べる22) ．

結晶構造解析法の展開に従って述べる．実空間と逆空間のまず，正規化された関数を用いているエントロピー S

インデックスをそれぞれ r と k として，まず，Jaynes

に対して，前節において紹介した正規化を不要とする I-

20)

divergence を導入する．

のエントロピー関数を導入する．

S(ρ̄, τ̄ ) = −ρ̄(r) ln

r

ρ̄(r)
,
τ̄ (r)

ρ(r)τ (r)ρ̄(r) = , τ̄ (r) = 
,
)

ρ(r

r
r  τ (r )よび τ̄ (r) は正規化された分布関数である．既知位相と未知の位相に関する拘束条件を以下のように与える．

位相回復においては，すべてが未知の位相であることから，未知の位相に関する構造因子を用いた逆空間拘束条件のみを用いる．GS ダイヤグラムに従い，失われている Fobs の位相として，Fcal の位相を適用することから，(4) 式から以下が得られる．

∂ |Fcal (k) − Fobs (k)|2 = 2(ρ(r) − ρ (r)),
∂ρ(r)

(23)

k

ただし，M1 および M2 は反射数を表わし，以下の Fcal は計算された構造因子を表わす．

Fcal (k) = V

exp(−i2πr · k),

(29)

(24)

k




ρ(r) 
+τ (r) −
ρ(r).
τ (r)
r
r
(28)

k

1  ||Fcal (k)| − |Fobs (k)||2
C2 =
,
M2
σ(k)2

ρ(r) ln(22)

布関数である．それらは電子密度関数に対応する．ρ̄(r) お

1  |Fcal (k) − Fobs (k)|2
,
M1
σ(k)2


r

ただし，ρ(r) は推定分布関数，τ (r) は事前仮説となる分C1 =I(ρ, τ ) =(21)ρ (r) は，|Fobs (k)| exp[iψ(k)] の逆フーリエ変換，ψ(k)は Fcal の位相とする．(29) 式を利用することで以下が得
られる．

(25)

ρ(r) = exp[ln τ (r) + C(ρ (r) − ρ(r))],

r

(30)

V はユニットセルの体積，Fobs は観測された構造因子，
σ(k)2 は Fobs の統計的な分散を表わす．構造因子の拘ただし C は正定数とする．これは，|ρ − τ | が十分小を仮定

束条件のもとでのエントロピーの最大化のために以下のに帰着される．つまり，位相回復と MEM 結晶構造解析は，

L(λ1 , λ2 ) を用いる．

一般化された情報量を通じて関係付けられることがわかる．

L(λ1 , λ2 ) = S(ρ, τ ) −λ1
λ2
C1 − C2 .
2
2

(26)

最大エントロピーとなる確率分布は，つぎのように求めら

5.

まとめ

失ったフーリエ位相を探すことで，直接観測が難しい対象の構造がわかることから，光源に対し広く応用できる位れる18) ．


ρ(r) = exp ln τ (r)+

相回復は有用な方法である．GS アルゴリズムと Fienup の

λ1 F0  Fobs (k) − Fcal (k)
M1
σ(k)2
k

× exp{−i2πk · r}
|Fobs (k)| exp[iψ(k)]−Fcal(k)

λ2 F0
+
M2
σ(k)2
k

× exp{−i2πk · r} ,
(27)

336

することに留意して，情報量 Dγ における γ → 0 の場合

HIO 等によって，有限回の繰り返し計算で得られる仮説実像が真か，あるいは真に十分近いかは，回折イメージング
における大きな問題である．80 年代初頭に Fienup によって位相回復アルゴリズムが展開されたと同時期に，国内においても位相回復アルゴリズムに関する研究は，精力的に進められた23)〜26) ．その後，位相回復における解の一意性など，根本の理論的問題が依然として残っている中で，本内容の後半において，位相回復を情報量最小化と関連させ，情報量を軸とした位相回復アルゴリズムの表現について述

計測と制御

第 50 巻

第5号

2011 年 5 月号

べた．これは，80 年代の位相回復のアルゴリズムに関連した研究とは異なる方向性を与えているものと著者は考えている．それに関連する研究として，与えられる Fobs がフーリエ変換の関係から量子ノイズによって離れている場合には，位相回復解の集合が関数空間上で球殻構造となること
や，繰り返し写像による統計的力学系の重要性を筆者は言及している27) ．数理を主体とした境界領域からのアプロー
チで，この問題に新たなブレークスルーが期待されている．
アルゴリズムや計算においては，境界領域からのこれまでとは違ったアプローチも含めて，位相回復解に関する基礎理論や実験データに対応できるアルゴリズムの開発などが，今後の重要な課題となる．
（2011 年 1 月 24 日受付）
参

考

文 献

1） Q. シェン，Q. ハオ，S. M. グリュナー：X 線回折像から，複雑な物質の構造を復元するために〜失われた位相を求めて〜(奥山健二訳)，パリティ，10 (2006)
2） J. Miao, P. Charalambous, J. Kirz and D. Sayre: Extending
the methodology of X-ray crystallography to allow imaging of micrometer-sized non-crystalline specimens, Nature,400, 342/344 (1999)
3） D. Sayre: Some implications of a theorem due to Shannon,
Acta Crystallogr., 5, p.843 (1952)
4） C.E. Shannon: A Mathematical Theory of Communication,
Bell Sys. Tech. Journal, 27, 456/477 (1948)
5） R.W. Gerchberg and W.O. Saxton: A practical algorithm
for the determination of phase from image and diﬀraction
plane pictures, Optic, 35, 237/246 (1972)
6） J.R. Fienup: Phase retrieval algorithms: a comparison, Applied Optics, 21, 2758/2769 (1982)
7） V. Elser: Phase retrieval by iterated projections, J. Opt.
Soc. Am. A, 20–1, 40/55 (2003)
8） G. Oszlányi and A. Süto: AB initio structure by charge
ﬂipping, Acta Cryst., A60, 134/141 (2003)
9） D.C. Youla and H. Webb: Image Restriction by the Methodof Convex Projections: Part 1 – Theory, IEEE Trans. Med.
Imaging., MI-1, 81/95 (1982)
10） R.W. Gerchberg: Super-resolution through Error Energy
Resolution, Opt. Acta., 21–9, 709/720 (1974)
11） A. Papoulis: A New Algorithm in Spectral Analysis andBand-limited Extrapolation, IEEE Trans. Circuits & Syst.,
CAS-22–9, 735/742 (1975)
12） H. Shioya and K. Gohara: Generalized phase retrieval algorithm based on information measures, Optics Communications, 266, 88/93 (2006)
13） A. Basu, I.R. Harris, N. Hjort and M.C. Jones: Robust andeﬃcient estimation by minimizing a density power divergence, Biometrika, 85, 549/559 (1998)14） 内田真人，塩谷浩之：不定性を用いた分布間情報量の拡張形式に関する検討，信学論，J87-A–4, 546/553 (2004) [英訳：Electronics and Communication in Japan: Part 3, Fundamental

計測と制御

第 50 巻 第 5 号

2011 年 5 月号

Electronic Science, 88–4, 35/42 (2005)]
15） I. Csiszár: Information-type measures of diﬀerence of probability distributions and indirect observation, Studia Sci.
Math. Hunger., 2, 299/318 (1967)16） I. Csiszár: Why least squares and maximum entropy? an axiomatic approach to inverse problems, Ann. Stat., 19,2033/2066 (1991)17） S. Kullback and R.A. Leibler: On Information and Suﬃciency, Ann. Math. Stat., 22, 79/86 (1951)
18） D.M. Collins: Electron density images from imperfect databy iterative entropy maximization, Nature (London), 298,
49 (1982)
19） M. Sakata and M. Sato: Accurate structure analysis by
the maximum-entropy method, Acta Cryst., A46, 263/270(1990)
20） E.T. Jaynes: Prior Probabilities, IEEE Trans. Syst. Cybern. SSC4, 227/241 (1968)
21） R.P. Millane: Phase retrieval in crystallography and optics,
J. Opt. Soc. Am. A, 7–3, 394/411 (1990)22） H. Shioya and K. Gohara: Maximum entropy method for
diﬀractive imaging, J. Opt. Soc. Am. A, 25–11, 2846/2850(2008)
23） J. Maeda and K. Murata: Retrieval of wave aberration frompoint spread function of optical transfer function data, Applied Optics, 20, 274/279 (1981)24） 前田純治，村田和美：光学における位相回復問題，光学，11–3,
230/240 (1982)
25） 中島伸治，朝倉利光：ヒルベルト変換による位相回復，光学，
11–3, 241/252 (1982)
26） H. Takajo and T. Takahashi: Least-squares phase estimation from the phase diﬀerence, J. Opt. Soc. Am. A, 5–3,
416/425 (1988)
27） H. Shioya and K. Gohara: Spherical shell structure of distribution of images reconstructed by diﬀractive imaging, J.
Opt. Soc. Am. A, 27–5, 1214/1218 (2010)［著者 紹 介］

しお

や

ひろ

ゆき

塩

谷

浩

之 君（正会員）
1990 年北海道大学理学部数学科卒業，北海道大学大学院工学研究科情報工学専攻修士課程修了，同大学院博士後期課程修了，博士 (工学)．北海道
大学工学研究科システム情報工学専攻助手，室蘭工業大学工学部助教授を経て，現在，同大学教授．
データ学習，ニューラルネット，強化学習，位相回復などに興味をもつ．電子情報通信学会，ACM,
OSA などの会員．

ごう

はら

かず

とし

郷

原

一

寿 君（正会員）（本号 p.331 参照）

337