# coding: utf-8
#
# Project: X-ray image reader
# https://github.com/silx-kit/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 .
#
"""
Authors:
........
* Henning O. Sorensen & Erik Knudsen:
Center for Fundamental Research: Metal Structures in Four Dimensions;
Risoe National Laboratory;
Frederiksborgvej 399;
DK-4000 Roskilde;
email:erik.knudsen@risoe.dk
* Jon Wright, Jérôme Kieffer & Gaël Goret:
European Synchrotron Radiation Facility;
Grenoble (France)
Supports Mar345 imaging plate and Mar555 flat panel
Documentation on the format is available from:
http://rayonix.com/site_media/downloads/mar345_formats.pdf
"""
# Get ready for python3:
from __future__ import with_statement, print_function, absolute_import
__authors__ = ["Henning O. Sorensen" , "Erik Knudsen", "Jon Wright", "Jérôme Kieffer"]
__date__ = "30/10/2015"
__status__ = "production"
__copyright__ = "2007-2009 Risoe National Laboratory; 2010-2015 ESRF"
__licence__ = "GPLv3+"
from .fabioimage import FabioImage
import numpy, struct, time, sys, traceback
import logging
logger = logging.getLogger("mar345image")
from .compression import compPCK, decPCK
class Mar345Image(FabioImage):
_need_a_real_file = True
def __init__(self, *args, **kwargs):
FabioImage.__init__(self, *args, **kwargs)
self.numhigh = None
self.numpixels = None
self.swap_needed = None
def read(self, fname, frame=None):
""" Read a mar345 image"""
self.filename = fname
f = self._open(self.filename, "rb")
self._readheader(f)
if 'compressed' in self.header['Format']:
self.data = decPCK(f, self.dim1, self.dim2, self.numhigh, swap_needed=self.swap_needed)
else:
logger.error("cannot handle these formats yet " + \
"due to lack of documentation")
return None
self.bytecode = numpy.uint32
f.close()
return self
def _readheader(self, infile=None):
""" Read a mar345 image header """
# clip was not used anywhere - commented out
# clip = '\x00'
# using a couple of local variables inside this function
f = infile
h = {}
# header is 4096 bytes long
l = f.read(64)
# the contents of the mar345 header is taken to be as
# described in
# http://www.mar-usa.com/support/downloads/mar345_formats.pdf
# the first 64 bytes are 4-byte integers (but in the CBFlib
# example image it seems to 128 bytes?)
# first 4-byte integer is a marker to check endianness
if struct.unpack(" 65535)[0]
nb_pix = pix_location.size
if nb_pix % 8 == 0:
tmp = numpy.zeros((nb_pix, 2), dtype="int32")
else:
tmp = numpy.zeros(((nb_pix // 8 + 1) * 8, 2), dtype="int32")
tmp[:nb_pix, 0] = pix_location + 1
tmp[:nb_pix, 1] = flt_data[pix_location]
if self.swap_needed:
tmp.byteswap(True)
return tmp.tostring()
def nb_overflow_pixels(self):
return (self.data > 65535).sum()
@staticmethod
def checkData(data=None):
if data is None:
return None
else:
# enforce square image
shape = data.shape
assert len(shape) == 2, "image has 2 dimensions"
mshape = max(shape)
z = numpy.zeros((mshape, mshape), dtype=int)
z[:shape[0], :shape[1]] = data
return z
mar345image = Mar345Image