schleyfox-rgeo 0.2.5

data/lib/rgeo/geographic/proj4_projector.rb

@@ -0,0 +1,98 @@
+# -----------------------------------------------------------------------------
+# 
+# Proj4 projection
+# 
+# -----------------------------------------------------------------------------
+# 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 Geographic
+    
+    
+    class Proj4Projector  # :nodoc:
+      
+      
+      def initialize(geography_factory_, projection_factory_)
+        @geography_factory = geography_factory_
+        @projection_factory = projection_factory_
+      end
+      
+      
+      def project(geometry_)
+        Feature.cast(geometry_, @projection_factory, :project)
+      end
+      
+      
+      def unproject(geometry_)
+        Feature.cast(geometry_, @geography_factory, :project)
+      end
+      
+      
+      def projection_factory
+        @projection_factory
+      end
+      
+      
+      def wraps?
+        false
+      end
+      
+      
+      def limits_window
+        nil
+      end
+      
+      
+      class << self
+        
+        
+        def create_from_existing_factory(geography_factory_, projection_factory_)
+          new(geography_factory_, projection_factory_)
+        end
+        
+        
+        def create_from_proj4(geography_factory_, proj4_, opts_={})
+          projection_factory_ = Cartesian.preferred_factory(:proj4 => proj4_, :coord_sys => opts_[:coord_sys], :srid => opts_[:srid], :buffer_resolution => opts_[:buffer_resolution], :lenient_multi_polygon_assertions => opts_[:lenient_multi_polygon_assertions], :has_z_coordinate => opts_[:has_z_coordinate], :has_m_coordinate => opts_[:has_m_coordinate])
+          new(geography_factory_, projection_factory_)
+        end
+        
+        
+      end
+      
+      
+    end
+    
+    
+  end
+  
+end

data/lib/rgeo/geographic/projected_feature_classes.rb

@@ -0,0 +1,213 @@
+# -----------------------------------------------------------------------------
+# 
+# Projtected geographic 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 Geographic
+    
+    
+    class ProjectedPointImpl  # :nodoc:
+      
+      
+      include Feature::Point
+      include ImplHelper::BasicGeometryMethods
+      include ImplHelper::BasicPointMethods
+      include ProjectedGeometryMethods
+      
+      
+      def _validate_geometry
+        @y = 85.0511287 if @y > 85.0511287
+        @y = -85.0511287 if @y < -85.0511287
+        super
+      end
+      
+      
+      def canonical_x
+        x_ = @x % 360.0
+        x_ -= 360.0 if x_ >= 180.0
+        x_
+      end
+      alias_method :canonical_longitude, :canonical_x
+      alias_method :canonical_lon, :canonical_x
+      
+      
+      def canonical_point
+        if @x >= -180.0 && @x < 180.0
+          self
+        else
+          PointImpl.new(@factory, canonical_x, @y)
+        end
+      end
+      
+      
+      alias_method :longitude, :x
+      alias_method :lon, :x
+      alias_method :latitude, :y
+      alias_method :lat, :y
+      
+      
+    end
+    
+    
+    class ProjectedLineStringImpl  # :nodoc:
+      
+      
+      include Feature::LineString
+      include ImplHelper::BasicGeometryMethods
+      include ImplHelper::BasicLineStringMethods
+      include ProjectedGeometryMethods
+      include ProjectedNCurveMethods
+      include ProjectedLineStringMethods
+      
+      
+    end
+    
+    
+    class ProjectedLinearRingImpl  # :nodoc:
+      
+      
+      include Feature::Line
+      include ImplHelper::BasicGeometryMethods
+      include ImplHelper::BasicLineStringMethods
+      include ImplHelper::BasicLinearRingMethods
+      include ProjectedGeometryMethods
+      include ProjectedNCurveMethods
+      include ProjectedLineStringMethods
+      
+      
+    end
+    
+    
+    class ProjectedLineImpl  # :nodoc:
+      
+      
+      include Feature::Line
+      include ImplHelper::BasicGeometryMethods
+      include ImplHelper::BasicLineStringMethods
+      include ImplHelper::BasicLineMethods
+      include ProjectedGeometryMethods
+      include ProjectedNCurveMethods
+      include ProjectedLineStringMethods
+      
+      
+    end
+    
+    
+    class ProjectedPolygonImpl  # :nodoc:
+      
+      
+      include Feature::Polygon
+      include ImplHelper::BasicGeometryMethods
+      include ImplHelper::BasicPolygonMethods
+      include ProjectedGeometryMethods
+      include ProjectedNSurfaceMethods
+      
+      
+      def _validate_geometry
+        super
+        unless projection
+          raise Error::InvalidGeometry, 'Polygon failed assertions'
+        end
+      end
+      
+      
+    end
+    
+    
+    class ProjectedGeometryCollectionImpl  # :nodoc:
+      
+      
+      include Feature::GeometryCollection
+      include ImplHelper::BasicGeometryMethods
+      include ImplHelper::BasicGeometryCollectionMethods
+      include ProjectedGeometryMethods
+      
+      
+    end
+    
+    
+    class ProjectedMultiPointImpl  # :nodoc:
+      
+      
+      include Feature::GeometryCollection
+      include ImplHelper::BasicGeometryMethods
+      include ImplHelper::BasicGeometryCollectionMethods
+      include ImplHelper::BasicMultiPointMethods
+      include ProjectedGeometryMethods
+      
+      
+    end
+    
+    
+    class ProjectedMultiLineStringImpl  # :nodoc:
+      
+      
+      include Feature::GeometryCollection
+      include ImplHelper::BasicGeometryMethods
+      include ImplHelper::BasicGeometryCollectionMethods
+      include ImplHelper::BasicMultiLineStringMethods
+      include ProjectedGeometryMethods
+      include ProjectedNCurveMethods
+      
+      
+    end
+    
+    
+    class ProjectedMultiPolygonImpl  # :nodoc:
+      
+      
+      include Feature::GeometryCollection
+      include ImplHelper::BasicGeometryMethods
+      include ImplHelper::BasicGeometryCollectionMethods
+      include ImplHelper::BasicMultiPolygonMethods
+      include ProjectedGeometryMethods
+      include ProjectedNSurfaceMethods
+      
+      
+      def _validate_geometry
+        super
+        unless projection
+          raise Error::InvalidGeometry, 'MultiPolygon failed assertions'
+        end
+      end
+      
+      
+    end
+    
+    
+  end
+  
+end

data/lib/rgeo/geographic/projected_feature_methods.rb

@@ -0,0 +1,228 @@
+# -----------------------------------------------------------------------------
+# 
+# Projected geographic common method definitions
+# 
+# -----------------------------------------------------------------------------
+# 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 Geographic
+    
+    
+    module ProjectedGeometryMethods  # :nodoc:
+      
+      
+      def srid
+        factory.srid
+      end
+      
+      
+      def projection
+        unless defined?(@projection)
+          @projection = factory.project(self)
+        end
+        @projection
+      end
+      
+      
+      def envelope
+        factory.unproject(projection.envelope)
+      end
+      
+      
+      def is_empty?
+        projection.is_empty?
+      end
+      
+      
+      def is_simple?
+        projection.is_simple?
+      end
+      
+      
+      def boundary
+        factory.unproject(projection.boundary)
+      end
+      
+      
+      def equals?(rhs_)
+        projection.equals?(Feature.cast(rhs_, factory).projection)
+      end
+      
+      
+      def disjoint?(rhs_)
+        projection.disjoint?(Feature.cast(rhs_, factory).projection)
+      end
+      
+      
+      def intersects?(rhs_)
+        projection.intersects?(Feature.cast(rhs_, factory).projection)
+      end
+      
+      
+      def touches?(rhs_)
+        projection.touches?(Feature.cast(rhs_, factory).projection)
+      end
+      
+      
+      def crosses?(rhs_)
+        projection.crosses?(Feature.cast(rhs_, factory).projection)
+      end
+      
+      
+      def within?(rhs_)
+        projection.within?(Feature.cast(rhs_, factory).projection)
+      end
+      
+      
+      def contains?(rhs_)
+        projection.contains?(Feature.cast(rhs_, factory).projection)
+      end
+      
+      
+      def overlaps?(rhs_)
+        projection.overlaps?(Feature.cast(rhs_, factory).projection)
+      end
+      
+      
+      def relate(rhs_, pattern_)
+        projection.relate(Feature.cast(rhs_, factory).projection, pattern_)
+      end
+      
+      
+      def distance(rhs_)
+        projection.distance(Feature.cast(rhs_, factory).projection)
+      end
+      
+      
+      def buffer(distance_)
+        factory.unproject(projection.buffer(distance_))
+      end
+      
+      
+      def convex_hull
+        factory.unproject(projection.convex_hull)
+      end
+      
+      
+      def intersection(rhs_)
+        factory.unproject(projection.intersection(Feature.cast(rhs_, factory).projection))
+      end
+      
+      
+      def union(rhs_)
+        factory.unproject(projection.union(Feature.cast(rhs_, factory).projection))
+      end
+      
+      
+      def difference(rhs_)
+        factory.unproject(projection.difference(Feature.cast(rhs_, factory).projection))
+      end
+      
+      
+      def sym_difference(rhs_)
+        factory.unproject(projection.sym_difference(Feature.cast(rhs_, factory).projection))
+      end
+      
+      
+    end
+    
+    
+    module ProjectedNCurveMethods  # :nodoc:
+      
+      
+      def length
+        projection.length
+      end
+      
+      
+    end
+    
+    
+    module ProjectedLineStringMethods  # :nodoc:
+      
+      
+      def _validate_geometry
+        size_ = @points.size
+        if size_ > 1
+          last_ = @points[0]
+          (1...size_).each do |i_|
+            p_ = @points[i_]
+            last_x_ = last_.x
+            p_x_ = p_.x
+            changed_ = true
+            if p_x_ < last_x_ - 180.0
+              p_x_ += 360.0 while p_x_ < last_x_ - 180.0
+            elsif p_x_ > last_x_ + 180.0
+              p_x_ +- 360.0 while p_x_ > last_x_ + 180.0
+            else
+              changed_ = false
+            end
+            if changed_
+              p_ = factory.point(p_x_, p_.y)
+              @points[i_] = p_
+            end
+            last_ = p_
+          end
+        end
+        super
+      end
+      
+      
+    end
+    
+    
+    module ProjectedNSurfaceMethods  # :nodoc:
+      
+      
+      def area
+        projection.area
+      end
+      
+      
+      def centroid
+        factory.unproject(projection.centroid)
+      end
+      
+      
+      def point_on_surface
+        factory.unproject(projection.point_on_surface)
+      end
+      
+      
+    end
+    
+    
+  end
+  
+end