「3次元再構成・逆投影」の版間の差分
(→重み付き逆投影法) |
(→ラドン変換による方法) |
||
(2人の利用者による、間の23版が非表示) | |||
行26: | 行26: | ||
== 3次元再構成の実際 == | == 3次元再構成の実際 == | ||
+ | ===フィルタ逆投影法=== | ||
+ | ====1軸回転==== | ||
+ | =====単粒子解析===== | ||
+ | [[Media:Input-BP-1.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <div>[[mrc2Dto3D]]はオプションm=1でフィルタ逆投影に設定できます。</div> | ||
+ | <div>makefileでの設定項目</div> | ||
+ | <pre> | ||
+ | IN_3D_TIFF_EXT=3d : 3次元画像をtiffにする | ||
+ | PRJ_MODE=0 : -mの値 | ||
+ | </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> | ||
+ | <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 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-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-2-2.png]]</td> | ||
+ | <td>[[ファイル:Outdata-BP-2.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <div>フィルタ逆投影の方がくっきりとした画像が得られます。</div> | ||
+ | <br> | ||
+ | |||
+ | =====電子線トモグラフィー===== | ||
+ | [[Media:Input-BP-3.zip|こちらのファイル]](makefile込み)を入力ファイルとします。±60度で、10度刻みの画像からの3次元再構成を示しています。 | ||
+ | <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 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-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-1.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | ====2軸回転==== | ||
+ | =====単粒子解析===== | ||
+ | [[Media:Input-BP-2.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-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 1</td> | ||
+ | <td>[[画像:Outdata-2BP-2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-2-2.png]]</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> | ||
+ | </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 1</td> | ||
+ | <td>[[画像:Outdata-T2BP-1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-1-2.png]]</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次元再構成法をするためのプログラムです。 | 重み付き逆投影法による3次元再構成法をするためのプログラムです。 | ||
行36: | 行220: | ||
$(TARGET).3dcounterが重みを表現しています。$(TARGET).3d及び$(TARGET).3dcounterがあれば、[[mrcImageAverage]]を使って、3次元画像を平均することができます。 | $(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> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | ===SIRT法=== | ||
+ | ====1軸回転==== | ||
+ | =====単粒子解析===== | ||
+ | [[Media:Input-BP-1.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <div>[[mrc2Dto3D]]はオプション-SIRTなどで設定できます。</div> | ||
+ | <div>makefileでの設定項目</div> | ||
+ | <pre> | ||
+ | #SIRT=-SIRT -max 30 -rms 0.0 : -maxで繰り返し回数、-SIRT無しのときはコメントアウト | ||
+ | </pre> | ||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>逆投影法 / 繰り返し回数</th> | ||
+ | <td>なし</td> | ||
+ | <td>3</td> | ||
+ | <td>5</td> | ||
+ | <td>10</td> | ||
+ | <td>30</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-BP-0-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-4-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-6-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-8-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-BP-0-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-4-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-6-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-8-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(zx平面)</td> | ||
+ | <td>[[画像:Outdata-BP-0-2.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-4-2.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-6-2.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-8-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-BP-0.gif]]</td> | ||
+ | <td>[[画像:Outdata-BP-1.gif]]</td> | ||
+ | <td>[[画像:Outdata-BP-4.gif]]</td> | ||
+ | <td>[[画像:Outdata-BP-6.gif]]</td> | ||
+ | <td>[[画像:Outdata-BP-8.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-BP-2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-7-0.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-9-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-BP-2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-7-1.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-9-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(zx平面)</td> | ||
+ | <td>[[画像:Outdata-BP-2-2.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-7-2.png]]</td> | ||
+ | <td>[[画像:Outdata-BP-9-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-BP-2.gif]]</td> | ||
+ | <td>[[画像:Outdata-BP-3.gif]]</td> | ||
+ | <td>[[画像:Outdata-BP-5.gif]]</td> | ||
+ | <td>[[画像:Outdata-BP-7.gif]]</td> | ||
+ | <td>[[画像:Outdata-BP-9.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | =====電子線トモグラフィー===== | ||
+ | [[Media:Input-BP-3.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>逆投影法 / 繰り返し回数</th> | ||
+ | <td>なし</td> | ||
+ | <td>3</td> | ||
+ | <td>5</td> | ||
+ | <td>10</td> | ||
+ | <td>30</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-TBP-0-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-4-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-6-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-8-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-TBP-0-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-4-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-6-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-8-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(zx平面)</td> | ||
+ | <td>[[画像:Outdata-TBP-0-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-2-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-4-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-6-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-8-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-TBP-0.gif]]</td> | ||
+ | <td>[[画像:Outdata-TBP-2.gif]]</td> | ||
+ | <td>[[画像:Outdata-TBP-4.gif]]</td> | ||
+ | <td>[[画像:Outdata-TBP-6.gif]]</td> | ||
+ | <td>[[画像:Outdata-TBP-8.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-TBP-1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-7-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-9-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-TBP-1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-7-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-9-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(zx平面)</td> | ||
+ | <td>[[画像:Outdata-TBP-1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-7-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TBP-9-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-TBP-1.gif]]</td> | ||
+ | <td>[[画像:Outdata-TBP-3.gif]]</td> | ||
+ | <td>[[画像:Outdata-TBP-5.gif]]</td> | ||
+ | <td>[[画像:Outdata-TBP-7.gif]]</td> | ||
+ | <td>[[画像:Outdata-TBP-9.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <div>回数を増やすと鮮明になりますが、yz平面はブレたままです。</div> | ||
+ | <br> | ||
+ | |||
+ | ====2軸回転==== | ||
+ | =====単粒子解析===== | ||
+ | [[Media:Input-BP-2.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>逆投影法 / 繰り返し回数</th> | ||
+ | <td>なし</td> | ||
+ | <td>3</td> | ||
+ | <td>5</td> | ||
+ | <td>10</td> | ||
+ | <td>30</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-2BP-0-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-4-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-6-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-8-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-2BP-0-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-4-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-6-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-8-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(zx平面)</td> | ||
+ | <td>[[画像:Outdata-2BP-0-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-4-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-6-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-8-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-2BP-0.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-1.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-4.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-6.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-8.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-2BP-d2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d7-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d9-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-2BP-d2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d7-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d9-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(zx平面)</td> | ||
+ | <td>[[画像:Outdata-2BP-d2-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d7-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d9-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-2BP-d2.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d3.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d5.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d7.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-d9.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-2BP-wd-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd10-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd30-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-2BP-wd-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd10-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd30-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-2BP-wd-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd10-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd30-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-2BP-wd.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd3.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd5.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd10.gif]]</td> | ||
+ | <td>[[画像:Outdata-2BP-wd30.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | =====電子線トモグラフィー===== | ||
+ | [[Media:Input-BP-4.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>逆投影法 / 繰り返し回数</th> | ||
+ | <td>なし</td> | ||
+ | <td>3</td> | ||
+ | <td>5</td> | ||
+ | <td>10</td> | ||
+ | <td>30</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-T2BP-0-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-4-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-6-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-8-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-T2BP-0-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-4-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-6-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-8-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(zx平面)</td> | ||
+ | <td>[[画像:Outdata-T2BP-0-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-2-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-4-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-6-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-8-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-T2BP-0.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-2.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-4.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-6.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-8.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-T2BP-d1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d7-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d9-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-T2BP-d1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d7-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d9-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(zx平面)</td> | ||
+ | <td>[[画像:Outdata-T2BP-d1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d7-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d9-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>フィルタ逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-T2BP-d1.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d3.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d5.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d7.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-d9.gif]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(xy平面)</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd10-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd30-0.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd10-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd30-1.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(yz平面)</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd10-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd30-2.png]]</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(全体)</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd3.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd5.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd10.gif]]</td> | ||
+ | <td>[[画像:Outdata-T2BP-wd30.gif]]</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | ===[[ラドン変換]]による方法=== | ||
+ | [[ラドン変換]]を用いて[[3次元再構成]]を行うには以下の[[コマンド]]を使用します。<br> | ||
<br> | <br> | ||
2次元画像のセット(複数枚) | 2次元画像のセット(複数枚) | ||
− | ↓ 各画像に[[mrcImageSinogramCreate]] | + | ↓ 各画像に[[mrcImageSinogramCreate]]にて[[シノグラム]](2次元ラドン空間)に変換する |
− | + | [[シノグラム]]のセット(複数枚) | |
↓ [[mrcRadon2Dto3D]]にて3次元ラドンにする | ↓ [[mrcRadon2Dto3D]]にて3次元ラドンにする | ||
行49: | 行788: | ||
3次元ラドンファイル | 3次元ラドンファイル | ||
− | ↓ [[mrcImageInverseRadonTransform]] | + | ↓ [[mrcImageInverseRadonTransform]]にて逆ラドン変換を行う |
3次元ファイル(再構成完了) | 3次元ファイル(再構成完了) | ||
+ | |||
+ | |||
+ | ====1軸回転==== | ||
+ | =====単粒子解析===== | ||
+ | [[Media:Input-BP-1.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <div>[[mrcImageInverseRadonTransform]]は-mなどで逆投影方法を設定できます。</div> | ||
+ | <div>makefileでの設定項目</div> | ||
+ | <pre> | ||
+ | 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のみ使用 | ||
+ | </pre> | ||
+ | <div>makefileでのコマンド</div> | ||
+ | <pre> | ||
+ | make Radon3D | ||
+ | make 3DTiff | ||
+ | </pre> | ||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>方法</th> | ||
+ | <th>対応オプション</th> | ||
+ | <th>xy平面</th> | ||
+ | <th>yz平面</th> | ||
+ | <th>zx平面</th> | ||
+ | <th>全体</th> | ||
+ | <th>補足</th> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影</td> | ||
+ | <td>-m 1</td> | ||
+ | <td>[[画像:Outdata-RBP-1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-1.gif]]</td> | ||
+ | <td></td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影</td> | ||
+ | <td>-m 2</td> | ||
+ | <td>[[画像:Outdata-RBP-2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-2-2.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-2.gif]]</td> | ||
+ | <td></td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(フィルタあり)</td> | ||
+ | <td>-m 3</td> | ||
+ | <td>[[画像:Outdata-RBP-3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-3.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>NearestCountRWeightedBackProjection</td> | ||
+ | <td>-m 4</td> | ||
+ | <td>[[画像:Outdata-RBP-4-0.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-4-1.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-4-2.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-4.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>WeightedCountRWeightedBackProjection</td> | ||
+ | <td>-m 5</td> | ||
+ | <td>[[画像:Outdata-RBP-5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-5.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1 -T 1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(フィルタあり)</td> | ||
+ | <td>-m 6</td> | ||
+ | <td>[[画像:Outdata-RBP-6-0.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-6-1.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-6-2.png]]</td> | ||
+ | <td>[[画像:Outdata-RBP-6.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</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> | ||
+ | <th>補足</th> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影</td> | ||
+ | <td>-m 1</td> | ||
+ | <td>[[画像:Outdata-TRBP-1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-1.gif]]</td> | ||
+ | <td></td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影</td> | ||
+ | <td>-m 2</td> | ||
+ | <td>[[画像:Outdata-TRBP-2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-2-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-2.gif]]</td> | ||
+ | <td></td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(フィルタあり)</td> | ||
+ | <td>-m 3</td> | ||
+ | <td>[[画像:Outdata-TRBP-3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-3.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>NearestCountRWeightedBackProjection</td> | ||
+ | <td>-m 4</td> | ||
+ | <td>[[画像:Outdata-TRBP-4-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-4-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-4-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-4.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>WeightedCountRWeightedBackProjection</td> | ||
+ | <td>-m 5</td> | ||
+ | <td>[[画像:Outdata-TRBP-5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-5.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1 -T 1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(フィルタあり)</td> | ||
+ | <td>-m 6</td> | ||
+ | <td>[[画像:Outdata-TRBP-6-0.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-6-1.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-6-2.png]]</td> | ||
+ | <td>[[画像:Outdata-TRBP-6.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> | ||
+ | |||
+ | ====2軸回転==== | ||
+ | =====単粒子解析===== | ||
+ | [[Media:Input-BP-2.zip|こちらのファイル]](makefile込み)を入力ファイルとします。 | ||
+ | <br> | ||
+ | |||
+ | <table border="1"> | ||
+ | <tr> | ||
+ | <th>方法</th> | ||
+ | <th>対応オプション</th> | ||
+ | <th>xy平面</th> | ||
+ | <th>yz平面</th> | ||
+ | <th>zx平面</th> | ||
+ | <th>全体</th> | ||
+ | <th>補足</th> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影</td> | ||
+ | <td>-m 1</td> | ||
+ | <td>[[画像:Outdata-2RBP-1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-1.gif]]</td> | ||
+ | <td></td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影</td> | ||
+ | <td>-m 2</td> | ||
+ | <td>[[画像:Outdata-2RBP-2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-2-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-2.gif]]</td> | ||
+ | <td></td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(フィルタあり)</td> | ||
+ | <td>-m 3</td> | ||
+ | <td>[[画像:Outdata-2RBP-3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-3.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>NearestCountRWeightedBackProjection</td> | ||
+ | <td>-m 4</td> | ||
+ | <td>[[画像:Outdata-2RBP-4-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-4-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-4-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-4.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>WeightedCountRWeightedBackProjection</td> | ||
+ | <td>-m 5</td> | ||
+ | <td>[[画像:Outdata-2RBP-5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-5.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1 -T 1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(フィルタあり)</td> | ||
+ | <td>-m 6</td> | ||
+ | <td>[[画像:Outdata-2RBP-6-0.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-6-1.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-6-2.png]]</td> | ||
+ | <td>[[画像:Outdata-2RBP-6.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</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> | ||
+ | <th>補足</th> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影</td> | ||
+ | <td>-m 1</td> | ||
+ | <td>[[画像:Outdata-T2RBP-1-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-1-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-1-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-1.gif]]</td> | ||
+ | <td></td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影</td> | ||
+ | <td>-m 2</td> | ||
+ | <td>[[画像:Outdata-T2RBP-2-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-2-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-2-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-2.gif]]</td> | ||
+ | <td></td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>重み付き逆投影(フィルタあり)</td> | ||
+ | <td>-m 3</td> | ||
+ | <td>[[画像:Outdata-T2RBP-3-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-3-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-3-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-3.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>NearestCountRWeightedBackProjection</td> | ||
+ | <td>-m 4</td> | ||
+ | <td>[[画像:Outdata-T2RBP-4-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-4-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-4-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-4.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>WeightedCountRWeightedBackProjection</td> | ||
+ | <td>-m 5</td> | ||
+ | <td>[[画像:Outdata-T2RBP-5-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-5-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-5-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-5.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1 -T 1</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td>単純逆投影(フィルタあり)</td> | ||
+ | <td>-m 6</td> | ||
+ | <td>[[画像:Outdata-T2RBP-6-0.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-6-1.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-6-2.png]]</td> | ||
+ | <td>[[画像:Outdata-T2RBP-6.gif]]</td> | ||
+ | <td>-Rmin 0.05 -Rmax 0.1</td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | <br> |
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込み)を入力ファイルとします。