nswtopo 2.0.0.pre.beta1

data/lib/nswtopo/helpers/colour.rb

@@ -0,0 +1,176 @@
+class Colour
+  Error = Class.new RuntimeError
+
+  COLOURS = YAML.load <<~YAML
+    aliceblue: [240, 248, 255]
+    antiquewhite: [250, 235, 215]
+    aqua: [0, 255, 255]
+    aquamarine: [127, 255, 212]
+    azure: [240, 255, 255]
+    beige: [245, 245, 220]
+    bisque: [255, 228, 196]
+    black: [0, 0, 0]
+    blanchedalmond: [255, 235, 205]
+    blue: [0, 0, 255]
+    blueviolet: [138, 43, 226]
+    brown: [165, 42, 42]
+    burlywood: [222, 184, 135]
+    cadetblue: [95, 158, 160]
+    chartreuse: [127, 255, 0]
+    chocolate: [210, 105, 30]
+    coral: [255, 127, 80]
+    cornflowerblue: [100, 149, 237]
+    cornsilk: [255, 248, 220]
+    crimson: [220, 20, 60]
+    cyan: [0, 255, 255]
+    darkblue: [0, 0, 139]
+    darkcyan: [0, 139, 139]
+    darkgoldenrod: [184, 134, 11]
+    darkgray: [169, 169, 169]
+    darkgreen: [0, 100, 0]
+    darkgrey: [169, 169, 169]
+    darkkhaki: [189, 183, 107]
+    darkmagenta: [139, 0, 139]
+    darkolivegreen: [85, 107, 47]
+    darkorange: [255, 140, 0]
+    darkorchid: [153, 50, 204]
+    darkred: [139, 0, 0]
+    darksalmon: [233, 150, 122]
+    darkseagreen: [143, 188, 143]
+    darkslateblue: [72, 61, 139]
+    darkslategray: [47, 79, 79]
+    darkslategrey: [47, 79, 79]
+    darkturquoise: [0, 206, 209]
+    darkviolet: [148, 0, 211]
+    deeppink: [255, 20, 147]
+    deepskyblue: [0, 191, 255]
+    dimgray: [105, 105, 105]
+    dimgrey: [105, 105, 105]
+    dodgerblue: [30, 144, 255]
+    firebrick: [178, 34, 34]
+    floralwhite: [255, 250, 240]
+    forestgreen: [34, 139, 34]
+    fuchsia: [255, 0, 255]
+    gainsboro: [220, 220, 220]
+    ghostwhite: [248, 248, 255]
+    gold: [255, 215, 0]
+    goldenrod: [218, 165, 32]
+    gray: [128, 128, 128]
+    grey: [128, 128, 128]
+    green: [0, 128, 0]
+    greenyellow: [173, 255, 47]
+    honeydew: [240, 255, 240]
+    hotpink: [255, 105, 180]
+    indianred: [205, 92, 92]
+    indigo: [75, 0, 130]
+    ivory: [255, 255, 240]
+    khaki: [240, 230, 140]
+    lavender: [230, 230, 250]
+    lavenderblush: [255, 240, 245]
+    lawngreen: [124, 252, 0]
+    lemonchiffon: [255, 250, 205]
+    lightblue: [173, 216, 230]
+    lightcoral: [240, 128, 128]
+    lightcyan: [224, 255, 255]
+    lightgoldenrodyellow: [250, 250, 210]
+    lightgray: [211, 211, 211]
+    lightgreen: [144, 238, 144]
+    lightgrey: [211, 211, 211]
+    lightpink: [255, 182, 193]
+    lightsalmon: [255, 160, 122]
+    lightseagreen: [32, 178, 170]
+    lightskyblue: [135, 206, 250]
+    lightslategray: [119, 136, 153]
+    lightslategrey: [119, 136, 153]
+    lightsteelblue: [176, 196, 222]
+    lightyellow: [255, 255, 224]
+    lime: [0, 255, 0]
+    limegreen: [50, 205, 50]
+    linen: [250, 240, 230]
+    magenta: [255, 0, 255]
+    maroon: [128, 0, 0]
+    mediumaquamarine: [102, 205, 170]
+    mediumblue: [0, 0, 205]
+    mediumorchid: [186, 85, 211]
+    mediumpurple: [147, 112, 219]
+    mediumseagreen: [60, 179, 113]
+    mediumslateblue: [123, 104, 238]
+    mediumspringgreen: [0, 250, 154]
+    mediumturquoise: [72, 209, 204]
+    mediumvioletred: [199, 21, 133]
+    midnightblue: [25, 25, 112]
+    mintcream: [245, 255, 250]
+    mistyrose: [255, 228, 225]
+    moccasin: [255, 228, 181]
+    navajowhite: [255, 222, 173]
+    navy: [0, 0, 128]
+    oldlace: [253, 245, 230]
+    olive: [128, 128, 0]
+    olivedrab: [107, 142, 35]
+    orange: [255, 165, 0]
+    orangered: [255, 69, 0]
+    orchid: [218, 112, 214]
+    palegoldenrod: [238, 232, 170]
+    palegreen: [152, 251, 152]
+    paleturquoise: [175, 238, 238]
+    palevioletred: [219, 112, 147]
+    papayawhip: [255, 239, 213]
+    peachpuff: [255, 218, 185]
+    peru: [205, 133, 63]
+    pink: [255, 192, 203]
+    plum: [221, 160, 221]
+    powderblue: [176, 224, 230]
+    purple: [128, 0, 128]
+    red: [255, 0, 0]
+    rosybrown: [188, 143, 143]
+    royalblue: [65, 105, 225]
+    saddlebrown: [139, 69, 19]
+    salmon: [250, 128, 114]
+    sandybrown: [244, 164, 96]
+    seagreen: [46, 139, 87]
+    seashell: [255, 245, 238]
+    sienna: [160, 82, 45]
+    silver: [192, 192, 192]
+    skyblue: [135, 206, 235]
+    slateblue: [106, 90, 205]
+    slategray: [112, 128, 144]
+    slategrey: [112, 128, 144]
+    snow: [255, 250, 250]
+    springgreen: [0, 255, 127]
+    steelblue: [70, 130, 180]
+    tan: [210, 180, 140]
+    teal: [0, 128, 128]
+    thistle: [216, 191, 216]
+    tomato: [255, 99, 71]
+    turquoise: [64, 224, 208]
+    violet: [238, 130, 238]
+    wheat: [245, 222, 179]
+    white: [255, 255, 255]
+    whitesmoke: [245, 245, 245]
+    yellow: [255, 255, 0]
+    yellowgreen: [154, 205, 50]
+  YAML
+
+  def initialize(string_or_array)
+    @triplet = case string_or_array
+    when Array then string_or_array.take(3).map(&:round)
+    when *COLOURS.keys
+      @name = string_or_array
+      COLOURS[string_or_array]
+    when /^#?([0-9A-F]{2})([0-9A-F]{2})([0-9A-F]{2})$/i
+      [$1, $2, $3].map { |hex| Integer("0x#{hex}") }
+    when /^rgb\((\d{1,3}),(\d{1,3}),(\d{1,3})\)$/
+      [$1, $2, $3].map(&:to_i)
+    end
+    raise Error, "invalid colour: #{string_or_array}" unless @triplet&.all?(0..255)
+  end
+  attr_reader :triplet
+
+  def mix(other, fraction)
+    Colour.new [triplet, other.triplet].along(fraction.to_f).map(&:to_i)
+  end
+
+  def to_s
+    @name || "#%.2X%.2X%.2X" % triplet
+  end
+end

data/lib/nswtopo/helpers/concurrently.rb

@@ -0,0 +1,27 @@
+module Concurrently
+  CORES = Etc.nprocessors rescue 1
+
+  def concurrently(threads = CORES, &block)
+    elements = Queue.new
+    threads.times.map do
+      Thread.new do
+        while element = elements.pop
+          block.call element
+        end
+      end
+    end.tap do
+      inject(elements, &:<<).close
+    end.each(&:join)
+    self
+  end
+
+  def concurrent_groups(threads = CORES, &block)
+    group_by.with_index do |item, index|
+      index % threads
+    end.values.map do |items|
+      Thread.new(items, &block)
+    end.each(&:join)
+  end
+end
+
+Enumerator.send :include, Concurrently

data/lib/nswtopo/helpers/dir.rb

@@ -0,0 +1,7 @@
+class Dir
+  def self.mktmppath
+    mktmpdir do |path|
+      yield Pathname.new(path)
+    end
+  end
+end

data/lib/nswtopo/helpers/hash.rb

@@ -0,0 +1,15 @@
+module HashHelpers
+  def deep_merge(other)
+    merge(other) do |key, old_value, new_value|
+      Hash === old_value ? Hash === new_value ? old_value.deep_merge(new_value) : new_value : new_value
+    end
+  end
+
+  def deep_merge!(other)
+    merge!(other) do |key, old_value, new_value|
+      Hash === old_value ? Hash === new_value ? old_value.deep_merge!(new_value) : new_value : new_value
+    end
+  end
+end
+
+Hash.send :include, HashHelpers

data/lib/nswtopo/helpers/tar_writer.rb

@@ -0,0 +1,11 @@
+module TarWriterHelpers
+  def add_entry(entry)
+    check_closed
+    @io.write entry.header
+    @io.write entry.read
+    @io.write ?\0 while @io.pos % 512 > 0
+    self
+  end
+end
+
+Gem::Package::TarWriter.send :include, TarWriterHelpers

data/lib/nswtopo/helpers.rb

@@ -0,0 +1,6 @@
+require_relative 'helpers/hash'
+require_relative 'helpers/array'
+require_relative 'helpers/dir'
+require_relative 'helpers/concurrently'
+require_relative 'helpers/tar_writer'
+require_relative 'helpers/colour'

data/lib/nswtopo/layer/arcgis_raster.rb

@@ -0,0 +1,73 @@
+module NSWTopo
+  module ArcGISRaster
+    include Raster, Log
+    CREATE = %w[url]
+
+    def get_raster(temp_dir)
+      raise "no resolution specified for #{@name}" unless Numeric === @resolution
+      txt_path = temp_dir / "mosaic.txt"
+      vrt_path = temp_dir / "mosaic.vrt"
+
+      ArcGISServer.start @url do |connection, service, projection|
+        local_bbox = @map.bounding_box
+        target_bbox = local_bbox.reproject_to projection
+        target_resolution = @resolution * Math::sqrt(target_bbox.first.area / local_bbox.first.area)
+
+        raise "not a tiled map or image server: #{@url}" unless tile_info = service["tileInfo"]
+        lods = tile_info["lods"]
+        origin = tile_info["origin"].values_at "x", "y"
+        tile_sizes = tile_info.values_at "cols", "rows"
+
+        lods.sort_by! do |lod|
+          -lod["resolution"]
+        end
+        lod = lods.find do |lod|
+          lod["resolution"] < target_resolution
+        end || lods.last
+        tile_level, tile_resolution = lod.values_at "level", "resolution"
+
+        tiles = target_bbox.coordinates.first.map do |corner|
+          corner.minus(origin)
+        end.transpose.map(&:minmax).zip(tile_sizes).map do |bound, tile_size|
+          bound / tile_resolution / tile_size
+        end.map do |min, max|
+          (min.floor..max.ceil).each_cons(2).to_a
+        end.inject(&:product).map do |cols, rows|
+          bounds = [cols, rows].zip(tile_sizes).map do |indices, tile_size|
+            indices.times(tile_size * tile_resolution)
+          end.transpose.map do |corner|
+            corner.plus(origin)
+          end.transpose
+
+          bbox = bounds.inject(&:product).values_at(0,2,3,1)
+          next unless target_bbox.first.clip(bbox)
+
+          row, col = rows[1].abs, cols[0]
+          rel_path = "tile/#{tile_level}/#{row}/#{col}"
+          jpg_path = temp_dir / "#{row}.#{col}" # could be png
+          tif_path = temp_dir / "#{row}.#{col}.tif"
+
+          ullr = bounds.inject(&:product).values_at(1,2).flatten
+          gdal_args = ["-a_srs", projection, "-a_ullr", *ullr, "-of", "GTiff", jpg_path, tif_path]
+
+          [rel_path, jpg_path, gdal_args, tif_path]
+        end.compact
+        tiles.each.with_index do |(rel_path, jpg_path, gdal_args, tif_path), index|
+          log_update "%s: retrieving tile %i of %i" % [@name, index + 1, tiles.length]
+          connection.get(rel_path, blankTile: true) do |response|
+            jpg_path.binwrite response.body
+          end
+        end
+      end.each do |rel_path, jpg_path, gdal_args, tif_path|
+        OS.gdal_translate *gdal_args
+      end.map(&:last).tap do |tif_paths|
+        txt_path.write tif_paths.join(?\n)
+        OS.gdalbuildvrt "-input_file_list", txt_path, vrt_path
+      end
+
+      OS.gdal_translate vrt_path, Pathname.pwd / "foo.tif"
+
+      return @resolution, vrt_path
+    end
+  end
+end

data/lib/nswtopo/layer/contour.rb

@@ -0,0 +1,233 @@
+module NSWTopo
+  module Contour
+    include Vector, DEM, Log
+    CREATE = %w[interval index smooth simplify thin density min-length no-depression knolls fill]
+    DEFAULTS = YAML.load <<~YAML
+      interval: 5
+      smooth: 0.2
+      density: 4.0
+      min-length: 2.0
+      knolls: 0.2
+      section: 100
+      stroke: "#805100"
+      stroke-width: 0.08
+      Depression:
+        symbolise:
+          interval: 2.0
+          line:
+            stroke-width: 0.12
+            y2: -0.3
+      labels:
+        font-size: 1.4
+        letter-spacing: 0.05
+        orientation: downhill
+        collate: true
+        min-radius: 5
+        max-turn: 20
+        sample: 10
+        minimum-area: 70
+        separation: 40
+        separation-all: 15
+        separation-along: 100
+    YAML
+
+    def margin
+      { mm: [3 * @smooth, 1].min }
+    end
+
+    def check_geos!
+      json = OS.ogr2ogr "-dialect", "SQLite", "-sql", "SELECT geos_version() AS version", "-f", "GeoJSON", "/vsistdout/", "/vsistdin/" do |stdin|
+        stdin.write GeoJSON::Collection.new.to_json
+      end
+      raise unless version = JSON.parse(json).dig("features", 0, "properties", "version")
+      raise unless (version.split(?-).first.split(?.).map(&:to_i) <=> [3, 3]) >= 0
+    rescue OS::Error, JSON::ParserError, RuntimeError
+      raise "contour thinning requires GDAL with SpatiaLite and GEOS support"
+    end
+
+    def get_features
+      @simplify ||= [0.5 * @interval / Math::tan(Math::PI * 85 / 180), 0.001 * 0.05 * @map.scale].min
+      @index ||= 10 * @interval
+      @params = {
+        "Index" => { "stroke-width" => 2 * @params["stroke-width"] },
+        "labels" => { "fill" => @fill || @params["stroke"] }
+      }.deep_merge(@params)
+
+      check_geos! if @thin
+      raise "%im index interval not a multiple of %im contour interval" % [@index, @interval] unless @index % @interval == 0
+
+      Dir.mktmppath do |temp_dir|
+        dem_path, blur_path = temp_dir / "dem.tif", temp_dir / "dem.blurred.tif"
+
+        if @smooth.zero?
+          get_dem temp_dir, blur_path
+        else
+          get_dem temp_dir, dem_path
+          blur_dem dem_path, blur_path
+        end
+
+        db_flags = @thin ? %w[-f SQLite -dsco SPATIALITE=YES] : ["-f", "ESRI Shapefile"]
+        db_path = temp_dir / "contour"
+
+        log_update "%s: generating contour lines" % @name
+        json = OS.gdal_contour "-q", "-a", "elevation", "-i", @interval, "-f", "GeoJSON", "-lco", "RFC7946=NO", blur_path, "/vsistdout/"
+        contours = GeoJSON::Collection.load json, @map.projection
+
+        if @no_depression.nil?
+          candidates = contours.select do |feature|
+            feature.coordinates.last == feature.coordinates.first &&
+            feature.coordinates.anticlockwise?
+          end.each do |feature|
+            feature["depression"] = 1
+          end
+          index = RTree.load(candidates, &:bounds)
+
+          contours.reject! do |feature|
+            next unless feature["depression"] == 1
+            index.search(feature.bounds).none? do |other|
+              next if other == feature
+              feature.coordinates.first.within?(other.coordinates) ||
+              other.coordinates.first.within?(feature.coordinates)
+            end
+          end
+        end
+
+        contours.reject! do |feature|
+          feature.coordinates.last == feature.coordinates.first &&
+          feature.bounds.all? { |min, max| max - min < @knolls * @map.scale / 1000.0 }
+        end
+
+        contours.each do |feature|
+          id, elevation, depression = feature.values_at "ID", "elevation", "depression"
+          feature.properties.replace("id" => id, "elevation" => elevation, "modulo" => elevation % @index, "depression" => depression || 0)
+        end
+
+        contours.reject! do |feature|
+          feature["elevation"].zero?
+        end
+
+        OS.ogr2ogr "-a_srs", @map.projection, "-nln", "contour", "-simplify", @simplify, *db_flags, db_path, "GeoJSON:/vsistdin/" do |stdin|
+          stdin.write contours.to_json
+        end
+
+        if @thin
+          slope_tif_path = temp_dir / "slope.tif"
+          slope_vrt_path = temp_dir / "slope.vrt"
+          min_length = @min_length * @map.scale / 1000.0
+
+          log_update "%s: generating slope masks" % @name
+          OS.gdaldem "slope", blur_path, slope_tif_path, "-compute_edges"
+          json = OS.gdalinfo "-json", slope_tif_path
+          width, height = JSON.parse(json)["size"]
+          srcwin = [ -2, -2, width + 4, height + 4 ]
+          OS.gdal_translate "-srcwin", *srcwin, "-a_nodata", "none", "-of", "VRT", slope_tif_path, slope_vrt_path
+
+          multiplier = @index / @interval
+          case multiplier
+          when  4 then [ [1,3], 2 ]
+          when  5 then [ [1,4], [2,3] ]
+          when  6 then [ [1,4], [2,5], 3 ]
+          when  7 then [ [2,5], [1,3,6], 4 ]
+          when  8 then [ [1,3,5,7], [2,6], 4 ]
+          when  9 then [ [1,4,7], [2,5,8], [3,6] ]
+          when 10 then [ [2,5,8], [1,4,6,9], [3,7] ]
+          else raise "contour thinning not available for specified index interval"
+          end.inject(multiplier) do |count, (*drop)|
+            angle = Math::atan(1000.0 * @index * @density / @map.scale / count) * 180.0 / Math::PI
+            mask_path = temp_dir / "mask.#{count}.sqlite"
+
+            OS.gdal_contour "-nln", "ring", "-a", "angle", "-fl", angle, *db_flags, slope_vrt_path, mask_path
+
+            OS.ogr2ogr "-update", "-nln", "mask", "-nlt", "MULTIPOLYGON", mask_path, mask_path, "-dialect", "SQLite", "-sql", <<~SQL
+              SELECT
+                ST_Buffer(ST_Buffer(ST_Polygonize(geometry), #{0.5 * min_length}, 6), #{-0.5 * min_length}, 6) AS geometry
+              FROM ring
+            SQL
+
+            drop.each do |index|
+              OS.ogr2ogr "-nln", "mask", "-update", "-append", "-explodecollections", "-q", db_path, mask_path, "-dialect", "SQLite", "-sql", <<~SQL
+                SELECT geometry, #{index * @interval} AS modulo
+                FROM mask
+              SQL
+            end
+
+            count - drop.count
+          end
+
+          log_update "%s: thinning contour lines" % @name
+          OS.ogr2ogr "-nln", "divided", "-update", "-explodecollections", db_path, db_path, "-dialect", "SQLite", "-sql", <<~SQL
+            WITH intersecting(contour, mask) AS (
+              SELECT contour.rowid, mask.rowid
+              FROM contour
+              INNER JOIN mask
+              ON
+                mask.modulo = contour.modulo AND
+                contour.rowid IN (
+                  SELECT rowid FROM SpatialIndex
+                  WHERE
+                    f_table_name = 'contour' AND
+                    search_frame = mask.geometry
+                ) AND
+                ST_Relate(contour.geometry, mask.geometry, 'T********')
+            )
+
+            SELECT contour.geometry, contour.id, contour.elevation, contour.modulo, contour.depression, 1 AS unmasked, 1 AS unaltered
+            FROM contour
+            LEFT JOIN intersecting ON intersecting.contour = contour.rowid
+            WHERE intersecting.contour IS NULL
+
+            UNION SELECT ExtractMultiLinestring(ST_Difference(contour.geometry, ST_Collect(mask.geometry))) AS geometry, contour.id, contour.elevation, contour.modulo, contour.depression, 1 AS unmasked, 0 AS unaltered
+            FROM contour
+            INNER JOIN intersecting ON intersecting.contour = contour.rowid
+            INNER JOIN mask ON intersecting.mask = mask.rowid
+            GROUP BY contour.rowid
+            HAVING min(ST_Relate(contour.geometry, mask.geometry, '**T******'))
+
+            UNION SELECT ExtractMultiLinestring(ST_Intersection(contour.geometry, ST_Collect(mask.geometry))) AS geometry, contour.id, contour.elevation, contour.modulo, contour.depression, 0 AS unmasked, 0 AS unaltered
+            FROM contour
+            INNER JOIN intersecting ON intersecting.contour = contour.rowid
+            INNER JOIN mask ON intersecting.mask = mask.rowid
+            GROUP BY contour.rowid
+          SQL
+
+          OS.ogr2ogr "-nln", "thinned", "-update", "-explodecollections", db_path, db_path, "-dialect", "SQLite", "-sql", <<~SQL
+            SELECT ST_LineMerge(ST_Collect(geometry)) AS geometry, id, elevation, modulo, depression, unaltered
+            FROM divided
+            WHERE unmasked OR ST_Length(geometry) < #{min_length}
+            GROUP BY id, elevation, modulo, unaltered
+          SQL
+
+          OS.ogr2ogr "-nln", "contour", "-update", "-overwrite", db_path, db_path, "-dialect", "SQLite", "-sql", <<~SQL
+            SELECT geometry, id, elevation, modulo, depression
+            FROM thinned
+            WHERE unaltered OR ST_Length(geometry) > #{min_length}
+          SQL
+        end
+
+        json = OS.ogr2ogr "-f", "GeoJSON", "-lco", "RFC7946=NO", "/vsistdout/", db_path, "contour"
+        GeoJSON::Collection.load(json, @map.projection).each do |feature|
+          elevation, modulo, depression = feature.values_at "elevation", "modulo", "depression"
+          category = modulo.zero? ? %w[Index] : %w[Standard]
+          category << "Depression" if depression == 1
+          feature.clear
+          feature["elevation"] = elevation
+          feature["category"] = category
+          feature["label"] = elevation.to_i.to_s if modulo.zero?
+        end
+      end
+    end
+
+    def to_s
+      elevations = features.map do |feature|
+        [feature["elevation"], feature["category"].include?("Index")]
+      end.uniq.sort_by(&:first)
+      range = elevations.map(&:first).minmax
+      interval, index = %i[itself last].map do |selector|
+        elevations.select(&selector).map(&:first).each_cons(2).map { |e0, e1| e1 - e0 }.min
+      end
+      [["%im intervals", interval], ["%im indices", index], ["%im-%im elevation", (range if range.all?)]].select(&:last).map do |label, value|
+        label % value
+      end.join(", ").prepend("%s: " % @name)
+    end
+  end
+end

data/lib/nswtopo/layer/control.rb

@@ -0,0 +1,94 @@
+module NSWTopo
+  module Control
+    include Vector
+    CREATE = %w[diameter spot font-size colour]
+    DEFAULTS = YAML.load <<~YAML
+      diameter: 7.0
+      colour: darkmagenta
+      stroke-width: 0.25
+      waterdrop:
+        stroke: blue
+      labels:
+        dupe: outline
+        outline:
+          stroke: white
+          fill: none
+          stroke-width: 0.25
+          stroke-opacity: 0.75
+        position: [ aboveright, belowright, aboveleft, belowleft, right, left, above, below ]
+        font-family: sans-serif
+        font-style: normal
+        stroke: none
+    YAML
+    SCALING_PARAMS = <<~YAML
+      fence: 2.0
+      control:
+        symbol:
+        - circle:
+            r: 1.0
+            fill: none
+      hashhouse:
+        symbol:
+          path:
+            d: M 0.0 -1.0 L -0.866 0.5 L 0.866 0.5 Z
+            fill: none
+      anc:
+        symbol:
+          path:
+            d: M 0.7071 0.7071 L -0.7071 0.7071 L -0.7071 -0.7071 L 0.7071 -0.7071 Z
+            fill: none
+      waterdrop:
+        symbol:
+          path:
+            d:
+              M 0 0 m -0.63954,0.063887 -0.0040064,0.32652 0.33453,0.034055 0,-0.38261 -0.33052,0.022034 z
+              M 0 0 m -0.0095612,-0.43108 0,0.11413
+              M 0 0 m 0.11225,-0.43108 0,0.11413
+              M 0 0 m -0.30901,0.04046 c 0.020447,0.0013422 0.029653,0.024004 0.049078,0.031048 0.018777,0.0053462 0.038045,0.0035119 0.05709,0.0020132 0.0063728,-0.0033554 0.042405,-0.030227 0.044414,-0.034299 0.01612,-0.013399 0.037913,-0.04687 0.037717,-0.065861 l 0.0040064,-0.10517 c 0.00020519,-0.0053688 -0.021845,-0.022347 -0.025354,-0.026427 -0.0042688,-0.0049656 0.0055968,-0.05284 -0.0026192,-0.051919 -0.0092762,-0.0035792 -0.013901,0.00044738 -0.02763,-0.0061064 -0.0033824,-0.010245 -0.0050752,-0.053478 0.0081307,-0.051037 0.022266,0.0013424 0.036219,0.0040488 0.061533,-0.0045856 0.015342,0 0.068088,-0.048747 0.083316,-0.048747 l 0.13477,0.00067107 c 0.0086143,0 0.025117,0.022723 0.029266,0.026872 0.019873,0.024221 0.046404,0.029805 0.082646,0.041671 0.0092704,0.0011184 0.041429,-0.0045408 0.044242,0.0020136 0.0063744,0.021228 0.0042104,0.0427 0.0024896,0.050818 -0.0023072,0.010894 -0.023914,-0.00015658 -0.025061,0.013936 -0.0009304,0.011453 -0.00032659,0.031954 -0.0034784,0.046116 -0.019142,0.019707 -0.042333,0.032122 -0.034371,0.040092 l -0.00052567,0.090013 c 0.0024592,0.02049 0.0088222,0.043385 0.016268,0.056554 0.020901,0.046083 0.064072,0.047163 0.06903,0.04844 0.0022744,0.00067108 0.1974,0.025825 0.30749,0.13121 0.0093678,0.00897 0.029965,0.026342 0.053084,0.097154 0.02234,0.068425 0.037538,0.2914 0.030048,0.29447 -0.06374,0.026123 -0.15033,0.053353 -0.2534,0.020043 -0.01286,-0.0041608 0.0013624,-0.088139 -0.0070112,-0.1282 -0.0085866,-0.041085 0.0002729,-0.12091 -0.08864,-0.13672 -0.037049,-0.0065768 -0.15685,0.014316 -0.24088,0.03205 -0.040012,0.0084555 -0.094156,0.013712 -0.12269,0.0080082 -0.017636,-0.0035344 -0.12587,-0.027022 -0.1458,-0.032006 -0.013224,-0.0033104 -0.042701,-0.0050104 -0.087246,0.0088806 -0.033064,0.010312 -0.050813,0.035506 -0.049812,0.026722 0.00073147,-0.00642 0.0012504,-0.371 -8.5944e-05,-0.37768 z
+              M 0 0 m 0.050734,-0.6446 c -0.020293,-0.00022369 -0.052848,0.029116 -0.09202,0.057281 -0.051683,0.032021 -0.11317,0.0094174 -0.16568,-0.00774 -0.049831,-0.019774 -0.12933,-0.0033552 -0.13608,0.058974 0.0013488,0.036319 0.014608,0.092861 0.058748,0.095185 0.074174,-0.010044 0.15057,-0.043112 0.22528,-0.018924 0.029827,0.013175 0.035194,0.033585 0.037678,0.033102 l 0.14636,0 c 0.002484,0.00044738 0.0078512,-0.019931 0.037678,-0.033102 0.074708,-0.024188 0.1511,0.0088806 0.22528,0.018924 0.044141,-0.0023264 0.057457,-0.058867 0.058805,-0.095185 -0.006752,-0.062332 -0.086307,-0.078753 -0.13614,-0.058974 -0.052507,0.017157 -0.114,0.039761 -0.16568,0.00774 -0.039878,-0.028673 -0.072888,-0.058556 -0.093094,-0.057225 -0.00036462,-2.4608e-05 -0.00072.451,-5.1456e-05 -0.0011296,-5.5928e-05 z
+              M 0 0 m -0.16679,-0.26528 c 0.23655,0.036829 0.43804,0.013466 0.43804,0.013466
+              M 0 0 m -0.17245,-0.2153 c 0.23655,0.036826 0.44371,0.017694 0.44371,0.017694
+            fill: none
+      labels:
+        margin: 1.4142
+        font-size: 1.5
+    YAML
+
+    def get_features
+      scaled_params = SCALING_PARAMS.gsub(/\-?\d\.\d+/) { |number| "%.5g" % (number.to_f * 0.5 * @diameter) }
+      scaled_params = YAML.load scaled_params
+      scaled_params["control"]["symbol"] << { "circle" => { "r" => 0.07, "stroke-width" => 0.14, "fill" => "none" } } if @spot
+      @params = scaled_params.deep_merge @params
+      @params["labels"]["font-size"] = @font_size if @font_size
+      @params["labels"]["fill"] = @params["stroke"] = @colour.to_s if @colour
+      points, controls = GPS.load(@path).points, GeoJSON::Collection.new
+      [["control",   /^(1?\d\d)W?$/ ],
+       ["hashhouse", /^(HH)$/       ],
+       ["anc",       /^(ANC)$/      ],
+       ["waterdrop", /^1?\d\dW$|^W$/],
+      ].each do |type, selector|
+        points.each do |point|
+          name = point["name"]
+          next unless name =~ selector
+          properties = [["category", [type, *$1]], ["label", $1]].select(&:last).to_h
+          controls.add_point point.coordinates, properties
+        end
+      end
+      controls
+    end
+
+    def to_s
+      counts = %w[control waterdrop hashhouse].map do |category|
+        waypoints = features.select do |feature|
+          feature["category"].any? category
+        end
+        next if waypoints.empty?
+        count = "%i %s%s" % [waypoints.length, category, waypoints.one? ? nil : ?s]
+        next count unless "control" == category
+        total = features.sum { |feature| feature["label"].to_i.floor(-1) }
+        count << " (%i points)" % total
+      end.compact
+      [@name, counts.join(", ")].join(": ")
+    end
+  end
+end