ffi-gdal 1.0.0.beta16 → 1.0.0

data/lib/gdal/{gridder.rb → extensions/gridder.rb}

@@ -1,11 +1,11 @@
 # frozen_string_literal: true
 
 require 'narray'
-require_relative '../gdal'
+require 'gdal'
+require 'gdal/options'
+require 'ogr'
 require_relative 'gridder_options'
 require_relative 'gridder/point_extracting'
-require_relative 'options'
-require_relative '../ogr'
 
 module GDAL
   # Somewhat analogous to the gdal_grid utility.
@@ -46,9 +46,9 @@ module GDAL
     # points, and writes out the newly gridder raster.
     def grid!
       dataset = build_dataset(@options.output_driver,
-        @dest_file_name,
-        @options.output_data_type,
-        @options.output_creation_options)
+                              @dest_file_name,
+                              @options.output_data_type,
+                              @options.output_creation_options)
 
       grid_and_write(dataset.raster_band(1), dataset.geo_transform)
 
@@ -74,10 +74,10 @@ module GDAL
     # @return [GDAL::Dataset]
     def build_dataset(driver, dest_file_name, data_type, creation_options = {})
       dataset = driver.create_dataset(dest_file_name,
-        output_width,
-        output_height,
-        data_type: data_type,
-        **creation_options)
+                                      output_width,
+                                      output_height,
+                                      data_type: data_type,
+                                      **creation_options)
 
       dataset.projection = build_output_spatial_reference
       dataset.geo_transform = build_output_geo_transform
@@ -146,12 +146,12 @@ module GDAL
       @y_max ||= @options.output_y_extent.fetch(:max) { @source_layer.extent.y_max }
     end
 
-    # @return [Fixnum]
+    # @return [Integer]
     def output_width
       @options.output_size[:width]
     end
 
-    # @return [Fixnum]
+    # @return [Integer]
     def output_height
       @options.output_size[:height]
     end
@@ -159,7 +159,7 @@ module GDAL
     # Figures out the proper block sizes to use for iterating over layer pixels,
     # gridding them, and writing them to the raster file.
     #
-    # @param raster_band_block_size [Fixnum]
+    # @param raster_band_block_size [Integer]
     def each_block(raster_band_block_size)
       data_type_size = @options.output_data_type_size
       block_size = build_block_sizes(raster_band_block_size, data_type_size)
@@ -202,7 +202,7 @@ module GDAL
         extents = { x_min: grid_x_min, x_max: grid_x_max, y_min: grid_y_min, y_max: grid_y_max }
 
         @options.grid.create(points, extents, data_ptr, output_size,
-          progress_arg, scaled_progress_ptr)
+                             progress_arg, scaled_progress_ptr)
 
         raster_band.raster_io('w', data_ptr, x_offset: x_offset, y_offset: y_offset,
                                              x_size: x_request, y_size: y_request,
@@ -210,9 +210,9 @@ module GDAL
       end
     end
 
-    # @param raster_band_block_size [Hash{x: Fixnum, y: Fixnum}]
-    # @param data_type_size [Fixnum]
-    # @return [Hash{x: Fixnum, y: Fixnum}]
+    # @param raster_band_block_size [Hash{x: Integer, y: Integer}]
+    # @param data_type_size [Integer]
+    # @return [Hash{x: Integer, y: Integer}]
     def build_block_sizes(raster_band_block_size, data_type_size)
       block_x_size = raster_band_block_size[:x]
       block_y_size = raster_band_block_size[:y]
@@ -241,8 +241,8 @@ module GDAL
     # blocks.
     #
     # @see http://gdal.sourcearchive.com/documentation/1.7.2/gdal_8h_904fbbb050e16c9d0ac028dc5113ef27.html
-    # @param block_number [Fixnum]
-    # @param block_count [Fixnum]
+    # @param block_number [Integer]
+    # @param block_count [Integer]
     # @return [FFI::Pointer]
     def build_scaled_progress_pointer(block_number, block_count)
       return unless @options.progress_formatter
@@ -257,9 +257,9 @@ module GDAL
 
     # Determines how large of a chunk of data to grid and rasterize.
     #
-    # @param block_size [Fixnum]
-    # @param raster_border [Fixnum]
-    # @return [Fixnum]
+    # @param block_size [Integer]
+    # @param raster_border [Integer]
+    # @return [Integer]
     def build_data_request_size(block_size, offset, raster_border)
       request = block_size
 
@@ -273,28 +273,28 @@ module GDAL
     # @param min [Float]
     # @param pixel_size [Float]
     # @param offset [Float]
-    # @param request_size [Fixnum]
-    # @return [Array<Fixnum, Fixnum>] The min and max values based on the given
+    # @param request_size [Integer]
+    # @return [Array<Integer, Integer>] The min and max values based on the given
     #   parameters.
     def build_grid_extents(min, pixel_size, offset, request_size)
-      grid_min = min + pixel_size * offset
-      grid_max = min + pixel_size * (offset + request_size)
+      grid_min = min + (pixel_size * offset)
+      grid_max = min + (pixel_size * (offset + request_size))
 
       [grid_min, grid_max]
     end
 
-    # @param block_x_size [Fixnum]
-    # @param block_y_size [Fixnum]
-    # @param raster_width [Fixnum]
-    # @param raster_height [Fixnum]
-    # @return [Fixnum] The total number of blocks that should be iterated
+    # @param block_x_size [Integer]
+    # @param block_y_size [Integer]
+    # @param raster_width [Integer]
+    # @param raster_height [Integer]
+    # @return [Integer] The total number of blocks that should be iterated
     #   through during the grid+rasterize process.
     def build_block_count(block_x_size, block_y_size, raster_width, raster_height)
       build_block_size(raster_width, block_x_size) * build_block_size(raster_height, block_y_size)
     end
 
-    # @param total_pixels [Fixnum] Number of pixels in the width or height.
-    # @param block_size [Fixnum] Size of the reported block.
+    # @param total_pixels [Integer] Number of pixels in the width or height.
+    # @param block_size [Integer] Size of the reported block.
     def build_block_size(total_pixels, block_size)
       (total_pixels + block_size - 1) / block_size
     end

data/lib/gdal/{gridder_options.rb → extensions/gridder_options.rb}

@@ -1,8 +1,10 @@
 # frozen_string_literal: true
 
-require_relative 'grid'
-require_relative 'options'
-require_relative '../ogr'
+require 'gdal/grid'
+require 'gdal/options'
+require 'ogr/exceptions'
+require 'ogr/geometry'
+require 'ogr/spatial_reference'
 
 module GDAL
   # Object to be used with a {GDAL::Gridder}.
@@ -140,7 +142,7 @@ module GDAL
     def input_clipping_geometry=(geometry)
       unless geometry.is_a?(OGR::Geometry)
         raise OGR::InvalidGeometry,
-          "Clipping geometry must be a OGR::Geometry type, but was a #{geometry.class}"
+              "Clipping geometry must be a OGR::Geometry type, but was a #{geometry.class}"
       end
 
       @input_clipping_geometry = geometry
@@ -151,14 +153,13 @@ module GDAL
       data_types = FFI::GDAL::GDAL::DataType.symbols
 
       unless data_types.include?(type)
-        raise GDAL::InvalidDataType,
-          "output_data_type must be one of #{data_types} but was #{type}"
+        raise GDAL::InvalidDataType, "output_data_type must be one of #{data_types} but was #{type}"
       end
 
       @grid.data_type = @output_data_type = type
     end
 
-    # @return [Fixnum]
+    # @return [Integer]
     def output_data_type_size
       GDAL::DataType.size(@output_data_type) / 8
     end
@@ -168,8 +169,7 @@ module GDAL
       driver_names = GDAL::Driver.short_names
 
       unless driver_names.include?(format)
-        raise GDAL::InvalidDriverName,
-          "output_form must be one of #{driver_names} but was #{format}"
+        raise GDAL::InvalidDriverName, "output_form must be one of #{driver_names} but was #{format}"
       end
 
       @output_format = format
@@ -183,14 +183,14 @@ module GDAL
       @output_driver ||= GDAL::Driver.by_name(@output_format)
     end
 
-    # @param min_max [Array<Fixnum>, Hash{min => Number, max => Number}]
+    # @param min_max [Array<Integer>, Hash{min => Number, max => Number}]
     def output_x_extent=(min_max)
       min, max = extract_min_max(min_max, :min, :max)
 
       @output_x_extent = { min: min, max: max }
     end
 
-    # @param min_max [Array<Fixnum>, Hash{min => Number, max => Number}]
+    # @param min_max [Array<Integer>, Hash{min => Number, max => Number}]
     def output_y_extent=(min_max)
       min, max = extract_min_max(min_max, :min, :max)
 
@@ -213,7 +213,7 @@ module GDAL
     def output_projection=(spatial_reference)
       unless spatial_reference.is_a?(OGR::SpatialReference)
         raise OGR::InvalidSpatialReference,
-          "output_projection must be an OGR::SpatialReference but was a #{spatial_reference.class}"
+              "output_projection must be an OGR::SpatialReference but was a #{spatial_reference.class}"
       end
 
       @output_projection = spatial_reference
@@ -250,8 +250,7 @@ module GDAL
     # @return [Array<Number>]
     def extract_min_max_from_array(content, min_name, max_name)
       unless content.length == 2
-        raise ArgumentError,
-          "Please supply only 2 elements, one for #{min_name}, one for #{max_name}"
+        raise ArgumentError, "Please supply only 2 elements, one for #{min_name}, one for #{max_name}"
       end
 
       [content[0], content[1]]

data/lib/gdal/{raster_attribute_table_mixins → extensions/raster_attribute_table}/extensions.rb

@@ -1,13 +1,13 @@
 # frozen_string_literal: true
 
-require 'json'
+require 'gdal/raster_attribute_table'
 
 module GDAL
-  module RasterAttributeTableMixins
+  class RasterAttributeTable
     module Extensions
-      # @param row [Fixnum]
-      # @param field [Fixnum]
-      # @param value [String, Float, Fixnum]
+      # @param row [Integer]
+      # @param field [Integer]
+      # @param value [String, Float, Integer]
       def set_value(row, field, value)
         case value.class
         when String then set_value_as_string(row, field, value)
@@ -19,7 +19,7 @@ module GDAL
 
       # Get +column_name+, +column_usage+, +column_type+ as a Hash.
       #
-      # @param index [Fixnum]
+      # @param index [Integer]
       # @return [Hash]
       def column(index)
         {
@@ -44,3 +44,5 @@ module GDAL
     end
   end
 end
+
+GDAL::RasterAttributeTable.include(GDAL::RasterAttributeTable::Extensions)

data/lib/gdal/{raster_band_mixins → extensions/raster_band}/algorithm_extensions.rb

@@ -1,9 +1,9 @@
 # frozen_string_literal: true
 
-require_relative '../../ogr'
+require 'gdal/raster_band'
 
 module GDAL
-  module RasterBandMixins
+  class RasterBand
     # RasterBand methods for doing warp-like things, but not using GDAL's Warp
     # API.
     module AlgorithmExtensions
@@ -51,14 +51,14 @@ module GDAL
       #   raster_band = dataset.raster_band(1)
       #
       #   raster_band.simple_erase! do |x, y, value|
-      #     if x < (raster_band.x_size / 2) && y < (raster_band.y_size / 2) && value < 0
+      #     if x < (raster_band.x_size / 2) && y < (raster_band.y_size / 2) && value.negative?
       #       false
       #     else
       #       true
       #     end
       #   end
       #
-      # @return [Fixnum] The number of pixels that were erased.
+      # @return [Integer] The number of pixels that were erased.
       def simple_erase!
         erase_value = no_data_value[:value]
         write_buffer = GDAL._buffer_from_data_type(data_type, x_size)
@@ -66,7 +66,7 @@ module GDAL
 
         unless erase_value
           raise GDAL::NoRasterEraseValue,
-            'Cannot erase values, RasterBand does not have a NODATA value set'
+                'Cannot erase values, RasterBand does not have a NODATA value set'
         end
 
         y_size.times do |line_num|
@@ -98,7 +98,7 @@ module GDAL
 
       # @param write_buffer [FFI::Buffer]
       # @param pixel_values [Array<Number>]
-      # @param line_number [Fixnum]
+      # @param line_number [Integer]
       def rewrite_pixel_row(write_buffer, pixel_values, line_number)
         GDAL._write_pointer(write_buffer, data_type, pixel_values)
         raster_io('w', write_buffer, x_size: x_size, y_size: 1, y_offset: line_number)
@@ -107,3 +107,5 @@ module GDAL
     end
   end
 end
+
+GDAL::RasterBand.include(GDAL::RasterBand::AlgorithmExtensions)

data/lib/gdal/{raster_band_mixins → extensions/raster_band}/coloring_extensions.rb

@@ -1,7 +1,7 @@
 # frozen_string_literal: true
 
 module GDAL
-  module RasterBandMixins
+  class RasterBand
     # RasterBand methods added for dealing with colorizing.
     module ColoringExtensions
       # Sets the band to be a Palette band, then applies an RGB color table using
@@ -9,7 +9,7 @@ module GDAL
       # on the number of colors given.  Note that this will overwrite any existing
       # color table that may be set on this band.
       #
-      # @param colors [Array<Fixnum, String>, Fixnum, String] Can be a single or
+      # @param colors [Array<Integer, String>, Integer, String] Can be a single or
       #   many colors, given as either R, G, and B integers (0-255) or as strings
       #   of Hex.
       #
@@ -45,7 +45,7 @@ module GDAL
           # TODO: Fix possible uninitialized rgb_array
           rgb_array = hex_to_rgb(color) unless color.is_a?(Array)
           table.add_color_entry(color_entry_index,
-            rgb_array[0], rgb_array[1], rgb_array[2], 255)
+                                rgb_array[0], rgb_array[1], rgb_array[2], 255)
         end
 
         self.color_table = table
@@ -54,7 +54,7 @@ module GDAL
       # Gets the colors from the associated ColorTable and returns an Array of
       # those, where each ColorEntry is [R, G, B, A].
       #
-      # @return [Array<Array<Fixnum>>]
+      # @return [Array<Array<Integer>>]
       def colors_as_rgb
         return [] unless color_table
 
@@ -84,3 +84,5 @@ module GDAL
     end
   end
 end
+
+GDAL::RasterBand.include(GDAL::RasterBand::ColoringExtensions)

data/lib/gdal/{raster_band_mixins → extensions/raster_band}/extensions.rb

@@ -1,9 +1,11 @@
 # frozen_string_literal: true
 
-require 'json'
+require 'gdal/raster_band'
+require 'numo/narray'
+require_relative 'io_extensions'
 
 module GDAL
-  module RasterBandMixins
+  class RasterBand
     module Extensions
       # @return [Enumerator]
       def each_overview
@@ -16,7 +18,7 @@ module GDAL
 
       # @return [Array]
       def overviews
-        eaach_overview.to_a
+        each_overview.to_a
       end
 
       # @return [Array]
@@ -44,7 +46,7 @@ module GDAL
       # Each pixel of the raster projected using the dataset's geo_transform.
       # The output NArray is a 3D array where the inner-most array is a the
       # lat an lon, those are contained in an array per pixel line, and finally
-      # the outter array contains each of the pixel lines.
+      # the outer array contains each of the pixel lines.
       #
       # @return [NArray]
       def projected_points
@@ -64,10 +66,12 @@ module GDAL
 
       # The total number of pixels in the raster band.
       #
-      # @return [Fixnum]
+      # @return [Integer]
       def pixel_count
         x_size * y_size
       end
     end
   end
 end
+
+GDAL::RasterBand.include(GDAL::RasterBand::Extensions)

data/lib/gdal/{raster_band_mixins → extensions/raster_band}/io_extensions.rb

@@ -1,7 +1,9 @@
 # frozen_string_literal: true
 
+require 'gdal/raster_band'
+
 module GDAL
-  module RasterBandMixins
+  class RasterBand
     # Methods for reading & writing RasterBands that didn't come from GDAL.
     module IOExtensions
       # Reads the raster line-by-line and returns as an NArray.  Will yield each
@@ -48,8 +50,8 @@ module GDAL
 
       # Convenience method for directly setting a single pixel value.
       #
-      # @param x [Fixnum] Pixel number in the row to set.
-      # @param y [Fixnum] Row number of the pixel to set.
+      # @param x [Integer] Pixel number in the row to set.
+      # @param y [Integer] Row number of the pixel to set.
       # @param new_value [Number] The value to set the pixel to.
       def set_pixel_value(x, y, new_value)
         data_pointer = GDAL._pointer_from_data_type(data_type)
@@ -60,8 +62,8 @@ module GDAL
 
       # Convenience method for directly getting a single pixel value.
       #
-      # @param x [Fixnum] Pixel number in the row to get.
-      # @param y [Fixnum] Row number of the pixel to get.
+      # @param x [Integer] Pixel number in the row to get.
+      # @param y [Integer] Row number of the pixel to get.
       # @return [Number]
       def pixel_value(x, y)
         output = raster_io('r', x_size: 1, y_size: 1, x_offset: x, y_offset: y, buffer_x_size: 1, buffer_y_size: 1)
@@ -74,8 +76,8 @@ module GDAL
       # x and y counts for when the total number of pixels and lines does not
       # divide evenly using GDAL's block count.
       #
-      # @return [Hash{x => Fixnum, x_remainder => Fixnum, y => Fixnum,
-      #   y_remainder => Fixnum}]
+      # @return [Hash{x => Integer, x_remainder => Integer, y => Integer,
+      #   y_remainder => Integer}]
       # @see http://www.gdal.org/classGDALRasterBand.html#a09e1d83971ddff0b43deffd54ef25eef
       def block_count
         x_size_plus_block_size = x_size + block_size[:x]
@@ -91,7 +93,7 @@ module GDAL
 
       # The buffer size to use for block-based IO, based on #block_size.
       #
-      # @return [Fixnum]
+      # @return [Integer]
       def block_buffer_size
         block_size[:x] * block_size[:y]
       end
@@ -120,10 +122,10 @@ module GDAL
       #
       # @yieldparam pixels [Array<Number>] An Array the same size as
       #   {#block_buffer_size} containing all pixel values in the current block.
-      # @yieldparam x_block_size [Fixnum] Instead of using only #{RasterBand#block_size},
+      # @yieldparam x_block_size [Integer] Instead of using only #{RasterBand#block_size},
       #   it will tell you the size of each block--handy for when the last block
       #   is smaller than the rest.
-      # @yieldparam y_block_size [Fixnum] Same as +x_block_siz+ but for y.
+      # @yieldparam y_block_size [Integer] Same as +x_block_siz+ but for y.
       # @return [Enumerator, nil]
       def read_blocks_by_block
         return enum_for(:read_blocks_by_block) unless block_given?
@@ -150,9 +152,9 @@ module GDAL
       # GDAL reports that the block size to read is 60, we still need to know
       # to read those last 5 lines when using block-related methods.
       #
-      # @param block_number [Fixnum] The number of the y block when iterating
+      # @param block_number [Integer] The number of the y block when iterating
       #   through y blocks.
-      # @return [Fixnum] The number of lines to read as part of the block.
+      # @return [Integer] The number of lines to read as part of the block.
       def calculate_y_block_size(block_number)
         # If it's the last block...
         if block_number == (block_count[:y] - 1)
@@ -167,9 +169,9 @@ module GDAL
       # and GDAL reports that the block size to read is 60, we still need to
       # know to read those last 5 pixels when using block-related methods.
       #
-      # @param block_number [Fixnum] The number of the x block when iterating
+      # @param block_number [Integer] The number of the x block when iterating
       #   through x blocks.
-      # @return [Fixnum] The number of pixels to read as part of the block.
+      # @return [Integer] The number of pixels to read as part of the block.
       def calculate_x_block_size(block_number)
         # If it's the last block...
         if block_number == (block_count[:x] - 1)
@@ -181,3 +183,5 @@ module GDAL
     end
   end
 end
+
+GDAL::RasterBand.include(GDAL::RasterBand::IOExtensions)

data/lib/gdal/{raster_band_classifier.rb → extensions/raster_band_classifier.rb}

@@ -1,6 +1,7 @@
 # frozen_string_literal: true
 
-require_relative '../gdal'
+require 'gdal'
+require 'numo/narray'
 
 module GDAL
   # Takes a list of Ranges of color values and remaps them.  Note that these
@@ -50,7 +51,7 @@ module GDAL
     # calculate evenly-weighted ranges.  If there are remainder values at the
     # max end of the values, those get lumped in with the last range.
     #
-    # @param range_count [Fixnum] The number of ranges to create.
+    # @param range_count [Integer] The number of ranges to create.
     # @return [Array<Hash>, nil]
     def equal_count_ranges(range_count)
       pixels = @raster_band.to_nna
@@ -103,9 +104,6 @@ module GDAL
         new_band_pixels[data_pixels & band_pixels.le(r[:range].max) & band_pixels.ge(r[:range].min)] = r[:map_to]
       end
 
-      in_range = @ranges.map { |r| band_pixels.le(r[:range].max) & band_pixels.ge(r[:range].min) }.reduce(:|)
-      new_band_pixels[data_pixels & ~in_range] = 0 if in_range
-
       mask_nan(new_band_pixels, data_pixels) if nodata_is_nan?
       @raster_band.write_xy_narray(new_band_pixels)
     end

data/lib/gdal/geo_transform.rb

@@ -1,13 +1,10 @@
 # frozen_string_literal: true
 
 require_relative '../gdal'
-require_relative 'geo_transform_mixins/extensions'
 
 module GDAL
   class GeoTransform
-    include GeoTransformMixins::Extensions
-
-    # @return [FFI::Pointer] A pointer that can be used to hold a GeoTransform.
+    # @return [FFI::MemoryPointer] A pointer that can be used to hold a GeoTransform.
     def self.new_pointer
       FFI::MemoryPointer.new(:double, 6)
     end
@@ -41,8 +38,6 @@ module GDAL
       self.pixel_height ||= 1.0
       self.x_rotation ||= 0.0
       self.y_rotation ||= 0.0
-
-      to_a
     end
 
     def null?

data/lib/gdal/grid.rb

@@ -25,10 +25,10 @@ module GDAL
     end
 
     # @param points [Array,NArray] An Array containing all x, y, and z points.
-    # @param extents [Hash{x_min: Fixnum, y_min: Fixnum, x_max: Fixnum, y_max: Fixnum}]
+    # @param extents [Hash{x_min: Integer, y_min: Integer, x_max: Integer, y_max: Integer}]
     # @param data_pointer [FFI::Pointer] Pointer that will contain the gridded
     #   data (after this method is done).
-    # @param output_size [Hash{x: Fixnum, y: Fixnum}] Overall dimensions of the
+    # @param output_size [Hash{x: Integer, y: Integer}] Overall dimensions of the
     #   area of the output raster to grid.
     # @param progress_block [Proc]
     # @param progress_arg [FFI::Pointer]
@@ -74,6 +74,7 @@ module GDAL
     # @param points [Array]
     def make_points_pointer(points)
       raise GDAL::Error, 'No points to make pointer for' if points.compact.empty?
+
       input_coordinates_ptr = FFI::MemoryPointer.new(:double, points.length)
       input_coordinates_ptr.write_array_of_double(points)
 

data/lib/gdal/grid_algorithms.rb

@@ -3,22 +3,22 @@
 module GDAL
   module GridAlgorithms
     autoload :InverseDistanceToAPower,
-      File.expand_path('grid_algorithms/inverse_distance_to_a_power', __dir__)
+             File.expand_path('grid_algorithms/inverse_distance_to_a_power', __dir__)
     autoload :MetricAverageDistance,
-      File.expand_path('grid_algorithms/metric_average_distance', __dir__)
+             File.expand_path('grid_algorithms/metric_average_distance', __dir__)
     autoload :MetricAverageDistancePts,
-      File.expand_path('grid_algorithms/metric_average_distance_pts', __dir__)
+             File.expand_path('grid_algorithms/metric_average_distance_pts', __dir__)
     autoload :MetricCount,
-      File.expand_path('grid_algorithms/metric_count', __dir__)
+             File.expand_path('grid_algorithms/metric_count', __dir__)
     autoload :MetricMaximum,
-      File.expand_path('grid_algorithms/metric_maximum', __dir__)
+             File.expand_path('grid_algorithms/metric_maximum', __dir__)
     autoload :MetricMinimum,
-      File.expand_path('grid_algorithms/metric_minimum', __dir__)
+             File.expand_path('grid_algorithms/metric_minimum', __dir__)
     autoload :MetricRange,
-      File.expand_path('grid_algorithms/metric_range', __dir__)
+             File.expand_path('grid_algorithms/metric_range', __dir__)
     autoload :MovingAverage,
-      File.expand_path('grid_algorithms/moving_average', __dir__)
+             File.expand_path('grid_algorithms/moving_average', __dir__)
     autoload :NearestNeighbor,
-      File.expand_path('grid_algorithms/nearest_neighbor', __dir__)
+             File.expand_path('grid_algorithms/nearest_neighbor', __dir__)
   end
 end