rgeo-dschee 0.5.4

data/test/coord_sys/proj4_srs_data_test.rb

@@ -0,0 +1,41 @@
+# -----------------------------------------------------------------------------
+#
+# Tests for OGC CS classes
+#
+# -----------------------------------------------------------------------------
+
+require "test/unit"
+require "rgeo"
+
+module RGeo
+  module Tests # :nodoc:
+    module CoordSys # :nodoc:
+      class TestProj4SRSData < ::Test::Unit::TestCase # :nodoc:
+        def test_epsg_4326
+          db_ = ::RGeo::CoordSys::SRSDatabase::Proj4Data.new("epsg")
+          entry_ = db_.get(4326)
+          allowed_vals_ = [
+            "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", # Proj 4.7
+            "+proj=longlat +datum=WGS84 +no_defs" # Proj 4.8
+          ]
+          assert(allowed_vals_.include?(entry_.proj4.original_str))
+          assert_equal("WGS 84", entry_.name)
+        end
+
+        def test_epsg_3857
+          db_ = ::RGeo::CoordSys::SRSDatabase::Proj4Data.new("epsg")
+          entry_ = db_.get(3857)
+          assert_equal("+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs", entry_.proj4.original_str)
+          assert_equal("WGS 84 / Pseudo-Mercator", entry_.name)
+        end
+
+        def test_nad83_4601
+          db_ = ::RGeo::CoordSys::SRSDatabase::Proj4Data.new("nad83")
+          entry_ = db_.get(4601)
+          assert_equal("+proj=lcc  +datum=NAD83 +lon_0=-120d50 +lat_1=48d44 +lat_2=47d30 +lat_0=47 +x_0=500000 +y_0=0 +no_defs", entry_.proj4.original_str)
+          assert_equal("4601: washington north: nad83", entry_.name)
+        end
+      end
+    end
+  end
+end if ::RGeo::CoordSys::Proj4.supported?

data/test/coord_sys/proj4_test.rb

@@ -0,0 +1,150 @@
+# -----------------------------------------------------------------------------
+#
+# Tests for proj4 wrapper
+#
+# -----------------------------------------------------------------------------
+
+require "test/unit"
+require "rgeo"
+
+module RGeo
+  module Tests # :nodoc:
+    module CoordSys # :nodoc:
+      class TestProj4 < ::Test::Unit::TestCase # :nodoc:
+        def test_proj4_version
+          assert_match(/^\d+\.\d+(\.\d+)?$/, ::RGeo::CoordSys::Proj4.version)
+        end
+
+        def test_create_wgs84
+          proj_ = ::RGeo::CoordSys::Proj4.create("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+          assert_equal(true, proj_.geographic?)
+          assert_equal("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", proj_.original_str)
+          assert_equal(" +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0", proj_.canonical_str)
+        end
+
+        def test_get_wgs84_geographic
+          proj_ = ::RGeo::CoordSys::Proj4.create("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+          proj2_ = proj_.get_geographic
+          assert_nil(proj2_.original_str)
+          assert_equal(true, proj2_.geographic?)
+          coords_ = ::RGeo::CoordSys::Proj4.transform_coords(proj_, proj2_, 1, 2, 0)
+          assert_equal([1, 2, 0], coords_)
+        end
+
+        def test_identity_transform
+          proj_ = ::RGeo::CoordSys::Proj4.create("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+          assert_equal([1, 2, 0], ::RGeo::CoordSys::Proj4.transform_coords(proj_, proj_, 1, 2, 0))
+          assert_equal([1, 2], ::RGeo::CoordSys::Proj4.transform_coords(proj_, proj_, 1, 2, nil))
+        end
+
+        def _project_merc(x_, y_)
+          [x_ * 6_378_137.0, ::Math.log(::Math.tan(::Math::PI / 4.0 + y_ / 2.0)) * 6_378_137.0]
+        end
+
+        def _unproject_merc(x_, y_)
+          [x_ / 6_378_137.0, (2.0 * ::Math.atan(::Math.exp(y_ / 6_378_137.0)) - ::Math::PI / 2.0)]
+        end
+
+        def _assert_close_enough(a_, b_)
+          delta_ = ::Math.sqrt(a_ * a_ + b_ * b_) * 0.00000001
+          delta_ = 0.000000000001 if delta_ < 0.000000000001
+          assert_in_delta(a_, b_, delta_)
+        end
+
+        def _assert_xy_close(xy1_, xy2_)
+          _assert_close_enough(xy1_[0], xy2_[0])
+          _assert_close_enough(xy1_[1], xy2_[1])
+        end
+
+        def test_simple_mercator_transform
+          geography_ = ::RGeo::CoordSys::Proj4.create("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", radians: true)
+          projection_ = ::RGeo::CoordSys::Proj4.create("+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs")
+          _assert_xy_close(_project_merc(0, 0), ::RGeo::CoordSys::Proj4.transform_coords(geography_, projection_, 0, 0, nil))
+          _assert_xy_close(_project_merc(0.01, 0.01), ::RGeo::CoordSys::Proj4.transform_coords(geography_, projection_, 0.01, 0.01, nil))
+          _assert_xy_close(_project_merc(1, 1), ::RGeo::CoordSys::Proj4.transform_coords(geography_, projection_, 1, 1, nil))
+          _assert_xy_close(_project_merc(-1, -1), ::RGeo::CoordSys::Proj4.transform_coords(geography_, projection_, -1, -1, nil))
+          _assert_xy_close(_unproject_merc(0, 0), ::RGeo::CoordSys::Proj4.transform_coords(projection_, geography_, 0, 0, nil))
+          _assert_xy_close(_unproject_merc(10_000, 10_000), ::RGeo::CoordSys::Proj4.transform_coords(projection_, geography_, 10_000, 10_000, nil))
+          _assert_xy_close(_unproject_merc(-20_000_000, -20_000_000), ::RGeo::CoordSys::Proj4.transform_coords(projection_, geography_, -20_000_000, -20_000_000, nil))
+        end
+
+        def test_equivalence
+          proj1_ = ::RGeo::CoordSys::Proj4.create("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+          proj2_ = ::RGeo::CoordSys::Proj4.create(" +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+          assert_equal(proj1_, proj2_)
+        end
+
+        def test_hashes_equal_for_equivalent_objects
+          proj1_ = ::RGeo::CoordSys::Proj4.create("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+          proj2_ = ::RGeo::CoordSys::Proj4.create(" +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+          assert_equal(proj1_.hash, proj2_.hash)
+        end
+
+        def test_point_projection_cast
+          geography_ = ::RGeo::Geos.factory(proj4: "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", srid: 4326)
+          projection_ = ::RGeo::Geos.factory(proj4: "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs", srid: 27_700)
+          proj_point_ = projection_.parse_wkt("POINT(473600.5000000000000000 186659.7999999999883585)")
+          geo_point_ = ::RGeo::Feature.cast(proj_point_, project: true, factory: geography_)
+          _assert_close_enough(-0.9393598527244420, geo_point_.x)
+          _assert_close_enough(51.5740106527552697, geo_point_.y)
+        end
+
+        def test_point_transform_lowlevel
+          geography_ = ::RGeo::Geos.factory(proj4: "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", srid: 4326)
+          projection_ = ::RGeo::Geos.factory(proj4: "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs", srid: 27_700)
+          proj_point_ = projection_.parse_wkt("POINT(473600.5000000000000000 186659.7999999999883585)")
+          geo_point_ = ::RGeo::CoordSys::Proj4.transform(projection_.proj4, proj_point_, geography_.proj4, geography_)
+          _assert_close_enough(-0.9393598527244420, geo_point_.x)
+          _assert_close_enough(51.5740106527552697, geo_point_.y)
+        end
+
+        def test_geocentric
+          obj1_ = ::RGeo::CoordSys::Proj4.create("+proj=geocent +ellps=WGS84")
+          assert_equal(true, obj1_.geocentric?)
+        end
+
+        def test_get_geographic
+          projection_ = ::RGeo::CoordSys::Proj4.create("+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs")
+          geographic_ = projection_.get_geographic
+          expected_ = ::RGeo::CoordSys::Proj4.create("+proj=latlong +a=6378137 +b=6378137 +nadgrids=@null")
+          assert_equal(expected_, geographic_)
+        end
+
+        def test_marshal_roundtrip
+          obj1_ = ::RGeo::CoordSys::Proj4.create("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+          dump_ = ::Marshal.dump(obj1_)
+          obj2_ = ::Marshal.load(dump_)
+          assert_equal(obj1_, obj2_)
+        end
+
+        def test_dup
+          obj1_ = ::RGeo::CoordSys::Proj4.create("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+          obj2_ = obj1_.dup
+          assert_equal(obj1_, obj2_)
+        end
+
+        def test_dup_of_get_geographic
+          obj1_ = ::RGeo::CoordSys::Proj4.create("+proj=latlong +datum=WGS84 +ellps=WGS84")
+          obj2_ = obj1_.get_geographic
+          obj3_ = obj2_.dup
+          assert_equal(obj1_, obj3_)
+        end
+
+        if ::RGeo.yaml_supported?
+
+          def test_yaml_roundtrip
+            obj1_ = ::RGeo::CoordSys::Proj4.create("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
+            dump_ = ::Psych.dump(obj1_)
+            obj2_ = ::Psych.load(dump_)
+            assert_equal(obj1_, obj2_)
+          end
+
+        end
+      end
+    end
+  end
+end if ::RGeo::CoordSys::Proj4.supported?
+
+unless ::RGeo::CoordSys::Proj4.supported?
+  puts "WARNING: Proj4 support not available. Related tests skipped."
+end

data/test/coord_sys/sr_org_test.rb

@@ -0,0 +1,32 @@
+# -----------------------------------------------------------------------------
+#
+# Tests for OGC CS classes
+#
+# -----------------------------------------------------------------------------
+
+require "test/unit"
+require "rgeo"
+
+module RGeo
+  module Tests # :nodoc:
+    module CoordSys # :nodoc:
+      class TestSrOrg < ::Test::Unit::TestCase # :nodoc:
+        def test_epsg_4326
+          db_ = ::RGeo::CoordSys::SRSDatabase::SrOrg.new("epsg")
+          entry_ = db_.get(4326)
+          assert_equal("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", entry_.proj4.original_str)
+          assert_kind_of(::RGeo::CoordSys::CS::GeographicCoordinateSystem, entry_.coord_sys)
+          assert_equal("WGS 84", entry_.name)
+        end
+
+        def test_epsg_3857
+          db_ = ::RGeo::CoordSys::SRSDatabase::SrOrg.new("epsg")
+          entry_ = db_.get(3857)
+          assert_equal("+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs", entry_.proj4.original_str)
+          assert_kind_of(::RGeo::CoordSys::CS::ProjectedCoordinateSystem, entry_.coord_sys)
+          assert_equal("Popular Visualisation CRS / Mercator", entry_.name)
+        end
+      end if false
+    end
+  end
+end

data/test/coord_sys/url_reader_test.rb

@@ -0,0 +1,42 @@
+# -----------------------------------------------------------------------------
+#
+# Tests for OGC CS classes
+#
+# -----------------------------------------------------------------------------
+
+require "test/unit"
+require "rgeo"
+
+module RGeo
+  module Tests # :nodoc:
+    module CoordSys # :nodoc:
+      class TestUrlReader < ::Test::Unit::TestCase # :nodoc:
+        def test_sr_org_epsg_4326_ogcwkt
+          db_ = ::RGeo::CoordSys::SRSDatabase::UrlReader.new
+          entry_ = db_.get("http://spatialreference.org/ref/epsg/4326/ogcwkt/")
+          assert_kind_of(::RGeo::CoordSys::CS::GeographicCoordinateSystem, entry_.coord_sys)
+          assert_equal("WGS 84", entry_.name)
+        end
+
+        def test_sr_org_epsg_4326_proj4
+          db_ = ::RGeo::CoordSys::SRSDatabase::UrlReader.new
+          entry_ = db_.get("http://spatialreference.org/ref/epsg/4326/proj4/")
+          assert_equal("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", entry_.proj4.original_str)
+        end
+
+        def test_sr_org_epsg_3857_ogcwkt
+          db_ = ::RGeo::CoordSys::SRSDatabase::UrlReader.new
+          entry_ = db_.get("http://spatialreference.org/ref/epsg/3857/ogcwkt/")
+          assert_kind_of(::RGeo::CoordSys::CS::ProjectedCoordinateSystem, entry_.coord_sys)
+          assert_equal("Popular Visualisation CRS / Mercator", entry_.name)
+        end
+
+        def test_sr_org_epsg_3857_proj4
+          db_ = ::RGeo::CoordSys::SRSDatabase::UrlReader.new
+          entry_ = db_.get("http://spatialreference.org/ref/epsg/3857/proj4/")
+          assert_equal("+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs", entry_.proj4.original_str)
+        end
+      end if false
+    end
+  end
+end

data/test/geos_capi/factory_test.rb

@@ -0,0 +1,31 @@
+# -----------------------------------------------------------------------------
+#
+# Tests for the GEOS factory
+#
+# -----------------------------------------------------------------------------
+
+require "test/unit"
+require "rgeo"
+
+require ::File.expand_path("../common/factory_tests.rb", ::File.dirname(__FILE__))
+
+module RGeo
+  module Tests # :nodoc:
+    module GeosCAPI # :nodoc:
+      class TestFactory < ::Test::Unit::TestCase # :nodoc:
+        def setup
+          @factory = ::RGeo::Geos.factory(srid: 1000)
+          @srid = 1000
+        end
+
+        include ::RGeo::Tests::Common::FactoryTests
+
+        def test_is_geos_factory
+          assert_equal(true, ::RGeo::Geos.is_geos?(@factory))
+          assert_equal(true, ::RGeo::Geos.is_capi_geos?(@factory))
+          assert_equal(false, ::RGeo::Geos.is_ffi_geos?(@factory))
+        end
+      end
+    end
+  end
+end if ::RGeo::Geos.capi_supported?

data/test/geos_capi/geometry_collection_test.rb

@@ -0,0 +1,24 @@
+# -----------------------------------------------------------------------------
+#
+# Tests for the GEOS geometry collection implementation
+#
+# -----------------------------------------------------------------------------
+
+require "test/unit"
+require "rgeo"
+
+require ::File.expand_path("../common/geometry_collection_tests.rb", ::File.dirname(__FILE__))
+
+module RGeo
+  module Tests # :nodoc:
+    module GeosCAPI # :nodoc:
+      class TestGeometryCollection < ::Test::Unit::TestCase # :nodoc:
+        def create_factory
+          ::RGeo::Geos.factory
+        end
+
+        include ::RGeo::Tests::Common::GeometryCollectionTests
+      end
+    end
+  end
+end if ::RGeo::Geos.capi_supported?

data/test/geos_capi/line_string_test.rb

@@ -0,0 +1,24 @@
+# -----------------------------------------------------------------------------
+#
+# Tests for the GEOS line string implementation
+#
+# -----------------------------------------------------------------------------
+
+require "test/unit"
+require "rgeo"
+
+require ::File.expand_path("../common/line_string_tests.rb", ::File.dirname(__FILE__))
+
+module RGeo
+  module Tests # :nodoc:
+    module GeosCAPI # :nodoc:
+      class TestLineString < ::Test::Unit::TestCase # :nodoc:
+        include ::RGeo::Tests::Common::LineStringTests
+
+        def setup
+          @factory = ::RGeo::Geos.factory
+        end
+      end
+    end
+  end
+end if ::RGeo::Geos.capi_supported?

data/test/geos_capi/misc_test.rb

@@ -0,0 +1,116 @@
+# -----------------------------------------------------------------------------
+#
+# Tests for miscellaneous GEOS stuff
+#
+# -----------------------------------------------------------------------------
+
+require "test/unit"
+require "rgeo"
+
+module RGeo
+  module Tests # :nodoc:
+    module GeosCAPI # :nodoc:
+      class TestMisc < ::Test::Unit::TestCase # :nodoc:
+        def setup
+          @factory = ::RGeo::Geos.factory(srid: 4326)
+        end
+
+        def test_marshal_dump_with_geos
+          @factory = ::RGeo::Geos.factory(
+            srid: 4326,
+            wkt_generator: :geos,
+            wkb_generator: :geos,
+            wkt_parser: :geos,
+            wkb_parser: :geos
+          )
+
+          dump = nil
+          assert_nothing_raised { dump = @factory.marshal_dump }
+          assert_equal({}, dump["wktg"])
+          assert_equal({}, dump["wkbg"])
+          assert_equal({}, dump["wktp"])
+          assert_equal({}, dump["wkbp"])
+        end
+
+        def test_encode_with_geos
+          @factory = ::RGeo::Geos.factory(
+            srid: 4326,
+            wkt_generator: :geos,
+            wkb_generator: :geos,
+            wkt_parser: :geos,
+            wkb_parser: :geos
+          )
+          coder = Psych::Coder.new("test")
+
+          assert_nothing_raised { @factory.encode_with(coder) }
+          assert_equal({}, coder["wkt_generator"])
+          assert_equal({}, coder["wkb_generator"])
+          assert_equal({}, coder["wkt_parser"])
+          assert_equal({}, coder["wkb_parser"])
+        end
+
+        def test_uninitialized
+          geom_ = ::RGeo::Geos::CAPIGeometryImpl.new
+          assert_equal(false, geom_.initialized?)
+          assert_nil(geom_.geometry_type)
+        end
+
+        def test_empty_geometries_equal
+          geom1_ = @factory.collection([])
+          geom2_ = @factory.line_string([])
+          assert(!geom1_.eql?(geom2_))
+          assert(geom1_.equals?(geom2_))
+        end
+
+        def test_prepare
+          p1_ = @factory.point(1, 2)
+          p2_ = @factory.point(3, 4)
+          p3_ = @factory.point(5, 2)
+          polygon_ = @factory.polygon(@factory.linear_ring([p1_, p2_, p3_, p1_]))
+          assert_equal(false, polygon_.prepared?)
+          polygon_.prepare!
+          assert_equal(true, polygon_.prepared?)
+        end
+
+        def test_auto_prepare
+          p1_ = @factory.point(1, 2)
+          p2_ = @factory.point(3, 4)
+          p3_ = @factory.point(5, 2)
+          polygon_ = @factory.polygon(@factory.linear_ring([p1_, p2_, p3_, p1_]))
+          assert_equal(false, polygon_.prepared?)
+          polygon_.intersects?(p1_)
+          assert_equal(false, polygon_.prepared?)
+          polygon_.intersects?(p2_)
+          assert_equal(true, polygon_.prepared?)
+
+          factory_no_auto_prepare_ = ::RGeo::Geos.factory(srid: 4326, auto_prepare: :disabled)
+          polygon2_ = factory_no_auto_prepare_.polygon(
+            factory_no_auto_prepare_.linear_ring([p1_, p2_, p3_, p1_]))
+          assert_equal(false, polygon2_.prepared?)
+          polygon2_.intersects?(p1_)
+          assert_equal(false, polygon2_.prepared?)
+          polygon2_.intersects?(p2_)
+          assert_equal(false, polygon2_.prepared?)
+        end
+
+        def test_gh_21
+          # Test for GH-21 (seg fault in rgeo_convert_to_geos_geometry)
+          # This seemed to fail under Ruby 1.8.7 only.
+          f_ = RGeo::Geographic.simple_mercator_factory
+          loc_ = f_.line_string([f_.point(-123, 37), f_.point(-122, 38)])
+          f2_ = f_.projection_factory
+          loc2_ = f2_.line_string([f2_.point(-123, 37), f2_.point(-122, 38)])
+          loc2_.intersection(loc_)
+        end
+
+        def test_geos_version
+          assert_match(/^\d+\.\d+(\.\d+)?$/, ::RGeo::Geos.version)
+        end
+      end
+    end
+  end
+end if ::RGeo::Geos.capi_supported?
+
+unless ::RGeo::Geos.capi_supported?
+  puts "WARNING: GEOS CAPI support not available. Related tests skipped."
+end