libgd-gis 0.4.0 → 0.4.2

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 74dd1161f2cc225b81f304e90a5e60c2bf249da1606c7a84798367f3d1d1c4bb
-  data.tar.gz: 51a6d88d8f49e6e14a363c68fbd9a0837aae91aef028b456f4df20865831fce8
+  metadata.gz: fa4dd3c210385e6b89998eeccbc88c4e82d5c45b6081be260b98843b4bd01124
+  data.tar.gz: c0e6f00da6d017e7d30f8ccc6c14b4a16ee3a5e7799d5cbf8c4a4b72b137bde6
 SHA512:
-  metadata.gz: 78987b17ccf96eb35607632956a97add64353d9500a5ec9b208028a7064a11501ee5d0b19ce0b1134e41062dd35387b891eb1d67316737a0f7c83279223d1e2c
-  data.tar.gz: e2a9b4e8149839b5b9d7b834757e304979a81e3879a4332f97159e75236be79ea4c4e179a54ed3efd94231c24d290edc83812793b2962ffdbd40ae12549bf40a
+  metadata.gz: fd5c30cbcdae3f214a9f0fb9605edd1244ca161774fc6e60fabcff70688b9b3712bc7444be326fd88c5517d470bfa8c06d8bc0522ca470edc86e5608235a26e9
+  data.tar.gz: 1cfd799c2c1abfb457ad761de7b293a97b5f880ddd55cd4e02298820c91f344977ee176abf8d7da9b77594247a1c8c1a1f34db951f69f15ff00ed7748d3c4da5

data/lib/gd/gis/bbox_resolver.rb

@@ -0,0 +1,27 @@
+module GD
+  module GIS
+    module BBoxResolver
+      def self.resolve(bbox)
+        case bbox
+        when Symbol, String
+          Extents.fetch(bbox)
+
+        when Array
+          validate!(bbox)
+          bbox.map(&:to_f)
+
+        else
+          raise ArgumentError,
+            "bbox must be Symbol, String or [min_lng, min_lat, max_lng, max_lat]"
+        end
+      end
+
+      def self.validate!(bbox)
+        unless bbox.is_a?(Array) && bbox.size == 4
+          raise ArgumentError,
+            "bbox must be [min_lng, min_lat, max_lng, max_lat]"
+        end
+      end
+    end
+  end
+end

data/lib/gd/gis/crs_normalizer.rb

@@ -69,6 +69,9 @@ module GD
           when EPSG3857
             mercator_to_wgs84(lon, lat)
 
+          when "EPSG:22195"
+            gk_to_wgs84(lon, lat)
+
           else
             raise ArgumentError, "Unsupported CRS: #{@crs}"
           end
@@ -97,6 +100,70 @@ module GD
           lat = ((2 * Math.atan(Math.exp(y / r))) - (Math::PI / 2)) * 180.0 / Math::PI
           [lon, lat]
         end
+
+        # Converts Gauss–Krüger (GK) projected coordinates to WGS84.
+        #
+        # This method converts easting/northing coordinates from
+        # Gauss–Krüger Argentina Zone 5 (EPSG:22195) into
+        # WGS84 longitude/latitude (degrees).
+        #
+        # The implementation is intended for cartographic rendering
+        # and visualization purposes, not for high-precision geodesy.
+        #
+        # @param easting [Numeric]
+        #   Easting value in meters.
+        #
+        # @param northing [Numeric]
+        #   Northing value in meters.
+        #
+        # @return [Array<Float>]
+        #   A `[longitude, latitude]` pair in decimal degrees (WGS84).
+        #
+        # @example Convert Gauss–Krüger coordinates
+        #   gk_to_wgs84(580_000, 6_176_000)
+        #   # => [longitude, latitude]
+        #
+        # @note
+        #   This method assumes:
+        #   - Central meridian: −60°
+        #   - False easting: 500,000 m
+        #   - WGS84-compatible ellipsoid
+        #
+        # @see https://epsg.io/22195
+
+        def gk_to_wgs84(easting, northing)
+          a = 6378137.0
+          f = 1 / 298.257223563
+          e2 = (2 * f) - (f * f)
+          lon0 = -60.0 * Math::PI / 180.0
+
+          x = easting - 500_000.0
+          y = northing - 10_000_000.0
+
+          m = y
+          mu = m / (a * (1 - (e2 / 4) - (3 * e2 * e2 / 64)))
+
+          e1 = (1 - Math.sqrt(1 - e2)) / (1 + Math.sqrt(1 - e2))
+
+          j1 = (3 * e1 / 2) - (27 * (e1**3) / 32)
+          j2 = (21 * (e1**2) / 16) - (55 * (e1**4) / 32)
+
+          fp = mu + (j1 * Math.sin(2 * mu)) + (j2 * Math.sin(4 * mu))
+
+          c1 = e2 * (Math.cos(fp)**2)
+          t1 = Math.tan(fp)**2
+          r1 = a * (1 - e2) / ((1 - (e2 * (Math.sin(fp)**2)))**1.5)
+          n1 = a / Math.sqrt(1 - (e2 * (Math.sin(fp)**2)))
+
+          d = x / n1
+
+          lat = fp - ((n1 * Math.tan(fp) / r1) *
+                    (((d**2) / 2) - ((5 + (3 * t1) + (10 * c1)) * (d**4) / 24)))
+
+          lon = lon0 + ((d - ((1 + (2 * t1) + c1) * (d**3) / 6)) / Math.cos(fp))
+
+          [lon * 180.0 / Math::PI, lat * 180.0 / Math::PI]
+        end
       end
     end
   end

data/lib/gd/gis/data/extents_global.json

@@ -0,0 +1,95 @@
+{
+  "world": [-180.0, -90.0, 180.0, 90.0],
+
+  "northern_hemisphere": [-180.0, 0.0, 180.0, 90.0],
+  "southern_hemisphere": [-180.0, -90.0, 180.0, 0.0],
+  "eastern_hemisphere": [0.0, -90.0, 180.0, 90.0],
+  "western_hemisphere": [-180.0, -90.0, 0.0, 90.0],
+
+  "africa": [-25.0, -35.0, 51.0, 38.0],
+  "antarctica": [-180.0, -90.0, 180.0, -60.0],
+  "asia": [26.0, -11.0, 169.0, 78.0],
+  "europe": [-31.0, 34.0, 39.0, 72.0],
+  "north_america": [-168.0, 7.0, -52.0, 83.0],
+  "south_america": [-82.0, -56.0, -34.0, 13.0],
+  "oceania": [110.0, -50.0, 180.0, 0.0],
+
+  "argentina": [-73.6, -55.1, -53.6, -21.7],
+  "bolivia": [-69.7, -22.9, -57.5, -9.7],
+  "brazil": [-74.0, -34.0, -34.0, 5.3],
+  "chile": [-75.6, -56.0, -66.4, -17.5],
+  "colombia": [-79.0, -4.2, -66.9, 13.5],
+  "ecuador": [-81.0, -5.0, -75.0, 1.5],
+  "guyana": [-61.4, 1.2, -56.5, 8.6],
+  "paraguay": [-62.6, -27.6, -54.3, -19.3],
+  "peru": [-81.4, -18.4, -68.7, -0.0],
+  "suriname": [-58.1, 1.8, -53.9, 6.0],
+  "uruguay": [-58.4, -35.0, -53.1, -30.1],
+  "venezuela": [-73.4, 0.6, -59.8, 12.2],
+
+  "canada": [-141.0, 41.0, -52.6, 83.1],
+  "united_states": [-125.0, 24.0, -66.9, 49.4],
+  "mexico": [-118.4, 14.5, -86.7, 32.7],
+  "greenland": [-73.0, 59.8, -12.2, 83.6],
+
+  "iceland": [-24.5, 63.3, -13.5, 66.6],
+  "ireland": [-10.5, 51.4, -6.0, 55.5],
+  "united_kingdom": [-8.6, 49.8, 1.8, 60.9],
+  "portugal": [-9.6, 36.9, -6.2, 42.2],
+  "spain": [-9.5, 35.9, 3.3, 43.8],
+  "france": [-5.2, 41.3, 9.6, 51.1],
+  "belgium": [2.5, 49.5, 6.4, 51.5],
+  "netherlands": [3.3, 50.7, 7.2, 53.7],
+  "germany": [5.9, 47.3, 15.0, 55.1],
+  "switzerland": [5.9, 45.8, 10.5, 47.8],
+  "italy": [6.6, 36.6, 18.5, 47.1],
+  "austria": [9.5, 46.4, 17.2, 49.0],
+  "poland": [14.1, 49.0, 24.1, 54.8],
+  "czechia": [12.1, 48.5, 18.9, 51.1],
+  "hungary": [16.1, 45.7, 22.9, 48.6],
+  "greece": [19.4, 34.8, 28.2, 41.8],
+  "norway": [4.9, 58.0, 31.3, 71.2],
+  "sweden": [11.1, 55.3, 24.2, 69.1],
+  "finland": [20.6, 59.8, 31.6, 70.1],
+  "denmark": [8.0, 54.5, 12.7, 57.8],
+
+  "turkey": [26.0, 36.0, 45.0, 42.1],
+  "russia": [19.0, 41.0, 180.0, 82.0],
+
+  "saudi_arabia": [34.5, 16.3, 55.7, 32.2],
+  "iran": [44.0, 25.0, 63.3, 39.8],
+  "iraq": [38.8, 29.0, 48.6, 37.4],
+  "israel": [34.2, 29.5, 35.9, 33.3],
+  "jordan": [34.9, 29.2, 39.3, 33.4],
+  "syria": [35.7, 32.3, 42.4, 37.3],
+
+  "india": [68.1, 6.5, 97.4, 35.5],
+  "pakistan": [60.9, 23.7, 77.8, 37.1],
+  "bangladesh": [88.0, 20.7, 92.7, 26.6],
+  "china": [73.5, 18.1, 135.1, 53.6],
+  "mongolia": [87.7, 41.6, 119.9, 52.1],
+  "japan": [129.4, 31.0, 145.8, 45.5],
+  "south_korea": [125.0, 33.1, 131.9, 38.6],
+  "north_korea": [124.0, 37.7, 130.7, 43.0],
+
+  "thailand": [97.3, 5.6, 105.6, 20.5],
+  "vietnam": [102.1, 8.2, 109.5, 23.4],
+  "malaysia": [99.6, 0.8, 119.3, 7.4],
+  "indonesia": [95.0, -11.0, 141.0, 6.0],
+  "philippines": [116.9, 4.6, 126.6, 21.2],
+
+  "egypt": [24.7, 22.0, 36.9, 31.7],
+  "libya": [9.3, 19.5, 25.1, 33.2],
+  "algeria": [-8.7, 19.0, 12.0, 37.1],
+  "morocco": [-13.2, 27.7, -1.0, 35.9],
+  "tunisia": [7.5, 30.2, 11.6, 37.3],
+
+  "nigeria": [2.7, 4.2, 14.7, 13.9],
+  "ethiopia": [32.9, 3.4, 47.9, 14.9],
+  "kenya": [33.9, -4.7, 41.9, 5.0],
+  "south_africa": [16.4, -34.8, 32.9, -22.1],
+
+  "australia": [112.9, -43.7, 153.6, -10.7],
+  "new_zealand": [166.4, -47.3, 178.6, -34.4],
+  "papua_new_guinea": [140.8, -11.7, 156.0, -1.4]
+}

data/lib/gd/gis/extents.rb

@@ -0,0 +1,39 @@
+require "json"
+
+module GD
+  module GIS
+    module Extents
+      DATA_PATH = File.expand_path(
+        "data/extents_global.json",
+        __dir__
+      )
+
+      @extents = nil
+
+      class << self
+        def fetch(name)
+          load_data!
+          @extents.fetch(name.to_s.downcase) do
+            raise ArgumentError, "Unknown extent: #{name}"
+          end
+        end
+
+        def [](name)
+          fetch(name)
+        end
+
+        def all
+          load_data!
+          @extents.keys.map(&:to_sym)
+        end
+
+        private
+
+        def load_data!
+          return if @extents
+          @extents = JSON.parse(File.read(DATA_PATH))
+        end
+      end
+    end
+  end
+end

data/lib/gd/gis/layer_points.rb

@@ -101,16 +101,15 @@ module GD
       #
       # @return [void]
       def render!(img, projector)
-        value = case @icon
-                when "numeric"
-                  @symbol
-                when "alphabetic"
-                  (@symbol + 96).chr
-                when "symbol"
-                  @symbol
-                else
-                  @icon
-                end
+        value =
+          case @icon
+          when "numeric", "symbol"
+            @symbol
+          when "alphabetic"
+            (@symbol + 96).chr
+          else
+            @icon
+          end
 
         if @icon.is_a?(GD::Image)
           w = @icon.width
@@ -176,7 +175,7 @@ module GD
         # baseline  = top_y + h = y + h/2
         text_x = (x - (w / 2.0)).round
         text_y = (y + (h / 2.0)).round
-        
+
         # 4) Draw number
         img.text(
           text,

data/lib/gd/gis/legend.rb

@@ -0,0 +1,38 @@
+# lib/gd/gis/legend.rb
+
+module GD
+  module GIS
+    LegendItem = Struct.new(:color, :label)
+
+    # Represents a map legend rendered as part of the final image.
+    #
+    # A Legend provides visual context for map elements by associating
+    # colors or symbols with human-readable labels.
+    #
+    # Legends are rendered server-side and embedded directly into the
+    # resulting map image, allowing the map to be self-explanatory
+    # without relying on external UI components.
+    #
+    # A Legend is typically created and configured via {Map#legend}
+    # and rendered automatically during the map rendering pipeline.
+    #
+    # @example Creating a legend
+    #   map.legend do |l|
+    #     l.add [76, 175, 80, 0], "Delivered"
+    #     l.add [255, 193, 7, 0], "In transit"
+    #     l.add [244, 67, 54, 0], "Delayed"
+    #
+    class Legend
+      attr_reader :items, :position
+
+      def initialize(position: :bottom_right)
+        @position = position
+        @items = []
+      end
+
+      def add(color, label)
+        @items << LegendItem.new(color, label)
+      end
+    end
+  end
+end

data/lib/gd/gis/map.rb

@@ -7,6 +7,7 @@ require_relative "layer_geojson"
 require_relative "layer_points"
 require_relative "layer_lines"
 require_relative "layer_polygons"
+require_relative "legend"
 
 LINE_GEOMS = %w[LineString MultiLineString].freeze
 POLY_GEOMS = %w[Polygon MultiPolygon].freeze
@@ -69,6 +70,9 @@ module GD
         crs: nil,
         fitted_bbox: false
       )
+        # resolve symbolic bbox
+        bbox = GD::GIS::BBoxResolver.resolve(bbox)
+
         # 1. Basic input validation
         raise ArgumentError, "bbox must be [min_lng, min_lat, max_lng, max_lat]" unless
           bbox.is_a?(Array) && bbox.size == 4
@@ -217,6 +221,110 @@ module GD
         @used_labels[key] = true
       end
 
+      def legend(position: :bottom_right)
+        @legend = Legend.new(position: position)
+        yield @legend
+      end
+
+      def legend_from_layers(position: :bottom_right)
+        @legend = Legend.new(position: position)
+
+        layers.each do |layer|
+          next unless layer.respond_to?(:color)
+
+          @legend.add(layer.color, layer.name)
+        end
+      end
+
+      def draw_legend
+        return unless @legend
+        return unless @image
+        return unless @style
+        return unless @style.global
+        return if @style.global[:label] == false
+
+        label_style = @style.global[:label] || {}
+
+        padding     = 10
+        box_size    = 12
+        line_height = 18
+        margin      = 15
+
+        # --- font (from style) -----------------------------------
+
+        font_path =
+          case label_style[:font]
+          when nil, "default"
+            GD::GIS::FontHelper.random
+          else
+            label_style[:font]
+          end
+
+        font_size  = label_style[:size] || 10
+        font_color = GD::Color.rgba(*(label_style[:color] || [0, 0, 0, 0]))
+
+        # --- measure text (CORRECT API) ---------------------------
+
+        text_widths = @legend.items.map do |i|
+          w, = @image.text_bbox(
+            i.label,
+            font: font_path,
+            size: font_size
+          )
+          w
+        end
+
+        width  = (text_widths.max || 0) + box_size + (padding * 3)
+        height = (@legend.items.size * line_height) + (padding * 2)
+
+        # --- position --------------------------------------------
+
+        x, y =
+          case @legend.position
+          when :bottom_right
+            [@image.width - width - margin, @image.height - height - margin]
+          when :bottom_left
+            [margin, @image.height - height - margin]
+          when :top_right
+            [@image.width - width - margin, margin]
+          else
+            [margin, margin]
+          end
+
+        # --- background ------------------------------------------
+
+        bg     = GD::Color.rgba(255, 255, 255, 80)
+        border = GD::Color.rgb(200, 200, 200)
+
+        @image.filled_rectangle(x, y, x + width, y + height, bg)
+        @image.rectangle(x, y, x + width, y + height, border)
+
+        # --- items -----------------------------------------------
+
+        @legend.items.each_with_index do |item, idx|
+          iy = y + padding + (idx * line_height)
+
+          # color box
+          @image.filled_rectangle(
+            x + padding,
+            iy,
+            x + padding + box_size,
+            iy + box_size,
+            GD::Color.rgba(*item.color)
+          )
+
+          # label text
+          @image.text_ft(
+            item.label,
+            x: x + padding + box_size + 8,
+            y: iy + box_size,
+            font: font_path,
+            size: font_size,
+            color: font_color
+          )
+        end
+      end
+
       # Loads features from a GeoJSON file.
       #
       # This method:
@@ -376,10 +484,10 @@ module GD
 
         @points_layers << GD::GIS::PointsLayer.new(
           [row],
-          lon: -> r { r[:lon] },
-          lat: -> r { r[:lat] },
+          lon: ->(r) { r[:lon] },
+          lat: ->(r) { r[:lat] },
           icon: icon || @style.point[:icon],
-          label: label ? -> r { r[:label] } : nil,
+          label: label ? ->(r) { r[:label] } : nil,
           font: font || @style.point[:font],
           size: size || @style.point[:size],
           color: color || @style.point[:color],
@@ -477,6 +585,8 @@ module GD
         @polygons_layers.each { |l| l.render!(@image, projection) }
         @lines_layers.each    { |l| l.render!(@image, projection) }
         @points_layers.each   { |l| l.render!(@image, projection) }
+
+        draw_legend
       end
 
       def render_viewport
@@ -536,6 +646,8 @@ module GD
         @polygons_layers.each { |l| l.render!(@image, projection) }
         @lines_layers.each    { |l| l.render!(@image, projection) }
         @points_layers.each   { |l| l.render!(@image, projection) }
+
+        draw_legend
       end
 
       # Saves the rendered image to disk.

data/lib/gd/gis.rb

@@ -32,6 +32,9 @@ require_relative "gis/font_helper"
 require_relative "gis/style"
 require_relative "gis/classifier"
 
+require_relative "gis/extents"
+require_relative "gis/bbox_resolver"
+
 require_relative "gis/feature"
 require_relative "gis/map"
 require_relative "gis/basemap"

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: libgd-gis
 version: !ruby/object:Gem::Version
-  version: 0.4.0
+  version: 0.4.2
 platform: ruby
 authors:
 - Germán Alberto Giménez Silva
@@ -40,9 +40,12 @@ files:
 - README.md
 - lib/gd/gis.rb
 - lib/gd/gis/basemap.rb
+- lib/gd/gis/bbox_resolver.rb
 - lib/gd/gis/classifier.rb
 - lib/gd/gis/color_helpers.rb
 - lib/gd/gis/crs_normalizer.rb
+- lib/gd/gis/data/extents_global.json
+- lib/gd/gis/extents.rb
 - lib/gd/gis/feature.rb
 - lib/gd/gis/font_helper.rb
 - lib/gd/gis/geometry.rb
@@ -50,8 +53,8 @@ files:
 - lib/gd/gis/layer_lines.rb
 - lib/gd/gis/layer_points.rb
 - lib/gd/gis/layer_polygons.rb
+- lib/gd/gis/legend.rb
 - lib/gd/gis/map.rb
-- lib/gd/gis/middleware.rb
 - lib/gd/gis/ontology.rb
 - lib/gd/gis/ontology.yml
 - lib/gd/gis/projection.rb
@@ -76,7 +79,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 4.0.5
+rubygems_version: 4.0.6
 specification_version: 4
 summary: Geospatial raster rendering for Ruby using libgd
 test_files: []

data/lib/gd/gis/middleware.rb

@@ -1,152 +0,0 @@
-# frozen_string_literal: true
-
-module GD
-  module GIS
-    # Coordinate Reference System (CRS) helpers and normalizers.
-    #
-    # This module defines commonly used CRS identifiers and
-    # utilities for normalizing coordinates into a single,
-    # consistent representation.
-    #
-    module CRS
-      # OGC CRS84 (longitude, latitude)
-      CRS84      = "urn:ogc:def:crs:OGC:1.3:CRS84"
-
-      # EPSG:4326 (latitude, longitude axis order)
-      EPSG4326   = "EPSG:4326"
-
-      # EPSG:3857 (Web Mercator)
-      EPSG3857   = "EPSG:3857"
-
-      # Gauss–Krüger Argentina, zone 5
-      #
-      # Note: This constant represents a *specific* GK zone
-      # and is not a generic Gauss–Krüger definition.
-      GK_ARGENTINA = "EPSG:22195"
-
-      # Normalizes coordinates from supported CRS definitions
-      # into CRS84 (longitude, latitude in degrees).
-      #
-      # Supported input CRS:
-      # - CRS84
-      # - EPSG:4326 (axis order normalization)
-      # - EPSG:3857 (Web Mercator)
-      # - EPSG:22195 (Gauss–Krüger Argentina, zone 5)
-      #
-      # All outputs are returned as:
-      #   [longitude, latitude] in degrees
-      #
-      # ⚠️ Projection conversions are intended for mapping
-      # and visualization, not for high-precision geodesy.
-      #
-      class Normalizer
-        # Creates a new CRS normalizer.
-        #
-        # @param crs [String, Symbol, nil]
-        #   CRS identifier; defaults to CRS84 if nil
-        def initialize(crs)
-          @crs = normalize_name(crs)
-        end
-
-        # Normalizes a coordinate pair into CRS84.
-        #
-        # @param lon [Numeric]
-        #   first coordinate (meaning depends on input CRS)
-        # @param lat [Numeric]
-        #   second coordinate (meaning depends on input CRS)
-        #
-        # @return [Array<Float>]
-        #   normalized [longitude, latitude] in degrees
-        #
-        # @raise [RuntimeError]
-        #   if the CRS is not supported
-        def normalize(lon, lat)
-          case @crs
-          when CRS84
-            [lon, lat]
-
-          when EPSG4326
-            # EPSG:4326 uses (lat, lon)
-            [lat, lon]
-
-          when GK_ARGENTINA
-            gk_to_wgs84(lon, lat)
-
-          when EPSG3857
-            mercator_to_wgs84(lon, lat)
-
-          else
-            raise "Unsupported CRS: #{@crs}"
-          end
-        end
-
-        private
-
-        # Normalizes a CRS name into a comparable string.
-        #
-        # @param name [Object]
-        # @return [String]
-        def normalize_name(name)
-          return CRS84 if name.nil?
-
-          name.to_s.strip
-        end
-
-        # Converts Web Mercator coordinates to WGS84.
-        #
-        # @param x [Numeric] X coordinate in meters
-        # @param y [Numeric] Y coordinate in meters
-        # @return [Array<Float>] [longitude, latitude] in degrees
-        def mercator_to_wgs84(x, y)
-          r = 6378137.0
-          lon = (x / r) * 180.0 / Math::PI
-          lat = ((2 * Math.atan(Math.exp(y / r))) - (Math::PI / 2)) * 180.0 / Math::PI
-          [lon, lat]
-        end
-
-        # Converts Gauss–Krüger Argentina (zone 5) coordinates to WGS84.
-        #
-        # This implementation provides sufficient accuracy for
-        # cartographic rendering and visualization.
-        #
-        # @param easting [Numeric] easting (meters)
-        # @param northing [Numeric] northing (meters)
-        # @return [Array<Float>] [longitude, latitude] in degrees
-        def gk_to_wgs84(easting, northing)
-          # Parameters for Argentina GK Zone 5
-          a = 6378137.0
-          f = 1 / 298.257223563
-          e2 = (2 * f) - (f * f)
-          lon0 = -60.0 * Math::PI / 180.0 # central meridian zone 5
-
-          x = easting - 500000.0
-          y = northing
-
-          m = y
-          mu = m / (a * (1 - (e2 / 4) - (3 * e2 * e2 / 64)))
-
-          e1 = (1 - Math.sqrt(1 - e2)) / (1 + Math.sqrt(1 - e2))
-
-          j1 = (3 * e1 / 2) - (27 * (e1**3) / 32)
-          j2 = (21 * (e1**2) / 16) - (55 * (e1**4) / 32)
-
-          fp = mu + (j1 * Math.sin(2 * mu)) + (j2 * Math.sin(4 * mu))
-
-          c1 = e2 * (Math.cos(fp)**2)
-          t1 = Math.tan(fp)**2
-          r1 = a * (1 - e2) / ((1 - (e2 * (Math.sin(fp)**2)))**1.5)
-          n1 = a / Math.sqrt(1 - (e2 * (Math.sin(fp)**2)))
-
-          d = x / n1
-
-          lat = fp - ((n1 * Math.tan(fp) / r1) *
-                     (((d**2) / 2) - ((5 + (3 * t1) + (10 * c1)) * (d**4) / 24)))
-
-          lon = lon0 + ((d - ((1 + (2 * t1) + c1) * (d**3) / 6)) / Math.cos(fp))
-
-          [lon * 180.0 / Math::PI, lat * 180.0 / Math::PI]
-        end
-      end
-    end
-  end
-end