libgd-gis 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 5e2c3fe258b7fdb52d1775d50e18095333bf692416c4b4fadf3659febf174604
+  data.tar.gz: 3e1ebd506ff4469f231ada7621a5a10a395464d2c373679ee81aacf4e6c1dd35
+SHA512:
+  metadata.gz: c9bf3495bb2918619807b7b26a3f26b8cb523e839ec0813b7fdb369f9b26ce78add87f0d322e13acff7c6e265093b7852374482afad538bd90010aba13f5b12e
+  data.tar.gz: bd5c40732c0126cc2b69c195cd675f4108d0dfc107da5acd8df00834e61d1d575a8bf1ce4b256a635e8fe80f6b07aafd8cc4fe9a45997805079a6361eb38171d

data/README.md

@@ -0,0 +1,127 @@
+# 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.
+
+---
+
+## 🌍 Example
+
+This map of Tanzania was generated entirely in Ruby from a GeoJSON dataset of hydroelectric power plants:
+
+![Tanzania hydro plants](examples/output/tanzania_hydro.png)
+
+No JavaScript.  
+No QGIS.  
+No Mapbox.  
+Just Ruby.
+
+---
+
+## What is this?
+
+`libgd-gis` is a **geospatial rendering engine** for Ruby built on top of [`ruby-libgd`](https://github.com/ggerman/ruby-libgd).
+
+It allows you to:
+
+- Load GeoJSON, CSV, or any dataset with coordinates  
+- Fetch real basemap tiles  
+- Reproject WGS84 (lat/lon) into Web Mercator  
+- Render points, icons, and layers onto a raster map  
+- Generate PNG maps or map tiles  
+
+This is the same type of pipeline used by professional GIS systems — implemented in Ruby.
+
+---
+
+## Installation
+
+### System dependency
+
+`libgd-gis` depends on **libgd**, via `ruby-libgd`.
+
+Install libgd first:
+
+**Ubuntu / Debian**
+```
+sudo apt install libgd-dev
+```
+
+**macOS**
+```
+brew install gd
+```
+
+---
+
+### Ruby gems
+
+```
+gem install ruby-libgd
+gem install libgd-gis
+```
+
+---
+
+## Quick Example
+
+Render hydroelectric plants from a GeoJSON file:
+
+```ruby
+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]
+
+map = GD::GIS::Map.new(
+  bbox: bbox,
+  zoom: 7,
+  basemap: :carto_light
+)
+
+map.add_points(
+  plants,
+  lon: ->(f) { f["geometry"]["coordinates"][0] },
+  lat: ->(f) { f["geometry"]["coordinates"][1] },
+  icon: "hydro.png"
+)
+
+map.render
+map.save("tanzania_hydro.png")
+```
+
+---
+
+## Features
+
+- Real basemap tiles  
+- WGS84 → Web Mercator projection  
+- GeoJSON point rendering  
+- CSV / JSON support  
+- Icon-based symbol layers  
+- Automatic bounding box fitting  
+- Raster output (PNG)  
+
+---
+
+## License
+
+MIT
+
+---
+
+## Author
+
+Germán Silva  
+https://github.com/ggerman  
+https://rubystacknews.com

data/examples/argentina.rb

@@ -0,0 +1,82 @@
+require "net/http"
+require "uri"
+require "fileutils"
+require "gd"
+
+ZOOM = 7
+TILE_SIZE = 256
+
+# Argentina bounding box (WGS84)
+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
+
+# Compute tile range for Argentina
+x_min = lon2tile(WEST,  ZOOM)
+x_max = lon2tile(EAST,  ZOOM)
+y_min = lat2tile(NORTH, ZOOM)
+y_max = lat2tile(SOUTH, ZOOM)
+
+puts "Argentina tiles at z=#{ZOOM}:"
+puts "x: #{x_min} .. #{x_max}"
+puts "y: #{y_min} .. #{y_max}"
+
+COLS = x_max - x_min + 1
+ROWS = y_max - y_min + 1
+
+WIDTH  = COLS * TILE_SIZE
+HEIGHT = ROWS * TILE_SIZE
+
+FileUtils.mkdir_p("tmp/tiles")
+FileUtils.mkdir_p("output")
+
+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 "Downloading tiles..."
+
+(x_min..x_max).each do |x|
+  (y_min..y_max).each do |y|
+    fetch_tile(ZOOM, x, y)
+  end
+end
+
+puts "Composing map..."
+
+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("tmp/tiles/#{ZOOM}_#{x}_#{y}.png")
+    map.copy(tile, cx * TILE_SIZE, cy * TILE_SIZE, 0, 0, TILE_SIZE, TILE_SIZE)
+  end
+end
+
+map.save("output/argentina.png")
+
+puts "Saved output/argentina.png"

data/examples/argentina_con_localidades.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
+
+# --------------------------------
+# Cargar mapa base
+# --------------------------------
+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
+
+# --------------------------------
+# Bounding box proyectado
+# --------------------------------
+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
+
+# --------------------------------
+# Cargar localidades
+# --------------------------------
+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
+
+# --------------------------------
+# Guardar resultado
+# --------------------------------
+img.save(OUTPUT)
+puts "Saved #{OUTPUT}"

data/examples/argentina_museos.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_museos.png")
+puts "Saved output/argentina_museos.png"

data/examples/hydro_plants.geojson

@@ -0,0 +1,25 @@
+{
+"type": "FeatureCollection",
+"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:OGC:1.3:CRS84" } },
+"features": [
+{ "type": "Feature", "properties": { "id": 0, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Kagera", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": "T" }, "geometry": { "type": "Point", "coordinates": [ 31.0933382, -1.196415 ] } },
+{ "type": "Feature", "properties": { "id": 7, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Nyumba ya Mungu", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 37.4403441, -3.8195803 ] } },
+{ "type": "Feature", "properties": { "id": 8, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Chutes Pangani I", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": "T" }, "geometry": { "type": "Point", "coordinates": [ 38.759691784504, -5.343275566027487 ] } },
+{ "type": "Feature", "properties": { "id": 9, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Chutes Pangani II", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 38.481225781243062, -5.269153170191239 ] } },
+{ "type": "Feature", "properties": { "id": 10, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "New Pangani", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 38.591509346890959, -5.378960936287418 ] } },
+{ "type": "Feature", "properties": { "id": 11, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Hale Tanesco", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": "T" }, "geometry": { "type": "Point", "coordinates": [ 38.635622773150118, -5.225224606056146 ] } },
+{ "type": "Feature", "properties": { "id": 13, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Mufindi", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": "T" }, "geometry": { "type": "Point", "coordinates": [ 35.5402297, -7.9807668 ] } },
+{ "type": "Feature", "properties": { "id": 14, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Mtera", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 35.891170578822909, -7.096171197430225 ] } },
+{ "type": "Feature", "properties": { "id": 15, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Kidata", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 36.970001220703125, -7.679999828338623 ] } },
+{ "type": "Feature", "properties": { "id": 16, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Bas Kihansi", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": "T" }, "geometry": { "type": "Point", "coordinates": [ 35.8792647, -8.6122779 ] } },
+{ "type": "Feature", "properties": { "id": 19, "status": "Existing", "source": "World Bank", "country": "Tanzania", "name": "Ruhudji", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 0, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 34.779586311911132, -9.462585584742733 ] } },
+{ "type": "Feature", "properties": { "id": 0, "status": "Proposed", "source": "Tanesco", "country": "Tanzania", "name": "Masigira", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 2015, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 35.057617146214611, -9.858152893694784 ] } },
+{ "type": "Feature", "properties": { "id": 0, "status": "Proposed", "source": "Tanesco", "country": "Tanzania", "name": "Ruhudji", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 2015, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 35.432706580089054, -9.2930181466573 ] } },
+{ "type": "Feature", "properties": { "id": 0, "status": "Proposed", "source": "Tanesco", "country": "Tanzania", "name": "Nakatuta", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 2015, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 35.472716119702326, -11.15346173867451 ] } },
+{ "type": "Feature", "properties": { "id": 0, "status": "Existing", "source": "Tanesco", "country": "Tanzania", "name": "Sunda Falls(MHP)", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 2015, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 36.89805596842519, -11.413523746160788 ] } },
+{ "type": "Feature", "properties": { "id": 0, "status": "Existing", "source": "Tanesco", "country": "Tanzania", "name": "Kihansi", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 2015, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 35.860183447113322, -8.140235053135704 ] } },
+{ "type": "Feature", "properties": { "id": 0, "status": "Existing", "source": "Tanesco", "country": "Tanzania", "name": "Tosamaganga", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 2015, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 35.385145143948549, -7.423783841805231 ] } },
+{ "type": "Feature", "properties": { "id": 0, "status": "Proposed", "source": "Tanesco", "country": "Tanzania", "name": "Mpanga", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 2015, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 35.752782896995242, -8.727883399619813 ] } },
+{ "type": "Feature", "properties": { "id": 0, "status": "Proposed", "source": "Tanesco", "country": "Tanzania", "name": "Rumakali", "type": "Hydro", "capacit_MW": 0, "operator": null, "year": 2015, "difference": null }, "geometry": { "type": "Point", "coordinates": [ 33.832324995558118, -8.942934675041157 ] } }
+]
+}