/usr/share/pyshared/fabio/fit2dspreadsheetimage.py is in python-fabio 0.0.8-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 | """
Read the fit2d ascii image output
+ Jon Wright, ESRF
"""
import numpy
from fabioimage import fabioimage
class fit2dspreadsheetimage(fabioimage):
"""
Read a fit2d ascii format
"""
def _readheader(self, infile):
"""
TODO : test for minimal attributes?
"""
line = infile.readline()
try:
items = line.split()
xdim = int(items[0])
ydim = int(items[1])
except:
raise
self.header['title'] = line
self.header['Dim_1'] = xdim
self.header['Dim_2'] = ydim
def read(self, fname):
"""
Read in header into self.header and
the data into self.data
"""
self.header = {}
self.resetvals()
infile = self._open(fname)
self._readheader(infile)
# Compute image size
try:
self.dim1 = int(self.header['Dim_1'])
self.dim2 = int(self.header['Dim_2'])
except:
raise Exception("file", str(fname) + \
"is corrupt, cannot read it")
bytecode = numpy.float32
self.bpp = len(numpy.array(0, bytecode).tostring())
#now read the data into the array
try:
vals = []
for line in infile.readlines():
try:
vals.append([float(x) for x in line.split()])
except:
pass
self.data = numpy.array(vals).astype(bytecode)
assert self.data.shape == (self.dim2, self.dim1)
except:
raise IOError, "Error reading ascii"
self.resetvals()
# ensure the PIL image is reset
self.pilimage = None
return self
if __name__ == "__main__":
import sys, time
start = time.time()
img = fit2dspreadsheetimage()
img.read(sys.argv[1])
print time.time() - start
print img.dim1, img.dim2, img.data.shape
from matplotlib.pylab import imshow, show
imshow(img.data.T)
show()
|