fwhm = eval(opts.fwhmStr) a = 1e4 # num.iinfo('uint16').max = 65535 if debug > 1: for g in gList: print g.grainSpots['predAngles'] spotParamsList = getSpotsFromGrains(gList, fwhm, a) if os.path.exists(synthStackName): os.remove(synthStackName) # writer = detector.ReadGE(None, omeMin, deltaOme, doFlip=False).getWriter(synthStackName) writer = detectorGeom.reader.getWriter(synthStackName) stack = xrdutil.makeSynthFrames(spotParamsList, detectorGeom, omegas, cutoffMult=cutoffMult, output=writer, debug=debug) writer.close() ''' # if a stack of sparse matrices is returned, can do things like: frame = stack[-1].todense() pw = detectorGeom.display(frame) detectorGeom.drawRings(pw, planeData) ''' ''' # readerTmp = reader.makeNew() readerTmp = detectorGeom.getNewReader(synthStackName, omeMin, deltaOme, doFlip=False, subtractDark=False) flatFrame = readerTmp.read(nframes=nOme, sumImg=num.maximum) pw = detectorGeom.display(flatFrame) '''
''' fwhm = eval(opts.fwhmStr) a = 1e4 # num.iinfo('uint16').max = 65535 if debug > 1: for g in gList: print g.grainSpots['predAngles'] spotParamsList = getSpotsFromGrains(gList, fwhm, a) if os.path.exists(synthStackName): os.remove(synthStackName) # writer = detector.ReadGE(None, omeMin, deltaOme, doFlip=False).getWriter(synthStackName) writer = detectorGeom.reader.getWriter(synthStackName) stack = xrdutil.makeSynthFrames(spotParamsList, detectorGeom, omegas, cutoffMult=3.0, # to keep things less computationally expensive output=writer, debug=debug) writer.close() ''' # if a stack of sparse matrices is returned, can do things like: frame = stack[-1].todense() pw = detectorGeom.display(frame) detectorGeom.drawRings(pw, planeData) ''' ''' # readerTmp = reader.makeNew() readerTmp = detectorGeom.getNewReader(synthStackName, omeMin, deltaOme, doFlip=False, subtractDark=False) flatFrame = readerTmp.read(nframes=nOme, sumImg=num.maximum) pw = detectorGeom.display(flatFrame) '''
''' fwhm = eval(opts.fwhmStr) a = 1e4 # num.iinfo('uint16').max = 65535 if debug > 1: for g in gList: print g.grainSpots['predAngles'] spotParamsList = getSpotsFromGrains(gList, fwhm, a) if os.path.exists(synthStackName): os.remove(synthStackName) # writer = detector.ReadGE(None, omeMin, deltaOme, doFlip=False).getWriter(synthStackName) writer = detectorGeom.reader.getWriter(synthStackName) stack = xrdutil.makeSynthFrames(spotParamsList, detectorGeom, omegas, cutoffMult=cutoffMult, output=writer, debug=debug) writer.close() ''' # if a stack of sparse matrices is returned, can do things like: frame = stack[-1].todense() pw = detectorGeom.display(frame) detectorGeom.drawRings(pw, planeData) ''' ''' # readerTmp = reader.makeNew() readerTmp = detectorGeom.getNewReader(synthStackName, omeMin, deltaOme, doFlip=False, subtractDark=False) flatFrame = readerTmp.read(nframes=nOme, sumImg=num.maximum) pw = detectorGeom.display(flatFrame) '''
fwhm = eval(opts.fwhmStr) a = 1e4 # num.iinfo('uint16').max = 65535 if debug > 1: for g in gList: print g.grainSpots['predAngles'] spotParamsList = getSpotsFromGrains(gList, fwhm, a) if os.path.exists(synthStackName): os.remove(synthStackName) # writer = detector.ReadGE(None, omeMin, deltaOme, doFlip=False).getWriter(synthStackName) writer = detectorGeom.reader.getWriter(synthStackName) stack = xrdutil.makeSynthFrames( spotParamsList, detectorGeom, omegas, cutoffMult=3.0, # to keep things less computationally expensive output=writer, debug=debug) writer.close() ''' # if a stack of sparse matrices is returned, can do things like: frame = stack[-1].todense() pw = detectorGeom.display(frame) detectorGeom.drawRings(pw, planeData) ''' ''' # readerTmp = reader.makeNew() readerTmp = detectorGeom.getNewReader(synthStackName, omeMin, deltaOme, doFlip=False, subtractDark=False) flatFrame = readerTmp.read(nframes=nOme, sumImg=num.maximum) pw = detectorGeom.display(flatFrame)