def solution_GaiaDR1_analysis(datadir, datafile, sollab): '''Comparison between the VLBI solution and Gaia DR1 ''' DiffData = sol_Gaia_diff_calc( "%s/%s" % (datadir, datafile), "%s/%s_Gaiadr1_dif.sou" % (datadir, datafile[:-4]), "%s_GaiaDR1" % sollab) # DiffData = [soucom, RAdeg, DCdeg, dRA, dRA_err, dDC, dDC_err, cov, # pos_sep, X_a, X_d, X, flgcom] # # IERS transformation parameter # print("# IERS transformation:") # # catalog_transfor(catdif, datadir, "%s_GaiaDR1" % sollab) # catalog_transfor(DiffData, # "%s/%s_Gaiadr1_dif.sou" % (datadir, datafile[:-4]), # label="%s_GaiaDR1" % sollab) # # datadir, "%s_GaiaDR1" % sollab) # VSH function parameters print("# VSH analysis:") # vsh_analysis(catdif, datadir, "%s_GaiaDR1" % sollab) vsh_analysis(DiffData, "%s/%s_Gaiadr1_dif.sou" % (datadir, datafile[:-4]), label="%s_GaiaDR1" % sollab)
def sou_position_com(cat1, cat2, comlabel): '''Compare the source position of two VLBI solutions ''' DiffData, dif_file = catlg_diff_calc(cat1, cat2, comlabel) # IERS transformation parameter print("# IERS transformation:") catalog_transfor(DiffData, dif_file, label=comlabel) print("# VSH analysis:") # VSH function parameters vsh_analysis(DiffData, dif_file, label=comlabel)
def solution_GaiaDR2_analysis(datadir, datafile, sollab): '''Comparison between the VLBI solution and Gaia DR1 ''' DiffData = sol_Gaia_diff_calc( "%s/%s" % (datadir, datafile), "%s/%s_Gaiadr2_dif.sou" % (datadir, datafile[:-4]), "%s_GaiaDR2" % sollab) # VSH function parameters print("# VSH analysis:") # vsh_analysis(cat_dif, datadir, "%s_GaiaDR1" % sollab) vsh_analysis(DiffData, "%s/%s_Gaiadr2_dif.sou" % (datadir, datafile[:-4]), label="%s_GaiaDR2" % sollab)
def solution_icrf2_analysis(datadir, datafile, sollab): '''Compare the VLBI solution with ICRF2 solution. The comparison is done by two methods: 1) Use the transformation equation recommended by IERS group; 2) Use degree-2 VSH analysis. Parameters ---------- datadir : string begining with '/' directory where the data are stored. Returns ---------- None ''' # catdif = "%s/%s_icrf2_dif.sou" % (datadir, datafile[:-4]) DiffData = sol_icrf2_diff_calc( "%s/%s" % (datadir, datafile), "%s/%s_icrf2_dif.sou" % (datadir, datafile[:-4]), "%s_icrf2" % sollab) # # IERS transformation parameters. # print("# IERS transformation:") # # catalog_transfor(catdif, datadir, "%s_icrf2" % sollab) # catalog_transfor(DiffData, # "%s/%s_icrf2_dif.sou" % (datadir, datafile[:-4]), # label="%s_icrf2" % sollab) # # datadir, "%s_icrf2" % sollab) # VSH function parameters. print("# VSH analysis:") # vsh_analysis(catdif, datadir, "%s_icrf2" % sollab) vsh_analysis(DiffData, "%s/%s_icrf2_dif.sou" % (datadir, datafile[:-4]), label="%s_icrf2" % sollab)
"# Created date: %s\n#" % (vlbi_cat, gaia_cat, time.strftime("%d/%m/%Y", time.localtime())), file=foffset) for (comsou_ivsi, comsou_iersi, ra_com_gi, dec_com_gi, sig_pos_max_com_gi, sig_pos_max_com_vi, num_ses_comi, drai, ddeci, dra_erri, ddec_erri, dra_ddec_covi, ang_sepi, Xai, Xdi, Xi, X2i) in zip(comsou_ivs, comsou_iers, ra_com_g, dec_com_g, sig_pos_max_com_g, sig_pos_max_com_v, num_ses_com, dra, ddec, dra_err, ddec_err, dra_ddec_cov, ang_sep, Xa, Xd, X, X2): print("%-8s %8s %14.9f %+14.9f %8.3f %8.3f %5d %+8.3f %+8.3f " "%8.3f %8.3f %+5.2f %8.3f %8.3f %8.3f %8.3f %8.3f" % (comsou_ivsi, comsou_iersi, ra_com_gi, dec_com_gi, sig_pos_max_com_gi, sig_pos_max_com_vi, num_ses_comi, drai, ddeci, dra_erri, ddec_erri, dra_ddec_covi, ang_sepi, Xai, Xdi, Xi, X2i), file=foffset) foffset.close() flg = np.empty_like(comsou_iers) pos_offset = [ comsou_iers, ra_com_g, dec_com_g, dra, dra_err, ddec, ddec_err, dra_ddec_cov, ang_sep, Xa, Xd, X, flg ] # VSH analysis vsh_analysis(pos_offset, vlbi_cat, label="opa-sx-180425-GA15") # --------------------------------- END --------------------------------
# O for non-icrf2 source tp.append('O') post_diff_plot(dRA, dRA_err, dDC, dDC_err, DC1n, np.array(tp)) fdif.close() # -------------------- MAIN ---------------------------------- # # Test # cat1 = "/home/nliu/solutions/test/a1/result_a1.sou" # cat2 = "/home/nliu/solutions/test/a2/result_a2.sou" # catdif = "/home/nliu/solutions/test/a1_a2_dif.sou" # catlg_diff_calc(cat1, cat2, catdif) # # cat_transfor(catdif) # vsh_analysis(catdif) # GA - Non-GA cat1 = "/home/nliu/solutions/test/GA/opa2017a_aprx.sou" cat2 = "/home/nliu/solutions/test/GA/opa2017a_ga.sou" catdif = "/home/nliu/solutions/test/GA/opa2017a_ga_dif.sou" catlg_diff_calc(cat1, cat2, catdif) cat_transfor(catdif) vsh_analysis(catdif) # --------------------------------- END -------------------------------