korba 0.5.0 → 0.7.0

data/lib/korba/tle.rb

@@ -5,7 +5,10 @@ module Korba
   class Tle
     include OrbitUtils
 
-    attr_reader :tle_json, :tle_string, :object_name, :object_id, :epoch, :mean_motion, :eccentricity, :inclination, :ra_of_asc_node, :arg_of_pericenter, :mean_anomaly
+    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
 
     def initialize(tle = nil, type: :string)
       return if tle.nil?
@@ -16,6 +19,8 @@ module Korba
       when :json
         initialize_from_json(tle)
       end
+      set_epoch_datetime_and_julian_date
+      set_element_set_record_values
     end
 
     def to_kep
@@ -30,8 +35,19 @@ module Korba
     end
 
     def to_car
-      kep = to_kep
-      kep.to_car
+      propagate_to(0)
+    end
+
+    def propagate_to(minutesAfterEpoch)
+      r = [0, 0, 0]
+      v = [0, 0, 0]
+
+      @element_set_record.error = 0
+      SGP4.sgp4(@element_set_record, minutesAfterEpoch, r, v)
+      @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])
     end
 
     private
@@ -45,13 +61,32 @@ 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
       epoch_year = line1_strings[3][0..1].to_i + 2000
-      epoch_day_of_year = line1_strings[3][2..].to_f
-      # Subtract 1 from epoch_day_of_year because January 1st is considered day 1
-      epoch_time = Time.new(epoch_year, 1, 1, 0, 0, 0, "+00:00") + (epoch_day_of_year - 1) * 24 * 60 * 60
+      @epoch_days = line1_strings[3][2..].to_f
+      # Subtract 1 from epoch_days because January 1st is considered day 1
+      epoch_time = Time.new(epoch_year, 1, 1, 0, 0, 0, "+00:00") + (epoch_days - 1) * 24 * 60 * 60
       @epoch = epoch_time.strftime("%Y-%m-%dT%H:%M:%S.%6N")
+
+      mean_motion_dot_sign = line1_strings[4][0] == "-" ? -1 : 1
+      mean_motion_dot_value_start = line1_strings[4].size == 10 ? 1 : 0
+      @mean_motion_dot = mean_motion_dot_sign * "0#{line1_strings[4][mean_motion_dot_value_start..]}".to_f
+
+      mean_motion_ddot_sign = line1_strings[5][0] == "-" ? -1 : 1
+      mean_motion_ddot_value_start = line1_strings[5].size == 8 ? 1 : 0
+      mean_motion_ddot_exponent = line1_strings[5][mean_motion_ddot_value_start + 5..mean_motion_ddot_value_start + 6].to_i
+      @mean_motion_ddot =
+        mean_motion_ddot_sign * "0.#{line1_strings[5][mean_motion_ddot_value_start..mean_motion_ddot_value_start + 4]}".to_f * 10 ** mean_motion_ddot_exponent
+
+      bstar_sign = line1_strings[6][0] == "-" ? -1 : 1
+      bstar_value_start = line1_strings[6].size == 8 ? 1 : 0
+      bstar_exponent = line1_strings[6][bstar_value_start + 5..bstar_value_start + 6].to_i
+      @bstar = bstar_sign * "0.#{line1_strings[6][bstar_value_start..bstar_value_start + 4]}".to_f * 10 ** bstar_exponent
+
+      @element_set_no = line1_strings[8][0..2].to_i
     end
 
     def parse_line2(line2_strings)
@@ -60,20 +95,62 @@ module Korba
       @eccentricity = "0.#{line2_strings[4]}".to_f
       @arg_of_pericenter = line2_strings[5].to_f
       @mean_anomaly = line2_strings[6].to_f
-      @mean_motion = line2_strings[7].to_f
+      @mean_motion = line2_strings[7][0..10].to_f
+      @revolution_number = line2_strings[7][11..15].to_i
     end
 
     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]
+      @mean_motion_dot = tle_json[:MEAN_MOTION_DOT]
+      @mean_motion_ddot = tle_json[:MEAN_MOTION_DDOT]
+      @bstar = tle_json[:BSTAR]
+      @element_set_no = tle_json[:ELEMENT_SET_NO]
       @mean_motion = tle_json[:MEAN_MOTION]
       @eccentricity = tle_json[:ECCENTRICITY]
       @inclination = tle_json[:INCLINATION]
       @ra_of_asc_node = tle_json[:RA_OF_ASC_NODE]
       @arg_of_pericenter = tle_json[:ARG_OF_PERICENTER]
       @mean_anomaly = tle_json[:MEAN_ANOMALY]
+      @revolution_number = tle_json[:REV_AT_EPOCH]
+    end
+
+    def set_epoch_datetime_and_julian_date
+      @epoch_datetime = Time.new(@epoch + " UTC")
+      @julian_date = SGP4.jday(
+        epoch_datetime.year, epoch_datetime.mon, epoch_datetime.day, epoch_datetime.hour, epoch_datetime.min, epoch_datetime.sec
+      )
+    end
+
+    def set_element_set_record_values
+      @element_set_record = ElsetRec.new
+      @element_set_record.whichconst = 3 # wgs84
+      @element_set_record.epochyr = epoch_datetime.year
+      @element_set_record.epochdays = epoch_days
+      @element_set_record.jdsatepoch = julian_date[0]
+      @element_set_record.jdsatepochF = julian_date[1]
+
+      @element_set_record.classification = classification_type
+      @element_set_record.elnum = element_set_no
+      @element_set_record.revnum = revolution_number
+      @element_set_record.satnum = satellite_number
+      @element_set_record.bstar = bstar
+      @element_set_record.inclo = deg_to_rad(inclination)
+      @element_set_record.nodeo = deg_to_rad(ra_of_asc_node)
+      @element_set_record.argpo = deg_to_rad(arg_of_pericenter)
+      @element_set_record.mo = deg_to_rad(mean_anomaly)
+      @element_set_record.ecco = eccentricity
+
+      xpdotp = 1440.0 / (2.0 * Math::PI)
+      @element_set_record.no_kozai = mean_motion / xpdotp
+      @element_set_record.ndot = mean_motion_dot / (xpdotp * 1440.0)
+      @element_set_record.nddot = mean_motion_ddot / (xpdotp * 1440.0 * 1440.0)
+      # update element_set_record
+      SGP4.sgp4init("a", @element_set_record)
     end
   end
 end

data/lib/korba/version.rb

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

data/lib/korba.rb

@@ -6,7 +6,9 @@ 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"
 
 module Korba
 end

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.5.0
+  version: 0.7.0
 platform: ruby
 authors:
 - kk0000-kk
-autorequire:
 bindir: exe
 cert_chain: []
-date: 2024-12-15 00:00:00.000000000 Z
+date: 1980-01-02 00:00:00.000000000 Z
 dependencies: []
 description: Orbital elements calculation library
 email:
@@ -25,11 +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: []
@@ -38,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
@@ -53,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