geodetic 0.1.0 → 0.3.0

data/docs/reference/serialization.md

@@ -28,7 +28,7 @@ All `to_s` methods accept an optional `precision` parameter controlling the numb
 Passing `0` returns integer values (no decimal point). MGRS and USNG are string-based and do not accept a precision parameter.
 
 ```ruby
-lla = Geodetic::Coordinates::LLA.new(lat: 47.6205, lng: -122.3493, alt: 184.0)
+lla = Geodetic::Coordinate::LLA.new(lat: 47.6205, lng: -122.3493, alt: 184.0)
 lla.to_s        # => "47.620500, -122.349300, 184.00"
 lla.to_s(3)     # => "47.620, -122.349, 184.00"
 lla.to_s(0)     # => "48, -122, 184"
@@ -41,14 +41,14 @@ lla.to_s(0)     # => "48, -122, 184"
 ### LLA
 
 ```ruby
-point = Geodetic::Coordinates::LLA.new(lat: 38.8977, lng: -77.0365, alt: 100.0)
+point = Geodetic::Coordinate::LLA.new(lat: 38.8977, lng: -77.0365, alt: 100.0)
 
 point.to_s    # => "38.897700, -77.036500, 100.00"
 point.to_s(2) # => "38.90, -77.04, 100.00"
 point.to_a    # => [38.8977, -77.0365, 100.0]
 
-Geodetic::Coordinates::LLA.from_string("38.8977, -77.0365, 100.0")
-Geodetic::Coordinates::LLA.from_array([38.8977, -77.0365, 100.0])
+Geodetic::Coordinate::LLA.from_string("38.8977, -77.0365, 100.0")
+Geodetic::Coordinate::LLA.from_array([38.8977, -77.0365, 100.0])
 ```
 
 **LLA-specific: DMS format**
@@ -59,7 +59,7 @@ LLA also supports degrees-minutes-seconds notation:
 point.to_dms
 # => "38 53' 51.72\" N, 77 2' 11.40\" W, 100.00 m"
 
-Geodetic::Coordinates::LLA.from_dms("38 53' 51.72\" N, 77 2' 11.40\" W, 100.0 m")
+Geodetic::Coordinate::LLA.from_dms("38 53' 51.72\" N, 77 2' 11.40\" W, 100.0 m")
 ```
 
 The DMS string format is: `DD MM' SS.ss" H, DDD MM' SS.ss" H, ALT m` where H is N/S for latitude and E/W for longitude. The altitude portion is optional in `from_dms` (defaults to 0.0).
@@ -67,26 +67,26 @@ The DMS string format is: `DD MM' SS.ss" H, DDD MM' SS.ss" H, ALT m` where H is
 ### ECEF
 
 ```ruby
-point = Geodetic::Coordinates::ECEF.new(x: 1130730.0, y: -4828583.0, z: 3991570.0)
+point = Geodetic::Coordinate::ECEF.new(x: 1130730.0, y: -4828583.0, z: 3991570.0)
 
 point.to_s    # => "1130730.00, -4828583.00, 3991570.00"
 point.to_s(0) # => "1130730, -4828583, 3991570"
 point.to_a    # => [1130730.0, -4828583.0, 3991570.0]
 
-Geodetic::Coordinates::ECEF.from_string("1130730.0, -4828583.0, 3991570.0")
-Geodetic::Coordinates::ECEF.from_array([1130730.0, -4828583.0, 3991570.0])
+Geodetic::Coordinate::ECEF.from_string("1130730.0, -4828583.0, 3991570.0")
+Geodetic::Coordinate::ECEF.from_array([1130730.0, -4828583.0, 3991570.0])
 ```
 
 ### UTM
 
 ```ruby
-point = Geodetic::Coordinates::UTM.new(easting: 323394.0, northing: 4307396.0, altitude: 100.0, zone: 18, hemisphere: 'N')
+point = Geodetic::Coordinate::UTM.new(easting: 323394.0, northing: 4307396.0, altitude: 100.0, zone: 18, hemisphere: 'N')
 
 point.to_s    # => "323394.00, 4307396.00, 100.00, 18, N"
 point.to_a    # => [323394.0, 4307396.0, 100.0, 18, "N"]
 
-Geodetic::Coordinates::UTM.from_string("323394.0, 4307396.0, 100.0, 18, N")
-Geodetic::Coordinates::UTM.from_array([323394.0, 4307396.0, 100.0, 18, "N"])
+Geodetic::Coordinate::UTM.from_string("323394.0, 4307396.0, 100.0, 18, N")
+Geodetic::Coordinate::UTM.from_array([323394.0, 4307396.0, 100.0, 18, "N"])
 ```
 
 Note: The array and string formats include all five components: easting, northing, altitude, zone number, and hemisphere.
@@ -94,67 +94,67 @@ Note: The array and string formats include all five components: easting, northin
 ### ENU
 
 ```ruby
-point = Geodetic::Coordinates::ENU.new(e: 100.0, n: 200.0, u: 50.0)
+point = Geodetic::Coordinate::ENU.new(e: 100.0, n: 200.0, u: 50.0)
 
 point.to_s    # => "100.00, 200.00, 50.00"
 point.to_a    # => [100.0, 200.0, 50.0]
 
-Geodetic::Coordinates::ENU.from_string("100.0, 200.0, 50.0")
-Geodetic::Coordinates::ENU.from_array([100.0, 200.0, 50.0])
+Geodetic::Coordinate::ENU.from_string("100.0, 200.0, 50.0")
+Geodetic::Coordinate::ENU.from_array([100.0, 200.0, 50.0])
 ```
 
 ### NED
 
 ```ruby
-point = Geodetic::Coordinates::NED.new(n: 200.0, e: 100.0, d: -50.0)
+point = Geodetic::Coordinate::NED.new(n: 200.0, e: 100.0, d: -50.0)
 
 point.to_s    # => "200.00, 100.00, -50.00"
 point.to_a    # => [200.0, 100.0, -50.0]
 
-Geodetic::Coordinates::NED.from_string("200.0, 100.0, -50.0")
-Geodetic::Coordinates::NED.from_array([200.0, 100.0, -50.0])
+Geodetic::Coordinate::NED.from_string("200.0, 100.0, -50.0")
+Geodetic::Coordinate::NED.from_array([200.0, 100.0, -50.0])
 ```
 
 ### WebMercator
 
 ```ruby
-point = Geodetic::Coordinates::WebMercator.new(x: -8575605.0, y: 4707175.0)
+point = Geodetic::Coordinate::WebMercator.new(x: -8575605.0, y: 4707175.0)
 
 point.to_s    # => "-8575605.00, 4707175.00"
 point.to_a    # => [-8575605.0, 4707175.0]
 
-Geodetic::Coordinates::WebMercator.from_string("-8575605.0, 4707175.0")
-Geodetic::Coordinates::WebMercator.from_array([-8575605.0, 4707175.0])
+Geodetic::Coordinate::WebMercator.from_string("-8575605.0, 4707175.0")
+Geodetic::Coordinate::WebMercator.from_array([-8575605.0, 4707175.0])
 ```
 
 ### UPS
 
 ```ruby
-point = Geodetic::Coordinates::UPS.new(easting: 2000000.0, northing: 2000000.0, hemisphere: 'N', zone: 'Y')
+point = Geodetic::Coordinate::UPS.new(easting: 2000000.0, northing: 2000000.0, hemisphere: 'N', zone: 'Y')
 
 point.to_s    # => "2000000.00, 2000000.00, N, Y"
 point.to_a    # => [2000000.0, 2000000.0, "N", "Y"]
 
-Geodetic::Coordinates::UPS.from_string("2000000.0, 2000000.0, N, Y")
-Geodetic::Coordinates::UPS.from_array([2000000.0, 2000000.0, "N", "Y"])
+Geodetic::Coordinate::UPS.from_string("2000000.0, 2000000.0, N, Y")
+Geodetic::Coordinate::UPS.from_array([2000000.0, 2000000.0, "N", "Y"])
 ```
 
 ### BNG
 
 ```ruby
-point = Geodetic::Coordinates::BNG.new(easting: 530000.0, northing: 180000.0)
+point = Geodetic::Coordinate::BNG.new(easting: 530000.0, northing: 180000.0)
 
 point.to_s    # => "530000.00, 180000.00"
 point.to_a    # => [530000.0, 180000.0]
 
-Geodetic::Coordinates::BNG.from_string("530000.0, 180000.0")
-Geodetic::Coordinates::BNG.from_array([530000.0, 180000.0])
+Geodetic::Coordinate::BNG.from_string("530000.0, 180000.0")
+Geodetic::Coordinate::BNG.from_array([530000.0, 180000.0])
 ```
 
 BNG also supports grid reference notation via the constructor:
 
 ```ruby
-point = Geodetic::Coordinates::BNG.new(grid_ref: "TQ 30 80")
+point = Geodetic::Coordinate::BNG.new(grid_ref: "TQ 30 80")
 point.to_grid_reference(6)  # => "TQ 300000 800000"
 point.to_grid_reference(0)  # => "TQ" (grid square only)
 ```
@@ -162,13 +162,13 @@ point.to_grid_reference(0)  # => "TQ" (grid square only)
 ### StatePlane
 
 ```ruby
-point = Geodetic::Coordinates::StatePlane.new(easting: 2000000.0, northing: 500000.0, zone_code: 'CA_I')
+point = Geodetic::Coordinate::StatePlane.new(easting: 2000000.0, northing: 500000.0, zone_code: 'CA_I')
 
 point.to_s    # => "2000000.00, 500000.00, CA_I"
 point.to_a    # => [2000000.0, 500000.0, "CA_I"]
 
-Geodetic::Coordinates::StatePlane.from_string("2000000.0, 500000.0, CA_I")
-Geodetic::Coordinates::StatePlane.from_array([2000000.0, 500000.0, "CA_I"])
+Geodetic::Coordinate::StatePlane.from_string("2000000.0, 500000.0, CA_I")
+Geodetic::Coordinate::StatePlane.from_array([2000000.0, 500000.0, "CA_I"])
 ```
 
 ---
@@ -180,12 +180,12 @@ MGRS and USNG use alphanumeric grid references rather than numeric arrays. They
 ### MGRS
 
 ```ruby
-mgrs = Geodetic::Coordinates::MGRS.new(mgrs_string: "18SUJ2034706880")
+mgrs = Geodetic::Coordinate::MGRS.new(mgrs_string: "18SUJ2034706880")
 mgrs.to_s    # => "18SUJ2034706880"
 mgrs.to_a    # => ["18S", "UJ", 20347.0, 6880.0, 5]
 
-Geodetic::Coordinates::MGRS.from_string("18SUJ2034706880")
-Geodetic::Coordinates::MGRS.from_array(["18S", "UJ", 20347, 6880, 5])
+Geodetic::Coordinate::MGRS.from_string("18SUJ2034706880")
+Geodetic::Coordinate::MGRS.from_array(["18S", "UJ", 20347, 6880, 5])
 ```
 
 The MGRS string format is: `{zone_number}{zone_letter}{square_id}{easting}{northing}` with no spaces. Precision varies based on the number of coordinate digits (0 to 5 pairs). The array format is `[grid_zone, square_id, easting, northing, precision]`.
@@ -193,12 +193,12 @@ The MGRS string format is: `{zone_number}{zone_letter}{square_id}{easting}{north
 ### USNG
 
 ```ruby
-usng = Geodetic::Coordinates::USNG.new(usng_string: "18S UJ 20347 06880")
+usng = Geodetic::Coordinate::USNG.new(usng_string: "18S UJ 20347 06880")
 usng.to_s    # => "18S UJ 20347 06880"
 usng.to_a    # => ["18S", "UJ", 20347.0, 6880.0, 5]
 
-Geodetic::Coordinates::USNG.from_string("18S UJ 20347 06880")
-Geodetic::Coordinates::USNG.from_array(["18S", "UJ", 20347, 6880, 5])
+Geodetic::Coordinate::USNG.from_string("18S UJ 20347 06880")
+Geodetic::Coordinate::USNG.from_array(["18S", "UJ", 20347, 6880, 5])
 ```
 
 USNG uses the same underlying format as MGRS but separates components with spaces for readability. Both spaced and non-spaced formats are accepted by `from_string`. The array format is `[grid_zone, square_id, easting, northing, precision]`. USNG also provides:
@@ -216,24 +216,24 @@ String and array serialization support full roundtrip fidelity:
 
 ```ruby
 # LLA roundtrip via string
-original = Geodetic::Coordinates::LLA.new(lat: 40.7128, lng: -74.0060, alt: 10.0)
-restored = Geodetic::Coordinates::LLA.from_string(original.to_s)
+original = Geodetic::Coordinate::LLA.new(lat: 40.7128, lng: -74.0060, alt: 10.0)
+restored = Geodetic::Coordinate::LLA.from_string(original.to_s)
 original == restored  # => true
 
 # ECEF roundtrip via array
-original = Geodetic::Coordinates::ECEF.new(x: 1334000.0, y: -4654000.0, z: 4138000.0)
-restored = Geodetic::Coordinates::ECEF.from_array(original.to_a)
+original = Geodetic::Coordinate::ECEF.new(x: 1334000.0, y: -4654000.0, z: 4138000.0)
+restored = Geodetic::Coordinate::ECEF.from_array(original.to_a)
 original == restored  # => true
 
 # UTM roundtrip via string
-original = Geodetic::Coordinates::UTM.new(easting: 583960.0, northing: 4507523.0, altitude: 10.0, zone: 18, hemisphere: 'N')
-restored = Geodetic::Coordinates::UTM.from_string(original.to_s)
+original = Geodetic::Coordinate::UTM.new(easting: 583960.0, northing: 4507523.0, altitude: 10.0, zone: 18, hemisphere: 'N')
+restored = Geodetic::Coordinate::UTM.from_string(original.to_s)
 original == restored  # => true
 
 # LLA roundtrip via DMS
-original = Geodetic::Coordinates::LLA.new(lat: 38.8977, lng: -77.0365, alt: 100.0)
+original = Geodetic::Coordinate::LLA.new(lat: 38.8977, lng: -77.0365, alt: 100.0)
 dms_string = original.to_dms
-restored = Geodetic::Coordinates::LLA.from_dms(dms_string)
+restored = Geodetic::Coordinate::LLA.from_dms(dms_string)
 # Note: DMS roundtrip has minor floating-point precision differences
 ```
 

data/examples/01_basic_conversions.rb

@@ -4,18 +4,18 @@
 # Shows all coordinate systems converting to/from each other
 
 require_relative '../lib/geodetic'
-require_relative '../lib/geodetic/coordinates/lla'
-require_relative '../lib/geodetic/coordinates/ecef'
-require_relative '../lib/geodetic/coordinates/enu'
-require_relative '../lib/geodetic/coordinates/ned'
-require_relative '../lib/geodetic/coordinates/utm'
+require_relative '../lib/geodetic/coordinate/lla'
+require_relative '../lib/geodetic/coordinate/ecef'
+require_relative '../lib/geodetic/coordinate/enu'
+require_relative '../lib/geodetic/coordinate/ned'
+require_relative '../lib/geodetic/coordinate/utm'
 
 include Geodetic
-LLA  = Coordinates::LLA
-ECEF = Coordinates::ECEF
-ENU  = Coordinates::ENU
-NED  = Coordinates::NED
-UTM  = Coordinates::UTM
+LLA  = Coordinate::LLA
+ECEF = Coordinate::ECEF
+ENU  = Coordinate::ENU
+NED  = Coordinate::NED
+UTM  = Coordinate::UTM
 
 puts "=== Orthogonal Coordinate System Conversions Demo ==="
 puts

data/examples/02_all_coordinate_systems.rb

@@ -3,33 +3,33 @@
 # Shows complete orthogonal conversions between all implemented coordinate systems
 
 require_relative '../lib/geodetic'
-require_relative '../lib/geodetic/coordinates/lla'
-require_relative '../lib/geodetic/coordinates/ecef'
-require_relative '../lib/geodetic/coordinates/utm'
-require_relative '../lib/geodetic/coordinates/enu'
-require_relative '../lib/geodetic/coordinates/ned'
-require_relative '../lib/geodetic/coordinates/mgrs'
-require_relative '../lib/geodetic/coordinates/web_mercator'
-require_relative '../lib/geodetic/coordinates/ups'
-require_relative '../lib/geodetic/coordinates/usng'
-require_relative '../lib/geodetic/coordinates/state_plane'
-require_relative '../lib/geodetic/coordinates/bng'
-require_relative '../lib/geodetic/coordinates/gh36'
+require_relative '../lib/geodetic/coordinate/lla'
+require_relative '../lib/geodetic/coordinate/ecef'
+require_relative '../lib/geodetic/coordinate/utm'
+require_relative '../lib/geodetic/coordinate/enu'
+require_relative '../lib/geodetic/coordinate/ned'
+require_relative '../lib/geodetic/coordinate/mgrs'
+require_relative '../lib/geodetic/coordinate/web_mercator'
+require_relative '../lib/geodetic/coordinate/ups'
+require_relative '../lib/geodetic/coordinate/usng'
+require_relative '../lib/geodetic/coordinate/state_plane'
+require_relative '../lib/geodetic/coordinate/bng'
+require_relative '../lib/geodetic/coordinate/gh36'
 require_relative '../lib/geodetic/geoid_height'
 
 include Geodetic
-LLA        = Coordinates::LLA
-ECEF       = Coordinates::ECEF
-UTM_Coord  = Coordinates::UTM
-ENU        = Coordinates::ENU
-NED        = Coordinates::NED
-MGRS       = Coordinates::MGRS
-USNG       = Coordinates::USNG
-WebMerc    = Coordinates::WebMercator
-UPS        = Coordinates::UPS
-StatePlane = Coordinates::StatePlane
-BNG        = Coordinates::BNG
-GH36       = Coordinates::GH36
+LLA        = Coordinate::LLA
+ECEF       = Coordinate::ECEF
+UTM_Coord  = Coordinate::UTM
+ENU        = Coordinate::ENU
+NED        = Coordinate::NED
+MGRS       = Coordinate::MGRS
+USNG       = Coordinate::USNG
+WebMerc    = Coordinate::WebMercator
+UPS        = Coordinate::UPS
+StatePlane = Coordinate::StatePlane
+BNG        = Coordinate::BNG
+GH36       = Coordinate::GH36
 
 def demo_coordinate_systems
   puts "=" * 80

data/lib/geodetic/areas/circle.rb

@@ -1,6 +1,6 @@
 # frozen_string_literal: true
 
-require_relative '../coordinates/lla'
+require_relative '../coordinate/lla'
 
 module Geodetic
   module Areas

data/lib/geodetic/areas/polygon.rb

@@ -1,6 +1,6 @@
 # frozen_string_literal: true
 
-require_relative '../coordinates/lla'
+require_relative '../coordinate/lla'
 
 module Geodetic
   module Areas
@@ -27,7 +27,7 @@ module Geodetic
         centroid_lng /= (6.0 * area)
         centroid_lat /= (6.0 * area)
 
-        @centroid = Coordinates::LLA.new(lat: centroid_lat, lng: centroid_lng, alt: 0.0)
+        @centroid = Coordinate::LLA.new(lat: centroid_lat, lng: centroid_lng, alt: 0.0)
       end
 
       def includes?(a_point)

data/lib/geodetic/areas/rectangle.rb

@@ -1,6 +1,6 @@
 # frozen_string_literal: true
 
-require_relative '../coordinates/lla'
+require_relative '../coordinate/lla'
 
 module Geodetic
   module Areas
@@ -15,13 +15,13 @@ module Geodetic
       #     se: LLA.new(lat: 40.0, lng: -74.0)
       #   )
       def initialize(nw:, se:)
-        @nw = nw.is_a?(Coordinates::LLA) ? nw : nw.to_lla
-        @se = se.is_a?(Coordinates::LLA) ? se : se.to_lla
+        @nw = nw.is_a?(Coordinate::LLA) ? nw : nw.to_lla
+        @se = se.is_a?(Coordinate::LLA) ? se : se.to_lla
 
         raise ArgumentError, "NW corner must have higher latitude than SE corner" if @nw.lat < @se.lat
         raise ArgumentError, "NW corner must have lower longitude than SE corner" if @nw.lng > @se.lng
 
-        @centroid = Coordinates::LLA.new(
+        @centroid = Coordinate::LLA.new(
           lat: (@nw.lat + @se.lat) / 2.0,
           lng: (@nw.lng + @se.lng) / 2.0,
           alt: 0.0
@@ -29,11 +29,11 @@ module Geodetic
       end
 
       def ne
-        Coordinates::LLA.new(lat: @nw.lat, lng: @se.lng, alt: 0.0)
+        Coordinate::LLA.new(lat: @nw.lat, lng: @se.lng, alt: 0.0)
       end
 
       def sw
-        Coordinates::LLA.new(lat: @se.lat, lng: @nw.lng, alt: 0.0)
+        Coordinate::LLA.new(lat: @se.lat, lng: @nw.lng, alt: 0.0)
       end
 
       def includes?(a_point)

data/lib/geodetic/{coordinates → coordinate}/bng.rb

@@ -5,7 +5,7 @@
 # Uses Transverse Mercator projection with specific parameters
 
 module Geodetic
-  module Coordinates
+  module Coordinate
     class BNG
       require_relative '../datum'
 
@@ -228,42 +228,6 @@ module Geodetic
         from_lla(lla_coord, datum)
       end
 
-      def to_gh36(datum = WGS84, precision: 10)
-        GH36.new(to_lla(datum), precision: precision)
-      end
-
-      def self.from_gh36(gh36_coord, datum = WGS84)
-        lla_coord = gh36_coord.to_lla(datum)
-        from_lla(lla_coord, datum)
-      end
-
-      def to_gh(datum = WGS84, precision: 12)
-        GH.new(to_lla(datum), precision: precision)
-      end
-
-      def self.from_gh(gh_coord, datum = WGS84)
-        lla_coord = gh_coord.to_lla(datum)
-        from_lla(lla_coord, datum)
-      end
-
-      def to_ham(datum = WGS84, precision: 6)
-        HAM.new(to_lla(datum), precision: precision)
-      end
-
-      def self.from_ham(ham_coord, datum = WGS84)
-        lla_coord = ham_coord.to_lla(datum)
-        from_lla(lla_coord, datum)
-      end
-
-      def to_olc(datum = WGS84, precision: 10)
-        OLC.new(to_lla(datum), precision: precision)
-      end
-
-      def self.from_olc(olc_coord, datum = WGS84)
-        lla_coord = olc_coord.to_lla(datum)
-        from_lla(lla_coord, datum)
-      end
-
       def ==(other)
         return false unless other.is_a?(BNG)
 
@@ -388,6 +352,8 @@ module Geodetic
 
         new(easting: easting, northing: northing)
       end
+
+      Coordinate.register_class(self)
     end
   end
 end

data/lib/geodetic/{coordinates → coordinate}/ecef.rb

@@ -3,7 +3,7 @@
 require_relative '../datum'
 
 module Geodetic
-  module Coordinates
+  module Coordinate
     class ECEF
       attr_reader :x, :y, :z
 
@@ -184,38 +184,6 @@ module Geodetic
         bng_coord.to_ecef(datum)
       end
 
-      def to_gh36(precision: 10)
-        GH36.new(to_lla, precision: precision)
-      end
-
-      def self.from_gh36(gh36_coord, datum = WGS84)
-        gh36_coord.to_ecef(datum)
-      end
-
-      def to_gh(precision: 12)
-        GH.new(to_lla, precision: precision)
-      end
-
-      def self.from_gh(gh_coord, datum = WGS84)
-        gh_coord.to_ecef(datum)
-      end
-
-      def to_ham(precision: 6)
-        HAM.new(to_lla, precision: precision)
-      end
-
-      def self.from_ham(ham_coord, datum = WGS84)
-        ham_coord.to_ecef(datum)
-      end
-
-      def to_olc(precision: 10)
-        OLC.new(to_lla, precision: precision)
-      end
-
-      def self.from_olc(olc_coord, datum = WGS84)
-        olc_coord.to_ecef(datum)
-      end
-
       def to_s(precision = 2)
         precision = precision.to_i
         if precision == 0
@@ -248,6 +216,8 @@ module Geodetic
         delta_x <= 1e-6 && delta_y <= 1e-6 && delta_z <= 1e-6
       end
 
+
+      Coordinate.register_class(self)
     end
   end
 end

data/lib/geodetic/{coordinates → coordinate}/enu.rb

@@ -3,7 +3,7 @@
 require_relative '../datum'
 
 module Geodetic
-  module Coordinates
+  module Coordinate
     class ENU
       attr_reader :e, :n, :u
       alias_method :east, :e
@@ -194,6 +194,33 @@ module Geodetic
         from_lla(lla, reference_lla)
       end
 
+      def to_georef(reference_lla, precision: 8)
+        GEOREF.new(to_lla(reference_lla), precision: precision)
+      end
+
+      def self.from_georef(georef_coord, reference_lla)
+        lla = georef_coord.to_lla
+        from_lla(lla, reference_lla)
+      end
+
+      def to_gars(reference_lla, precision: 7)
+        GARS.new(to_lla(reference_lla), precision: precision)
+      end
+
+      def self.from_gars(gars_coord, reference_lla)
+        lla = gars_coord.to_lla
+        from_lla(lla, reference_lla)
+      end
+
+      def to_h3(reference_lla, precision: 7)
+        H3.new(to_lla(reference_lla), precision: precision)
+      end
+
+      def self.from_h3(h3_coord, reference_lla)
+        lla = h3_coord.to_lla
+        from_lla(lla, reference_lla)
+      end
+
       def to_s(precision = 2)
         precision = precision.to_i
         if precision == 0
@@ -266,6 +293,8 @@ module Geodetic
       def horizontal_distance_to_origin
         Math.sqrt(@e**2 + @n**2)
       end
+
+      Coordinate.register_class(self)
     end
   end
 end