korba 0.8.0 → 0.9.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 5dc3a4abbdd189f4388673e3960a615940a4feb68c359cce37a29ff47fa7efcb
-  data.tar.gz: c312b2be706e321242d0f4c762134ad82ea681348e5631ec36a13ed91c36db11
+  metadata.gz: e2c09e63837c73c4e90d0154d64f821f9ce483dff5658cc4260010ccdc34d867
+  data.tar.gz: a5cff68159d77759daa51b9ab3b6fe338bbe5f267273074b81f81da3c4abf6ee
 SHA512:
-  metadata.gz: f0cc15b9fa16cec125d3902e3b52e1336acf4246c72860049955c32456d1eaa0e56a2fd99953477d36141da23b3b8696ed0b8a611c80bc51e1249368c7f07f84
-  data.tar.gz: ffcbf5ed21c86f069d83276227c9db29e5b59d257d1630f5a1fb15eaff0550b881b1e7753fcc551ce6c647c49faf2bc84ee52efb015fbd9ad8efa9355b2da6ad
+  metadata.gz: 89c6ea54a2fc3e6499bc93b3cc6079269e3f8b31b3655339cd7d8ea9a82165bb46507fe8db71fbba2342559d466ff6df32a3afe435a679ab77631046674bb089
+  data.tar.gz: 750c1c7f2b926cc1b070550bfd6e4f18f359b0ea9658a92873a2ef3ade1ee30b63a20374490f22088d6fdf14d98fcf471cfb883900ed0476240f8a9ccfe96366

data/CHANGELOG.md

@@ -1,3 +1,12 @@
+## [0.9.1] - 2026-05-06
+
+- FIX: TLE epoch handling: convert to datetime and adjust SGP4 propagation epoch
+- feat: Implement Orbit Model
+
+## [0.9.0] - 2026-05-04
+
+- Enable orbit propagation with Keplerian elements
+
 ## [0.8.0] - 2026-04-19
 
 - Enable orbit propagation using 4th-order Runge-Kutta

data/lib/korba/converter.rb

@@ -0,0 +1,29 @@
+# frozen_string_literal: true
+
+module Korba
+  module Converter
+    def tle_to_cartesian(tle)
+      return nil if tle.nil?
+
+      tle.to_car
+    end
+
+    def tle_to_keplerian(tle)
+      return nil if tle.nil?
+
+      tle.to_kep
+    end
+
+    def keplerian_to_cartesian(keplerian)
+      return nil if keplerian.nil?
+
+      keplerian.to_car
+    end
+
+    def cartesian_to_keplerian(cartesian)
+      return nil if cartesian.nil?
+
+      cartesian.to_kep
+    end
+  end
+end

data/lib/korba/orbit.rb

@@ -0,0 +1,55 @@
+# frozen_string_literal: true
+
+module Korba
+  class Orbit
+    include Converter
+    attr_reader :epoch, :name, :tle
+
+    class << self
+      def from_tle(tle = nil, type: :string)
+        tle = Tle.new(tle, type:)
+        new(tle:, epoch: tle.epoch_datetime, name: tle.object_name)
+      end
+
+      def from_keplerian(keplerian)
+        new(keplerian:, epoch: keplerian.epoch, name: keplerian.object_name)
+      end
+
+      def from_cartesian(cartesian)
+        new(cartesian:, epoch: cartesian.epoch, name: cartesian.object_name)
+      end
+    end
+
+    def initialize(keplerian: nil, cartesian: nil, tle: nil, epoch: nil, name: nil)
+      @keplerian = keplerian
+      @cartesian = cartesian
+      @tle = tle
+      @epoch = epoch
+      @name = name
+    end
+
+    def cartesian
+      @cartesian ||= tle_to_cartesian(tle) || keplerian_to_cartesian(keplerian)
+    end
+
+    def keplerian
+      @keplerian ||= tle_to_keplerian(tle) || cartesian_to_keplerian(cartesian)
+    end
+
+    def propagate(type:, seconds_after_epoch:, disable_j2: false)
+      case type.to_sym
+      in :sgp4
+        propagated_cartesian = tle.propagate_to(seconds_after_epoch / 60.0)
+        return Orbit.from_cartesian(propagated_cartesian)
+      in :rk4
+        propagator = Propagator::Rk4.new(cartesian, disable_j2:)
+        propagated_cartesian = propagator.propagate(seconds_after_epoch)
+        return Orbit.from_cartesian(propagated_cartesian)
+      in :kepler
+        propagator = Propagator::Kepler.new(keplerian, disable_j2:)
+        propagated_keplerian = propagator.propagate(seconds_after_epoch)
+        return Orbit.from_keplerian(propagated_keplerian)
+      end
+    end
+  end
+end

data/lib/korba/orbit_utils.rb

@@ -7,6 +7,10 @@ module Korba
       (Constant::GME / (mean_motion * 2.0 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
     end
 
+    def period
+      2.0 * Math::PI * Math.sqrt(semi_major_axis ** 3 / Constant::GME)
+    end
+
     def height_at_perigee
       semi_major_axis * (1 - eccentricity) - Constant::EARTH_RADIUS
     end
@@ -50,10 +54,7 @@ module Korba
     end
 
     def normalize_rad(rad)
-      rad = rad + 2.0 * Math::PI if rad < 0
-      normalize_rad = rad > 2.0 * Math::PI ? rad - 2.0 * Math::PI : rad
-      normalize_rad(normalize_rad) if normalize_rad != rad
-      normalize_rad
+      rad % (2.0 * Math::PI)
     end
 
     def rad_to_deg(rad)
@@ -62,10 +63,7 @@ module Korba
     end
 
     def normalize_deg(deg)
-      deg = deg + 360.0 if deg < 0
-      normalized_deg = deg > 360.0 ? deg - 360.0 : deg
-      normalize_deg(normalized_deg) if normalized_deg != deg
-      normalized_deg
+      deg % 360.0
     end
   end
 end

data/lib/korba/propagator/kepler.rb

@@ -0,0 +1,54 @@
+module Korba
+  module Propagator
+    class Kepler
+      include OrbitUtils
+
+      def initialize(initial_kep, disable_j2: false)
+        @initial_kep = initial_kep
+        @disable_j2 = disable_j2
+      end
+
+      def propagate(seconds_after_epoch)
+        mean_motion = Math.sqrt(Constant::GME / @initial_kep.semi_major_axis ** 3)
+        mean_anomaly = rad_to_deg(deg_to_rad(@initial_kep.mean_anomaly) + mean_motion * seconds_after_epoch)
+
+        Kep.new(
+          object_name: @initial_kep.object_name,
+          epoch: @initial_kep.epoch + seconds_after_epoch,
+          semi_major_axis: @initial_kep.semi_major_axis,
+          eccentricity: @initial_kep.eccentricity,
+          inclination: @initial_kep.inclination,
+          ra_of_asc_node: ra_of_asc_node(seconds_after_epoch),
+          arg_of_pericenter: arg_of_pericenter(seconds_after_epoch),
+          mean_anomaly: mean_anomaly,
+        )
+      end
+
+      private
+
+      def ra_of_asc_node(seconds_after_epoch)
+        return @initial_kep.ra_of_asc_node if @disable_j2
+
+        rad_to_deg(deg_to_rad(@initial_kep.ra_of_asc_node) + ra_of_asc_node_dot_radian * seconds_after_epoch)
+      end
+
+      def ra_of_asc_node_dot_radian
+        -(3 * Math::PI * Constant::J2 *
+          (Constant::EARTH_RADIUS / @initial_kep.semi_major_axis * (1 - @initial_kep.eccentricity ** 2)) ** 2 *
+          Math.cos(deg_to_rad(@initial_kep.inclination))) / @initial_kep.period
+      end
+
+      def arg_of_pericenter(seconds_after_epoch)
+        return @initial_kep.arg_of_pericenter if @disable_j2
+
+        rad_to_deg(deg_to_rad(@initial_kep.arg_of_pericenter) + arg_of_pericenter_dot_radian * seconds_after_epoch)
+      end
+
+      def arg_of_pericenter_dot_radian
+        (3 * Math::PI / 2.0 * Constant::J2 *
+         (Constant::EARTH_RADIUS / @initial_kep.semi_major_axis * (1 - @initial_kep.eccentricity ** 2)) ** 2 *
+         (5 * (Math.cos(deg_to_rad(@initial_kep.inclination)) ** 2) - 1)) / @initial_kep.period
+      end
+    end
+  end
+end

data/lib/korba/tle.rb

@@ -25,7 +25,7 @@ module Korba
 
     def to_kep
       Kep.new(object_name:,
-              epoch:,
+              epoch: epoch_datetime,
               semi_major_axis:,
               eccentricity:,
               inclination:,
@@ -47,7 +47,7 @@ module Korba
       @sgp4Error = @element_set_record.error
       r = r.map { _1 * 1000 }
       v = v.map { _1 * 1000 }
-      Car.new(object_name:, epoch:, x: r[0], y: r[1], z: r[2], vx: v[0], vy: v[1], vz: v[2])
+      Car.new(object_name:, epoch: epoch_datetime + minutesAfterEpoch * 60.0, x: r[0], y: r[1], z: r[2], vx: v[0], vy: v[1], vz: v[2])
     end
 
     private

data/lib/korba/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module Korba
-  VERSION = "0.8.0"
+  VERSION = "0.9.1"
 end

data/lib/korba.rb

@@ -6,10 +6,13 @@ require_relative "korba/tle"
 require_relative "korba/kep"
 require_relative "korba/car"
 require_relative "korba/orbit_utils"
+require_relative "korba/converter"
+require_relative "korba/orbit"
 require_relative "korba/keplers_equation"
 require_relative "korba/sgp4/elset_rec"
 require_relative "korba/propagator/sgp4"
 require_relative "korba/propagator/rk4"
+require_relative "korba/propagator/kepler"
 
 module Korba
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: korba
 version: !ruby/object:Gem::Version
-  version: 0.8.0
+  version: 0.9.1
 platform: ruby
 authors:
 - kk0000-kk
@@ -23,9 +23,12 @@ files:
 - lib/korba.rb
 - lib/korba/car.rb
 - lib/korba/constant.rb
+- lib/korba/converter.rb
 - lib/korba/kep.rb
 - lib/korba/keplers_equation.rb
+- lib/korba/orbit.rb
 - lib/korba/orbit_utils.rb
+- lib/korba/propagator/kepler.rb
 - lib/korba/propagator/rk4.rb
 - lib/korba/propagator/sgp4.rb
 - lib/korba/sgp4/elset_rec.rb
@@ -34,7 +37,8 @@ files:
 - mise.toml
 - sig/korba.rbs
 homepage: https://github.com/kk0000-kk/korba
-licenses: []
+licenses:
+- MIT
 metadata:
   allowed_push_host: https://rubygems.org
   homepage_uri: https://github.com/kk0000-kk/korba