Source code for G2img_CBF

# -*- coding: utf-8 -*-
########### SVN repository information ###################
# $Date: 2020-10-14 19:56:09 +0000 (Wed, 14 Oct 2020) $
# $Author: vondreele $
# $Revision: 4594 $
# $URL: https://subversion.xray.aps.anl.gov/pyGSAS/trunk/imports/G2img_CBF.py $
# $Id: G2img_CBF.py 4594 2020-10-14 19:56:09Z vondreele $
########### SVN repository information ###################
'''
*Module G2img_CBF: .cbf cif image file*
---------------------------------------

'''

from __future__ import division, print_function
import time
import GSASIIobj as G2obj
import GSASIIpath
import struct as st
import numpy as np
GSASIIpath.SetVersionNumber("$Revision: 4594 $")
[docs]class CBF_ReaderClass(G2obj.ImportImage): '''Routine to read a Read cif image data .cbf file. This is used by Pilatus. ''' def __init__(self): super(self.__class__,self).__init__( # fancy way to self-reference extensionlist=('.cbf',), strictExtension=True, formatName = 'CBF image', longFormatName = 'CIF Binary Data Format image file (NB: Slow!)' )
[docs] def ContentsValidator(self, filename): '''no test used at this time ''' return True
[docs] def Reader(self,filename, ParentFrame=None, **unused): '''Read using Bob's routine :func:`GetCbfData` ''' self.Comments,self.Data,self.Npix,self.Image = GetCbfData(self,filename) if self.Npix == 0 or not self.Comments: return False self.LoadImage(ParentFrame,filename) return True
[docs]def GetCbfData(self,filename): 'Read cif binarydetector data cbf file' import unpack_cbf as cbf if GSASIIpath.GetConfigValue('debug'): print ('Read cif binary detector data cbf file: '+filename) File = open(filename,'rb') sizexy = [0,0] pixSize = [172,172] #Pixium4700? cent = [0,0] wave = 1.54187 #default <CuKa> det2theta = 0.0 dist = 1000. byteOrd = '<' stream = File.read() if 'bytes' in str(type(stream)): stream = stream.decode('latin-1') starter = '\x0c\x1a\x04\xd5' imageBeg = stream.find(starter)+4 head = stream[:imageBeg] term = '\r\n' #CR-LF if term in head: pass else: term = '\n' #LF only head = head.split(term) for line in head: fields = line.split() if 'Wavelength' in line: wave = float(fields[2]) elif 'Detector_distance' in line: dist = float(fields[2])*1000. elif 'Pixel_size' in line: pixSize = [float(fields[2])*1.e6,float(fields[5])*1.e6] elif 'Beam_xy' in line: cent = [float(fields[2].strip('(').strip(',')),float(fields[3].strip(')'))] elif 'X-Binary-Size:' in line: compImageSize = int(fields[1]) elif 'BIG_ENDIAN' in line: byteOrd = '>' elif 'Fastest-Dimension' in line: sizexy[0] = int(fields[1]) elif 'Second-Dimension' in line: sizexy[1] = int(fields[1]) elif 'Number-of-Elements' in line: Npix = int(fields[1]) elif 'Detector_2theta' in line: det2theta = float(fields[2]) nxy = sizexy[0]*sizexy[1] cent = [cent[0]*pixSize[0]/1000.,cent[1]*pixSize[1]/1000.] File.seek(0) img = File.read()[imageBeg:imageBeg+compImageSize] File.close() nimg = len(img) image = np.zeros(nxy,dtype=np.int32) time0 = time.time() if 'bytes' in str(type(img)): img = np.frombuffer(img,dtype=np.uint8) image = cbf.unpack_cbf3(nimg,img,nxy,image) else: image = cbf.unpack_cbf(nimg,img,nxy,image) image = np.reshape(image,(sizexy[1],sizexy[0])) print ('import time: %.3f'%(time.time()-time0)) data = {'pixelSize':pixSize,'wavelength':wave,'distance':dist,'center':cent,'size':sizexy,'det2theta':det2theta} Npix = sizexy[0]*sizexy[1] return head,data,Npix,image