korba 0.5.0 → 0.6.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, :object_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
@@ -34,6 +39,18 @@ module Korba
       kep.to_car
     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
 
     def initialize_from_string(tle_string)
@@ -46,12 +63,31 @@ module Korba
 
     def parse_line1(line1_strings)
       # TODO: object_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 +96,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]
+      @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.6.0"
 end

data/lib/korba.rb

@@ -7,6 +7,8 @@ require_relative "korba/kep"
 require_relative "korba/car"
 require_relative "korba/orbit_utils"
 require_relative "korba/kepler_equation_fucntion"
+require_relative "korba/sgp4/elset_rec"
+require_relative "korba/sgp4/sgp4"
 
 module Korba
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: korba
 version: !ruby/object:Gem::Version
-  version: 0.5.0
+  version: 0.6.0
 platform: ruby
 authors:
 - kk0000-kk
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2024-12-15 00:00:00.000000000 Z
+date: 2024-12-26 00:00:00.000000000 Z
 dependencies: []
 description: Orbital elements calculation library
 email:
@@ -28,6 +28,8 @@ files:
 - lib/korba/kepler_equation_fucntion.rb
 - lib/korba/newton_function.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
 - sig/korba.rbs