View source for Electron Tomography

'''Electron Tomography''' is one of 3D Reconstruction  using electron microscope. It is the method that reconstruct 3D image by electron micrographs which is projected the same photo field for several angle, then it creates Tomogram by using computer.

　In the following, describe its method in order.

== [[take tilt-series images]] ==
<div>　When sample is taken by this method, image that is tilted at a regular interval around an axis is made. As tilt angle increases, the photo field become wider, thus objects in image come around to the center gradually.</div>
<br>

	<table> 
		<tr> 
			<td>[[File:Input-Tomogram1.png]]<br>
				<p align="Center">1 Axis Tilt(Center is 0°'s image)<br></p>
			</td> 
		</tr> 
	</table> 
<br>

<div>When you will take 2 Axis Tilt 's image, rotate the sample 90 degree, and take as well as 1 Axis Tilt.</div>
	<table> 
		<tr> 
			<td>[[File:Input-Tomogram2.png]]<br>
				<p align="Center">2 Axis Tilt(Center is 0°'s image)(as 90° rotation image for 1 Axis Tilt's image)<br></p>
			</td> 
		</tr> 
	</table> 
<br>

== [[Image Correction]] ==
　Tilt-series image's defocus is larger than normal electron photograph's one, because it is taken as whole field is underfocus. In this case, CTF Correction is not required, because outside data of 1st thon ring don't have almost information if you use electron gun as LaB6. But, in the case of using field-emission electron gun (FEG), note that you might get a incorrect data.<br>
<br>

== [[Rough Alignment]] ==
　Align roughly each tilt image by using the correlation among images around center.

<div>[[Media:1WDC-Tom-2dSet.zip|Input file]](rotation at range: ±60°, interval: 2°, around Y-axis)<br>
	<table> 
		<tr> 
			<td><p align="Center">[[File:1WDC-Tom-2dSet.png]]<br>
				10° interval<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

=== Preprocess(Windowing) ===
<div>　By Windowing, remove noise of edge of image.</div>
<br>

<div>Windowing of this time is x: 80%, y: 90%</div>
<div>Command(In the case of 1st image): [[mrcImageWindowing]] -i Set1000.roi -o Set1000.mask -W 0.2 0.2 0.1 0.1</div>
<br>

<div>When you use [[Media:Makefile-Tomogram.zip|This Makefile]], it performs for all roi file.</div>
<div>Setting data</div>
<pre>
# For Windowing
WIN_X=0.2
WIN_X_MAX=0.2
WIN_Y=0.1
WIN_Y_MAX=0.1
WIN_MODE=18
</pre>

<div>Command</div>
<pre>
make Windowing
</pre>
<br>

<div>[[Media:1WDC-Tom-mask-2dSet.zip|Output file]](.mask)<br>
	<table> 
		<tr> 
			<td><p align="Center">[[File:1WDC-Tom-mask-2dSet.png]]<br>
				10° interval<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

===Alignment===
<div>　Because objects move parallelly depending on tilt angle, align about it. Now, align among the images whose angle is close to each other so as not to be affected by tilt.<br>
e.g. case of images whose interval is 2 degree<br>
　Align 2° image to 0° image.<br>
　Align 4° image to '''Aligned''' 2° image.<br>
　Align 6° image to '''Aligned''' 4° image.<br>
　.<br>
　.<br>
　(Align "-" as well as "+")</div>
<br>

====case of  using [[mrcImageCorrelation]]====
<div>　Perform alignment by using [[Media:Makefile-Tomogram.zip|This Makefile]]. Even if you don't use the Makefile, you perform command sequentially, and you can align. But, Input file list for [[mrc2Dto3D]] is required to be created.</div>
<br>

=====Example 1=====
<div>[[Media:1WDC-Tom-mask-2dSet.zip|Input file]]<br>
	<table> 
		<tr> 
			<td><p align="Center">[[File:1WDC-Tom-mask-2dSet.png]]<br>
				10° interval<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>Command</div>
<pre>
make CorFit1
</pre>
<br>

<div>[[Media:1WDC-Tom-Fit-2d.zip|Output file]](.fit)</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:1WDC-Tom-Fit-2d.png]]<br>
				10° interval<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

== [[Fine Alignment]] ==
<div>　Align the position and the angle of tilt axis as possible.</div>
<br>

===Example 1(Alignment of tilt axis)===
<div>　Calculate tilt of axis by using [[mrcImageTiltAxisSearch]].</div>
<br>

<div>[[Media:1WDC-Tom-Tilt-2dSet.zip|Input file]](1 axis tilt data whose axis is tilted by 10 degree)<br>
	<table> 
		<tr> 
			<td><p align="Center">[[File:1WDC-Tom-Tilt-2dSet.png]]<br>
				10° interval<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>[[Media:Makefile-Tomogram1.zip|Used Makefile]]'s Setting</div>
<pre>
# 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
TILTN=10
TILTITER=1
TILTSCALE=5
</pre>
<br>

<div>Command</div>
<pre>
make TiltFit
</pre>
<br>

<div>Output of axis tilt(It is stored at .tiltinfofile)</div>
<pre>
9.895
</pre>
<br>

<div>difference about 3D Reconstruction</div>
	<table> 
		<tr> 
			<td>Before correction</td> 
		</tr> 
		<tr> 
			<td>[[File:Outdata-Tomogram4.png]]</td> 
			<td>[[File:Outdata1-Tomogram4.png]]</td> 
			<td>[[File:Outdata2-Tomogram4.png]]</td> 
		</tr> 
		<tr> 
			<td>After correction</td> 
		</tr> 
		<tr> 
			<td>[[File:Outdata-Tomogram5.png]]</td> 
			<td>[[File:Outdata1-Tomogram5.png]]</td> 
			<td>[[File:Outdata2-Tomogram5.png]]</td> 
		</tr> 
		<tr> 
			<td>xy-plane</td> 
			<td>yz-plane</td> 
			<td>zx-plane</td> 
		</tr> 
	</table> 
<br>

===Example2(Repeatedly Alignment of tilt axis)===
<div>Repeat to use [[mrcImageTiltAxisSearch]], referring tilt of output axis.</div>
<br>

<div>[[Media:1WDC-Tom-Tilt1-2dSet.zip|Input file]](1 axis tilt data whose axis is tilted by 1 degree)<br>
	<table> 
		<tr> 
			<td><p align="Center">[[File:1WDC-Tom-Tilt1-2dSet.png]]<br>
				10° interval<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>[[Media:Makefile-Tomogram2.zip|Used Makefile]]'s Setting</div>
<pre>
# For mrcImageTiltAxisSearch
IN_TILT_EXT=roi
# RotMode
ROTMODE=ZEYR
# Rot1
ROT1MIN=1
ROT1MAX=1
ROT1D=1
# 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
</pre>
<br>

<div>Command</div>
<pre>
make TiltFit
</pre>
<br>

<div>transition of tilt of axis</div>
<pre>
0.084
0.161
0.230
0.294
0.345

...

0.701
0.700
0.685
0.706
0.701
</pre>
<div>The tilt of axis approaches gradually 1.</div>
<br>

== [[3D Reconstruction]] ==
<div>　２次元画像のセットから３次元画像を再構成します.</div>
<br>

<div>[[Media:1WDC-shift1.3d|元の3Ddata]]'s image</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:1WDC-shift1-3d.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:1WDC-shift1-3d2.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:1WDC-shift1-3d3.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">0 (0, 0, 0)<br> 
				4 (31, 26, 26)<br> 
				0.0116844<br> 
				0.147033<br> 
				0.000287174<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

===[[mrc2Dto3D]]===
====実行例1(1軸傾斜での再構成)====
<div>　連続傾斜像の撮影で軸にズレがない場合は角度刻みをそのままInput fileリストにSettingします.</div>
<br>

<div>[[Media:1WDC-Tom-Fit-2d.zip|Input file]]<br>
	<table> 
		<tr> 
			<td><p align="Center">[[File:1WDC-Tom-Fit-2d.png]]<br>
				10°刻みで表示<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>コマンド: [[mrc2Dto3D]] -i Input.3dlst -o Input.3d -m 1</div>
<br>

<div>[[Media:Makefile-Tomogram.zip|こちらのMakefile]]を使用すれば,以下のコマンドでも実行できます.</div>
<div>Setting内容</div>
<pre>
# For Reconstruction
IN_3D_EXT=fit
# RotMode
ROTMODE=YOYS
# Rot1
ROT1MIN=-60
ROT1MAX=60
ROT1D=2
# Rot2
ROT2MIN=0
ROT2MAX=0
# Rot3
ROT3MIN=0
ROT3MAX=0
</pre>
<div>コマンド</div>
<pre>
make 3DList
make Input.3d
</pre>
<br>

<div>Output file(Input.3d)</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Outdata-mrc2Dto3D.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata1-mrc2Dto3D.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata2-mrc2Dto3D.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-0.00437076 (39, 34, 36)<br> 
				0.00799233 (37, 34, 36)<br> 
				3.45089e-06<br> 
				0.000342836<br> 
				6.69601e-07<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

====実行例2(2軸傾斜での再構成)====
<div>２軸で傾斜した画像にて３次元再構成を行います.</div>
<div>[[Media:1WDC-Tom1-2dSet.zip|Input file]](X, Y-axis周りに±60°の範囲でそれぞれ10°刻みで傾斜)<br>
	<table> 
		<tr> 
			<td><p align="Center">[[File:1WDC-Tom1-2dSet.png]]<br>
				1行目がX-axis傾斜,２行目がY-axis傾斜<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>[[Media:Makefile-Tomogram6.zip|こちらのMakefile]]を使って３次元再構成を行います.</div>
<div>コマンド</div>
<pre>
make all
</pre>
<br>

<div>途中で使用する範囲を決めるために[[Display2]]が起動しますが,今回は全ての範囲を選択してROIInformationを保存します.<br>
今回OutputしたInformationfile</div>
<pre>
DataA_006.mrcsmth-0000.roi Rect        0.000000       0.000000       63.000000       0.000000       63.000000       63.000000        0.000000       63.000000
</pre>
<br>

<div>Output file</div>
	<table> 
		<tr>
			<td><p align="Center">[[File:Outdata-Tomogram2.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata1-Tomogram2.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata2-Tomogram2.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-9.32199 (29, 37, 34)<br> 
				21.995 (40, 35, 32)<br> 
				0.0280746<br> 
				 1.38771<br> 
				0.00271036<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

====実行例3(Doubleを使用した2軸傾斜での再構成)====
<div>[[mrc2Dto3D]]のOption -Doubleを使用して３次元再構成を行います.</div>
<br>

<div>変更箇所</div>
<pre>
.roilst.mrc3d:
#	mrc2Dto3D -I $*.roilst -o $*.mrc3d -single 0 -InterpolationMode 2 -m 1 
	mrc2Dto3D -I $*.roilst -o $*.mrc3d -Double -InterpolationMode 2 -m 1 
</pre>
<br>

<div>Output file</div>
	<table> 
		<tr>
			<td><p align="Center">[[File:Outdata-Tomogram6.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata1-Tomogram6.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata2-Tomogram6.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-10360.9 (41, 24, 0)<br> 
				8132.05 (40, 35, 31)<br> 
				86.0574<br> 
				673.594<br> 
				1.31561<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

===ラドン変換を使った３次元再構成===
<div>　[[mrc2Dto3D]]を使用するほかにラドン変換から３次元像を作る方法があります.ラドン変換は以下の流れで行います.
<br>

位置合わせ済みの2Dリスト<br>
↓[[mrcImageSinogramCreate]]<br>
Sinogramリスト<br>
↓[[mrcRadon2Dto3D]]<br>
3Dラドンfile<br>
↓[[mrcImageInverseRadonTransform]]<br>
3Dfile(完了)</div>
<br>

====実行例1====
<div>[[Media:Makefile-Tomogram7.zip|こちらのMakefile]]を使い,[[Media:1WDC-Tom2-2dSet.zip|実行例2でCreateされた .roilstfile]]について実行します.</div>
<br>

<div>コマンド</div>
<pre>
make Radon3D
</pre>

<div>Output file</div>
	<table> 
		<tr> 
			<td>xy-plane</td> 
			<td>yz-plane</td> 
			<td>zx-plane</td> 
		</tr> 
		<tr> 
			<td>[[File:Outdata-Tomogram3.png]]</td> 
			<td>[[File:Outdata1-Tomogram3.png]]</td> 
			<td>[[File:Outdata2-Tomogram3.png]]</td> 
		</tr> 
	</table> 
	<table> 
		<tr>
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">1.84923e+06 (0, 59, 28)<br> 
				1.17991e+07 (34, 31, 30)<br> 
				3.63631e+06<br> 
				1.64358e+06<br> 
				3536.42<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

====実行例2====
<div>　[[Media:Makefile-Tomogram7.zip|こちらのMakefile]]から逆ラドン変換のOptionを変更して実行します.</div>
<br>

<div>変更内容</div>
<pre>
### RadonTransform
RBP_MODE=4
RBP_OPTION=-Rmin 0.05 -Rmax 0.1
</pre>

<div>コマンド</div>
<pre>
make Radon3D
</pre>

<div>Output file</div>
	<table> 
		<tr> 
			<td>xy-plane</td> 
			<td>yz-plane</td> 
			<td>zx-plane</td> 
		</tr> 
		<tr> 
			<td>[[File:Outdata-Tomogram7.png]]</td> 
			<td>[[File:Outdata1-Tomogram7.png]]</td> 
			<td>[[File:Outdata2-Tomogram7.png]]</td> 
		</tr> 
	</table> 
	<table> 
		<tr>
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-344789 (25, 25, 0)<br> 
				631546 (37, 34, 28)<br> 
				4381.1<br> 
				76358.4<br> 
				164.297<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

== [[Problems of Electron Tomography image]] ==
===ミッシングエリア===
<div>　全方位からのProjection像を得られないために,ミッシングエリア(１軸の場合はウェッジ,２軸の場合はピラミッド)とよばれるInformationがない領域があります.そのため,方向依存でのボケを生じます.</div>
<br>

<div>[[:Media:Makefile-Tomogram5.zip|使用MakeFile]].</div>
<div>[[:Media:Input-ellipsoidal.3d|元の３次元像]](Target.ini3dにリネーム)</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Input-ellipsoidal.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input1-ellipsoidal.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input2-ellipsoidal.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">0 (0, 0, 0)<br> 
				2 (29, 29, 14)<br> 
				0.108302<br> 
				0.364429<br> 
				0.000711774<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>１軸回転でProjectionして３次元再構成した場合(Y-axis: -60° ~ 60°: 10°刻み)</div>
<div>変更箇所</div>
<pre>
TILTAXISNUMBER=1	# Single: 1 Double: 2 
</pre>
<br>

<div>コマンド</div>
<pre>
make Target.ini2d
make TestData2DSet
make all
</pre>
<br>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Output-ellipsoidal.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Output1-ellipsoidal.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Output2-ellipsoidal.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-0.00330818 (52, 36, 32)<br> 
				0.00425516 (13, 37, 32)<br> 
				2.08481e-05<br> 
				0.000467848<br> 
				9.13766e-07<br></p> 
			</td> 
		</tr> 
	</table> 
<div></div>
<br>

<div>２軸回転でProjectionして３次元再構成した場合(X, Y-axis: 各-60° ~ 60°: 10°刻み)</div>
<div>変更箇所</div>
<pre>
TILTAXISNUMBER=2	# Single: 1 Double: 2 
</pre>
<div>コマンド</div>
<pre>
make Target.ini2d1
make TestData2DSetDouble
make all
</pre>
<br>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Output3-ellipsoidal.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Output4-ellipsoidal.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Output5-ellipsoidal.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-0.00248567 (12, 38, 32)<br> 
				0.00344027 (13, 37, 32)<br> 
				2.07529e-05<br> 
				0.000401237<br> 
				7.83666e-07<br></p> 
			</td> 
		</tr> 
	</table> 
<div></div>
<br>

===画像の位置,フォーカスの問題===
<div>　試料を傾けながら撮影するため回転軸からの距離に応じてレンズとの距離が変化していきます.これによって生じる問題をいかに補正するかも一つのポイントとなります.</div>
<br>

====parallel transformの補正====
<div>　傾斜する角度が大きくなるほど,視野が広くなるため試料は傾斜軸の方へ近づいていきます.これを補正するために[[mrcImageCorrelation]]などを使ってどれだけparallel transformが必要かを計算します.また,試料を傾斜する毎に形も変わっていきますので比較はできるだけ近い画像同士で行います.</div>
<br>

====軸の決定====
<div>　撮影した画像は必ずしもSettingした軸を中心に傾斜しているとは限りません.これは電子顕微鏡自体の性能や試料の設置方法などによって多少のズレが生じてしまうからです.傾斜軸のズレは[[mrcImageTiltAxisSearch]]などを使って計算します.</div>
<br>

== トモグラフに使われる画像処理 ==
<div>　再構成した３次元画像のAnalysisを行う前処理として以下のような処理を施します.</div>
<br>

=== Smoothing ===
<div>　再構成した３次元像のNoiseを除去するためSmoothingを行います.</div>
<br>

====実行例1====
<div>[[Media:Input-Tomogram3D.mrc|Input file]]</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Input-Tomogram3D.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input1-Tomogram3D.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input2-Tomogram3D.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-0.00268221 (76, 92, 66)<br> 
				0.00642324 (89, 42, 71)<br> 
				2.0386e-06<br> 
				0.000254286<br> 
				1.63933e-07<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>コマンド: [[mrcImageSmoothing]] -i Input-Tomogram3D.mrc -o Input-TomogramSmth.mrc -m 1</div>
<br>

<div>[[Media:Input-TomogramSmth.mrc|Output file]]</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Input-TomogramSmth.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input1-TomogramSmth.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input2-TomogramSmth.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-0.00136636 (75, 42, 133)<br> 
				0.00292404 (94, 65, 68)<br> 
				-9.48921e-07<br> 
				0.00013513<br> 
				8.71152e-08<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

=== セグメンテーション ===
<div>　２次元画像や３次元画像を必要な部分ごとに切り出します.</div>
<br>

====３次元画像の切り出し====
<div>　再構成された３次元画像から粒子と考えられる部分などを切り出します.切り出したdataは３次元画像のAnalysisや平均化するときに使用することができます.</div>
<br>

=====実行例1=====
<div>切り出したい部分のcoordinates,Widthが分かっている場合は[[mrcImageCenterGet]]を使います.</div>
<br>

<div>[[Media:Input-TomogramSmth.mrc|Input file]]</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Input-TomogramSmth.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input1-TomogramSmth.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input2-TomogramSmth.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-0.00136636 (75, 42, 133)<br> 
				0.00292404 (94, 65, 68)<br> 
				-9.48921e-07<br> 
				0.00013513<br> 
				8.71152e-08<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>コマンド: [[mrcImageCenterGet]] -i Input-TomogramSmth.mrc -o Input-TomogramSub.mrc -Cx 67 -Cy 67 -Cz 67 -Nx 27 -Ny 27 -Nz 27</div>
<br>

<div>[[Media:Input-TomogramSub.mrc|Output file]]</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Input-TomogramSub.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input1-TomogramSub.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input2-TomogramSub.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-0.0012017 (5, 9, 13)<br> 
				0.00263957 (13, 6, 15)<br> 
				2.22923e-05<br> 
				0.000388094<br> 
				2.76625e-06<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

====２次元画像の切り出し====
<div>　視野全体ではなく,一部だけ再構成するために傾斜画像から部分的に２次元画像を切り出す場合もあります.視野全体の再構成よりもコストが削減されるので,部分的なdataだけあれば十分な場合はこの手法を用います.試料を傾けながら撮影しているので,画像から切り出す位置や角度を画像によって変える必要があります.切り出すときには[[mrcImageCenterGet]]や[[mrcImageROI]],[[mrcImageROIs]]などを使用します.切り出し範囲を全て手動で決めるのはかなりの時間を要しますので,基準(0°)'s imageでの切り出し範囲を手動で決めると残りの切り出し位置を傾斜角度に応じて自動で計算して切り出しを行う方法をとります.</div>
操作動画: [[:Media:Tomography4-1.mov|(.mov)]] [[:Media:Tomography4-1.mp4|(.mp4)]]<br>
<br>

<div>以下の流れで３次元再構成まで行います.</div>
<pre>
傾斜画像(2D画像複数)(前処理済)

↓	切り出し範囲のSetting(Display2: 0°'s imageのみ)

0°の切り出しInformationfile(ROIInformation)

↓	各傾斜画像の切り出し範囲を計算する(現在はMakefile内で処理)

各傾斜画像の切り出しInformation

↓	画像の切り出し(mrcImageROIs)

切り出された複数's imagefile(切り出し数 × 傾斜数分)

↓	位置合わせ(mrcImageCorrelarion + Makefileの処理)

位置合わせ済みの切り出しInformation(傾斜数分)

↓	再度画像の切り出し(mrcImageROIs) (※ 切り出し画像をシフトするのではなく,切り出し位置を再Settingして切り出します)

位置合わせ済み's image(ROI)file(切り出し数 × 傾斜数分)

↓	角度InformationfileのCreate(Makefile内の処理)

角度Informationfile(切り出し数分)

↓	３次元再構成(mrc2Dto3D)

3次元画像(切り出し数分)
</pre>

=====実行例1(1軸傾斜のサブトモグラム)=====
<div>[[Media:Input-Tomogram2DSet1.zip|Input file]]'s image(一部)(y-axis回転: -60° ~ 60°: 2°刻み)</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Input-Tomogram2DSet1.png]]<br>
				-60°(p1_001.mrc)<br></p> 
			</td> 
			<td><p align="Center">[[File:Input1-Tomogram2DSet1.png]]<br>
				0°(p1_031.mrc)<br></p> 
			</td> 
			<td><p align="Center">[[File:Input2-Tomogram2DSet1.png]]<br>
				60°(p1_061.mrc)<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>[[:Media:Makefile-Tomogram4.zip|こちらのMakefile]].</div>
<div>コマンド</div>
<pre>
make all
</pre>

<div>コマンドをInput後に[[Display2]](2Dビューワー)が開きますので,切り出したい部分(複数可)を選択しInformationfile(.roiinfo)をCreateします.</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Tomogram1.png]]</p></td> 
			<td><p align="Center">[[File:Tomogram2.png]]</p></td> 
		</tr> 
		<tr> 
			<td><p align="Center">切り出し部分決めてEdit->OKで決定します.</p></td> 
			<td><p align="Center">複数切り出したいときはROI->MultiROIを選択します.</p> </td> 
		</tr> 
	</table> 
	<table> 
		<tr> 
			<td><p align="Center">[[File:Tomogram3.png]]</p></td> 
		</tr> 
		<tr> 
			<td>今回はInformationだけCreateすれば,切り出しは自動で行われます.<br>
			ROI InformationウィンドウにてInfoFileName右側のsaveボタンを押すとfileがCreateされます.<br>
			(またはExtractModeをInformationにしてExtractボタンを押します.)</td> 
		</tr> 
	</table> 
<br>

<div>今回OutputしたInformationfile</div>
<pre>
p1_031-0000.roi	Rect	20	30	60	30	60	70	20	70
p1_031-0001.roi	Rect	15	75	55	75	55	115	15	115
p1_031-0002.roi	Rect	25	110	65	110	65	150	25	150
p1_031-0003.roi	Rect	45	0	85	0	85	40	45	40
p1_031-0004.roi	Rect	60	50	100	50	100	90	60	90
p1_031-0005.roi	Rect	55	80	95	80	95	120	55	120
p1_031-0006.roi	Rect	55	120	95	120	95	160	55	160
p1_031-0007.roi	Rect	85	15	125	15	125	55	85	55
p1_031-0008.roi	Rect	85	51	125	51	125	91	85	91
p1_031-0009.roi	Rect	85	91	125	91	125	131	85	131
p1_031-0010.roi	Rect	93	130	133	130	133	170	93	170
p1_031-0011.roi	Rect	115	25	155	25	155	65	115	65
p1_031-0012.roi	Rect	120	60	160	60	160	100	120	100
p1_031-0013.roi	Rect	123	100	163	100	163	140	123	140
p1_031-0014.roi	Rect	115	140	155	140	155	180	115	180
p1_031-0015.roi	Rect	152	152	192	152	192	192	152	192
</pre>
<br>

<div>1枚's imageについて[[Display2]]からInformationfile(.roiinfo)をCreateすると,残り's imageについては傾きから切り出す範囲を計算し,自動で切り出しが行われます.</div>
<br>
<div>自動で切り出された画像(縦: 選択範囲毎, 横: file(角度)毎)(.pad)</div>
	<table> 
		<tr> 
			<td></td> 
			<td>[[File:Tomogram-Segmentation.png]]</td>
			<td></td> 
		</tr> 
		<tr> 
			<td>-60°</td> 
			<td><p align="Center">0°</p></td> 
			<td>60°</td> 
		</tr> 
	</table> 
<br>

<div>切り出し処理を終えた後はそれぞれ's imageに対して３次元再構成が行われます.</div>
<br>

<div>再構成されたそれぞれの３次元画像(.mrc3d)</div>
	<table> 
		<tr> 
			<td>xy-plane</td> 
			<td>[[File:Outdata-Tomogram.png]]</td>
		</tr> 
		<tr> 
			<td>yz-plane</td> 
			<td>[[File:Outdata1-Tomogram.png]]</td>
		</tr> 
		<tr> 
			<td>zx-plane</td> 
			<td>[[File:Outdata2-Tomogram.png]]</td>
		</tr> 
	</table>
<br>

=====実行例2(2軸傾斜のサブトモグラム)=====
<div>　２軸で傾斜した画像から切り出し及び再構成を行ってみます.</div>
<br>

<div>[[Media:Input-Tomogram2DSet2.zip|Input file]]'s image(一部)(x-axis, y-axis回転: -60° ~ 60°: 10°刻み)</div>
	<table> 
		<tr> 
			<td>x-axis回転</td>
			<td><p align="Center">[[File:Input-Tomogram2DSet2.png]]<br>
				-60°<br></p> 
			</td> 
			<td><p align="Center">[[File:Input1-Tomogram2DSet2.png]]<br>
				0°<br></p> 
			</td> 
			<td><p align="Center">[[File:Input2-Tomogram2DSet2.png]]<br>
				60°<br></p> 
			</td> 
		</tr> 
		<tr> 
			<td>y-axis回転</td>
			<td><p align="Center">[[File:Input3-Tomogram2DSet2.png]]<br>
				-60°<br></p> 
			</td> 
			<td><p align="Center">[[File:Input4-Tomogram2DSet2.png]]<br>
				0°<br></p> 
			</td> 
			<td><p align="Center">[[File:Input5-Tomogram2DSet2.png]]<br>
				60°<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>[[:Media:Makefile-Tomogram5.zip|こちらのMakefile]].</div>
<br>

<div>コマンド</div>
<pre>
make all
</pre>
<br>

<div>今回はこの領域でROIfileの切り出しを行います.</div>
<div>[[File:Tomogram2-Segmanetation.png]]</div>
<div>今回OutputしたInformationfile</div>
<pre>
DataA_006-0000.roi Rect        0.722314      137.147614       33.000000      106.000000       49.053597      122.636041       16.775911      153.783654
DataA_006-0001.roi Rect       15.000000       79.000000       42.000000       79.000000       42.000000      107.000000       15.000000      107.000000
DataA_006-0002.roi Rect       13.000000       39.000000       58.000000       39.000000       58.000000       64.000000       13.000000       64.000000
DataA_006-0003.roi Rect       41.000000      129.000000       84.000000      129.000000       84.000000      154.000000       41.000000      154.000000
DataA_006-0004.roi Rect       51.000000       85.000000       82.000000       85.000000       82.000000      119.000000       51.000000      119.000000
DataA_006-0005.roi Rect       47.269722       69.494676       78.455043       40.316766       95.000000       58.000000       63.814679       87.177909
DataA_006-0006.roi Rect       66.000000        4.000000       92.000000        4.000000       92.000000       44.000000       66.000000       44.000000
DataA_006-0007.roi Rect       77.597890      150.922358       95.484138      132.052220      120.886248      156.129861      103.000000      175.000000
DataA_006-0008.roi Rect       93.000000       92.000000      115.000000       92.000000      115.000000      132.000000       93.000000      132.000000
DataA_006-0009.roi Rect      100.000000       51.000000      121.000000       51.000000      121.000000       93.000000      100.000000       93.000000
DataA_006-0010.roi Rect       98.589377       19.087800      119.760822       10.100593      135.000000       46.000000      113.828555       54.987208
DataA_006-0011.roi Rect      122.000000      141.000000      146.000000      141.000000      146.000000      184.000000      122.000000      184.000000
DataA_006-0012.roi Rect      136.883039       95.213292      160.000000      104.000000      144.116961      145.786708      121.000000      137.000000
DataA_006-0013.roi Rect      127.000000       73.000000      146.061183       60.635990      169.000000       96.000000      149.938817      108.364010
DataA_006-0014.roi Rect      140.000000       22.000000      177.361289       43.286947      166.470292       62.402025      129.109004       41.115078
DataA_006-0015.roi Rect      163.000000      152.000000      183.000000      152.000000      183.000000      190.000000      163.000000      190.000000

</pre>
<br>

<div>自動で切り出された画像(縦: 選択範囲毎, 横: file(角度)毎)(.pad)</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Tomogram2-Segmanetation1.png]]<br>
				DataA(x-axis回転)<br></p> 
			</td> 
			<td><p align="Center">[[File:Tomogram2-Segmanetation2.png]]<br>
				DataB(y-axis回転)<br></p> 
			</td> 
		</tr> 
	</table> 
<div>切り出しサイズがバラバラで分かりにくくなっていますが,いずれの傾斜軸でもROIfileが自動で切り出されています.</div>
<br>

<div>再構成されたそれぞれの３次元画像(.mrc3d)</div>
	<table> 
		<tr> 
			<td>xy-plane</td> 
			<td>[[File:Outdata-Tomogram1.png]]</td>
		</tr> 
		<tr> 
			<td>yz-plane</td> 
			<td>[[File:Outdata1-Tomogram1.png]]</td>
		</tr> 
		<tr> 
			<td>zx-plane</td> 
			<td>[[File:Outdata2-Tomogram1.png]]</td>
		</tr> 
	</table>
<br>

== [[Interpretation of 3D Reconstruction]] ==
　また,多様な構造を含むため,３次元画像から興味あるセグメントを切り出すなどの作業が必要です.そのための支援ソフトウェアが必要です.

=== [[Averaging of Subtomograms]]===
<div>　セグメンテーションによって分割された３次元data同士が別方向を向いている同じ粒子であるとき,それを平均化することでミッシングエリアを補間し,より精密な３次元像にすることができます.</div>
<br>

<div>分割したdata1</div>
	<table> 
		<tr> 
			<td>[[File:Input-SubTomogram.png]]</td> 
			<td>[[File:Input1-SubTomogram.png]]</td> 
			<td>[[File:Input2-SubTomogram.png]]</td> 
		</tr> 
	</table> 
<div>分割したdata2(90°向きが違う)</div>
	<table> 
		<tr> 
			<td>[[File:Input-SubTomogram1.png]]</td> 
			<td>[[File:Input1-SubTomogram2.png]]</td> 
			<td>[[File:Input2-SubTomogram3.png]]</td> 
		</tr> 
	</table> 
<div>平均したdata</div>
	<table> 
		<tr> 
			<td>[[File:Output-SubTomogram1.png]]</td> 
			<td>[[File:Output1-SubTomogram2.png]]</td> 
			<td>[[File:Output2-SubTomogram3.png]]</td> 
		</tr> 
		<tr> 
			<td>xy-plane</td> 
			<td>yz-plane</td> 
			<td>zx-plane</td> 
		</tr> 
	</table> 
<br>

====実行例1====
<div>[[Media:Input-Mean-SubTomogram.zip|Input file]]</div>
	<table> 
		<tr> 
			<td>xy-plane</td> 
			<td>[[File:Input-Mean-SubTomogram.png]]</td> 
		</tr> 
		<tr> 
			<td>yz-plane</td> 
			<td>[[File:Input1-Mean-SubTomogram.png]]</td> 
		</tr> 
		<tr> 
			<td>zx-plane</td> 
			<td>[[File:Input2-Mean-SubTomogram.png]]</td> 
		</tr> 
	</table> 
<br>

<div>３次元像を平均化する前に[[mrcImageAutoRotationCorrelation3D]]を使って,向きを合わせます.</div>
<br>

<div>[[Media:Input-TomogramSub.mrc|Reference file]](Input fileの真ん中のdataを使用)</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Input-TomogramSub.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input1-TomogramSub.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input2-TomogramSub.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-0.0012017 (5, 9, 13)<br> 
				0.00263957 (13, 6, 15)<br> 
				2.22923e-05<br> 
				0.000388094<br> 
				2.76625e-06<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>コマンド: [[mrcImageAutoRotationCorrelation3D]] -i c1-1.mrc -r c2-2.mrc -fit f1-1.mrc -EA YOYS -Rot1 0 0 30 -Rot2 0 330 30 -Rot3 0 330 30 -Xrange 0 0 -Yrange 0 0  -Zrange 0 0  </div>
<br>

<div>[[Media:Input-Fit-SubTomogram.zip|補正後のdata]](-fitのOutput file)</div>
	<table> 
		<tr> 
			<td>xy-plane</td> 
			<td>[[File:Input-Fit-SubTomogram.png]]</td> 
		</tr> 
		<tr> 
			<td>yz-plane</td> 
			<td>[[File:Input1-Fit-SubTomogram.png]]</td> 
		</tr> 
		<tr> 
			<td>zx-plane</td> 
			<td>[[File:Input2-Fit-SubTomogram.png]]</td> 
		</tr> 
	</table> 
<br>

<div>補正fileに対し,[[mrcImageAverage]]にて平均化を行います.</div>
<div>Output file</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Outdata-TomogramSub.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata1-TomogramSub.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata2-TomogramSub.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-0.000536652 (6, 13, 12)<br> 
				0.00177258 (13, 17, 13)<br> 
				3.22171e-05<br> 
				0.000228739<br> 
				1.6304e-06<br></p> 
			</td> 
		</tr> 
	</table> 
<div>平均化することにより３次元像の質(特にyz-plane)が良くなりました.</div>
<br>

====実行例2====
<div>[[Media:Outdata-Tomogram.zip|Input file]]</div>
	<table> 
		<tr> 
			<td>xy-plane</td> 
			<td>[[File:Outdata-Tomogram.png]]</td> 
		</tr> 
		<tr> 
			<td>yz-plane</td> 
			<td>[[File:Outdata1-Tomogram.png]]</td> 
		</tr> 
		<tr> 
			<td>zx-plane</td> 
			<td>[[File:Outdata2-Tomogram.png]]</td> 
		</tr> 
	</table> 
<br>

<div>[[Media:Input-TomogramSub1.mrc|Reference file]](Input fileの先頭のdataを使用)</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Input-TomogramSub1.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input1-TomogramSub1.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Input2-TomogramSub1.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-5.47817 (34, 19, 20)<br> 
				16.0715 (31, 22, 19)<br> 
				0.142234<br> 
				1.22547<br> 
				0.0048441<br></p> 
			</td> 
		</tr> 
	</table> 
<br>

<div>同様にまず向きを合わせます.</div>
<div>コマンド: [[mrcImageAutoRotationCorrelation3D]] -i p1-0001.mrc3d -r Target.mrc3dref -fit p1-0001.3dfit -EA XEYR -Rot1 0 355 5 -Rot2 0 355 5 -Rot3 0 0 30 -M 18</div>
<br>
<div>[[Media:Input-Fit-SubTomogram1.zip|補正後のdata]](-fitのOutput file)</div>
	<table> 
		<tr> 
			<td>xy-plane</td> 
			<td>[[File:Input-Fit-SubTomogram1.png]]</td> 
		</tr> 
		<tr> 
			<td>yz-plane</td> 
			<td>[[File:Input1-Fit-SubTomogram1.png]]</td> 
		</tr> 
		<tr> 
			<td>zx-plane</td> 
			<td>[[File:Input2-Fit-SubTomogram1.png]]</td> 
		</tr> 
	</table> 
<br>

<div>Reference fileに向きを合っているものと合っていないものがあります.<br>
これはそれぞれのサブトモグラムの３次元像にてミッシングエリアが異なるので,どちらかにdataがない角度では相関値が大Widthに下がってしまうためです.</div>
<div>[[File:Fig-SubTomogram.png]]</div>
<br>

<div>補正fileに対し,[[mrcImageAverage]]にて平均化を行います.</div>
<div>Output file</div>
	<table> 
		<tr> 
			<td><p align="Center">[[File:Outdata-TomogramSub1.png]]<br>
				xy-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata1-TomogramSub1.png]]<br>
				yz-plane<br></p> 
			</td> 
			<td><p align="Center">[[File:Outdata2-TomogramSub1.png]]<br>
				zx-plane<br></p> 
			</td> 
			<td><p align="left">Min<br> 
				Max<br> 
				Mean<br> 
				SD<br> 
				SE<br></p> 
			</td> 
			<td><p align="left">-2.32263 (22, 29, 20)<br> 
				5.20816 (19, 24, 22)<br> 
				0.0844992<br> 
				0.49404<br> 
				0.00195286<br></p> 
			</td> 
		</tr> 
	</table> 
<br>