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

提供: Eospedia
移動: 案内検索
行71: 行71:
  
 
===位置合わせ===
 
===位置合わせ===
====[[mrcImageAutoRotationCorrelation]]を使用した場合====
+
====[[mrcImageCorrelation]]を使用した場合====
<div> [[Media:Makefile-Tomogram.zip|こちらのMakefile]]を使用して位置合わせを実行します。Makefileを使用せずにコマンドを逐次入力しても位置合わせは可能ですが、-Oの出力ファイルを[[mrc2Dto3D]]用のフォーマットに変換する必要があります。</div>
+
<div> [[Media:Makefile-Tomogram.zip|こちらのMakefile]]を使用して位置合わせを実行します。Makefileを使用せずにコマンドを逐次入力しても位置合わせは可能ですが、[[mrc2Dto3D]]用の入力ファイルリストを作成する必要があります。</div>
 
<br>
 
<br>
  
 
=====実行例1=====
 
=====実行例1=====
<div>[[Media:1WDC-win-2dSet.zip|入力ファイル]](全ファイル)<br>
+
<div>[[Media:1WDC-Tom-Mask-2dSet.zip|入力ファイル]](全ファイル)<br>
 
<table>  
 
<table>  
 
<tr>  
 
<tr>  
<td><p align="Center">[[画像:1WDC-win-2dSet.png]]<br>
+
<td><p align="Center">[[画像:1WDC-Tom-mask-2dSet.png]]<br>
[[mrcImageMake2DArrayImage]](N=12)で表示<br></p>  
+
10°刻みで表示<br></p>  
 
</td>  
 
</td>  
 
</tr>  
 
</tr>  
 
</table>  
 
</table>  
 
<br>
 
<br>
 
<div>[[Media:1WDC-shift1.ref2d|参照ファイル]]の画像</div>
 
<table>
 
<tr>
 
<td><p align="Center">[[画像:1WDC-shift1-2d.png]]<br>
 
[[mrcImageMake2DArrayImage]](N=12)で表示<br></p>
 
</td>
 
</tr>
 
</table>
 
<br>
 
 
<div>設定内容</div>
 
<pre>
 
### For mrcImageAutoRotationCorrelation
 
# mode
 
COR_MODE=18
 
# StackOrder
 
STC_MODE=1
 
</pre>
 
  
 
<div>コマンド</div>
 
<div>コマンド</div>
 
<pre>
 
<pre>
make RotFit
+
make CorFit1
</pre>
+
<br>
+
 
+
<div>実行結果</div>
+
<div>[[Media:1WDC-corfit.corinfo.zip|角度情報]]</div>
+
<pre>
+
Set1002.fit YOYS 60.000002 0.000000 0.000000 0.097586
+
Set1000.fit YOYS 0.000000 0.000000 0.000000 0.103920
+
Set1001.fit YOYS 30.000001 0.000000 0.000000 0.094755
+
Set1005.fit YOYS 149.999997 0.000000 0.000000 0.120492
+
Set1004.fit YOYS 120.000003 0.000000 0.000000 0.115170
+
Set1003.fit YOYS 90.000003 0.000000 0.000000 0.096009
+
Set1006.fit YOYS 180.000005 0.000000 0.000000 0.103048
+
Set1008.fit YOYS 240.000007 0.000000 0.000000 0.096656
+
Set1007.fit YOYS 209.999999 0.000000 0.000000 0.099589
+
Set1009.fit YOYS 270.000001 0.000000 0.000000 0.101168
+
Set1010.fit YOYS 299.999995 0.000000 0.000000 0.115695
+
Set1011.fit YOYS 329.999989 0.000000 0.000000 0.122864
+
 
</pre>
 
</pre>
<div>今回は入力ファイルと参照ファイルがほぼ同じであるので、角度情報は一致します。</div>
 
 
<br>
 
<br>
  
<div>[[Media:1WDC-corfit-2DSet.zip|出力ファイル]]</div>
+
<div>[[Media:1WDC-Tom-Fit-2DSet.zip|出力ファイル]]</div>
 
<table>  
 
<table>  
 
<tr>  
 
<tr>  
<td><p align="Center">[[画像:1WDC-corfit-2d.png]]<br>
+
<td><p align="Center">[[画像:1WDC-Tom-Fit-2d.png]]<br>
[[mrcImageMake2DArrayImage]](N=12)で表示<br></p>  
+
10°刻みで表示<br></p>  
 
</td>  
 
</td>  
 
</tr>  
 
</tr>  
行152: 行114:
 
<br>
 
<br>
  
<div>[[Media:1WDC-win-2dSet.zip|入力ファイル]]<br>
+
<div>[[Media:1WDC-Tom-Fit-2dSet.zip|入力ファイル]]<br>
 
<table>  
 
<table>  
 
<tr>  
 
<tr>  
<td><p align="Center">[[画像:1WDC-win-2dSet.png]]<br>
+
<td><p align="Center">[[画像:1WDC-Tom-Fit-2dSet.png]]<br>
[[mrcImageMake2DArrayImage]](N=12)で表示<br></p>  
+
10°刻みで表示<br></p>  
 
</td>  
 
</td>  
 
</tr>  
 
</tr>  
行163: 行125:
  
 
<div>コマンド: [[mrc2Dto3D]] -i Input.3dlst -o Input.3d -m 1</div>
 
<div>コマンド: [[mrc2Dto3D]] -i Input.3dlst -o Input.3d -m 1</div>
<br>
 
 
<div>Input.3dlstのファイル</div>
 
<pre>
 
Set1000.fit YOYS 0 0 0
 
Set1001.fit YOYS 30 0 0
 
Set1002.fit YOYS 60 0 0
 
Set1003.fit YOYS 90 0 0
 
Set1004.fit YOYS 120 0 0
 
Set1005.fit YOYS 150 0 0
 
Set1006.fit YOYS 180 0 0
 
Set1007.fit YOYS 210 0 0
 
Set1008.fit YOYS 240 0 0
 
Set1009.fit YOYS 270 0 0
 
Set1010.fit YOYS 300 0 0
 
Set1011.fit YOYS 330 0 0
 
</pre>
 
 
<br>
 
<br>
  
行193: 行138:
 
<td>[[画像:Outdata1-mrc2Dto3D.png]]</td>  
 
<td>[[画像:Outdata1-mrc2Dto3D.png]]</td>  
 
<td>[[画像:Outdata2-mrc2Dto3D.png]]</td>  
 
<td>[[画像:Outdata2-mrc2Dto3D.png]]</td>  
</tr>
 
<tr>
 
<td>[[画像:Outdata3-mrc2Dto3D.png]]</td>
 
<td>[[画像:Outdata4-mrc2Dto3D.png]]</td>
 
<td>[[mrcImageLowPassFilter]]を掛けた場合</td>
 
 
</tr>  
 
</tr>  
 
</table>  
 
</table>  
行208: 行148:
 
標準誤差<br></p>  
 
標準誤差<br></p>  
 
</td>  
 
</td>  
<td><p align="left">-0.00811993 (43, 34, 36)<br>  
+
<td><p align="left">-0.00437076 (39, 34, 36)<br>  
0.0128274 (23, 25, 33)<br>  
+
0.00799233 (37, 34, 36)<br>  
1.79873e-06<br>  
+
3.45089e-06<br>  
0.00120493<br>  
+
0.000342836<br>  
2.35338e-06<br></p>  
+
6.69601e-07<br></p>  
 
</td>  
 
</td>  
 
</tr>  
 
</tr>  

2014年4月22日 (火) 08:59時点における版

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

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

連続傾斜像の撮影

 同一視野で撮影すると試料を特定の軸に対して一定間隔で回転した画像になります。


1WDC-Tom-2dSet.png

1軸回転

Input1-Tomogram.png

2軸回転


画像の補正

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

ラフ・アラインメント

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

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

1WDC-Tom-2dSet.png
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


出力ファイル(全ファイル)

1WDC-Tom-mask-2dSet.png
10°刻みで表示


位置合わせ

mrcImageCorrelationを使用した場合

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


実行例1
入力ファイル(全ファイル)

1WDC-Tom-mask-2dSet.png
10°刻みで表示


コマンド
make CorFit1


1WDC-Tom-Fit-2d.png
10°刻みで表示


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

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

3次元再構成

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

mrc2Dto3Dで実行

実行例1

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


入力ファイル

ファイル:1WDC-Tom-Fit-2dSet.png
10°刻みで表示


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


出力ファイル
xy平面 yz平面 zx平面
Outdata-mrc2Dto3D.png Outdata1-mrc2Dto3D.png Outdata2-mrc2Dto3D.png

最小

最大
平均値
標準偏差

標準誤差

-0.00437076 (39, 34, 36)

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

6.69601e-07


実行例2(2軸回転)

2軸で回転した画像にて3次元再構成を行います。
入力ファイル(全ファイル)

1WDC-double-2dSet.png
mrcImageMake2DArrayImage(N=12)で表示


1WDC-shift1-3d.png
xy平面

1WDC-shift1-3d2.png
yz平面

最小

最大
平均値
標準偏差

標準誤差

0 (0, 0, 0)

4 (31, 26, 26)
0.0116844
0.147033

0.000287174


今回はこちらのMakefileを使って3次元像を作成します。



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

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

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

平滑化

セグメンテーション

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

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

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