「電子線トモグラフィー」の版間の差分

2014年4月24日 (木) 01:06時点における版

電子線トモグラフィーとは、電子顕微鏡を用いた３次元再構成法の一つで、同一視野を様々な方向からの投影された電子顕微鏡像をコンピュータの中で三次元像に再構成し、コンピュータを使って断層像（トモグラム）を作成する手法です。

　下記にその手順を順に述べていきます。

1 連続傾斜像の撮影
2 画像の補正
3 ラフ・アラインメント
- 3.1 前処理(窓関数)
- 3.2 位置合わせ
  - 3.2.1 mrcImageCorrelationを使用した場合
    - 3.2.1.1 実行例1
4 ファイン・アラインメント
5 ３次元再構成
- 5.1 mrc2Dto3Dで実行
  - 5.1.1 実行例1
  - 5.1.2 実行例2(２軸回転)
- 5.2 ラドン変換を使った３次元再構成で実行
6 電子線トモグラフィー画像のもつ問題点
- 6.1 ミッシングエリア
7 トモグラフに使われる画像処理
- 7.1 平滑化
- 7.2 セグメンテーション
8 トモグラム画像の解釈 ==

連続傾斜像の撮影

　この手法で撮影すると試料を特定の軸に対して一定間隔で回転した画像になります。また、回転角度が大きくなるにつれて視野は広くなるので、映っている物体は回転軸の中心に集まってきます。

１軸回転

２軸回転

画像の補正

　傾斜画像では、視野全体がアンダーフォーカスとなるように、通常の電子顕微鏡撮影よりも大きなデフォカース値をとる場合が多いです。その場合、LaB6などの電子銃では第一トーンリングより外側の情報がほとんど無いために、CTF補正等は必要ありませんが、電界放出銃を利用した場合には、場合によっては間違った情報を与える場合があるので注意が必要です。

ラフ・アラインメント

　傾斜画像の中心付近の画像の相関から、それぞれの傾斜画像の大まかな位置を合わせます。

入力ファイル(Y軸周りに±60°の範囲で2°刻みで回転)

10°刻みで表示

前処理(窓関数)

　窓関数を掛けて画像端のノイズを除去します。

今回は横80%、縦90%の窓関数を掛けます。

コマンド(1枚目の場合): mrcImageWindowing -i Set1000.roi -o Set1000.mask -W 0.2 0.2 0.1 0.1

こちらのMakefileを使用すると全てのroiファイルについて処理を行います。

設定内容

# For Windowing
WIN_X=0.2
WIN_X_MAX=0.2
WIN_Y=0.1
WIN_Y_MAX=0.1
WIN_MODE=18

コマンド

make Windowing

出力ファイル

10°刻みで表示

位置合わせ

　回転角度に応じて対象が平行移動しているので、位置を合わせます。このとき回転による影響を受けにくくするために隣同士の角度で位置合わせを行います。

例. 2°刻みで撮影した画像の場合
　2°画像の位置を0°画像の位置に合わせる
　4°画像の位置を位置補正後の2°画像の位置に合わせる
　6°画像の位置を位置補正後の4°画像の位置に合わせる
　.
　.

　(+が終わったら-についても同様に位置を合わせる)

mrcImageCorrelationを使用した場合

　こちらのMakefileを使用して位置合わせを実行します。Makefileを使用せずにコマンドを逐次入力しても位置合わせは可能ですが、mrc2Dto3D用の入力ファイルリストを作成する必要があります。

実行例1

入力ファイル

10°刻みで表示

コマンド

make CorFit1

出力ファイル

10°刻みで表示

ファイン・アラインメント

　傾斜画像の軸の位置、角度を出来る限り一致させます。

３次元再構成

　２次元画像のセットから３次元画像を再構成します。

元の3Dデータの画像

xy平面

yz平面

zx平面

最小

最大
平均値
標準偏差

標準誤差

0 (0, 0, 0)

4 (31, 26, 26)
0.0116844
0.147033

0.000287174

mrc2Dto3Dで実行

実行例1

　連続傾斜像の撮影で軸にズレがない場合は角度刻みをそのまま入力ファイルリストに設定します。

入力ファイル

10°刻みで表示

コマンド: mrc2Dto3D -i Input.3dlst -o Input.3d -m 1

こちらのMakefileを使用すれば、以下のコマンドでも実行できます。

make 3DList
make Input.3d

出力ファイル

xy平面

yz平面

zx平面

最小

最大
平均値
標準偏差

標準誤差

-0.00437076 (39, 34, 36)

0.00799233 (37, 34, 36)
3.45089e-06
0.000342836

6.69601e-07

実行例2(２軸回転)

２軸で回転した画像にて３次元再構成を行います。

入力ファイル(X, Y軸周りに±60°の範囲でそれぞれ10°刻みで回転)

X軸は30°、Y軸は10°刻みで表示

こちらのMakefileを使って、補正およびアラインメント、３次元再構成を行います。

コマンド

make All

出力ファイル

xy平面

yz平面

zx平面

最小

最大
平均値
標準偏差

標準誤差

-0.0012516 (30, 25, 51)

0.0027149 (31, 31, 38)
4.84671e-06
0.000196091

3.82989e-07

ラドン変換を使った３次元再構成で実行

　mrc2Dto3Dを使用するほかにラドン変換から３次元像を作る方法があります。ラドン変換は以下の流れで行います。

位置合わせ済みの2Dリスト
↓mrcImageSinogramCreate
シノグラムリスト
↓mrcRadon2Dto3D
3Dラドンファイル
↓mrcImageInverseRadonTransform

3Dファイル(完了)

こちらのMakefileを使い、mrc2Dto3Dの実行例2と同じ入力ファイルで実行します。

コマンド

make All

出力ファイル

xy平面	yz平面	zx平面

最小

最大
平均値
標準偏差

標準誤差

393.858 (36, 9, 39)

1427.57 (39, 39, 40)
612.613
138.895

0.194111

電子線トモグラフィー画像のもつ問題点

ミッシングエリア

　全方位からの投影像を得られないために、ミッシングエリアとよばれる情報がない領域があります。そのため、方向依存でのボケを生じます。

元の３次元像

xy平面

yz平面

zx平面

最小

最大
平均値
標準偏差

標準誤差

0 (0, 0, 0)

2 (29, 29, 14)
0.108302
0.364429

0.000711774

１軸回転で投影して３次元再構成した場合

xy平面

yz平面

zx平面

最小

最大
平均値
標準偏差

標準誤差

-0.00231156 (14, 23, 26)

0.00260412 (24, 21, 37)
1.72166e-05
0.000406119

7.93202e-07

Y軸: -60° ~ 60°: 10°刻み

２軸回転で投影して３次元再構成した場合

xy平面

yz平面

zx平面

最小

最大
平均値
標準偏差

標準誤差

-0.00164345 (19, 37, 18)

0.00201538 (48, 23, 32)
2.15537e-05
0.000351101

6.85745e-07

X, Y軸: 各-60° ~ 60°: 10°刻み

トモグラフに使われる画像処理

平滑化

セグメンテーション

トモグラム画像の解釈 ==

　また、多様な構造を含むため、３次元画像から興味あるセグメントを切り出すなどの作業が必要です。そのための支援ソフトウェアが必要です。

== サブトモグラムの平均化==

@@ 行301: / 行301: @@
 == [[電子線トモグラフィー画像のもつ問題点]] ==
-　全方位からの投影像を得られないために、ミッシングエリアとよばれる情報がない領域があります。そのため、方向依存でのボケを生じます。
+===ミッシングエリア===
+<div>　全方位からの投影像を得られないために、ミッシングエリアとよばれる情報がない領域があります。そのため、方向依存でのボケを生じます。</div>
+<br>
+<div>[[Media:Input-ellipsoidal.3d|元の３次元像]]</div>
+	<table>
+		<tr>
+			<td><p align="Center">[[画像:Input-ellipsoidal.png]]<br>
+				xy平面<br></p>
+			</td>
+			<td><p align="Center">[[画像:Input1-ellipsoidal.png]]<br>
+				yz平面<br></p>
+			</td>
+			<td><p align="Center">[[画像:Input2-ellipsoidal.png]]<br>
+				zx平面<br></p>
+			</td>
+			<td><p align="left">最小<br>
+				最大<br>
+				平均値<br>
+				標準偏差<br>
+				標準誤差<br></p>
+			</td>
+			<td><p align="left">0 (0, 0, 0)<br>
+(29, 29, 14)<br>
+.108302<br>
+.364429<br>
+.000711774<br></p>
+			</td>
+		</tr>
+	</table>
+<br>
+<div>１軸回転で投影して３次元再構成した場合</div>
+	<table>
+		<tr>
+			<td><p align="Center">[[画像:Output-ellipsoidal.png]]<br>
+				xy平面<br></p>
+			</td>
+			<td><p align="Center">[[画像:Output1-ellipsoidal.png]]<br>
+				yz平面<br></p>
+			</td>
+			<td><p align="Center">[[画像:Output2-ellipsoidal.png]]<br>
+				zx平面<br></p>
+			</td>
+			<td><p align="left">最小<br>
+				最大<br>
+				平均値<br>
+				標準偏差<br>
+				標準誤差<br></p>
+			</td>
+			<td><p align="left">-0.00231156 (14, 23, 26)<br>
+.00260412 (24, 21, 37)<br>
+.72166e-05<br>
+.000406119<br>
+.93202e-07<br></p>
+			</td>
+		</tr>
+	</table>
+<div>Y軸: -60° ~ 60°: 10°刻み</div>
+<br>
+<div>２軸回転で投影して３次元再構成した場合</div>
+	<table>
+		<tr>
+			<td><p align="Center">[[画像:Output3-ellipsoidal.png]]<br>
+				xy平面<br></p>
+			</td>
+			<td><p align="Center">[[画像:Output4-ellipsoidal.png]]<br>
+				yz平面<br></p>
+			</td>
+			<td><p align="Center">[[画像:Output5-ellipsoidal.png]]<br>
+				zx平面<br></p>
+			</td>
+			<td><p align="left">最小<br>
+				最大<br>
+				平均値<br>
+				標準偏差<br>
+				標準誤差<br></p>
+			</td>
+			<td><p align="left">-0.00164345 (19, 37, 18)<br>
+.00201538 (48, 23, 32)<br>
+.15537e-05<br>
+.000351101<br>
+.85745e-07<br></p>
+			</td>
+		</tr>
+	</table>
+<div>X, Y軸: 各-60° ~ 60°: 10°刻み</div>
+<br>
 == トモグラフに使われる画像処理 ==

「電子線トモグラフィー」の版間の差分

2014年4月24日 (木) 01:06時点における版

目次

連続傾斜像の撮影

画像の補正

ラフ・アラインメント

前処理(窓関数)

位置合わせ

mrcImageCorrelationを使用した場合

実行例1

ファイン・アラインメント

３次元再構成

mrc2Dto3Dで実行

実行例1

実行例2(２軸回転)

ラドン変換を使った３次元再構成で実行

電子線トモグラフィー画像のもつ問題点

ミッシングエリア

トモグラフに使われる画像処理

平滑化

セグメンテーション

トモグラム画像の解釈 ==

案内メニュー

表示

個人用ツール

案内

検索

ツール