georuby-ext 0.0.1 → 0.0.2

data/spec/lib/georuby/multi_polygon_spec.rb

@@ -0,0 +1,69 @@
+require "spec_helper"
+
+describe  GeoRuby::SimpleFeatures::MultiPolygon do  
+  
+  describe "#to_wgs84" do
+    let(:multi_polygon_google) { geometry "MULTIPOLYGON(((0.0 -7.08115455161362e-10,0.0 111325.142866385,111319.490793272 -7.08115455161362e-10,0.0 -7.08115455161362e-10)))", 900913 }
+    let(:multi_polygon_wgs84) { geometry "MULTIPOLYGON(((0 0,0 1,1 0,0 0)))" } 
+    
+    it "should return a polygon in wgs84 coordinates" do
+      [multi_polygon_google.to_wgs84.points, multi_polygon_wgs84.points].transpose.each{ |p1,p2|
+        p1.x.should be_within(0.0000001).of(p2.x) 
+        p1.y.should be_within(0.0000001).of(p2.y)
+      }
+    end
+
+    it "should return a multi_polygon with wgs84 srid" do
+      multi_polygon_google.to_wgs84.srid.should == multi_polygon_wgs84.srid
+    end
+  end
+
+  describe "#polygons" do
+    let(:geo_multi_polygon) { geometry "MULTIPOLYGON(((0 0,0 1,1 1,1 0)), ((0 0,0 1,1 1,1 0)))" } 
+    let(:geo_polygon){ polygon("(0 0,0 1,1 1,1 0)") } 
+
+    it "should return an array of polygons" do
+      geo_multi_polygon.polygons.should == [geo_polygon, geo_polygon]
+    end
+  end
+
+  describe "#difference" do    
+    let(:g_multi_polygon) { geometry("MULTIPOLYGON (((0 0, 0 5, 5 5, 5 0, 0 0)), ((6 1, 6 2, 7 2, 7 1, 6 1)))") } 
+    let(:empty_geometry_collection) { GeoRuby::SimpleFeatures::GeometryCollection.from_geometries([], 4326) }
+
+    it "should return an empty geometry collection if multi polygons are the same" do
+      g_multi_polygon.difference(g_multi_polygon).should == empty_geometry_collection
+    end
+
+    it "should return an rgeo polygon if there's only one polygon left' " do
+      g_multi_polygon2 = geometry("MULTIPOLYGON (((0 0, 0 5, 5 5, 5 0, 0 0)))") 
+      g_multi_polygon.difference(g_multi_polygon2).should == geometry("POLYGON ((6 1, 6 2, 7 2, 7 1, 6 1))") 
+    end
+
+    it "should return a multi polygon if there's more than one polygon left" do
+      g_multi_polygon2 = geometry("MULTIPOLYGON (((0.25 0.25, 0.25 0.5, 0.75 0.5, 0.75 0.25, 0.25 0.25)))") 
+      g_multi_polygon.difference(g_multi_polygon2).should == geometry("MULTIPOLYGON(((0 0, 0 5, 5 5, 5 0, 0 0),(0.25 0.25,0.75 0.25,0.75 0.5,0.25 0.5,0.25 0.25)),((6 1, 6 2, 7 2, 7 1, 6 1)))")
+    end
+    
+    it "should return the same polygon if the polygon to substract isn't inside the polygon" do
+      g_multi_polygon2 = geometry("MULTIPOLYGON (((0 0, 0 -1, -1 -1, -1 0, 0 0)))") 
+      g_multi_polygon.difference(g_multi_polygon2).should == geometry("MULTIPOLYGON (((0 0, 0 5, 5 5, 5 0, 0 0)), ((6 1, 6 2, 7 2, 7 1, 6 1)))")
+    end
+
+    it "should return a multi polygon with 2 holes" do
+      g_multi_polygon2 = geometry("MULTIPOLYGON ( ((1 1, 2 1, 2 2, 1 2, 1 1)), ((3 3, 4 3, 4 4, 3 4, 3 3)) )") 
+      g_multi_polygon.difference(g_multi_polygon2).should == geometry("MULTIPOLYGON ( ((0 0, 0 5, 5 5, 5 0, 0 0), (1 1, 2 1, 2 2, 1 2, 1 1), (3 3, 4 3, 4 4, 3 4, 3 3)), ((6 1, 6 2, 7 2, 7 1, 6 1)) )")
+    end
+
+  end
+
+  describe "#to_rgeo" do
+    let(:g_multi_polygon) { geometry("MULTIPOLYGON(((0.0 0.0,0.0 1.0,1.0 1.0,1.0 0.0,0.0 0.0),(0.25 0.25,0.75 0.25,0.75 0.5,0.25 0.5,0.25 0.25)),((0.0 2.0,0.0 3.0,3.0 3.0,2.0 2.0,0.0 2.0)))") } 
+    let(:rg_multi_polygon) { rgeometry("MULTIPOLYGON(((0.0 0.0,0.0 1.0,1.0 1.0,1.0 0.0,0.0 0.0),(0.25 0.25,0.75 0.25,0.75 0.5,0.25 0.5,0.25 0.25)),((0.0 2.0,0.0 3.0,3.0 3.0,2.0 2.0,0.0 2.0)))") }   
+
+    it "should return an rgeo multi polygons" do     
+      g_multi_polygon.to_rgeo.should == rg_multi_polygon
+    end
+  end
+  
+end

data/spec/lib/georuby/point_spec.rb

@@ -0,0 +1,248 @@
+require 'spec_helper'
+
+describe GeoRuby::SimpleFeatures::Point do
+
+  subject { point 1.0, 2.0 }
+
+  describe "#==" do
+
+    let(:other) { mock :other, :lat => nil, :lng => nil }
+    
+    it "should return false if other is nil" do
+      subject.should_not == nil
+    end
+
+    it "should return true if spherical_distance to other is smaller than 10e-3 (< 1m)" do
+      subject.stub :spherical_distance => 10e-4
+      subject.should == other
+    end
+
+    it "should return false if spherical distance to other is greater or equals to 10e-3 (>= 1m)" do
+      subject.stub :spherical_distance => 10e-3
+      subject.should_not == other
+    end
+
+    it "should match the same point into 2 different srid" do
+      subject.to_google.should == subject
+    end
+
+  end
+
+  describe "to_rgeo" do
+
+    it "should have the same information" do
+      subject.to_rgeo.should have_same(:x, :y, :srid).than(subject)
+    end
+
+  end
+
+  describe ".centroid" do
+
+    def centroid(points)
+      points = points(points) if String === points
+      GeoRuby::SimpleFeatures::Point.centroid points
+    end
+    
+    it "should return nil if no points in input" do
+      centroid([]).should be_nil
+    end
+    
+    it "should return point in input if only one" do
+      centroid("0 0").should == point(0,0)
+    end
+
+    it "should return middle if two points in input" do
+      centroid("0 0,1 0").should == point(0.5,0)
+    end
+    
+    it "should return centroid of given points" do
+      centroid("0 0,0 2,2 2,2 0").should == point(1,1)
+    end
+
+    it "should have the same srid than points" do
+      google_point = point.to_google
+      centroid([google_point]).srid.should == google_point.srid
+    end
+    
+  end
+
+  describe "#to_proj4" do
+    
+    it "should return a Proj4::Point" do
+      subject.to_proj4.should be_instance_of(Proj4::Point)
+    end
+
+    it "should use the given ratio (if specified)" do
+      subject.to_proj4(10).x.should == subject.x * 10
+      subject.to_proj4(10).y.should == subject.y * 10
+    end
+
+    context "when Point is wgs84" do
+
+      it "should transform x to radians (using ratio #{Proj4::DEG_TO_RAD})" do
+        subject.to_proj4.x.should == subject.x * Proj4::DEG_TO_RAD
+      end
+
+      it "should transform y to radians (using ratio #{Proj4::DEG_TO_RAD})" do
+        subject.to_proj4.y.should == subject.y * Proj4::DEG_TO_RAD
+      end
+
+    end
+
+    it "should have the same z" do
+      subject.to_proj4.z.should == nil
+    end
+
+  end
+
+  describe ".from_pro4j" do
+
+    let(:proj4_point) { subject.to_proj4 }
+    let(:srid) { subject.srid }
+
+    def from_pro4j(point, srid = srid)
+      GeoRuby::SimpleFeatures::Point.from_pro4j point, srid
+    end
+
+    it "should transform x to degres" do
+      from_pro4j(proj4_point).x.should == proj4_point.x * Proj4::RAD_TO_DEG
+    end
+
+    it "should transform y to degres" do
+      from_pro4j(proj4_point).y.should == proj4_point.y * Proj4::RAD_TO_DEG
+    end
+
+    it "should have the specified srid" do
+      from_pro4j(proj4_point, srid).srid.should == srid
+    end
+    
+  end
+
+  it "should be used as Hash key" do
+    hash = { point(0,0) => :found }
+    hash[point(0,0)].should == :found
+  end
+
+  describe "#eql?" do
+    
+    it "should compare x, y, z and srid" do
+      point(0,0).should eq(point(0,0))
+    end
+
+  end
+
+  describe "hash" do
+
+    it "should use x, y, z and srid" do
+      point(0,0).hash.should == point(0,0).hash
+    end
+    
+  end
+
+  describe "to_s" do
+    
+    it "should contain 'y,x'" do
+      point(1,2).to_s.should == "2,1"
+    end
+
+  end
+
+  describe "#bounding_box" do
+
+    it "should be itself twice" do
+      subject.bounding_box.should == [subject, subject]
+    end
+
+    it "should return 3 points with_z" do
+      subject.with_z = true
+      subject.bounding_box.should == [subject] * 3
+    end
+    
+  end
+
+  describe "#projection" do
+
+    let(:projection) { mock :projection }
+    
+    it "should return projection associated to srid" do
+      Proj4::Projection.should_receive(:for_srid).with(subject.srid).and_return(projection)
+      subject.projection.should == projection
+    end
+
+  end
+
+  describe "#project_to" do
+
+    let(:tour_eiffel_in_wgs84) { point 2.2946, 48.8580, 4326 }
+    let(:tour_eiffel_in_google) { point 255433.703574246, 6250801.22232508, 900913 }
+
+    it "should return the Point when the target srid is the same" do
+      subject.project_to(subject.srid).should == subject
+    end
+
+    it "should return a Point with the target srid" do
+      subject.project_to(900913).srid.should == 900913
+    end
+
+    it "should return a Point projected into the target srid" do
+      tour_eiffel_in_google.project_to(tour_eiffel_in_wgs84.srid).should be_close_to tour_eiffel_in_wgs84
+      tour_eiffel_in_wgs84.project_to(tour_eiffel_in_google.srid).should be_close_to tour_eiffel_in_google
+    end
+
+  end
+
+  describe ".bounds" do
+
+    let(:some_points) { points("0 0,1 1,1 0,0 1") }
+
+    it "should return an Envelope with the same srid than points" do
+      GeoRuby::SimpleFeatures::Point.bounds(some_points).srid.should == GeoRuby::SimpleFeatures::Point.srid!(some_points)
+    end
+
+    it "should return '0 0,1 1' for '0 0,1 1,1 0,0 1'" do
+      bounds = GeoRuby::SimpleFeatures::Point.bounds(points('0 0,1 1,1 0,0 1'))
+      bounds.lower_corner.should == point(0,0)
+      bounds.upper_corner.should == point(1,1)
+    end
+
+  end
+
+  describe "#endpoint" do
+
+    let(:distance) { 10000 }
+    
+    it "should return a point at the given distance" do
+      subject.endpoint(rand(360), distance).spherical_distance(subject).should be_within(0.001).of(distance)
+    end
+
+    it "should not change longitude when heading is 0 (north)" do
+      subject.endpoint(0, distance).should have_same(:lng).than(subject)
+    end
+
+    it "should not change latitude when heading is 90 (east)" do
+      subject.endpoint(90, distance).lat.should be_within(0.0001).of(subject.lat)
+    end
+
+  end
+
+  describe "#change" do
+    
+    it "should change x if specified" do
+      subject.change(:x => 42).x.should == 42
+    end
+
+    it "should change y if specified" do
+      subject.change(:y => 42).y.should == 42
+    end
+
+    it "should change the srid if specified" do
+      subject.change(:srid => 1).srid.should == 1
+    end
+
+    it "should not change unspecified attributes" do
+      subject.change(:x => 1).should have_same(:y, :srid, :with_z, :with_m).than(subject)
+    end
+
+  end
+
+end

data/spec/lib/georuby/polygon_spec.rb

@@ -0,0 +1,175 @@
+require 'spec_helper'
+
+describe GeoRuby::SimpleFeatures::Polygon do    
+
+  describe "circle" do
+    
+    let(:center) { point 1,1 }
+    let(:radius) { 1000 }
+    let(:sides) { 8 }
+    
+    def circle(arguments = {})
+      arguments = { :center => center, :radius => radius, :sides => sides }.merge(arguments)
+      GeoRuby::SimpleFeatures::Polygon.circle *arguments.values_at(:center, :radius, :sides)
+    end
+
+    it "should create a square" do
+      p1 = center.endpoint 0, radius
+      p2 = center.endpoint 90, radius
+      p3 = center.endpoint 180, radius
+      p4 = center.endpoint 270, radius
+
+      circle(:sides => 4).should == polygon(p1, p2, p3, p4, p1)
+    end
+    
+    it "should have the given side count" do
+      circle(:sides => 16).side_count.should == 16
+    end
+    
+    it "should be closed" do
+      circle(:sides => 16).rings.each do |ring|
+        ring.is_closed.should be_true
+      end
+    end
+    
+    it "should have all its points as the radius distance of the center" do
+      circle.points.all? do |point| 
+        point.spherical_distance(center).should be_within(0.0001).of(radius)
+      end  
+    end
+
+    it "should have the same distance between all its points" do
+      average_distance = circle.perimeter / circle.side_count
+      circle.points.each_cons(2) do |point, next_point|
+        point.spherical_distance(next_point).should be_within(0.01).of(average_distance)
+      end
+    end
+  end
+
+  describe "#project_to" do
+
+    let(:target_srid) { 900913 }
+    subject { polygon("(0 0,1 1,1 0)") }
+
+    it "should project all rings into wgs84" do
+      subject.project_to(target_srid).rings.each_with_index do |ring, index|
+        ring.should == subject[index].project_to(target_srid)
+      end
+    end
+
+    it "should have the target srid" do
+      subject.project_to(target_srid).srid.should == target_srid
+    end
+
+    it "should not change other attributes" do
+      subject.project_to(target_srid).should have_same(:with_z, :with_m).than(subject)
+    end
+
+  end
+
+  describe "to_rgeo" do   
+
+    it "should return a rgeo polygon" do
+      rgeo_polygon =  rgeometry("POLYGON((0 0, 0 5, 5 5, 5 0, 0 0))")    
+      geo_polygon = geometry("POLYGON((0 0, 0 5, 5 5, 5 0, 0 0))") 
+      geo_polygon.to_rgeo.should == rgeo_polygon
+    end
+
+    it "should return a rgeo polygon with an hole" do
+      rgeo_polygon =  rgeometry("POLYGON((0 0, 0 5, 5 5, 5 0, 0 0), (1 1, 2 1, 2 2, 1 2, 1 1))")    
+      geo_polygon = geometry("POLYGON((0 0, 0 5, 5 5, 5 0, 0 0), (1 1, 2 1, 2 2, 1 2, 1 1) )") 
+      geo_polygon.to_rgeo.should == rgeo_polygon
+    end
+
+  end
+  
+  describe "centroid" do
+    let(:georuby_polygon){ geometry("POLYGON((0 0, 0 2, 2 2, 2 0, 0 0))") }      
+    
+    it "should return centroid for a polygon" do
+      georuby_polygon.centroid.should == point(1,1)
+    end
+    
+  end
+
+
+  describe "union" do
+     
+    let(:georuby_polygon) { geometry("POLYGON((0 0, 0 2, 2 2, 2 0, 0 0))") }            
+
+    it "should return the same polygon than in input" do   
+      GeoRuby::SimpleFeatures::Polygon.union([georuby_polygon, georuby_polygon]).should == georuby_polygon
+    end
+    
+    it "should return union of 2 polygons which have a common side" do
+      result =  geometry("POLYGON((0.0 0.0,0.0 2.0,0.0 4.0,2.0 4.0,2.0 2.0,2.0 0.0,0.0 0.0))")
+      georuby_polygon2 = geometry("POLYGON((0 2, 0 4, 2 4, 2 2, 0 2))")
+      GeoRuby::SimpleFeatures::Polygon.union([georuby_polygon, georuby_polygon2]).should == result
+    end
+
+    it "should return 2 polygons which have an area in common" do
+      result =  geometry("POLYGON((0 0,0 1,0 2,0 4,2 4,2 2,2 1,2 0,0 0))")
+      georuby_polygon2 = geometry("POLYGON((0 1, 0 4, 2 4, 2 1, 0 1))")
+      GeoRuby::SimpleFeatures::Polygon.union([georuby_polygon, georuby_polygon2]).should == result
+    end
+
+    it "should return 2 polygons which have no intersection" do
+      result =  geometry("MULTIPOLYGON( ((0 0, 0 2, 2 2, 2 0, 0 0)), ((3 3, 3 4, 4 4, 4 3, 3 3)) )")
+      georuby_polygon2 = geometry("POLYGON((3 3, 4 3, 4 4, 3 4, 3 3))")
+      GeoRuby::SimpleFeatures::Polygon.union([georuby_polygon, georuby_polygon2]).should == result
+    end
+    
+  end
+
+  describe "intersect" do
+    let(:georuby_polygon){ polygon("(0 0,0 2,2 2,2 0,0 0)")  }      
+    let(:georuby_polygon2){ polygon("(0 0,0 4,2 4,2 0,0 0)") }   
+    
+    it "should return intersect of polygons" do
+      test = GeoRuby::SimpleFeatures::Polygon.intersection([georuby_polygon, georuby_polygon2])
+      GeoRuby::SimpleFeatures::Polygon.intersection([georuby_polygon, georuby_polygon2]).text_representation.should == georuby_polygon.text_representation
+    end
+  end
+
+  describe "#difference" do
+    let(:georuby_polygon) {  polygon("(0 0,0 2,2 2,2 0,0 0)")  }      
+    let(:georuby_polygon2){  polygon("(0 0,0 1,2 1,2 0,0 0)") }    
+    
+    it "should be empty geometry collection if polygons are the same" do
+      georuby_polygon.difference(georuby_polygon).should == GeoRuby::SimpleFeatures::GeometryCollection.from_geometries([], 4326)
+    end
+
+    it "should return  polygon('(0 1,0 2,2 2,2 1,0 1)' as result" do
+      result =  polygon("(0 1,0 2,2 2,2 1,0 1)")
+      georuby_polygon.difference(georuby_polygon2).should == result
+    end
+
+    it "should be empty geometry collection if polygon is smaller than the other" do
+      georuby_polygon2.difference(georuby_polygon).should == GeoRuby::SimpleFeatures::GeometryCollection.from_geometries([], 4326)
+    end
+
+    it "should return a result if we make difference with a multi polygon" do
+      geo_polygon = geometry("POLYGON ((0 0, 0 5, 5 5, 5 0, 0 0))")
+      geo_multi_polygon = geometry("MULTIPOLYGON ( ((1 1, 2 1, 2 2, 1 2, 1 1)), ((3 3, 4 3, 4 4, 3 4, 3 3)) )")
+      geo_polygon.difference(geo_multi_polygon).should ==  geometry("POLYGON ( ((0 0, 0 5, 5 5, 5 0, 0 0), (1 1, 2 1, 2 2, 1 2, 1 1), (3 3, 4 3, 4 4, 3 4, 3 3)) )")
+    end
+
+  end
+
+  describe "#==" do
+    
+    it "should be true when points are same" do
+      geometry("POLYGON((0 0,0 1,2 1,2 0,0 0))").should == geometry("POLYGON((0 0,0 1,2 1,2 0,0 0))")
+    end
+
+  end
+
+  describe "side_count" do
+    
+    it "should return 3 for a triangle" do
+      polygon("(0 0,1 0,0 1,0 0)").side_count.should == 3
+    end
+
+  end
+  
+end