turf-ruby 0.2.0 → 0.3.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: f8ed0e8a70920277e4e8507100d1f14084e75d0904c03feea99dc8196530561a
-  data.tar.gz: baba690e3110b85c5ece7eed4f5bd3229b5ecfbbeeca7ad4a7ca8e05d182763f
+  metadata.gz: b858fdec044cd9bcb71a29727849dee8ff75f74cf4bf2f97163fafb08be006ab
+  data.tar.gz: 718bd40cab7ece92b2f95d080d7596bfbac2de785b3ba129078c01a2085cc252
 SHA512:
-  metadata.gz: 9e331328f338fff791262493ae0f63da0baea68b5ec457c911ce62711803fc5f91b989bb5a2abc5191ff354f25306c664cb1b00a75c857484da518da4367a475
-  data.tar.gz: 7d1814173daf223ec7d8f9bbd4ae080c18c85bcfe4c0d307ada85abde67772aaa5d58bf87024fb5757825b5faa7544ab908ff12bc67c3cc5302ecf5e5281a0d6
+  metadata.gz: c115975844db3097b16c9ba05b601104a9f8de6ff4216797dab70bfd539b282493d9988835074ed2ba36f43ae806f8d101278e021d4615dcb7153c8cadec95fc
+  data.tar.gz: f155f788ea240f39d93603f6cd33e6712ccd7acfb71ca63d6052ad3f1df9d14d2b52d32ec8395f24ee6f922920dca57961d874d9fc41ebc566d4d3fc340c43fe

data/.gitignore

@@ -6,3 +6,4 @@
 /pkg/
 /spec/reports/
 /tmp/
+.ruby-version

data/Gemfile

@@ -6,4 +6,5 @@ source "https://rubygems.org"
 gemspec
 
 gem "minitest", "~> 5.0"
+gem "minitest-focus"
 gem "rake", "~> 12.0"

data/Gemfile.lock

@@ -1,12 +1,14 @@
 PATH
   remote: .
   specs:
-    turf-ruby (0.2.0)
+    turf-ruby (0.3.0)
 
 GEM
   remote: https://rubygems.org/
   specs:
     minitest (5.14.0)
+    minitest-focus (1.2.1)
+      minitest (>= 4, < 6)
     rake (12.3.3)
 
 PLATFORMS
@@ -14,6 +16,7 @@ PLATFORMS
 
 DEPENDENCIES
   minitest (~> 5.0)
+  minitest-focus
   rake (~> 12.0)
   turf-ruby!
 

data/lib/turf/along.rb

@@ -0,0 +1,27 @@
+# frozen_string_literal: true
+
+module Turf
+  def self.along(line, distance, **options)
+    travelled = 0
+
+    geom = get_geom line
+    coords = geom[:coordinates]
+
+    coords.each_with_index do |coord, i|
+      break if distance >= travelled && i == coords.length - 1
+
+      if travelled >= distance
+        overshot = distance - travelled
+        return point(coord) if overshot.zero?
+
+        direction = bearing(coord, coords[i - 1]) - 180
+        interpolated = destination(coord, overshot, direction, options)
+        return interpolated
+      else
+        travelled += distance(coords[i], coords[i + 1], options)
+      end
+    end
+
+    point(coords.last)
+  end
+end

data/lib/turf/bearing.rb

@@ -0,0 +1,25 @@
+# frozen_string_literal: true
+
+module Turf
+  def self.bearing(from, to, **options)
+    return calculate_final_bearing(from, to) if options[:final] == true
+
+    coordinates1 = get_coord from
+    coordinates2 = get_coord to
+
+    lon1 = degrees_to_radians(coordinates1[0])
+    lon2 = degrees_to_radians(coordinates2[0])
+    lat1 = degrees_to_radians(coordinates1[1])
+    lat2 = degrees_to_radians(coordinates2[1])
+    a = Math.sin(lon2 - lon1) * Math.cos(lat2)
+    b = Math.cos(lat1) * Math.sin(lat2) -
+        Math.sin(lat1) * Math.cos(lat2) * Math.cos(lon2 - lon1)
+
+    radians_to_degrees(Math.atan2(a, b))
+  end
+
+  def self.calculate_final_bearing(from, to)
+    bear = bearing(to, from)
+    (bear + 180) % 360
+  end
+end

data/lib/turf/boolean_point_in_polygon.rb

@@ -0,0 +1,78 @@
+# frozen_string_literal: true
+
+# https://github.com/Turfjs/turf/blob/master/packages/turf-boolean-point-in-polygon/index.ts
+def in_bbox(point, bbox)
+  bbox[0] <= point[0] &&
+    bbox[1] <= point[1] &&
+    bbox[2] >= point[0] &&
+    bbox[3] >= point[1]
+end
+
+def in_ring(point, ring, ignore_boundary)
+  is_inside = false
+  is_ring_valid = ring[0][0] == ring[ring.size - 1][0]
+  is_ring_valid &&= ring[0][1] == ring[ring.size - 1][1]
+  if is_ring_valid
+    ring = ring.slice(0, ring.size - 1)
+  end
+  ring.each_with_index do |ring_pt, ring_pt_index|
+    ring_pt2 = ring[(ring_pt_index + 1) % ring.size]
+
+    xi = ring_pt[0]
+    yi = ring_pt[1]
+    xj = ring_pt2[0]
+    yj = ring_pt2[1]
+
+    on_boundary = (
+      point[1] * (xi - xj) + yi * (xj - point[0]) + yj * (point[0] - xi)
+    ).zero?
+    on_boundary &&= ((xi - point[0]) * (xj - point[0]) <= 0)
+    on_boundary &&= ((yi - point[1]) * (yj - point[1]) <= 0)
+    if on_boundary
+      return !ignore_boundary
+    end
+
+    intersect = ((yi > point[1]) != (yj > point[1])) &&
+                (point[0] < (xj - xi) * (point[1] - yi) / (yj - yi) + xi)
+    if intersect
+      is_inside = !is_inside
+    end
+  end
+  is_inside
+end
+
+module Turf
+  def self.boolean_point_in_polygon(point, polygon, ignore_boundary: false)
+    polygon = deep_symbolize_keys(polygon)
+    pt = get_coord(point)
+    geom = get_geom(polygon)
+    type = geom.fetch(:type)
+    bbox = polygon[:bbox]
+    polys = geom.fetch(:coordinates)
+
+    # Quick elimination if point is not inside bbox
+    return false if bbox && !in_bbox(pt, bbox)
+
+    # normalize to multipolygon
+    polys = [polys] if type == "Polygon"
+
+    inside_poly = false
+    polys.each do |poly|
+      # check if it is in the outer ring first
+      next unless in_ring(pt, poly[0], ignore_boundary)
+
+      in_hole = false
+
+      # check for the point in any of the holes
+      poly.slice(1, poly.size - 1).each do |hole|
+        if in_ring(pt, hole, !ignore_boundary)
+          in_hole = true
+        end
+      end
+      if !in_hole
+        inside_poly = true
+      end
+    end
+    inside_poly
+  end
+end

data/lib/turf/destination.rb

@@ -0,0 +1,22 @@
+# frozen_string_literal: true
+
+module Turf
+  def self.destination(origin, distance, bearing, **options)
+    coordinates1 = get_coord origin
+    longitude1 = degrees_to_radians coordinates1[0]
+    latitude1 = degrees_to_radians coordinates1[1]
+    bearing_radians = degrees_to_radians bearing
+    radians = length_to_radians distance, options[:units]
+
+    latitude2 = Math.asin(Math.sin(latitude1) * Math.cos(radians) +
+      Math.cos(latitude1) * Math.sin(radians) * Math.cos(bearing_radians))
+    longitude2 = longitude1 + Math.atan2(
+      Math.sin(bearing_radians) * Math.sin(radians) * Math.cos(latitude1),
+      Math.cos(radians) - Math.sin(latitude1) * Math.sin(latitude2),
+    )
+    lng = radians_to_degrees(longitude2)
+    lat = radians_to_degrees(latitude2)
+
+    point([lng, lat], options[:properties])
+  end
+end

data/lib/turf/helper.rb

@@ -33,6 +33,15 @@ module Turf
     feat
   end
 
+  def self.feature_collection(features, options = {})
+    fc = { type: "FeatureCollection" }
+    fc[:id] = options[:id] if options[:id]
+    fc[:bbox] = options[:bbox] if options[:bbox]
+    fc[:features] = features
+
+    fc
+  end
+
   def self.line_string(coordinates, properties = nil, options = {})
     if coordinates.size < 2
       raise Error, "coordinates must be an array of two or more positions"
@@ -75,6 +84,11 @@ module Turf
     feature(geom, properties, options)
   end
 
+  def self.radians_to_degrees(radians)
+    degrees = radians.remainder(2 * Math::PI)
+    degrees * 180 / Math::PI
+  end
+
   def self.degrees_to_radians(degrees)
     radians = degrees.remainder(360)
     radians * Math::PI / 180

data/lib/turf/invariant.rb

@@ -26,4 +26,10 @@ module Turf
 
     raise Error, "coord must be GeoJSON Point or an Array of numbers"
   end
+
+  def self.get_geom(geojson)
+    return geojson[:geometry] if geojson[:type] == "Feature"
+
+    geojson
+  end
 end

data/lib/turf/version.rb

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

data/lib/turf_ruby.rb

@@ -1,6 +1,10 @@
 # frozen_string_literal: true
 
 require "turf"
+require "turf/along"
+require "turf/bearing"
+require "turf/destination"
 require "turf/helper"
 require "turf/measurement"
 require "turf/invariant"
+require "turf/boolean_point_in_polygon"

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: turf-ruby
 version: !ruby/object:Gem::Version
-  version: 0.2.0
+  version: 0.3.0
 platform: ruby
 authors:
 - Rafael Santos
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2020-03-10 00:00:00.000000000 Z
+date: 2020-10-05 00:00:00.000000000 Z
 dependencies: []
 description: Turf Ruby is a Ruby library for spatial analysis. It includes traditional
   spatial operations, helper functions for creating GeoJSON data, and data classification
@@ -32,6 +32,10 @@ files:
 - bin/console
 - bin/setup
 - lib/turf.rb
+- lib/turf/along.rb
+- lib/turf/bearing.rb
+- lib/turf/boolean_point_in_polygon.rb
+- lib/turf/destination.rb
 - lib/turf/helper.rb
 - lib/turf/invariant.rb
 - lib/turf/measurement.rb