korba 0.1.0 → 0.2.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 9653687d464e692461e7f9d1a02bf22c969b9cbd410be77270f2774f04efde50
-  data.tar.gz: d5d043ecefc767e6ba7c158158065cfc42fc3bcda73b47e6f79315c478733b4a
+  metadata.gz: dc5cb24ea2424abd3408fb27d0bd8c414b60e02fc3923116f4b4a3f278be049c
+  data.tar.gz: 6a67a53fa7eda0a04d51bf34732a5fd82b5d97a57963e17acd28fbd64715fbf1
 SHA512:
-  metadata.gz: 64c57ea852bd747be82d9c0d42a963c1990167d92316447b6ccf338be2fc8e4c9775fb2155f1372939a1aa5d5a2c554a9e4abb213367e13afd51df5bd51b22a2
-  data.tar.gz: 90e986f11ccfa6f197df247339d39e1c928c788e36c0796027a571e09718c83ac75c4f210104ba5b209bc28c3eb183a021cb00e7c2caf08439082a1468dca97b
+  metadata.gz: e3d5ee6cf6802c4c47b04986a7ae9ebf1bc279bb76d9f8e83bebc4bdabe0521be616c662a0358195b373e22a4974c1bca2e316225e4fa79e0fd63c01e35b6920
+  data.tar.gz: c9e5748953f5d1ad950ca294c7ecf82d099f83ebbe31753453e57930dcc9aaf4ff26ffad97677efb577ad13ba429d287917ff6a47a2dad280a5d2962b5a8d8a6

data/CHANGELOG.md

@@ -1,4 +1,7 @@
-## [Unreleased]
+## [0.2.0] - 2024-12-08
+
+- Enable to create tle object from string
+- Enable to calculate height
 
 ## [0.1.0] - 2024-11-30
 

data/lib/korba/constant.rb

@@ -2,7 +2,9 @@
 
 module Korba
   class Constant
-    # 3.986004418(8)×10^14 m3/s2
+    # geocentric gravitational constant in m3/s2
     GME = 3.9860044188e14
+    # nominal equatorial Earth radius in m
+    EARTH_RADIUS = 6378137
   end
 end

data/lib/korba/kep.rb

@@ -1,7 +1,10 @@
 # frozen_string_literal: true
+require_relative "orbit_utils"
 
 module Korba
   class Kep
+    include Korba::OrbitUtils
+
     attr_reader :object_name, :epoch, :semi_major_axis, :eccentricity, :inclination, :ra_of_asc_node, :arg_of_pericenter, :mean_anomaly
 
     def initialize(object_name:, epoch:, semi_major_axis:, eccentricity:, inclination:, ra_of_asc_node:, arg_of_pericenter:, mean_anomaly:)

data/lib/korba/orbit_utils.rb

@@ -0,0 +1,18 @@
+# frozen_string_literal: true
+
+module Korba
+  module OrbitUtils
+    def semi_major_axis
+      # a = (μ / n^2)^(1/3) m
+      (Korba::Constant::GME / (mean_motion * 2 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
+    end
+
+    def height_at_perigee
+      semi_major_axis * (1 - eccentricity) - Korba::Constant::EARTH_RADIUS
+    end
+
+    def height_at_apogee
+      semi_major_axis * (1 + eccentricity) - Korba::Constant::EARTH_RADIUS
+    end
+  end
+end

data/lib/korba/tle.rb

@@ -1,24 +1,19 @@
 # frozen_string_literal: true
+require_relative "orbit_utils"
+
 module Korba
   class Tle
-    attr_reader :tle_json, :object_name, :object_id, :epoch, :mean_motion, :eccentricity, :inclination, :ra_of_asc_node, :arg_of_pericenter, :mean_anomaly
+    include Korba::OrbitUtils
 
-    def initialize(tle_json)
-      @tle_json = tle_json
-      @object_name = tle_json[:OBJECT_NAME]
-      @object_id = tle_json[:OBJECT_ID]
-      @epoch = tle_json[:EPOCH]
-      @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]
-    end
+    attr_reader :tle_json, :object_name, :object_id, :epoch, :mean_motion, :eccentricity, :inclination, :ra_of_asc_node, :arg_of_pericenter, :mean_anomaly
 
-    def semi_major_axis
-      # a = (μ / n^2)^(1/3) m
-      (Korba::Constant::GME / (mean_motion * 2 * Math::PI / 86400.0) ** 2.0) ** (1.0 / 3.0)
+    def initialize(tle, type: :string)
+      case type
+      when :string
+        initialize_from_string(tle)
+      when :json
+        initialize_from_json(tle)
+      end
     end
 
     def to_kep
@@ -31,5 +26,46 @@ module Korba
                      arg_of_pericenter:,
                      mean_anomaly:)
     end
+
+    private
+
+    def initialize_from_string(tle_string)
+      lines = tle_string.split(/\R/)
+      @object_name = lines.shift if lines.size > 2
+      parse_line1(lines[0].split(" "))
+      parse_line2(lines[1].split(" "))
+    end
+
+    def parse_line1(line1_strings)
+      # TODO: object_id
+      # 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 = epoch_time.strftime("%Y-%m-%dT%H:%M:%S.%6N")
+    end
+
+    def parse_line2(line2_strings)
+      @inclination = line2_strings[2].to_f
+      @ra_of_asc_node = line2_strings[3].to_f
+      @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
+    end
+
+    def initialize_from_json(tle_json)
+      @tle_json = tle_json
+      @object_name = tle_json[:OBJECT_NAME]
+      @object_id = tle_json[:OBJECT_ID]
+      @epoch = tle_json[:EPOCH]
+      @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]
+    end
   end
 end

data/lib/korba/version.rb

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

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: korba
 version: !ruby/object:Gem::Version
-  version: 0.1.0
+  version: 0.2.0
 platform: ruby
 authors:
 - kk0000-kk
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2024-12-07 00:00:00.000000000 Z
+date: 2024-12-08 00:00:00.000000000 Z
 dependencies: []
 description: Orbital elements calculation library
 email:
@@ -24,6 +24,7 @@ files:
 - lib/korba.rb
 - lib/korba/constant.rb
 - lib/korba/kep.rb
+- lib/korba/orbit_utils.rb
 - lib/korba/tle.rb
 - lib/korba/version.rb
 - sig/korba.rbs