rgeo 0.5.1 → 0.5.2

data/ext/geos_c_impl/geometry_collection.h

@@ -0,0 +1,46 @@
+/*
+  Geometry collection methods for GEOS wrapper
+*/
+
+
+#ifndef RGEO_GEOS_GEOMETRY_COLLECTION_INCLUDED
+#define RGEO_GEOS_GEOMETRY_COLLECTION_INCLUDED
+
+#include <ruby.h>
+#include <geos_c.h>
+
+#include "factory.h"
+
+RGEO_BEGIN_C
+
+
+/*
+  Initializes the geometry collection module. This should be called after
+  the geometry module is initialized.
+*/
+void rgeo_init_geos_geometry_collection(RGeo_Globals* globals);
+
+/*
+  Comopares the contents of two geometry collections. Does not test the
+  types of the collections themselves, but tests the types, values, and
+  contents of all the contents. The two given geometries MUST be
+  collection types-- i.e. GeometryCollection, MultiPoint, MultiLineString,
+  or MultiPolygon.
+  Returns Qtrue if the contents of the two geometry collections are equal,
+  Qfalse if they are inequal, or Qnil if an error occurs.
+*/
+VALUE rgeo_geos_geometry_collections_eql(GEOSContextHandle_t context, const GEOSGeometry* geom1, const GEOSGeometry* geom2, char check_z);
+
+/*
+  A tool for building up hash values.
+  You must pass in the context, a geos geometry, and a seed hash.
+  Returns an updated hash.
+  This call is useful in sequence, and should be bracketed by calls to
+  rb_hash_start and rb_hash_end.
+*/
+st_index_t rgeo_geos_geometry_collection_hash(GEOSContextHandle_t context, const GEOSGeometry* geom, st_index_t hash);
+
+
+RGEO_END_C
+
+#endif

data/ext/geos_c_impl/line_string.c

@@ -0,0 +1,675 @@
+/*
+  Line string methods for GEOS wrapper
+*/
+
+
+#include "preface.h"
+
+#ifdef RGEO_GEOS_SUPPORTED
+
+#include <string.h>
+#include <ruby.h>
+#include <geos_c.h>
+
+#include "factory.h"
+#include "geometry.h"
+#include "point.h"
+#include "line_string.h"
+
+#include "coordinates.h"
+
+RGEO_BEGIN_C
+
+
+static VALUE method_line_string_geometry_type(VALUE self)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  if (self_data->geom) {
+    result = RGEO_FACTORY_DATA_PTR(self_data->factory)->globals->feature_line_string;
+  }
+  return result;
+}
+
+
+static VALUE method_linear_ring_geometry_type(VALUE self)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  if (self_data->geom) {
+    result = RGEO_FACTORY_DATA_PTR(self_data->factory)->globals->feature_linear_ring;
+  }
+  return result;
+}
+
+
+static VALUE method_line_geometry_type(VALUE self)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  if (self_data->geom) {
+    result = RGEO_FACTORY_DATA_PTR(self_data->factory)->globals->feature_line;
+  }
+  return result;
+}
+
+
+static VALUE method_line_string_length(VALUE self)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+  const GEOSGeometry* self_geom;
+  double len;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  self_geom = self_data->geom;
+  if (self_geom) {
+    if (GEOSLength_r(self_data->geos_context, self_geom, &len)) {
+      result = rb_float_new(len);
+    }
+  }
+  return result;
+}
+
+
+static VALUE method_line_string_num_points(VALUE self)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+  const GEOSGeometry* self_geom;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  self_geom = self_data->geom;
+  if (self_geom) {
+    result = INT2NUM(GEOSGetNumCoordinates_r(self_data->geos_context, self_geom));
+  }
+  return result;
+}
+
+
+static void d(VALUE v) {
+    ID sym_puts = rb_intern("puts");
+    ID sym_inspect = rb_intern("inspect");
+    rb_funcall(rb_mKernel, sym_puts, 1,
+        rb_funcall(v, sym_inspect, 0));
+}
+
+static VALUE method_line_string_coordinates(VALUE self)
+{
+  VALUE result;
+  VALUE point_tuple;
+  RGeo_GeometryData* self_data;
+  const GEOSGeometry* self_geom;
+  const GEOSCoordSequence* coord_sequence;
+  int zCoordinate;
+
+  GEOSContextHandle_t context;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  self_geom = self_data->geom;
+
+  if (self_geom) {
+    zCoordinate = RGEO_FACTORY_DATA_PTR(self_data->factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M;
+    context = self_data->geos_context;
+    coord_sequence = GEOSGeom_getCoordSeq_r(context, self_geom);
+    if(coord_sequence) {
+      result = extract_points_from_coordinate_sequence(context, coord_sequence, zCoordinate);
+    }
+  }
+  return result;
+}
+
+
+
+static VALUE get_point_from_coordseq(VALUE self, const GEOSCoordSequence* coord_seq, unsigned int i, char has_z)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+  GEOSContextHandle_t self_context;
+  double x, y, z;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  self_context = self_data->geos_context;
+  if (GEOSCoordSeq_getX_r(self_context, coord_seq, i, &x)) {
+    if (GEOSCoordSeq_getY_r(self_context, coord_seq, i, &y)) {
+      if (has_z) {
+        if (!GEOSCoordSeq_getZ_r(self_context, coord_seq, i, &z)) {
+          z = 0.0;
+        }
+      }
+      else {
+        z = 0.0;
+      }
+      result = rgeo_create_geos_point(self_data->factory, x, y, z);
+    }
+  }
+  return result;
+}
+
+
+static VALUE method_line_string_point_n(VALUE self, VALUE n)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+  const GEOSGeometry* self_geom;
+  GEOSContextHandle_t self_context;
+  const GEOSCoordSequence* coord_seq;
+  char has_z;
+  int si;
+  unsigned int i;
+  unsigned int size;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  self_geom = self_data->geom;
+  if (self_geom) {
+    self_context = self_data->geos_context;
+    coord_seq = GEOSGeom_getCoordSeq_r(self_context, self_geom);
+    if (coord_seq) {
+      has_z = (char)(RGEO_FACTORY_DATA_PTR(self_data->factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
+      si = NUM2INT(n);
+      if (si >= 0) {
+        i = si;
+        if (GEOSCoordSeq_getSize_r(self_context, coord_seq, &size)) {
+          if (i < size) {
+            result = get_point_from_coordseq(self, coord_seq, i, has_z);
+          }
+        }
+      }
+    }
+  }
+  return result;
+}
+
+
+static VALUE method_line_string_points(VALUE self)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+  const GEOSGeometry* self_geom;
+  GEOSContextHandle_t self_context;
+  const GEOSCoordSequence* coord_seq;
+  char has_z;
+  unsigned int size;
+  double x;
+  double y;
+  double z;
+  unsigned int i;
+  VALUE point;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  self_geom = self_data->geom;
+  if (self_geom) {
+    self_context = self_data->geos_context;
+    coord_seq = GEOSGeom_getCoordSeq_r(self_context, self_geom);
+    if (coord_seq) {
+      has_z = (char)(RGEO_FACTORY_DATA_PTR(self_data->factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
+      if (GEOSCoordSeq_getSize_r(self_context, coord_seq, &size)) {
+        result = rb_ary_new2(size);
+        for (i=0; i<size; ++i) {
+          point = get_point_from_coordseq(self, coord_seq, i, has_z);
+          if (!NIL_P(point)) {
+            rb_ary_store(result, i, point);
+          }
+        }
+      }
+    }
+  }
+  return result;
+}
+
+
+static VALUE method_line_string_start_point(VALUE self)
+{
+  return method_line_string_point_n(self, INT2NUM(0));
+}
+
+
+static VALUE method_line_string_end_point(VALUE self)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+  const GEOSGeometry* self_geom;
+  unsigned int n;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  self_geom = self_data->geom;
+  if (self_geom) {
+    n = GEOSGetNumCoordinates_r(self_data->geos_context, self_geom);
+    if (n > 0) {
+      result = method_line_string_point_n(self, INT2NUM(n-1));
+    }
+  }
+  return result;
+}
+
+
+static VALUE method_line_string_is_closed(VALUE self)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+  const GEOSGeometry* self_geom;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  self_geom = self_data->geom;
+  if (self_geom) {
+    result = rgeo_is_geos_line_string_closed(self_data->geos_context, self_geom);
+  }
+  return result;
+}
+
+
+static VALUE method_line_string_is_ring(VALUE self)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+  const GEOSGeometry* self_geom;
+  char val;
+
+  result = Qnil;
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  self_geom = self_data->geom;
+  if (self_geom) {
+    val = GEOSisRing_r(self_data->geos_context, self_geom);
+    if (val == 0) {
+      result = Qfalse;
+    }
+    else if (val == 1) {
+      result = Qtrue;
+    }
+  }
+  return result;
+}
+
+
+static VALUE method_line_string_eql(VALUE self, VALUE rhs)
+{
+  VALUE result;
+  RGeo_GeometryData* self_data;
+
+  result = rgeo_geos_klasses_and_factories_eql(self, rhs);
+  if (RTEST(result)) {
+    self_data = RGEO_GEOMETRY_DATA_PTR(self);
+    result = rgeo_geos_coordseqs_eql(self_data->geos_context, self_data->geom, RGEO_GEOMETRY_DATA_PTR(rhs)->geom, RGEO_FACTORY_DATA_PTR(self_data->factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
+  }
+  return result;
+}
+
+
+static VALUE method_line_string_hash(VALUE self)
+{
+  st_index_t hash;
+  RGeo_GeometryData* self_data;
+  VALUE factory;
+
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  factory = self_data->factory;
+  hash = rb_hash_start(0);
+  hash = rgeo_geos_objbase_hash(factory,
+    RGEO_FACTORY_DATA_PTR(factory)->globals->feature_line_string, hash);
+  hash = rgeo_geos_coordseq_hash(self_data->geos_context, self_data->geom, hash);
+  return LONG2FIX(rb_hash_end(hash));
+}
+
+
+static VALUE method_linear_ring_hash(VALUE self)
+{
+  st_index_t hash;
+  RGeo_GeometryData* self_data;
+  VALUE factory;
+
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  factory = self_data->factory;
+  hash = rb_hash_start(0);
+  hash = rgeo_geos_objbase_hash(factory,
+    RGEO_FACTORY_DATA_PTR(factory)->globals->feature_linear_ring, hash);
+  hash = rgeo_geos_coordseq_hash(self_data->geos_context, self_data->geom, hash);
+  return LONG2FIX(rb_hash_end(hash));
+}
+
+
+static VALUE method_line_hash(VALUE self)
+{
+  st_index_t hash;
+  RGeo_GeometryData* self_data;
+  VALUE factory;
+
+  self_data = RGEO_GEOMETRY_DATA_PTR(self);
+  factory = self_data->factory;
+  hash = rb_hash_start(0);
+  hash = rgeo_geos_objbase_hash(factory,
+    RGEO_FACTORY_DATA_PTR(factory)->globals->feature_line, hash);
+  hash = rgeo_geos_coordseq_hash(self_data->geos_context, self_data->geom, hash);
+  return LONG2FIX(rb_hash_end(hash));
+}
+
+
+static GEOSCoordSequence* coord_seq_from_array(VALUE factory, VALUE array, char close)
+{
+  RGeo_FactoryData* factory_data;
+  VALUE point_type;
+  unsigned int len;
+  char has_z;
+  unsigned int dims;
+  double* coords;
+  GEOSContextHandle_t context;
+  unsigned int i;
+  char good;
+  const GEOSGeometry* entry_geom;
+  const GEOSCoordSequence* entry_cs;
+  double x;
+  GEOSCoordSequence* coord_seq;
+
+  Check_Type(array, T_ARRAY);
+  factory_data = RGEO_FACTORY_DATA_PTR(factory);
+  point_type = factory_data->globals->feature_point;
+  len = (unsigned int)RARRAY_LEN(array);
+  has_z = (char)(RGEO_FACTORY_DATA_PTR(factory)->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
+  dims = has_z ? 3 : 2;
+  coords = ALLOC_N(double, len == 0 ? 1 : len * dims);
+  if (!coords) {
+    return NULL;
+  }
+  context = factory_data->geos_context;
+  for (i=0; i<len; ++i) {
+    good = 0;
+    entry_geom = rgeo_convert_to_geos_geometry(factory, rb_ary_entry(array, i), point_type);
+    if (entry_geom) {
+      entry_cs = GEOSGeom_getCoordSeq_r(context, entry_geom);
+      if (entry_cs) {
+        if (GEOSCoordSeq_getX_r(context, entry_cs, 0, &x)) {
+          coords[i*dims] = x;
+          if (GEOSCoordSeq_getY_r(context, entry_cs, 0, &x)) {
+            coords[i*dims+1] = x;
+            good = 1;
+            if (has_z) {
+              if (GEOSCoordSeq_getZ_r(context, entry_cs, 0, &x)) {
+                coords[i*dims+2] = x;
+              }
+              else {
+                good = 0;
+              }
+            }
+          }
+        }
+      }
+    }
+    if (!good) {
+      free(coords);
+      return NULL;
+    }
+  }
+  if (len > 0 && close) {
+    if (coords[0] == coords[(len-1)*dims] && coords[1] == coords[(len-1)*dims+1]) {
+      close = 0;
+    }
+  }
+  else {
+    close = 0;
+  }
+  coord_seq = GEOSCoordSeq_create_r(context, len + close, 3);
+  if (coord_seq) {
+    for (i=0; i<len; ++i) {
+      GEOSCoordSeq_setX_r(context, coord_seq, i, coords[i*dims]);
+      GEOSCoordSeq_setY_r(context, coord_seq, i, coords[i*dims+1]);
+      GEOSCoordSeq_setZ_r(context, coord_seq, i, has_z ? coords[i*dims+2] : 0);
+    }
+    if (close) {
+      GEOSCoordSeq_setX_r(context, coord_seq, len, coords[0]);
+      GEOSCoordSeq_setY_r(context, coord_seq, len, coords[1]);
+      GEOSCoordSeq_setZ_r(context, coord_seq, len, has_z ? coords[2] : 0);
+    }
+  }
+  free(coords);
+  return coord_seq;
+}
+
+
+static VALUE cmethod_create_line_string(VALUE module, VALUE factory, VALUE array)
+{
+  VALUE result;
+  GEOSCoordSequence* coord_seq;
+  RGeo_FactoryData* factory_data;
+  GEOSGeometry* geom;
+
+  result = Qnil;
+  coord_seq = coord_seq_from_array(factory, array, 0);
+  if (coord_seq) {
+    factory_data = RGEO_FACTORY_DATA_PTR(factory);
+    geom = GEOSGeom_createLineString_r(factory_data->geos_context, coord_seq);
+    if (geom) {
+      result = rgeo_wrap_geos_geometry(factory, geom, factory_data->globals->geos_line_string);
+    }
+  }
+  return result;
+}
+
+
+static VALUE cmethod_create_linear_ring(VALUE module, VALUE factory, VALUE array)
+{
+  VALUE result;
+  GEOSCoordSequence* coord_seq;
+  RGeo_FactoryData* factory_data;
+  GEOSGeometry* geom;
+
+  result = Qnil;
+  coord_seq = coord_seq_from_array(factory, array, 1);
+  if (coord_seq) {
+    factory_data = RGEO_FACTORY_DATA_PTR(factory);
+    geom = GEOSGeom_createLinearRing_r(factory_data->geos_context, coord_seq);
+    if (geom) {
+      result = rgeo_wrap_geos_geometry(factory, geom, factory_data->globals->geos_linear_ring);
+    }
+  }
+  return result;
+}
+
+
+static void populate_geom_into_coord_seq(GEOSContextHandle_t context, const GEOSGeometry* geom, GEOSCoordSequence* coord_seq, unsigned int i, char has_z)
+{
+  const GEOSCoordSequence* cs;
+  double x;
+
+  cs = GEOSGeom_getCoordSeq_r(context, geom);
+  x = 0;
+  if (cs) {
+    GEOSCoordSeq_getX_r(context, cs, 0, &x);
+  }
+  GEOSCoordSeq_setX_r(context, coord_seq, i, x);
+  x = 0;
+  if (cs) {
+    GEOSCoordSeq_getY_r(context, cs, 0, &x);
+  }
+  GEOSCoordSeq_setY_r(context, coord_seq, i, x);
+  x = 0;
+  if (has_z && cs) {
+    GEOSCoordSeq_getZ_r(context, cs, 0, &x);
+  }
+  GEOSCoordSeq_setZ_r(context, coord_seq, i, x);
+}
+
+
+static VALUE cmethod_create_line(VALUE module, VALUE factory, VALUE start, VALUE end)
+{
+  VALUE result;
+  RGeo_FactoryData* factory_data;
+  char has_z;
+  VALUE point_type;
+  GEOSContextHandle_t context;
+  const GEOSGeometry* start_geom;
+  const GEOSGeometry* end_geom;
+  GEOSCoordSequence* coord_seq;
+  GEOSGeometry* geom;
+
+  result = Qnil;
+  factory_data = RGEO_FACTORY_DATA_PTR(factory);
+  has_z = (char)(factory_data->flags & RGEO_FACTORYFLAGS_SUPPORTS_Z_OR_M);
+  point_type = factory_data->globals->feature_point;
+  context = factory_data->geos_context;
+
+  start_geom = rgeo_convert_to_geos_geometry(factory, start, point_type);
+  if (start_geom) {
+    end_geom = rgeo_convert_to_geos_geometry(factory, end, point_type);
+    if (end_geom) {
+      coord_seq = GEOSCoordSeq_create_r(context, 2, 3);
+      if (coord_seq) {
+        populate_geom_into_coord_seq(context, start_geom, coord_seq, 0, has_z);
+        populate_geom_into_coord_seq(context, end_geom, coord_seq, 1, has_z);
+        geom = GEOSGeom_createLineString_r(context, coord_seq);
+        if (geom) {
+          result = rgeo_wrap_geos_geometry(factory, geom, factory_data->globals->geos_line);
+        }
+      }
+    }
+  }
+
+  return result;
+}
+
+
+static VALUE impl_copy_from(VALUE klass, VALUE factory, VALUE original, char subtype)
+{
+  VALUE result;
+  const GEOSGeometry* original_geom;
+  GEOSContextHandle_t context;
+  const GEOSCoordSequence* original_coord_seq;
+  GEOSCoordSequence* coord_seq;
+  GEOSGeometry* geom;
+
+  result = Qnil;
+  original_geom = RGEO_GEOMETRY_DATA_PTR(original)->geom;
+  if (original_geom) {
+    context = RGEO_FACTORY_DATA_PTR(factory)->geos_context;
+    if (subtype == 1 && GEOSGetNumCoordinates_r(context, original_geom) != 2) {
+      original_geom = NULL;
+    }
+    if (original_geom) {
+      original_coord_seq = GEOSGeom_getCoordSeq_r(context, original_geom);
+      if (original_coord_seq) {
+        coord_seq = GEOSCoordSeq_clone_r(context, original_coord_seq);
+        if (coord_seq) {
+          geom = subtype == 2 ? GEOSGeom_createLinearRing_r(context, coord_seq) : GEOSGeom_createLineString_r(context, coord_seq);
+          if (geom) {
+            result = rgeo_wrap_geos_geometry(factory, geom, klass);
+          }
+        }
+      }
+    }
+  }
+  return result;
+}
+
+
+static VALUE cmethod_line_string_copy_from(VALUE klass, VALUE factory, VALUE original)
+{
+  return impl_copy_from(klass, factory, original, 0);
+}
+
+
+static VALUE cmethod_line_copy_from(VALUE klass, VALUE factory, VALUE original)
+{
+  return impl_copy_from(klass, factory, original, 1);
+}
+
+
+static VALUE cmethod_linear_ring_copy_from(VALUE klass, VALUE factory, VALUE original)
+{
+  return impl_copy_from(klass, factory, original, 2);
+}
+
+
+void rgeo_init_geos_line_string(RGeo_Globals* globals)
+{
+  VALUE geos_line_string_methods;
+  VALUE geos_linear_ring_methods;
+  VALUE geos_line_methods;
+
+  // Class methods for CAPILineStringImpl
+  rb_define_module_function(globals->geos_line_string, "create", cmethod_create_line_string, 2);
+  rb_define_module_function(globals->geos_line_string, "_copy_from", cmethod_line_string_copy_from, 2);
+
+  // Class methods for CAPILinearRingImpl
+  rb_define_module_function(globals->geos_linear_ring, "create", cmethod_create_linear_ring, 2);
+  rb_define_module_function(globals->geos_linear_ring, "_copy_from", cmethod_linear_ring_copy_from, 2);
+
+  // Class methods for CAPILineImpl
+  rb_define_module_function(globals->geos_line, "create", cmethod_create_line, 3);
+  rb_define_module_function(globals->geos_line, "_copy_from", cmethod_line_copy_from, 2);
+
+  // CAPILineStringMethods module
+  geos_line_string_methods = rb_define_module_under(globals->geos_module, "CAPILineStringMethods");
+  rb_define_method(geos_line_string_methods, "rep_equals?", method_line_string_eql, 1);
+  rb_define_method(geos_line_string_methods, "eql?", method_line_string_eql, 1);
+  rb_define_method(geos_line_string_methods, "hash", method_line_string_hash, 0);
+  rb_define_method(geos_line_string_methods, "geometry_type", method_line_string_geometry_type, 0);
+  rb_define_method(geos_line_string_methods, "length", method_line_string_length, 0);
+  rb_define_method(geos_line_string_methods, "num_points", method_line_string_num_points, 0);
+  rb_define_method(geos_line_string_methods, "point_n", method_line_string_point_n, 1);
+  rb_define_method(geos_line_string_methods, "points", method_line_string_points, 0);
+  rb_define_method(geos_line_string_methods, "start_point", method_line_string_start_point, 0);
+  rb_define_method(geos_line_string_methods, "end_point", method_line_string_end_point, 0);
+  rb_define_method(geos_line_string_methods, "is_closed?", method_line_string_is_closed, 0);
+  rb_define_method(geos_line_string_methods, "is_ring?", method_line_string_is_ring, 0);
+  rb_define_method(geos_line_string_methods, "coordinates", method_line_string_coordinates, 0);
+
+  // CAPILinearRingMethods module
+  geos_linear_ring_methods = rb_define_module_under(globals->geos_module, "CAPILinearRingMethods");
+  rb_define_method(geos_linear_ring_methods, "geometry_type", method_linear_ring_geometry_type, 0);
+  rb_define_method(geos_linear_ring_methods, "hash", method_linear_ring_hash, 0);
+
+  // CAPILineMethods module
+  geos_line_methods = rb_define_module_under(globals->geos_module, "CAPILineMethods");
+  rb_define_method(geos_line_methods, "geometry_type", method_line_geometry_type, 0);
+  rb_define_method(geos_line_methods, "hash", method_line_hash, 0);
+}
+
+
+VALUE rgeo_is_geos_line_string_closed(GEOSContextHandle_t context, const GEOSGeometry* geom)
+{
+  VALUE result;
+  unsigned int n;
+  double x1, x2, y1, y2, z1, z2;
+  const GEOSCoordSequence* coord_seq;
+
+  result = Qnil;
+  n = GEOSGetNumCoordinates_r(context, geom);
+  if (n > 0) {
+    coord_seq = GEOSGeom_getCoordSeq_r(context, geom);
+    if (GEOSCoordSeq_getX_r(context, coord_seq, 0, &x1)) {
+      if (GEOSCoordSeq_getX_r(context, coord_seq, n-1, &x2)) {
+        if (x1 == x2) {
+          if (GEOSCoordSeq_getY_r(context, coord_seq, 0, &y1)) {
+            if (GEOSCoordSeq_getY_r(context, coord_seq, n-1, &y2)) {
+              result = y1 == y2 ? Qtrue : Qfalse;
+            }
+          }
+        }
+        else {
+          result = Qfalse;
+        }
+      }
+    }
+  }
+  return result;
+}
+
+
+RGEO_END_C
+
+#endif