korba 0.6.0 → 0.7.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 2e1fb6aee8c669efdda379b50b45af24fc6fabf6514972cbfc93a8156c54f533
-  data.tar.gz: cde0a77ead9cecb840aa1fad00693c5a5c61428ae6cd6f905f2faf05ec1973d7
+  metadata.gz: d032139625c6722891cef771fc76ffb90e8b082f89146ebb860a5f8034004b15
+  data.tar.gz: 354e2f0f7c382059bf7a424ffd952c32b56b6ae616b5c7bd42f50a148d2f42c3
 SHA512:
-  metadata.gz: 8981644e546a11c656d83c751aa9e9ce19046df56258287a1f9219401377feef7def725de5d4cfa75481dd156179cd93404cc89d20aa80b0761ac0fba42c4d93
-  data.tar.gz: 17306bb686d434c116d531681a6dcabd4998825e24ec042e65f1d08b84350e80a2645603ece086edf2c8d124811f641f0706011e5d32b29da687482c3ee61c5b
+  metadata.gz: 6f73a3a19efff8c89299173d67dbfa1aaf754ce1c0828148f7bd3933cd57d8de8ba74352cc12849fcd405882c964f345353e1f6d4204b4e451d00a8a6ebf5337
+  data.tar.gz: 293456b12a49007bb1f1c1acc6201fcf99678e5427b1c81942d9f1e53ab5c9cc8716ffdeaaa562846c45ef2b245400f7d73758a05494da91df1528da9e06cba0

data/CHANGELOG.md

@@ -1,3 +1,12 @@
+## [0.7.1] - 2026-04-05
+
+- FIX: Fixed calculation when orbital inclination and eccentricity are 0
+
+## [0.7.0] - 2026-04-01
+
+- FIX: Use SGP4 for TLE to Cartesian conversion
+- FIX: Remove Newton's method library and implement Kepler equation processing
+
 ## [0.6.0] - 2024-12-27
 
 - Enable use SGP4

data/README.md

@@ -70,6 +70,8 @@ kep.to_car
   - https://www.coronasha.co.jp/np/isbn/9784339012316/
 - 人工衛星の軌道 概論
   - https://www.coronasha.co.jp/np/isbn/9784339046403/
+- SGP4
+  - https://github.com/aholinch/sgp4
 
 ## Development
 

data/lib/korba/kep.rb

@@ -19,13 +19,23 @@ module Korba
     end
 
     def to_car
-      vector_n = Vector[
-        Math.sin(deg_to_rad(inclination)) * Math.sin(deg_to_rad(ra_of_asc_node)),
-        -Math.sin(deg_to_rad(inclination)) * Math.cos(deg_to_rad(ra_of_asc_node)),
-        Math.cos(deg_to_rad(inclination))
-      ]
-      vector_omega = Vector[Math.cos(deg_to_rad(ra_of_asc_node)), Math.sin(deg_to_rad(ra_of_asc_node)), 0]
-      vector_m = vector_n.cross(vector_omega)
+      # 特異点対応: inclinationが0（または極めて小さい）場合の処理
+      is_equatorial = inclination.abs < 1e-10
+
+      if is_equatorial
+        # 赤道面内の場合、x軸を基準にする
+        vector_omega = Vector[1.0, 0.0, 0.0]
+        vector_m = Vector[0.0, 1.0, 0.0]
+        vector_n = Vector[0.0, 0.0, 1.0]
+      else
+        vector_n = Vector[
+          Math.sin(deg_to_rad(inclination)) * Math.sin(deg_to_rad(ra_of_asc_node)),
+          -Math.sin(deg_to_rad(inclination)) * Math.cos(deg_to_rad(ra_of_asc_node)),
+          Math.cos(deg_to_rad(inclination))
+        ]
+        vector_omega = Vector[Math.cos(deg_to_rad(ra_of_asc_node)), Math.sin(deg_to_rad(ra_of_asc_node)), 0.0]
+        vector_m = vector_n.cross(vector_omega)
+      end
 
       r_angle_factor = deg_to_rad(arg_of_pericenter + true_anomaly)
       vector_r = distance * (Math.cos(r_angle_factor) * vector_omega + Math.sin(r_angle_factor) * vector_m)

data/lib/korba/keplers_equation.rb

@@ -0,0 +1,46 @@
+# frozen_string_literal: true
+
+module Korba
+  class KeplersEquation
+    include OrbitUtils
+
+    attr_reader :eccentricity, :mean_anomaly
+
+    def initialize(eccentricity:, mean_anomaly:)
+      super()
+      @eccentricity = eccentricity
+      @mean_anomaly = mean_anomaly
+    end
+
+    MAX_ITERATIONS = 1000
+    def solve(tolerance = 1e-10)
+      m_rad = deg_to_rad(mean_anomaly)
+
+      if eccentricity < 0.8
+        en = m_rad
+      else
+        en = Math::PI
+      end
+
+      MAX_ITERATIONS.times do
+        # f(E) = E - e*sin(E) - M
+        # f'(E) = 1 - e*cos(E)
+
+        sin_en = Math.sin(en.to_f)
+        cos_en = Math.cos(en.to_f)
+
+        f = en - eccentricity * sin_en - m_rad
+        f_dot = 1 - eccentricity * cos_en
+
+        # 更新
+        delta = f / f_dot
+        en = en - delta
+
+        # 収束判定
+        break if delta.abs < tolerance
+      end
+
+      en
+    end
+  end
+end

data/lib/korba/orbit_utils.rb

@@ -4,7 +4,7 @@ module Korba
   module OrbitUtils
     def semi_major_axis
       # a = (μ / n^2)^(1/3) m
-      (Constant::GME / (mean_motion * 2 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
+      (Constant::GME / (mean_motion * 2.0 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
     end
 
     def height_at_perigee
@@ -16,15 +16,17 @@ module Korba
     end
 
     def eccentric_anomaly
-      f = KeplerEquationFunction.new(eccentricity:, mean_anomaly:)
-      x = [deg_to_rad(mean_anomaly)]
-      nlsolve(f, x)
-      rad_to_deg(x[0])
+      kepler = KeplersEquation.new(eccentricity:, mean_anomaly:)
+      rad_to_deg(kepler.solve)
     end
 
     def true_anomaly
-      factor = (Math.cos(deg_to_rad(eccentric_anomaly)) - eccentricity) / (1 - eccentricity * Math.cos(deg_to_rad(eccentric_anomaly)))
-      rad_to_deg(Math.acos(factor))
+      e_rad = deg_to_rad(eccentric_anomaly)
+      y = Math.sqrt(1 - eccentricity ** 2) * Math.sin(e_rad)
+      x = Math.cos(e_rad) - eccentricity
+
+      # atan2(y, x) で正確な位相（-pi から pi）を求める
+      rad_to_deg(Math.atan2(y, x))
     end
 
     def distance
@@ -32,10 +34,12 @@ module Korba
     end
 
     def velocity
-      Math.sqrt(Constant::GME * (2 / distance - 1 / semi_major_axis))
+      Math.sqrt(Constant::GME * (2.0 / distance - 1.0 / semi_major_axis))
     end
 
     def path_angle
+      return 0.0 if eccentricity < 1e-15
+
       factor = Math.sqrt(Constant::GME * semi_major_axis * (1 - eccentricity ** 2)) / (distance * velocity)
       rad_to_deg(Math.acos(factor))
     end
@@ -46,8 +50,8 @@ module Korba
     end
 
     def normalize_rad(rad)
-      rad = rad + 2 * Math::PI if rad < 0
-      normalize_rad = rad > 2 * Math::PI ? rad - 2 * Math::PI : 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
     end
@@ -58,8 +62,8 @@ module Korba
     end
 
     def normalize_deg(deg)
-      deg = deg + 360 if deg < 0
-      normalized_deg = deg > 360 ? deg - 360 : 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
     end

data/lib/korba/tle.rb

@@ -5,7 +5,7 @@ module Korba
   class Tle
     include OrbitUtils
 
-    attr_reader :tle_json, :tle_string, :object_id, :object_name, :epoch_datetime, :julian_date,
+    attr_reader :tle_json, :tle_string, :sattelite_id, :object_name, :epoch_datetime, :julian_date,
                 :satellite_number, :classification_type, :epoch, :mean_motion_dot, :mean_motion_ddot, :bstar, :element_set_no,
                 :inclination, :ra_of_asc_node, :eccentricity, :arg_of_pericenter, :mean_anomaly, :mean_motion, :revolution_number,
                 :element_set_record, :epoch_days
@@ -35,8 +35,7 @@ module Korba
     end
 
     def to_car
-      kep = to_kep
-      kep.to_car
+      propagate_to(0)
     end
 
     def propagate_to(minutesAfterEpoch)
@@ -62,7 +61,7 @@ module Korba
     end
 
     def parse_line1(line1_strings)
-      # TODO: object_id
+      # TODO: sattelite_id
       @satellite_number = line1_strings[1][0..4].to_i
       @classification_type = line1_strings[1][-1]
       # For now, only supports years after 2000
@@ -103,7 +102,7 @@ module Korba
     def initialize_from_json(tle_json)
       @tle_json = tle_json
       @object_name = tle_json[:OBJECT_NAME]
-      @object_id = tle_json[:OBJECT_ID]
+      @sattelite_id = tle_json[:OBJECT_ID]
       @satellite_number = tle_json[:NORAD_CAT_ID]
       @classification_type = tle_json[:CLASSIFICATION_TYPE]
       @epoch = tle_json[:EPOCH]

data/lib/korba/version.rb

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

data/lib/korba.rb

@@ -6,7 +6,7 @@ require_relative "korba/tle"
 require_relative "korba/kep"
 require_relative "korba/car"
 require_relative "korba/orbit_utils"
-require_relative "korba/kepler_equation_fucntion"
+require_relative "korba/keplers_equation"
 require_relative "korba/sgp4/elset_rec"
 require_relative "korba/sgp4/sgp4"
 

data/mise.toml

@@ -0,0 +1,2 @@
+[tools]
+ruby = "3.4.5"

metadata

@@ -1,14 +1,13 @@
 --- !ruby/object:Gem::Specification
 name: korba
 version: !ruby/object:Gem::Version
-  version: 0.6.0
+  version: 0.7.1
 platform: ruby
 authors:
 - kk0000-kk
-autorequire:
 bindir: exe
 cert_chain: []
-date: 2024-12-26 00:00:00.000000000 Z
+date: 1980-01-02 00:00:00.000000000 Z
 dependencies: []
 description: Orbital elements calculation library
 email:
@@ -25,13 +24,13 @@ files:
 - lib/korba/car.rb
 - lib/korba/constant.rb
 - lib/korba/kep.rb
-- lib/korba/kepler_equation_fucntion.rb
-- lib/korba/newton_function.rb
+- lib/korba/keplers_equation.rb
 - lib/korba/orbit_utils.rb
 - lib/korba/sgp4/elset_rec.rb
 - lib/korba/sgp4/sgp4.rb
 - lib/korba/tle.rb
 - lib/korba/version.rb
+- mise.toml
 - sig/korba.rbs
 homepage: https://github.com/kk0000-kk/korba
 licenses: []
@@ -40,7 +39,6 @@ metadata:
   homepage_uri: https://github.com/kk0000-kk/korba
   source_code_uri: https://github.com/kk0000-kk/korba
   changelog_uri: https://github.com/kk0000-kk/korba/blob/main/CHANGELOG.md
-post_install_message:
 rdoc_options: []
 require_paths:
 - lib
@@ -55,8 +53,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.5.23
-signing_key:
+rubygems_version: 3.7.2
 specification_version: 4
 summary: Orbital elements calculation library
 test_files: []

data/lib/korba/kepler_equation_fucntion.rb

@@ -1,20 +0,0 @@
-# frozen_string_literal: true
-require_relative "newton_function"
-
-module Korba
-  class KeplerEquationFunction < NewtonFunction
-    include OrbitUtils
-
-    attr_reader :eccentricity, :mean_anomaly
-
-    def initialize(eccentricity:, mean_anomaly:)
-      super()
-      @eccentricity = eccentricity
-      @mean_anomaly = mean_anomaly
-    end
-
-    def values(x)
-      [x[0] - eccentricity * Math.sin(x[0]) - deg_to_rad(mean_anomaly)]
-    end
-  end
-end

data/lib/korba/newton_function.rb

@@ -1,25 +0,0 @@
-# frozen_string_literal: true
-require "bigdecimal/newton"
-include Newton
-
-module Korba
-  class NewtonFunction
-    def initialize()
-      @zero = BigDecimal("0.0")
-      @one = BigDecimal("1.0")
-      @two = BigDecimal("2.0")
-      @ten = BigDecimal("10.0")
-      @eps = BigDecimal("1.0e-14")
-    end
-
-    def zero; @zero; end
-    def one; @one; end
-    def two; @two; end
-    def ten; @ten; end
-    def eps; @eps; end
-
-    def values(x)
-      raise NotImplementedError
-    end
-  end
-end