/usr/lib/python2.7/dist-packages/fabio/pixiimage.py is in python-fabio 0.3.0+dfsg-1build1.
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 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 | # coding: utf-8
#
# Project: X-ray image reader
# https://github.com/kif/fabio
#
#
# Copyright (C) European Synchrotron Radiation Facility, Grenoble, France
#
# Principal author: Jérôme Kieffer (Jerome.Kieffer@ESRF.eu)
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
#
"""
Author: Jon Wright, ESRF.
"""
# Get ready for python3:
from __future__ import with_statement, print_function
__authors__ = ["Jon Wright", "Jérôme Kieffer"]
__contact__ = "wright@esrf.fr"
__license__ = "GPLv3+"
__copyright__ = "European Synchrotron Radiation Facility, Grenoble, France"
__date__ = "30/10/2015"
import numpy
import sys
import os
from .fabioimage import FabioImage
from .fabioutils import previous_filename, next_filename
class PixiImage(FabioImage):
_need_a_seek_to_read = True
def _readheader(self, infile):
infile.seek(0)
self.header = self.check_header()
byt = infile.read(4)
framesize = numpy.fromstring(byt, numpy.int32)
if framesize == 243722:
# life is good
width = 476
height = 512
offset = 24
self.header['framesize'] = framesize
self.header['width'] = width
self.header['height'] = height
self.header['offset'] = offset
else:
print("Pixiimage, bad framesize: %s" % framesize)
raise
def read(self, fname, frame=None):
if frame is None:
frame = 0
self.header = self.check_header()
self.resetvals()
infile = self._open(fname, "rb")
self.sequencefilename = fname
self._readheader(infile)
self.nframes = os.path.getsize(fname) / 487448
self._readframe(infile, frame)
infile.close()
return self
def _makeframename(self):
self.filename = "%s$%04d" % (self.sequencefilename,
self.currentframe)
def _readframe(self, filepointer, img_num):
if (img_num > self.nframes or img_num < 0):
raise Exception("Bad image number")
imgstart = self.header['offset'] + img_num * (512 * 476 * 2 + 24)
filepointer.seek(imgstart, 0)
self.data = numpy.fromstring(filepointer.read(512 * 476 * 2),
numpy.uint16)
self.data.shape = self.header['height'], self.header['width']
self.dim2, self.dim1 = self.data.shape
self.currentframe = int(img_num)
self._makeframename()
def write(self, fname, force_type=numpy.uint16):
raise Exception("Write is not implemented")
def getframe(self, num):
"""
Returns a frame as a new FabioImage object
"""
if num < 0 or num > self.nframes:
raise Exception("Requested frame number is out of range")
# Do a deep copy of the header to make a new one
newheader = {}
for k in self.header.keys():
newheader[k] = self.header[k]
frame = pixiimage(header=newheader)
frame.nframes = self.nframes
frame.sequencefilename = self.sequencefilename
infile = frame._open(self.sequencefilename, "rb")
frame._readframe(infile, num)
infile.close()
return frame
def next(self):
"""
Get the next image in a series as a fabio image
"""
if self.currentframe < (self.nframes - 1) and self.nframes > 1:
return self.getframe(self.currentframe + 1)
else:
newobj = pixiimage()
newobj.read(next_filename(
self.sequencefilename))
return newobj
def previous(self):
"""
Get the previous image in a series as a fabio image
"""
if self.currentframe > 0:
return self.getframe(self.currentframe - 1)
else:
newobj = pixiimage()
newobj.read(previous_filename(
self.sequencefilename))
return newobj
pixiimage = PixiImage
def demo(fname):
i = PixiImage()
i.read(fname)
import pylab
pylab.imshow(numpy.log(i.data))
print("%s\t%s\t%s\t%s" % (i.filename, i.data.max(), i.data.min(), i.data.mean()))
pylab.title(i.filename)
pylab.show()
while 1:
i = i.next()
pylab.imshow(numpy.log(i.data))
pylab.title(i.filename)
pylab.show()
raw_input()
if __name__ == "__main__":
import sys
demo(sys.argv[1])
|