diagnosis.py 文件源码

def ensure_no_stuckbits(F, args, fn): 
    bias = re.split('[\[ \] \: \,]', F[0].header['BIASSEC'])[1:-1]
    biassec = [int(t)-((i+1)%2) for i,t in enumerate(bias)]
    data = F[0].data[biassec[2]:biassec[3], biassec[0]:biassec[1]]
    mode_value = int(mode(data.ravel()).mode)
    missing_values = []
    for i in np.arange(mode_value-8, mode_value+9):
        if (data == i).sum() == 0:
            missing_values.append(i)
    for missing in missing_values:
        args.log.warning('The value %i is not represented '
                         'in the overscan region for %s' %(missing, fn))
    if len(missing_values):
        return False
    else:
        return True