rgeo 2.3.0 → 3.0.0.pre.rc.1

data/lib/rgeo/cartesian/sweepline_intersector.rb

@@ -0,0 +1,149 @@
+# frozen_string_literal: true
+
+# -----------------------------------------------------------------------------
+#
+# Simple Sweepline Intersector Class
+#
+# -----------------------------------------------------------------------------
+
+module RGeo
+  module Cartesian
+    # Implements a Sweepline intersector to find all intersections
+    # in a group of segments. The idea is to use a horizontal line starting
+    # at y = +Infinity that sweeps down to y = -Infinity and every time it hits
+    # a new line, it will check if it intersects with any of the segments
+    # the line currently intersects at that y value.
+    # This is a more simplistic implementation that uses an array to hold
+    # observed segments instead of a sorted BST, so performance may be significantly
+    # worse in the case of lots of segments overlapping in y-ranges.
+    class SweeplineIntersector
+      Event = Struct.new(:point, :segment, :is_start)
+      Intersection = Struct.new(:point, :s1, :s2)
+
+      def initialize(segments)
+        @segments = segments
+      end
+      attr_reader :segments
+
+      # Returns the "proper" intersections from the list of segments.
+      #
+      # This will only return intersections that are not the start/end or
+      # end/start of the 2 segments. This could be useful for finding intersections
+      # in a ring for example, because knowing that segments are connected in a linestring
+      # is not always helpful, but those are reported by default.
+      #
+      # Note: This is not the true definition of a proper intersection. A
+      # truly proper intersection does not include colinear intersections and
+      # the intersection must lie in the interior of both segments.
+      #
+      # @return [Array<RGeo::Cartesian::SweeplineIntersector::Intersection>]
+      def proper_intersections
+        return @proper_intersections if defined?(@proper_intersections)
+
+        @proper_intersections = []
+        intersections.each do |intersection|
+          s1 = intersection.s1
+          s2 = intersection.s2
+          pt = intersection.point
+
+          unless (pt == s1.s && pt == s2.e) || (pt == s1.e && pt == s2.s)
+            @proper_intersections << intersection
+          end
+        end
+        @proper_intersections
+      end
+
+      # Computes the intersections of the input segments.
+      #
+      # Creates an event queue from the +events+ and adds segments to the
+      # +observed_segments+ array while their ending event has not been popped
+      # from the queue.
+      #
+      # Compares the new segment from the +is_start+ event to each observed segment
+      # then adds it to +observed_segments+. Records any intersections in to the
+      # returned array.
+      #
+      # @return [Array<RGeo::Cartesian::SweeplineIntersector::Intersection>]
+      def intersections
+        return @intersections if defined?(@intersections)
+
+        @intersections = []
+        observed_segments = Set.new
+        events.each do |e|
+          seg = e.segment
+
+          if e.is_start
+            observed_segments.each do |oseg|
+              int_pt = seg.segment_intersection(oseg)
+              if int_pt
+                intersect = Intersection.new(int_pt, seg, oseg)
+                @intersections << intersect
+              end
+            end
+            observed_segments << seg
+          else
+            observed_segments.delete(seg)
+          end
+        end
+        @intersections
+      end
+
+      # Returns an ordered array of events from the input segments. Events
+      # are the start and endpoints of each segment with an is_start tag to
+      # indicate if this is the starting or ending event for that segment.
+      #
+      # Ordering is done by greatest-y -> smallest-x -> is_start = true.
+      #
+      # @return [Array]
+      def events
+        return @events if defined?(@events)
+
+        @events = []
+        segments.each do |segment|
+          event_pair = create_event_pair(segment)
+          @events.concat(event_pair)
+        end
+
+        @events.sort! do |a, b|
+          if a.point == b.point
+            if a.is_start
+              -1
+            else
+              1
+            end
+          elsif a.point.y == b.point.y
+            a.point.x <=> b.point.x
+          else
+            b.point.y <=> a.point.y
+          end
+        end
+        @events
+      end
+
+      private
+
+      # Creates a pair of events from a segment
+      #
+      # @param segment [Segment]
+      #
+      # @return [Array]
+      def create_event_pair(segment)
+        s = segment.s
+        e = segment.e
+
+        s_event = Event.new(s, segment)
+        e_event = Event.new(e, segment)
+
+        if s.y > e.y || (s.y == e.y && s.x < e.x)
+          s_event.is_start = true
+          e_event.is_start = false
+        else
+          s_event.is_start = false
+          e_event.is_start = true
+        end
+
+        [s_event, e_event]
+      end
+    end
+  end
+end

data/lib/rgeo/cartesian/valid_op.rb

@@ -0,0 +1,71 @@
+# frozen_string_literal: true
+
+module RGeo
+  module Cartesian
+    module ValidOp
+      include ImplHelper::ValidOp
+
+      def validity_helper
+        ValidOpHelpers
+      end
+    end
+
+    module ValidOpHelpers
+      include ImplHelper::ValidOpHelpers
+
+      module_function(*ImplHelper::ValidOpHelpers.singleton_methods) # rubocop:disable Style/AccessModifierDeclarations
+
+      module_function
+
+      # Checks that there are no invalid intersections between the components
+      # of a polygon.
+      #
+      # @param [RGeo::Feature::Polygon] poly
+      #
+      # @return [String] invalid_reason
+      def check_consistent_area(poly)
+        # Get set of unique coords
+        pts = poly.exterior_ring.coordinates.to_set
+        poly.interior_rings.each do |ring|
+          pts += ring.coordinates
+        end
+        num_points = pts.size
+
+        # if additional nodes were added, there must be an intersection
+        # through a boundary.
+        if poly.send(:graph).incident_edges.size > num_points
+          return Error::SELF_INTERSECTION
+        end
+
+        rings = [poly.exterior_ring] + poly.interior_rings
+        return Error::SELF_INTERSECTION if rings.uniq.size != rings.size
+
+        nil
+      end
+
+      # Checks that the interior of a polygon is connected.
+      #
+      # Process to do this is to walk around an interior cycle of the
+      # exterior shell in the polygon's geometry graph. It will keep track
+      # of all the nodes it visited and if that set is a superset of the
+      # coordinates in the exterior_ring, the interior is connected, otherwise
+      # it is split somewhere.
+      #
+      # @param [RGeo::Feature::Polygon] poly
+      #
+      # @return [String] invalid_reason
+      def check_connected_interiors(poly)
+        exterior_coords = poly.exterior_ring.coordinates.to_set
+
+        visited = Set.new
+        poly.send(:graph).geom_edges.first.exterior_edge.and_connected do |hedge|
+          visited << hedge.origin.coordinates
+        end
+
+        return Error::DISCONNECTED_INTERIOR unless exterior_coords.subset?(visited)
+
+        nil
+      end
+    end
+  end
+end

data/lib/rgeo/cartesian.rb

@@ -8,8 +8,11 @@
 
 require_relative "cartesian/calculations"
 require_relative "cartesian/feature_methods"
+require_relative "cartesian/valid_op"
 require_relative "cartesian/feature_classes"
 require_relative "cartesian/factory"
 require_relative "cartesian/interface"
 require_relative "cartesian/bounding_box"
 require_relative "cartesian/analysis"
+require_relative "cartesian/sweepline_intersector"
+require_relative "cartesian/planar_graph"

data/lib/rgeo/coord_sys/cs/wkt_parser.rb

@@ -24,7 +24,7 @@ module RGeo
             return value
           end
           unless @cur_token.is_a?(TypeString)
-            raise Error::ParseError("Found token #{@cur_token} when we expected a value")
+            raise Error::ParseError, "Found token #{@cur_token} when we expected a value"
           end
           type = @cur_token
           next_token
@@ -48,7 +48,7 @@ module RGeo
           when "TOWGS84"
             bursa_wolf_params = args.find_all(Numeric)
             unless bursa_wolf_params.size == 7
-              raise Error::ParseError("Expected 7 Bursa Wolf parameters but found #{bursa_wolf_params.size}")
+              raise Error::ParseError, "Expected 7 Bursa Wolf parameters but found #{bursa_wolf_params.size}"
             end
             obj = WGS84ConversionInfo.create(*bursa_wolf_params)
           when "UNIT"
@@ -87,7 +87,7 @@ module RGeo
             unit = args.find_first(Unit)
             axes = args.find_all(AxisInfo)
             unless axes.size > 0
-              raise Error::ParseError("Expected at least one AXIS in a LOCAL_CS")
+              raise Error::ParseError, "Expected at least one AXIS in a LOCAL_CS"
             end
             obj = LocalCoordinateSystem.create(name, local_datum, unit, axes, *args.create_optionals)
           when "GEOCCS"
@@ -97,7 +97,7 @@ module RGeo
             linear_unit = args.find_first(LinearUnit)
             axes = args.find_all(AxisInfo)
             unless axes.size == 0 || axes.size == 3
-              raise Error::ParseError("GEOCCS must contain either 0 or 3 AXIS parameters")
+              raise Error::ParseError, "GEOCCS must contain either 0 or 3 AXIS parameters"
             end
             obj = GeocentricCoordinateSystem.create(name, horizontal_datum, prime_meridian, linear_unit, axes[0], axes[1], axes[2], *args.create_optionals)
           when "VERT_CS"
@@ -113,7 +113,7 @@ module RGeo
             angular_unit = args.find_first(AngularUnit)
             axes = args.find_all(AxisInfo)
             unless axes.size == 0 || axes.size == 2
-              raise Error::ParseError("GEOGCS must contain either 0 or 2 AXIS parameters")
+              raise Error::ParseError, "GEOGCS must contain either 0 or 2 AXIS parameters"
             end
             obj = GeographicCoordinateSystem.create(name, angular_unit, horizontal_datum, prime_meridian, axes[0], axes[1], *args.create_optionals)
           when "PROJCS"
@@ -125,7 +125,7 @@ module RGeo
             linear_unit = args.find_first(LinearUnit)
             axes = args.find_all(AxisInfo)
             unless axes.size == 0 || axes.size == 2
-              raise Error::ParseError("PROJCS must contain either 0 or 2 AXIS parameters")
+              raise Error::ParseError, "PROJCS must contain either 0 or 2 AXIS parameters"
             end
             obj = ProjectedCoordinateSystem.create(name, geographic_coordinate_system, projection, linear_unit, axes[0], axes[1], *args.create_optionals)
           else

data/lib/rgeo/error.rb

@@ -29,5 +29,20 @@ module RGeo
     # Parsing failed
     class ParseError < RGeoError
     end
+
+    # Standard error messages from
+    # https://github.com/locationtech/jts/blob/0afbfb1956ec24912a8b4dc4edff0f1200442857/modules/core/src/main/java/org/locationtech/jts/operation/valid/TopologyValidationError.java#L98-L110
+    TOPOLOGY_VALIDATION_ERR = "Topology Validation Error"
+    REPEATED_POINT = "Repeated Point"
+    HOLE_OUTSIDE_SHELL = "Hole lies outside shell"
+    NESTED_HOLES = "Holes are nested"
+    DISCONNECTED_INTERIOR = "Interior is disconnected"
+    SELF_INTERSECTION = "Self-intersection"
+    RING_SELF_INTERSECTION = "Ring Self-intersection"
+    NESTED_SHELLS = "Nested shells"
+    DUPLICATE_RINGS = "Duplicate Rings"
+    TOO_FEW_POINTS = "Too few distinct points in geometry component"
+    INVALID_COORDINATE = "Invalid Coordinate"
+    UNCLOSED_RING = "Ring is not closed"
   end
 end

data/lib/rgeo/feature/curve.rb

@@ -90,8 +90,13 @@ module RGeo
       # Returns a boolean value. Note that this is different from the SFS
       # specification, which stipulates an integer return value.
 
+      def closed?
+        raise Error::UnsupportedOperation, "Method Curve#closed? not defined."
+      end
+
       def is_closed?
-        raise Error::UnsupportedOperation, "Method Curve#is_closed? not defined."
+        warn "The is_closed? method is deprecated, please use the closed? counterpart, will be removed in v3" unless ENV["RGEO_SILENCE_DEPRECATION"]
+        closed?
       end
 
       # === SFS 1.1 Description
@@ -105,8 +110,13 @@ module RGeo
       # Returns a boolean value. Note that this is different from the SFS
       # specification, which stipulates an integer return value.
 
+      def ring?
+        raise Error::UnsupportedOperation, "Method Curve#ring? not defined."
+      end
+
       def is_ring?
-        raise Error::UnsupportedOperation, "Method Curve#is_ring? not defined."
+        warn "The is_ring? method is deprecated, please use the ring? counterpart, will be removed in v3" unless ENV["RGEO_SILENCE_DEPRECATION"]
+        ring?
       end
     end
   end

data/lib/rgeo/feature/geometry.rb

@@ -90,7 +90,7 @@ module RGeo
       # operations on them.)
 
       def factory
-        raise Error::UnsupportedOperation, "Method Geometry#factory not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#factory not defined."
       end
 
       # === SFS 1.1 Description
@@ -105,7 +105,7 @@ module RGeo
       # point geometries, 1 for curves, and 2 for surfaces.
 
       def dimension
-        raise Error::UnsupportedOperation, "Method Geometry#dimension not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#dimension not defined."
       end
 
       # === SFS 1.1 Description
@@ -122,7 +122,7 @@ module RGeo
       # call the +type_name+ method of the returned module.
 
       def geometry_type
-        raise Error::UnsupportedOperation, "Method Geometry#geometry_type not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#geometry_type not defined."
       end
 
       # === SFS 1.1 Description
@@ -137,7 +137,7 @@ module RGeo
       # stored in either the same or some other datastore.
 
       def srid
-        raise Error::UnsupportedOperation, "Method Geometry#srid not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#srid not defined."
       end
 
       # === SFS 1.1 Description
@@ -151,7 +151,7 @@ module RGeo
       # Returns an object that supports the Geometry interface.
 
       def envelope
-        raise Error::UnsupportedOperation, "Method Geometry#envelope not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#envelope not defined."
       end
 
       # === SFS 1.1 Description
@@ -164,7 +164,7 @@ module RGeo
       # Returns an ASCII string.
 
       def as_text
-        raise Error::UnsupportedOperation, "Method Geometry#as_text not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#as_text not defined."
       end
 
       # === SFS 1.1 Description
@@ -177,7 +177,7 @@ module RGeo
       # Returns a binary string.
 
       def as_binary
-        raise Error::UnsupportedOperation, "Method Geometry#as_binary not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#as_binary not defined."
       end
 
       # === SFS 1.1 Description
@@ -191,8 +191,13 @@ module RGeo
       # Returns a boolean value. Note that this is different from the SFS
       # specification, which stipulates an integer return value.
 
+      def empty?
+        raise Error::UnsupportedOperation, "Method #{self.class}#empty? not defined."
+      end
+
       def is_empty?
-        raise Error::UnsupportedOperation, "Method Geometry#is_empty? not defined."
+        warn "The is_empty? method is deprecated, please use the empty? counterpart, will be removed in v3" unless ENV["RGEO_SILENCE_DEPRECATION"]
+        empty?
       end
 
       # === SFS 1.1 Description
@@ -208,8 +213,13 @@ module RGeo
       # Returns a boolean value. Note that this is different from the SFS
       # specification, which stipulates an integer return value.
 
+      def simple?
+        raise Error::UnsupportedOperation, "Method #{self.class}#simple? not defined."
+      end
+
       def is_simple?
-        raise Error::UnsupportedOperation, "Method Geometry#is_simple? not defined."
+        warn "The is_simple? method is deprecated, please use the simple? counterpart, will be removed in v3" unless ENV["RGEO_SILENCE_DEPRECATION"]
+        simple?
       end
 
       # === SFS 1.1 Description
@@ -224,7 +234,7 @@ module RGeo
       # Returns an object that supports the Geometry interface.
 
       def boundary
-        raise Error::UnsupportedOperation, "Method Geometry#boundary not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#boundary not defined."
       end
 
       # === SFS 1.1 Description
@@ -243,7 +253,7 @@ module RGeo
       # from different factories is undefined.
 
       def equals?(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#equals? not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#equals? not defined."
       end
 
       # === SFS 1.1 Description
@@ -262,7 +272,7 @@ module RGeo
       # from different factories is undefined.
 
       def disjoint?(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#disjoint? not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#disjoint? not defined."
       end
 
       # === SFS 1.1 Description
@@ -281,7 +291,7 @@ module RGeo
       # from different factories is undefined.
 
       def intersects?(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#intersects? not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#intersects? not defined."
       end
 
       # === SFS 1.1 Description
@@ -300,7 +310,7 @@ module RGeo
       # from different factories is undefined.
 
       def touches?(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#touches? not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#touches? not defined."
       end
 
       # === SFS 1.1 Description
@@ -319,7 +329,7 @@ module RGeo
       # from different factories is undefined.
 
       def crosses?(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#crosses? not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#crosses? not defined."
       end
 
       # === SFS 1.1 Description
@@ -338,7 +348,7 @@ module RGeo
       # from different factories is undefined.
 
       def within?(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#within? not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#within? not defined."
       end
 
       # === SFS 1.1 Description
@@ -357,7 +367,7 @@ module RGeo
       # from different factories is undefined.
 
       def contains?(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#contains? not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#contains? not defined."
       end
 
       # === SFS 1.1 Description
@@ -376,7 +386,7 @@ module RGeo
       # from different factories is undefined.
 
       def overlaps?(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#overlaps? not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#overlaps? not defined."
       end
 
       # === SFS 1.1 Description
@@ -402,7 +412,7 @@ module RGeo
       # from different factories is undefined.
 
       def relate?(another_geometry, _intersection_pattern_matrix_)
-        raise Error::UnsupportedOperation, "Method Geometry#relate not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#relate not defined."
       end
 
       # === SFS 1.1 Description
@@ -421,7 +431,7 @@ module RGeo
       # distance between objects from different factories is undefined.
 
       def distance(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#distance not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#distance not defined."
       end
 
       # === SFS 1.1 Description
@@ -436,7 +446,7 @@ module RGeo
       # Returns an object that supports the Geometry interface.
 
       def buffer(_distance_)
-        raise Error::UnsupportedOperation, "Method Geometry#buffer not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#buffer not defined."
       end
 
       # === SFS 1.1 Description
@@ -449,7 +459,7 @@ module RGeo
       # Returns an object that supports the Geometry interface.
 
       def convex_hull
-        raise Error::UnsupportedOperation, "Method Geometry#convex_hull not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#convex_hull not defined."
       end
 
       # === SFS 1.1 Description
@@ -467,7 +477,7 @@ module RGeo
       # operations on objects from different factories is undefined.
 
       def intersection(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#intersection not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#intersection not defined."
       end
 
       # === SFS 1.1 Description
@@ -485,7 +495,7 @@ module RGeo
       # operations on objects from different factories is undefined.
 
       def union(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#union not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#union not defined."
       end
 
       # === SFS 1.1 Description
@@ -503,7 +513,7 @@ module RGeo
       # operations on objects from different factories is undefined.
 
       def difference(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#difference not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#difference not defined."
       end
 
       # === SFS 1.1 Description
@@ -521,7 +531,7 @@ module RGeo
       # operations on objects from different factories is undefined.
 
       def sym_difference(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#sym_difference not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#sym_difference not defined."
       end
 
       # Returns true if this geometric object is representationally
@@ -533,7 +543,7 @@ module RGeo
       # from different factories is undefined.
 
       def rep_equals?(another_geometry)
-        raise Error::UnsupportedOperation, "Method Geometry#rep_equals? not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#rep_equals? not defined."
       end
 
       # Unions a collection of Geometry or a single Geometry
@@ -550,7 +560,7 @@ module RGeo
       # returns nil.
       #
       def unary_union
-        raise Error::UnsupportedOperation, "Method Geometry#unary_union not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#unary_union not defined."
       end
 
       # This method should behave almost the same as the rep_equals?

data/lib/rgeo/feature/geometry_collection.rb

@@ -44,7 +44,7 @@ module RGeo
       # Returns an integer.
 
       def num_geometries
-        raise Error::UnsupportedOperation, "Method GeometryCollection#num_geometries not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#num_geometries not defined."
       end
 
       # === SFS 1.1 Description
@@ -59,7 +59,7 @@ module RGeo
       # in that it does not support negative indexes.
 
       def geometry_n(n)
-        raise Error::UnsupportedOperation, "Method GeometryCollection#geometry_n not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#geometry_n not defined."
       end
 
       # Alias of the num_geometries method.
@@ -79,13 +79,13 @@ module RGeo
       # returns nil, where [-1] returns the last element of the collection.
 
       def [](n)
-        raise Error::UnsupportedOperation, "Method GeometryCollection#[] not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#[] not defined."
       end
 
       # Nodes the linework in a list of Geometries
       #
       def node
-        raise Error::UnsupportedOperation, "Method GeometryCollection#node not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#node not defined."
       end
 
       # Iterates over the geometries of this GeometryCollection.
@@ -96,7 +96,15 @@ module RGeo
       # include the Enumerable mixin.
 
       def each(&block)
-        raise Error::UnsupportedOperation, "Method GeometryCollection#each not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#each not defined."
+      end
+
+      # Gives a point that is guaranteed to be within the geometry.
+      #
+      # Extends OGC SFS 1.1 and follows PostGIS standards.
+      # @see https://postgis.net/docs/ST_PointOnSurface.html
+      def point_on_surface
+        raise Error::UnsupportedOperation, "Method #{self.class}#each not defined."
       end
     end
   end

data/lib/rgeo/feature/line_string.rb

@@ -34,7 +34,7 @@ module RGeo
       # Returns an integer.
 
       def num_points
-        raise Error::UnsupportedOperation, "Method LineString#num_points not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#num_points not defined."
       end
 
       # === SFS 1.1 Description
@@ -48,14 +48,14 @@ module RGeo
       # Does not support negative indexes.
 
       def point_n(n)
-        raise Error::UnsupportedOperation, "Method LineString#point_n not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#point_n not defined."
       end
 
       # Returns the constituent points as an array of objects that
       # support the Point interface.
 
       def points
-        raise Error::UnsupportedOperation, "Method LineString#points not defined."
+        raise Error::UnsupportedOperation, "Method #{self.class}#points not defined."
       end
     end
   end