geo_graf 1.0.2

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: b2efd6937e6852e7b2fc4d422629843b005acbf8
+  data.tar.gz: 3a832e48ff664b270f751f33fd070098189301cd
+SHA512:
+  metadata.gz: 47a2b6cfa2e69e6a60985f654fc52f2fea73db5da92d1364500c99af7233ffea763044271009f1bb3481f351add9b69fe5745ef58920a6c54327f9dd22d3bc04
+  data.tar.gz: c2f0a08bf4bb9f17a7bfb8a46ff4f5a4e59eaf840762e721349b4a377c9fb9f8d2c947b4d5aed429b28980915785557aff8dd9c8990bcd0fe44603f58bb4040e

data/.gitignore

@@ -0,0 +1,2 @@
+tmp/
+*.gem

data/.travis.yml

@@ -0,0 +1,6 @@
+before_install: sudo apt-get install libgeos++-dev
+
+language: ruby
+rvm:
+  - 2.0.0
+  - ruby-head

data/Gemfile

@@ -0,0 +1,4 @@
+source 'https://rubygems.org'
+
+# Specify your gem's dependencies in wimdu-geograf.gemspec
+gemspec

data/Gemfile.lock

@@ -0,0 +1,59 @@
+PATH
+  remote: .
+  specs:
+    geo_graf (1.0.2)
+      rgeo
+
+GEM
+  remote: https://rubygems.org/
+  specs:
+    celluloid (0.15.2)
+      timers (~> 1.1.0)
+    coderay (1.0.9)
+    diff-lcs (1.2.4)
+    ffi (1.9.0)
+    formatador (0.2.4)
+    guard (2.0.3)
+      formatador (>= 0.2.4)
+      listen (~> 2.0)
+      lumberjack (~> 1.0)
+      pry (>= 0.9.12)
+      thor (>= 0.18.1)
+    guard-rspec (3.1.0)
+      guard (>= 1.8)
+      rspec (~> 2.13)
+    listen (2.0.3)
+      celluloid (>= 0.15.2)
+      rb-fsevent (>= 0.9.3)
+      rb-inotify (>= 0.9)
+    lumberjack (1.0.4)
+    method_source (0.8.2)
+    pry (0.9.12.2)
+      coderay (~> 1.0.5)
+      method_source (~> 0.8)
+      slop (~> 3.4)
+    rake (10.1.0)
+    rb-fsevent (0.9.3)
+    rb-inotify (0.9.2)
+      ffi (>= 0.5.0)
+    rgeo (0.3.20)
+    rspec (2.14.1)
+      rspec-core (~> 2.14.0)
+      rspec-expectations (~> 2.14.0)
+      rspec-mocks (~> 2.14.0)
+    rspec-core (2.14.5)
+    rspec-expectations (2.14.3)
+      diff-lcs (>= 1.1.3, < 2.0)
+    rspec-mocks (2.14.3)
+    slop (3.4.6)
+    thor (0.18.1)
+    timers (1.1.0)
+
+PLATFORMS
+  ruby
+
+DEPENDENCIES
+  bundler (~> 1.3)
+  geo_graf!
+  guard-rspec
+  rake

data/Guardfile

@@ -0,0 +1,4 @@
+guard 'rspec', keep_failed: true, cli: '--format documentation' do
+  watch(%r{^spec/.+_spec\.rb$})
+  watch(%r{^lib/(.+)\.rb$})     { |m| "spec/#{m[1]}_spec.rb" }
+end

data/LICENSE.txt

@@ -0,0 +1,22 @@
+Copyright (c) 2013 Wimdu GmbH
+
+MIT License
+
+Permission is hereby granted, free of charge, to any person obtaining
+a copy of this software and associated documentation files (the
+"Software"), to deal in the Software without restriction, including
+without limitation the rights to use, copy, modify, merge, publish,
+distribute, sublicense, and/or sell copies of the Software, and to
+permit persons to whom the Software is furnished to do so, subject to
+the following conditions:
+
+The above copyright notice and this permission notice shall be
+included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
+EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
+MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
+NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
+LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
+OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
+WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

data/README.md

@@ -0,0 +1,90 @@
+# GeoGraf
+
+[![Build Status](https://travis-ci.org/wimdu/geo_graf.png?branch=master)](https://travis-ci.org/wimdu/geo_graf)
+
+## Overview
+
+This gem computes the relation between multiple polygons (do they overlap and
+by how much?), such as Wimdu's GeoLocations.
+
+GeoGraf uses [RGeo](http://dazuma.github.io/rgeo/rdoc/) to perform the
+calculations. It uses the
+[Mercator projection](http://en.wikipedia.org/wiki/Mercator_projection), which
+projects the surface of earth onto a flat surface. This is inaccurate, as
+areas appear bigger when they are further away from the equator (have a look at
+Greenland and the DR Congo on Google Maps, they both are roughly 2 million
+square kilometers big). However, the assumption is that the error of areas that
+relate to each other is rather small, as the overlapping area is at the same
+latitude. And we're not interested in the correct area calculation, we just need
+to know if they overlap, which is the bigger one and roughly how much they
+overlap.
+
+The input is expected to be an array of polygon descriptions. Each polygon
+description has to be a hash with an `:id` key and `:polygon_coords`, which
+itself is an array of coordinates describing the polygon. Each coordinate is
+an array of it's latitude and longitude (in this order).
+
+Example:
+
+```ruby
+require 'geo_graf'
+
+GeoGraf.intersections_for(
+  [
+    {
+      id: 1,
+      polygon_coords: [[-2, -2], [2, -2], [2, 2], [-2, 2], [-2, -2]]
+    },
+    {
+      id: 2,
+      polygon_coords: [[-3, -3], [-1, -3], [-1, -1], [-3, -1], [-3, -3]]
+    }
+  ]
+)
+```
+
+The output is an array of relation descriptions. Each description is a hash in
+the following format:
+  * `:id`: The ID of the polygon. It's always the ID of the smaller polygon
+    (by area).
+  * `:contained_area_percentage`: How much of the area of the smaller polygon
+    is contained in the larger one.
+  * `:container_id`: The ID of the larger polygon.
+
+If two polygons don't overlap, no relation exists and there will be no
+description returned for these two.
+
+For the example above, the following result will be returned:
+
+```ruby
+[
+  {
+    id: 2,
+    contained_area_percentage: 25,
+    container_id: 1
+  }
+]
+```
+
+This means that 25% of the area of polygon 2 overlap with the area of polygon
+1.
+
+## Installation
+
+GeoGraf depends on RGeo with [GEOS](http://trac.osgeo.org/geos/) >= 3.3.3. You
+need to install it before installing RGeo.
+
+### MacOs with homebrew
+
+```sh
+brew install geos
+```
+
+### Ubuntu/Debian
+
+```sh
+sudo apt-get install libgeos++-dev
+```
+
+## License
+Copyright (c) 2013 Wimdu GmbH (MIT License). See LICENSE.txt for details.

data/Rakefile

@@ -0,0 +1,9 @@
+require 'bundler/gem_tasks'
+require 'bundler/setup'
+require 'rspec/core/rake_task'
+
+RSpec::Core::RakeTask.new do |t|
+  t.rspec_opts = ['--color', '--format doc', '--order rand']
+end
+
+task default: :spec

data/geo_graf.gemspec

@@ -0,0 +1,20 @@
+Gem::Specification.new do |spec|
+  spec.name          = 'geo_graf'
+  spec.version       = '1.0.2'
+  spec.date          = Time.now.strftime('%Y-%m-%d')
+  spec.authors       = ['Marek Nowak', 'Johannes Barre']
+  spec.email         = ['marek.nowak@wimdu.com', 'johannes.barre@wimdu.com']
+  spec.description   = %q{Calculates the relations between overlapping polygons}
+  spec.summary       = %q{Calculates the relations between overlapping polygons}
+  spec.homepage      = 'https://github.com/wimdu/geo_graf'
+
+  spec.files         = `git ls-files`.split($/)
+  spec.executables   = spec.files.grep(%r{^bin/}) { |f| File.basename(f) }
+  spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
+  spec.require_paths = ['lib']
+
+  spec.add_dependency 'rgeo'
+  spec.add_development_dependency 'bundler', '~> 1.3'
+  spec.add_development_dependency 'guard-rspec'
+  spec.add_development_dependency 'rake'
+end

data/lib/geo_graf.rb

@@ -0,0 +1,10 @@
+require_relative 'geo_graf/intersection_calculator'
+
+module GeoGraf
+  class << self
+
+    def intersections_for(input_geodata)
+      IntersectionCalculator.new(input_geodata).intersections
+    end
+  end
+end

data/lib/geo_graf/intersection_calculator.rb

@@ -0,0 +1,45 @@
+require 'rgeo'
+require 'geo_graf/polygon'
+
+module GeoGraf
+  class GeosNotInstalledError < StandardError; end
+
+  class IntersectionCalculator
+    def initialize(input_geodata)
+      @geodata ||= input_geodata
+        .map { |geodatum| {id: geodatum[:id], polygon: polygon_from_coords(geodatum[:polygon_coords])} }
+        .sort_by { |geodatum| geodatum[:polygon].area } # always have smaller first
+    end
+
+    def intersections
+      intersections = []
+
+      geodata.combination(2) do |smaller, bigger|
+        intersection_area = smaller[:polygon].intersection_area(bigger[:polygon])
+
+        unless intersection_area.zero?
+          intersections << {
+            id: smaller[:id],
+            contained_area_percentage: (intersection_area.to_f / smaller[:polygon].area * 100).round,
+            container_id: bigger[:id]
+          }
+        end
+      end
+
+      intersections
+    end
+
+    private
+
+    attr_reader :geodata
+
+    def polygon_from_coords(coords)
+      Polygon.new(coords, rgeo_factory)
+    end
+
+    def rgeo_factory
+      raise(GeosNotInstalledError, 'The Geos library needs to be installed (see README.md).') unless RGeo::Geos.supported?
+      @factory ||= RGeo::Geographic.simple_mercator_factory
+    end
+  end
+end

data/lib/geo_graf/polygon.rb

@@ -0,0 +1,83 @@
+module GeoGraf
+  class Polygon
+    def initialize(coordinates, factory)
+      @rgeo_polygons = [create_polygon(coordinates, factory)]
+
+      create_shadow_polygon_if_required(coordinates, factory)
+    end
+
+    def intersection_area(other)
+      intersection = intersection(other)
+
+      if intersection.respond_to?(:area)
+        intersection.area
+      elsif intersection.respond_to?(:map)
+        sum_of_areas_in(intersection)
+      else
+        0.0
+      end
+    end
+
+    def area
+      rgeo_polygons.first.area
+    end
+
+    protected
+
+    attr_reader :rgeo_polygons
+
+    private
+
+    def intersection(other)
+      intersection = nil
+
+      # If you wonder why for a seemingly simple intersection of two
+      # polygons we are actually intersecting multiple rgeo_polygons
+      # with other's rgeo_polygons, check out the comment above
+      # #create_shadow_polygon_if_required. Peace.
+
+      rgeo_polygons.each do |polygon|
+        other.rgeo_polygons.each do |other_polygon|
+          intersection = polygon.intersection(other_polygon)
+          return intersection if intersection && !intersection.is_empty?
+        end
+      end 
+      intersection
+    end
+
+    def create_polygon(coordinates, factory, longitude_offset: 0)
+      points = coordinates.map { |c| factory.point(c[1] + longitude_offset, c[0]) }
+      ring = factory.linear_ring(points)
+      factory.polygon(ring)
+    end
+
+    def sum_of_areas_in(collection)
+      collection
+      .select { |i| i.respond_to?(:area) }
+      .map(&:area)
+      .inject(0, :+)
+    end
+
+    # Ah, the shadow polygon... We use it because RGeo is kind of retarded
+    # and doesn't deal with situations in which shapes wrap around
+    # the 180th meridian correctly. Therefore we create a fake shape with
+    # an offset of either 360 or -360 degrees to ensure that we catch
+    # the intersection if there is one. And, the name is awesome, right?
+
+    def create_shadow_polygon_if_required(coordinates, factory)
+      if line_180(factory).intersects?(rgeo_polygons.first)
+        rgeo_polygons << create_polygon(coordinates, factory, longitude_offset: -360)
+      elsif line_minus_180(factory).intersects?(rgeo_polygons.first)
+        rgeo_polygons << create_polygon(coordinates, factory, longitude_offset: 360)
+      end
+    end
+
+    def line_180(factory)
+      factory.line(factory.point(180, 90), factory.point(180, -90))
+    end
+
+    def line_minus_180(factory)
+      factory.line(factory.point(-180, 90), factory.point(-180, -90))
+    end
+  end
+end

data/spec/geo_graf/intersection_calculator_spec.rb

@@ -0,0 +1,215 @@
+$: << '../lib'
+
+require 'geo_graf/intersection_calculator'
+
+describe GeoGraf::IntersectionCalculator do
+  describe "#intersections" do
+    subject(:intersections) { described_class.new(input).intersections }
+
+    context "given no shape" do
+      let(:input) { [] }
+
+      it "returns no intersections" do
+        expect(intersections).to be_empty
+      end
+    end
+
+    context "given only one shape" do
+      let(:input) {
+        [{
+          id: 1,
+          polygon_coords: [[-2, -2], [2, -2], [-2, 1], [-2, -2]]
+        }]
+      }
+
+      it "returns no intersections" do
+        expect(intersections).to be_empty
+      end
+    end
+
+    context "given two nonintersecting shapes" do
+      let(:input) {
+        [
+          {
+            id: 1,
+            polygon_coords: [[-2, -2], [2, -2], [-2, 1], [-2, -2]]
+          },
+          {
+            id: 2,
+            polygon_coords: [[100, 100], [102, 100], [100, 101], [100, 100]]
+          }
+        ]
+      }
+
+      it "returns no intersections" do
+        expect(intersections).to be_empty
+      end
+    end
+
+    context "given two identical shapes" do
+      let(:input) {
+        [
+          {
+            id: 1,
+            polygon_coords: [[-2, -2], [2, -2], [-2, 1], [-2, -2]]
+          },
+          {
+            id: 2,
+            polygon_coords: [[-2, -2], [2, -2], [-2, 1], [-2, -2]]
+          }
+        ]
+      }
+
+      it "returns an intersection of the second one containing the first one 100%" do
+        expect(intersections).to match_array([{id: 1, contained_area_percentage: 100, container_id: 2}])
+      end
+    end
+
+    context "given one square inside the other" do
+      let(:input) {
+        [
+          {
+            id: 1,
+            polygon_coords: [[-2, -2], [2, -2], [2, 2], [-2, 2], [-2, -2]]
+          },
+          {
+            id: 2,
+            polygon_coords: [[-2, -2], [0, -2], [0, 0], [-2, 0], [-2, -2]]
+          }
+        ]
+      }
+
+      it "returns an intersection of the first one containing the second one 100%" do
+        expect(intersections).to match_array([{id: 2, contained_area_percentage: 100, container_id: 1}])
+      end
+    end
+
+    context "given one square intersecting the other 25%" do
+      let(:input) {
+        [
+          {
+            id: 1,
+            polygon_coords: [[-2, -2], [2, -2], [2, 2], [-2, 2], [-2, -2]]
+          },
+          {
+            id: 2,
+            polygon_coords: [[-3, -3], [-1, -3], [-1, -1], [-3, -1], [-3, -3]]
+          }
+        ]
+      }
+
+      it "returns an intersection of the first one containing the second one 25%" do
+        expect(intersections).to match_array([{id: 2, contained_area_percentage: 25, container_id: 1}])
+      end
+    end
+
+    context "given one square intersecting the other square of the same size 50%" do
+      let(:input) {
+        [
+          {
+            id: 1,
+            polygon_coords: [[-2, -2], [2, -2], [2, 2], [-2, 2], [-2, -2]]
+          },
+          {
+            id: 2,
+            polygon_coords: [[-4, -2], [0, -2], [0, 2], [-4, 2], [-4, -2]]
+          }
+        ]
+      }
+
+      it "returns an intersection of the second one containing the first one 50%" do
+        expect(intersections).to match_array([{id: 1, contained_area_percentage: 50, container_id: 2}])
+      end
+    end
+
+    context "given multiple shapes of which some intersect" do
+      let(:input) {
+        [
+          {
+            id: 1,
+            polygon_coords: [[-2, -2], [2, -2], [2, 2], [-2, 2], [-2, -2]]
+          },
+          {
+            id: 2,
+            polygon_coords: [[-1, 1], [-4, 1], [-4, -2], [-1, -2], [-1, 1]]
+          },
+          {
+            id: 3,
+            polygon_coords: [[-3, -3], [-1, -3], [-1, -1], [-3, -1], [-3, -3]]
+          },
+          {
+            id: 4,
+            polygon_coords: [[-1, 8], [-3, 8], [-3, 6], [-1, 6], [-1, 8]]
+          },
+          {
+            id: 5,
+            polygon_coords: [[11, 6], [11, 0], [2, 6], [11, 6]]
+          },
+          {
+            id: 6,
+            polygon_coords: [[5, 4], [8, 4], [8, 2], [5, 2], [5, 4]]
+          },
+          {
+            id: 7,
+            polygon_coords: [[3, -3], [1, -3], [1, -1], [3, -1], [3, -3]]
+          }
+        ]
+      }
+
+      it "returns 4 intersections" do
+        expect(intersections).to match_array([
+          {id: 2, contained_area_percentage: 33, container_id: 1},
+          {id: 3, contained_area_percentage: 50, container_id: 2},
+          {id: 3, contained_area_percentage: 25, container_id: 1},
+          {id: 6, contained_area_percentage: 49, container_id: 5}, # 49% because of it's not a plane but a sphere
+          {id: 7, contained_area_percentage: 25, container_id: 1}
+        ])
+      end
+    end
+
+    context "given multiple shapes on the both sides of 180 degrees line of which some intersect" do
+      let(:input) {
+        [
+          {
+            id: 1,
+            polygon_coords: [[-2, 178], [2, 178], [2, -178], [-2, -178], [-2, 178]]
+          },
+          {
+            id: 2,
+            polygon_coords: [[-1, -179], [-4, -179], [-4, 178], [-1, 178], [-1, -179]]
+          },
+          {
+            id: 3,
+            polygon_coords: [[-3, 177], [-1, 177], [-1, 179], [-3, 179], [-3, 177]]
+          },
+          {
+            id: 4,
+            polygon_coords: [[-1, -172], [-3, -172], [-3, -174], [-1, -174], [-1, -172]]
+          },
+          {
+            id: 5,
+            polygon_coords: [[11, -174], [11, 180], [2, -174], [11, -174]]
+          },
+          {
+            id: 6,
+            polygon_coords: [[5, -176], [8, -176], [8, -178], [5, -178], [5, -176]]
+          },
+          {
+            id: 7,
+            polygon_coords: [[3, 177], [1, 177], [1, 179], [3, 179], [3, 177]]
+          }
+        ]
+      }
+
+      it "returns 4 intersections" do
+        expect(intersections).to match_array([
+          {id: 2, contained_area_percentage: 33, container_id: 1},
+          {id: 3, contained_area_percentage: 50, container_id: 2},
+          {id: 3, contained_area_percentage: 25, container_id: 1},
+          {id: 6, contained_area_percentage: 49, container_id: 5}, # 49% because of it's not a plane but a sphere
+          {id: 7, contained_area_percentage: 25, container_id: 1}
+        ])
+      end
+    end
+  end
+end

data/spec/geo_graf/polygon_spec.rb

@@ -0,0 +1,146 @@
+$: << '../lib'
+
+require 'geo_graf/polygon'
+
+describe GeoGraf::Polygon do
+  let(:polygon) { described_class.new(input, factory) }
+  let(:factory) { double("Factory", polygon: rgeo_polygon, point: rgeo_point, linear_ring: rgeo_linear_ring, line: line_180) }
+  let(:rgeo_polygon) { double("RGeoPolygon") } 
+  let(:rgeo_point) { double("RGeoPoint") }
+  let(:rgeo_linear_ring) { double("LinearRing") }
+  let(:line_180) { double("Line", intersects?: false) }
+  let(:input) { [[-2, -2], [2, -2], [-2, 1], [-2, -2]] }
+
+  describe ".new" do
+    subject(:create_new_polygon) { polygon }
+
+    context "when the polygon doesn't span across 180th meridian" do
+      it "creates a point for each coordinate" do
+        factory.should_receive(:point).exactly(8).times # 4 + 4 used for creation of 180 and -180 lines
+
+        create_new_polygon
+      end
+
+      it "creates a linear ring out of the points" do
+        factory.should_receive(:linear_ring).with([rgeo_point] * 4)
+
+        create_new_polygon
+      end
+
+      it "creates the polygon using the linear ring" do
+        factory.should_receive(:polygon).with(rgeo_linear_ring)
+
+        create_new_polygon
+      end
+    end
+
+    context "when the polygon spans across 180th meridian" do
+      let(:line_180) { double("Line", intersects?: true) }
+
+      it "creates two polygons using the linear rings" do
+        factory.should_receive(:polygon).with(rgeo_linear_ring).twice
+
+        create_new_polygon
+      end
+    end
+  end
+
+  describe "#area" do
+    subject(:get_area) { polygon.area }
+
+    it "delegates to the internal RGeo polygon" do
+      rgeo_polygon.should_receive(:area)
+
+      get_area
+    end
+  end
+
+  describe "#intersection_area" do
+    let(:other_polygon) { double('OtherPolygon', rgeo_polygons: [others_rgeo_polygon_1]) }
+    let(:others_rgeo_polygon_1) { double('OtherRgeoPolygon1') }
+    let(:intersection) { double('Intersection', is_empty?: false, area: 1234.5) }
+    let(:no_intersection) { double('NoIntersection', is_empty?: true) }
+
+    context "when the polygon doesn't span across 180th meridian" do
+      context "when it intersects the other polygon" do
+        before do
+          rgeo_polygon.stub(:intersection).with(others_rgeo_polygon_1).and_return(intersection)
+        end
+
+        it "returns the area of the intersection" do
+          expect(polygon.intersection_area(other_polygon)).to eq(1234.5)
+        end
+
+        context "when the intersection spans over multiple polygons and points" do
+          let(:intersection) { [double('Point'), double('Polygon', area: 123.67), double('Polygon', area: 3.12)] }
+
+          before do
+            intersection.stub(:is_empty?).and_return(false)
+          end
+
+          it "returns the sum of areas of the intersection" do
+            expect(polygon.intersection_area(other_polygon)).to eq(126.79)
+          end
+        end
+      end
+
+      context "when it doesn't intersect the other polygon" do
+        before do
+          rgeo_polygon.stub(:intersection).with(others_rgeo_polygon_1).and_return(no_intersection)
+        end
+
+        it "returns 0.0" do
+          expect(polygon.intersection_area(other_polygon)).to eq(0.0)
+        end
+      end
+    end
+
+    context "when the polygon spans across 180th meridian" do
+      let(:other_polygon) {
+        double('OtherPolygon', rgeo_polygons: [others_rgeo_polygon_1, others_rgeo_polygon_2])
+      }
+      let(:others_rgeo_polygon_2) { double('OtherRgeoPolygon2') }
+      let(:shadow_rgeo_polygon) { double("ShadowRgeoPolygon") } 
+      let(:line_180) { double("Line", intersects?: true) }
+
+      before do
+        factory.stub(:polygon).and_return(rgeo_polygon, shadow_rgeo_polygon)
+        rgeo_polygon.stub(:intersection).with(others_rgeo_polygon_1).and_return(no_intersection)
+        rgeo_polygon.stub(:intersection).with(others_rgeo_polygon_2).and_return(no_intersection)
+        shadow_rgeo_polygon.stub(:intersection).with(others_rgeo_polygon_2).and_return(no_intersection)
+      end
+
+      context "when it intersects the other polygon" do
+        before do
+          shadow_rgeo_polygon.stub(:intersection).with(others_rgeo_polygon_1).and_return(intersection)
+        end
+
+        it "returns the area of the intersection" do
+          expect(polygon.intersection_area(other_polygon)).to eq(1234.5)
+        end
+
+        context "when the intersection spans over multiple polygon" do
+          let(:intersection) { [double('Point'), double('Polygon', area: 123.67), double('Polygon', area: 3.12)] }
+
+          before do
+            intersection.stub(:is_empty?).and_return(false)
+          end
+
+          it "returns the sum of areas of the intersection" do
+            expect(polygon.intersection_area(other_polygon)).to eq(126.79)
+          end
+        end
+      end
+
+      context "when it doesn't intersect the other polygon" do
+        before do
+          shadow_rgeo_polygon.stub(:intersection).with(others_rgeo_polygon_1).and_return(no_intersection)
+        end
+
+        it "returns 0.0" do
+          expect(polygon.intersection_area(other_polygon)).to eq(0.0)
+        end
+      end
+    end
+  end
+end

data/spec/geo_graf_spec.rb

@@ -0,0 +1,22 @@
+$: << '../lib'
+
+require 'geo_graf'
+
+describe GeoGraf do
+  describe "#intersections_for" do
+    subject(:intersections_for_input) { described_class.intersections_for(input) }
+    
+    let(:calculator) { double('Calculator') }
+    let(:input) { double('Input') }
+    let(:intersection_calculator_output) { double('Output') }
+
+    before do
+      allow(GeoGraf::IntersectionCalculator).to receive(:new).with(input).and_return(calculator)
+      allow(calculator).to receive(:intersections).and_return(intersection_calculator_output)
+    end
+
+    it "runs intersections on a new instance of IntersectionCalculator" do
+      expect(intersections_for_input).to eq(intersection_calculator_output)
+    end
+  end
+end

metadata

@@ -0,0 +1,119 @@
+--- !ruby/object:Gem::Specification
+name: geo_graf
+version: !ruby/object:Gem::Version
+  version: 1.0.2
+platform: ruby
+authors:
+- Marek Nowak
+- Johannes Barre
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2013-10-11 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: rgeo
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '1.3'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ~>
+      - !ruby/object:Gem::Version
+        version: '1.3'
+- !ruby/object:Gem::Dependency
+  name: guard-rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - '>='
+      - !ruby/object:Gem::Version
+        version: '0'
+description: Calculates the relations between overlapping polygons
+email:
+- marek.nowak@wimdu.com
+- johannes.barre@wimdu.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- .gitignore
+- .travis.yml
+- Gemfile
+- Gemfile.lock
+- Guardfile
+- LICENSE.txt
+- README.md
+- Rakefile
+- geo_graf.gemspec
+- lib/geo_graf.rb
+- lib/geo_graf/intersection_calculator.rb
+- lib/geo_graf/polygon.rb
+- spec/geo_graf/intersection_calculator_spec.rb
+- spec/geo_graf/polygon_spec.rb
+- spec/geo_graf_spec.rb
+homepage: https://github.com/wimdu/geo_graf
+licenses: []
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - '>='
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.0.3
+signing_key: 
+specification_version: 4
+summary: Calculates the relations between overlapping polygons
+test_files:
+- spec/geo_graf/intersection_calculator_spec.rb
+- spec/geo_graf/polygon_spec.rb
+- spec/geo_graf_spec.rb