Bug #1470
closedFrealign prep dtype float32 error
0%
Description
Hey Jim,
Niko and I were trying to do a frealign job today and came across a bug with prepRefineFrealign. Here is the error we are getting:
Traceback (most recent call last): File "/home/lfisher/myami/appion/bin/prepRefineFrealign.py", line 192, in ? app.start() File "/home/lfisher/myami/appion/appionlib/apPrepRefine.py", line 336, in start self.__processStack(prepdata) File "/home/lfisher/myami/appion/appionlib/apPrepRefine.py", line 297, in __processStack self.convertToRefineStack() File "/home/lfisher/myami/appion/bin/prepRefineFrealign.py", line 110, in convertToRefineStack self.ImagicStackToFrealignMrcStack() File "/home/lfisher/myami/appion/bin/prepRefineFrealign.py", line 94, in ImagicStackToFrealignMrcStack apIMAGIC.convertImagicStackToMrcStack(stackroot,stackbaseroot+'.mrc') File "/home/lfisher/myami/appion/appionlib/apIMAGIC.py", line 410, in convertImagicStackToMrcStack imagic2mrc.imagic_to_mrc(infile,outfile) File "/home/lfisher/myami/pyami/imagic2mrc.py", line 22, in imagic_to_mrc print pyami.mrc.numpy2mrc[mrc_header['dtype']] KeyError: dtype('float32')
Using Python 2.4.3, imagic.py code_to_dtype is setting dtype to float32 instead of numpy.float32, which causes the keyerror in imagic2mrc.py.
Anchi provided assistance and when we changed the code from:
code_to_dtype = { numpy.frombuffer('REAL', numpy.int32)[0]: numpy.dtype(numpy.float32), # ... }
to:
code_to_dtype = { numpy.frombuffer('REAL', numpy.int32)[0]: numpy.float32, # ... }
it solved the issue.
Updated by Lauren Fisher about 13 years ago
- Category set to Python scripting
- Assignee set to Jim Pulokas
Updated by Anchi Cheng almost 13 years ago
Got the same error. Somehow imagic_header is not numpy. As a workaround, I commented out the offending line in imagic2mrc.py in my sandbox and it worked.
Jim. please still investigate.
Updated by Anchi Cheng almost 13 years ago
- Target version set to Appion/Leginon 2.2.0
- Affected Version changed from Appion/Leginon 2.1.0 to Appion/Leginon 2.2.0 (trunk)
Updated by Anchi Cheng almost 13 years ago
- Status changed from New to In Test
- Assignee changed from Jim Pulokas to Lauren Fisher
Jim fixed the bug (my workaround didn't work) in r16262
Updated by Amber Herold over 12 years ago
Can this be closed? It is a very small low risk change and the code has been in for 3 months and it sounds like you may have run Jims fix Anchi?
Updated by Anchi Cheng over 12 years ago
I did test it, when I updated the issue. Just wanted to make sure the reporter has a chance to test.
Updated by Lauren Fisher over 12 years ago
- Status changed from In Test to Closed