ppe-georuby 1.7.2

data/spec/geo_ruby/simple_features/line_string_spec.rb

@@ -0,0 +1,245 @@
+require File.expand_path(File.dirname(__FILE__) + '/../../spec_helper')
+
+module LineSpecHelper
+  def mock_points(num)
+#    @point = mock(Point, :x => 1.0, :y => 2.0)
+    Array.new(num) { |i| mock_point(i,i) }
+  end
+
+  def mock_point(x=1,y=2)
+    mock(Point, :x => x, :y => y, :text_representation => "#{x} #{y}")
+  end
+end
+
+describe LineString do
+  include LineSpecHelper
+
+  before(:each) do
+    @line = LineString.from_points([mock(Point)])
+  end
+
+  it "should instantiate" do
+    violated unless @line
+  end
+
+  it "should have binary_geometry_type 2" do
+    @line.binary_geometry_type.should eql(2)
+  end
+
+  it "should have text_geometry_type" do
+    @line.text_geometry_type.should eql("LINESTRING")
+  end
+
+  it "should have a points array" do
+    @line.points.should be_instance_of(Array)
+  end
+
+  describe "tu converted" do
+
+    it "should concat points" do
+      line_string = LineString::new
+      line_string.concat([Point.from_x_y(12.4,45.3),Point.from_x_y(45.4,41.6)])
+
+      line_string.length.should eql(2)
+      line_string[0].should == Point.from_x_y(12.4,45.3)
+
+      point=Point.from_x_y(123,45.8777)
+      line_string[0]=point
+      line_string[0].should == point
+
+      points=[Point.from_x_y(123,45.8777),Point.from_x_y(45.4,41.6)]
+      line_string.each_index {|i| line_string[i].should == points[i] }
+
+      point=Point.from_x_y(22.4,13.56)
+      line_string << point
+      line_string.length.should eql(3)
+      line_string[2].should eql(point)
+    end
+
+    it "should create" do
+      line_string = LineString.from_points([Point.from_x_y(12.4,-45.3),Point.from_x_y(45.4,41.6)],123)
+      line_string.class.should eql(LineString)
+      line_string.length.should eql(2)
+      line_string[0].should == Point.from_x_y(12.4,-45.3)
+      line_string[1].should == Point.from_x_y(45.4,41.6)
+
+      line_string = LineString.from_coordinates([[12.4,-45.3],[45.4,41.6],[4.456,1.0698]],123)
+      line_string.class.should eql(LineString)
+      line_string.length.should eql(3)
+      line_string[0].should == Point.from_x_y(12.4,-45.3)
+      line_string[1].should == Point.from_x_y(45.4,41.6)
+
+      line_string = LineString.from_coordinates([[12.4,-45.3,123],[45.4,41.6,333],[4.456,1.0698,987]],123,true)
+      line_string.class.should eql(LineString)
+      line_string.length.should eql(3)
+      line_string[0].should == Point.from_x_y_z(12.4,-45.3,123,123)
+
+      line_string = LineString.from_coordinates([[12.4,-45.3,123],[45.4,41.6,333],[4.456,1.0698,987]],123,true)
+      line_string.class.should eql(LineString)
+      line_string.length.should eql(3)
+      line_string[0].should == Point.from_x_y_z(12.4,-45.3,123,123)
+    end
+
+    it "should bbox it" do
+      bbox = LineString.from_coordinates([[12.4,-45.3,123],[45.4,41.6,333],[4.456,1.0698,987]],123,true).bounding_box
+      bbox.length.should eql(2)
+      bbox[0].should == Point.from_x_y_z(4.456,-45.3,123)
+      bbox[1].should == Point.from_x_y_z(45.4,41.6,987)
+    end
+
+    it "should test_line_string_equal" do
+      line_string1 = LineString.from_coordinates([[12.4,-45.3],[45.4,41.6],[4.456,1.0698]],123)
+      line_string2 = LineString.from_coordinates([[12.4,-45.3],[45.4,41.6]],123)
+      point = Point.from_x_y(12.4,-45.3,123)
+
+      line_string1.should == LineString.from_coordinates([[12.4,-45.3],[45.4,41.6],[4.456,1.0698]],123)
+      line_string1.should_not == line_string2
+      line_string1.should_not == point
+    end
+
+    it "should test_line_string_binary" do
+      line_string = LineString.from_coordinates([[12.4,-45.3],[45.4,41.6]],256)
+      line_string.as_hex_ewkb.should eql("01020000200001000002000000CDCCCCCCCCCC28406666666666A646C03333333333B34640CDCCCCCCCCCC4440")
+
+      line_string = LineString.from_coordinates([[12.4,-45.3,35.3],[45.4,41.6,12.3]],256,true)
+      line_string.as_hex_ewkb.should eql("01020000A00001000002000000CDCCCCCCCCCC28406666666666A646C06666666666A641403333333333B34640CDCCCCCCCCCC44409A99999999992840")
+
+      line_string = LineString.from_coordinates([[12.4,-45.3,35.3,45.1],[45.4,41.6,12.3,40.23]],256,true,true)
+      line_string.as_hex_ewkb.should eql("01020000E00001000002000000CDCCCCCCCCCC28406666666666A646C06666666666A64140CDCCCCCCCC8C46403333333333B34640CDCCCCCCCCCC44409A999999999928403D0AD7A3701D4440")
+    end
+
+    it "test_line_string_text" do
+      line_string = LineString.from_coordinates([[12.4,-45.3],[45.4,41.6]],256)
+      line_string.as_ewkt.should eql("SRID=256;LINESTRING(12.4 -45.3,45.4 41.6)")
+
+      line_string = LineString.from_coordinates([[12.4,-45.3,35.3],[45.4,41.6,12.3]],256,true)
+      line_string.as_ewkt.should eql("SRID=256;LINESTRING(12.4 -45.3 35.3,45.4 41.6 12.3)")
+
+      line_string = LineString.from_coordinates([[12.4,-45.3,35.3],[45.4,41.6,12.3]],256,false,true)
+      line_string.as_ewkt.should eql("SRID=256;LINESTRINGM(12.4 -45.3 35.3,45.4 41.6 12.3)")
+
+      line_string = LineString.from_coordinates([[12.4,-45.3,35.3,25.2],[45.4,41.6,12.3,13.75]],256,true,true)
+      line_string.as_ewkt.should eql("SRID=256;LINESTRING(12.4 -45.3 35.3 25.2,45.4 41.6 12.3 13.75)")
+    end
+
+    it "should test_linear_ring_creation" do
+      #testing just the constructor helpers since the rest is the same as for line_string
+      linear_ring = LinearRing.from_coordinates([[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]],345)
+      linear_ring.class.should eql(LinearRing)
+      linear_ring.length.should eql(4)
+      linear_ring.should be_closed
+      linear_ring[1].should == Point.from_x_y(45.4,41.6,345)
+    end
+  end
+
+  describe "> Coordinates" do
+
+    before(:each) do
+      Point.should_receive(:from_coordinates).
+        exactly(4).with(anything, 4326, false, false).and_return(mock_point)
+      @line = LineString.from_coordinates([1.2,2.5,2.2,4.5])
+    end
+
+    it "should instantiate from coordinates" do
+      @line.points.length.should eql(4)
+    end
+
+  end
+
+  describe "> Instantiated" do
+
+    before(:each) do
+      @line = LineString.from_points(mock_points(7))
+    end
+
+    it "should be closed if the last point equals the first" do
+      @line.push(@line.first)
+      @line.should be_closed
+      @line.length.should eql(8)
+    end
+
+    it "should print the text representation" do
+      @line.text_representation.should eql("0 0,1 1,2 2,3 3,4 4,5 5,6 6")
+    end
+
+    it "should print the georss_simple_representation" do
+      @line.georss_simple_representation({:geom_attr => nil}).
+        should eql("<georss:line>0 0 1 1 2 2 3 3 4 4 5 5 6 6</georss:line>\n")
+    end
+
+    it "should map the georss_poslist" do
+      @line.georss_poslist.should eql("0 0 1 1 2 2 3 3 4 4 5 5 6 6")
+    end
+
+    it "should print the kml_representation" do
+      @line.kml_representation.should
+        eql("<LineString>\n<coordinates>0,0 1,1 2,2 3,3 4,4 5,5 6,6</coordinates>\n</LineString>\n")
+    end
+
+    it "should print the kml_poslist without reverse or z" do
+      @line.kml_poslist({}).should eql("0,0 1,1 2,2 3,3 4,4 5,5 6,6")
+    end
+
+    it "should print the kml_poslist reverse" do
+      @line.kml_poslist({:reverse => true}).should eql("6,6 5,5 4,4 3,3 2,2 1,1 0,0")
+    end
+  end
+
+  describe "> Distances..." do
+    before(:each) do
+      @p1 = mock(Point)
+      @p2 = mock(Point)
+      @p3 = mock(Point)
+      @line = LineString.from_points([@p1,@p2,@p3])
+    end
+
+    it "should print the length with haversine" do
+      @p1.should_receive(:spherical_distance).with(@p2).and_return(10)
+      @p2.should_receive(:spherical_distance).with(@p3).and_return(10)
+      @line.spherical_distance.should eql(20)
+    end
+
+    it "should print lenght as euclidian" do
+      @p1.should_receive(:euclidian_distance).with(@p2).and_return(10)
+      @p2.should_receive(:euclidian_distance).with(@p3).and_return(10)
+      @line.euclidian_distance.should eql(20)
+    end
+  end
+
+  describe "Simplify" do
+
+    before do
+      @line = LineString.from_coordinates([[6,0],[4,1],[3,4],[4,6],[5,8],[5,9],[4,10],[6,15]], 4326)
+    end
+
+    it "should simplify a simple linestring" do
+      line =  LineString.from_coordinates([[12,15],[14,17],[17, 20]], 4326)
+      line.simplify.should have(2).points
+    end
+
+    it "should simplify a harder linestring" do
+      @line.simplify(6).should have(6).points
+      @line.simplify(9).should have(4).points
+      @line.simplify(10).should have(3).points
+    end
+
+    it "should barely touch it" do
+      @line.simplify(1).should have(7).points
+    end
+
+    it "should simplify to five points" do
+      @line.simplify(7).should have(5).points
+    end
+
+    it "should flatten it" do
+      @line.simplify(11).should have(2).points
+    end
+
+    it "should be the first and last in a flatten" do
+      line = @line.simplify(11)
+      line[0].should be_a_point(6, 0)
+      line[1].should be_a_point(6, 15)
+    end
+
+  end
+end

data/spec/geo_ruby/simple_features/linear_ring_spec.rb

@@ -0,0 +1,14 @@
+require File.expand_path(File.dirname(__FILE__) + '/../../spec_helper')
+
+describe LinearRing do
+#  before(:each) do
+#    @po = Polygon.new
+#  end
+
+  it "should instantiate" do
+    lr = LinearRing.new([1.2,2.5,2.2,4.5])
+    lr.should be_instance_of(LinearRing)
+  end
+
+
+end

data/spec/geo_ruby/simple_features/multi_line_string_spec.rb

@@ -0,0 +1,54 @@
+require File.expand_path(File.dirname(__FILE__) + '/../../spec_helper')
+
+describe MultiLineString do
+
+ it "test_multi_line_string_creation" do
+    multi_line_string1 = MultiLineString.from_line_strings([LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012]],256),LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012],[45.123,123.3]],256)],256)
+    multi_line_string1.should be_instance_of(MultiLineString)
+    multi_line_string1.length.should eql(2)
+    multi_line_string1[0].should == LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012]],256)
+
+    multi_line_string2= MultiLineString.from_coordinates([[[1.5,45.2],[-54.12312,-0.012]],[[1.5,45.2],[-54.12312,-0.012],[45.123,123.3]]],256);
+    multi_line_string1.should be_instance_of(MultiLineString)
+    multi_line_string1.length.should eql(2)
+    multi_line_string2[0].should == LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012]],256)
+
+    multi_line_string2.should == multi_line_string2
+  end
+
+  it "test_multi_line_string_binary" do
+    multi_line_string = MultiLineString.from_line_strings([LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012]],256),LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012],[45.123,123.3]],256)],256)
+    multi_line_string.as_hex_ewkb.should eql("01050000200001000002000000010200000002000000000000000000F83F9A99999999994640E4BD6A65C20F4BC0FA7E6ABC749388BF010200000003000000000000000000F83F9A99999999994640E4BD6A65C20F4BC0FA7E6ABC749388BF39B4C876BE8F46403333333333D35E40")
+
+    multi_line_string = MultiLineString.from_line_strings([LineString.from_coordinates([[1.5,45.2,1.3,1.2],[-54.12312,-0.012,1.2,4.5]],256,true,true),LineString.from_coordinates([[1.5,45.2,5.1,-4.5],[-54.12312,-0.012,-6.8,3.4],[45.123,123.3,4.5,-5.3]],256,true,true)],256,true,true)
+    multi_line_string.as_hex_ewkb.should eql("0105000020000100000200000001020000C002000000000000000000F83F9A99999999994640CDCCCCCCCCCCF43F333333333333F33FE4BD6A65C20F4BC0FA7E6ABC749388BF333333333333F33F000000000000124001020000C003000000000000000000F83F9A99999999994640666666666666144000000000000012C0E4BD6A65C20F4BC0FA7E6ABC749388BF3333333333331BC03333333333330B4039B4C876BE8F46403333333333D35E40000000000000124033333333333315C0")
+  end
+
+  it "test_multi_line_string_text" do
+    multi_line_string = MultiLineString.from_line_strings([LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012]],256),LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012],[45.123,123.3]],256)],256)
+    multi_line_string.as_ewkt.should eql("SRID=256;MULTILINESTRING((1.5 45.2,-54.12312 -0.012),(1.5 45.2,-54.12312 -0.012,45.123 123.3))")
+
+    multi_line_string = MultiLineString.from_line_strings([LineString.from_coordinates([[1.5,45.2,1.3,1.2],[-54.12312,-0.012,1.2,4.5]],256,true,true),LineString.from_coordinates([[1.5,45.2,5.1,-4.5],[-54.12312,-0.012,-6.8,3.4],[45.123,123.3,4.5,-5.3]],256,true,true)],256,true,true)
+    multi_line_string.as_ewkt.should eql("SRID=256;MULTILINESTRING((1.5 45.2 1.3 1.2,-54.12312 -0.012 1.2 4.5),(1.5 45.2 5.1 -4.5,-54.12312 -0.012 -6.8 3.4,45.123 123.3 4.5 -5.3))")
+  end
+
+  describe "Some More" do
+    before do
+      @mls = MultiLineString.from_line_strings([LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012]],256),LineString.from_coordinates([[1.5,45.2],[-54.12312,-0.012],[45.123,123.3]],256)],256)
+    end
+
+    it "should have a accessor to all points" do
+      @mls.points.should be_instance_of(Array)
+    end
+
+    it "should flatten the array" do
+      @mls.should have(5).points
+    end
+
+    it "should simplify to linestring" do
+      ls = @mls.to_line_string
+      ls.should be_instance_of(LineString)
+      ls.should have(5).points
+    end
+  end
+end

data/spec/geo_ruby/simple_features/multi_point_spec.rb

@@ -0,0 +1,35 @@
+require File.expand_path(File.dirname(__FILE__) + '/../../spec_helper')
+
+describe MultiPoint do
+
+  it "test_multi_point_creation" do
+    multi_point = MultiPoint.from_coordinates([[12.4,-123.3],[-65.1,123.4],[123.55555555,123]],444)
+    multi_point.should be_instance_of(MultiPoint)
+    multi_point.length.should eql(3)
+    multi_point[0].should == Point.from_x_y(12.4,-123.3,444)
+    multi_point[2].should == Point.from_x_y(123.55555555,123,444)
+  end
+
+  it "test_multi_point_binary" do
+    multi_point = MultiPoint.from_coordinates([[12.4,-123.3],[-65.1,123.4],[123.55555555,123]],444)
+    multi_point.as_hex_ewkb.should eql("0104000020BC010000030000000101000000CDCCCCCCCCCC28403333333333D35EC0010100000066666666664650C09A99999999D95E4001010000001F97DD388EE35E400000000000C05E40")
+
+    multi_point = MultiPoint.from_coordinates([[12.4,-123.3,4.5],[-65.1,123.4,1.2],[123.55555555,123,2.3]],444,true)
+    multi_point.as_hex_ewkb.should eql("01040000A0BC010000030000000101000080CDCCCCCCCCCC28403333333333D35EC00000000000001240010100008066666666664650C09A99999999D95E40333333333333F33F01010000801F97DD388EE35E400000000000C05E406666666666660240")
+  end
+
+  it "test_multi_point_text" do
+    multi_point = MultiPoint.from_coordinates([[12.4,-123.3],[-65.1,123.4],[123.55555555,123]],444)
+    multi_point.as_ewkt.should eql("SRID=444;MULTIPOINT((12.4 -123.3),(-65.1 123.4),(123.55555555 123))")
+
+    multi_point = MultiPoint.from_coordinates([[12.4,-123.3,4.5],[-65.1,123.4,6.7],[123.55555555,123,7.8]],444,true)
+    multi_point.as_ewkt.should eql("SRID=444;MULTIPOINT((12.4 -123.3 4.5),(-65.1 123.4 6.7),(123.55555555 123 7.8))")
+  end
+
+  it "should respond to points" do
+    mp = MultiPoint.from_coordinates([[12.4,-123.3],[-65.1,123.4],[123.55555555,123]],444)
+    mp.should have(3).geometries
+    mp.should have(3).points
+  end
+
+end

data/spec/geo_ruby/simple_features/multi_polygon_spec.rb

@@ -0,0 +1,50 @@
+require File.expand_path(File.dirname(__FILE__) + '/../../spec_helper')
+
+describe MultiPolygon do
+
+   it "test_multi_polygon_creation" do
+    multi_polygon1 = MultiPolygon.from_polygons([Polygon.from_coordinates([[[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]],[[2.4,5.3],[5.4,1.4263],[14.46,1.06],[2.4,5.3]]],256),Polygon.from_coordinates([[[0,0],[4,0],[4,4],[0,4],[0,0]],[[1,1],[3,1],[3,3],[1,3],[1,1]]],256)],256)
+    multi_polygon1.should be_instance_of(MultiPolygon)
+    multi_polygon1.length.should eql(2)
+    multi_polygon1[0].should == Polygon.from_coordinates([[[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]],[[2.4,5.3],[5.4,1.4263],[14.46,1.06],[2.4,5.3]]],256)
+
+    multi_polygon2 = MultiPolygon.from_coordinates([[[[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]],[[2.4,5.3],[5.4,1.4263],[14.46,1.06],[2.4,5.3]]],[[[0,0],[4,0],[4,4],[0,4],[0,0]],[[1,1],[3,1],[3,3],[1,3],[1,1]]]],256)
+    multi_polygon1.should be_instance_of(MultiPolygon)
+    multi_polygon1.length.should eql(2)
+
+    multi_polygon2[0].should == Polygon.from_coordinates([[[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]],[[2.4,5.3],[5.4,1.4263],[14.46,1.06],[2.4,5.3]]],256)
+    multi_polygon1.should == multi_polygon2
+  end
+
+  it "test_multi_polygon_binary" do
+    multi_polygon = MultiPolygon.from_polygons([Polygon.from_coordinates([[[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]],[[2.4,5.3],[5.4,1.4263],[14.46,1.06],[2.4,5.3]]],256),Polygon.from_coordinates([[[0,0],[4,0],[4,4],[0,4],[0,0]],[[1,1],[3,1],[3,3],[1,3],[1,1]]],256)],256)
+    multi_polygon.as_hex_ewkb.should eql("0106000020000100000200000001030000000200000004000000CDCCCCCCCCCC28406666666666A646C03333333333B34640CDCCCCCCCCCC44406DE7FBA9F1D211403D2CD49AE61DF13FCDCCCCCCCCCC28406666666666A646C004000000333333333333034033333333333315409A999999999915408A8EE4F21FD2F63FEC51B81E85EB2C40F6285C8FC2F5F03F3333333333330340333333333333154001030000000200000005000000000000000000000000000000000000000000000000001040000000000000000000000000000010400000000000001040000000000000000000000000000010400000000000000000000000000000000005000000000000000000F03F000000000000F03F0000000000000840000000000000F03F00000000000008400000000000000840000000000000F03F0000000000000840000000000000F03F000000000000F03F")
+
+    multi_polygon = MultiPolygon.from_polygons([Polygon.from_coordinates([[[12.4,-45.3,1.2],[45.4,41.6,1.2],[4.456,1.0698,1.2],[12.4,-45.3,1.2]],[[2.4,5.3,1.2],[5.4,1.4263,1.2],[14.46,1.06,1.2],[2.4,5.3,1.2]]],256,false,true),Polygon.from_coordinates([[[0,0,1.2],[4,0,1.2],[4,4,2.3],[0,4,1.2],[0,0,1.2]],[[1,1,2.2],[3,1,3.3],[3,3,1.1],[1,3,2.4],[1,1,2.2]]],256,false,true)],256,false,true)
+    multi_polygon.as_hex_ewkb.should eql("0106000020000100000200000001030000400200000004000000CDCCCCCCCCCC28406666666666A646C0333333333333F33F3333333333B34640CDCCCCCCCCCC4440333333333333F33F6DE7FBA9F1D211403D2CD49AE61DF13F333333333333F33FCDCCCCCCCCCC28406666666666A646C0333333333333F33F0400000033333333333303403333333333331540333333333333F33F9A999999999915408A8EE4F21FD2F63F333333333333F33FEC51B81E85EB2C40F6285C8FC2F5F03F333333333333F33F33333333333303403333333333331540333333333333F33F0103000040020000000500000000000000000000000000000000000000333333333333F33F00000000000010400000000000000000333333333333F33F00000000000010400000000000001040666666666666024000000000000000000000000000001040333333333333F33F00000000000000000000000000000000333333333333F33F05000000000000000000F03F000000000000F03F9A999999999901400000000000000840000000000000F03F6666666666660A40000000000000084000000000000008409A9999999999F13F000000000000F03F00000000000008403333333333330340000000000000F03F000000000000F03F9A99999999990140")
+  end
+
+  it "test_multi_polygon_text" do
+    multi_polygon = MultiPolygon.from_polygons([Polygon.from_coordinates([[[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]],[[2.4,5.3],[5.4,1.4263],[14.46,1.06],[2.4,5.3]]],256),Polygon.from_coordinates([[[0,0],[4,0],[4,4],[0,4],[0,0]],[[1,1],[3,1],[3,3],[1,3],[1,1]]],256)],256)
+    multi_polygon.as_ewkt.should eql("SRID=256;MULTIPOLYGON(((12.4 -45.3,45.4 41.6,4.456 1.0698,12.4 -45.3),(2.4 5.3,5.4 1.4263,14.46 1.06,2.4 5.3)),((0 0,4 0,4 4,0 4,0 0),(1 1,3 1,3 3,1 3,1 1)))")
+
+    multi_polygon = MultiPolygon.from_polygons([Polygon.from_coordinates([[[12.4,-45.3,2],[45.4,41.6,3],[4.456,1.0698,4],[12.4,-45.3,2]],[[2.4,5.3,1],[5.4,1.4263,3.44],[14.46,1.06,4.5],[2.4,5.3,1]]],4326,true),Polygon.from_coordinates([[[0,0,5.6],[4,0,5.4],[4,4,1],[0,4,23],[0,0,5.6]],[[1,1,2.3],[3,1,4],[3,3,5],[1,3,6],[1,1,2.3]]],4326,true)],4326,true)
+    multi_polygon.as_ewkt.should eql("SRID=4326;MULTIPOLYGON(((12.4 -45.3 2,45.4 41.6 3,4.456 1.0698 4,12.4 -45.3 2),(2.4 5.3 1,5.4 1.4263 3.44,14.46 1.06 4.5,2.4 5.3 1)),((0 0 5.6,4 0 5.4,4 4 1,0 4 23,0 0 5.6),(1 1 2.3,3 1 4,3 3 5,1 3 6,1 1 2.3)))")
+  end
+
+  describe "Counting" do
+
+    before do
+      @mp = MultiPolygon.from_polygons([Polygon.from_coordinates([[[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]],[[2.4,5.3],[5.4,1.4263],[14.46,1.06],[2.4,5.3]]],256),Polygon.from_coordinates([[[0,0],[4,0],[4,4],[0,4],[0,0]],[[1,1],[3,1],[3,3],[1,3],[1,1]]],256)],256)
+    end
+
+    it "should have a points method" do
+      @mp.points[0].should be_instance_of(Point)
+    end
+
+    it "should flatten it right" do
+      @mp.should have(18).points
+    end
+
+  end
+end

data/spec/geo_ruby/simple_features/point_spec.rb

@@ -0,0 +1,356 @@
+# -*- coding: utf-8 -*-
+require File.expand_path(File.dirname(__FILE__) + '/../../spec_helper')
+
+describe Point do
+  before(:each) do
+    @point = Point.new(4326)
+  end
+
+  it "should instantiatember" do
+    violated unless @point
+    @point.should be_instance_of(Point)
+  end
+
+  it "should have a nice matcher" do
+    @point.should be_a_point
+  end
+
+  it "should have a very nice matcher" do
+    @point.should be_a_point(0.0, 0.0)
+  end
+
+  it "should have a very nice matcher" do
+    Point.from_x_y_z_m(1,2,3.33,"t").should be_a_point(1, 2, 3.33, "t")
+  end
+
+  it "should have a dumb matcher" do
+    Point.should be_geometric
+  end
+
+  it "should be subclassable" do
+    place = Class.new(Point)
+    p = place.from_x_y(1,2)
+    p.should be_a place
+  end
+
+  it "should have binary_geometry_type 2" do
+    @point.binary_geometry_type.should eql(1)
+  end
+
+  it "should have the correct srid" do
+    @point.srid.should eql(4326)
+  end
+
+  it "should not have z or m" do
+    @point.with_z.should be_false
+    @point.should_not be_with_z
+    @point.with_m.should be_false
+    @point.should_not be_with_m
+  end
+
+  it "should set params to 0.0" do
+    @point.x.should eql(0.0)
+    @point.y.should eql(0.0)
+    @point.z.should eql(0.0)
+    @point.m.should eql(0.0)
+  end
+
+  it "should compare ok" do
+    point1= Point::new
+    point1.set_x_y(1.5,45.4)
+    point2= Point::new
+    point2.set_x_y(1.5,45.4)
+    point3= Point::new
+    point3.set_x_y(4.5,12.3)
+    point4= Point::new
+    point4.set_x_y_z(1.5,45.4,423)
+    point5= Point::new
+    point5.set_x_y(1.5,45.4)
+    point5.m=15
+    geometry= Geometry::new
+
+    point1.should == point2
+    point1.should_not == point3
+    point1.should_not == point4
+    point1.should_not == point5
+    point1.should_not == geometry
+  end
+
+  describe "> Instantiation" do
+
+    it "should instantiate a 2d point" do
+      point = Point.from_x_y(10,20,123)
+      point.x.should eql(10)
+      point.y.should eql(20)
+      point.srid.should eql(123)
+      point.z.should eql(0.0)
+    end
+
+    it "should instantiate a 3d point" do
+      point = Point.from_x_y_z(-10,-20,-30)
+      point.x.should eql(-10)
+      point.y.should eql(-20)
+      point.z.should eql(-30)
+    end
+
+    it "should instantiate a 3d(m) point" do
+      point = Point.from_x_y_m(10,20,30)
+      point.x.should eql(10)
+      point.y.should eql(20)
+      point.m.should eql(30)
+      point.z.should eql(0.0)
+    end
+
+    it "should instantiate a 4d point" do
+      point = Point.from_x_y_z_m(10,20,30,40,123)
+      point.x.should eql(10)
+      point.y.should eql(20)
+      point.z.should eql(30)
+      point.m.should eql(40)
+      point.srid.should eql(123)
+    end
+
+    it "should instantiate a point from polar coordinates" do
+      point = Point.from_r_t(1.4142,45)
+      point.y.should be_close(1, 0.00001)
+      point.x.should be_close(1, 0.00001)
+    end
+
+    it "should instantiate from coordinates x,y" do
+      point = Point.from_coordinates([1.6,2.8],123)
+      point.x.should eql(1.6)
+      point.y.should eql(2.8)
+      point.should_not be_with_z
+      point.z.should eql(0.0)
+      point.srid.should eql(123)
+    end
+
+    it "should instantiate from coordinates x,y,z" do
+      point = Point.from_coordinates([1.6,2.8,3.4],123, true)
+      point.x.should eql(1.6)
+      point.y.should eql(2.8)
+      point.z.should eql(3.4)
+      point.should be_with_z
+      point.srid.should eql(123)
+    end
+
+    it "should instantiate from coordinates x,y,z,m" do
+      point = Point.from_coordinates([1.6,2.8,3.4,15],123, true, true)
+      point.x.should eql(1.6)
+      point.y.should eql(2.8)
+      point.z.should eql(3.4)
+      point.m.should eql(15)
+      point.should be_with_z
+      point.should be_with_m
+      point.srid.should eql(123)
+    end
+
+    it "should have a bbox" do
+      bbox = Point.from_x_y_z_m(-1.6,2.8,-3.4,15,123).bounding_box
+      bbox.length.should eql(2)
+      bbox[0].should == Point.from_x_y_z(-1.6,2.8,-3.4)
+      bbox[1].should == Point.from_x_y_z(-1.6,2.8,-3.4)
+    end
+
+    it "should parse lat long" do
+      Point.from_latlong("-20° 47' 26.37","-20° 47' 26.37").x.should be_close(-20.790658, 0.00001)
+      Point.from_latlong("20° 47' 26.378","20° 47' 26.378").y.should be_close(20.790658, 0.00001)
+    end
+
+    it "should parse lat long w/o sec" do
+      Point.from_latlong("-20°47′26″","-20°47′26″").x.should be_close(-20.790555, 0.00001)
+      Point.from_latlong("20°47′26″","20°47′26″").y.should be_close(20.790555, 0.00001)
+    end
+
+    it "should accept with W or S notation" do
+      Point.from_latlong("20° 47' 26.37 W","20° 47' 26.37 S").x.should be_close(-20.790658, 0.00001)
+      Point.from_latlong("20° 47' 26.37 W","20° 47' 26.37 S").y.should be_close(-20.790658, 0.00001)
+    end
+
+    it "should instantiate a point from positive degrees" do
+      point = Point.from_latlong('47`20 06.09E','22`50 77.35N')
+      point.y.should be_close(22.8548194, 0.000001)
+      point.x.should be_close(47.335025, 0.000001)
+    end
+
+    it "should instantiate a point from negative degrees" do
+      point = Point.from_latlong('47`20 06.09W','22`50 77.35S')
+      point.y.should be_close(-22.8548194, 0.000001)
+      point.x.should be_close(-47.335025, 0.000001)
+    end
+
+    it "should print out nicely" do
+      Point.from_x_y(47.88, -20.1).as_latlong.should eql("47°52′48″, -20°06′00″")
+    end
+
+    it "should print out nicely" do
+      Point.from_x_y(-20.78, 20.78).as_latlong(:full => true).should eql("-20°46′48.00″, 20°46′48.00″")
+    end
+
+    it "should print out nicely" do
+      Point.from_x_y(47.11, -20.2).as_latlong(:full => true).should eql("47°06′36.00″, -20°11′60.00″")
+    end
+
+    it "should print out nicely" do
+      Point.from_x_y(47.11, -20.2).as_latlong(:coord => true).should eql("47°06′36″N, 20°11′60″W")
+    end
+
+    it "should print out nicely" do
+      Point.from_x_y(-47.11, 20.2).as_latlong(:full => true,:coord => true).should eql("47°06′36.00″S, 20°11′60.00″E")
+    end
+
+  end
+
+  describe " > Distance & Bearing" do
+
+    before(:each) do
+      @p1 = Point.from_x_y(1,1)
+      @p2 = Point.from_x_y(2,2)
+    end
+
+    it "and a 3th grade child should calculate euclidian distance" do
+      @p1.euclidian_distance(@p2).
+        should be_close(1.4142135623731, 0.00000001)
+    end
+
+    it "should calculate spherical distance" do
+      @p1.spherical_distance(@p2).
+        should be_close(157225.358003181,0.00000001)
+    end
+
+    it "should calculate ellipsoidal distance" do
+      @p1.ellipsoidal_distance(@p2).
+        should be_close(156876.149400742, 0.00000001)
+    end
+
+    describe "Orthogonal Distance" do
+      before do
+        @line = LineString.from_coordinates([[0,0],[1,3]], 4326)
+        @line2 = LineString.from_coordinates([[1,1],[1,2]], 4326)
+      end
+
+      it "should calcula orthogonal distance from a line (90 deg)" do
+        @p1.orthogonal_distance(@line).should be_close(1.414, 0.001)
+      end
+
+      it "should calcula orthogonal distance very close..." do
+        @p1.orthogonal_distance(@line2).should be_zero
+      end
+
+      it "should calcula orthogonal distance from a line (90 deg)" do
+        @p2.orthogonal_distance(@line).should be_close(2.828, 0.001)
+      end
+
+      it "should calcula orthogonal distance from a line (0 deg)" do
+        @p2.orthogonal_distance(@line2).should be_close(1.0, 0.1)
+      end
+
+      it "should calcula orthogonal distance from a line (0 deg)" do
+        @p2.orthogonal_distance(@line2).should be_close(1.0, 0.1)
+      end
+
+    end
+
+    it "should calculate the bearing from apoint to another in degrees" do
+      @p1.bearing_to(@p2).should be_close(45.0, 0.01)
+    end
+
+    it "should calculate the bearing from apoint to another in degrees" do
+      p3 = Point.from_x_y(1,-1)
+      @p1.bearing_to(p3).should be_close(180.0, 0.01)
+    end
+
+    it "should calculate the bearing from apoint to another in degrees" do
+      p3 = Point.from_x_y(-1,-1)
+      @p1.bearing_to(p3).should be_close(225.0, 0.01)
+    end
+
+    it "should calculate the bearing from apoint to another in degrees" do
+      p3 = Point.from_x_y(-1,1)
+      @p1.bearing_to(p3).should be_close(270.0, 0.01)
+    end
+
+    it "should calculate the bearing from apoint to another in degrees" do
+      p3 = Point.from_x_y(2,-1)
+      @p1.bearing_to(p3).should be_close(153.4349488, 0.0001)
+    end
+
+    it "should calculate a clone point bearing to 0" do
+      @p1.bearing_to(@p1).should eql(0)
+    end
+
+    it "should calculate the bearing from apoint to another in degrees" do
+      @p1.bearing_text(@p2).should eql(:ne)
+    end
+
+    it "should calculate the bearing from apoint to another in degrees" do
+      p3 = Point.from_x_y(-1,1)
+      @p1.bearing_text(p3).should eql(:w)
+    end
+
+  end
+
+  describe "> Export Formats" do
+
+    before(:each) do
+      @point = Point.from_x_y( -11.2431, 32.3141 )
+    end
+
+    it "should print nicely" do
+      @point.text_representation.should eql("-11.2431 32.3141")
+    end
+
+    it "should printoout as binary" do
+      Point.from_x_y(12.4,45.3,123).as_hex_ewkb.should eql("01010000207B000000CDCCCCCCCCCC28406666666666A64640")
+      point = Point.from_x_y_z_m(12.4,45.3,-3.5,15,123)
+      point.as_hex_ewkb.should eql("01010000E07B000000CDCCCCCCCCCC28406666666666A646400000000000000CC00000000000002E40")
+      point.as_hex_wkb.should eql("0101000000CDCCCCCCCCCC28406666666666A64640")
+    end
+
+    it "should printoout as text" do
+      Point.from_x_y(12.4,45.3,123).as_ewkt.should eql("SRID=123;POINT(12.4 45.3)")
+      point = Point.from_x_y_z(12.4,45.3,-3.5,123)
+      point.as_ewkt.should eql("SRID=123;POINT(12.4 45.3 -3.5)")
+      point.as_wkt.should eql("POINT(12.4 45.3)")
+      point.as_ewkt(false,true).should eql("POINT(12.4 45.3 -3.5)")
+      point = Point.from_x_y_m(12.4,45.3,-3.5,123)
+      point.as_ewkt.should eql("SRID=123;POINTM(12.4 45.3 -3.5)")
+      point.as_ewkt(true,true,false).should eql("SRID=123;POINT(12.4 45.3)")
+    end
+
+    it "should have a nice bounding box" do
+      @point.should have(2).bounding_box
+      @point.bounding_box.each do |point|
+        point.x.should eql(@point.x)
+        point.y.should eql(@point.y)
+      end
+    end
+
+    it "should print as kml too" do
+      @point.kml_representation.should eql("<Point>\n<coordinates>-11.2431,32.3141</coordinates>\n</Point>\n")
+    end
+
+    it "should print as georss" do
+      @point.georss_simple_representation(:georss_ns => 'hey').should eql("<hey:point>32.3141 -11.2431</hey:point>\n")
+    end
+
+    it "should print r (polar coords)" do
+      @point.r.should be_close(34.214154, 0.00001)
+    end
+
+    it "should print theta as degrees" do
+      @point.theta_deg.should be_close(289.184406352127, 0.0001)
+    end
+
+    it "should print theta as radians" do
+      @point.theta_rad.should be_close(5.04722003626982, 0.0001)
+    end
+
+    it "should output as polar" do
+      @point.as_polar.should be_instance_of(Array)
+      @point.should have(2).as_polar #.length.should eql(2)
+    end
+
+  end
+
+end