check_result new

import pandas as pd
import os

'''
modify time:2022-01-19 15:10:28
editor:zhaoyuanyuan
'''
# start_time = time.time()
# yesterday = datetime.today()+timedelta(-1)
# yesterday_format = yesterday.strftime('%Y_%m_%d')
yesterday_format = '2022_03_01'
fitspath = "/data2/data/" + yesterday_format + "/"
# fitspath = "/data2/data/2022_03_01/"
wcspath = '/data1/result/'+yesterday_format+ "/"
print(fitspath)


if (os.path.exists(fitspath)):
    print(',,,,')
    os.system('ls ' + fitspath + '*fits >fits.dat')
    os.system('sort fits.dat -k1.84 >sfits.dat')
    os.system('ls '+wcspath+'*.wcs >wcs.dat')
    ls_fits = []
    ls_obstime =[]
    ls_exptime = []
    blist = [line.rstrip() for line in open('sfits.dat')]
    for j in blist:
        try:
            w = j[-56:-49]
            w1 = int(w)
            tm = j[-32:-18]
            expt = float(j[-44:-37])
            ls_fits.append(w1)
            ls_obstime.append(tm)
            ls_exptime.append(expt)

        except Exception:
            pass
        continue
    fits_obstdict = dict(zip(ls_fits,ls_obstime))
    fits_exptdict = dict(zip(ls_fits,ls_exptime))
    print(len(ls_fits),'len fits')
    ls_obj = list(set(ls_fits))
    ls_lenobj = []
    for fnum in ls_obj:
        lenfnum = ls_fits.count(fnum)
        ls_lenobj.append(lenfnum)

    ls_wcs = []
    alist = [line.rstrip() for line in open('wcs.dat')]

    for j in alist:
        try:
            a = j[-55:-48]
            wcsnum1 = int(a)
            ls_wcs.append(wcsnum1)
        except Exception:
            pass
        continue
    print(len(ls_wcs), 'len wcs')
    wcsnum = list(set(ls_wcs))
    lenobjwcs = []
    lenobjfit = []
    lenwcsprobality=[]
    objobstime =[]
    objexptime = []
    for target in ls_obj:
        # if target in wcsnum:
        lenwcs =  ls_wcs.count(target)
        lenfit = ls_fits.count(target)
        probality = float('%.2f'%(lenwcs/lenfit))
        obstime = fits_obstdict[target]
        exptime = fits_exptdict[target]
        objobstime.append(obstime)
        objexptime.append(exptime)
        lenobjfit.append(lenfit)
        lenobjwcs.append(lenwcs)
        lenwcsprobality.append(probality)
    df = pd.DataFrame({'OBJ':ls_obj,'obs_time':objobstime,'exp_time':objexptime,'len_fit':lenobjfit,
                       'len_wcs':lenobjwcs,'probablity':lenwcsprobality})

    df.to_csv('wcsresult.xlsx')