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 



112
113
114
115
116
117
 
# File 'lib/geospatial/polygon.rb', line 112

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

#area(radius = 1.0) ⇒ Object

Parameters:

  • radius (Float) (defaults to: 1.0) —

    The radius of the sphere on which to compute the area.



81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
 
# File 'lib/geospatial/polygon.rb', line 81

def area(radius = 1.0)
	if @points.size > 2
		area = 0.0
		
		self.edges.each do |p1, p2|
			r1 = (p2[0] - p1[0]) * D2R
			r2 = 2 + Math::sin(p1[1] * D2R) + Math::sin(p2[1] * D2R)
			
			area += r1 * r2
		end
		
		return (area * radius * radius / 2.0).abs
	else
		return 0.0
	end
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 



155
156
157
158
159
160
161
162
163
164
165
166
167
 
# File 'lib/geospatial/polygon.rb', line 155

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) 


178
179
180
 
# File 'lib/geospatial/polygon.rb', line 178

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

#include_point?(point) ⇒ Boolean

Returns:

  • (Boolean) 


141
142
143
144
145
 
# File 'lib/geospatial/polygon.rb', line 141

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

#intersect?(other) ⇒ Boolean

Returns:

  • (Boolean) 


169
170
171
172
173
174
175
176
 
# File 'lib/geospatial/polygon.rb', line 169

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) 


147
148
149
150
151
152
153
 
# File 'lib/geospatial/polygon.rb', line 147

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 



98
99
100
101
102
103
104
105
106
107
108
109
110
 
# File 'lib/geospatial/polygon.rb', line 98

def simplify(minimum_distance = 1)
	simplified_points = @points.first(1)
	
	@points.each_with_index do |point, index|
		next_point = @points[(index+1) % @points.size]
		
		if yield(simplified_points.last, point, next_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).



122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
 
# File 'lib/geospatial/polygon.rb', line 122

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