xpcs_timepixel_debug.py 文件源码

def get_timepixel_image( x,y,t, det_shape = [256, 256], delta_time = None   ):
    '''give x,y, t data to get image in a period of delta_time (in second)'''
    t0 = t.min() *6.1
    tm = t.max() *6.1

    if delta_time is not None:
        delta_time *=1e12
        if delta_time > tm:
            delta_time = tm            
    else:
        delta_time = tm
    #print( delta_time)
    t_ = t[t<delta_time]
    x_ = x[:len(t_)]
    y_ = y[:len(t_)]

    img = np.zeros( det_shape, dtype= np.int32 )
    pixlist = x_*det_shape[0] + y_ 
    his = np.histogram( pixlist, bins= np.arange( det_shape[0]*det_shape[1] +1) )[0] 
    np.ravel( img )[:] = his
    print( 'The max photon count is %d.'%img.max())
    return img