charta 0.1.12 → 0.1.17

checksums.yaml

@@ -1,7 +1,7 @@
 ---
-SHA1:
-  metadata.gz: c2df82681d22428e8a351eb7021b92633c58734d
-  data.tar.gz: 882e807081a9b46d51ce58e67b823268f98921a7
+SHA256:
+  metadata.gz: 3bdf87fdffb17cf24b38fb76b86d84593bc3f6ca0d0d3e8ffcf5df3d76d3db1f
+  data.tar.gz: bda9a9fd027f0bfdb514a7de6e92429977dc55e5fe7db38c27f6a59a871df875
 SHA512:
-  metadata.gz: d1bda79d73e0533c5f768266ea2543d5af2776d5435efaf5f81dab0f519ae5fb39864a2d20ae3ef0d4b80300e3af5b965f46dc767b057f1db8034bb0e7276722
-  data.tar.gz: 2de5ce38da7e6a7258336fb96a0b05c9a24f7e7d0d8d531a232d82abc279eddd0c111bdaed944120e60c8c1c85fa2b06114e7c91781d06c7509ddabd264677fa
+  metadata.gz: a808aff9f3da978cd293d0084e9c2a61c7a414fe259ca71821ed195afc50f57e99b7d352e0651723dfa7379d70881dfbf01c4d2959e5602c121e94a78bf14b90
+  data.tar.gz: 54a3e3e081602b978bd4b60d9c8bbbf894d22e14b501715f34278ba48c4ee76489b29b92f2f51419ccb311cd046f9871865240565385b6f5848a63a2b257b09c

data/.gitignore

@@ -9,3 +9,4 @@
 /spec/reports/
 /tmp/
 .byebug_history
+.ruby-version

data/.gitlab-ci.yml

@@ -0,0 +1,13 @@
+stages:
+  - test
+
+test:
+  image: registry.gitlab.com/ekylibre/docker-base-images/ruby2.3:1
+  before_script:
+    - bundle install --path vendor/bundle
+  cache:
+    key: bundle
+    paths:
+      - vendor/bundle
+  script:
+    - bundle exec rake test

data/lib/charta.rb

@@ -1,6 +1,10 @@
 # Gathers geomatic calculations
 # Completes RGeo
+require 'bigdecimal'
+require 'bigdecimal/util'
 require 'rgeo'
+require 'charta/coordinates'
+require 'charta/ewkt_serializer'
 require 'charta/geometry'
 require 'charta/geometry_collection'
 require 'charta/point'

data/lib/charta/bounding_box.rb

@@ -24,5 +24,9 @@ module Charta
     def to_a
       [[@y_min, @x_min], [@y_max, @x_max]]
     end
+
+    def to_bbox_string
+      "#{@x_min}, #{@y_min}, #{x_max}, #{y_max}"
+    end
   end
 end

data/lib/charta/coordinates.rb

@@ -0,0 +1,65 @@
+module Charta
+  module Coordinates
+    class << self
+
+      # Force coordinates to 2D
+      def flatten(hash)
+        map_coordinates(hash) { |position| position[0..1] }
+      end
+
+      def map_coordinates(hash, &block)
+        case hash['type']
+          when 'FeatureCollection'
+            map_feature_collection_coordinates hash, &block
+          when 'Feature'
+            map_feature_coordinates hash, &block
+          else
+            map_geometry_coordinates hash, &block
+        end
+      end
+
+      def normalize_4326_geometry(json)
+        map_coordinates json do |(x, y)|
+          [((x + 180.to_d) % 360.to_d) - 180.to_d, ((y + 90.to_d) % 180.to_d) - 90.to_d]
+        end
+      end
+
+      private
+
+        def map_feature_collection_coordinates(hash, &block)
+          hash.merge 'features' => hash['features'].map { |feature| map_feature_coordinates feature, &block }
+        end
+
+        def map_feature_coordinates(hash, &block)
+          hash.merge 'geometry' => map_geometry_coordinates(hash['geometry'], &block)
+        end
+
+        def map_geometry_coordinates(hash, &block)
+          if hash['type'] == 'GeometryCollection'
+            map_geometry_collection_coordinates hash, &block
+          else
+            coordinates = hash['coordinates']
+            mapped =
+              case hash['type']
+                when 'Point' then
+                  block.call coordinates
+                when 'MultiPoint', 'LineString'
+                  coordinates.map(&block)
+                when 'MultiLineString', 'Polygon'
+                  coordinates.map { |line| line.map(&block) }
+                when 'MultiPolygon'
+                  coordinates.map { |poly| poly.map { |line| line.map(&block) } }
+                else
+                  raise StandardError, "Cannot handle: #{hash['type'].inspect}. In #{hash.inspect}"
+              end
+
+            hash.merge 'coordinates' => mapped
+          end
+        end
+
+        def map_geometry_collection_coordinates(hash, &block)
+          hash.merge 'geometries' => hash['geometries'].map { |geometry| map_geometry_coordinates(geometry, &block) }
+        end
+    end
+  end
+end

data/lib/charta/ewkt_serializer.rb

@@ -0,0 +1,92 @@
+module Charta
+  module EwktSerializer
+
+    class << self
+      def object_to_ewkt(hash)
+        type = hash[:type] || hash['type']
+        send("#{type.gsub(/(.)([A-Z])/, '\1_\2').downcase}_to_ewkt", hash)
+      end
+
+      private
+
+        def feature_collection_to_ewkt(hash)
+          return 'GEOMETRYCOLLECTION EMPTY' if hash['features'].nil?
+          'GEOMETRYCOLLECTION(' + hash['features'].collect do |feature|
+            object_to_ewkt(feature)
+          end.join(', ') + ')'
+        end
+
+        def geometry_collection_to_ewkt(hash)
+          return 'GEOMETRYCOLLECTION EMPTY' if hash['geometries'].nil?
+          'GEOMETRYCOLLECTION(' + hash['geometries'].collect do |feature|
+            object_to_ewkt(feature)
+          end.join(', ') + ')'
+        end
+
+        def feature_to_ewkt(hash)
+          object_to_ewkt(hash['geometry'])
+        end
+
+        def point_to_ewkt(hash)
+          return 'POINT EMPTY' if hash['coordinates'].nil?
+          'POINT(' + hash['coordinates'].join(' ') + ')'
+        end
+
+        def line_string_to_ewkt(hash)
+          return 'LINESTRING EMPTY' if hash['coordinates'].nil?
+          'LINESTRING(' + hash['coordinates'].collect do |point|
+            point.join(' ')
+          end.join(', ') + ')'
+        end
+
+        def polygon_to_ewkt(hash)
+          return 'POLYGON EMPTY' if hash['coordinates'].nil?
+          'POLYGON(' + hash['coordinates'].collect do |hole|
+            '(' + hole.collect do |point|
+              point.join(' ')
+            end.join(', ') + ')'
+          end.join(', ') + ')'
+        end
+
+        def multi_point_to_ewkt(hash)
+          return 'MULTIPOINT EMPTY' if hash['coordinates'].nil?
+          'MULTIPOINT(' + hash['coordinates'].collect do |point|
+            '(' + point.join(' ') + ')'
+          end.join(', ') + ')'
+        end
+
+        def multi_line_string_to_ewkt(hash)
+          return 'MULTILINESTRING EMPTY' if hash['coordinates'].nil?
+          'MULTILINESTRING(' + hash['coordinates'].collect do |line|
+            '(' + line.collect do |point|
+              point.join(' ')
+            end.join(', ') + ')'
+          end.join(', ') + ')'
+        end
+
+        def multipolygon_to_ewkt(hash)
+          return 'MULTIPOLYGON EMPTY' if hash['coordinates'].nil?
+          'MULTIPOLYGON(' + hash['coordinates'].collect do |polygon|
+            '(' + polygon.collect do |hole|
+              '(' + hole.collect do |point|
+                point.join(' ')
+              end.join(', ') + ')'
+            end.join(', ') + ')'
+          end.join(', ') + ')'
+        end
+
+        # for PostGIS ST_ASGeoJSON compatibility
+        def multi_polygon_to_ewkt(hash)
+          return 'MULTIPOLYGON EMPTY' if hash['coordinates'].nil?
+          'MULTIPOLYGON(' + hash['coordinates'].collect do |polygon|
+            '(' + polygon.collect do |hole|
+              '(' + hole.collect do |point|
+                point.join(' ')
+              end.join(', ') + ')'
+            end.join(', ') + ')'
+          end.join(', ') + ')'
+        end
+    end
+
+  end
+end

data/lib/charta/geo_json.rb

@@ -7,12 +7,14 @@ module Charta
 
     def initialize(data, srid = :WGS84)
       srid ||= :WGS84
-      @json = self.class.flatten(data.is_a?(Hash) ? data : JSON.parse(data))
+      @json = Coordinates.flatten(data.is_a?(Hash) ? data : JSON.parse(data))
       lsrid = @json['crs']['properties']['name'] if @json.is_a?(Hash) &&
                                                     @json['crs'].is_a?(Hash) &&
                                                     @json['crs']['properties'].is_a?(Hash)
       lsrid ||= srid
       @srid = ::Charta.find_srid(lsrid)
+
+      @json = Coordinates.normalize_4326_geometry(@json) if @srid.to_i == 4326
     end
 
     def geom
@@ -24,7 +26,7 @@ module Charta
     end
 
     def to_ewkt
-      "SRID=#{@srid};" + self.class.object_to_ewkt(@json)
+      "SRID=#{srid};" + EwktSerializer.object_to_ewkt(@json)
     end
 
     def valid?
@@ -42,131 +44,14 @@ module Charta
         false
       end
 
-      # Force coordinates to 2D
       def flatten(hash)
-        if hash['type'] == 'FeatureCollection'
-          flatten_feature_collection(hash)
-        elsif hash['type'] == 'Feature'
-          flatten_feature(hash)
-        else
-          flatten_geometry(hash)
-        end
-      end
-
-      def flatten_feature_collection(hash)
-        hash.merge('features' => hash['features'].map { |f| flatten_feature(f) })
-      end
-
-      def flatten_feature(hash)
-        hash.merge('geometry' => flatten_geometry(hash['geometry']))
-      end
-
-      def flatten_geometry(hash)
-        coordinates = hash['coordinates']
-        flattened =
-          case hash['type']
-          when 'Point' then
-            flatten_position(coordinates)
-          when 'MultiPoint', 'LineString'
-            coordinates.map { |p| flatten_position(p) }
-          when 'MultiLineString', 'Polygon'
-            coordinates.map { |l| l.map { |p| flatten_position(p) } }
-          when 'MultiPolygon'
-            coordinates.map { |m| m.map { |l| l.map { |p| flatten_position(p) } } }
-          when 'GeometryCollection' then
-            return hash.merge('geometries' => hash['geometries'].map { |g| flatten_geometry(g) })
-          else
-            raise StandardError, "Cannot handle: #{hash['type'].inspect}. In #{hash.inspect}"
-          end
-
-        hash.merge('coordinates' => flattened)
-      end
-
-      def flatten_position(position)
-        position[0..1]
-      end
-
-      def object_to_ewkt(hash)
-        type = hash[:type] || hash['type']
-        send("#{type.gsub(/(.)([A-Z])/, '\1_\2').downcase}_to_ewkt", hash)
-      end
-
-      def feature_collection_to_ewkt(hash)
-        return 'GEOMETRYCOLLECTION EMPTY' if hash['features'].nil?
-        'GEOMETRYCOLLECTION(' + hash['features'].collect do |feature|
-          object_to_ewkt(feature)
-        end.join(', ') + ')'
-      end
-
-      def geometry_collection_to_ewkt(hash)
-        return 'GEOMETRYCOLLECTION EMPTY' if hash['geometries'].nil?
-        'GEOMETRYCOLLECTION(' + hash['geometries'].collect do |feature|
-          object_to_ewkt(feature)
-        end.join(', ') + ')'
-      end
-
-      def feature_to_ewkt(hash)
-        object_to_ewkt(hash['geometry'])
-      end
-
-      def point_to_ewkt(hash)
-        return 'POINT EMPTY' if hash['coordinates'].nil?
-        'POINT(' + hash['coordinates'].join(' ') + ')'
-      end
-
-      def line_string_to_ewkt(hash)
-        return 'LINESTRING EMPTY' if hash['coordinates'].nil?
-        'LINESTRING(' + hash['coordinates'].collect do |point|
-          point.join(' ')
-        end.join(', ') + ')'
-      end
-
-      def polygon_to_ewkt(hash)
-        return 'POLYGON EMPTY' if hash['coordinates'].nil?
-        'POLYGON(' + hash['coordinates'].collect do |hole|
-          '(' + hole.collect do |point|
-            point.join(' ')
-          end.join(', ') + ')'
-        end.join(', ') + ')'
-      end
-
-      def multi_point_to_ewkt(hash)
-        return 'MULTIPOINT EMPTY' if hash['coordinates'].nil?
-        'MULTIPOINT(' + hash['coordinates'].collect do |point|
-          '(' + point.join(' ') + ')'
-        end.join(', ') + ')'
-      end
-
-      def multi_line_string_to_ewkt(hash)
-        return 'MULTILINESTRING EMPTY' if hash['coordinates'].nil?
-        'MULTILINESTRING(' + hash['coordinates'].collect do |line|
-          '(' + line.collect do |point|
-            point.join(' ')
-          end.join(', ') + ')'
-        end.join(', ') + ')'
-      end
-
-      def multipolygon_to_ewkt(hash)
-        return 'MULTIPOLYGON EMPTY' if hash['coordinates'].nil?
-        'MULTIPOLYGON(' + hash['coordinates'].collect do |polygon|
-          '(' + polygon.collect do |hole|
-            '(' + hole.collect do |point|
-              point.join(' ')
-            end.join(', ') + ')'
-          end.join(', ') + ')'
-        end.join(', ') + ')'
+        Coordinates.flatten hash
       end
 
-      # for PostGIS ST_ASGeoJSON compatibility
-      def multi_polygon_to_ewkt(hash)
-        return 'MULTIPOLYGON EMPTY' if hash['coordinates'].nil?
-        'MULTIPOLYGON(' + hash['coordinates'].collect do |polygon|
-          '(' + polygon.collect do |hole|
-            '(' + hole.collect do |point|
-              point.join(' ')
-            end.join(', ') + ')'
-          end.join(', ') + ')'
-        end.join(', ') + ')'
+      %i[object_to_ewkt feature_collection_to_ewkt geometry_collection_to_ewkt feature_to_ewkt point_to_ewkt line_string_to_ewkt polygon_to_ewkt multi_point_to_ewkt multi_line_string_to_ewkt multipolygon_to_ewkt multi_polygon_to_ewkt].each do |m|
+        define_method m do |*args|
+          EwktSerializer.send m, *args
+        end
       end
     end
   end

data/lib/charta/geometry.rb

@@ -43,6 +43,7 @@ module Charta
     def to_text
       feature.as_text.match(/\ASRID=.*;(.*)/)[1]
     end
+
     alias as_text to_text
     alias to_wkt to_text
 
@@ -50,6 +51,7 @@ module Charta
     def to_ewkt
       Charta.generate_ewkt(feature).to_s
     end
+
     alias to_s to_ewkt
 
     def ewkt
@@ -62,6 +64,7 @@ module Charta
       generator = RGeo::WKRep::WKBGenerator.new(tag_format: :ewkbt, emit_ewkbt_srid: true)
       generator.generate(feature)
     end
+
     alias to_ewkb to_binary
 
     # Pas bien compris le fonctionnement
@@ -85,6 +88,7 @@ module Charta
     def to_geojson
       to_json_object.to_json
     end
+
     alias to_json to_geojson
 
     # Returns object in JSON (Hash)
@@ -110,12 +114,24 @@ module Charta
 
     # Returns true if Geometry is a Surface
     def surface?
-      [RGeo::Feature::Polygon, RGeo::Feature::MultiPolygon].include? feature.geometry_type
+      if collection?
+        feature.any? { |geometry| Charta.new_geometry(geometry).surface? }
+      else
+        [RGeo::Feature::Polygon, RGeo::Feature::MultiPolygon].include? feature.geometry_type
+      end
     end
 
     # Returns area in unit corresponding to the SRS
     def area
-      surface? ? feature.area : 0
+      if surface?
+        if collection?
+          feature.sum { |geometry| Charta.new_geometry(geometry).area }
+        else
+          feature.area
+        end
+      else
+        0
+      end
     end
 
     # Returns true if this Geometry is an empty geometrycollection, polygon,
@@ -123,6 +139,7 @@ module Charta
     def empty?
       feature.is_empty?
     end
+
     alias blank? empty?
 
     # Computes the geometric center of a geometry, or equivalently, the center
@@ -142,10 +159,16 @@ module Charta
 
     def convert_to(new_type)
       case new_type
-      when type then self
-      when :multi_point then flatten_multi(:point)
-      when :multi_line_string then flatten_multi(:line_string)
-      when :multi_polygon then flatten_multi(:polygon)
+        when type then
+          self
+        when :multi_point then
+          flatten_multi(:point)
+        when :multi_line_string then
+          flatten_multi(:line_string)
+        when :multi_polygon then
+          flatten_multi(:polygon)
+        else
+          self
       end
     end
 
@@ -196,6 +219,7 @@ module Charta
       other_geometry = Charta.new_geometry(other).transform(srid)
       feature.union(other_geometry.feature)
     end
+
     alias + merge
 
     def intersection(other)
@@ -203,10 +227,16 @@ module Charta
       feature.intersection(other_geometry.feature)
     end
 
+    def intersects?(other)
+      other_geometry = Charta.new_geometry(other).transform(srid)
+      feature.intersects?(other_geometry.feature)
+    end
+
     def difference(other)
       other_geometry = Charta.new_geometry(other).transform(srid)
       feature.difference(other_geometry.feature)
     end
+
     alias - difference
 
     def bounding_box
@@ -288,52 +318,52 @@ module Charta
 
       private
 
-      def geos_factory(srid)
-        RGeo::Geos.factory(
-          srid: srid,
-          wkt_generator: {
-            type_format: :ewkt,
-            emit_ewkt_srid: true,
-            convert_case: :upper
-          },
-          wkt_parser: {
-            support_ewkt: true
-          },
-          wkb_generator:  {
-            type_format: :ewkb,
-            emit_ewkb_srid: true,
-            hex_format: true
-          },
-          wkb_parser: {
-            support_ewkb: true
-          }
-        )
-      end
+        def geos_factory(srid)
+          RGeo::Geos.factory(
+            srid: srid,
+            wkt_generator: {
+              type_format: :ewkt,
+              emit_ewkt_srid: true,
+              convert_case: :upper
+            },
+            wkt_parser: {
+              support_ewkt: true
+            },
+            wkb_generator: {
+              type_format: :ewkb,
+              emit_ewkb_srid: true,
+              hex_format: true
+            },
+            wkb_parser: {
+              support_ewkb: true
+            }
+          )
+        end
 
-      def projected_factory(srid)
-        proj4 = '+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs'
-        RGeo::Geographic.projected_factory(
-          srid: srid,
-          wkt_generator: {
-            type_format: :ewkt,
-            emit_ewkt_srid: true,
-            convert_case: :upper
-          },
-          wkt_parser: {
-            support_ewkt: true
-          },
-          wkb_generator:  {
-            type_format: :ewkb,
-            emit_ewkb_srid: true,
-            hex_format: true
-          },
-          wkb_parser: {
-            support_ewkb: true
-          },
-          projection_srid: 6933,
-          projection_proj4: proj4
-        )
-      end
+        def projected_factory(srid)
+          proj4 = '+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs'
+          RGeo::Geographic.projected_factory(
+            srid: srid,
+            wkt_generator: {
+              type_format: :ewkt,
+              emit_ewkt_srid: true,
+              convert_case: :upper
+            },
+            wkt_parser: {
+              support_ewkt: true
+            },
+            wkb_generator: {
+              type_format: :ewkb,
+              emit_ewkb_srid: true,
+              hex_format: true
+            },
+            wkb_parser: {
+              support_ewkb: true
+            },
+            projection_srid: 6933,
+            projection_proj4: proj4
+          )
+        end
     end
   end
 end

data/lib/charta/point.rb

@@ -10,5 +10,10 @@ module Charta
       feature.y
     end
     alias latitude y
+
+    def distance(point)
+      raise ArgumentError, "wrong type: Charta::Point required" if point.class.name != "Charta::Point"
+      to_rgeo.distance(point.to_rgeo)
+    end
   end
 end

data/lib/charta/polygon.rb

@@ -8,5 +8,10 @@ module Charta
       end
       @exterior_ring
     end
+
+    def distance(point)
+      polygon_centroid = Charta.new_point(*centroid, 4326)
+      polygon_centroid.distance(point)
+    end
   end
 end

data/lib/charta/version.rb

@@ -1,3 +1,3 @@
 module Charta
-  VERSION = '0.1.12'.freeze
+  VERSION = '0.1.17'.freeze
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: charta
 version: !ruby/object:Gem::Version
-  version: 0.1.12
+  version: 0.1.17
 platform: ruby
 authors:
 - Brice TEXIER
-autorequire: 
+autorequire:
 bindir: bin
 cert_chain: []
-date: 2019-07-01 00:00:00.000000000 Z
+date: 2020-09-08 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: nokogiri
@@ -122,7 +122,7 @@ dependencies:
     - - ">="
       - !ruby/object:Gem::Version
         version: '0'
-description: 
+description:
 email:
 - brice@ekylibre.com
 executables: []
@@ -130,6 +130,7 @@ extensions: []
 extra_rdoc_files: []
 files:
 - ".gitignore"
+- ".gitlab-ci.yml"
 - ".travis.yml"
 - CODE_OF_CONDUCT.md
 - Gemfile
@@ -139,6 +140,8 @@ files:
 - charta.gemspec
 - lib/charta.rb
 - lib/charta/bounding_box.rb
+- lib/charta/coordinates.rb
+- lib/charta/ewkt_serializer.rb
 - lib/charta/geo_json.rb
 - lib/charta/geojson_import.rb
 - lib/charta/geometry.rb
@@ -156,7 +159,7 @@ homepage: https://github.com/ekylibre/charta
 licenses:
 - MIT
 metadata: {}
-post_install_message: 
+post_install_message:
 rdoc_options: []
 require_paths:
 - lib
@@ -171,9 +174,8 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubyforge_project: 
-rubygems_version: 2.6.14
-signing_key: 
+rubygems_version: 3.0.3
+signing_key:
 specification_version: 4
 summary: Simple tool over geos and co
 test_files: []