jorgepiloto

from astropy.coordinates import SkyCoord, Angle, Distance
from astropy.table import Table
from astropy import units as u


def main():

    # Read stars data into an astropy table
    # https://www.astronexus.com/downloads/catalogs/hygdata_v3.csv.gz
    stars_table = Table.read("hygdata_v3.csv")

jorgepiloto

CLA

This is just a test for a CLA.

jorgepiloto

""" Validate poliastro frames against Orekit ones """

from orekit.pyhelpers import setup_orekit_curdir
from org.hipparchus.geometry.euclidean.threed import Vector3D
from org.orekit.bodies import CelestialBodyFactory
from org.orekit.frames import FramesFactory, Transform, UpdatableFrame
from org.orekit.time import AbsoluteDate
from org.orekit.utils import PVCoordinatesProvider, TimeStampedPVCoordinates

import orekit

jorgepiloto

%----------------------------------------
%---------- Spacecraft
%----------------------------------------
% Create spacecraft and define its orbit
Create Spacecraft ss_0;
GMAT ss_0.DateFormat = TAIModJulian;
GMAT ss_0.Epoch = '21545';
GMAT ss_0.CoordinateSystem = EarthICRF;
GMAT ss_0.DisplayStateType = Cartesian;
GMAT ss_0.X = 7200;

jorgepiloto

""" Modeling Hohmann impulsive transfer with Orekit Python wrapper """

import numpy as np

# Setup the virtual machine and data location
import orekit
from orekit.pyhelpers import setup_orekit_curdir

vm = orekit.initVM()

jorgepiloto

# -*- coding: utf-8 -*-
"""
Introductory example - Plotting sin
===================================

This is a general example demonstrating a Matplotlib plot output, embedded
rST, the use of math notation and cross-linking to other examples. It would be
useful to compare the :download:`source Python file <plot_0_sin.py>` with the
output below.

jorgepiloto

from astropy.io import ascii                                                    
from astropy import units as u                                                  
from poliastro.atmosphere.models import COESA62                                 
from prettytable import PrettyTable                                             
import numpy as np                                                              
                                                                                 
coesa62 = COESA62()                                                             
                                                                                 
t = PrettyTable()                                                               
t.align='l'                                                                     

jorgepiloto

from astropy import units as u
from poliastro.twobody.angles import nu_to_M
from poliastro.core.angles import nu_to_M as nu_to_M_core

nu = 150 * u.deg
ecc = 0.08 * u.one

# If making use of radians in high-level function
M = nu_to_M(nu.to(u.rad), ecc)
print("\nHigh level function radian input results")

jorgepiloto

import matplotlib.pyplot as plt
from astropy import units as u
from astropy.time import Time

from poliastro.bodies import Earth, Mars
from poliastro.maneuver import Maneuver
from poliastro.plotting.static import StaticOrbitPlotter
from poliastro.twobody import Orbit

# Departure and time of flight for the mission

jorgepiloto

import matplotlib.pyplot as plt
from astropy import units as u
from astropy.time import Time

from poliastro.bodies import Earth
from poliastro.frames import GeocentricSolarEcliptic
from poliastro.plotting.static import StaticOrbitPlotter
from poliastro.twobody import Orbit

EPOCH = Time("2018-02-19")