ruby-geometry 0.0.2 → 0.0.3

data/lib/algorithms/point_in_polygon.rb

@@ -0,0 +1,64 @@
+module Geometry
+  class PointInPolygon < Struct.new(:point, :polygon)
+    extend ActiveSupport::Memoizable
+    
+    def inside?
+      point_location == :inside
+    end
+
+    def outside?
+      point_location == :outside
+    end
+
+    def on_the_boundary?
+      point_location == :on_the_boundary
+    end
+
+    def point_location
+      return :outside unless bounding_box.contains?(point)
+      return :on_the_boundary if point_is_vertex? || point_on_edge?
+
+      intersection_count(choose_good_ray).odd? ? :inside : :outside
+    end
+
+    delegate :vertices, :edges, :bounding_box, :to => :polygon
+    memoize :point_location, :edges, :bounding_box
+
+    private
+
+    def point_is_vertex?
+      vertices.any? { |vertex| vertex == point }
+    end
+
+    def point_on_edge?
+      edges.any? { |edge| edge.contains_point?(point) }
+    end
+
+    def choose_good_ray
+      ray = random_ray
+      while ! good_ray?(ray) do
+        ray = random_ray
+      end
+      ray
+    end
+
+    def good_ray?(ray)
+      edges.none? { |edge| edge.parallel_to?(ray) } && vertices.none? { |vertex| ray.contains_point?(vertex) }
+    end
+    
+    def intersection_count(ray)
+      edges.select { |edge| edge.intersects_with?(ray) }.size
+    end
+    
+    def random_ray
+      random_direction = rand * (2 * Math::PI)
+      
+      ray_endpoint = Point sufficient_ray_radius * Math.cos(random_direction), sufficient_ray_radius * Math.sin(random_direction)
+      Segment point, ray_endpoint
+    end
+
+    def sufficient_ray_radius
+      @sufficient_ray_radius ||= bounding_box.diagonal.length * 2
+    end
+  end
+end

data/lib/bounding_box.rb

@@ -0,0 +1,11 @@
+module Geometry
+  class BoundingBox < Struct.new(:leftbottom, :righttop)
+    def diagonal
+      Segment leftbottom, righttop
+    end
+
+    def contains?(point)
+      point.x.between?(leftbottom.x, righttop.x) && point.y.between?(leftbottom.y, righttop.y)
+    end
+  end
+end

data/lib/geometry.rb

@@ -1,8 +1,13 @@
+require "active_support/core_ext/module/delegation"
+require "active_support/memoizable"
+
 require 'point'
 require 'segment'
 require 'vector'
 require 'polygon'
 require 'line'
+require 'bounding_box'
+require 'algorithms/point_in_polygon'
 
 module Geometry
   include Math

data/lib/point.rb

@@ -8,4 +8,8 @@ module Geometry
       x === another_point.x && y === another_point.y
     end
   end
-end
+end
+
+def Point(x, y)
+  Geometry::Point.new(x, y)
+end

data/lib/polygon.rb

@@ -9,5 +9,21 @@ module Geometry
 
       edges << Segment.new(vertices.last, vertices.first)
     end
+
+    def bounding_box
+      leftbottom = Point vertices.map(&:x).min, vertices.map(&:y).min
+      righttop = Point vertices.map(&:x).max, vertices.map(&:y).max
+      
+      BoundingBox.new leftbottom, righttop
+    end
+
+    def contains?(point)
+      point_in_polygon = PointInPolygon.new(point, self)
+      point_in_polygon.inside? || point_in_polygon.on_the_boundary?
+    end
   end
-end
+end
+
+def Polygon(vertices)
+  Geometry::Polygon.new(vertices)
+end

data/lib/segment.rb

@@ -102,4 +102,8 @@ module Geometry
         self.to_vector.cross_product(vector_to_second_endpoint) <= 0
     end
   end
-end
+end
+
+def Segment(point1, point2)
+  Geometry::Segment.new point1, point2
+end

data/test/bounding_box/contains_test.rb

@@ -0,0 +1,29 @@
+require 'test/unit'
+require 'geometry'
+
+class ContainsTest < Test::Unit::TestCase
+  include Geometry
+
+  def bounding_box
+    BoundingBox.new Point(-1, -1), Point(1, 1)
+  end
+
+  def test_contains
+    assert bounding_box.contains?(Point(0, 0))
+  end
+
+  def test_does_not_contain
+    assert ! bounding_box.contains?(Point(2, 2))
+    assert ! bounding_box.contains?(Point(-2, -2))
+    assert ! bounding_box.contains?(Point(1, 2))
+    assert ! bounding_box.contains?(Point(0, -1.1))
+  end
+
+  def test_on_the_edge
+    assert bounding_box.contains?(Point(1, 0))
+  end
+
+  def test_at_vertex
+    assert bounding_box.contains?(Point(1, 1))
+  end
+end

data/test/polygon/bounding_box_test.rb

@@ -0,0 +1,43 @@
+require 'test/unit'
+require 'geometry'
+
+class BoundingBoxTest < Test::Unit::TestCase
+  include Geometry
+
+  def test_rectangle
+    polygon = Polygon [
+                       Point(-1, -1),
+                       Point(1, -1),
+                       Point(1, 1),
+                       Point(-1, 1)
+                      ]
+    assert polygon.bounding_box.leftbottom == Point(-1, -1)
+    assert polygon.bounding_box.righttop == Point(1, 1)
+  end
+
+  def test_sloped_rectangle
+    polygon = Polygon [
+                       Point(-1, 0),
+                       Point(0, -1),
+                       Point(1, 0),
+                       Point(0, 1)
+                      ]
+    
+    assert polygon.bounding_box.leftbottom == Point(-1, -1)
+    assert polygon.bounding_box.righttop == Point(1, 1)
+  end
+
+  def test_nonconvex
+    polygon = Polygon [
+                       Point(0, 2),
+                       Point(2, 0),
+                       Point(5, 1),
+                       Point(3, 8),
+                       Point(3, 1),
+                       Point(1, 4)
+                      ]
+
+    assert polygon.bounding_box.leftbottom == Point(0, 0)
+    assert polygon.bounding_box.righttop == Point(5, 8)
+  end
+end

data/test/polygon/contains_point_test.rb

@@ -0,0 +1,73 @@
+require 'test/unit'
+require 'geometry'
+
+class ContainsPointTest < Test::Unit::TestCase
+  include Geometry
+
+  def test_convex
+    rectangle = Polygon.new [
+                             Point(0, 0),
+                             Point(1, 0),
+                             Point(1, 1),
+                             Point(0, 1)
+                            ]
+    
+
+    inner = Point(0.5, 0.5)
+    assert rectangle.contains?(inner)
+
+    outer = Point(1.5, 1.5)
+    assert ! rectangle.contains?(outer)
+
+    on_edge = Point(0.5, 1)
+    assert rectangle.contains?(on_edge)
+
+    at_vertex = Point(1, 1)
+    assert rectangle.contains?(at_vertex)
+  end
+  
+
+  # +--------+
+  # |        |
+  # |    +---+
+  # |    |
+  # |    +---+
+  # |        |
+  # +--------+ 
+
+  def test_nonconvex
+    nonconvex_polygon = Polygon.new [
+                                     Point(0, 0),
+                                     Point(0, 6),
+                                     Point(4, 6),
+                                     Point(4, 4),
+                                     Point(2, 4),
+                                     Point(2, 2),
+                                     Point(4, 2),
+                                     Point(4, 0)
+                                    ]
+
+    inner_points = [
+                    Point(1, 5),
+                    Point(3, 5),
+                    Point(1, 3),
+                    Point(1, 1),
+                    Point(3, 1)
+                   ]
+
+    outer_points = [
+                    Point(7, 5),
+                    Point(5, 3),
+                    Point(7, 3),
+                    Point(7, 1)
+                   ]
+
+    inner_points.each do |inner_point|
+      assert nonconvex_polygon.contains?(inner_point)
+    end
+
+    outer_points.each do |outer_point|
+      assert ! nonconvex_polygon.contains?(outer_point)
+    end
+  end
+end

metadata

@@ -1,125 +1,115 @@
---- !ruby/object:Gem::Specification 
+--- !ruby/object:Gem::Specification
 name: ruby-geometry
-version: !ruby/object:Gem::Version 
-  prerelease: false
-  segments: 
-  - 0
-  - 0
-  - 2
-  version: 0.0.2
+version: !ruby/object:Gem::Version
+  version: 0.0.3
+  prerelease: 
 platform: ruby
-authors: 
+authors:
 - Daniel Vartanov
 autorequire: 
 bindir: bin
 cert_chain: []
-
-date: 2010-10-17 00:00:00 -04:00
-default_executable: 
+date: 2011-12-04 00:00:00.000000000Z
 dependencies: []
-
 description: 
 email: daniel.vartanov@gmail.com
 executables: []
-
 extensions: []
-
 extra_rdoc_files: []
-
-files: 
-- lib/point.rb
-- lib/line.rb
+files:
+- lib/algorithms/point_in_polygon.rb
+- lib/bounding_box.rb
 - lib/comparison_with_precision.rb
-- lib/polygon.rb
-- lib/vector.rb
 - lib/geometry.rb
+- lib/line.rb
+- lib/point.rb
+- lib/polygon.rb
 - lib/segment.rb
+- lib/vector.rb
+- test/bounding_box/contains_test.rb
 - test/geometry/distance_test.rb
+- test/line/angle_to_test.rb
+- test/line/horizontal_test.rb
+- test/line/initialize_test.rb
+- test/line/intersect_x_test.rb
+- test/line/parallel_to_test.rb
+- test/line/slope_test.rb
+- test/line/vertical_test.rb
+- test/line/x_intercept_test.rb
+- test/line/y_intercept_test.rb
+- test/point/equals_test.rb
+- test/point/initialize_test.rb
+- test/polygon/bounding_box_test.rb
+- test/polygon/contains_point_test.rb
 - test/polygon/edges_test.rb
+- test/segment/bounds_test.rb
 - test/segment/contains_point_test.rb
-- test/segment/overlaps_test.rb
-- test/segment/to_vector_test.rb
-- test/segment/length_test.rb
+- test/segment/initialize_test.rb
 - test/segment/intersection_point_with_test.rb
-- test/segment/parallel_to_test.rb
 - test/segment/intersects_with_test.rb
-- test/segment/bounds_test.rb
-- test/segment/initialize_test.rb
-- test/line/x_intercept_test.rb
-- test/line/y_intercept_test.rb
-- test/line/parallel_to_test.rb
-- test/line/angle_to_test.rb
-- test/line/initialize_test.rb
-- test/line/vertical_test.rb
-- test/line/slope_test.rb
-- test/line/intersect_x_test.rb
-- test/line/horizontal_test.rb
-- test/vector/equals_test.rb
-- test/vector/modulus_test.rb
+- test/segment/length_test.rb
+- test/segment/overlaps_test.rb
+- test/segment/parallel_to_test.rb
+- test/segment/to_vector_test.rb
 - test/vector/arithmetics_test.rb
+- test/vector/collinear_with_test.rb
 - test/vector/cross_product_test.rb
+- test/vector/equals_test.rb
+- test/vector/modulus_test.rb
 - test/vector/scalar_product_test.rb
-- test/vector/collinear_with_test.rb
-- test/point/equals_test.rb
-- test/point/initialize_test.rb
-has_rdoc: true
 homepage: http://github.com/DanielVartanov/ruby-geometry
 licenses: []
-
 post_install_message: 
 rdoc_options: []
-
-require_paths: 
+require_paths:
 - lib
-required_ruby_version: !ruby/object:Gem::Requirement 
+required_ruby_version: !ruby/object:Gem::Requirement
   none: false
-  requirements: 
-  - - ">="
-    - !ruby/object:Gem::Version 
-      segments: 
-      - 0
-      version: "0"
-required_rubygems_version: !ruby/object:Gem::Requirement 
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
   none: false
-  requirements: 
-  - - ">="
-    - !ruby/object:Gem::Version 
-      segments: 
-      - 0
-      version: "0"
+  requirements:
+  - - ! '>='
+    - !ruby/object:Gem::Version
+      version: '0'
 requirements: []
-
 rubyforge_project: 
-rubygems_version: 1.3.7
+rubygems_version: 1.8.10
 signing_key: 
 specification_version: 3
 summary: Implementation of basic 2D geometry algorithms in Ruby
-test_files: 
+test_files:
+- test/bounding_box/contains_test.rb
 - test/geometry/distance_test.rb
+- test/line/angle_to_test.rb
+- test/line/horizontal_test.rb
+- test/line/initialize_test.rb
+- test/line/intersect_x_test.rb
+- test/line/parallel_to_test.rb
+- test/line/slope_test.rb
+- test/line/vertical_test.rb
+- test/line/x_intercept_test.rb
+- test/line/y_intercept_test.rb
+- test/point/equals_test.rb
+- test/point/initialize_test.rb
+- test/polygon/bounding_box_test.rb
+- test/polygon/contains_point_test.rb
 - test/polygon/edges_test.rb
+- test/segment/bounds_test.rb
 - test/segment/contains_point_test.rb
-- test/segment/overlaps_test.rb
-- test/segment/to_vector_test.rb
-- test/segment/length_test.rb
+- test/segment/initialize_test.rb
 - test/segment/intersection_point_with_test.rb
-- test/segment/parallel_to_test.rb
 - test/segment/intersects_with_test.rb
-- test/segment/bounds_test.rb
-- test/segment/initialize_test.rb
-- test/line/x_intercept_test.rb
-- test/line/y_intercept_test.rb
-- test/line/parallel_to_test.rb
-- test/line/angle_to_test.rb
-- test/line/initialize_test.rb
-- test/line/vertical_test.rb
-- test/line/slope_test.rb
-- test/line/intersect_x_test.rb
-- test/line/horizontal_test.rb
-- test/vector/equals_test.rb
-- test/vector/modulus_test.rb
+- test/segment/length_test.rb
+- test/segment/overlaps_test.rb
+- test/segment/parallel_to_test.rb
+- test/segment/to_vector_test.rb
 - test/vector/arithmetics_test.rb
+- test/vector/collinear_with_test.rb
 - test/vector/cross_product_test.rb
+- test/vector/equals_test.rb
+- test/vector/modulus_test.rb
 - test/vector/scalar_product_test.rb
-- test/vector/collinear_with_test.rb
-- test/point/equals_test.rb
-- test/point/initialize_test.rb