「電子線トモグラフィー」の版間の差分
行71: | 行71: | ||
===位置合わせ=== | ===位置合わせ=== | ||
− | <div> | + | <div> 回転角度に応じて対象が平行移動しているので、位置を合わせます。このとき回転による影響を受けにくくするために隣同士の角度で位置合わせを行います。<br> |
+ | 例. 2°刻みで撮影した画像の場合<br> | ||
+ | 2°画像の位置を0°画像の位置に合わせる<br> | ||
+ | 4°画像の位置を'''位置補正後の'''2°画像の位置に合わせる<br> | ||
+ | 6°画像の位置を'''位置補正後の'''4°画像の位置に合わせる<br> | ||
+ | .<br> | ||
+ | .<br> | ||
+ | (+が終わったら-についても同様に位置を合わせる)</div> | ||
<br> | <br> | ||
行106: | 行113: | ||
== [[ファイン・アラインメント]] == | == [[ファイン・アラインメント]] == | ||
− | <div> | + | <div> 傾斜画像の軸の位置、角度を出来る限り一致させます。</div> |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
<br> | <br> | ||
2014年4月30日 (水) 00:00時点における版
電子線トモグラフィーとは、電子顕微鏡を用いた3次元再構成法の一つで、同一視野を様々な方向からの投影された電子顕微鏡像をコンピュータの中で三次元像に再構成し、コンピュータを使って断層像(トモグラム)を作成する手法です。
下記にその手順を順に述べていきます。
目次
連続傾斜像の撮影
1軸回転 2軸回転 |
画像の補正
傾斜画像では、視野全体がアンダーフォーカスとなるように、通常の電子顕微鏡撮影よりも大きなデフォカース値をとる場合が多いです。その場合、LaB6などの電子銃では第一トーンリングより外側の情報がほとんど無いために、CTF補正等は必要ありませんが、電界放出銃を利用した場合には、場合によっては間違った情報を与える場合があるので注意が必要です。
ラフ・アラインメント
傾斜画像の中心付近の画像の相関から、それぞれの傾斜画像の大まかな位置を合わせます。
前処理(窓関数)
# 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
位置合わせ
例. 2°刻みで撮影した画像の場合
2°画像の位置を0°画像の位置に合わせる
4°画像の位置を位置補正後の2°画像の位置に合わせる
6°画像の位置を位置補正後の4°画像の位置に合わせる
.
.
mrcImageCorrelationを使用した場合
実行例1
make CorFit1
ファイン・アラインメント
実行例1
# RotMode ROTMODE=ZOYS # Rot1 ROT1MIN=10 ROT1MAX=10 ROT1D=10 # Rot2 ROT2MIN=-60 ROT2MAX=60 ROT2D=2 # Rot3 ROT3MIN=0 ROT3MAX=0 ROT3D=30 ### For mrcImageMove SHIFT2MAX=10 SHIFT3MAX=0 ### For mrcImageTiltAxisSearch TILTMIN=0 TILTMAX=20 TILTD=2
make TiltFit
9.895
補正前 | ||
補正後 | ||
xy平面 | yz平面 | zx平面 |
実行例2(繰り返し)
# For mrcImageTiltAxisSearch IN_TILT_EXT=roi # RotMode ROTMODE=ZOYR # Rot1 ROT1MIN=10 ROT1MAX=10 ROT1D=10 # Rot2 ROT2MIN=-60 ROT2MAX=60 ROT2D=10 # Rot3 ROT3MIN=0 ROT3MAX=0 ROT3D=10 ### For mrcImageTiltAxisSearch TILTMIN=-10 TILTMAX=10 TILTN=10 TILTITER=100 TILTSCALE=5
make TiltFit
0.084 0.161 0.230 0.294 0.345 -中略- 0.701 0.700 0.685 0.706 0.701
3次元再構成
最小 最大 |
0 (0, 0, 0) 4 (31, 26, 26) |
mrc2Dto3Dで実行
実行例1
make 3DList make Input.3d
最小 最大 |
-0.00437076 (39, 34, 36) 0.00799233 (37, 34, 36) |
実行例2(2軸回転)
make Windowing make CorFit1 make CorFit2 make 3D
最小 最大 |
-0.0012516 (30, 25, 51) 0.0027149 (31, 31, 38) |
ラドン変換を使った3次元再構成で実行
位置合わせ済みの2Dリスト
↓mrcImageSinogramCreate
シノグラムリスト
↓mrcRadon2Dto3D
3Dラドンファイル
↓mrcImageInverseRadonTransform
make Windowing make CorFit1 make CorFit2 make Radon3D
xy平面 | yz平面 | zx平面 |
最小 最大 |
248.518 (3, 0, 15) 1509.49 (31, 24, 39) |
電子線トモグラフィー画像のもつ問題点
ミッシングエリア
最小 最大 |
0 (0, 0, 0) 2 (29, 29, 14) |
# Initial Reference Model INITIAL=Input-ellipsoidal # For mrcImageCorrelation IN_CORR_EXT=roi IN_CORR_LIST=ROIs ### For mrcImageMove SHIFT2MAX=0 SHIFT3MAX=0
make TestData make CorFit1 make CorFit2 make 3D
# Rot2 ROT2MIN=-60 ROT2MAX=60 ROT2D=10 # Rot3 ROT3MIN=0 ROT3MAX=0 ROT3D=10
最小 最大 |
-0.00196191 (61, 38, 32) 0.00269849 (24, 21, 37) |
# Rot2 ROT2MIN=-60 ROT2MAX=60 ROT2D=10 # Rot3 ROT3MIN=-60 ROT3MAX=60 ROT3D=10
最小 最大 |
-0.00143731 (50, 32, 45) 0.00230465 (53, 31, 37) |
画像の位置、フォーカスの問題
平行移動の補正
軸の決定
トモグラフに使われる画像処理
平滑化
実行例1
最小 最大 |
-0.00268221 (76, 92, 66) 0.00642324 (89, 42, 71) |
最小 最大 |
-0.00136636 (75, 42, 133) 0.00292404 (94, 65, 68) |
セグメンテーション
実行例1
最小 最大 |
-0.00136636 (75, 42, 133) 0.00292404 (94, 65, 68) |
最小 最大 |
-0.0012017 (5, 9, 13) 0.00263957 (13, 6, 15) |
トモグラム画像の解釈
また、多様な構造を含むため、3次元画像から興味あるセグメントを切り出すなどの作業が必要です。そのための支援ソフトウェアが必要です。
サブトモグラムの平均化
xy平面 | yz平面 | zx平面 |
実行例1
xy平面 |
yz平面 |
zx平面 |
最小 最大 |
-0.0012017 (5, 9, 13) 0.00263957 (13, 6, 15) |
xy平面 |
yz平面 |
zx平面 |
最小 最大 |
-0.000536652 (6, 13, 12) 0.00177258 (13, 17, 13) |