Content-Type: text/html; charset=UTF-8 Double Star Calculator - Orbital Elements

Double Star Calculator - Orbit

- Revision:   $Rev: 1216 $
- Report created:   2025-07-11 19:33:41

Orbital Elements for SHJ345AB
WDS: SHJ345AB WDS ID/Discoverer
P: 3500.000 Orbital period (years)
a: 14.558 Semi-major axis (arcseconds)
i: 34.3 Inclination (degrees)
Ω: 289.1 Longitude of ascending node (degrees)
T: 2020.450 Epoch of periastron passage (year)
e: 0.913 Eccentricity
ω: 158.4 Longitude of periastron (degrees)
Start-Year: 275.0 Year
End-Year: 3775.0 Year
Epoch: 2025.0 Besselian Epoch

Further data from the WDS Sixth Catalog of Orbits for SHJ345AB
Position: 222634.30-164431.9 Epoch-2000 right ascension and declination
WDS-ID: 22266-1645 WDS designation
Discoverer Code: SHJ345AB Discover designation and components
ADS number: 15934 Aitken Double Star catalog
HD number: 212697 HD catalog number
Hipparcos number: 110778 Hipparcos catalog number
V1: 6.29 Magnitude of the primary
V2: 6.39 Magnitude of the secondary
Equinox: Equinox
Date: 2019 Date of last observation
Orbit grade: 4 Orbit grade, ranging from 1 'definitive' to 5 'indeterminate'
Flags: n Flags

Double Star Orbit for SHJ345AB
--> -->
 
 
ValueError		 Python 3.12.3: /usr/bin/python3
Fri Jul 11 19:33:41 2025

A problem occurred in a Python script. Here is the sequence of function calls leading up to the error, in the order they occurred.

 /var/www/cgi-bin/ds_orbit.py in <module>
    891 
    892 # Smooth curve if necessary
=>  893 x_smooth, y_smooth = Q_get_smoothed_orbit(x_list, y_list)
    894 plt.plot(x_smooth, y_smooth, label=f'Orbit of {wds}', color='#528B8B')
    895 
x_smooth undefined, y_smooth undefined, Q_get_smoothed_orbit = <function Q_get_smoothed_orbit>, x_list = [27.241884381189195, 27.242359032356944, 27.242821142532048, 27.243270711406115, 27.243707738660703, 27.244132223967316, 27.24454416698742, 27.244943567372424, 27.245330424763686, 27.245704738792533, 27.246066509080222, 27.24641573523795, 27.24675241686688, 27.24707655355811, 27.24738814489267, 27.24768719044155, 27.247973689765665, 27.248247642415862, 27.248509047932938, 27.24875790584762, ...], y_list = [0.491385783608894, 0.4961011579616068, 0.5008163039346676, 0.5055312193556755, 0.5102459020520589, 0.5149603498509727, 0.5196745605793816, 0.5243885320639813, 0.5291022621312235, 0.5338157486073118, 0.5385289893181765, 0.5432419820894564, 0.5479547247465858, 0.5526672151146131, 0.5573794510183779, 0.5620914302823902, 0.5668031507308255, 0.5715146101875705, 0.5762258064761748, 0.5809367374199021, ...]
 /var/www/cgi-bin/ds_orbit.py in Q_get_smoothed_orbit(x_list=[27.241884381189195, 27.242359032356944, 27.242821142532048, 27.243270711406115, 27.243707738660703, 27.244132223967316, 27.24454416698742, 27.244943567372424, 27.245330424763686, 27.245704738792533, 27.246066509080222, 27.24641573523795, 27.24675241686688, 27.24707655355811, 27.24738814489267, 27.24768719044155, 27.247973689765665, 27.248247642415862, 27.248509047932938, 27.24875790584762, ...], y_list=[0.491385783608894, 0.4961011579616068, 0.5008163039346676, 0.5055312193556755, 0.5102459020520589, 0.5149603498509727, 0.5196745605793816, 0.5243885320639813, 0.5291022621312235, 0.5338157486073118, 0.5385289893181765, 0.5432419820894564, 0.5479547247465858, 0.5526672151146131, 0.5573794510183779, 0.5620914302823902, 0.5668031507308255, 0.5715146101875705, 0.5762258064761748, 0.5809367374199021, ...])
    272     # Nur wenn ausreichend Punkte vorhanden sind
    273     if len(x) >= 4:
=>  274         tck, u = splprep([x, y], s=0.0, per=True) # per=True für geschlossene Kurve (Orbit)
    275         u_fine = np.linspace(0, 1, 300)           # Feinere Parameterwerte
    276         return splev(u_fine, tck)
tck undefined, u undefined, splprep = <function splprep>, x = array([-1.24439731, -1.24985272, -1.24929712, ..., -1.19193555, -1.21533291, -1.23287335]), y = array([-0.07838142, -0.18168147, -0.28410942, ..., 0.23179831, 0.12885897, 0.02530025]), s undefined, per undefined
 /usr/lib/python3/dist-packages/scipy/interpolate/_fitpack_py.py in splprep(x=[array([-1.24439731, -1.24985272, -1.24929712, ..., -1.19193555, -1.21533291, -1.23287335]), array([-0.07838142, -0.18168147, -0.28410942, ..., 0.23179831, 0.12885897, 0.02530025])], w=None, u=None, ub=None, ue=None, k=3, task=0, s=0.0, t=None, full_output=0, nest=None, per=True, quiet=1)
    153     """
    154 
=>  155     res = _impl.splprep(x, w, u, ub, ue, k, task, s, t, full_output, nest, per,
=>  156                         quiet)
    157     return res
res undefined, global _impl = <module 'scipy.interpolate._fitpack_impl' from '...ist-packages/scipy/interpolate/_fitpack_impl.py'>, _impl.splprep = <function splprep>, x = [array([-1.24439731, -1.24985272, -1.24929712, ..., -1.19193555, -1.21533291, -1.23287335]), array([-0.07838142, -0.18168147, -0.28410942, ..., 0.23179831, 0.12885897, 0.02530025])], w = None, u = None, ub = None, ue = None, k = 3, task = 0, s = 0.0, t = None, full_output = 0, nest = None, per = True, quiet = 1
 /usr/lib/python3/dist-packages/scipy/interpolate/_fitpack_impl.py in splprep(x=array([[-1.24439731, -1.24985272, -1.24929712, .... 0.23179831, 0.12885897, -0.07838142]]), w=array([1., 1., 1., ..., 1., 1., 1.]), u=array([0. , 0.00161448, 0.00321314, ..., 0.99508629, 0.99673388, 1. ]), ub=0, ue=1, k=3, task=0, s=0.0, t=array([], dtype=float64), full_output=0, nest=3507, per=True, quiet=1)
    172     wrk = _parcur_cache['wrk']
    173     iwrk = _parcur_cache['iwrk']
=>  174     t, c, o = _fitpack._parcur(ravel(transpose(x)), w, u, ub, ue, k,
=>  175                                task, ipar, s, t, nest, wrk, iwrk, per)
    176     _parcur_cache['u'] = o['u']
t = array([], dtype=float64), c undefined, o undefined, global _fitpack = <module 'scipy.interpolate._fitpack' from '/usr/...polate/_fitpack.cpython-312-x86_64-linux-gnu.so'>, _fitpack._parcur = <built-in function _parcur>, global ravel = <function ravel>, global transpose = <function transpose>, x = array([[-1.24439731, -1.24985272, -1.24929712, .... 0.23179831, 0.12885897, -0.07838142]]), w = array([1., 1., 1., ..., 1., 1., 1.]), u = array([0. , 0.00161448, 0.00321314, ..., 0.99508629, 0.99673388, 1. ]), ub = 0, ue = 1, k = 3, task = 0, ipar = False, s = 0.0, nest = 3507, wrk = array([], dtype=float64), iwrk = array([], dtype=int32), per = True

ValueError: Invalid inputs.
      add_note = <built-in method add_note of ValueError object>
      args = ('Invalid inputs.',)
      with_traceback = <built-in method with_traceback of ValueError object>