spatial_features 0.0.1

Sign up to get free protection for your applications and to get access to all the features.
Files changed (51) hide show
  1. checksums.yaml +7 -0
  2. data/MIT-LICENSE +20 -0
  3. data/README.rdoc +3 -0
  4. data/Rakefile +19 -0
  5. data/lib/spatial_features/caching.rb +66 -0
  6. data/lib/spatial_features/controller_helpers/spatial_extensions.rb +35 -0
  7. data/lib/spatial_features/has_spatial_features.rb +158 -0
  8. data/lib/spatial_features/import/arcgis_kmz_features.rb +104 -0
  9. data/lib/spatial_features/models/feature.rb +59 -0
  10. data/lib/spatial_features/models/spatial_cache.rb +3 -0
  11. data/lib/spatial_features/models/spatial_proximity.rb +4 -0
  12. data/lib/spatial_features/venn_polygons.rb +56 -0
  13. data/lib/spatial_features/version.rb +3 -0
  14. data/lib/spatial_features.rb +21 -0
  15. data/lib/tasks/spatial_features_tasks.rake +4 -0
  16. data/test/dummy/README.rdoc +28 -0
  17. data/test/dummy/Rakefile +6 -0
  18. data/test/dummy/app/assets/javascripts/application.js +13 -0
  19. data/test/dummy/app/assets/stylesheets/application.css +15 -0
  20. data/test/dummy/app/controllers/application_controller.rb +5 -0
  21. data/test/dummy/app/helpers/application_helper.rb +2 -0
  22. data/test/dummy/app/views/layouts/application.html.erb +14 -0
  23. data/test/dummy/bin/bundle +3 -0
  24. data/test/dummy/bin/rails +4 -0
  25. data/test/dummy/bin/rake +4 -0
  26. data/test/dummy/config/application.rb +23 -0
  27. data/test/dummy/config/boot.rb +5 -0
  28. data/test/dummy/config/database.yml +25 -0
  29. data/test/dummy/config/environment.rb +5 -0
  30. data/test/dummy/config/environments/development.rb +37 -0
  31. data/test/dummy/config/environments/production.rb +82 -0
  32. data/test/dummy/config/environments/test.rb +39 -0
  33. data/test/dummy/config/initializers/assets.rb +8 -0
  34. data/test/dummy/config/initializers/backtrace_silencers.rb +7 -0
  35. data/test/dummy/config/initializers/cookies_serializer.rb +3 -0
  36. data/test/dummy/config/initializers/filter_parameter_logging.rb +4 -0
  37. data/test/dummy/config/initializers/inflections.rb +16 -0
  38. data/test/dummy/config/initializers/mime_types.rb +4 -0
  39. data/test/dummy/config/initializers/session_store.rb +3 -0
  40. data/test/dummy/config/initializers/wrap_parameters.rb +14 -0
  41. data/test/dummy/config/locales/en.yml +23 -0
  42. data/test/dummy/config/routes.rb +56 -0
  43. data/test/dummy/config/secrets.yml +22 -0
  44. data/test/dummy/config.ru +4 -0
  45. data/test/dummy/public/404.html +67 -0
  46. data/test/dummy/public/422.html +67 -0
  47. data/test/dummy/public/500.html +66 -0
  48. data/test/dummy/public/favicon.ico +0 -0
  49. data/test/spatial_features_test.rb +7 -0
  50. data/test/test_helper.rb +15 -0
  51. metadata +171 -0
checksums.yaml ADDED
@@ -0,0 +1,7 @@
1
+ ---
2
+ SHA1:
3
+ metadata.gz: c42265b78f3d41490bf8e89691d57f1cab48d9b2
4
+ data.tar.gz: bdc3542a36c9584b6f18bceddda2b7e175997a2e
5
+ SHA512:
6
+ metadata.gz: a17091a2662a74ac313c47059ac13f63a220a45fbce3c11aedffc7ccc6c74c8425b07aec89ccc0a1b00772c5405bcaa5432ac39bd20aa8398d93ac405d7407eb
7
+ data.tar.gz: 7a15bb57318905523faa487ff24ac34b48a331e893d4088dfa77aff0c8dda933604599428f2030263d9b67bf287cfd4e316486b007c893978c51ab9cd263ea0d
data/MIT-LICENSE ADDED
@@ -0,0 +1,20 @@
1
+ Copyright 2014 YOURNAME
2
+
3
+ Permission is hereby granted, free of charge, to any person obtaining
4
+ a copy of this software and associated documentation files (the
5
+ "Software"), to deal in the Software without restriction, including
6
+ without limitation the rights to use, copy, modify, merge, publish,
7
+ distribute, sublicense, and/or sell copies of the Software, and to
8
+ permit persons to whom the Software is furnished to do so, subject to
9
+ the following conditions:
10
+
11
+ The above copyright notice and this permission notice shall be
12
+ included in all copies or substantial portions of the Software.
13
+
14
+ THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
15
+ EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
16
+ MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
17
+ NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
18
+ LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
19
+ OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
20
+ WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
data/README.rdoc ADDED
@@ -0,0 +1,3 @@
1
+ = Spatial Features
2
+
3
+ Adds spatial methods to a model.
data/Rakefile ADDED
@@ -0,0 +1,19 @@
1
+ begin
2
+ require 'bundler/setup'
3
+ rescue LoadError
4
+ puts 'You must `gem install bundler` and `bundle install` to run rake tasks'
5
+ end
6
+
7
+ Bundler::GemHelper.install_tasks
8
+
9
+ require 'rake/testtask'
10
+
11
+ Rake::TestTask.new(:test) do |t|
12
+ t.libs << 'lib'
13
+ t.libs << 'test'
14
+ t.pattern = 'test/**/*_test.rb'
15
+ t.verbose = false
16
+ end
17
+
18
+
19
+ task default: :test
@@ -0,0 +1,66 @@
1
+ module SpatialFeatures
2
+ mattr_accessor :default_cache_buffer_in_meters
3
+ self.default_cache_buffer_in_meters = 100
4
+
5
+ # Create or update the spatial cache of a spatial class in relation to another
6
+ # NOTE: Arguments are order independent, so their names do not reflect the _a _b
7
+ # naming scheme used in other cache methods
8
+ def self.cache_proximity(klass, clazz)
9
+ clear_cache(klass, clazz)
10
+
11
+ klass.find_each do |record|
12
+ create_spatial_proximities(record, clazz)
13
+ create_spatial_cache(record, clazz)
14
+ end
15
+
16
+ clazz.find_each do |record|
17
+ create_spatial_cache(record, klass)
18
+ end
19
+ end
20
+
21
+ # Create or update the spatial cache of a single record in relation to another spatial class
22
+ def self.cache_record_proximity(record, klass)
23
+ clear_record_cache(record, klass)
24
+ create_spatial_proximities(record, klass)
25
+ create_spatial_cache(record, klass)
26
+ end
27
+
28
+ # Delete all cache entries relating klass to clazz
29
+ def self.clear_cache(klass = nil, clazz = nil)
30
+ if klass.blank? && clazz.blank?
31
+ SpatialCache.delete_all
32
+ SpatialProximity.delete_all
33
+ else
34
+ SpatialCache.where(:spatial_model_type => klass, :intersection_model_type => clazz).delete_all
35
+ SpatialCache.where(:spatial_model_type => clazz, :intersection_model_type => klass).delete_all
36
+ SpatialProximity.where(:model_a_type => klass, :model_b_type => clazz).delete_all
37
+ SpatialProximity.where(:model_a_type => clazz, :model_b_type => klass).delete_all
38
+ end
39
+ end
40
+
41
+ def self.clear_record_cache(record, klass)
42
+ record.spatial_cache.where(:intersection_model_type => klass.name).delete_all
43
+ SpatialProximity.where(:model_a_type => record.class.name, :model_a_id => record.id, :model_b_type => klass.name).delete_all
44
+ SpatialProximity.where(:model_b_type => record.class.name, :model_b_id => record.id, :model_a_type => klass.name).delete_all
45
+ end
46
+
47
+ def self.create_spatial_proximities(record, klass)
48
+ record_is_a = record.class.name < klass.name
49
+
50
+ scope = klass.within_buffer(record, default_cache_buffer_in_meters, :intersection_area => true, :distance => true, :cache => false)
51
+ scope.find_each do |klass_record|
52
+ SpatialProximity.create!(
53
+ :model_a => record_is_a ? record : klass_record,
54
+ :model_b => record_is_a ? klass_record : record,
55
+ :distance_in_meters => klass_record.distance_in_meters,
56
+ :intersection_area_in_square_meters => klass_record.intersection_area_in_square_meters)
57
+ end
58
+ end
59
+
60
+ def self.create_spatial_cache(model, klass)
61
+ SpatialCache.create!(
62
+ :spatial_model => model,
63
+ :intersection_model_type => klass.name,
64
+ :cache_distance => default_cache_buffer_in_meters)
65
+ end
66
+ end
@@ -0,0 +1,35 @@
1
+ module SpatialExtensions
2
+ private
3
+
4
+ def abstract_proximity_action(scope, target, distance, &block)
5
+ @nearby_records = scope_for_search(scope).within_buffer(target, distance, :distance => true, :intersection_area => true).order('distance_in_meters ASC, intersection_area_in_square_meters DESC, id ASC')
6
+ @target = target
7
+
8
+ yield if block_given?
9
+
10
+ respond_to do |format|
11
+ format.html { render :template => 'shared/spatial/feature_proximity', :layout => false }
12
+ format.kml { render :template => 'shared/spatial/feature_proximity' }
13
+ end
14
+ end
15
+
16
+ def abstract_venn_polygons_action(scope, target, &block)
17
+ @venn_polygons = SpatialFeatures.venn_polygons(scope_for_search(scope).intersecting(target), target.class.where(:id => target), :simplified => false)
18
+ @klass = klass_for_search(scope)
19
+ @target = target
20
+
21
+ yield if block_given?
22
+
23
+ respond_to do |format|
24
+ format.kml { render :template => 'shared/spatial/feature_venn_polygons' }
25
+ end
26
+ end
27
+
28
+ def klass_for_search(scope_or_class)
29
+ scope_or_class.is_a?(ActiveRecord::Relation) ? scope_or_class.klass : scope_or_class
30
+ end
31
+
32
+ def scope_for_search(scope)
33
+ params.key?(:ids) ? scope.where(:id => params[:ids]) : scope
34
+ end
35
+ end
@@ -0,0 +1,158 @@
1
+ module SpatialFeatures
2
+ module ActMethod
3
+ def has_spatial_features(options = {})
4
+ has_many :features, :as => :spatial_model, :dependent => :delete_all
5
+ scope :with_features, lambda { where(:id => Feature.select(:spatial_model_id).where(:spatial_model_type => name)) }
6
+ scope :without_features, lambda { where.not(:id => Feature.select(:spatial_model_id).where(:spatial_model_type => name)) }
7
+
8
+ has_many :spatial_cache, :as => :spatial_model
9
+
10
+ extend SpatialFeatures::ClassMethods
11
+ include SpatialFeatures::InstanceMethods
12
+ end
13
+ end
14
+
15
+ module ClassMethods
16
+ # Add methods to generate cache keys for a record or all records of this class
17
+ # NOTE: features are never updated, only deleted and created, therefore we can
18
+ # tell if they have changed by finding the maximum id and count instead of needing timestamps
19
+ def features_cache_key
20
+ "#{name}/#{Feature.where(:spatial_model_type => self).maximum(:id)}-#{Feature.where(:spatial_model_type => self).count}"
21
+ end
22
+
23
+ def intersecting(other, options = {})
24
+ within_buffer(other, 0, options)
25
+ end
26
+
27
+ def within_buffer(other, buffer_in_meters = 0, options = {})
28
+ raise "Can't intersect with #{other} because it does not implement has_features" unless other.has_spatial_features?
29
+
30
+ if other.spatial_cache_for?(self, buffer_in_meters) && options[:cache] != false # CACHED
31
+ scope = cached_spatial_join(other)
32
+ .select("#{table_name}.*, spatial_proximities.distance_in_meters, spatial_proximities.intersection_area_in_square_meters")
33
+
34
+ scope = scope.where("spatial_proximities.distance_in_meters <= ?", buffer_in_meters) if buffer_in_meters
35
+ else # NON-CACHED
36
+ scope = joins_features_for(other)
37
+ .select("#{table_name}.*")
38
+ .group("#{table_name}.#{primary_key}")
39
+
40
+ scope = scope.where('ST_DWithin(features_for.geog_lowres, features_for_other.geog_lowres, ?)', buffer_in_meters) if buffer_in_meters
41
+ scope = scope.select("MIN(ST_Distance(features_for.geog_lowres, features_for_other.geog_lowres)) AS distance_in_meters") if options[:distance]
42
+ scope = scope.select("SUM(ST_Area(ST_Intersection(features_for.geog_lowres, features_for_other.geog_lowres))) AS intersection_area_in_square_meters") if options[:intersection_area]
43
+ end
44
+
45
+ return scope
46
+ end
47
+
48
+ def polygons
49
+ Feature.polygons.where(:spatial_model_type => self.class)
50
+ end
51
+
52
+ def lines
53
+ Feature.lines.where(:spatial_model_type => self.class)
54
+ end
55
+
56
+ def points
57
+ Feature.points.where(:spatial_model_type => self.class)
58
+ end
59
+
60
+ def cached_spatial_join(other)
61
+ raise "Cannot use cached spatial join for the same class" if other.class.name == self.name
62
+
63
+ other_column = other.class.name < self.name ? :model_a : :model_b
64
+ self_column = other_column == :model_a ? :model_b : :model_a
65
+
66
+ joins("INNER JOIN spatial_proximities ON spatial_proximities.#{self_column}_type = '#{self}' AND spatial_proximities.#{self_column}_id = #{table_name}.id AND spatial_proximities.#{other_column}_type = '#{other.class}' AND spatial_proximities.#{other_column}_id = '#{other.id}'")
67
+ end
68
+
69
+ def joins_features_for(other, table_alias = 'features_for')
70
+ joins_features(table_alias)
71
+ .joins(%Q(INNER JOIN features "#{table_alias}_other" ON "#{table_alias}_other".spatial_model_type = '#{other.class.name}' AND "#{table_alias}_other".spatial_model_id = #{other.id}))
72
+ end
73
+
74
+ def joins_features(table_alias = 'features_for')
75
+ joins(%Q(INNER JOIN features "#{table_alias}" ON "#{table_alias}".spatial_model_type = '#{name}' AND "#{table_alias}".spatial_model_id = #{table_name}.id))
76
+ end
77
+ end
78
+
79
+ module InstanceMethods
80
+ def has_spatial_features?
81
+ true
82
+ end
83
+
84
+ def features_cache_key
85
+ "#{self.class.name}/#{self.id}-#{features.maximum(:id)}-#{features.size}"
86
+ end
87
+
88
+ def polygons?
89
+ !features.polygons.empty?
90
+ end
91
+
92
+ def lines?
93
+ !features.lines.empty?
94
+ end
95
+
96
+ def points?
97
+ !features.points.empty?
98
+ end
99
+
100
+ def features?
101
+ features.present?
102
+ end
103
+
104
+ # Returns true if the model stores a hash of the features so we don't need to process the features if they haven't changed
105
+ def has_spatial_features_hash?
106
+ respond_to?(:features_hash)
107
+ end
108
+
109
+ def intersects?(other)
110
+ self.class.intersecting(other).exists?(self)
111
+ end
112
+
113
+ def total_intersection_area_in_square_meters(klass, options = {})
114
+ self.class
115
+ .select(%Q(ST_Area(ST_Intersection(ST_Union(features_for.geog_lowres::geometry), ST_Union(features_for_other.geog_lowres::geometry))::geography) AS intersection_area_in_square_meters))
116
+ .joins(%Q(INNER JOIN features "features_for" ON "features_for".spatial_model_type = '#{self.class}' AND "features_for".spatial_model_id = #{self.class.table_name}.id))
117
+ .joins(%Q(INNER JOIN features "features_for_other" ON "features_for_other".spatial_model_type = '#{klass}'))
118
+ .where(:id => self.id)
119
+ .where('ST_DWithin(features_for.geog_lowres, features_for_other.geog_lowres, 0)')
120
+ .group("#{self.class.table_name}.id")
121
+ .first
122
+ .try(:intersection_area_in_square_meters) || 0
123
+ end
124
+
125
+ def total_intersection_area_in_hectares(klass)
126
+ Formatters::HECTARES.call(total_intersection_area_in_square_meters(klass))
127
+ end
128
+
129
+ def total_intersection_area_percentage(klass)
130
+ return 0.0 unless features_area_in_square_meters > 0
131
+
132
+ ((total_intersection_area_in_square_meters(klass) / features_area_in_square_meters) * 100).round(1)
133
+ end
134
+
135
+ def features_area_in_square_meters
136
+ @features_area_in_square_meters ||= features.sum('ST_Area(features.geog_lowres)')
137
+ end
138
+
139
+ def features_area_in_hectares
140
+ Formatters::HECTARES.call(features_area_in_square_meters)
141
+ end
142
+
143
+ def spatial_cache_for(klass)
144
+ spatial_cache.where(:intersection_model_type => klass).first
145
+ end
146
+
147
+ def spatial_cache_for?(klass, buffer_in_meters)
148
+ if cache = spatial_cache_for(klass)
149
+ return cache.cache_distance.nil? if buffer_in_meters.nil? # cache must be total if no buffer_in_meters
150
+ return true if cache.cache_distance.nil? # always good if cache is total
151
+
152
+ return buffer_in_meters <= cache.cache_distance
153
+ else
154
+ return false
155
+ end
156
+ end
157
+ end
158
+ end
@@ -0,0 +1,104 @@
1
+ module ArcGISKmzFeatures
2
+ require 'open-uri'
3
+ require 'digest/md5'
4
+
5
+ def update_features!
6
+ @feature_error_messages = []
7
+ kml_array = []
8
+ cache_kml = ''
9
+
10
+ Array(arcgis_kmz_url).each do |url|
11
+ kml_array << open_kmz_url(url)
12
+ cache_kml << kml_array.last.to_s
13
+ end
14
+
15
+ if has_spatial_features_hash?
16
+ new_features_hash = Digest::MD5.hexdigest(cache_kml) if cache_kml.present?
17
+
18
+ if new_features_hash != self.features_hash
19
+ replace_features(kml_array)
20
+ update_attributes(:features_hash => new_features_hash)
21
+ else
22
+ return false
23
+ end
24
+ else
25
+ replace_features(kml_array)
26
+ end
27
+
28
+ return true
29
+ end
30
+
31
+ private
32
+
33
+ def replace_features(kml_array)
34
+ new_features = []
35
+ kml_array.each {|kml| new_features.concat build_features(kml) }
36
+
37
+ ActiveRecord::Base.transaction do
38
+ self.features.destroy_all
39
+ new_features.each(&:save)
40
+
41
+ @feature_error_messages.concat new_features.collect {|feature| "Feature #{feature.name}: #{feature.errors.full_messages.to_sentence}" if feature.errors.present? }.compact.flatten
42
+ if @feature_error_messages.present?
43
+ raise UpdateError, "Error updating #{self.class} #{self.id}. #{@feature_error_messages.to_sentence}"
44
+ end
45
+ end
46
+ end
47
+
48
+ def build_features(kml)
49
+ new_type_features = []
50
+
51
+ extract_kml_features(kml) do |feature_type, feature, name, metadata|
52
+ begin
53
+ new_type_features << build_feature(feature_type, name, metadata, build_geom(feature))
54
+ rescue => e
55
+ @feature_error_messages << e.message
56
+ end
57
+ end
58
+
59
+ return new_type_features
60
+ end
61
+
62
+ # Use ST_Force_2D to discard z-coordinates that cause failures later in the process
63
+ def build_geom(feature)
64
+ if make_valid?
65
+ geom = ActiveRecord::Base.connection.select_value("SELECT ST_CollectionExtract(ST_MakeValid(ST_Force_2D(ST_GeomFromKML('#{feature}'))),3 )")
66
+ else
67
+ geom = ActiveRecord::Base.connection.select_value("SELECT ST_Force_2D(ST_GeomFromKML('#{feature}'))")
68
+ end
69
+ end
70
+
71
+ def extract_kml_features(kml, &block)
72
+ Nokogiri::XML(kml).css('Placemark').each do |placemark|
73
+ name = placemark.css('name').text
74
+ metadata = Hash[Nokogiri::XML(placemark.css('description').text).css('td').collect(&:text).each_slice(2).to_a]
75
+
76
+ {'Polygon' => 'POLYGON', 'LineString' => 'LINE', 'Point' => 'POINT'}.each do |kml_type, sql_type|
77
+ placemark.css(kml_type).each do |feature|
78
+ yield sql_type, feature, name, metadata
79
+ end
80
+ end
81
+ end
82
+ end
83
+
84
+ def build_feature(feature_type, name, metadata, geom)
85
+ Feature.new(spatial_model: self, name: name, metadata: metadata, feature_type: feature_type, geog: geom)
86
+ end
87
+
88
+ def open_kmz_url(url)
89
+ Zip::InputStream.open(open(url)) do |io|
90
+ while (entry = io.get_next_entry)
91
+ return io.read if entry.name.downcase == 'doc.kml'
92
+ end
93
+ end
94
+
95
+ return nil
96
+ end
97
+
98
+ # Can be overridden to use PostGIS to force geometry to be valid
99
+ def make_valid?
100
+ false
101
+ end
102
+
103
+ class UpdateError < StandardError; end
104
+ end
@@ -0,0 +1,59 @@
1
+ class Feature < ActiveRecord::Base
2
+ belongs_to :spatial_model, :polymorphic => :true
3
+
4
+ self.rgeo_factory_generator = RGeo::Geos.factory_generator
5
+
6
+ before_validation :sanitize_feature_type
7
+ validates_presence_of :geog
8
+ validate :geometry_is_valid
9
+ validates_inclusion_of :feature_type, :in => ['polygon', 'point', 'line']
10
+ after_save :cache_derivatives
11
+
12
+ store :metadata
13
+
14
+ def self.polygons
15
+ where(:feature_type => 'polygon')
16
+ end
17
+
18
+ def self.lines
19
+ where(:feature_type => 'line')
20
+ end
21
+
22
+ def self.points
23
+ where(:feature_type => 'point')
24
+ end
25
+
26
+ def self.for_kml(options = {})
27
+ if options[:simplified]
28
+ select("features.name, features.kml, features.metadata").where("features.name IS NULL OR features.name NOT IN ('s', 't')")
29
+ else
30
+ select("features.name, ST_AsKML(features.geog, 6) AS kml, features.metadata")
31
+ end
32
+ end
33
+
34
+ def self.invalid
35
+ select('features.*, ST_IsValidReason(geog::geometry) AS invalid_geometry_message').where.not('ST_IsValid(geog::geometry)')
36
+ end
37
+
38
+ def envelope(buffer_in_meters = 0)
39
+ self.class.select("ST_Envelope(ST_Buffer(features.geog, #{buffer_in_meters})::geometry) AS result").where(:id => id).first.result.exterior_ring.points.values_at(0,2)
40
+ end
41
+
42
+ private
43
+
44
+ def cache_derivatives
45
+ self.class.connection.execute "UPDATE features SET geog_lowres = ST_SimplifyPreserveTopology(geog::geometry, 0.0001) WHERE id = #{self.id}"
46
+ self.class.connection.execute "UPDATE features SET kml = ST_AsKML(geog_lowres::geometry, 5) WHERE id = #{self.id}"
47
+ end
48
+
49
+ def geometry_is_valid
50
+ if geog?
51
+ instance = self.class.unscoped.invalid.from("(SELECT ST_GeometryFromText('#{self.geog}') AS geog) #{self.class.table_name}").to_a.first
52
+ errors.add :geog, instance.invalid_geometry_message if instance
53
+ end
54
+ end
55
+
56
+ def sanitize_feature_type
57
+ self.feature_type = self.feature_type.to_s.strip.downcase
58
+ end
59
+ end
@@ -0,0 +1,3 @@
1
+ class SpatialCache < ActiveRecord::Base
2
+ belongs_to :spatial_model, :polymorphic => true
3
+ end
@@ -0,0 +1,4 @@
1
+ class SpatialProximity < ActiveRecord::Base
2
+ belongs_to :model_a, :polymorphic => true
3
+ belongs_to :model_b, :polymorphic => true
4
+ end
@@ -0,0 +1,56 @@
1
+ module SpatialFeatures
2
+ # Splits overlapping features into separate polygons at their areas of overlap, and returns an array of objects
3
+ # with kml for the overlapping area and a list of the record ids whose kml overlapped within that area
4
+ def self.venn_polygons(*scopes)
5
+ options = scopes.extract_options!
6
+ column = options[:simplified] ? 'geog_lowres' : 'geog'
7
+ scope = scopes.collect do |scope|
8
+ scope.klass.from(scope, scope.klass.table_name).joins(:features).where('features.feature_type = ?', 'polygon').select("features.#{column}::geometry AS the_geom").to_sql
9
+ end.join(' UNION ')
10
+
11
+ sql = "
12
+ SELECT scope.id, scope.type, ST_AsKML(geom) AS kml FROM ST_Dump((
13
+ SELECT ST_Polygonize(the_geom) AS the_geom FROM (
14
+
15
+ SELECT ST_Union(the_geom) AS the_geom FROM (
16
+
17
+ -- Handle Multigeometry
18
+ SELECT ST_ExteriorRing((ST_DumpRings(the_geom)).geom) AS the_geom
19
+ FROM (#{scope}) AS scope
20
+
21
+ ) AS exterior_lines
22
+
23
+ ) AS noded_lines
24
+ WHERE NOT ST_IsEmpty(the_geom) -- Ignore empty geometry from ST_Union if there are no polygons because polygonize will explode
25
+
26
+ )) AS venn_polygons
27
+ "
28
+
29
+ # Join with the original polygons so we can determine which original polygons each venn polygon came from
30
+ scope = scopes.collect do |scope|
31
+ scope.klass.from(scope, scope.klass.table_name).joins(:features).where('features.feature_type = ?', 'polygon').select("#{scope.klass.table_name}.id, features.spatial_model_type AS type, features.#{column}").to_sql
32
+ end.join(' UNION ')
33
+ sql <<
34
+ "INNER JOIN (#{scope}) AS scope
35
+ ON ST_Covers(scope.#{column}, ST_PointOnSurface(venn_polygons.geom)) -- Shrink the venn polygons so they don't share edges with the original polygons which could cause varying results due to tiny inaccuracy"
36
+
37
+ # Eager load the records for each venn polygon
38
+ eager_load_hash = Hash.new {|hash, key| hash[key] = []}
39
+ polygons = ActiveRecord::Base.connection.select_all(sql)
40
+ polygons.group_by{|row| row['type']}.each do |record_type, rows|
41
+ rows.each do |row|
42
+ eager_load_hash[record_type] << row['id']
43
+ end
44
+ end
45
+ eager_load_hash.each do |record_type, ids|
46
+ eager_load_hash[record_type] = record_type.constantize.find(ids)
47
+ end
48
+
49
+ # Instantiate objects to hold the kml and records for each venn polygon
50
+ polygons.group_by{|row| row['kml']}.collect do |kml, rows|
51
+ # Uniq on row id in case a single record had self intersecting multi geometry, which would cause it to appear duplicated on a single venn polygon
52
+ records = rows.uniq{|row| row.values_at('id', 'type') }.collect{|row| eager_load_hash.fetch(row['type']).detect{|record| record.id == row['id'].to_i } }
53
+ OpenStruct.new(:kml => kml, :records => records)
54
+ end
55
+ end
56
+ end
@@ -0,0 +1,3 @@
1
+ module SpatialFeatures
2
+ VERSION = "0.0.1"
3
+ end
@@ -0,0 +1,21 @@
1
+ # GEMS
2
+ require "activerecord-postgis-adapter"
3
+
4
+ # LIB
5
+ require 'spatial_features/caching'
6
+ require 'spatial_features/venn_polygons'
7
+ require 'spatial_features/has_spatial_features'
8
+
9
+ require 'spatial_features/import/arcgis_kmz_features'
10
+
11
+ require 'spatial_features/controller_helpers/spatial_extensions'
12
+
13
+ require 'spatial_features/models/feature'
14
+ require 'spatial_features/models/spatial_cache'
15
+ require 'spatial_features/models/spatial_proximity'
16
+
17
+ module SpatialFeatures
18
+ end
19
+
20
+ # Load the act method
21
+ ActiveRecord::Base.send :extend, SpatialFeatures::ActMethod
@@ -0,0 +1,4 @@
1
+ # desc "Explaining what the task does"
2
+ # task :spatial_features do
3
+ # # Task goes here
4
+ # end
@@ -0,0 +1,28 @@
1
+ == README
2
+
3
+ This README would normally document whatever steps are necessary to get the
4
+ application up and running.
5
+
6
+ Things you may want to cover:
7
+
8
+ * Ruby version
9
+
10
+ * System dependencies
11
+
12
+ * Configuration
13
+
14
+ * Database creation
15
+
16
+ * Database initialization
17
+
18
+ * How to run the test suite
19
+
20
+ * Services (job queues, cache servers, search engines, etc.)
21
+
22
+ * Deployment instructions
23
+
24
+ * ...
25
+
26
+
27
+ Please feel free to use a different markup language if you do not plan to run
28
+ <tt>rake doc:app</tt>.
@@ -0,0 +1,6 @@
1
+ # Add your own tasks in files placed in lib/tasks ending in .rake,
2
+ # for example lib/tasks/capistrano.rake, and they will automatically be available to Rake.
3
+
4
+ require File.expand_path('../config/application', __FILE__)
5
+
6
+ Rails.application.load_tasks
@@ -0,0 +1,13 @@
1
+ // This is a manifest file that'll be compiled into application.js, which will include all the files
2
+ // listed below.
3
+ //
4
+ // Any JavaScript/Coffee file within this directory, lib/assets/javascripts, vendor/assets/javascripts,
5
+ // or vendor/assets/javascripts of plugins, if any, can be referenced here using a relative path.
6
+ //
7
+ // It's not advisable to add code directly here, but if you do, it'll appear at the bottom of the
8
+ // compiled file.
9
+ //
10
+ // Read Sprockets README (https://github.com/sstephenson/sprockets#sprockets-directives) for details
11
+ // about supported directives.
12
+ //
13
+ //= require_tree .
@@ -0,0 +1,15 @@
1
+ /*
2
+ * This is a manifest file that'll be compiled into application.css, which will include all the files
3
+ * listed below.
4
+ *
5
+ * Any CSS and SCSS file within this directory, lib/assets/stylesheets, vendor/assets/stylesheets,
6
+ * or vendor/assets/stylesheets of plugins, if any, can be referenced here using a relative path.
7
+ *
8
+ * You're free to add application-wide styles to this file and they'll appear at the bottom of the
9
+ * compiled file so the styles you add here take precedence over styles defined in any styles
10
+ * defined in the other CSS/SCSS files in this directory. It is generally better to create a new
11
+ * file per style scope.
12
+ *
13
+ *= require_tree .
14
+ *= require_self
15
+ */
@@ -0,0 +1,5 @@
1
+ class ApplicationController < ActionController::Base
2
+ # Prevent CSRF attacks by raising an exception.
3
+ # For APIs, you may want to use :null_session instead.
4
+ protect_from_forgery with: :exception
5
+ end
@@ -0,0 +1,2 @@
1
+ module ApplicationHelper
2
+ end
@@ -0,0 +1,14 @@
1
+ <!DOCTYPE html>
2
+ <html>
3
+ <head>
4
+ <title>Dummy</title>
5
+ <%= stylesheet_link_tag 'application', media: 'all', 'data-turbolinks-track' => true %>
6
+ <%= javascript_include_tag 'application', 'data-turbolinks-track' => true %>
7
+ <%= csrf_meta_tags %>
8
+ </head>
9
+ <body>
10
+
11
+ <%= yield %>
12
+
13
+ </body>
14
+ </html>