korba 0.6.0 → 0.7.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 2e1fb6aee8c669efdda379b50b45af24fc6fabf6514972cbfc93a8156c54f533
-  data.tar.gz: cde0a77ead9cecb840aa1fad00693c5a5c61428ae6cd6f905f2faf05ec1973d7
+  metadata.gz: a8a9dc1d59d577592a4ef91601bb34a4921e2bb1eb731c0fe866698759717d09
+  data.tar.gz: bb572412946fd8f08a00ab844c4a4975d08953871984a142225ffe09b786417d
 SHA512:
-  metadata.gz: 8981644e546a11c656d83c751aa9e9ce19046df56258287a1f9219401377feef7def725de5d4cfa75481dd156179cd93404cc89d20aa80b0761ac0fba42c4d93
-  data.tar.gz: 17306bb686d434c116d531681a6dcabd4998825e24ec042e65f1d08b84350e80a2645603ece086edf2c8d124811f641f0706011e5d32b29da687482c3ee61c5b
+  metadata.gz: a6f13b11b06a7dc3159b5219204a3c83791a5708f7fc8d64c0d72734bd751c655b8438ef7d29fbb2599f94dfa6657e9c282235dbf86239fab49ad12a2eafff9d
+  data.tar.gz: bcc2a2d8b20d9accbe8b6079334ce6c4610a22e78363d7c9e6243f39a1608c46a90adc0b7d0e52c53adb9fb7c81b7a4dc60f1a291e523e3166bc54b68f9a6897

data/CHANGELOG.md

@@ -1,3 +1,8 @@
+## [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/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

@@ -16,10 +16,8 @@ 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

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.0"
 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.0
 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