schleyfox-rgeo 0.2.5

data/lib/rgeo/cartesian/calculations.rb

@@ -0,0 +1,161 @@
+# -----------------------------------------------------------------------------
+# 
+# Core calculations in the plane
+# 
+# -----------------------------------------------------------------------------
+# 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 Cartesian
+    
+    
+    # Represents a line segment in the plane.
+    
+    class Segment  # :nodoc:
+      
+      def initialize(start_, end_)
+        @s = start_
+        @e = end_
+        @sx = @s.x
+        @sy = @s.y
+        @ex = @e.x
+        @ey = @e.y
+        @dx = @ex - @sx
+        @dy = @ey - @sy
+        @lensq = @dx * @dx + @dy * @dy
+      end
+      
+      
+      attr_reader :s
+      attr_reader :e
+      attr_reader :dx
+      attr_reader :dy
+      
+      
+      def to_s
+        "#{@s} - #{@e}"
+      end
+      
+      
+      def eql?(rhs_)
+        rhs_.kind_of?(Segment) && @s == rhs_.s && @e == rhs_.e
+      end
+      alias_method :==, :eql?
+      
+      
+      def degenerate?
+        @lensq == 0
+      end
+      
+      
+      # Returns a negative value if the point is to the left,
+      # a positive value if the point is to the right, or
+      # 0 if the point is collinear to the segment.
+      
+      def side(p_)
+        px_ = p_.x
+        py_ = p_.y
+        (@sx - px_) * (@ey - py_) - (@sy - py_) * (@ex - px_)
+      end
+      
+      
+      def tproj(p_)
+        if @lensq == 0
+          nil
+        else
+          (@dx * (p_.x - @sx) + @dy * (p_.y - @sy)) / @lensq
+        end
+      end
+      
+      
+      def contains_point?(p_)
+        if side(p_) == 0
+          t_ = tproj(p_)
+          t_ && t_ >= 0.0 && t_ <= 1.0
+        else
+          false
+        end
+      end
+      
+      
+      def intersects_segment?(seg_)
+        s2_ = seg_.s
+        # Handle degenerate cases
+        if seg_.degenerate?
+          if @lensq == 0
+            return @s == s2_
+          else
+            return contains_point?(s2_)
+          end
+        elsif @lensq == 0
+          return seg_.contains_point?(@s)
+        end
+        # Both segments have nonzero length.
+        sx2_ = s2_.x
+        sy2_ = s2_.y
+        dx2_ = seg_.dx
+        dy2_ = seg_.dy
+        denom_ = @dx*dy2_ - @dy*dx2_
+        if denom_ == 0
+          # Segments are parallel. Make sure they are collinear.
+          return false unless side(s2_) == 0
+          # 1-D check.
+          ts_ = (@dx * (sx2_ - @sx) + @dy * (sy2_ - @sy)) / @lensq
+          te_ = (@dx * (sx2_ + dx2_ - @sx) + @dy * (sy2_ + dy2_ - @sy)) / @lensq
+          if ts_ < te_
+            te_ >= 0.0 && ts_ <= 1.0
+          else
+            ts_ >= 0.0 && te_ <= 1.0
+          end
+        else
+          # Segments are not parallel. Check the intersection of their
+          # containing lines.
+          t_ = (dy2_ * (sx2_ - @sx) + dx2_ * (@sy - sy2_)) / denom_
+          return false if t_ < 0.0 || t_ > 1.0
+          t2_ = (@dy * (sx2_ - @sx) + @dx * (@sy - sy2_)) / denom_
+          t2_ >= 0.0 && t2_ <= 1.0
+        end
+      end
+      
+      
+      def length
+        ::Math.sqrt(@lensq)
+      end
+      
+      
+    end
+    
+    
+  end
+  
+end

data/lib/rgeo/cartesian/factory.rb

@@ -0,0 +1,209 @@
+# -----------------------------------------------------------------------------
+# 
+# Geographic data 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 Cartesian
+    
+    
+    # This class implements the factory for the simple cartesian
+    # implementation.
+    
+    class Factory
+      
+      include Feature::Factory::Instance
+      
+      
+      # Create a new simple cartesian factory.
+      # 
+      # See ::RGeo::Cartesian::simple_factory for a list of supported options.
+      
+      def initialize(opts_={})
+        @has_z = opts_[:has_z_coordinate] ? true : false
+        @has_m = opts_[:has_m_coordinate] ? true : false
+        @proj4 = opts_[:proj4]
+        if CoordSys::Proj4.supported?
+          if @proj4.kind_of?(::String) || @proj4.kind_of?(::Hash)
+            @proj4 = CoordSys::Proj4.create(@proj4)
+          end
+        else
+          @proj4 = nil
+        end
+        srid_ = opts_[:srid]
+        @coord_sys = opts_[:coord_sys]
+        if @coord_sys.kind_of?(::String)
+          @coord_sys = CoordSys::CS.create_from_wkt(@coord_sys) rescue nil
+        end
+        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
+        @srid = srid_.to_i
+      end
+      
+      
+      # Equivalence test.
+      
+      def eql?(rhs_)
+        rhs_.is_a?(self.class) && @srid == rhs_.srid && @has_z == rhs_.property(:has_z_coordinate) && @has_m == rhs_.property(:has_m_coordinate)
+      end
+      alias_method :==, :eql?
+      
+      
+      # Returns the SRID.
+      
+      def srid
+        @srid
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#property
+      
+      def property(name_)
+        case name_
+        when :has_z_coordinate
+          @has_z
+        when :has_m_coordinate
+          @has_m
+        when :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_, *extra_)
+        PointImpl.new(self, x_, y_, *extra_) rescue nil
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#line_string
+      
+      def line_string(points_)
+        LineStringImpl.new(self, points_) rescue nil
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#line
+      
+      def line(start_, end_)
+        LineImpl.new(self, start_, end_) rescue nil
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#linear_ring
+      
+      def linear_ring(points_)
+        LinearRingImpl.new(self, points_) rescue nil
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#polygon
+      
+      def polygon(outer_ring_, inner_rings_=nil)
+        PolygonImpl.new(self, outer_ring_, inner_rings_) rescue nil
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#collection
+      
+      def collection(elems_)
+        GeometryCollectionImpl.new(self, elems_) rescue nil
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#multi_point
+      
+      def multi_point(elems_)
+        MultiPointImpl.new(self, elems_) rescue nil
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#multi_line_string
+      
+      def multi_line_string(elems_)
+        MultiLineStringImpl.new(self, elems_) rescue nil
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#multi_polygon
+      
+      def multi_polygon(elems_)
+        MultiPolygonImpl.new(self, elems_) rescue nil
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#proj4
+      
+      def proj4
+        @proj4
+      end
+      
+      
+      # See ::RGeo::Feature::Factory#coord_sys
+      
+      def coord_sys
+        @coord_sys
+      end
+      
+      
+    end
+    
+    
+  end
+  
+end

data/lib/rgeo/cartesian/feature_classes.rb

@@ -0,0 +1,173 @@
+# -----------------------------------------------------------------------------
+# 
+# Spherical geography feature classes
+# 
+# -----------------------------------------------------------------------------
+# 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 Cartesian
+    
+    
+    class PointImpl  # :nodoc:
+      
+      
+      include ::RGeo::Feature::Point
+      include ::RGeo::ImplHelper::BasicGeometryMethods
+      include ::RGeo::ImplHelper::BasicPointMethods
+      include ::RGeo::Cartesian::GeometryMethods
+      
+      
+      def distance(rhs_)
+        rhs_ = ::RGeo::Feature.cast(rhs_, @factory)
+        case rhs_
+        when PointImpl
+          dx_ = @x - rhs_.x
+          dy_ = @y - rhs_.y
+          ::Math.sqrt(dx_ * dx_ + dy_ * dy_)
+        else
+          super
+        end
+      end
+      
+      
+    end
+    
+    
+    class LineStringImpl  # :nodoc:
+      
+      
+      include ::RGeo::Feature::LineString
+      include ::RGeo::ImplHelper::BasicGeometryMethods
+      include ::RGeo::ImplHelper::BasicLineStringMethods
+      include ::RGeo::Cartesian::GeometryMethods
+      include ::RGeo::Cartesian::LineStringMethods
+      
+      
+    end
+    
+    
+    class LineImpl  # :nodoc:
+      
+      
+      include ::RGeo::Feature::Line
+      include ::RGeo::ImplHelper::BasicGeometryMethods
+      include ::RGeo::ImplHelper::BasicLineStringMethods
+      include ::RGeo::ImplHelper::BasicLineMethods
+      include ::RGeo::Cartesian::GeometryMethods
+      include ::RGeo::Cartesian::LineStringMethods
+      
+      
+    end
+    
+    
+    class LinearRingImpl  # :nodoc:
+      
+      
+      include ::RGeo::Feature::Line
+      include ::RGeo::ImplHelper::BasicGeometryMethods
+      include ::RGeo::ImplHelper::BasicLineStringMethods
+      include ::RGeo::ImplHelper::BasicLinearRingMethods
+      include ::RGeo::Cartesian::GeometryMethods
+      include ::RGeo::Cartesian::LineStringMethods
+      
+      
+    end
+    
+    
+    class PolygonImpl  # :nodoc:
+      
+      
+      include ::RGeo::Feature::Polygon
+      include ::RGeo::ImplHelper::BasicGeometryMethods
+      include ::RGeo::ImplHelper::BasicPolygonMethods
+      include ::RGeo::Cartesian::GeometryMethods
+      
+      
+    end
+    
+    
+    class GeometryCollectionImpl  # :nodoc:
+      
+      
+      include ::RGeo::Feature::GeometryCollection
+      include ::RGeo::ImplHelper::BasicGeometryMethods
+      include ::RGeo::ImplHelper::BasicGeometryCollectionMethods
+      include ::RGeo::Cartesian::GeometryMethods
+      
+      
+    end
+    
+    
+    class MultiPointImpl  # :nodoc:
+      
+      
+      include ::RGeo::Feature::GeometryCollection
+      include ::RGeo::ImplHelper::BasicGeometryMethods
+      include ::RGeo::ImplHelper::BasicGeometryCollectionMethods
+      include ::RGeo::ImplHelper::BasicMultiPointMethods
+      include ::RGeo::Cartesian::GeometryMethods
+      
+      
+    end
+    
+    
+    class MultiLineStringImpl  # :nodoc:
+      
+      
+      include ::RGeo::Feature::GeometryCollection
+      include ::RGeo::ImplHelper::BasicGeometryMethods
+      include ::RGeo::ImplHelper::BasicGeometryCollectionMethods
+      include ::RGeo::ImplHelper::BasicMultiLineStringMethods
+      include ::RGeo::Cartesian::GeometryMethods
+      
+      
+    end
+    
+    
+    class MultiPolygonImpl  # :nodoc:
+      
+      
+      include ::RGeo::Feature::GeometryCollection
+      include ::RGeo::ImplHelper::BasicGeometryMethods
+      include ::RGeo::ImplHelper::BasicGeometryCollectionMethods
+      include ::RGeo::ImplHelper::BasicMultiPolygonMethods
+      include ::RGeo::Cartesian::GeometryMethods
+      
+      
+    end
+    
+    
+  end
+  
+end