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 -------------------------------