Commit c1e6c18b authored by Samuel  HANOT's avatar Samuel HANOT

use 0.5 FSC criterion to get the resolution

parent d8a02c9f
......@@ -2,7 +2,7 @@ set terminal pdfcairo size 8.5cm, 5.25cm font 'helvetica,12' enhanced
unset key
set output 'plots/fscs.pdf'
p for [f in system('ls -v [0-9]*/*threshold.fsc')] f t f, 0.143
p for [f in system('ls -v [0-9]*/*threshold.fsc')] f t f, 0.5
set output 'plots/resolution.pdf'
resolution=system('cat resolution.txt')
set logscale x 10
......
......@@ -28,8 +28,8 @@ ${bindir}/gmconvert/gmconvert -igmm ${n}/${n}.gmm -imap threshold.map -omap ${n}
e2proc3d.py ${map_name} ${n}/${n}.gmm.mrc.fsc --calcfsc ${n}/${n}.gmm.mrc
e2proc3d.py threshold.map ${n}/${n}.gmm.mrc.threshold.fsc --calcfsc ${n}/${n}.gmm.mrc
e2proc3d.py ${n}/${n}.gmm.mrc tmp.map --calcsf ${n}/${n}.gmm.mrc.sf
awk '$2<=0.143 && y>0.143{print 1.0/(x+(0.143-y)*($1-x)/($2-y))} {x=$1;y=$2}' < ${n}/${n}.gmm.mrc.fsc > ${n}/${n}.gmm.resolution
awk '$2<=0.143 && y>0.143{print 1.0/(x+(0.143-y)*($1-x)/($2-y))} {x=$1;y=$2}' < ${n}/${n}.gmm.mrc.threshold.fsc > ${n}/${n}.gmm.threshold.resolution
awk '$2<=0.5 && y>0.5{print 1.0/(x+(0.5-y)*($1-x)/($2-y))} {x=$1;y=$2}' < ${n}/${n}.gmm.mrc.fsc > ${n}/${n}.gmm.resolution
awk '$2<=0.5 && y>0.5{print 1.0/(x+(0.5-y)*($1-x)/($2-y))} {x=$1;y=$2}' < ${n}/${n}.gmm.mrc.threshold.fsc > ${n}/${n}.gmm.threshold.resolution
rm -f $n/sub*
rm -rf {.,$n}/tmp*
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment