Project

General

Profile

Actions

Bug #2053

closed

ACE2 not accurate

Added by Gabriel Lander about 12 years ago. Updated about 10 years ago.

Status:
Closed
Priority:
Normal
Category:
-
Target version:
-
Start date:
09/18/2012
Due date:
% Done:

0%

Estimated time:
Affected Version:
Appion/Leginon 3.0.0
Show in known bugs:
No
Workaround:

Description

We finally got redux running and can look at the plots from ACE2, and we see surprisingly poor estimations. I am attaching a two 1D plots for comparison of results from ACE2 and from CTFfind, and it is clear that ACE2 is pretty far off. These plots are for what I would consider a good image with Thon rings out to nearly 8 Angstroms.
Additionally, we've noticed that vast majority of images output "!!! WARNING: Image is possibly over-focused" when the 1/30A - 1/10A and 5 peak confidences are being calculated.
Any ideas why we might be having these issues?
running myami r17119

Here's the output for the ACE2 run for this image, but I can send the MRC file if you need it for debugging:

Starting image 1 ( skip:0, remain:0 ) id:696246, file: 12sep10a_01gr_082sq_02hl_02en ... Pixel size: 1.15450433333
... Ace2 executable: /archive/glander/myami/appion/bin/ace2.exe
... running ace2 at Thu Sep 13 08:08:13 2012
/archive/glander/myami/appion/bin/ace2.exe -i /leginon4/glander/12sep10a/rawdata/12sep10a_00001gr_00082sq_v02_00002hl_00002en.mrc -b 2 -c 2.2 -k 120.0 -a 1.15450433333 -e 10.0,0.001 -r 0.0

Program: /archive/glander/myami/appion/bin/ace2.exe in source ace2.m, svn revision 17012,
which was compiled on Aug 9 2012 at 10:54:00.

Reading Image /leginon4/glander/12sep10a/rawdata/12sep10a_00001gr_00082sq_v02_00002hl_00002en.mrc: DONE in 0.12 seconds
Processing image with 4096 4096 pixels
APIX: 1.154504 KV: 120.000000 CS: 2.200000
Binning image by 2 pixels... DONE in 0.05 seconds
Blurring 100 edge pixels... DONE in 0.02 seconds
Generating power spectrum... DONE in 4.31 seconds
post-binning by 2...
Blurring edges by 10.0 pixels...
Finding edges for ellipse fitting... DONE in 0.25 seconds (12495 edges)
Using RANSAC to find ellipse parameters...
RANSAC SUCCESS, 608/12495 (4.87%-99.00%)
DONE in 1.23 seconds
--
... Committing ctf parameters for 12sep10a_01gr_082sq_02hl_02en to database
... Reading image...
[CTF param] def1: 2.42e-06 | def2: 2.47e-06 | angle: 35.2 | ampcontr 0.27 | defratio 1.021
[CTF stats] conf: 0.978 | conf_d: 0.978
Final params: def1: 2.42e-06 | def2: 2.47e-06 | angle: 35.2 | defratio 1.02
... Defocus Astig Percent Diff 1.02 -- 2.419e-06, 2.469e-06
Computing power spectra.........
... Compute PSD with fieldsize 2048 and 9 images complete in 10.34 sec
... Preform a rotational average and remove spikes...
PART 1: SETUP PARAMETERS AND ELLIPTICAL AVERAGE
... Number of available peaks is 250
... Determine and subtract noise model
PART 2: BACKGROUND NOISE SUBTRACTION
... starting peak extension
Peak Extension Complete in 0.98 msec
... Constraint violation: 0.014 < -0.001
... Constraint violation: -0.009 < -0.001
... Constraint violation: -0.028 < -0.001
... Constraint violation: -0.006 < -0.001
Noise Model Complete in 86.04 msec
... Constraint violation: -0.005 < -0.001
-

... Saving 1D graph to file 12sep10a_01gr_082sq_02hl_02en-plots.png
PART 8: NORMALIZE THE 2D IMAGE
PART 9: DRAW THE 2D POWERSPEC IMAGE
... Saving 2D powerspectra to file: 12sep10a_01gr_082sq_02hl_02en-powerspec.jpg
[CTF run] method: ace2 | runname acetwo1
[CTF param] def1: 2.42e-06 | def2: 2.47e-06 | angle: 35.2 | ampcontr 0.27 | defratio 1.021
[CTF stats] conf_30-10: 0.363 | conf_5peak: 0.509 | res_0.8: 100.0A | res_0.5 100.0A
[CTF stats] conf: 0.978 | conf_d: 0.978

SUMMARY: pyace2
------------------------------------------
PEAKS: 0 peaks of 0
AVG PEAKS: 0.0 +/- 0.0 peaks
(- ESTIMATE: 0 to 0 total peaks )
-----------------------------------------

TIME: 24.83 sec
AVG TIME: 26.67 +/- 32.91 sec
------------------------------------------

Files

acetwo.png (474 KB) acetwo.png Gabriel Lander, 09/18/2012 12:45 PM
ctffind.png (582 KB) ctffind.png Gabriel Lander, 09/18/2012 12:45 PM

Related issues 2 (1 open1 closed)

Related to Appion - Bug #2016: How to get CTF data from databaseClosedNeil Voss08/23/2012

Actions
Related to Appion - Bug #2068: Resolve CTF problem effort by NRAMM local teamNew10/01/2012

Actions
Actions

Also available in: Atom PDF