libgd-gis 0.1.2 → 0.2.0

data/examples/a-world_peaks.rb

@@ -1,42 +0,0 @@
-require "json"
-require "gd/gis"
-
-# ---------------------------
-# Bounding box mundial
-# ---------------------------
-WORLD = [-170, -60, -30, 75]
-
-# ---------------------------
-# Crear mapa
-# ---------------------------
-map = GD::GIS::Map.new(
-  bbox: WORLD,
-  zoom: 4,
-  basemap: :carto_light
-)
-
-# Cargar datos
-# ---------------------------
-peaks = JSON.parse(File.read("picks.json"))
-
-# ---------------------------
-# Agregar capa de puntos
-# ---------------------------
-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]
-)
-
-# ---------------------------
-# Renderizar y guardar
-# ---------------------------
-map.render
-map.save("output/world_peaks-s.png")
-
-puts "Saved output/world_peaks.png"

data/examples/argentina_con_localidades.rb

@@ -1,89 +0,0 @@
-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

@@ -1,70 +0,0 @@
-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

@@ -1,25 +0,0 @@
-{
-"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 ] } }
-]
-}

data/examples/montanas.rb

@@ -1,118 +0,0 @@
-require "net/http"
-require "uri"
-require "fileutils"
-require "json"
-require "gd"
-
-# =========================
-# Parámetros
-# =========================
-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")
-
-# =========================
-# Web Mercator helpers
-# =========================
-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
-
-# =========================
-# Descargar tile
-# =========================
-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
-
-# =========================
-# Construir mapa base
-# =========================
-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
-
-# =========================
-# Cargar picos
-# =========================
-peaks = JSON.parse(File.read("picks.json"))
-
-red   = GD::Color.rgb(220, 40, 40)
-black = GD::Color.rgb(0, 0, 0)
-
-# =========================
-# Dibujar picos
-# =========================
-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
-
-  # Marcador
-  map.filled_circle(x, y, 7, red)
-
-  # Texto con FreeType
-  map.text(
-    name,
-    x: x + 10,
-    y: y - 5,
-    size: 10,
-    color: black,
-    font: FONT
-  )
-end
-
-# =========================
-# Guardar
-# =========================
-map.save("output/world_peaks.png")
-puts "Saved output/world_peaks.png"

data/examples/picks.json

@@ -1,50 +0,0 @@
-[
-  {"name": "Mount Everest", "elevation_m": 8849, "latitude": 27.9881, "longitude": 86.9250, "country": "Nepal/China"},
-  {"name": "K2", "elevation_m": 8611, "latitude": 35.8825, "longitude": 76.5133, "country": "Pakistan/China"},
-  {"name": "Kangchenjunga", "elevation_m": 8586, "latitude": 27.7025, "longitude": 88.1475, "country": "Nepal/India"},
-  {"name": "Lhotse", "elevation_m": 8516, "latitude": 27.9617, "longitude": 86.9333, "country": "Nepal/China"},
-  {"name": "Makalu", "elevation_m": 8485, "latitude": 27.8867, "longitude": 87.0867, "country": "Nepal/China"},
-  {"name": "Cho Oyu", "elevation_m": 8188, "latitude": 28.0942, "longitude": 86.6608, "country": "Nepal/China"},
-  {"name": "Dhaulagiri I", "elevation_m": 8167, "latitude": 28.6967, "longitude": 83.4867, "country": "Nepal"},
-  {"name": "Manaslu", "elevation_m": 8163, "latitude": 28.5497, "longitude": 84.5617, "country": "Nepal"},
-  {"name": "Nanga Parbat", "elevation_m": 8126, "latitude": 35.2375, "longitude": 74.5892, "country": "Pakistan"},
-  {"name": "Annapurna I", "elevation_m": 8091, "latitude": 28.5961, "longitude": 83.8203, "country": "Nepal"},
-  {"name": "Gasherbrum I", "elevation_m": 8080, "latitude": 35.7211, "longitude": 76.6764, "country": "Pakistan/China"},
-  {"name": "Broad Peak", "elevation_m": 8051, "latitude": 35.8128, "longitude": 76.5647, "country": "Pakistan/China"},
-  {"name": "Gasherbrum II", "elevation_m": 8035, "latitude": 35.7533, "longitude": 76.6333, "country": "Pakistan/China"},
-  {"name": "Shishapangma", "elevation_m": 8027, "latitude": 28.3486, "longitude": 85.7781, "country": "China"},
-  {"name": "Aconcagua", "elevation_m": 6961, "latitude": -32.6532, "longitude": -70.0109, "country": "Argentina"},
-  {"name": "Ojos del Salado", "elevation_m": 6893, "latitude": -27.1192, "longitude": -68.5422, "country": "Chile/Argentina"},
-  {"name": "Monte Pissis", "elevation_m": 6795, "latitude": -27.7500, "longitude": -68.7700, "country": "Argentina"},
-  {"name": "Huascarán", "elevation_m": 6768, "latitude": -9.1211, "longitude": -77.6153, "country": "Peru"},
-  {"name": "Cerro Bonete", "elevation_m": 6759, "latitude": -28.5400, "longitude": -67.7300, "country": "Argentina"},
-  {"name": "Volcán Llullaillaco", "elevation_m": 6739, "latitude": -24.7214, "longitude": -68.5328, "country": "Chile/Argentina"},
-  {"name": "Chimborazo", "elevation_m": 6263, "latitude": -1.4690, "longitude": -78.8170, "country": "Ecuador"},
-  {"name": "Denali", "elevation_m": 6190, "latitude": 63.0692, "longitude": -151.0070, "country": "United States"},
-  {"name": "Mount Logan", "elevation_m": 5959, "latitude": 60.5639, "longitude": -140.3964, "country": "Canada"},
-  {"name": "Pico de Orizaba", "elevation_m": 5636, "latitude": 19.0300, "longitude": -97.2667, "country": "Mexico"},
-  {"name": "Sierra Nevada de Santa Marta", "elevation_m": 5700, "latitude": 10.8500, "longitude": -73.5000, "country": "Colombia"},
-  {"name": "Volcán Nevado del Huila", "elevation_m": 5364, "latitude": 2.9428, "longitude": -76.0267, "country": "Colombia"},
-  {"name": "Volcán Tajumulco", "elevation_m": 4203, "latitude": 15.0333, "longitude": -91.9000, "country": "Guatemala"},
-  {"name": "Pico Duarte", "elevation_m": 3098, "latitude": 19.0833, "longitude": -70.7667, "country": "Dominican Republic"},
-  {"name": "Cerro El Pital", "elevation_m": 2730, "latitude": 14.2956, "longitude": -88.3075, "country": "El Salvador/Honduras"},
-  {"name": "Mount Vinson", "elevation_m": 4892, "latitude": -78.5250, "longitude": -85.6167, "country": "Antarctica"},
-  {"name": "Mount Kilimanjaro", "elevation_m": 5895, "latitude": -3.0674, "longitude": 37.3556, "country": "Tanzania"},
-  {"name": "Mount Kenya", "elevation_m": 5199, "latitude": -0.1517, "longitude": 37.3081, "country": "Kenya"},
-  {"name": "Mount Stanley", "elevation_m": 5109, "latitude": 0.3322, "longitude": 29.8739, "country": "Uganda/Democratic Republic of the Congo"},
-  {"name": "Mount Elbrus", "elevation_m": 5642, "latitude": 43.3547, "longitude": 42.4483, "country": "Russia"},
-  {"name": "Mont Blanc", "elevation_m": 4810, "latitude": 45.8326, "longitude": 6.8652, "country": "France/Italy"},
-  {"name": "Mount Etna", "elevation_m": 3357, "latitude": 37.7510, "longitude": 14.9937, "country": "Italy"},
-  {"name": "Mount Olympus", "elevation_m": 2917, "latitude": 40.0833, "longitude": 22.3333, "country": "Greece"},
-  {"name": "Mount Ararat", "elevation_m": 5137, "latitude": 39.7022, "longitude": 44.2931, "country": "Turkey"},
-  {"name": "Mount Damavand", "elevation_m": 5610, "latitude": 35.9522, "longitude": 52.1083, "country": "Iran"},
-  {"name": "Jabal An-Nabi Shu'ayb", "elevation_m": 3666, "latitude": 15.4900, "longitude": 44.0583, "country": "Yemen"},
-  {"name": "Jabal al-Lawz", "elevation_m": 2543, "latitude": 29.5833, "longitude": 38.6333, "country": "Saudi Arabia"},
-  {"name": "Mount Fuji", "elevation_m": 3776, "latitude": 35.3606, "longitude": 138.7274, "country": "Japan"},
-  {"name": "Mount Apo", "elevation_m": 2954, "latitude": 7.1031, "longitude": 125.2333, "country": "Philippines"},
-  {"name": "Gunung Puncak Jaya", "elevation_m": 4884, "latitude": -4.0825, "longitude": 137.1517, "country": "Indonesia"},
-  {"name": "Gunung Rinjani", "elevation_m": 3726, "latitude": -8.4083, "longitude": 116.4736, "country": "Indonesia"},
-  {"name": "Mount Kosciuszko", "elevation_m": 2228, "latitude": -36.4567, "longitude": 148.2622, "country": "Australia"},
-  {"name": "Mount Wilhelm", "elevation_m": 4509, "latitude": -5.7900, "longitude": 144.8700, "country": "Papua New Guinea"},
-  {"name": "Mount Cook (Aoraki)", "elevation_m": 3724, "latitude": -43.5942, "longitude": 170.1419, "country": "New Zealand"}
-]

data/examples/render_tile.rb

@@ -1,31 +0,0 @@
-require "net/http"
-require "fileutils"
-require_relative "../lib/libgd_gis"
-
-def render
-  tmp = Tempfile.new(["tile", ".png"])
-  tmp.binmode
-
-  uri = URI(@source)
-  req = Net::HTTP::Get.new(uri)
-  req["User-Agent"] = "ruby-libgd-gis/0.1 (https://github.com/ggerman/ruby-libgd-gis)"
-
-  Net::HTTP.start(uri.host, uri.port, use_ssl: true) do |http|
-    http.request(req) do |res|
-      res.read_body { |chunk| tmp.write(chunk) }
-    end
-  end
-
-  tmp.flush
-  @image = GD::Image.open(tmp.path)
-ensure
-  tmp.close
-end
-
-FileUtils.mkdir_p("output")
-
-tile = LibGD::GIS::Tile.osm(z: 4, x: 8, y: 5)
-tile.render
-tile.save("output/tile.png")
-
-puts "Tile written to output/tile.png"