Class: Geospatial::Polygon

Inherits:
Object
  • Object
show all
Defined in:
lib/geospatial/polygon.rb

Instance Attribute Summary collapse

Class Method Summary collapse

Instance Method Summary collapse

Constructor Details

#initialize(points, bounding_box = nil) ⇒ Polygon

Returns a new instance of Polygon.



41
42
43
44
 
# File 'lib/geospatial/polygon.rb', line 41

def initialize(points, bounding_box = nil)
  @points = points
  @bounding_box = bounding_box
end

Instance Attribute Details

#points ⇒ Object (readonly)

Returns the value of attribute points.



46
47
48
 
# File 'lib/geospatial/polygon.rb', line 46

def points
  @points
end

Class Method Details

.[](*points) ⇒ Object 



25
26
27
 
# File 'lib/geospatial/polygon.rb', line 25

def self.[] *points
  self.new(points)
end

.dump(polygon) ⇒ Object 



35
36
37
38
39
 
# File 'lib/geospatial/polygon.rb', line 35

def self.dump(polygon)
  if polygon
    JSON.dump(polygon.points.map(&:to_a))
  end
end

.is_left(p0, p1, p2) ⇒ Object 



92
93
94
95
96
97
 
# File 'lib/geospatial/polygon.rb', line 92

def self.is_left(p0, p1, p2)
  a = p1 - p0
  b = p2 - p0
  
  return (a[0] * b[1]) - (b[0] * a[1])
end

.load(data) ⇒ Object 



29
30
31
32
33
 
# File 'lib/geospatial/polygon.rb', line 29

def self.load(data)
  if data
    self.new(JSON.parse(data).map{|point| Vector.elements(point)})
  end
end

Instance Method Details

#[](index) ⇒ Object 



48
49
50
51
52
53
 
# File 'lib/geospatial/polygon.rb', line 48

def [] index
  a = @points[index.floor]
  b = @points[index.ceil % @points.size]
  
  return a + (b - a) * (index % 1.0)
end

#bounding_box ⇒ Object 



59
60
61
 
# File 'lib/geospatial/polygon.rb', line 59

def bounding_box
  @bounding_box ||= Box.enclosing_points(@points).freeze
end

#edge_intersection(a, b) ⇒ Object 



135
136
137
138
139
140
141
142
143
144
145
146
147
 
# File 'lib/geospatial/polygon.rb', line 135

def edge_intersection(a, b)
  line = Line.new(a, b)
  
  edges.each_with_index do |(pa, pb), i|
    edge = Line.new(pa, pb)
    
    if line.intersect?(edge)
      return i
    end
  end
  
  return nil
end

#edges ⇒ Object 



70
71
72
73
74
75
76
77
78
 
# File 'lib/geospatial/polygon.rb', line 70

def edges
  return to_enum(:edges) unless block_given?
  
  size = @points.size
  
  @points.each_with_index do |point, index|
    yield point, @points[(index+1)%size]
  end
end

#freeze ⇒ Object 



63
64
65
66
67
68
 
# File 'lib/geospatial/polygon.rb', line 63

def freeze
  @points.freeze
  bounding_box.freeze
  
  super
end

#include?(other) ⇒ Boolean

Returns:

  • (Boolean) 


158
159
160
 
# File 'lib/geospatial/polygon.rb', line 158

def include?(other)
  other.corners.all?{|corner| self.include_point?(corner)}
end

#include_point?(point) ⇒ Boolean

Returns:

  • (Boolean) 


121
122
123
124
125
 
# File 'lib/geospatial/polygon.rb', line 121

def include_point?(point)
  return false unless bounding_box.include_point?(point)
  
  self.winding_number(point) == 1
end

#intersect?(other) ⇒ Boolean

Returns:

  • (Boolean) 


149
150
151
152
153
154
155
156
 
# File 'lib/geospatial/polygon.rb', line 149

def intersect?(other)
  case other
  when Box
    intersect_with_box?(other)
  when Circle
    intersect_with_circle?(other)
  end
end

#intersect_with_box?(other) ⇒ Boolean

Returns:

  • (Boolean) 


127
128
129
130
131
132
133
 
# File 'lib/geospatial/polygon.rb', line 127

def intersect_with_box?(other)
  return true if @points.any?{|point| other.include_point?(point)}
  
  return true if other.corners.any?{|corner| self.include_point?(corner)}
  
  return false
end

#simplify(minimum_distance = 1) ⇒ Object 



80
81
82
83
84
85
86
87
88
89
90
 
# File 'lib/geospatial/polygon.rb', line 80

def simplify(minimum_distance = 1)
  simplified_points = @points.first(1)
  
  @points.each do |point|
    if yield(simplified_points.last, point)
      simplified_points << point
    end
  end
  
  self.class.new(simplified_points, bounding_box)
end

#to_s ⇒ Object 



55
56
57
 
# File 'lib/geospatial/polygon.rb', line 55

def to_s
  "#{self.class}#{@points.inspect}"
end

#winding_number(p) ⇒ Number

Test a 2D point for inclusion in the polygon.

Parameters:

  • p (Vector) —

    The point to test.

Returns:

  • (Number) —

    The number of times the polygon winds around the point (0 if outside).



102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
 
# File 'lib/geospatial/polygon.rb', line 102

def winding_number(p)
  count = 0
  
  edges.each do |pa, pb|
    if pa[1] <= p[1] 
      if pb[1] >= p[1] and Polygon.is_left(pa, pb, p) > 0
        count += 1
      end
    else
      if pb[1] <= p[1] and Polygon.is_left(pa, pb, p) < 0
        count -= 1
      end
    end
    
  end
  
  return count
end