「3次元再構成・逆投影」の版間の差分
(→電子線トモグラフィー) |
(→ラドン変換による方法) |
||
(2人の利用者による、間の7版が非表示) | |||
行38: | 行38: | ||
PRJ_MODE=0 : -mの値 | PRJ_MODE=0 : -mの値 | ||
</pre> | </pre> | ||
+ | <br> | ||
+ | |||
+ | 2軸で逆投影する場合は下記を変更します。<br> | ||
+ | <pre> | ||
+ | # mrc2Dto3D -I $*.3dlst -o $*.3d -InterpolationMode $(IN_MODE) -m $(PRJ_MODE) -single $(AXIS1) $(SIRT); | ||
+ | mrc2Dto3D -I $*.3dlst -o $*.3d -InterpolationMode $(IN_MODE) -m $(PRJ_MODE) -Double $(SIRT); | ||
+ | </pre> | ||
+ | <br> | ||
+ | |||
<div>makefileでのコマンド</div> | <div>makefileでのコマンド</div> | ||
<pre> | <pre> | ||
行63: | 行72: | ||
<td>フィルタ逆投影</td> | <td>フィルタ逆投影</td> | ||
<td>-m 1</td> | <td>-m 1</td> | ||
+ | <td>[[画像:Outdata-BP-2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-2-2.png]]</td> | ||
+ | <td>[[ファイル:Outdata-BP-2.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(2軸)</td> | ||
+ | <td>-m 1 -Double</td> | ||
<td>[[画像:Outdata-BP-2-0.png]]</td> | <td>[[画像:Outdata-BP-2-0.png]]</td> | ||
<td>[[画像:Outdata-BP-2-1.png]]</td> | <td>[[画像:Outdata-BP-2-1.png]]</td> | ||
行73: | 行90: | ||
=====電子線トモグラフィー===== | =====電子線トモグラフィー===== | ||
− | [[Media:Input-BP-3.zip|こちらのファイル]](makefile込み) | + | [[Media:Input-BP-3.zip|こちらのファイル]](makefile込み)を入力ファイルとします。±60度で、10度刻みの画像からの3次元再構成を示しています。 |
<br> | <br> | ||
行96: | 行113: | ||
<td>フィルタ逆投影</td> | <td>フィルタ逆投影</td> | ||
<td>-m 1</td> | <td>-m 1</td> | ||
+ | <td>[[画像:Outdata-TBP-1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-1.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(2軸)</td> | ||
+ | <td>-m 1 -Double</td> | ||
<td>[[画像:Outdata-TBP-1-0.png]]</td> | <td>[[画像:Outdata-TBP-1-0.png]]</td> | ||
<td>[[画像:Outdata-TBP-1-1.png]]</td> | <td>[[画像:Outdata-TBP-1-1.png]]</td> | ||
行133: | 行158: | ||
<td>[[画像:Outdata-2BP-2-2.png]]</td> | <td>[[画像:Outdata-2BP-2-2.png]]</td> | ||
<td>[[画像:Outdata-2BP-2.gif]]</td> | <td>[[画像:Outdata-2BP-2.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(2軸)</td> | ||
+ | <td>-m 1 -Double</td> | ||
+ | <td>[[画像:Outdata-2BP-d2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d2-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d2.gif]]</td> | ||
</tr> | </tr> | ||
</table> | </table> | ||
行165: | 行198: | ||
<td>[[画像:Outdata-T2BP-1-2.png]]</td> | <td>[[画像:Outdata-T2BP-1-2.png]]</td> | ||
<td>[[画像:Outdata-T2BP-1.gif]]</td> | <td>[[画像:Outdata-T2BP-1.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(2軸)</td> | ||
+ | <td>-m 1 Double</td> | ||
+ | <td>[[画像:Outdata-T2BP-d1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d1.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | ===重み付き逆投影法=== | ||
+ | 重み付き逆投影法による3次元再構成法をするためのプログラムです。 | ||
+ | |||
+ | <pre> | ||
+ | .3dlst.3d: | ||
+ | mrc2Dto3D -I $(TARGET).3dlst -o $(TARGET).3d -InterpolationMode 2 -Double -DoubleCounter $(TARGET).3dcounter -CounterThreshold 0.5 -m 1 -WeightMode 2 | ||
+ | </pre> | ||
+ | |||
+ | $(TARGET).3dcounterが重みを表現しています。$(TARGET).3d及び$(TARGET).3dcounterがあれば、[[mrcImageAverage]]を使って、3次元画像を平均することができます。 | ||
+ | |||
+ | ====1軸回転==== | ||
+ | =====単粒子解析===== | ||
+ | [[Media:Input-BP-1.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <div>[[mrc2Dto3D]]はオプションm=2で重み付き逆投影に設定できます。</div> | ||
+ | <div>makefileでの設定項目</div> | ||
+ | <pre> | ||
+ | IN_3D_TIFF_EXT=3d : 3次元画像をtiffにする | ||
+ | PRJ_MODE=2 : -mの値 | ||
+ | </pre> | ||
+ | <div>makefileでのコマンド</div> | ||
+ | <pre> | ||
+ | make 3D | ||
+ | make 3DTiff | ||
+ | </pre> | ||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>方法</th> | ||
+ | <th>対応オプション</th> | ||
+ | <th>xy平面</th> | ||
+ | <th>yz平面</th> | ||
+ | <th>zx平面</th> | ||
+ | <th>全体</th> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影</td> | ||
+ | <td>-m 0</td> | ||
+ | <td>[[画像:Outdata-BP-0-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-0-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-0-2.png]]</td> | ||
+ | <td>[[ファイル:Outdata-BP-0.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影</td> | ||
+ | <td>-m 2</td> | ||
+ | <td>[[画像:Outdata-BP-w-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-w-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-w-2.png]]</td> | ||
+ | <td>[[ファイル:Outdata-BP-w.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | =====電子線トモグラフィー===== | ||
+ | [[Media:Input-BP-3.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>方法</th> | ||
+ | <th>対応オプション</th> | ||
+ | <th>xy平面</th> | ||
+ | <th>yz平面</th> | ||
+ | <th>zx平面</th> | ||
+ | <th>全体</th> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影</td> | ||
+ | <td>-m 0</td> | ||
+ | <td>[[画像:Outdata-TBP-0-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-0-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-0-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-0.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影</td> | ||
+ | <td>-m 2</td> | ||
+ | <td>[[画像:Outdata-TBP-w-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-w-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-w-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-w.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | ====2軸回転==== | ||
+ | =====単粒子解析===== | ||
+ | [[Media:Input-BP-2.zip|こちらのファイル]](makefile込み)を入力ファイルとします。<br> | ||
+ | 2軸で逆投影する場合は下記を変更します。<br> | ||
+ | <pre> | ||
+ | # mrc2Dto3D -I $*.3dlst -o $*.3d -InterpolationMode $(IN_MODE) -m $(PRJ_MODE) -single $(AXIS1) $(SIRT); | ||
+ | mrc2Dto3D -I $*.3dlst -o $*.3d -InterpolationMode $(IN_MODE) -m $(PRJ_MODE) -Double $(SIRT); | ||
+ | </pre> | ||
+ | <br> | ||
+ | |||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>方法</th> | ||
+ | <th>対応オプション</th> | ||
+ | <th>xy平面</th> | ||
+ | <th>yz平面</th> | ||
+ | <th>zx平面</th> | ||
+ | <th>全体</th> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影</td> | ||
+ | <td>-m 0</td> | ||
+ | <td>[[画像:Outdata-2BP-0-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-0-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-0-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-0.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影</td> | ||
+ | <td>-m 2</td> | ||
+ | <td>[[画像:Outdata-2BP-w-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-w-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-w-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-w.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(2軸)</td> | ||
+ | <td>-m 2 -Double</td> | ||
+ | <td>[[画像:Outdata-2BP-wd-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | =====電子線トモグラフィー===== | ||
+ | [[Media:Input-BP-4.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>方法</th> | ||
+ | <th>対応オプション</th> | ||
+ | <th>xy平面</th> | ||
+ | <th>yz平面</th> | ||
+ | <th>zx平面</th> | ||
+ | <th>全体</th> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影</td> | ||
+ | <td>-m 0</td> | ||
+ | <td>[[画像:Outdata-T2BP-0-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-0-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-0-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-0.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影</td> | ||
+ | <td>-m 2</td> | ||
+ | <td>[[画像:Outdata-T2BP-w-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-w-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-w-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-w.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(2軸)</td> | ||
+ | <td>-m 2 -Double</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd.gif]]</td> | ||
</tr> | </tr> | ||
</table> | </table> | ||
行385: | 行598: | ||
<tr> | <tr> | ||
<td>フィルタ逆投影(xy平面)</td> | <td>フィルタ逆投影(xy平面)</td> | ||
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d2-0.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d3-0.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d5-0.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d7-0.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d9-0.png]]</td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
<td>フィルタ逆投影(yz平面)</td> | <td>フィルタ逆投影(yz平面)</td> | ||
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d2-1.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d3-1.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d5-1.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d7-1.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d9-1.png]]</td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
<td>フィルタ逆投影(zx平面)</td> | <td>フィルタ逆投影(zx平面)</td> | ||
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d2-2.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d3-2.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d5-2.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d7-2.png]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d9-2.png]]</td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
<td>フィルタ逆投影(全体)</td> | <td>フィルタ逆投影(全体)</td> | ||
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d2.gif]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d3.gif]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d5.gif]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d7.gif]]</td> |
− | <td>[[画像:Outdata-2BP- | + | <td>[[画像:Outdata-2BP-d9.gif]]</td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
行497: | 行710: | ||
<tr> | <tr> | ||
<td>フィルタ逆投影(xy平面)</td> | <td>フィルタ逆投影(xy平面)</td> | ||
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d1-0.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d3-0.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d5-0.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d7-0.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d9-0.png]]</td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
<td>フィルタ逆投影(yz平面)</td> | <td>フィルタ逆投影(yz平面)</td> | ||
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d1-1.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d3-1.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d5-1.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d7-1.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d9-1.png]]</td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
<td>フィルタ逆投影(zx平面)</td> | <td>フィルタ逆投影(zx平面)</td> | ||
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d1-2.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d3-2.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d5-2.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d7-2.png]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d9-2.png]]</td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
<td>フィルタ逆投影(全体)</td> | <td>フィルタ逆投影(全体)</td> | ||
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d1.gif]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d3.gif]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d5.gif]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d7.gif]]</td> |
− | <td>[[画像:Outdata-T2BP- | + | <td>[[画像:Outdata-T2BP-d9.gif]]</td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
行559: | 行772: | ||
<td>[[画像:Outdata-T2BP-wd30.gif]]</td> | <td>[[画像:Outdata-T2BP-wd30.gif]]</td> | ||
</tr> | </tr> | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
</table> | </table> | ||
<br> | <br> | ||
行739: | 行780: | ||
2次元画像のセット(複数枚) | 2次元画像のセット(複数枚) | ||
− | ↓ 各画像に[[mrcImageSinogramCreate]] | + | ↓ 各画像に[[mrcImageSinogramCreate]]にて[[シノグラム]](2次元ラドン空間)に変換する |
− | + | [[シノグラム]]のセット(複数枚) | |
↓ [[mrcRadon2Dto3D]]にて3次元ラドンにする | ↓ [[mrcRadon2Dto3D]]にて3次元ラドンにする | ||
行747: | 行788: | ||
3次元ラドンファイル | 3次元ラドンファイル | ||
− | ↓ [[mrcImageInverseRadonTransform]] | + | ↓ [[mrcImageInverseRadonTransform]]にて逆ラドン変換を行う |
3次元ファイル(再構成完了) | 3次元ファイル(再構成完了) |
2014年7月1日 (火) 06:05時点における最新版
3次元再構成法の概要
3次元再構成法にはいくつかの方法があります。(Eos)としたものは、Eosで対応している手法です。
1.単純逆投影法(Simple Backprojection) (Eos) 単純に逆投影を実施するものです。
2.重み付き逆投影法(Weighted Backprojection) ()
3.フィルター付逆投影法(Eos)
投影像の投影角に基づいて、フーリエ空間でフィルタをかけます。 ローフィルタというフィルタが1軸投影(トモグラフィー)などの場合にはよく用いられます。 2、3の区別は非常に難しいですが、原理的にはまったく違いがないと考えるべきでしょう.
4.繰り返し法(Iterative Method)
4.1 SIRT法(Simultaneous Iterative Reconstruction Technique)
4.2 ART法(Algebric Iterative Reconstruction Technique)(未対応)
5.ラドン変換による方法 (Eos)
演習では、重み付き逆投影法を用います。
3次元再構成の実際
フィルタ逆投影法
1軸回転
単粒子解析
こちらのファイル(makefile込み)を入力ファイルとします。
IN_3D_TIFF_EXT=3d : 3次元画像をtiffにする PRJ_MODE=0 : -mの値
2軸で逆投影する場合は下記を変更します。
# mrc2Dto3D -I $*.3dlst -o $*.3d -InterpolationMode $(IN_MODE) -m $(PRJ_MODE) -single $(AXIS1) $(SIRT); mrc2Dto3D -I $*.3dlst -o $*.3d -InterpolationMode $(IN_MODE) -m $(PRJ_MODE) -Double $(SIRT);
make 3D make 3DTiff
方法 | 対応オプション | xy平面 | yz平面 | zx平面 | 全体 |
---|---|---|---|---|---|
単純逆投影 | -m 0 | ||||
フィルタ逆投影 | -m 1 | ||||
フィルタ逆投影(2軸) | -m 1 -Double |
電子線トモグラフィー
こちらのファイル(makefile込み)を入力ファイルとします。±60度で、10度刻みの画像からの3次元再構成を示しています。
方法 | 対応オプション | xy平面 | yz平面 | zx平面 | 全体 |
---|---|---|---|---|---|
単純逆投影 | -m 0 | ||||
フィルタ逆投影 | -m 1 | ||||
フィルタ逆投影(2軸) | -m 1 -Double |
2軸回転
単粒子解析
こちらのファイル(makefile込み)を入力ファイルとします。
方法 | 対応オプション | xy平面 | yz平面 | zx平面 | 全体 |
---|---|---|---|---|---|
単純逆投影 | -m 0 | ||||
フィルタ逆投影 | -m 1 | ||||
フィルタ逆投影(2軸) | -m 1 -Double |
電子線トモグラフィー
こちらのファイル(makefile込み)を入力ファイルとします。
方法 | 対応オプション | xy平面 | yz平面 | zx平面 | 全体 |
---|---|---|---|---|---|
単純逆投影 | -m 0 | ||||
フィルタ逆投影 | -m 1 | ||||
フィルタ逆投影(2軸) | -m 1 Double |
重み付き逆投影法
重み付き逆投影法による3次元再構成法をするためのプログラムです。
.3dlst.3d: mrc2Dto3D -I $(TARGET).3dlst -o $(TARGET).3d -InterpolationMode 2 -Double -DoubleCounter $(TARGET).3dcounter -CounterThreshold 0.5 -m 1 -WeightMode 2
$(TARGET).3dcounterが重みを表現しています。$(TARGET).3d及び$(TARGET).3dcounterがあれば、mrcImageAverageを使って、3次元画像を平均することができます。
1軸回転
単粒子解析
こちらのファイル(makefile込み)を入力ファイルとします。
IN_3D_TIFF_EXT=3d : 3次元画像をtiffにする PRJ_MODE=2 : -mの値
make 3D make 3DTiff
方法 | 対応オプション | xy平面 | yz平面 | zx平面 | 全体 |
---|---|---|---|---|---|
単純逆投影 | -m 0 | ||||
重み付き逆投影 | -m 2 |
電子線トモグラフィー
こちらのファイル(makefile込み)を入力ファイルとします。
方法 | 対応オプション | xy平面 | yz平面 | zx平面 | 全体 |
---|---|---|---|---|---|
単純逆投影 | -m 0 | ||||
重み付き逆投影 | -m 2 |
2軸回転
単粒子解析
こちらのファイル(makefile込み)を入力ファイルとします。
2軸で逆投影する場合は下記を変更します。
# mrc2Dto3D -I $*.3dlst -o $*.3d -InterpolationMode $(IN_MODE) -m $(PRJ_MODE) -single $(AXIS1) $(SIRT); mrc2Dto3D -I $*.3dlst -o $*.3d -InterpolationMode $(IN_MODE) -m $(PRJ_MODE) -Double $(SIRT);
方法 | 対応オプション | xy平面 | yz平面 | zx平面 | 全体 |
---|---|---|---|---|---|
単純逆投影 | -m 0 | ||||
重み付き逆投影 | -m 2 | ||||
重み付き逆投影(2軸) | -m 2 -Double |
電子線トモグラフィー
こちらのファイル(makefile込み)を入力ファイルとします。
方法 | 対応オプション | xy平面 | yz平面 | zx平面 | 全体 |
---|---|---|---|---|---|
単純逆投影 | -m 0 | ||||
重み付き逆投影 | -m 2 | ||||
重み付き逆投影(2軸) | -m 2 -Double |
SIRT法
1軸回転
単粒子解析
こちらのファイル(makefile込み)を入力ファイルとします。
#SIRT=-SIRT -max 30 -rms 0.0 : -maxで繰り返し回数、-SIRT無しのときはコメントアウト
逆投影法 / 繰り返し回数 | なし | 3 | 5 | 10 | 30 |
---|---|---|---|---|---|
単純逆投影(xy平面) | |||||
単純逆投影(yz平面) | |||||
単純逆投影(zx平面) | |||||
単純逆投影(全体) | |||||
フィルタ逆投影(xy平面) | |||||
フィルタ逆投影(yz平面) | |||||
フィルタ逆投影(zx平面) | |||||
フィルタ逆投影(全体) |
電子線トモグラフィー
こちらのファイル(makefile込み)を入力ファイルとします。
逆投影法 / 繰り返し回数 | なし | 3 | 5 | 10 | 30 |
---|---|---|---|---|---|
単純逆投影(xy平面) | |||||
単純逆投影(yz平面) | |||||
単純逆投影(zx平面) | |||||
単純逆投影(全体) | |||||
フィルタ逆投影(xy平面) | |||||
フィルタ逆投影(yz平面) | |||||
フィルタ逆投影(zx平面) | |||||
フィルタ逆投影(全体) |
2軸回転
単粒子解析
こちらのファイル(makefile込み)を入力ファイルとします。
逆投影法 / 繰り返し回数 | なし | 3 | 5 | 10 | 30 |
---|---|---|---|---|---|
単純逆投影(xy平面) | |||||
単純逆投影(yz平面) | |||||
単純逆投影(zx平面) | |||||
単純逆投影(全体) | |||||
フィルタ逆投影(xy平面) | |||||
フィルタ逆投影(yz平面) | |||||
フィルタ逆投影(zx平面) | |||||
フィルタ逆投影(全体) | |||||
重み付き逆投影(xy平面) | |||||
重み付き逆投影(yz平面) | |||||
重み付き逆投影(yz平面) | |||||
重み付き逆投影(全体) |
電子線トモグラフィー
こちらのファイル(makefile込み)を入力ファイルとします。
逆投影法 / 繰り返し回数 | なし | 3 | 5 | 10 | 30 |
---|---|---|---|---|---|
単純逆投影(xy平面) | |||||
単純逆投影(yz平面) | |||||
単純逆投影(zx平面) | |||||
単純逆投影(全体) | |||||
フィルタ逆投影(xy平面) | |||||
フィルタ逆投影(yz平面) | |||||
フィルタ逆投影(zx平面) | |||||
フィルタ逆投影(全体) | |||||
重み付き逆投影(xy平面) | |||||
重み付き逆投影(yz平面) | |||||
重み付き逆投影(yz平面) | |||||
重み付き逆投影(全体) |
ラドン変換による方法
ラドン変換を用いて3次元再構成を行うには以下のコマンドを使用します。
2次元画像のセット(複数枚)
↓ 各画像にmrcImageSinogramCreateにてシノグラム(2次元ラドン空間)に変換する
シノグラムのセット(複数枚)
↓ mrcRadon2Dto3Dにて3次元ラドンにする
3次元ラドンファイル
↓ mrcImageInverseRadonTransformにて逆ラドン変換を行う
3次元ファイル(再構成完了)
1軸回転
単粒子解析
こちらのファイル(makefile込み)を入力ファイルとします。
IN_3D_TIFF_EXT=3dr : 3次元画像をtiffにする RBP_MODE=1 : -mの値 #R_MAX=-Rmax 0.1 : ハイパスフィルタ閾値 #R_MIN=-Rmin 0.05 : ローパスフィルタ閾値 #THRES=-T 10 : m=5のみ使用
make Radon3D make 3DTiff
電子線トモグラフィー
こちらのファイル(makefile込み)を入力ファイルとします。
2軸回転
単粒子解析
こちらのファイル(makefile込み)を入力ファイルとします。
電子線トモグラフィー
こちらのファイル(makefile込み)を入力ファイルとします。