Project

General

Profile

Actions

Bug #3406

open

CTF display crashes on close to focus images

Added by Gabriel Lander about 9 years ago. Updated about 9 years ago.

Status:
Assigned
Priority:
Normal
Assignee:
Category:
Image Processing
Target version:
-
Start date:
08/21/2015
Due date:
% Done:

0%

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

Description

We're running CTFFind on negative stain data that was collected very close to focus (~0.4 microns), and we get the attached error. We don't have this problem if we run myami/3.0

Beginning Main Loop

Starting image 1 ( skip:344, remain:67 ) id:3719578, file: 15aug17j_033sq_07hl_04en
... Pixel size: 2.05044
!!! WARNING: no best CTF value for image 15aug17j_033sq_07hl_04en
... Xmag=76081, dstep=1.56e-05, mpix=2.05e-10
!!! WARNING: Defocus minimum is less than zero
Defocus search range: 1000 A to 30000 A (0.10 to 3.00 um)
... running ctf estimation at Fri Aug 21 16:23:46 2015
/opt/applications/frealign/8.09/bin/ctffind64.exe
15aug17j_033sq_07hl_04en.mrc
15aug17j_033sq_07hl_04en-pow.mrc
2.2,120.0,0.15,76081.2313455,15.6,2
512,80.0,6.0,1000.0,30000.0,1000.0,500.0
... ctf estimation completed in 33.9 sec
... nominal=5.0e-07, bestdef=5.0e-07,
... def_1=6.1e-07, def_2=5.8e-07, astig_angle=-7.6, cross_corr=0.273,

... reading MRC: 15aug17j_033sq_07hl_04en-pow.mrc size:(512, 512) dtype:float32
... writing JPEG: 15aug17j_033sq_07hl_04en-pow.jpg ==== Committing data to database ====
... Committing ctf parameters for 15aug17j_033sq_07hl_04en to database
!!! WARNING: |def1| > |def2|, flipping defocus axes
... Reading image...
[CTF param] def1: 5.81e-07 | def2: 6.10e-07 | angle: 82.4 | ampcontr 0.15 | defratio 1.050
Final params: def1: 5.81e-07 | def2: 6.10e-07 | angle: 82.4 | defratio 1.05
... Defocus Astig Percent Diff 2.43 -- 5.809e-07, 6.099e-07
/opt/applications/python/lib/python2.6/site-packages/numpy-1.9.1-py2.6-linux-x86_64.egg/numpy/core/numeric.py:1846: DeprecationWarning: using a non-integer number instead of an integer wi
ll result in an error in the future
res = empty((N,)+dimensions, dtype=dtype)
/opt/applications/python/lib/python2.6/site-packages/numpy-1.9.1-py2.6-linux-x86_64.egg/numpy/core/numeric.py:1849: DeprecationWarning: using a non-integer number instead of an integer wi
ll result in an error in the future
tmp.shape = (1,)*i + (dim,)+(1,)*(N-i-1)
/opt/applications/python/lib/python2.6/site-packages/numpy-1.9.1-py2.6-linux-x86_64.egg/numpy/core/numeric.py:1851: DeprecationWarning: using a non-integer number instead of an integer wi
ll result in an error in the future
val = zeros(newdim, dtype)
Computing power spectra in 2048x2048 blocks./opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctfpower.py:101: DeprecationWarning: using a non-integer number inste
ad of an integer will result in an error in the future
cutout = image[x1:x2, y1:y2]
........
... Computing median of power spectra series
... Compute PSD with fieldsize 2048 and 9 images complete in 22.18 sec
!!! WARNING: Requested resolution (3.000) is not available (4.101)
... Preform a rotational average and remove spikes...
PART 1: SETUP PARAMETERS AND ELLIPTICAL AVERAGE
... Number of available peaks is 24
/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctftools.py:490: DeprecationWarning: using a non-integer number instead of an integer will result in an error in
the future
xdataint = xdataint[:outercutsize]
... Determine and subtract noise model
PART 2: BACKGROUND NOISE SUBTRACTION
... starting peak extension
Peak Extension Complete in 0.88 msec

Found 0 ctf values
Traceback (most recent call last):
File "/opt/applications/myami/trunk/bin/ctfestimate.py", line 445, in <module>
imgLoop.run()
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/appionLoop2.py", line 93, in run
self.loopCommitToDatabase(imgdata)
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/appionLoop2.py", line 143, in loopCommitToDatabase
return self.commitToDatabase(imgdata)
File "/opt/applications/myami/trunk/bin/ctfestimate.py", line 345, in commitToDatabase
ctfinsert.validateAndInsertCTFData(imgdata, self.ctfvalues, self.ctfrun, self.params['rundir'])
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctfinsert.py", line 45, in validateAndInsertCTFData
ctfvalues = runCTFdisplayTools(imgdata, ctfvalues, opimagedir, fftpath, fftfreq)
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctfinsert.py", line 83, in runCTFdisplayTools
ctfdisplaydict = ctfdisplay.makeCtfImages(imgdata, ctfvalues, fftpath, fftfreq)
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctfdisplay.py", line 1257, in makeCtfImages
ctfdisplaydict = a.CTFpowerspec(imgdata, ctfdata, fftpath, fftfreq, twod=twod)
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctfdisplay.py", line 1139, in CTFpowerspec
normpowerspec = self.normalizeCtf(powerspec, twod=twod)
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctfdisplay.py", line 210, in normalizeCtf
fitvalleydata[startIndex:endIndex], "below")
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctfnoise.py", line 550, in modelCTFNoise
filterctfdata, contraintFunction)
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctfnoise.py", line 399, in getAllInitialParameters
fitparams, value = self.fitLinear(xdata, ctfdata, contraintFunction)
File "/opt/applications/myami/trunk/lib/python2.6/site-packages/appionlib/apCtf/ctfnoise.py", line 183, in fitLinear
z = numpy.polyfit(xdata, ctfdata, 1)
File "/opt/applications/python/lib/python2.6/site-packages/numpy-1.9.1-py2.6-linux-x86_64.egg/numpy/lib/polynomial.py", line 552, in polyfit
raise TypeError("expected non-empty vector for x")
TypeError: expected non-empty vector for x

Actions

Also available in: Atom PDF