Changeset 2773


Ignore:
Timestamp:
Apr 5, 2017 8:21:48 AM (5 years ago)
Author:
vondreele
Message:

avoid DetDep? missing error from old gpx files
speedups in image importers

Location:
trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/GSASIIimage.py

    r2760 r2773  
    380380    x0 = dist*tand(tilt)
    381381    phi = data['rotation']
    382     dep = data['DetDepth']
     382    dep = data.get('DetDepth',0.)
    383383    azmthoff = data['azmthOff']
    384384    dx = np.array(x-cent[0],dtype=np.float32)
  • trunk/imports/G2img_1TIF.py

    r2772 r2773  
    208208                File.seek(8)
    209209                print 'Read GE-detector tiff file: ',filename
    210                 image = np.fromfile(File,dtype=np.int16,count=2*Npix)[:Npix]
     210                image = np.array(np.frombuffer(File.read(2*Npix),dtype=np.int16),dtype=np.int32)
     211#                image = np.fromfile(File,dtype=np.int16,count=2*Npix)[:Npix]
    211212#                image = np.array(ar.array('H',File.read(2*Npix)),dtype=np.int32)
    212213            elif IFD[258][2][0] == 32:
     
    251252                print 'Read APS PE-detector tiff file: ',filename
    252253                if dataType == 5:
    253                     image = np.fromfile(File,dtype=np.float32,count=4*Npix)[:Npix]
     254                    image = np.array(np.frombuffer(File.read(4*Npix),dtype=np.float32),dtype=np.int32)  #fastest
     255#                    image = np.fromfile(File,dtype=np.float32,count=4*Npix)[:Npix]
    254256#                    image = np.array(ar.array('f',File.read(4*Npix)),dtype=np.int32)
    255257                else:
     
    272274            File.seek(4096)
    273275            print 'Read MAR CCD tiff file: ',filename
    274             image = np.fromfile(File,dtype=np.int16,count=2*Npix)[:Npix]
     276            image = np.array(np.frombuffer(File.read(2*Npix),dtype=np.uint16),dtype=np.int32)
     277#            image = np.fromfile(File,dtype=np.int16,count=2*Npix)[:Npix]
    275278#            image = np.array(ar.array('H',File.read(2*Npix)),dtype=np.int32)
    276279        elif IFD[273][2][0] == 512:
  • trunk/imports/G2img_GE.py

    r2539 r2773  
    117117    import struct as st
    118118    import cPickle
     119    import time
    119120    more = False
     121    time0 = time.time()
    120122    File = open(filename,'rb')
    121123    if filename.split('.')[-1] in ['sum',]:
     
    124126        Npix = sizexy[0]*sizexy[1]
    125127        image = np.array(np.frombuffer(File.read(4*Npix),dtype=np.float32),dtype=np.int32)
    126 #        image = np.array(ar.array('f',File.read(4*Npix)),dtype=np.int32)
    127128    elif filename.split('.')[-1] in ['avg','cor']:
    128129        File.seek(0,2)
     
    134135        Npix = sizexy[0]*sizexy[1]
    135136        image = np.array(np.frombuffer(File.read(2*Npix),dtype=np.int16),dtype=np.int32)
    136 #        image = np.array(ar.array('H',File.read(2*Npix)),dtype=np.int32)
    137137    else:
    138138        head = ['GE detector raw data',]
     
    154154        File.seek(pos)
    155155        image = np.array(np.frombuffer(File.read(2*Npix),dtype=np.int16),dtype=np.int32)
    156 #        image = np.array(ar.array('H',File.read(2*Npix)),dtype=np.int32)
    157156        if len(image) != sizexy[1]*sizexy[0]:
    158157            print('not enough images while reading GE file: '+filename+'image #'+str(imagenum))
     
    164163                try:
    165164                    image += np.array(np.frombuffer(File.read(2*Npix),dtype=np.int16),dtype=np.int32)
    166 #                    image += np.array(ar.array('H',File.read(2*Npix)),dtype=np.int32)
    167165                except ValueError:
    168166                    break
     
    183181    data = {'pixelSize':[200,200],'wavelength':0.15,'distance':250.0,'center':[204.8,204.8],'size':sizexy}
    184182    File.close()
     183    print 'Image read time %.2fs'%(time.time()-time0)
    185184    if GSASIIpath.GetConfigValue('debug'):
    186185        print 'Read GE file: '+filename+' image #'+'%04d'%(imagenum)
Note: See TracChangeset for help on using the changeset viewer.