geodetic 0.5.0 → 0.5.2

data/examples/geodetic_demo.geojson

@@ -0,0 +1,305 @@
+{
+  "type": "FeatureCollection",
+  "features": [
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "Point",
+        "coordinates": [
+          -122.3493,
+          47.6205
+        ]
+      },
+      "properties": {
+        "name": "Seattle",
+        "type": "city"
+      }
+    },
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "Point",
+        "coordinates": [
+          -122.6784,
+          45.5152
+        ]
+      },
+      "properties": {
+        "name": "Portland",
+        "type": "city"
+      }
+    },
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "Point",
+        "coordinates": [
+          -122.4194,
+          37.7749
+        ]
+      },
+      "properties": {
+        "name": "San Francisco",
+        "type": "city"
+      }
+    },
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "Point",
+        "coordinates": [
+          -118.2437,
+          34.0522
+        ]
+      },
+      "properties": {
+        "name": "Los Angeles",
+        "type": "city"
+      }
+    },
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "Point",
+        "coordinates": [
+          -74.006,
+          40.7128
+        ]
+      },
+      "properties": {
+        "name": "New York",
+        "type": "city"
+      }
+    },
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "LineString",
+        "coordinates": [
+          [
+            -122.3493,
+            47.6205
+          ],
+          [
+            -122.6784,
+            45.5152
+          ],
+          [
+            -122.4194,
+            37.7749
+          ],
+          [
+            -118.2437,
+            34.0522
+          ]
+        ]
+      },
+      "properties": {
+        "name": "West Coast Route",
+        "mode": "driving"
+      }
+    },
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "Point",
+        "coordinates": [
+          2.3522,
+          48.8566
+        ]
+      },
+      "properties": {}
+    },
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "Point",
+        "coordinates": [
+          -0.1278,
+          51.5074
+        ]
+      },
+      "properties": {}
+    },
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "Polygon",
+        "coordinates": [
+          [
+            [
+              -123.0,
+              48.0
+            ],
+            [
+              -121.0,
+              48.0
+            ],
+            [
+              -121.0,
+              46.0
+            ],
+            [
+              -123.0,
+              46.0
+            ],
+            [
+              -123.0,
+              48.0
+            ]
+          ]
+        ]
+      },
+      "properties": {
+        "name": "Pacific NW",
+        "region": true
+      }
+    },
+    {
+      "type": "Feature",
+      "geometry": {
+        "type": "Polygon",
+        "coordinates": [
+          [
+            [
+              -122.3493,
+              47.710441151017214
+            ],
+            [
+              -122.32330338127136,
+              47.708710028492575
+            ],
+            [
+              -122.29831079292721,
+              47.70358352014371
+            ],
+            [
+              -122.27528691572327,
+              47.69525958316883
+            ],
+            [
+              -122.25511936150583,
+              47.684059519927274
+            ],
+            [
+              -122.23858413342101,
+              47.67041541322468
+            ],
+            [
+              -122.22631565848305,
+              47.654853266080394
+            ],
+            [
+              -122.21878256554382,
+              47.63797253321124
+            ],
+            [
+              -122.21627011382594,
+              47.62042286978415
+            ],
+            [
+              -122.21886987894307,
+              47.602879023582524
+            ],
+            [
+              -122.2264769927538,
+              47.586014856783315
+            ],
+            [
+              -122.23879492708923,
+              47.57047750255195
+            ],
+            [
+              -122.25534752329217,
+              47.5568626411414
+            ],
+            [
+              -122.27549771006404,
+              47.54569182325386
+            ],
+            [
+              -122.2984721281491,
+              47.53739267934147
+            ],
+            [
+              -122.32339069534318,
+              47.532282737212164
+            ],
+            [
+              -122.3493,
+              47.53055743197494
+            ],
+            [
+              -122.37520930465682,
+              47.532282737212164
+            ],
+            [
+              -122.4001278718509,
+              47.53739267934147
+            ],
+            [
+              -122.42310228993595,
+              47.54569182325386
+            ],
+            [
+              -122.44325247670783,
+              47.55686264114141
+            ],
+            [
+              -122.45980507291077,
+              47.57047750255195
+            ],
+            [
+              -122.4721230072462,
+              47.58601485678331
+            ],
+            [
+              -122.47973012105693,
+              47.60287902358252
+            ],
+            [
+              -122.48232988617406,
+              47.62042286978415
+            ],
+            [
+              -122.47981743445618,
+              47.63797253321123
+            ],
+            [
+              -122.47228434151695,
+              47.6548532660804
+            ],
+            [
+              -122.46001586657898,
+              47.670415413224674
+            ],
+            [
+              -122.44348063849417,
+              47.684059519927274
+            ],
+            [
+              -122.42331308427673,
+              47.69525958316883
+            ],
+            [
+              -122.40028920707279,
+              47.70358352014371
+            ],
+            [
+              -122.37529661872864,
+              47.708710028492575
+            ],
+            [
+              -122.3493,
+              47.710441151017214
+            ]
+          ]
+        ]
+      },
+      "properties": {
+        "name": "Seattle Metro",
+        "radius_km": 10
+      }
+    }
+  ]
+}

data/lib/geodetic/areas/regular_polygon.rb

@@ -31,29 +31,9 @@ module Geodetic
 
         n.times.map do |i|
           angle = @bearing + step * i
-          destination(@center, @radius, angle)
+          Vector.new(distance: @radius, bearing: angle).destination_from(@center)
         end
       end
-
-      # Compute a destination point given start, distance (meters), and bearing (degrees).
-      # Uses a simplified flat-earth projection adequate for polygon-scale distances.
-      def destination(origin, distance_m, bearing_deg)
-        bearing_rad = bearing_deg * Geodetic::RAD_PER_DEG
-        lat_rad = origin.lat * Geodetic::RAD_PER_DEG
-
-        # Approximate meters per degree at this latitude
-        m_per_deg_lat = 111_320.0
-        m_per_deg_lng = 111_320.0 * Math.cos(lat_rad)
-
-        dlat = distance_m * Math.cos(bearing_rad) / m_per_deg_lat
-        dlng = distance_m * Math.sin(bearing_rad) / m_per_deg_lng
-
-        Coordinate::LLA.new(
-          lat: origin.lat + dlat,
-          lng: origin.lng + dlng,
-          alt: origin.alt
-        )
-      end
     end
   end
 end

data/lib/geodetic/coordinate/bng.rb

@@ -138,15 +138,15 @@ module Geodetic
         end
       end
 
-      def self.from_lla(lla_coord, datum = WGS84)
+      def self.from_lla(lla, datum = WGS84)
         # Convert to OSGB36 if needed
         if datum == WGS84
           # Simplified WGS84 to OSGB36 transformation
-          lat = lla_coord.lat + 0.00015
-          lng = lla_coord.lng - 0.00045
-          osgb36_lla = LLA.new(lat: lat, lng: lng, alt: lla_coord.alt)
+          lat = lla.lat + 0.00015
+          lng = lla.lng - 0.00045
+          osgb36_lla = LLA.new(lat: lat, lng: lng, alt: lla.alt)
         else
-          osgb36_lla = lla_coord
+          osgb36_lla = lla
         end
 
         # Convert to BNG using Transverse Mercator
@@ -158,8 +158,8 @@ module Geodetic
         to_lla(datum).to_ecef(datum)
       end
 
-      def self.from_ecef(ecef_coord, datum = WGS84)
-        lla_coord = ecef_coord.to_lla(datum)
+      def self.from_ecef(ecef, datum = WGS84)
+        lla_coord = ecef.to_lla(datum)
         from_lla(lla_coord, datum)
       end
 
@@ -167,8 +167,8 @@ module Geodetic
         to_lla(datum).to_utm(datum)
       end
 
-      def self.from_utm(utm_coord, datum = WGS84)
-        lla_coord = utm_coord.to_lla(datum)
+      def self.from_utm(utm, datum = WGS84)
+        lla_coord = utm.to_lla(datum)
         from_lla(lla_coord, datum)
       end
 
@@ -176,8 +176,8 @@ module Geodetic
         to_lla(datum).to_enu(reference_lla)
       end
 
-      def self.from_enu(enu_coord, reference_lla, datum = WGS84)
-        lla_coord = enu_coord.to_lla(reference_lla)
+      def self.from_enu(enu, reference_lla, datum = WGS84)
+        lla_coord = enu.to_lla(reference_lla)
         from_lla(lla_coord, datum)
       end
 
@@ -185,8 +185,8 @@ module Geodetic
         to_lla(datum).to_ned(reference_lla)
       end
 
-      def self.from_ned(ned_coord, reference_lla, datum = WGS84)
-        lla_coord = ned_coord.to_lla(reference_lla)
+      def self.from_ned(ned, reference_lla, datum = WGS84)
+        lla_coord = ned.to_lla(reference_lla)
         from_lla(lla_coord, datum)
       end
 
@@ -194,8 +194,8 @@ module Geodetic
         MGRS.from_lla(to_lla(datum), datum, precision)
       end
 
-      def self.from_mgrs(mgrs_coord, datum = WGS84)
-        lla_coord = mgrs_coord.to_lla(datum)
+      def self.from_mgrs(mgrs, datum = WGS84)
+        lla_coord = mgrs.to_lla(datum)
         from_lla(lla_coord, datum)
       end
 
@@ -203,8 +203,8 @@ module Geodetic
         USNG.from_lla(to_lla(datum), datum, precision)
       end
 
-      def self.from_usng(usng_coord, datum = WGS84)
-        lla_coord = usng_coord.to_lla(datum)
+      def self.from_usng(usng, datum = WGS84)
+        lla_coord = usng.to_lla(datum)
         from_lla(lla_coord, datum)
       end
 
@@ -212,8 +212,8 @@ module Geodetic
         WebMercator.from_lla(to_lla(datum), datum)
       end
 
-      def self.from_web_mercator(web_mercator_coord, datum = WGS84)
-        lla_coord = web_mercator_coord.to_lla(datum)
+      def self.from_web_mercator(web_mercator, datum = WGS84)
+        lla_coord = web_mercator.to_lla(datum)
         from_lla(lla_coord, datum)
       end
 
@@ -221,8 +221,8 @@ module Geodetic
         UPS.from_lla(to_lla(datum), datum)
       end
 
-      def self.from_ups(ups_coord, datum = WGS84)
-        lla_coord = ups_coord.to_lla(datum)
+      def self.from_ups(ups, datum = WGS84)
+        lla_coord = ups.to_lla(datum)
         from_lla(lla_coord, datum)
       end
 
@@ -230,8 +230,8 @@ module Geodetic
         StatePlane.from_lla(to_lla(datum), zone_code, datum)
       end
 
-      def self.from_state_plane(state_plane_coord, datum = WGS84)
-        lla_coord = state_plane_coord.to_lla(datum)
+      def self.from_state_plane(state_plane, datum = WGS84)
+        lla_coord = state_plane.to_lla(datum)
         from_lla(lla_coord, datum)
       end
 
@@ -312,13 +312,13 @@ module Geodetic
         LLA.new(lat: lat_deg, lng: lng_deg, alt: 0.0)
       end
 
-      def self.transverse_mercator_forward(lla_coord, datum)
+      def self.transverse_mercator_forward(lla, datum)
         # Forward Transverse Mercator projection from OSGB36 lat/lng to BNG
         a = datum.a
         e2 = datum.e2
 
-        lat = lla_coord.lat * RAD_PER_DEG
-        lng = lla_coord.lng * RAD_PER_DEG
+        lat = lla.lat * RAD_PER_DEG
+        lng = lla.lng * RAD_PER_DEG
         lat0_rad = ORIGIN_LATITUDE * RAD_PER_DEG
         lng0_rad = ORIGIN_LONGITUDE * RAD_PER_DEG
         k0 = SCALE_FACTOR

data/lib/geodetic/coordinate/ecef.rb

@@ -140,48 +140,48 @@ module Geodetic
         MGRS.from_lla(to_lla(datum), datum, precision)
       end
 
-      def self.from_mgrs(mgrs_coord, datum = WGS84)
-        mgrs_coord.to_ecef(datum)
+      def self.from_mgrs(mgrs, datum = WGS84)
+        mgrs.to_ecef(datum)
       end
 
       def to_usng(datum = WGS84, precision = 5)
         USNG.from_lla(to_lla(datum), datum, precision)
       end
 
-      def self.from_usng(usng_coord, datum = WGS84)
-        usng_coord.to_ecef(datum)
+      def self.from_usng(usng, datum = WGS84)
+        usng.to_ecef(datum)
       end
 
       def to_web_mercator(datum = WGS84)
         WebMercator.from_lla(to_lla(datum), datum)
       end
 
-      def self.from_web_mercator(wm_coord, datum = WGS84)
-        wm_coord.to_ecef(datum)
+      def self.from_web_mercator(web_mercator, datum = WGS84)
+        web_mercator.to_ecef(datum)
       end
 
       def to_ups(datum = WGS84)
         UPS.from_lla(to_lla(datum), datum)
       end
 
-      def self.from_ups(ups_coord, datum = WGS84)
-        ups_coord.to_ecef(datum)
+      def self.from_ups(ups, datum = WGS84)
+        ups.to_ecef(datum)
       end
 
       def to_state_plane(zone_code, datum = WGS84)
         StatePlane.from_lla(to_lla(datum), zone_code, datum)
       end
 
-      def self.from_state_plane(sp_coord, datum = WGS84)
-        sp_coord.to_ecef(datum)
+      def self.from_state_plane(state_plane, datum = WGS84)
+        state_plane.to_ecef(datum)
       end
 
       def to_bng(datum = WGS84)
         BNG.from_lla(to_lla(datum))
       end
 
-      def self.from_bng(bng_coord, datum = WGS84)
-        bng_coord.to_ecef(datum)
+      def self.from_bng(bng, datum = WGS84)
+        bng.to_ecef(datum)
       end
 
       def to_s(precision = 2)

data/lib/geodetic/coordinate/enu.rb

@@ -108,8 +108,8 @@ module Geodetic
         MGRS.from_lla(to_lla(reference_lla), datum, precision)
       end
 
-      def self.from_mgrs(mgrs_coord, reference_lla, datum = WGS84)
-        lla = mgrs_coord.to_lla(datum)
+      def self.from_mgrs(mgrs, reference_lla, datum = WGS84)
+        lla = mgrs.to_lla(datum)
         from_lla(lla, reference_lla)
       end
 
@@ -117,8 +117,8 @@ module Geodetic
         USNG.from_lla(to_lla(reference_lla), datum, precision)
       end
 
-      def self.from_usng(usng_coord, reference_lla, datum = WGS84)
-        lla = usng_coord.to_lla(datum)
+      def self.from_usng(usng, reference_lla, datum = WGS84)
+        lla = usng.to_lla(datum)
         from_lla(lla, reference_lla)
       end
 
@@ -126,8 +126,8 @@ module Geodetic
         WebMercator.from_lla(to_lla(reference_lla), datum)
       end
 
-      def self.from_web_mercator(wm_coord, reference_lla, datum = WGS84)
-        lla = wm_coord.to_lla(datum)
+      def self.from_web_mercator(web_mercator, reference_lla, datum = WGS84)
+        lla = web_mercator.to_lla(datum)
         from_lla(lla, reference_lla)
       end
 
@@ -135,8 +135,8 @@ module Geodetic
         UPS.from_lla(to_lla(reference_lla), datum)
       end
 
-      def self.from_ups(ups_coord, reference_lla, datum = WGS84)
-        lla = ups_coord.to_lla(datum)
+      def self.from_ups(ups, reference_lla, datum = WGS84)
+        lla = ups.to_lla(datum)
         from_lla(lla, reference_lla)
       end
 
@@ -144,8 +144,8 @@ module Geodetic
         StatePlane.from_lla(to_lla(reference_lla), zone_code, datum)
       end
 
-      def self.from_state_plane(sp_coord, reference_lla, datum = WGS84)
-        lla = sp_coord.to_lla(datum)
+      def self.from_state_plane(state_plane, reference_lla, datum = WGS84)
+        lla = state_plane.to_lla(datum)
         from_lla(lla, reference_lla)
       end
 
@@ -153,8 +153,8 @@ module Geodetic
         BNG.from_lla(to_lla(reference_lla))
       end
 
-      def self.from_bng(bng_coord, reference_lla)
-        lla = bng_coord.to_lla
+      def self.from_bng(bng, reference_lla)
+        lla = bng.to_lla
         from_lla(lla, reference_lla)
       end
 
@@ -162,8 +162,8 @@ module Geodetic
         GH36.new(to_lla(reference_lla), precision: precision)
       end
 
-      def self.from_gh36(gh36_coord, reference_lla)
-        lla = gh36_coord.to_lla
+      def self.from_gh36(gh36, reference_lla)
+        lla = gh36.to_lla
         from_lla(lla, reference_lla)
       end
 
@@ -171,8 +171,8 @@ module Geodetic
         GH.new(to_lla(reference_lla), precision: precision)
       end
 
-      def self.from_gh(gh_coord, reference_lla)
-        lla = gh_coord.to_lla
+      def self.from_gh(gh, reference_lla)
+        lla = gh.to_lla
         from_lla(lla, reference_lla)
       end
 
@@ -180,8 +180,8 @@ module Geodetic
         HAM.new(to_lla(reference_lla), precision: precision)
       end
 
-      def self.from_ham(ham_coord, reference_lla)
-        lla = ham_coord.to_lla
+      def self.from_ham(ham, reference_lla)
+        lla = ham.to_lla
         from_lla(lla, reference_lla)
       end
 
@@ -189,8 +189,8 @@ module Geodetic
         OLC.new(to_lla(reference_lla), precision: precision)
       end
 
-      def self.from_olc(olc_coord, reference_lla)
-        lla = olc_coord.to_lla
+      def self.from_olc(olc, reference_lla)
+        lla = olc.to_lla
         from_lla(lla, reference_lla)
       end
 
@@ -198,8 +198,8 @@ module Geodetic
         GEOREF.new(to_lla(reference_lla), precision: precision)
       end
 
-      def self.from_georef(georef_coord, reference_lla)
-        lla = georef_coord.to_lla
+      def self.from_georef(georef, reference_lla)
+        lla = georef.to_lla
         from_lla(lla, reference_lla)
       end
 
@@ -207,8 +207,8 @@ module Geodetic
         GARS.new(to_lla(reference_lla), precision: precision)
       end
 
-      def self.from_gars(gars_coord, reference_lla)
-        lla = gars_coord.to_lla
+      def self.from_gars(gars, reference_lla)
+        lla = gars.to_lla
         from_lla(lla, reference_lla)
       end
 
@@ -216,8 +216,8 @@ module Geodetic
         H3.new(to_lla(reference_lla), precision: precision)
       end
 
-      def self.from_h3(h3_coord, reference_lla)
-        lla = h3_coord.to_lla
+      def self.from_h3(h3, reference_lla)
+        lla = h3.to_lla
         from_lla(lla, reference_lla)
       end