Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

My GSoC 2018 final project PR: add Gauss method and least squares for Earth- and Sun-orbiting bodies from ra/dec observations #141

Open
wants to merge 144 commits into
base: master
from

Conversation

Projects
None yet
2 participants
@PerezHz
Copy link

PerezHz commented Aug 14, 2018

This is the PR associated to my GSoC 2018 project. In this blog post all the details about my GSoC project are described.

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Aug 14, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Sep 3, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Sep 3, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Sep 3, 2018

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Sep 3, 2018

@PerezHz PerezHz force-pushed the Jorge_dev branch from b969ad8 to a88eeb5 Feb 26, 2019

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Feb 26, 2019

mu = mu_Sun # cts.GM_sun.to(uts.Unit("au3 / day2")).value
r1, r2, r3, v2, D, R, rho1, rho2, rho3, tau1, tau3, f1, g1, f3, g3, Ea_hc_pos, rho_1_sr, rho_2_sr, rho_3_sr, obs_t = gauss_estimate_mpc(mpc_object_data, mpc_observatories_data, inds, r2_root_ind=r2_root_ind)
# Apply refinement to Gauss' method, `refiters` iterations
for i in range(0,refiters):

This comment has been minimized.

@codacy-bot

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Feb 26, 2019

# compute residuals vector for ra/dec observations with pre-computed observed radec values vector
def radec_res_vec_rov_mpc(x, inds, mpc_object_data, mpc_observatories_data, rov):
"""Compute vector of observed minus computed (O-C) residuals for ra/dec
MPC-formatted observations of minor planets (asteroids, comets, etc.), with

This comment has been minimized.

# compute residuals vector for ra/dec observations with pre-computed observed radec values vector
def radec_res_vec_rov_mpc(x, inds, mpc_object_data, mpc_observatories_data, rov):
"""Compute vector of observed minus computed (O-C) residuals for ra/dec
MPC-formatted observations of minor planets (asteroids, comets, etc.), with

This comment has been minimized.

@codacy-bot

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Feb 26, 2019

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Feb 26, 2019

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Feb 26, 2019

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Feb 26, 2019

from astropy.time import Time
from datetime import datetime, timedelta
import matplotlib.pyplot as plt
from mpl_toolkits import mplot3d

This comment has been minimized.

@aerospaceresearch aerospaceresearch deleted a comment from codacy-bot Feb 26, 2019

M0 = np.mod(M,2*np.pi)
E0 = M0 + np.sign(np.sin(M0))*0.85*e #Murray-Dermotts' initial estimate
# successive approximations via Newtons' method
for i in range(0,4):

This comment has been minimized.

@codacy-bot
ra_res_vec = np.rad2deg(rv_star[0::2])*(3600.0)
dec_res_vec = np.rad2deg(rv_star[1::2])*(3600.0)

f, axarr = plt.subplots(2, sharex=True)

This comment has been minimized.

@codacy-bot
# compute residuals vector for ra/dec observations with pre-computed observed radec values vector
def radec_res_vec_rov_mpc(x, inds, mpc_object_data, mpc_observatories_data, rov):
"""Compute vector of observed minus computed (O-C) residuals for ra/dec
MPC-formatted observations of minor planets (asteroids, comets, etc.), with

This comment has been minimized.

@codacy-bot
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
You can’t perform that action at this time.