georuby 1.9.3

data/spec/geo_ruby/geojson_spec.rb

@@ -0,0 +1,147 @@
+require File.expand_path(File.dirname(__FILE__) + '/../spec_helper')
+
+DATA_DIR = File.dirname(__FILE__) + '/../data/geojson/'
+
+# All geojson test examples are from the GeoJSON spec unless otherwise
+# specified
+#
+# TODO Refactor comon test approaches into methods
+# TODO Add use of contexts?
+describe GeojsonParser do
+
+  it "should create a specified Point" do
+    point_json = %{ { "type": "Point", "coordinates": [100.0, 0.0] } }
+    point = Geometry.from_geojson(point_json)
+    point.class.should eql(Point)
+    point_hash = JSON.parse(point_json)
+    point.to_coordinates.should eql(point_hash['coordinates'])
+  end
+
+  it "should create a specified LineString" do
+    ls_json = %{ { "type": "LineString", "coordinates": [ [100.0, 0.0], [101.0, 1.0] ]} }
+    line_string = Geometry.from_geojson(ls_json)
+    line_string.class.should eql(LineString)
+    ls_hash = JSON.parse(ls_json)
+    line_string.to_coordinates.should eql(ls_hash['coordinates'])
+  end
+
+  it "should create a specified Polygon" do
+    poly_json = <<-EOJ
+      { "type": "Polygon",
+        "coordinates": [
+          [ [100.0, 0.0], [101.0, 0.0], [101.0, 1.0], [100.0, 1.0], [100.0, 0.0] ],
+          [ [100.2, 0.2], [100.8, 0.2], [100.8, 0.8], [100.2, 0.8], [100.2, 0.2] ]
+          ]
+      }
+    EOJ
+    polygon = Geometry.from_geojson(poly_json)
+    polygon.class.should eql(Polygon)
+    polygon.rings.size.should eql(2)
+    poly_hash = JSON.parse(poly_json)
+    polygon.to_coordinates.should eql(poly_hash['coordinates'])
+  end
+
+  it "should create a specified MultiPoint" do
+    mp_json = <<-EOJ
+      { "type": "MultiPoint",
+        "coordinates": [ [100.0, 0.0], [101.0, 1.0] ]
+        }
+    EOJ
+    multi_point = Geometry.from_geojson(mp_json)
+    multi_point.class.should eql(MultiPoint)
+    mp_hash = JSON.parse(mp_json)
+    multi_point.to_coordinates.should eql(mp_hash['coordinates'])
+  end
+
+  it "should create a specified MultiLineString" do
+    mls_json = <<-EOJ
+      { "type": "MultiLineString",
+        "coordinates": [
+            [ [100.0, 0.0], [101.0, 1.0] ],
+            [ [102.0, 2.0], [103.0, 3.0] ]
+          ]
+      }
+    EOJ
+    multi_ls = Geometry.from_geojson(mls_json)
+    multi_ls.class.should eql(MultiLineString)
+    mls_hash = JSON.parse(mls_json)
+    multi_ls.to_coordinates.should eql(mls_hash['coordinates'])
+  end
+
+  it "should create a specifiead MultiPolygon" do
+    mpoly_json = <<-EOJ
+      { "type": "MultiPolygon",
+        "coordinates": [
+          [[[102.0, 2.0], [103.0, 2.0], [103.0, 3.0], [102.0, 3.0], [102.0, 2.0]]],
+          [[[100.0, 0.0], [101.0, 0.0], [101.0, 1.0], [100.0, 1.0], [100.0, 0.0]],
+           [[100.2, 0.2], [100.8, 0.2], [100.8, 0.8], [100.2, 0.8], [100.2, 0.2]]]
+          ]
+        }
+    EOJ
+    mpoly = Geometry.from_geojson(mpoly_json)
+    mpoly.class.should eql(MultiPolygon)
+    mpoly_hash = JSON.parse(mpoly_json)
+    mpoly.to_coordinates.should eql(mpoly_hash['coordinates'])
+  end
+
+  it "should create a specified GeometryCollection" do
+    gcol_json = <<-EOJ
+      { "type": "GeometryCollection",
+        "geometries": [
+          { "type": "Point",
+            "coordinates": [100.0, 0.0]
+            },
+          { "type": "LineString",
+            "coordinates": [ [101.0, 0.0], [102.0, 1.0] ]
+            }
+        ]
+      }
+    EOJ
+    gcol = Geometry.from_geojson(gcol_json)
+    gcol.class.should eql(GeometryCollection)
+    gcol_hash = JSON.parse(gcol_json)
+    gcol.geometries.each_with_index do |g,i|
+      gh = gcol_hash['geometries'][i]
+      g.class.should eql(GeoRuby::SimpleFeatures.const_get(gh['type']))
+      g.to_coordinates.should eql(gh['coordinates'])
+    end
+  end
+
+  # Feature GeoJSON test example from wikipedia entry
+  it "should create a specified Feature" do
+    feature_json = <<-EOJ
+      {
+        "type":"Feature",
+        "id":"OpenLayers.Feature.Vector_314",
+        "properties":{"prop0": "value0"},
+        "geometry":{
+          "type":"Point",
+          "coordinates":[97.03125, 39.7265625]
+        }
+      }
+    EOJ
+    f = Geometry.from_geojson(feature_json)
+    f.class.should eql(GeojsonFeature)
+    feature_hash = JSON.parse(feature_json)
+    f.id.should eql(feature_hash['id'])
+    f.properties.should eql(feature_hash['properties'])
+    f.geometry.class.should eql(GeoRuby::SimpleFeatures.const_get(feature_hash['geometry']['type']))
+    f.geometry.to_coordinates.should eql(feature_hash['geometry']['coordinates'])
+  end
+
+  it "should create a specified FeatureCollection" do
+    fcol_json = File.read(DATA_DIR + 'feature_collection.json')
+    fcol = Geometry.from_geojson(fcol_json)
+    fcol.class.should eql(GeojsonFeatureCollection)
+    fcol_hash = JSON.parse(fcol_json)
+    fcol.features.each_with_index do |f,i|
+      fgh = fcol_hash['features'][i]['geometry']
+      fg = f.geometry
+      f.properties.should eql(fcol_hash['features'][i]['properties'])
+      fg.class.should eql(GeoRuby::SimpleFeatures.const_get(fgh['type']))
+      fg.to_coordinates.should eql(fgh['coordinates'])
+    end
+  end
+
+end
+

data/spec/geo_ruby/georss.rb

@@ -0,0 +1,218 @@
+require File.expand_path(File.dirname(__FILE__) + '/../../spec_helper')
+
+describe GeorssParser do
+  it "test_point_creation" do
+    point = Point.from_x_y(3,4)
+
+    point.as_georss(:dialect => :simple, :elev => 45.7, :featuretypetag => "hoyoyo").gsub("\n","").should eql("<georss:point featuretypetag=\"hoyoyo\" elev=\"45.7\">4 3</georss:point>")
+    point.as_georss(:dialect => :w3cgeo).gsub("\n","").should eql("<geo:lat>4</geo:lat><geo:long>3</geo:long>")
+    point.as_georss(:dialect => :gml).gsub("\n","").should eql("<georss:where><gml:Point><gml:pos>4 3</gml:pos></gml:Point></georss:where>")
+
+    point.as_kml(:id => "HOYOYO-42").gsub("\n","").should eql("<Point id=\"HOYOYO-42\"><coordinates>3,4</coordinates></Point>")
+  end
+
+  it "test_line_string" do
+    ls = LineString.from_points([Point.from_lon_lat_z(12.4,-45.3,56),Point.from_lon_lat_z(45.4,41.6,45)],123,true)
+
+    ls.as_georss.gsub("\n","").should eql("<georss:line>-45.3 12.4 41.6 45.4</georss:line>")
+    ls.as_georss(:dialect => :w3cgeo).gsub("\n","").should eql("<geo:lat>-45.3</geo:lat><geo:long>12.4</geo:long>")
+    ls.as_georss(:dialect => :gml).gsub("\n","").should eql("<georss:where><gml:LineString><gml:posList>-45.3 12.4 41.6 45.4</gml:posList></gml:LineString></georss:where>")
+
+    ls.as_kml(:extrude => 1, :altitude_mode => "absolute").gsub("\n","").should eql("<LineString><extrude>1</extrude><altitudeMode>absolute</altitudeMode><coordinates>12.4,-45.3,56 45.4,41.6,45</coordinates></LineString>")
+  end
+
+  it "test_polygon" do
+    linear_ring1 = LinearRing.from_coordinates([[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]],256)
+    linear_ring2 = LinearRing.from_coordinates([[2.4,5.3],[5.4,1.4263],[14.46,1.06],[2.4,5.3]],256)
+    polygon = Polygon.from_linear_rings([linear_ring1,linear_ring2],256)
+
+    polygon.as_georss(:georss_ns => "hoyoyo").gsub("\n","").should eql("<hoyoyo:polygon>-45.3 12.4 41.6 45.4 1.0698 4.456 -45.3 12.4</hoyoyo:polygon>")
+    polygon.as_georss(:dialect => :w3cgeo, :w3cgeo_ns => "bouyoul").gsub("\n","").should eql("<bouyoul:lat>-45.3</bouyoul:lat><bouyoul:long>12.4</bouyoul:long>")
+    polygon.as_georss(:dialect => :gml).gsub("\n","").should eql("<georss:where><gml:Polygon><gml:exterior><gml:LinearRing><gml:posList>-45.3 12.4 41.6 45.4 1.0698 4.456 -45.3 12.4</gml:posList></gml:LinearRing></gml:exterior></gml:Polygon></georss:where>")
+
+    polygon.as_kml.gsub("\n","").should eql("<Polygon><outerBoundaryIs><LinearRing><coordinates>12.4,-45.3 45.4,41.6 4.456,1.0698 12.4,-45.3</coordinates></LinearRing></outerBoundaryIs><innerBoundaryIs><LinearRing><coordinates>2.4,5.3 5.4,1.4263 14.46,1.06 2.4,5.3</coordinates></LinearRing></innerBoundaryIs></Polygon>")
+  end
+
+  it "test_geometry_collection" do
+    gc = GeometryCollection.from_geometries([Point.from_x_y(4.67,45.4,256),LineString.from_coordinates([[5.7,12.45],[67.55,54]],256)],256)
+
+    #only the first geometry is output
+    gc.as_georss(:dialect => :simple,:floor => 4).gsub("\n","").should eql("<georss:point floor=\"4\">45.4 4.67</georss:point>")
+    gc.as_georss(:dialect => :w3cgeo).gsub("\n","").should eql("<geo:lat>45.4</geo:lat><geo:long>4.67</geo:long>")
+    gc.as_georss(:dialect => :gml).gsub("\n","").should eql("<georss:where><gml:Point><gml:pos>45.4 4.67</gml:pos></gml:Point></georss:where>")
+
+    gc.as_kml(:id => "HOYOYO-42").gsub("\n","").should eql("<MultiGeometry id=\"HOYOYO-42\"><Point><coordinates>4.67,45.4</coordinates></Point><LineString><coordinates>5.7,12.45 67.55,54</coordinates></LineString></MultiGeometry>")
+  end
+
+  it "test_envelope" do
+    linear_ring1 = LinearRing.from_coordinates([[12,-45,5],[45,41,6],[4,1,8],[12.4,-45,3]],256,true)
+    linear_ring2 = LinearRing.from_coordinates([[2,5,9],[5.4,1,-5.4],[14,1,34],[2,5,3]],256,true)
+    polygon = Polygon.from_linear_rings([linear_ring1,linear_ring2],256,true)
+
+    e = polygon.envelope
+
+    e.as_georss(:dialect => :simple).gsub("\n","").should eql("<georss:box>-45 4 41 45</georss:box>")
+    #center
+    e.as_georss(:dialect => :w3cgeo).gsub("\n","").should eql("<geo:lat>-2</geo:lat><geo:long>24</geo:long>")
+    e.as_georss(:dialect => :gml).gsub("\n","").should eql("<georss:where><gml:Envelope><gml:LowerCorner>-45 4</gml:LowerCorner><gml:UpperCorner>41 45</gml:UpperCorner></gml:Envelope></georss:where>")
+
+    e.as_kml.gsub("\n","").should eql("<LatLonAltBox><north>41</north><south>-45</south><east>45</east><west>4</west><minAltitude>-5.4</minAltitude><maxAltitude>34</maxAltitude></LatLonAltBox>")
+  end
+
+  it "test_point_georss_read" do
+    #W3CGeo
+    str = "   <geo:lat >12.3</geo:lat >\n\t  <geo:long>   4.56</geo:long> "
+    geom = Geometry.from_georss(str)
+    geom.class.should eql(Point)
+    geom.lat.should eql(12.3)
+    geom.lon.should eql(4.56)
+
+    str = " <geo:Point> \n \t  <geo:long>   4.56</geo:long> \n\t  <geo:lat >12.3</geo:lat > </geo:Point>  "
+    geom = Geometry.from_georss(str)
+    geom.class.should eql(Point)
+    geom.lat.should eql(12.3)
+    geom.lon.should eql(4.56)
+
+    #gml
+    str = " <georss:where> \t\r  <gml:Point  > \t <gml:pos> 4 \t 3 </gml:pos> </gml:Point> </georss:where>"
+    geom = Geometry.from_georss(str)
+    geom.class.should eql(Point)
+    geom.lat.should eql(4.0)
+    geom.lon.should eql(3.0)
+
+    #simple
+    str = "<georss:point > 4 \r\t  3 \t</georss:point >"
+    geom  = Geometry.from_georss(str)
+    geom.class.should eql(Point)
+    geom.lat.should eql(4.0)
+    geom.lon.should eql(3.0)
+
+    #simple with tags
+    str = "<georss:point featuretypetag=\"hoyoyo\"  elev=\"45.7\" \n floor=\"2\" relationshiptag=\"puyopuyo\" radius=\"42\" > 4 \n 3 \t</georss:point >"
+    geom,tags = Geometry.from_georss_with_tags(str)
+    geom.class.should eql(Point)
+    geom.lat.should eql(4.0)
+    geom.lon.should eql(3.0)
+    tags.featuretypetag.should eql("hoyoyo")
+    tags.elev.should eql(45.7)
+    tags.relationshiptag.should eql("puyopuyo")
+    tags.floor.should eql(2)
+    tags.radius.should eql(42.0)
+  end
+
+  it "test_line_string_georss_read" do
+    ls = LineString.from_points([Point.from_lon_lat(12.4,-45.3),Point.from_lon_lat(45.4,41.6)])
+
+    str = "<georss:line > -45.3 12.4 \n \r41.6\t 45.4</georss:line>"
+    geom  = Geometry.from_georss(str)
+    geom.class.should eql(LineString)
+    ls.should == geom
+
+    str = "<georss:where><gml:LineString><gml:posList>-45.3 12.4 41.6 45.4</gml:posList></gml:LineString></georss:where>"
+    geom  = Geometry.from_georss(str)
+    geom.class.should eql(LineString)
+    ls.should == geom
+  end
+
+  it "test_polygon_georss_read" do
+    linear_ring = LinearRing.from_coordinates([[12.4,-45.3],[45.4,41.6],[4.456,1.0698],[12.4,-45.3]])
+    polygon = Polygon.from_linear_rings([linear_ring])
+
+    str = "<hoyoyo:polygon featuretypetag=\"42\"  > -45.3 12.4 41.6 \n\r 45.4 1.0698 \r 4.456 -45.3 12.4 </hoyoyo:polygon>"
+    geom = Geometry.from_georss(str)
+    geom.class.should eql(Polygon)
+    polygon.should == geom
+
+    str = "<georss:where>\r\t \n  <gml:Polygon><gml:exterior>   <gml:LinearRing><gml:posList> -45.3 \n\r 12.4 41.6 \n\t 45.4 1.0698 4.456 -45.3 12.4</gml:posList></gml:LinearRing></gml:exterior></gml:Polygon></georss:where>"
+    geom = Geometry.from_georss(str)
+    geom.class.should eql(Polygon)
+    polygon.should == geom
+  end
+
+  it "test_envelope_georss_read" do
+    e = Envelope.from_coordinates([[4.456,-45.3],[45.4,41.6]])
+
+    str = "<georss:box  >-45.3 4.456 \n41.6 45.4</georss:box>"
+    geom = Geometry.from_georss(str)
+    geom.class.should eql(Envelope)
+    geom.should == e
+
+    str = "<georss:where><gml:Envelope><gml:lowerCorner>-45.3 \n 4.456</gml:lowerCorner><gml:upperCorner>41.6 \t\n 45.4</gml:upperCorner></gml:Envelope></georss:where>"
+    geom = Geometry.from_georss(str)
+    geom.class.should eql(Envelope)
+    geom.should == e
+  end
+
+  it "test_kml_read" do
+    g = Geometry.from_kml("<Point><coordinates>45,12,25</coordinates></Point>")
+    g.should be_a Point
+    g.should == Point.from_x_y_z(45,12,25)
+
+    g = Geometry.from_kml("<LineString>
+      <extrude>1</extrude>
+      <tessellate>1</tessellate>
+      <coordinates>
+        -122.364383,37.824664,0 -122.364152,37.824322,0
+      </coordinates>
+    </LineString>")
+    g.should be_a LineString
+    g.length.should eql(2)
+    g.should == LineString.from_points([Point.from_x_y_z(-122.364383,37.824664,0),Point.from_x_y_z(-122.364152,37.824322,0)],4326,true)
+
+    g = Geometry.from_kml("<Polygon>
+      <extrude>1</extrude>
+      <altitudeMode>relativeToGround</altitudeMode>
+      <outerBoundaryIs>
+        <LinearRing>
+          <coordinates>
+            -122.366278,37.818844,30
+            -122.365248,37.819267,30
+            -122.365640,37.819861,30
+            -122.366669,37.819429,30
+            -122.366278,37.818844,30
+          </coordinates>
+        </LinearRing>
+      </outerBoundaryIs>
+      <innerBoundaryIs>
+        <LinearRing>
+          <coordinates>
+            -122.366212,37.818977,30
+            -122.365424,37.819294,30
+            -122.365704,37.819731,30
+            -122.366488,37.819402,30
+            -122.366212,37.818977,30
+          </coordinates>
+        </LinearRing>
+      </innerBoundaryIs>
+      <innerBoundaryIs>
+        <LinearRing>
+          <coordinates>
+            -122.366212,37.818977,30
+            -122.365424,37.819294,30
+            -122.365704,37.819731,30
+            -122.366488,37.819402,30
+            -122.366212,37.818977,30
+          </coordinates>
+        </LinearRing>
+      </innerBoundaryIs>
+    </Polygon>")
+    g.should be_a Polygon
+    g.length.should eql(3)
+  end
+
+  it "test_to_kml_for_point_does_not_raise_type_error_if_geom_data_not_provided" do
+    point = Point.from_coordinates([1.6,2.8],123)
+    lambda { point.kml_representation }.should_not raise_error(TypeError)
+  end
+
+  it "test_to_kml_for_polygon_does_not_raise_type_error_if_geom_data_not_provided" do
+    polygon =  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)
+    lambda { polygon.kml_representation }.should_not raise_error(TypeError)
+  end
+
+  it "test_to_kml_for_line_string_does_not_raise_type_error_if_geom_data_not_provided" do
+    ls = LineString.from_coordinates([[5.7,12.45],[67.55,54]],256)
+    lambda { ls.kml_representation }.should_not raise_error(TypeError)
+  end
+
+end

data/spec/geo_ruby/georss_spec.rb

@@ -0,0 +1,14 @@
+require File.expand_path(File.dirname(__FILE__) + '/../spec_helper')
+
+RSS_DATA_DIR = File.dirname(__FILE__) + '/../data/georss/'
+
+describe GeorssParser do
+
+  it "should parse an rss file" do
+    geo = GeorssParser.new.parse(File.read(RSS_DATA_DIR + "/w3c.xml"))
+    geo.should be_a Point
+  end
+
+
+
+end

data/spec/geo_ruby/gpx4r/gpx_spec.rb

@@ -0,0 +1,106 @@
+require File.expand_path(File.dirname(__FILE__) + '/../../spec_helper')
+
+include GeoRuby::Gpx4r
+include GeoRuby::SimpleFeatures
+
+describe Gpx4r do
+
+  it "should add gpx extension and raise if doesn't exists" do
+    lambda do
+      File.should_receive(:exists?).with("short.gpx").and_return(false)
+      GpxFile.open('short').should be_true
+    end.should raise_error MalformedGpxException
+  end
+
+  describe "Waypoints" do
+
+    before(:all) do
+      @gpxfile = GpxFile.open(File.dirname(__FILE__) + '/../../data/gpx/short.gpx', :with_z => true, :with_m => true)
+      @gpxfile2 = GpxFile.open(File.dirname(__FILE__) + '/../../data/gpx/fells_loop', :with_z => true, :with_m => true)
+      @gpxfile3 = GpxFile.open(File.dirname(__FILE__) + '/../../data/gpx/tracktreks.gpx', :with_z => true)
+    end
+
+    it "should open and parse" do
+      @gpxfile.record_count.should eql(2724)
+    end
+
+    it "should open and parse no trkpt one" do
+      @gpxfile2.record_count.should eql(86)
+    end
+
+    it "should open and parse 3" do
+      @gpxfile3.record_count.should eql(225)
+    end
+
+    it "should read X and Y" do
+      @gpxfile[0].x.should be_within(0.0001).of(9.093942)
+      @gpxfile[0].y.should be_within(0.0001).of(48.731813)
+    end
+
+    it "should read Z and M" do
+      @gpxfile[0].z.should eql(468.0)
+      @gpxfile[0].m.should eql("2008-09-07T17:36:57Z")
+    end
+
+    it "should read X and Y 2" do
+      @gpxfile2[0].x.should be_within(0.0001).of(-71.119277)
+      @gpxfile2[0].y.should be_within(0.0001).of(42.438878)
+    end
+
+    it "should read Z and M 2" do
+      @gpxfile2[0].z.should eql(44.586548)
+      @gpxfile2[0].m.should eql("2001-11-28T21:05:28Z")
+    end
+
+    it "should read X and Y 3" do
+      @gpxfile3[0].x.should be_within(0.0001).of(-149.8358011)
+      @gpxfile3[0].y.should be_within(0.0001).of(-17.5326508)
+    end
+
+    it "should read Z and M 3" do
+      @gpxfile3[0].z.should eql(88.5787354)
+      @gpxfile3[0].m.should eql(0.0)
+    end
+
+    it "should return it as a linestring" do
+      @gpxfile.as_line_string.should be_instance_of LineString
+      @gpxfile.as_polyline.should be_instance_of LineString
+    end
+
+    it "should return it as a linestring 3" do
+      @gpxfile3.as_line_string.should be_instance_of LineString
+      @gpxfile3.as_polyline.should be_instance_of LineString
+    end
+
+    it "should return a envelope" do
+      @gpxfile.envelope.should be_instance_of Envelope
+      @gpxfile.envelope.lower_corner.x.should be_within(0.001).of(9.08128)
+      @gpxfile.envelope.lower_corner.y.should be_within(0.001).of(48.7169)
+    end
+
+    it "should return a envelope 3" do
+      @gpxfile3.envelope.should be_instance_of Envelope
+      @gpxfile3.envelope.lower_corner.x.should be_within(0.001).of(-149.8422613)
+      @gpxfile3.envelope.lower_corner.y.should be_within(0.001).of(-17.547636)
+    end
+
+    it "should return it as a polygon" do
+      [@gpxfile, @gpxfile2, @gpxfile3].each do |g|
+        g.as_polygon.should be_instance_of Polygon
+        g.as_polygon[0].should be_instance_of LinearRing
+        g.as_polygon[0].should be_closed
+        g.as_polygon[1].should be_nil
+      end
+    end
+
+    it "should close the polygon" do
+      se = Point.from_x_y(-44, -23)
+      sw = Point.from_x_y(-42, -22)
+      nw = Point.from_x_y(-42, -25)
+      ne = Point.from_x_y(-44, -21)
+      @gpxfile.instance_variable_set(:@points, [se,sw,nw,ne])
+      @gpxfile.as_polygon.should == Polygon.from_points([[se,sw,nw,ne,se]])
+    end
+  end
+
+end