Changeset 2297 for trunk/imports/G2img_GE.py
- Timestamp:
- May 27, 2016 11:25:21 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/imports/G2img_GE.py
r2296 r2297 60 60 if imagenum is None: imagenum = 1 61 61 self.Comments,self.Data,self.Npix,self.Image,more = \ 62 GetGEsumData( filename,imagenum=imagenum)62 GetGEsumData(self,filename,imagenum=imagenum) 63 63 if self.Npix == 0 or not self.Comments: 64 64 return False … … 82 82 extensionlist=('.ge1','.ge2','.ge3','.ge4'), 83 83 strictExtension=True, 84 formatName = ' GE multi-image summed',84 formatName = 'sum GE multi-image', 85 85 longFormatName = 'sum of GE multi-image file' 86 86 ) … … 104 104 if imagenum is None: imagenum = 1 105 105 self.Comments,self.Data,self.Npix,self.Image,more = \ 106 GetGEsumData( filename,imagenum=imagenum,sum=True)106 GetGEsumData(self,filename,imagenum=imagenum,sum=True) 107 107 if self.Npix == 0 or not self.Comments: 108 108 return False … … 112 112 return True 113 113 114 def GetGEsumData( filename,imagenum=1,sum=False):114 def GetGEsumData(self,filename,imagenum=1,sum=False): 115 115 '''Read G.E. detector images from various files as produced at 1-ID and 116 116 with Detector Pool detector. Also sums multiple image files if desired … … 118 118 import struct as st 119 119 import array as ar 120 import cPickle 120 121 more = False 121 122 File = open(filename,'rb') … … 154 155 head += ['file: '+filename+' image #'+str(imagenum),] 155 156 if sum: #will ignore imagenum 157 print 'Frames to read %d'%(nframes), 156 158 while nframes > 1: #OK, this will sum the frames. 157 image += np.array(ar.array('H',File.read(2*Npix)),dtype=np.int32) 159 try: 160 image += np.array(ar.array('H',File.read(2*Npix)),dtype=np.int32) 161 except ValueError: 162 break 158 163 nframes -= 1 164 print '%d'%(nframes), 165 print '' 159 166 more = False 167 filename = os.path.splitext(filename)[0]+'.G2img' 168 File = open(filename,'wb') 169 Data = {'pixelSize':[200,200],'wavelength':0.15,'distance':250.0,'center':[204.8,204.8],'size':sizexy} 170 image = np.reshape(image,(sizexy[1],sizexy[0])) 171 cPickle.dump([head,Data,Npix,image],File,1) 172 File.close() 173 self.sumfile = filename 174 self.formatName = 'Summed GSAS-II image' 175 sum = False 160 176 image = np.reshape(image,(sizexy[1],sizexy[0])) 161 177 data = {'pixelSize':[200,200],'wavelength':0.15,'distance':250.0,'center':[204.8,204.8],'size':sizexy}
Note: See TracChangeset
for help on using the changeset viewer.