from read_GaiaDR2 import read_gaiadr2_iers_position
from cross_match import list_crossmatch, pos_max_calc, \
overall_err_calc, postional_difference_calc
from VSH_analysis import vsh_analysis
# ----------------------------- FUNCTIONS -----------------------------
# ------------------------------- MAINS -------------------------------
# Load VLBI solution
# vlbi_cat = ("/Users/Neo/Astronomy/Works/201711_GDR2_ICRF3/data/"
# "opa-sx-180425-GA15.cat")
vlbi_cat = "../data/opa-sx-180425-GA15.cat"
[
ivs_name_v, iers_name_v, ra_v, dec_v, ra_error_v, dec_error_v,
ra_dec_corr_v, num_ses, num_obs
] = read_cat(vlbi_cat)
# ellipe semi-major axis
sig_pos_max_v = pos_max_calc(ra_error_v, dec_error_v, ra_dec_corr_v)
# overall formal uncertainty
overall_err_v = overall_err_calc(ra_error_v, dec_error_v, ra_dec_corr_v)
# Load Gaia DR2 data
gaia_cat = "/Users/Neo/Astronomy/Data/catalogs/Gaia_DR2/gaiadr2_iers.fits"
[iers_name_g, ra_g, ra_error_g, dec_g, dec_error_g,
ra_dec_corr_g] = read_gaiadr2_iers_position(gaia_cat)
# ellipe semi-major axis
sig_pos_max_g = pos_max_calc(ra_error_g, dec_error_g, ra_dec_corr_g)
# --------------------------------- MAIN -------------------------------
# ICRF1 catalog
[icrf_name_i1, iers_name_i1, RA_i1, Dec_i1, e_RA_i1, e_DE_i1, corr_i1,
_] = read_icrf1_pos()
# ICRF2 catalog
[
icrf_name_i2, ivs_name_i2, iers_name_i2, RA_i2, Dec_i2, e_RA_i2, e_DE_i2,
corr_i2, _, _
] = read_icrf2()
# our solutions
# OPA-A
[ivs_name_a, iers_name_a, RA_a, Dec_a, e_RA_a, e_DE_a, corr_a, _,
_] = read_cat(
"/Users/Neo/Astronomy/Data/VLBISolutions/opa/gaia/opa-sx-180425-noGA/"
"opa-sx-180425-noGA.cat")
# OPA-B
[ivs_name_b, iers_name_b, RA_b, Dec_b, e_RA_b, e_DE_b, corr_b, _,
_] = read_cat(
"/Users/Neo/Astronomy/Data/VLBISolutions/opa/gaia/opa-sx-180425-GA00/"
"opa-sx-180425-GA00.cat")
# OPA-C
[ivs_name_c, iers_name_c, RA_c, Dec_c, e_RA_c, e_DE_c, corr_c, _,
_] = read_cat(
"/Users/Neo/Astronomy/Data/VLBISolutions/opa/gaia/opa-sx-180425-GA15/"
"opa-sx-180425-GA15.cat")
# Gaia DR1 catalog
# GaiaDR2-GA15b
gv = np.array([g1[7], g2[7], g3[7]])
glide_decomposed_plot(gv, "GaiaDR2-GA15b")
# GaiaDR2-noGAb
gv = np.array([g1[8], g2[8], g3[8]])
glide_decomposed_plot(gv, "GaiaDR2-noGAb")
# Read Catalog
# gaiadr2-timespan-noGA
[
ivs_name_nga, iers_name_nga, RA_nga, Dec_nga, RAc_err_nga, Dec_err_nga,
corr_nga, num_ses_nga, num_obs_nga
] = read_cat(
"/Users/Neo/Astronomy/Data/VLBISolutions/vlbi2_server/GalacticAberration/"
"gaiadr2-timespan-noGA/gaiadr2-timespan-noGA.cat")
# gaiadr2-timespan-noGAb
[
ivs_name_ngab, iers_name_ngab, RA_ngab, Dec_ngab, RAc_err_ngab,
Dec_err_ngab, corr_ngab, num_ses_ngab, num_obs_ngab
] = read_cat(
"/Users/Neo/Astronomy/Data/VLBISolutions/vlbi2_server/GalacticAberration/"
"gaiadr2-timespan-noGAb/gaiadr2-timespan-noGAb.cat")
# gaiadr2-timespan-GA00
[
ivs_name_g00, iers_name_g00, RA_g00, Dec_g00, RAc_err_g00, Dec_err_g00,
corr_g00, num_ses_g00, num_obs_g00
] = read_cat(