libgd-gis 0.1.2 → 0.2.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: c2942d8d8fd3401c3c6c5bd47b8ea648aa63bc415e9621fb4c2f40873a80c7c7
-  data.tar.gz: d4a50c693ad4ab58e72365f8062d2420a8f683ccf692b97f4a7103f3a8ebe959
+  metadata.gz: 3c053cbc0c95a2ab748a9e2dca5bcfdfe2340569ebdb05c383a42ee6194c8df4
+  data.tar.gz: 81d4c110d6c26041c1b451499cec406f2d0b2326aef27fe430f6142c3df29daf
 SHA512:
-  metadata.gz: 7fc569940f1fea02dbe396b32a05286c68a07712b0a3c81665c1a5d3ac512dc1be2b76c349986b1d0a042afde4e900907aa51efb66b1cd2e43782b440aef6fc4
-  data.tar.gz: 3a062854e5e353a6d77632a162e26425071f808d6b8d84d08a60dc7e7f713edffbcc9ecff98351058cc868c49557a86b66e4da6a098502378310e4d05d2561d7
+  metadata.gz: 0e9accc97be10902e35d833587de5a0923c7a0e88b88944c9c3bc78cdda3c8cab0cd742c95563a6bfa177e48759a7232e636ed11ca1dc20ecd781deab8c859c6
+  data.tar.gz: f0d1a6d0a854cd046f29ccaffa125f48b86d147b965c43501c6cedeefb2b3ebf27962e055f6f266714c24b4b3e21cb57606611a594c2bfab30d0ec617d64da58

data/README.md

@@ -1,24 +1,56 @@
-# libgd-gis
-
-**A geospatial raster engine for Ruby**
-
-Render real maps, GeoJSON layers, and geospatial data directly in Ruby using a native raster backend powered by **libgd**.
-
-This project brings back something Ruby lost over time:  
-the ability to generate **maps, tiles, and geospatial visualizations natively**, without external tools like QGIS, ImageMagick, or Mapbox.
+# LibGD-GIS
+
+<p align="center">
+  <a href="https://rubystacknews.com/2026/01/07/ruby-can-now-draw-maps-and-i-started-with-ice-cream/">
+    <img src="https://img.shields.io/badge/RubyStackNews-CC342D?style=for-the-badge&logo=ruby&logoColor=white" />
+  </a>
+  <a href="https://x.com/ruby_stack_news">
+    <img src="https://img.shields.io/badge/Twitter%20@RubyStackNews-1DA1F2?style=for-the-badge&logo=twitter&logoColor=white" />
+  </a>
+  <a href="https://www.linkedin.com/in/germ%C3%A1n-silva-56a12622/">
+    <img src="https://img.shields.io/badge/Germán%20Silva-0A66C2?style=for-the-badge&logo=linkedin&logoColor=white" />
+  </a>
+</p>
+
+<p align="center">
+  <a href="https://rubygems.org/gems/libgd-gis">
+    <img src="https://img.shields.io/badge/RubyGems-libgd--gis-CC342D?style=for-the-badge&logo=rubygems&logoColor=white" />
+  </a>
+  <a href="https://github.com/ggerman/libgd-gis">
+    <img src="https://img.shields.io/badge/GitHub-libgd--gis-181717?style=for-the-badge&logo=github&logoColor=white" />
+  </a>
+  <a href="https://github.com/ggerman/ruby-libgd">
+    <img src="https://img.shields.io/badge/Engine-ruby--libgd-CC342D?style=for-the-badge&logo=ruby&logoColor=white" />
+  </a>
+</p>
+
+<p align="right">
+  <img src="docs/images/logo-gis.png" width="160" />
+</p>
+
+![CI](https://github.com/ggerman/libgd-gis/actions/workflows/ci.yml/badge.svg)
+
+| Examples | Examples | Examples |
+| :----: | :----: | :--: |
+| <img src="docs/examples/america.png" height="250"> | <img src="docs/examples/argentina_museum.png" height="250"> | <img src="docs/examples/museos_parana.png" height="250"> |
+| <img src="docs/examples/asia.png" height="250"> | <img src="docs/examples/europe.png" height="250"> | <img src="docs/examples/icecream_parana.png" height="250"> |
+| <img src="docs/examples/argentina_cities.png" height="250"> | <img src="docs/examples/tanzania_hydro.png" height="250"> | <img src="docs/examples/parana_polygon.png" height="250"> |
+| <img src="docs/examples/parana_carto_dark.png" height="250"> | <img src="docs/examples/ramirez_avenue.png" height="250"> | |
 
 ---
 
-## 🌍 Example
+## A geospatial raster engine for Ruby.
 
-This map of Tanzania was generated entirely in Ruby from a GeoJSON dataset of hydroelectric power plants:
+libgd-gis allows Ruby to render real maps, GeoJSON layers, vector features, and geospatial tiles using a native raster backend powered by **libgd**.
 
-![Tanzania hydro plants](examples/output/tanzania_hydro.png)
+It restores something Ruby lost over time:
+ the ability to generate **maps, tiles, and GIS-grade visualizations natively**, without relying on external tools like QGIS, Mapnik, ImageMagick, or Mapbox.
 
-No JavaScript.  
-No QGIS.  
-No Mapbox.  
-Just Ruby.
+Built on top of **ruby-libgd**, this project turns Ruby into a **map rendering engine**, capable of producing spatial graphics, tiled maps, and geospatial outputs directly inside Ruby processes.
+
+- No external renderers.
+-  No shelling out.
+-  Just Ruby, raster, and GIS.
 
 ---
 
@@ -75,29 +107,46 @@ Render hydroelectric plants from a GeoJSON file:
 require "json"
 require "gd/gis"
 
-geo = JSON.parse(File.read("hydro_plants.geojson"))
-plants = geo["features"]
-
-lons = plants.map { |f| f["geometry"]["coordinates"][0] }
-lats = plants.map { |f| f["geometry"]["coordinates"][1] }
-
-bbox = [lons.min, lats.min, lons.max, lats.max]
+# ---------------------------
+# Bounding box mundial
+# ---------------------------
+AMERICA = [-170, -60, -30, 75]
 
+# ---------------------------
+# Crear mapa
+# ---------------------------
 map = GD::GIS::Map.new(
-  bbox: bbox,
-  zoom: 7,
+  bbox: AMERICA,
+  zoom: 4,
   basemap: :carto_light
 )
 
+# Cargar datos
+# ---------------------------
+peaks = JSON.parse(File.read("picks.json"))
+
+# ---------------------------
+# Agregar capa de puntos
+# ---------------------------
 map.add_points(
-  plants,
-  lon: ->(f) { f["geometry"]["coordinates"][0] },
-  lat: ->(f) { f["geometry"]["coordinates"][1] },
-  icon: "hydro.png"
+  peaks,
+  lon: ->(p) { p["longitude"] },
+  lat: ->(p) { p["latitude"] },
+  icon: "peak.png",
+  label: ->(p) { p["name"] },
+  font: "./fonts/DejaVuSans.ttf",
+  size: 10,
+  color: [0,0,0]
 )
 
+# ---------------------------
+# Renderizar y guardar
+# ---------------------------
 map.render
-map.save("tanzania_hydro.png")
+map.save("output/america.png")
+
+puts "Saved output/america.png"
+
 ```
 
 ---
@@ -122,6 +171,6 @@ MIT
 
 ## Author
 
-Germán Silva  
-https://github.com/ggerman  
+Germán Silva
+https://github.com/ggerman
 https://rubystacknews.com

data/examples/argentina_cities.rb

@@ -0,0 +1,89 @@
+require "json"
+require "gd"
+
+# --------------------------------
+# Input / output
+# --------------------------------
+BASE_MAP = "output/argentina.png"
+DATA     = "localidades.json"
+OUTPUT   = "output/argentina_localidades.png"
+
+# --------------------------------
+# Bounding box de Argentina (WGS84)
+# --------------------------------
+ARG_MIN_LON = -73.6
+ARG_MAX_LON = -53.6
+ARG_MIN_LAT = -55.1
+ARG_MAX_LAT = -21.8
+
+R = 6378137.0
+
+def mercator_x(lon)
+  lon * Math::PI / 180.0 * R
+end
+
+def mercator_y(lat)
+  Math.log(Math.tan(Math::PI/4 + lat * Math::PI / 360.0)) * R
+end
+
+# --------------------------------
+# Load base map
+# --------------------------------
+img = GD::Image.open(BASE_MAP)
+
+width  = img.width
+height = img.height
+
+flag = GD::Image.open("flag.png")
+flag_w = flag.width
+flag_h = flag.height
+
+# --------------------------------
+# Projected bounding box
+# --------------------------------
+min_x = mercator_x(ARG_MIN_LON)
+max_x = mercator_x(ARG_MAX_LON)
+min_y = mercator_y(ARG_MIN_LAT)
+max_y = mercator_y(ARG_MAX_LAT)
+
+def 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
+
+# --------------------------------
+# Load locations
+# --------------------------------
+data = JSON.parse(File.read(DATA))
+
+red = GD::Color.rgb(220, 40, 40)
+
+puts "Plotting #{data["localidades"].size} localidades..."
+
+data["localidades"].each do |loc|
+  lon = loc["centroide"]["lon"]
+  lat = loc["centroide"]["lat"]
+
+  x, y = lonlat_to_pixel(lon, lat, min_x, max_x, min_y, max_y, width, height)
+
+  # dibujar un punto
+    img.copy(
+        flag,
+        x - flag_w / 2,
+        y - flag_h / 2,
+        0, 0,
+        flag_w,
+        flag_h
+    )
+end
+
+# --------------------------------
+# Save result
+# --------------------------------
+img.save(OUTPUT)
+puts "Saved #{OUTPUT}"

data/examples/argentina_museums.rb

@@ -0,0 +1,70 @@
+require "csv"
+require "gd"
+require "pry"
+
+ZOOM = 7
+TILE_SIZE = 256
+
+# Bounding box Argentina (igual que en argentina_tiles.rb)
+WEST  = -73.6
+EAST  = -53.6
+NORTH = -21.8
+SOUTH = -55.1
+
+def lon2tile(lon, z)
+  ((lon + 180.0) / 360.0 * (2 ** z)).floor
+end
+
+def lat2tile(lat, z)
+  rad = lat * Math::PI / 180
+  ((1 - Math.log(Math.tan(rad) + 1 / Math.cos(rad)) / Math::PI) / 2 * (2 ** z)).floor
+end
+
+def lon2px(lon, z)
+  ((lon + 180.0) / 360.0 * (256 * 2**z))
+end
+
+def lat2px(lat, z)
+  rad = lat * Math::PI / 180
+  ((1 - Math.log(Math.tan(rad) + 1 / Math.cos(rad)) / Math::PI) / 2 * (256 * 2**z))
+end
+
+# Same tile bbox used to build argentina.png
+x_min = lon2tile(WEST,  ZOOM)
+y_min = lat2tile(NORTH, ZOOM)
+
+offset_x = x_min * TILE_SIZE
+offset_y = y_min * TILE_SIZE
+
+# Load map
+img = GD::Image.open("output/argentina.png")
+
+# Load icon
+icon = GD::Image.open("museos.png")
+iw = icon.width
+ih = icon.height
+
+puts "Plotting museums..."
+
+CSV.foreach("museos_datosabiertos.csv", headers: true, encoding: "bom|utf-8") do |row|
+  next if row["Latitud"].nil? || row["Longitud"].nil?
+
+  lat = row["Latitud"].to_f
+  lon = row["Longitud"].to_f
+
+  # global pixel in Web Mercator
+  px = lon2px(lon, ZOOM)
+  py = lat2px(lat, ZOOM)
+
+  # local pixel inside argentina.png
+  x = px - offset_x
+  y = py - offset_y
+
+  # skip if outside raster
+  next if x < 0 || y < 0 || x >= img.width || y >= img.height
+
+  img.copy(icon, x.to_i - iw/2, y.to_i - ih/2, 0, 0, iw, ih)
+end
+
+img.save("output/argentina_museums.png")
+puts "Saved output/argentina_museums.png"

data/examples/asia-peaks.rb

@@ -0,0 +1,42 @@
+require "json"
+require "gd/gis"
+
+# ---------------------------
+# Bounding box mundial
+# ---------------------------
+ASIA = [25.0, -10.0, 150.0, 60.0]
+
+# ---------------------------
+# Create map
+# ---------------------------
+map = GD::GIS::Map.new(
+  bbox: ASIA,
+  zoom: 4,
+  basemap: :carto_light
+)
+
+# Cargar datos
+# ---------------------------
+peaks = JSON.parse(File.read("picks.json"))
+
+# ---------------------------
+# Add points layer
+# ---------------------------
+map.add_points(
+  peaks,
+  lon: ->(p) { p["longitude"] },
+  lat: ->(p) { p["latitude"] },
+  icon: "peak.png",
+  label: ->(p) { p["name"] },
+  font: "./fonts/DejaVuSans.ttf",
+  size: 10,
+  color: [0,0,0]
+)
+
+# ---------------------------
+# Render and save
+# ---------------------------
+map.render
+map.save("output/asia.png")
+
+puts "Saved output/asia.png"

data/examples/avenue_ramirez.rb

@@ -0,0 +1,44 @@
+require "json"
+require "gd/gis"
+
+# Paraná bbox
+BBOX = [-60.556640625, -31.8402326679, -60.46875, -31.6907818061]
+
+map = GD::GIS::Map.new(
+  bbox: BBOX,
+  zoom: 14,
+  basemap: :esri_satellite
+)
+
+geo = JSON.parse(File.read("ramirez_full.geojson"))
+features = geo["features"]
+
+map.add_lines(
+  features,
+  stroke: [0, 0, 0, 90],       # borde
+  fill:   [0, 0, 0, 90],   # calzada
+  width:  50                   # ancho en metros
+)
+
+map.add_lines(
+  features,
+  stroke: [255, 165, 0, 90],
+  width:  2
+)
+
+# Label
+map.add_points(
+  [{ lon: -60.5205, lat: -31.76, name: "Av. Ramírez" }],
+  lon: ->(p){ p[:lon] },
+  lat: ->(p){ p[:lat] },
+  label: ->(p){ p[:name] },
+  font: "./fonts/DejaVuSans.ttf",
+  size: 16,
+  color: [0,0,0,160],
+  icon: "mark.png"
+)
+
+map.render
+map.save("ramirez_gis.png")
+
+puts "Generated ramirez_gis.png"

data/examples/europe-peaks.rb

@@ -0,0 +1,29 @@
+require "json"
+require "gd/gis"
+
+EUROPE = [-15.0, 34.0, 40.0, 72.0]
+
+map = GD::GIS::Map.new(
+  bbox: EUROPE,
+  zoom: 4,
+  basemap: :carto_light
+)
+
+
+peaks = JSON.parse(File.read("picks.json"))
+
+map.add_points(
+  peaks,
+  lon: ->(p) { p["longitude"] },
+  lat: ->(p) { p["latitude"] },
+  icon: "peak.png",
+  label: ->(p) { p["name"] },
+  font: "./fonts/DejaVuSans.ttf",
+  size: 10,
+  color: [0,0,0]
+)
+
+map.render
+map.save("output/europe.png")
+
+puts "Saved output/europe.png"

data/examples/ice-cream-shops.rb

@@ -0,0 +1,37 @@
+require "json"
+require "gd/gis"
+
+CITY = [-60.56, -31.76, -60.50, -31.71]
+
+map = GD::GIS::Map.new(
+  bbox: CITY,
+  zoom: 15,
+  basemap: :carto_light
+)
+
+# --------------------------------
+# Cargar datos
+# --------------------------------
+shops = JSON.parse(File.read("heladerias.json"))
+
+# --------------------------------
+# Ice cream shops layer
+# --------------------------------
+map.add_points(
+  shops,
+  lon: ->(s) { s["lon"] },
+  lat: ->(s) { s["lat"] },
+  icon: "ice-cream.png",
+  label: ->(s) { s["name"] },
+  font: "./fonts/DejaVuSans.ttf",
+  size: 11,
+  color: [40,40,40]
+)
+
+# --------------------------------
+# Render
+# --------------------------------
+map.render
+map.save("output/heladerias_parana.png")
+
+puts "Saved output/heladerias_parana.png"

data/examples/mountains.rb

@@ -0,0 +1,98 @@
+require "net/http"
+require "uri"
+require "fileutils"
+require "json"
+require "gd"
+
+ZOOM = 4
+TILE_SIZE = 256
+
+X_MIN = 0
+Y_MIN = 0
+X_MAX = (2 ** ZOOM) - 1
+Y_MAX = (2 ** ZOOM) - 1
+
+COLS = X_MAX - X_MIN + 1
+ROWS = Y_MAX - Y_MIN + 1
+
+WIDTH  = COLS * TILE_SIZE
+HEIGHT = ROWS * TILE_SIZE
+
+FONT = "./fonts/DejaVuSans-Bold.ttf"
+
+FileUtils.mkdir_p("tmp/tiles")
+FileUtils.mkdir_p("output")
+
+def lon2px(lon, z)
+  ((lon + 180.0) / 360.0 * (256 * 2**z))
+end
+
+def lat2px(lat, z)
+  rad = lat * Math::PI / 180
+  ((1 - Math.log(Math.tan(rad) + 1 / Math.cos(rad)) / Math::PI) / 2 * (256 * 2**z))
+end
+
+def fetch_tile(z, x, y)
+  path = "tmp/tiles/#{z}_#{x}_#{y}.png"
+  return path if File.exist?(path)
+
+  url = URI("https://basemaps.cartocdn.com/light_all/#{z}/#{x}/#{y}.png")
+  puts "Downloading #{url}"
+
+  Net::HTTP.start(url.host, url.port, use_ssl: true) do |http|
+    req = Net::HTTP::Get.new(url)
+    req["User-Agent"] = "ruby-libgd-gis"
+    res = http.request(req)
+    raise "Tile failed #{z}/#{x}/#{y}" unless res.code == "200"
+    File.binwrite(path, res.body)
+  end
+
+  path
+end
+
+puts "Building world map #{WIDTH}x#{HEIGHT}..."
+
+map = GD::Image.new(WIDTH, HEIGHT)
+
+(X_MIN..X_MAX).each_with_index do |x, cx|
+  (Y_MIN..Y_MAX).each_with_index do |y, cy|
+    tile = GD::Image.open(fetch_tile(ZOOM, x, y))
+    map.copy(tile, cx * TILE_SIZE, cy * TILE_SIZE, 0, 0, TILE_SIZE, TILE_SIZE)
+  end
+end
+
+peaks = JSON.parse(File.read("picks.json"))
+
+red   = GD::Color.rgb(220, 40, 40)
+black = GD::Color.rgb(0, 0, 0)
+
+puts "Plotting peaks..."
+
+peaks.each do |p|
+  lon  = p["longitude"]
+  lat  = p["latitude"]
+  name = p["name"]
+
+  px = lon2px(lon, ZOOM)
+  py = lat2px(lat, ZOOM)
+
+  x = px.to_i
+  y = py.to_i
+
+  map.filled_circle(x, y, 7, red)
+
+  map.text(
+    name,
+    x: x + 10,
+    y: y - 5,
+    size: 10,
+    color: black,
+    font: FONT
+  )
+end
+
+# =========================
+# Save
+# =========================
+map.save("output/world_peaks.png")
+puts "Saved output/world_peaks.png"

data/examples/parana-city.rb

@@ -0,0 +1,29 @@
+require "gd/gis"
+
+BBOX = [-60.53, -31.77, -60.51, -31.74]
+
+map = GD::GIS::Map.new(
+  bbox: BBOX,
+  zoom: 15,
+  basemap: :carto_dark
+)
+
+polygon = [
+  [-60.525, -31.750],
+  [-60.520, -31.750],
+  [-60.520, -31.755],
+  [-60.525, -31.755],
+  [-60.525, -31.750]
+]
+
+map.add_polygons(
+  [polygon],
+  fill:   [255, 140, 0, 120],
+  stroke: [0, 0, 0, 200],
+  width:  2
+)
+
+map.render
+map.save("carto_dark.png")
+
+puts "Generated carto_dark.png"