Example #1
0
# Apply IERS corrections.
x += dx
y += dy
# SOFA CIO locator. */
s = erfa.s06(tt1, tt2, x, y)

# Populate the context.
astrom = erfa.apci(tt1, tt2, pvb, pvh[0], x, y, s)

# Carry out the transformation and report the results.
ri, di = erfa.atciq(rc, dc, pr, pd, px, rv, *astrom)
reprd ( "ICRS -> CIRS (JPL, IERS):", ri, di )

# The same but with Saturn then Jupiter then Sun light deflection.
b0 = erfa.LDBODY((0.00028574, 3e-10,
                  ((-7.8101442680818964, -5.6095668114887358, -1.9807981923749924),
                   (0.0030723248971152, -0.0040699547707598, -0.0018133584165345))))
b1 = erfa.LDBODY((0.00095435, 3e-9,
                  ((0.7380987962351833, .6365869247538951, 1.9693136030111202),
                   (-0.0075581692172088, 0.0012691372216750, 0.0007279990012801))))
b2 = erfa.LDBODY((1.0, 6e-6,
                  ((-0.0007121743770509, -0.0023047830339257, -0.0010586596574639),
                   (0.0000062923521264, -0.0000003308883872, -0.0000002964866231))))
b= [b0, b1, b2]
ri, di = erfa.atciqn(rc, dc, pr, pd, px, rv, b, *astrom)
reprd ( "ICRS -> CIRS (+ planets):", ri, di )

# CIRS to ICRS (astrometric).
rca, dca = erfa.aticqn(ri, di, b, *astrom)
reprd ( "CIRS -> astrometric:", rca, dca )
Example #2
0
    (
        0.00095435,
        3e-9,
        np.array(
            (
                (0.7380987962351833, 0.6365869247538951, 1.9693136030111202),
                (-0.0075581692172088, 0.0012691372216750, 0.0007279990012801),
            )
        ),
    )
)
b2 = erfa.LDBODY(
    (
        1.0,
        6e-6,
        np.array(
            (
                (-0.0007121743770509, -0.0023047830339257, -0.0010586596574639),
                (0.0000062923521264, -0.0000003308883872, -0.0000002964866231),
            )
        ),
    )
)
b = [[b0, b1, b2]]
ri, di = erfa.atciqn(rc, dc, pr, pd, px, rv, astrom, b)
reprd("ICRS -> CIRS (+ planets):", ri, di)

# CIRS to ICRS (astrometric).
rca, dca = erfa.aticqn(ri, di, astrom, b)
reprd("CIRS -> astrometric:", rca, dca)
Example #3
0
# SOFA CIO locator. */
s = erfa.s06(tt1, tt2, x, y)

# Populate the context.
astrom = erfa.apci(tt1, tt2, pvb, pvh[0], x, y, s)

# Carry out the transformation and report the results.
ri, di = erfa.atciq(rc, dc, pr, pd, px, rv, *astrom)
reprd("ICRS -> CIRS (JPL, IERS):", ri, di)

# The same but with Saturn then Jupiter then Sun light deflection.
b0 = erfa.LDBODY(
    (0.00028574, 3e-10,
     ((-7.8101442680818964, -5.6095668114887358, -1.9807981923749924),
      (0.0030723248971152, -0.0040699547707598, -0.0018133584165345))))
b1 = erfa.LDBODY(
    (0.00095435, 3e-9,
     ((0.7380987962351833, .6365869247538951, 1.9693136030111202),
      (-0.0075581692172088, 0.0012691372216750, 0.0007279990012801))))
b2 = erfa.LDBODY(
    (1.0, 6e-6,
     ((-0.0007121743770509, -0.0023047830339257, -0.0010586596574639),
      (0.0000062923521264, -0.0000003308883872, -0.0000002964866231))))
b = [b0, b1, b2]
ri, di = erfa.atciqn(rc, dc, pr, pd, px, rv, b, *astrom)
reprd("ICRS -> CIRS (+ planets):", ri, di)

# CIRS to ICRS (astrometric).
rca, dca = erfa.aticqn(ri, di, b, *astrom)
reprd("CIRS -> astrometric:", rca, dca)