libgd-gis 0.1.0

data/lib/gd/gis/basemap.rb

@@ -0,0 +1,61 @@
+require "net/http"
+require "fileutils"
+
+module GD
+  module GIS
+    class Basemap
+      TILE_SIZE = 256
+
+      def initialize(zoom, bbox, provider=:carto_light)
+        @zoom = zoom
+        @bbox = bbox
+        @provider = provider
+      end
+
+      def url(z,x,y)
+        "https://basemaps.cartocdn.com/light_all/#{z}/#{x}/#{y}.png"
+      end
+
+      def lon2tile(lon)
+        ((lon + 180.0) / 360.0 * (2 ** @zoom)).floor
+      end
+
+      def lat2tile(lat)
+        rad = lat * Math::PI / 180
+        ((1 - Math.log(Math.tan(rad) + 1 / Math.cos(rad)) / Math::PI) / 2 * (2 ** @zoom)).floor
+      end
+
+      def fetch_tiles
+        west, south, east, north = @bbox
+
+        x_min = lon2tile(west)
+        x_max = lon2tile(east)
+        y_min = lat2tile(north)
+        y_max = lat2tile(south)
+
+        @x_min = x_min
+        @y_min = y_min
+
+        FileUtils.mkdir_p("tmp/tiles")
+
+        tiles = []
+
+        (x_min..x_max).each do |x|
+          (y_min..y_max).each do |y|
+            path = "tmp/tiles/#{@zoom}_#{x}_#{y}.png"
+            unless File.exist?(path)
+              uri = URI(url(@zoom,x,y))
+              Net::HTTP.start(uri.host, uri.port, use_ssl: true) do |http|
+                res = http.get(uri.path)
+                File.binwrite(path, res.body)
+              end
+            end
+            tiles << [x,y,path]
+          end
+        end
+
+        [tiles, x_min, y_min]
+      end
+    end
+  end
+end

data/lib/gd/gis/layer_points.rb

@@ -0,0 +1,25 @@
+module GD
+  module GIS
+    class PointsLayer
+      def initialize(data, lon:, lat:, icon:)
+        @data = data
+        @lon = lon
+        @lat = lat
+        @icon = GD::Image.open(icon)
+      end
+
+      def render!(img, projection)
+        w = @icon.width
+        h = @icon.height
+
+        @data.each do |row|
+          lon = @lon.call(row)
+          lat = @lat.call(row)
+
+          x,y = projection.call(lon,lat)
+          img.copy(@icon, x - w/2, y - h/2, 0,0,w,h)
+        end
+      end
+    end
+  end
+end

data/lib/gd/gis/map.rb

@@ -0,0 +1,108 @@
+require "gd"
+require_relative "basemap"
+require_relative "projection"
+require_relative "layer_points"
+
+module GD
+  module GIS
+    class Map
+      def initialize(bbox:, zoom:, basemap:)
+        @bbox = bbox
+        @zoom = zoom
+        @basemap = Basemap.new(zoom, bbox, basemap)
+        @layers = []
+      end
+
+      def add_points(data, lon:, lat:, icon:)
+        @layers << PointsLayer.new(data, lon: lon, lat: lat, icon: icon)
+      end
+
+      def render
+        tiles, x_min, y_min = @basemap.fetch_tiles
+
+        cols = tiles.map{|t| t[0]}.uniq.size
+        rows = tiles.map{|t| t[1]}.uniq.size
+
+        width = cols * 256
+        height = rows * 256
+
+        @img = GD::Image.new(width,height)
+
+        tiles.each do |x,y,path|
+          tile = GD::Image.open(path)
+          cx = x - x_min
+          cy = y - y_min
+          @img.copy(tile, cx*256, cy*256, 0,0,256,256)
+        end
+
+        # projection bbox
+        min_x = Projection.mercator_x(@bbox[0])
+        max_x = Projection.mercator_x(@bbox[2])
+        min_y = Projection.mercator_y(@bbox[1])
+        max_y = Projection.mercator_y(@bbox[3])
+
+        projector = ->(lon,lat) {
+          Projection.lonlat_to_pixel(lon,lat,min_x,max_x,min_y,max_y,width,height)
+        }
+
+        @layers.each { |l| l.render!(@img, projector) }
+      end
+
+      def tile2lon(x, z)
+        x.to_f / (2**z) * 360.0 - 180.0
+      end
+
+      def tile2lat(y, z)
+        n = Math::PI - (2.0 * Math::PI * y.to_f) / (2**z)
+        180.0 / Math::PI * Math.atan(0.5 * (Math.exp(n) - Math.exp(-n)))
+      end
+
+      def tile_bbox(z,x,y)
+        west  = tile2lon(x, z)
+        east  = tile2lon(x+1, z)
+        north = tile2lat(y, z)
+        south = tile2lat(y+1, z)
+        [west, south, east, north]
+      end
+
+      def intersect_bbox(a,b)
+        [
+          [a[0], b[0]].max,
+          [a[1], b[1]].max,
+          [a[2], b[2]].min,
+          [a[3], b[3]].min
+        ]
+      end
+
+      def tile(z,x,y)
+        tb = tile_bbox(z,x,y)
+        ib = intersect_bbox(@bbox, tb)
+
+        # Si no hay intersección, devolvemos tile vacío
+        if ib[0] >= ib[2] || ib[1] >= ib[3]
+          return GD::Image.new(256,256)
+        end
+
+        submap = GD::GIS::Map.new(
+          bbox: ib,
+          zoom: z,
+          basemap: :carto_light
+        )
+
+        # copiar capas
+        @layers.each { |l| submap.instance_variable_get(:@layers) << l }
+
+        submap.render
+
+        # recortar al tamaño exacto del tile
+        img = GD::Image.new(256,256)
+        img.copy(submap.instance_variable_get(:@img), 0,0,0,0,256,256)
+        img
+      end
+
+      def save(path)
+        @img.save(path)
+      end
+    end
+  end
+end

data/lib/gd/gis/projection.rb

@@ -0,0 +1,25 @@
+module GD
+  module GIS
+    module Projection
+      R = 6378137.0
+
+      def self.mercator_x(lon)
+        lon * Math::PI / 180.0 * R
+      end
+
+      def self.mercator_y(lat)
+        Math.log(Math.tan(Math::PI/4 + lat * Math::PI / 360.0)) * R
+      end
+
+      def self.lonlat_to_pixel(lon, lat, min_x, max_x, min_y, max_y, width, height)
+        x = mercator_x(lon)
+        y = mercator_y(lat)
+
+        px = (x - min_x) / (max_x - min_x) * width
+        py = height - (y - min_y) / (max_y - min_y) * height
+
+        [px.to_i, py.to_i]
+      end
+    end
+  end
+end

data/lib/gd/gis.rb

@@ -0,0 +1,3 @@
+# lib/gd/gis.rb
+require "gd"
+require_relative "gis/map"

data/lib/libgd_gis.rb

@@ -0,0 +1,44 @@
+require "libgd_gis"
+
+require "open-uri"
+require "tempfile"
+require "gd"
+
+module LibGD
+  module GIS
+    class Tile
+        attr_reader :z, :x, :y, :image
+
+        def self.osm(z:, x:, y:)
+        new(
+            z: z,
+            x: x,
+            y: y,
+            source: "https://api.maptiler.com/maps/basic/#{z}/#{x}/#{y}.png?key=GetYourOwnKey"
+        )
+        end
+
+        def initialize(z:, x:, y:, source:)
+            @z = z
+            @x = x
+            @y = y
+            @source = source
+        end
+
+        def render
+            tmp = Tempfile.new(["tile", ".png"])
+            tmp.binmode
+            tmp.write URI.open(@source).read
+            tmp.flush
+
+            @image = GD::Image.open(tmp.path)
+        ensure
+            tmp.close
+        end
+
+        def save(path)
+            @image.save(path)
+        end
+    end
+  end
+end

metadata

@@ -0,0 +1,281 @@
+--- !ruby/object:Gem::Specification
+name: libgd-gis
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- Germán Alberto Giménez Silva
+autorequire:
+bindir: bin
+cert_chain: []
+date: 2026-01-06 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: ruby-libgd
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+description: A native GIS raster engine for Ruby built on libgd. Render maps, GeoJSON,
+  heatmaps and tiles.
+email:
+- ggerman@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- README.md
+- examples/argentina.png
+- examples/argentina.rb
+- examples/argentina_con_localidades.rb
+- examples/argentina_museos.rb
+- examples/flag.png
+- examples/hydro.png
+- examples/hydro_plants.geojson
+- examples/localidades.json
+- examples/museos.png
+- examples/museos_datosabiertos.csv
+- examples/output/argentina.png
+- examples/output/argentina_localidades.png
+- examples/output/argentina_museos.png
+- examples/output/tanzania_hydro.png
+- examples/render_tile.rb
+- examples/tanzania.rb
+- examples/target.rb
+- examples/tmp/tiles/7_27_43.png
+- examples/tmp/tiles/7_27_44.png
+- examples/tmp/tiles/7_27_45.png
+- examples/tmp/tiles/7_27_46.png
+- examples/tmp/tiles/7_27_47.png
+- examples/tmp/tiles/7_27_48.png
+- examples/tmp/tiles/7_27_49.png
+- examples/tmp/tiles/7_27_50.png
+- examples/tmp/tiles/7_28_43.png
+- examples/tmp/tiles/7_28_44.png
+- examples/tmp/tiles/7_28_45.png
+- examples/tmp/tiles/7_28_46.png
+- examples/tmp/tiles/7_28_47.png
+- examples/tmp/tiles/7_28_48.png
+- examples/tmp/tiles/7_28_49.png
+- examples/tmp/tiles/7_28_50.png
+- examples/tmp/tiles/7_29_43.png
+- examples/tmp/tiles/7_29_44.png
+- examples/tmp/tiles/7_29_45.png
+- examples/tmp/tiles/7_29_46.png
+- examples/tmp/tiles/7_29_47.png
+- examples/tmp/tiles/7_29_48.png
+- examples/tmp/tiles/7_29_49.png
+- examples/tmp/tiles/7_29_50.png
+- examples/tmp/tiles/7_30_43.png
+- examples/tmp/tiles/7_30_44.png
+- examples/tmp/tiles/7_30_45.png
+- examples/tmp/tiles/7_30_46.png
+- examples/tmp/tiles/7_30_47.png
+- examples/tmp/tiles/7_30_48.png
+- examples/tmp/tiles/7_30_49.png
+- examples/tmp/tiles/7_30_50.png
+- examples/tmp/tiles/7_31_43.png
+- examples/tmp/tiles/7_31_44.png
+- examples/tmp/tiles/7_31_45.png
+- examples/tmp/tiles/7_31_46.png
+- examples/tmp/tiles/7_31_47.png
+- examples/tmp/tiles/7_31_48.png
+- examples/tmp/tiles/7_31_49.png
+- examples/tmp/tiles/7_31_50.png
+- examples/tmp/tiles/7_32_43.png
+- examples/tmp/tiles/7_32_44.png
+- examples/tmp/tiles/7_32_45.png
+- examples/tmp/tiles/7_32_46.png
+- examples/tmp/tiles/7_32_47.png
+- examples/tmp/tiles/7_32_48.png
+- examples/tmp/tiles/7_32_49.png
+- examples/tmp/tiles/7_32_50.png
+- examples/tmp/tiles/7_37_71.png
+- examples/tmp/tiles/7_37_72.png
+- examples/tmp/tiles/7_37_73.png
+- examples/tmp/tiles/7_37_74.png
+- examples/tmp/tiles/7_37_75.png
+- examples/tmp/tiles/7_37_76.png
+- examples/tmp/tiles/7_37_77.png
+- examples/tmp/tiles/7_37_78.png
+- examples/tmp/tiles/7_37_79.png
+- examples/tmp/tiles/7_37_80.png
+- examples/tmp/tiles/7_37_81.png
+- examples/tmp/tiles/7_37_82.png
+- examples/tmp/tiles/7_37_83.png
+- examples/tmp/tiles/7_37_84.png
+- examples/tmp/tiles/7_37_85.png
+- examples/tmp/tiles/7_37_86.png
+- examples/tmp/tiles/7_37_87.png
+- examples/tmp/tiles/7_38_71.png
+- examples/tmp/tiles/7_38_72.png
+- examples/tmp/tiles/7_38_73.png
+- examples/tmp/tiles/7_38_74.png
+- examples/tmp/tiles/7_38_75.png
+- examples/tmp/tiles/7_38_76.png
+- examples/tmp/tiles/7_38_77.png
+- examples/tmp/tiles/7_38_78.png
+- examples/tmp/tiles/7_38_79.png
+- examples/tmp/tiles/7_38_80.png
+- examples/tmp/tiles/7_38_81.png
+- examples/tmp/tiles/7_38_82.png
+- examples/tmp/tiles/7_38_83.png
+- examples/tmp/tiles/7_38_84.png
+- examples/tmp/tiles/7_38_85.png
+- examples/tmp/tiles/7_38_86.png
+- examples/tmp/tiles/7_38_87.png
+- examples/tmp/tiles/7_39_71.png
+- examples/tmp/tiles/7_39_72.png
+- examples/tmp/tiles/7_39_73.png
+- examples/tmp/tiles/7_39_74.png
+- examples/tmp/tiles/7_39_75.png
+- examples/tmp/tiles/7_39_76.png
+- examples/tmp/tiles/7_39_77.png
+- examples/tmp/tiles/7_39_78.png
+- examples/tmp/tiles/7_39_79.png
+- examples/tmp/tiles/7_39_80.png
+- examples/tmp/tiles/7_39_81.png
+- examples/tmp/tiles/7_39_82.png
+- examples/tmp/tiles/7_39_83.png
+- examples/tmp/tiles/7_39_84.png
+- examples/tmp/tiles/7_39_85.png
+- examples/tmp/tiles/7_39_86.png
+- examples/tmp/tiles/7_39_87.png
+- examples/tmp/tiles/7_40_71.png
+- examples/tmp/tiles/7_40_72.png
+- examples/tmp/tiles/7_40_73.png
+- examples/tmp/tiles/7_40_74.png
+- examples/tmp/tiles/7_40_75.png
+- examples/tmp/tiles/7_40_76.png
+- examples/tmp/tiles/7_40_77.png
+- examples/tmp/tiles/7_40_78.png
+- examples/tmp/tiles/7_40_79.png
+- examples/tmp/tiles/7_40_80.png
+- examples/tmp/tiles/7_40_81.png
+- examples/tmp/tiles/7_40_82.png
+- examples/tmp/tiles/7_40_83.png
+- examples/tmp/tiles/7_40_84.png
+- examples/tmp/tiles/7_40_85.png
+- examples/tmp/tiles/7_40_86.png
+- examples/tmp/tiles/7_40_87.png
+- examples/tmp/tiles/7_41_71.png
+- examples/tmp/tiles/7_41_72.png
+- examples/tmp/tiles/7_41_73.png
+- examples/tmp/tiles/7_41_74.png
+- examples/tmp/tiles/7_41_75.png
+- examples/tmp/tiles/7_41_76.png
+- examples/tmp/tiles/7_41_77.png
+- examples/tmp/tiles/7_41_78.png
+- examples/tmp/tiles/7_41_79.png
+- examples/tmp/tiles/7_41_80.png
+- examples/tmp/tiles/7_41_81.png
+- examples/tmp/tiles/7_41_82.png
+- examples/tmp/tiles/7_41_83.png
+- examples/tmp/tiles/7_41_84.png
+- examples/tmp/tiles/7_41_85.png
+- examples/tmp/tiles/7_41_86.png
+- examples/tmp/tiles/7_41_87.png
+- examples/tmp/tiles/7_42_71.png
+- examples/tmp/tiles/7_42_72.png
+- examples/tmp/tiles/7_42_73.png
+- examples/tmp/tiles/7_42_74.png
+- examples/tmp/tiles/7_42_75.png
+- examples/tmp/tiles/7_42_76.png
+- examples/tmp/tiles/7_42_77.png
+- examples/tmp/tiles/7_42_78.png
+- examples/tmp/tiles/7_42_79.png
+- examples/tmp/tiles/7_42_80.png
+- examples/tmp/tiles/7_42_81.png
+- examples/tmp/tiles/7_42_82.png
+- examples/tmp/tiles/7_42_83.png
+- examples/tmp/tiles/7_42_84.png
+- examples/tmp/tiles/7_42_85.png
+- examples/tmp/tiles/7_42_86.png
+- examples/tmp/tiles/7_42_87.png
+- examples/tmp/tiles/7_43_71.png
+- examples/tmp/tiles/7_43_72.png
+- examples/tmp/tiles/7_43_73.png
+- examples/tmp/tiles/7_43_74.png
+- examples/tmp/tiles/7_43_75.png
+- examples/tmp/tiles/7_43_76.png
+- examples/tmp/tiles/7_43_77.png
+- examples/tmp/tiles/7_43_78.png
+- examples/tmp/tiles/7_43_79.png
+- examples/tmp/tiles/7_43_80.png
+- examples/tmp/tiles/7_43_81.png
+- examples/tmp/tiles/7_43_82.png
+- examples/tmp/tiles/7_43_83.png
+- examples/tmp/tiles/7_43_84.png
+- examples/tmp/tiles/7_43_85.png
+- examples/tmp/tiles/7_43_86.png
+- examples/tmp/tiles/7_43_87.png
+- examples/tmp/tiles/7_44_71.png
+- examples/tmp/tiles/7_44_72.png
+- examples/tmp/tiles/7_44_73.png
+- examples/tmp/tiles/7_44_74.png
+- examples/tmp/tiles/7_44_75.png
+- examples/tmp/tiles/7_44_76.png
+- examples/tmp/tiles/7_44_77.png
+- examples/tmp/tiles/7_44_78.png
+- examples/tmp/tiles/7_44_79.png
+- examples/tmp/tiles/7_44_80.png
+- examples/tmp/tiles/7_44_81.png
+- examples/tmp/tiles/7_44_82.png
+- examples/tmp/tiles/7_44_83.png
+- examples/tmp/tiles/7_44_84.png
+- examples/tmp/tiles/7_44_85.png
+- examples/tmp/tiles/7_44_86.png
+- examples/tmp/tiles/7_44_87.png
+- examples/tmp/tiles/7_75_64.png
+- examples/tmp/tiles/7_75_65.png
+- examples/tmp/tiles/7_75_66.png
+- examples/tmp/tiles/7_75_67.png
+- examples/tmp/tiles/7_75_68.png
+- examples/tmp/tiles/7_76_64.png
+- examples/tmp/tiles/7_76_65.png
+- examples/tmp/tiles/7_76_66.png
+- examples/tmp/tiles/7_76_67.png
+- examples/tmp/tiles/7_76_68.png
+- examples/tmp/tiles/7_77_64.png
+- examples/tmp/tiles/7_77_65.png
+- examples/tmp/tiles/7_77_66.png
+- examples/tmp/tiles/7_77_67.png
+- examples/tmp/tiles/7_77_68.png
+- lib/gd/gis.rb
+- lib/gd/gis/basemap.rb
+- lib/gd/gis/layer_points.rb
+- lib/gd/gis/map.rb
+- lib/gd/gis/projection.rb
+- lib/libgd_gis.rb
+homepage: https://github.com/ggerman/libgd-gis
+licenses:
+- MIT
+metadata: {}
+post_install_message:
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubygems_version: 3.5.22
+signing_key:
+specification_version: 4
+summary: Geospatial raster rendering for Ruby using libgd
+test_files: []