Пример #1
0
rca, dca, eo = erfa.atic13 ( ri, di, tt1, tt2)
#
reprd ( "CIRS -> astrometric:", rca, dca )

#ICRS (astrometric) to CIRS (geocentric observer).
ri, di, eo = erfa.atci13 ( rca, dca, 0.0, 0.0, 0.0, 0.0, tt1, tt2)
reprd ( "astrometric -> CIRS:", ri, di );

# Apparent place.
ra = erfa.anp ( ri - eo )
da = di
reprd ( "geocentric apparent:", ra, da )

# CIRS to topocentric.
aot, zot, hot, dot, rot = erfa.atio13(ri, di, utc1, utc2, dut1,
                                      elong, phi, hm, xp, yp,
                                      0.0, 0.0, 0.0, 0.0)
reprd( "CIRS -> topocentric:", rot, dot )

# CIRS to topocentric.
aob, zob, hob, dob, rob = erfa.atio13(ri, di, utc1, utc2, dut1,
                                      elong, phi, hm, xp, yp,
                                      phpa, tc, rh, wl)
reprd( "CIRS -> observed:", rob, dob )

# ICRS to observed.
aob, zob, hob, dob, rob, eo = erfa.atco13 ( rc, dc, pr, pd, px, rv,
                                            utc1, utc2, dut1,
                                            elong, phi, hm, xp, yp,
                                            phpa, tc, rh, wl)
reprd ( "ICRS -> observed:", rob, dob )
Пример #2
0
rca, dca, eo = erfa.atic13(ri, di, tt1, tt2)
#
reprd("CIRS -> astrometric:", rca, dca)

# ICRS (astrometric) to CIRS (geocentric observer).
ri, di, eo = erfa.atci13(rca, dca, np.array([0.0]), np.array([0.0]), np.array([0.0]), np.array([0.0]), tt1, tt2)
reprd("astrometric -> CIRS:", ri, di)

# Apparent place.
ra = erfa.anp(ri - eo)
da = di
reprd("geocentric apparent:", ra, da)

# CIRS to topocentric.
aot, zot, hot, dot, rot = erfa.atio13(
    ri, di, utc1, utc2, dut1, elong, phi, hm, xp, yp, np.array([0.0]), np.array([0.0]), np.array([0.0]), np.array([0.0])
)
reprd("CIRS -> topocentric:", rot, dot)

# CIRS to topocentric.
aob, zob, hob, dob, rob = erfa.atio13(ri, di, utc1, utc2, dut1, elong, phi, hm, xp, yp, phpa, tc, rh, wl)
reprd("CIRS -> observed:", rob, dob)

# ICRS to observed.
aob, zob, hob, dob, rob, eo = erfa.atco13(
    rc, dc, pr, pd, px, rv, utc1, utc2, dut1, elong, phi, hm, xp, yp, phpa, tc, rh, wl
)
reprd("ICRS -> observed:", rob, dob)

# ICRS to CIRS using some user-supplied parameters.
# SOFA heliocentric Earth ephemeris.
Пример #3
0
# CIRS to ICRS (astrometric).
rca, dca, eo = erfa.atic13(ri, di, tt1, tt2)
#
reprd("CIRS -> astrometric:", rca, dca)

#ICRS (astrometric) to CIRS (geocentric observer).
ri, di, eo = erfa.atci13(rca, dca, 0.0, 0.0, 0.0, 0.0, tt1, tt2)
reprd("astrometric -> CIRS:", ri, di)

# Apparent place.
ra = erfa.anp(ri - eo)
da = di
reprd("geocentric apparent:", ra, da)

# CIRS to topocentric.
aot, zot, hot, dot, rot = erfa.atio13(ri, di, utc1, utc2, dut1, elong, phi, hm,
                                      xp, yp, 0.0, 0.0, 0.0, 0.0)
reprd("CIRS -> topocentric:", rot, dot)

# CIRS to topocentric.
aob, zob, hob, dob, rob = erfa.atio13(ri, di, utc1, utc2, dut1, elong, phi, hm,
                                      xp, yp, phpa, tc, rh, wl)
reprd("CIRS -> observed:", rob, dob)

# ICRS to observed.
aob, zob, hob, dob, rob, eo = erfa.atco13(rc, dc, pr, pd, px, rv, utc1, utc2,
                                          dut1, elong, phi, hm, xp, yp, phpa,
                                          tc, rh, wl)
reprd("ICRS -> observed:", rob, dob)

# ICRS to CIRS using some user-supplied parameters.
# SOFA heliocentric Earth ephemeris.