korba 0.4.0 → 0.6.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: e02130eb6662759d8869d606d53c5d7ecc64c26ed59a5605487fdf5e8d785f01
-  data.tar.gz: 90d8fe4257842deffd1f45ca5fab34013f9decb275ff03f9d59c369f3182e11e
+  metadata.gz: 2e1fb6aee8c669efdda379b50b45af24fc6fabf6514972cbfc93a8156c54f533
+  data.tar.gz: cde0a77ead9cecb840aa1fad00693c5a5c61428ae6cd6f905f2faf05ec1973d7
 SHA512:
-  metadata.gz: 95b4ade626d6622599c3d7e6b533db46a49b820c81cdbc3a8f57838901220950dfdb2643de622d8a414c63a4a24fa15831a8dc489c821ec15c1d10335a929d6c
-  data.tar.gz: 84b650670d1d10fc81b0cd9e9b5005a2dba965bb119df8913a7e38d0997b0c0f31f53ba105d9c3d13198ae4b768df5c7f23b5eb8e3d4694ea86bc7bf0b3e1ff6
+  metadata.gz: 8981644e546a11c656d83c751aa9e9ce19046df56258287a1f9219401377feef7def725de5d4cfa75481dd156179cd93404cc89d20aa80b0761ac0fba42c4d93
+  data.tar.gz: 17306bb686d434c116d531681a6dcabd4998825e24ec042e65f1d08b84350e80a2645603ece086edf2c8d124811f641f0706011e5d32b29da687482c3ee61c5b

data/CHANGELOG.md

@@ -1,3 +1,11 @@
+## [0.6.0] - 2024-12-27
+
+- Enable use SGP4
+
+## [0.5.0] - 2024-12-15
+
+- Enable transformation to Cartesian elements
+
 ## [0.4.0] - 2024-12-09
 
 - Enable save tle_string

data/README.md

@@ -38,7 +38,7 @@ orbit.height_at_perigee
 orbit.height_at_apogee
 # => 420372.03623388056
 
-orbit.to_kep
+kep = orbit.to_kep
 # =>
 # #<Korba::Kep:xxxxxxxxxxxxxxxx
 #  @arg_of_pericenter=314.0303,
@@ -49,8 +49,28 @@ orbit.to_kep
 #  @object_name="ISS (ZARYA)",
 #  @ra_of_asc_node=174.9565,
 #  @semi_major_axis=6793877.649839985>
+
+kep.to_car
+# =>
+# #<Korba::Car:xxxxxxxxxxxxxxxx
+#  @epoch="2024-12-07T20:37:24.085055",
+#  @object_name="ISS (ZARYA)",
+#  @vx=6150.772410883998,
+#  @vy=2489.3298780751356,
+#  @vz=-3816.0301666253677,
+#  @x=4019753.8621700387,
+#  @y=-3623966.518673545,
+#  @z=4114361.6934797494>
+
 ```
 
+## References
+
+- 宇宙工学シリーズ 3 人工衛星と宇宙探査機
+  - https://www.coronasha.co.jp/np/isbn/9784339012316/
+- 人工衛星の軌道 概論
+  - https://www.coronasha.co.jp/np/isbn/9784339046403/
+
 ## Development
 
 After checking out the repo, run `bin/setup` to install dependencies. Then, run `rake spec` to run the tests. You can also run `bin/console` for an interactive prompt that will allow you to experiment.

data/lib/korba/car.rb

@@ -0,0 +1,18 @@
+# frozen_string_literal: true
+
+module Korba
+  class Car
+    attr_reader :object_name, :epoch, :x, :y, :z, :vx, :vy, :vz
+
+    def initialize(object_name:, epoch:, x:, y:, z:, vx:, vy:, vz:)
+      @object_name = object_name
+      @epoch = epoch
+      @x = x
+      @y = y
+      @z = z
+      @vx = vx
+      @vy = vy
+      @vz = vz
+    end
+  end
+end

data/lib/korba/kep.rb

@@ -1,9 +1,9 @@
 # frozen_string_literal: true
-require_relative "orbit_utils"
+require "matrix"
 
 module Korba
   class Kep
-    include Korba::OrbitUtils
+    include OrbitUtils
 
     attr_reader :object_name, :epoch, :semi_major_axis, :eccentricity, :inclination, :ra_of_asc_node, :arg_of_pericenter, :mean_anomaly
 
@@ -17,5 +17,22 @@ module Korba
       @arg_of_pericenter = arg_of_pericenter
       @mean_anomaly = mean_anomaly
     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)
+
+      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)
+      v_angle_factor = deg_to_rad(arg_of_pericenter + true_anomaly - path_angle)
+      vector_v = velocity * (-Math.sin(v_angle_factor) * vector_omega + (Math.cos(v_angle_factor)) * vector_m)
+
+      Car.new(object_name:, epoch:, x: vector_r[0], y: vector_r[1], z: vector_r[2], vx: vector_v[0], vy: vector_v[1], vz: vector_v[2])
+    end
   end
 end

data/lib/korba/kepler_equation_fucntion.rb

@@ -0,0 +1,20 @@
+# 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

@@ -0,0 +1,25 @@
+# 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

data/lib/korba/orbit_utils.rb

@@ -4,15 +4,64 @@ 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)
+      (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
+      semi_major_axis * (1 - eccentricity) - Constant::EARTH_RADIUS
     end
 
     def height_at_apogee
-      semi_major_axis * (1 + eccentricity) - Korba::Constant::EARTH_RADIUS
+      semi_major_axis * (1 + eccentricity) - Constant::EARTH_RADIUS
+    end
+
+    def eccentric_anomaly
+      f = KeplerEquationFunction.new(eccentricity:, mean_anomaly:)
+      x = [deg_to_rad(mean_anomaly)]
+      nlsolve(f, x)
+      rad_to_deg(x[0])
+    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))
+    end
+
+    def distance
+      semi_major_axis * (1 - eccentricity * Math.cos(deg_to_rad(eccentric_anomaly)))
+    end
+
+    def velocity
+      Math.sqrt(Constant::GME * (2 / distance - 1 / semi_major_axis))
+    end
+
+    def path_angle
+      factor = Math.sqrt(Constant::GME * semi_major_axis * (1 - eccentricity ** 2)) / (distance * velocity)
+      rad_to_deg(Math.acos(factor))
+    end
+
+    def deg_to_rad(deg)
+      rad = deg * Math::PI / 180.0
+      normalize_rad(rad)
+    end
+
+    def normalize_rad(rad)
+      rad = rad + 2 * Math::PI if rad < 0
+      normalize_rad = rad > 2 * Math::PI ? rad - 2 * Math::PI : rad
+      normalize_rad(normalize_rad) if normalize_rad != rad
+      normalize_rad
+    end
+
+    def rad_to_deg(rad)
+      deg = rad * 180.0 / Math::PI
+      normalize_deg(deg)
+    end
+
+    def normalize_deg(deg)
+      deg = deg + 360 if deg < 0
+      normalized_deg = deg > 360 ? deg - 360 : deg
+      normalize_deg(normalized_deg) if normalized_deg != deg
+      normalized_deg
     end
   end
 end

data/lib/korba/sgp4/elset_rec.rb

@@ -0,0 +1,432 @@
+module Korba
+  #
+  # This class implements the elsetrec data type from Vallado's SGP4 code.
+  #
+  # From SGP4.h
+  # #define SGP4Version  "SGP4 Version 2016-03-09"
+  #
+  # @author aholinch
+  #
+  #
+  class ElsetRec
+    attr_accessor :whichconst
+    attr_accessor :satnum
+    attr_accessor :epochyr
+    attr_accessor :epochtynumrev
+    attr_accessor :error
+    attr_accessor :operationmode
+    attr_accessor :init
+    attr_accessor :method
+    attr_accessor :a
+    attr_accessor :altp
+    attr_accessor :alta
+    attr_accessor :epochdays
+    attr_accessor :jdsatepoch
+    attr_accessor :jdsatepochF
+    attr_accessor :nddot
+    attr_accessor :ndot
+    attr_accessor :bstar
+    attr_accessor :rcse
+    attr_accessor :inclo
+    attr_accessor :nodeo
+    attr_accessor :ecco
+    attr_accessor :argpo
+    attr_accessor :mo
+    attr_accessor :no_kozai
+    attr_accessor :classification
+    attr_accessor :intldesg
+    attr_accessor :ephtype
+    attr_accessor :elnum
+    attr_accessor :revnum
+    attr_accessor :no_unkozai
+    attr_accessor :am
+    attr_accessor :em
+    attr_accessor :im
+    attr_accessor :Om
+    attr_accessor :om
+    attr_accessor :mm
+    attr_accessor :nm
+    attr_accessor :t
+    attr_accessor :tumin
+    attr_accessor :mu
+    attr_accessor :radiusearthkm
+    attr_accessor :xke
+    attr_accessor :j2
+    attr_accessor :j3
+    attr_accessor :j4
+    attr_accessor :j3oj2
+    attr_accessor :dia_mm
+    attr_accessor :period_sec
+    attr_accessor :active
+    attr_accessor :not_orbital
+    attr_accessor :rcs_m2
+    attr_accessor :ep
+    attr_accessor :inclp
+    attr_accessor :nodep
+    attr_accessor :argpp
+    attr_accessor :mp
+    attr_accessor :isimp
+    attr_accessor :aycof
+    attr_accessor :con41
+    attr_accessor :cc1
+    attr_accessor :cc4
+    attr_accessor :cc5
+    attr_accessor :d2
+    attr_accessor :d3
+    attr_accessor :d4
+    attr_accessor :delmo
+    attr_accessor :eta
+    attr_accessor :argpdot
+    attr_accessor :omgcof
+    attr_accessor :sinmao
+    attr_accessor :t2cof
+    attr_accessor :t3cof
+    attr_accessor :t4cof
+    attr_accessor :t5cof
+    attr_accessor :x1mth2
+    attr_accessor :x7thm1
+    attr_accessor :mdot
+    attr_accessor :nodedot
+    attr_accessor :xlcof
+    attr_accessor :xmcof
+    attr_accessor :nodecf
+    attr_accessor :irez
+    attr_accessor :d2201
+    attr_accessor :d2211
+    attr_accessor :d3210
+    attr_accessor :d3222
+    attr_accessor :d4410
+    attr_accessor :d4422
+    attr_accessor :d5220
+    attr_accessor :d5232
+    attr_accessor :d5421
+    attr_accessor :d5433
+    attr_accessor :dedt
+    attr_accessor :del1
+    attr_accessor :del2
+    attr_accessor :del3
+    attr_accessor :didt
+    attr_accessor :dmdt
+    attr_accessor :dnodt
+    attr_accessor :domdt
+    attr_accessor :e3
+    attr_accessor :ee2
+    attr_accessor :peo
+    attr_accessor :pgho
+    attr_accessor :pho
+    attr_accessor :pinco
+    attr_accessor :plo
+    attr_accessor :se2
+    attr_accessor :se3
+    attr_accessor :sgh2
+    attr_accessor :sgh3
+    attr_accessor :sgh4
+    attr_accessor :sh2
+    attr_accessor :sh3
+    attr_accessor :si2
+    attr_accessor :si3
+    attr_accessor :sl2
+    attr_accessor :sl3
+    attr_accessor :sl4
+    attr_accessor :gsto
+    attr_accessor :xfact
+    attr_accessor :xgh2
+    attr_accessor :xgh3
+    attr_accessor :xgh4
+    attr_accessor :xh2
+    attr_accessor :xh3
+    attr_accessor :xi2
+    attr_accessor :xi3
+    attr_accessor :xl2
+    attr_accessor :xl3
+    attr_accessor :xl4
+    attr_accessor :xlamo
+    attr_accessor :zmol
+    attr_accessor :zmos
+    attr_accessor :atime
+    attr_accessor :xli
+    attr_accessor :xni
+    attr_accessor :snodm
+    attr_accessor :cnodm
+    attr_accessor :sinim
+    attr_accessor :cosim
+    attr_accessor :sinomm
+    attr_accessor :cosomm
+    attr_accessor :day
+    attr_accessor :emsq
+    attr_accessor :gam
+    attr_accessor :rtemsq
+    attr_accessor :s1
+    attr_accessor :s2
+    attr_accessor :s3
+    attr_accessor :s4
+    attr_accessor :s5
+    attr_accessor :s6
+    attr_accessor :s7
+    attr_accessor :ss1
+    attr_accessor :ss2
+    attr_accessor :ss3
+    attr_accessor :ss4
+    attr_accessor :ss5
+    attr_accessor :ss6
+    attr_accessor :ss7
+    attr_accessor :sz1
+    attr_accessor :sz2
+    attr_accessor :sz3
+    attr_accessor :sz11
+    attr_accessor :sz12
+    attr_accessor :sz13
+    attr_accessor :sz21
+    attr_accessor :sz22
+    attr_accessor :sz23
+    attr_accessor :sz31
+    attr_accessor :sz32
+    attr_accessor :sz33
+    attr_accessor :z1
+    attr_accessor :z2
+    attr_accessor :z3
+    attr_accessor :z11
+    attr_accessor :z12
+    attr_accessor :z13
+    attr_accessor :z21
+    attr_accessor :z22
+    attr_accessor :z23
+    attr_accessor :z31
+    attr_accessor :z32
+    attr_accessor :z33
+    attr_accessor :argpm
+    attr_accessor :inclm
+    attr_accessor :nodem
+    attr_accessor :dndt
+    attr_accessor :eccsq
+    attr_accessor :ainv
+    attr_accessor :ao
+    attr_accessor :con42
+    attr_accessor :cosio
+    attr_accessor :cosio2
+    attr_accessor :omeosq
+    attr_accessor :posq
+    attr_accessor :rp
+    attr_accessor :rteosq
+    attr_accessor :sinio
+
+    def initialize
+      @whichconst = 2    #SGP4.wgs72
+      @satnum = 0
+      @epochyr = 0
+      @epochtynumrev = 0
+      @error = 0
+      @operationmode = 0
+      @init = 0
+      @method = 0
+      @a = 0.0
+      @altp = 0.0
+      @alta = 0.0
+      @epochdays = 0.0
+      @jdsatepoch = 0.0
+      @jdsatepochF = 0.0
+      @nddot = 0.0
+      @ndot = 0.0
+      @bstar = 0.0
+      @rcse = 0.0
+      @inclo = 0.0
+      @nodeo = 0.0
+      @ecco = 0.0
+      @argpo = 0.0
+      @mo = 0.0
+      @no_kozai = 0.0
+
+      # sgp4fix add new variables from tle
+      @classification = "U"
+      @intldesg = ""
+      @ephtype = 0
+      @elnum = 0
+      @revnum = 0
+
+      # sgp4fix add unkozai'd variable
+      @no_unkozai = 0.0
+      # sgp4fix add singly averaged variables
+      @am = 0.0
+      @em = 0.0
+      @im = 0.0
+      @Om = 0.0
+      @om = 0.0
+      @mm = 0.0
+      @nm = 0.0
+      @t = 0.0
+
+      # sgp4fix add constant parameters to eliminate mutliple calls during execution
+      @tumin = 0.0
+      @mu = 0.0
+      @radiusearthkm = 0.0
+      @xke = 0.0
+      @j2 = 0.0
+      @j3 = 0.0
+      @j4 = 0.0
+      @j3oj2 = 0.0
+
+      #       Additional elements to capture relevant TLE and object information:
+      @dia_mm = 0 # RSO dia in mm
+      @period_sec = 0.0 # Period in seconds
+      @active = 0 # "Active S/C" flag (0=n, 1=y)
+      @not_orbital = 0 # "Orbiting S/C" flag (0=n, 1=y)
+      @rcs_m2 = 0.0 # "RCS (m^2)" storage
+
+      # temporary variables because the original authors call the same method with different variables
+      @ep = 0.0
+      @inclp = 0.0
+      @nodep = 0.0
+      @argpp = 0.0
+      @mp = 0.0
+
+      @isimp = 0
+      @aycof = 0.0
+      @con41 = 0.0
+      @cc1 = 0.0
+      @cc4 = 0.0
+      @cc5 = 0.0
+      @d2 = 0.0
+      @d3 = 0.0
+      @d4 = 0.0
+      @delmo = 0.0
+      @eta = 0.0
+      @argpdot = 0.0
+      @omgcof = 0.0
+      @sinmao = 0.0
+      @t2cof = 0.0
+      @t3cof = 0.0
+      @t4cof = 0.0
+      @t5cof = 0.0
+      @x1mth2 = 0.0
+      @x7thm1 = 0.0
+      @mdot = 0.0
+      @nodedot = 0.0
+      @xlcof = 0.0
+      @xmcof = 0.0
+      @nodecf = 0.0
+
+      # deep space
+      @irez = 0
+      @d2201 = 0.0
+      @d2211 = 0.0
+      @d3210 = 0.0
+      @d3222 = 0.0
+      @d4410 = 0.0
+      @d4422 = 0.0
+      @d5220 = 0.0
+      @d5232 = 0.0
+      @d5421 = 0.0
+      @d5433 = 0.0
+      @dedt = 0.0
+      @del1 = 0.0
+      @del2 = 0.0
+      @del3 = 0.0
+      @didt = 0.0
+      @dmdt = 0.0
+      @dnodt = 0.0
+      @domdt = 0.0
+      @e3 = 0.0
+      @ee2 = 0.0
+      @peo = 0.0
+      @pgho = 0.0
+      @pho = 0.0
+      @pinco = 0.0
+      @plo = 0.0
+      @se2 = 0.0
+      @se3 = 0.0
+      @sgh2 = 0.0
+      @sgh3 = 0.0
+      @sgh4 = 0.0
+      @sh2 = 0.0
+      @sh3 = 0.0
+      @si2 = 0.0
+      @si3 = 0.0
+      @sl2 = 0.0
+      @sl3 = 0.0
+      @sl4 = 0.0
+      @gsto = 0.0
+      @xfact = 0.0
+      @xgh2 = 0.0
+      @xgh3 = 0.0
+      @xgh4 = 0.0
+      @xh2 = 0.0
+      @xh3 = 0.0
+      @xi2 = 0.0
+      @xi3 = 0.0
+      @xl2 = 0.0
+      @xl3 = 0.0
+      @xl4 = 0.0
+      @xlamo = 0.0
+      @zmol = 0.0
+      @zmos = 0.0
+      @atime = 0.0
+      @xli = 0.0
+      @xni = 0.0
+      @snodm = 0.0
+      @cnodm = 0.0
+      @sinim = 0.0
+      @cosim = 0.0
+      @sinomm = 0.0
+      @cosomm = 0.0
+      @day = 0.0
+      @emsq = 0.0
+      @gam = 0.0
+      @rtemsq = 0.0
+      @s1 = 0.0
+      @s2 = 0.0
+      @s3 = 0.0
+      @s4 = 0.0
+      @s5 = 0.0
+      @s6 = 0.0
+      @s7 = 0.0
+      @ss1 = 0.0
+      @ss2 = 0.0
+      @ss3 = 0.0
+      @ss4 = 0.0
+      @ss5 = 0.0
+      @ss6 = 0.0
+      @ss7 = 0.0
+      @sz1 = 0.0
+      @sz2 = 0.0
+      @sz3 = 0.0
+      @sz11 = 0.0
+      @sz12 = 0.0
+      @sz13 = 0.0
+      @sz21 = 0.0
+      @sz22 = 0.0
+      @sz23 = 0.0
+      @sz31 = 0.0
+      @sz32 = 0.0
+      @sz33 = 0.0
+      @z1 = 0.0
+      @z2 = 0.0
+      @z3 = 0.0
+      @z11 = 0.0
+      @z12 = 0.0
+      @z13 = 0.0
+      @z21 = 0.0
+      @z22 = 0.0
+      @z23 = 0.0
+      @z31 = 0.0
+      @z32 = 0.0
+      @z33 = 0.0
+      @argpm = 0.0
+      @inclm = 0.0
+      @nodem = 0.0
+      @dndt = 0.0
+      @eccsq = 0.0
+
+      # for initl
+      @ainv = 0.0
+      @ao = 0.0
+      @con42 = 0.0
+      @cosio = 0.0
+      @cosio2 = 0.0
+      @omeosq = 0.0
+      @posq = 0.0
+      @rp = 0.0
+      @rteosq = 0.0
+      @sinio = 0.0
+    end
+  end
+end