schleyfox-rgeo 0.2.5

data/lib/rgeo/geos/impl_additions.rb

@@ -0,0 +1,76 @@
+# -----------------------------------------------------------------------------
+# 
+# GEOS implementation additions written in Ruby
+# 
+# -----------------------------------------------------------------------------
+# Copyright 2010 Daniel Azuma
+# 
+# All rights reserved.
+# 
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions are met:
+# 
+# * Redistributions of source code must retain the above copyright notice,
+#   this list of conditions and the following disclaimer.
+# * Redistributions in binary form must reproduce the above copyright notice,
+#   this list of conditions and the following disclaimer in the documentation
+#   and/or other materials provided with the distribution.
+# * Neither the name of the copyright holder, nor the names of any other
+#   contributors to this software, may be used to endorse or promote products
+#   derived from this software without specific prior written permission.
+# 
+# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+# AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+# IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+# LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR # CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+# SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+# CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+# ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+# POSSIBILITY OF SUCH DAMAGE.
+# -----------------------------------------------------------------------------
+;
+
+
+module RGeo
+  
+  module Geos
+    
+    
+    class GeometryImpl  # :nodoc:
+      
+      include Feature::Instance
+      
+      def inspect
+        "#<#{self.class}:0x#{object_id.to_s(16)} #{as_text.inspect}>"
+      end
+      
+    end
+    
+    
+    class Factory
+      
+      
+      # :stopdoc:
+      if ::RGeo::Geos.supported?
+        IMPL_CLASSES = {
+          Feature::Point => PointImpl,
+          Feature::LineString => LineStringImpl,
+          Feature::LinearRing => LinearRingImpl,
+          Feature::Line => LineImpl,
+          Feature::GeometryCollection => GeometryCollectionImpl,
+          Feature::MultiPoint => MultiPointImpl,
+          Feature::MultiLineString => MultiLineStringImpl,
+          Feature::MultiPolygon => MultiPolygonImpl,
+        }
+      end
+      # :startdoc:
+      
+      
+    end
+    
+    
+  end
+  
+end

data/lib/rgeo/geos/interface.rb

@@ -0,0 +1,139 @@
+# -----------------------------------------------------------------------------
+# 
+# GEOS toplevel interface
+# 
+# -----------------------------------------------------------------------------
+# Copyright 2010 Daniel Azuma
+# 
+# All rights reserved.
+# 
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions are met:
+# 
+# * Redistributions of source code must retain the above copyright notice,
+#   this list of conditions and the following disclaimer.
+# * Redistributions in binary form must reproduce the above copyright notice,
+#   this list of conditions and the following disclaimer in the documentation
+#   and/or other materials provided with the distribution.
+# * Neither the name of the copyright holder, nor the names of any other
+#   contributors to this software, may be used to endorse or promote products
+#   derived from this software without specific prior written permission.
+# 
+# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+# AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+# IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+# LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+# CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+# SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+# CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+# ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+# POSSIBILITY OF SUCH DAMAGE.
+# -----------------------------------------------------------------------------
+;
+
+
+module RGeo
+  
+  module Geos
+    
+    @supported = nil
+    
+    class << self
+      
+      
+      # Returns true if GEOS implementation is supported.
+      # If this returns false, GEOS features are not available.
+      
+      def supported?
+        @supported.nil? ? (@supported = Factory.respond_to?(:_create)) : @supported
+      end
+      
+      
+      # Returns true if the given feature is a GEOS feature, or if the given
+      # factory is a GEOS factory.
+      
+      def is_geos?(object_)
+        supported? && (Factory === object_ || GeometryImpl === object_ || ZMFactory === object_ || ZMGeometryImpl === object_)
+      end
+      
+      
+      # Returns a factory for the GEOS implementation.
+      # Returns nil if the GEOS implementation is not supported.
+      # 
+      # Note that GEOS does not natively support 4-dimensional data
+      # (i.e. both z and m values). However, RGeo's GEOS wrapper does
+      # provide a 4-dimensional factory that utilizes an extra native
+      # GEOS object to handle the extra coordinate. Hence, a factory
+      # configured with both Z and M support will work, but will be
+      # slower than a 2-dimensional or 3-dimensional factory.
+      # 
+      # Options include:
+      # 
+      # [<tt>:lenient_multi_polygon_assertions</tt>]
+      #   If set to true, assertion checking on MultiPolygon is disabled.
+      #   This may speed up creation of MultiPolygon objects, at the
+      #   expense of not doing the proper checking for OGC MultiPolygon
+      #   compliance. See RGeo::Feature::MultiPolygon for details on
+      #   the MultiPolygon assertions. Default is false.
+      # [<tt>:buffer_resolution</tt>]
+      #   The resolution of buffers around geometries created by this
+      #   factory. This controls the number of line segments used to
+      #   approximate curves. The default is 1, which causes, for
+      #   example, the buffer around a point to be approximated by a
+      #   4-sided polygon. A resolution of 2 would cause that buffer
+      #   to be approximated by an 8-sided polygon. The exact behavior
+      #   for different kinds of buffers is defined by GEOS.
+      # [<tt>:srid</tt>]
+      #   Set the SRID returned by geometries created by this factory.
+      #   Default is 0.
+      # [<tt>:proj4</tt>]
+      #   The coordinate system in Proj4 format, either as a
+      #   CoordSys::Proj4 object or as a string or hash representing the
+      #   proj4 format. Optional.
+      # [<tt>:coord_sys</tt>]
+      #   The coordinate system in OGC form, either as a subclass of
+      #   CoordSys::CS::CoordinateSystem, or as a string in WKT format.
+      #   Optional.
+      # [<tt>:srs_database</tt>]
+      #   Optional. If provided, the value should be an implementation of
+      #   CoordSys::SRSDatabase::Interface. If both this and an SRID are
+      #   provided, they are used to look up the proj4 and coord_sys
+      #   objects from a spatial reference system database.
+      # [<tt>:has_z_coordinate</tt>]
+      #   Support <tt>z_coordinate</tt>. Default is false.
+      # [<tt>:has_m_coordinate</tt>]
+      #   Support <tt>m_coordinate</tt>. Default is false.
+      
+      def factory(opts_={})
+        if supported?
+          if opts_[:has_z_coordinate] && opts_[:has_m_coordinate]
+            ZMFactory.new(opts_)
+          else
+            Factory.create(opts_)
+          end
+        else
+          nil
+        end
+      end
+      
+      
+      # Returns a Feature::FactoryGenerator that creates Geos-backed
+      # factories. The given options are used as the default options.
+      # 
+      # A common case for this is to provide the <tt>:srs_database</tt>
+      # as a default. Then, the factory generator need only be passed
+      # an SRID and it will automatically fetch the appropriate Proj4
+      # and CoordSys objects.
+      
+      def factory_generator(defaults_={})
+        ::Proc.new{ |c_| factory(defaults_.merge(c_)) }
+      end
+      
+      
+    end
+    
+  end
+  
+end

data/lib/rgeo/geos/zm_factory.rb

@@ -0,0 +1,275 @@
+# -----------------------------------------------------------------------------
+# 
+# GEOS zm factory implementation
+# 
+# -----------------------------------------------------------------------------
+# Copyright 2010 Daniel Azuma
+# 
+# All rights reserved.
+# 
+# Redistribution and use in source and binary forms, with or without
+# modification, are permitted provided that the following conditions are met:
+# 
+# * Redistributions of source code must retain the above copyright notice,
+#   this list of conditions and the following disclaimer.
+# * Redistributions in binary form must reproduce the above copyright notice,
+#   this list of conditions and the following disclaimer in the documentation
+#   and/or other materials provided with the distribution.
+# * Neither the name of the copyright holder, nor the names of any other
+#   contributors to this software, may be used to endorse or promote products
+#   derived from this software without specific prior written permission.
+# 
+# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+# AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+# IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
+# ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
+# LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
+# CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
+# SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
+# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
+# CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
+# ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
+# POSSIBILITY OF SUCH DAMAGE.
+# -----------------------------------------------------------------------------
+;
+
+
+module RGeo
+  
+  module Geos
+    
+    
+    # A factory for Geos that handles both Z and M.
+    
+    class ZMFactory
+      
+      include Feature::Factory::Instance
+      
+      
+      class << self
+        
+        
+        # Create a new factory. Returns nil if the GEOS implementation is
+        # not supported.
+        
+        def create(opts_={})
+          return nil unless Geos.supported?
+          new(opts_)
+        end
+        
+        
+      end
+      
+      
+      def initialize(opts_={})  # :nodoc:
+        proj4_ = opts_[:proj4]
+        coord_sys_ = opts_[:coord_sys]
+        srid_ = opts_[:srid]
+        if (!proj4_ || !coord_sys_) && srid_ && (db_ = opts_[:srs_database])
+          entry_ = db_.get(srid_.to_i)
+          if entry_
+            proj4_ ||= entry_.proj4
+            coord_sys_ ||= entry_.coord_sys
+          end
+        end
+        srid_ ||= coord_sys_.authority_code if coord_sys_
+        config_ = {
+          :lenient_multi_polygon_assertions => opts_[:lenient_multi_polygon_assertions],
+          :buffer_resolution => opts_[:buffer_resolution],
+          :srid => srid_.to_i, :proj4 => proj4_, :coord_sys => coord_sys_,
+        }
+        @zfactory = Factory.create(config_.merge(:has_z_coordinate => true))
+        @mfactory = Factory.create(config_.merge(:has_m_coordinate => true))
+      end
+      
+      
+      # Returns the SRID of geometries created by this factory.
+      
+      def srid
+        @zfactory._srid
+      end
+      
+      
+      # Returns the resolution used by buffer calculations on geometries
+      # created by this factory
+      
+      def buffer_resolution
+        @zfactory._buffer_resolution
+      end
+      
+      
+      # Returns true if this factory is lenient with MultiPolygon assertions
+      
+      def lenient_multi_polygon_assertions?
+        @zfactory.lenient_multi_polygon_assertions?
+      end
+      
+      
+      # Returns the z-only factory corresponding to this factory.
+      
+      def z_factory
+        @zfactory
+      end
+      
+      
+      # Returns the m-only factory corresponding to this factory.
+      
+      def m_factory
+        @mfactory
+      end
+      
+      
+      # Factory equivalence test.
+      
+      def eql?(rhs_)
+        rhs_.is_a?(ZMFactory) && rhs_.z_factory == @zfactory
+      end
+      alias_method :==, :eql?
+      
+      
+      # See ::RGeo::Feature::Factory#property
+      
+      def property(name_)
+        case name_
+        when :has_z_coordinate, :has_m_coordinate, :is_cartesian
+          true
+        else
+          nil
+        end
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#parse_wkt
+      
+      def parse_wkt(str_)
+        WKRep::WKTParser.new(self).parse(str_)
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#parse_wkb
+      
+      def parse_wkb(str_)
+        WKRep::WKBParser.new(self).parse(str_)
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#point
+      
+      def point(x_, y_, z_=0, m_=0)
+        ZMPointImpl.create(self, @zfactory.point(x_, y_, z_), @mfactory.point(x_, y_, m_))
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#line_string
+      
+      def line_string(points_)
+        ZMLineStringImpl.create(self, @zfactory.line_string(points_), @mfactory.line_string(points_))
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#line
+      
+      def line(start_, end_)
+        ZMLineStringImpl.create(self, @zfactory.line(start_, end_), @mfactory.line(start_, end_))
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#linear_ring
+      
+      def linear_ring(points_)
+        ZMLineStringImpl.create(self, @zfactory.linear_ring(points_), @mfactory.linear_ring(points_))
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#polygon
+      
+      def polygon(outer_ring_, inner_rings_=nil)
+        ZMPolygonImpl.create(self, @zfactory.polygon(outer_ring_, inner_rings_), @mfactory.polygon(outer_ring_, inner_rings_))
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#collection
+      
+      def collection(elems_)
+        ZMGeometryCollectionImpl.create(self, @zfactory.collection(elems_), @mfactory.collection(elems_))
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#multi_point
+      
+      def multi_point(elems_)
+        ZMGeometryCollectionImpl.create(self, @zfactory.multi_point(elems_), @mfactory.multi_point(elems_))
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#multi_line_string
+      
+      def multi_line_string(elems_)
+        ZMMultiLineStringImpl.create(self, @zfactory.multi_line_string(elems_), @mfactory.multi_line_string(elems_))
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#multi_polygon
+      
+      def multi_polygon(elems_)
+        ZMMultiPolygonImpl.create(self, @zfactory.multi_polygon(elems_), @mfactory.multi_polygon(elems_))
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#proj4
+      
+      def proj4
+        @zfactory.proj4
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#coord_sys
+      
+      def coord_sys
+        @zfactory.coord_sys
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#override_cast
+      
+      def override_cast(original_, ntype_, flags_)
+        return nil unless Geos.supported?
+        keep_subtype_ = flags_[:keep_subtype]
+        force_new_ = flags_[:force_new]
+        project_ = flags_[:project]
+        type_ = original_.geometry_type
+        ntype_ = type_ if keep_subtype_ && type_.include?(ntype_)
+        case original_
+        when ZMGeometryImpl
+          # Optimization if we're just changing factories, but to
+          # another ZM factory.
+          if original_.factory != self && ntype_ == type_ &&
+              (!project_ || original_.factory.proj4 == @proj4)
+          then
+            zresult_ = original_.z_geometry.dup
+            zresult_._set_factory(@zfactory)
+            mresult_ = original_.m_geometry.dup
+            mresult_._set_factory(@mfactory)
+            return original_.class.create(self, zresult_, mresult_)
+          end
+          # LineString conversion optimization.
+          if (original_.factory != self || ntype_ != type_) &&
+              (!project_ || original_.factory.proj4 == @proj4) &&
+              type_.subtype_of?(Feature::LineString) && ntype_.subtype_of?(Feature::LineString)
+          then
+            klass_ = Factory::IMPL_CLASSES[ntype_]
+            zresult_ = klass_._copy_from(@zfactory, original_.z_geometry)
+            mresult_ = klass_._copy_from(@mfactory, original_.m_geometry)
+            return ZMLineStringImpl.create(self, zresult_, mresult_)
+          end
+        end
+        false
+      end
+      
+      
+    end
+    
+    
+  end
+  
+end