panomosity 0.1.36 → 0.1.41

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: ed12e60dac2b5b41da899dadf401e2c51eaa498508fd78a6d0280d3bb262208f
-  data.tar.gz: 511594cc4b28962db894a9ad399a5e8b81f90808ff1a1ba32de4fab6938a4760
+  metadata.gz: 43cc742ee7f3a7e491f8e2ffe3a2f47096c649949e4c9a7cc144c7a2baea819c
+  data.tar.gz: 2754ed973a428206f5d2fa4a2327a77fa43a5fedb7151be499f7aaa029d0ca50
 SHA512:
-  metadata.gz: ab3adea770a98b961b048e456dc8a1b88f69cdd11503f4fa4bef201c436d357b1021d0a7b8e16bd48202669f0b0785955be780312b08dcceb224d2bd9f4847b2
-  data.tar.gz: 29ed19c8aabeb5d6c2d9ca784dc7dea1c7f3f7540537735a72bf5e21af55fde8a296ad48cec08b0c0aa36b6d28c74ff0f653ac42d8456d5e9d35b6d379c50972
+  metadata.gz: b1a7250d3b3c71680547decb6062086889c5dc60b03ba32a0bff7331264696fb435131fb9239a8b9c37e0aec0e90e6e5dcc56156eb2b7424ec7013cbe364928d
+  data.tar.gz: f430213bf4698399eae2607930ba443c529da150df03034f2e3100e68a45d43d18d228a8e64fc74f63561518a89e26db04b653d9a0e751de70a38c059f887c38

data/Gemfile.lock

@@ -1,31 +1,38 @@
 PATH
   remote: .
   specs:
-    panomosity (0.1.36)
+    panomosity (0.1.39)
+      write_xlsx
 
 GEM
   remote: https://rubygems.org/
   specs:
     coderay (1.1.2)
-    diff-lcs (1.3)
+    diff-lcs (1.4.4)
     method_source (0.9.0)
     pry (0.11.3)
       coderay (~> 1.1.0)
       method_source (~> 0.9.0)
     rake (10.5.0)
-    rspec (3.8.0)
-      rspec-core (~> 3.8.0)
-      rspec-expectations (~> 3.8.0)
-      rspec-mocks (~> 3.8.0)
-    rspec-core (3.8.0)
-      rspec-support (~> 3.8.0)
-    rspec-expectations (3.8.1)
+    rspec (3.9.0)
+      rspec-core (~> 3.9.0)
+      rspec-expectations (~> 3.9.0)
+      rspec-mocks (~> 3.9.0)
+    rspec-core (3.9.2)
+      rspec-support (~> 3.9.3)
+    rspec-expectations (3.9.2)
       diff-lcs (>= 1.2.0, < 2.0)
-      rspec-support (~> 3.8.0)
-    rspec-mocks (3.8.0)
+      rspec-support (~> 3.9.0)
+    rspec-mocks (3.9.1)
       diff-lcs (>= 1.2.0, < 2.0)
-      rspec-support (~> 3.8.0)
-    rspec-support (3.8.0)
+      rspec-support (~> 3.9.0)
+    rspec-support (3.9.3)
+    rubyzip (2.3.0)
+    write_xlsx (0.85.9)
+      rubyzip (>= 1.0.0)
+      zip-zip
+    zip-zip (0.3)
+      rubyzip (>= 1.0.0)
 
 PLATFORMS
   ruby
@@ -38,4 +45,4 @@ DEPENDENCIES
   rspec (~> 3.0)
 
 BUNDLED WITH
-   1.16.1
+   1.17.3

data/lib/panomosity.rb

@@ -1,5 +1,7 @@
 require 'panomosity/control_point'
+require 'panomosity/generalized_neighborhood'
 require 'panomosity/image'
+require 'panomosity/measure'
 require 'panomosity/neighborhood'
 require 'panomosity/neighborhood_group'
 require 'panomosity/optimisation_variable'
@@ -9,10 +11,12 @@ require 'panomosity/panorama'
 require 'panomosity/panorama_variable'
 require 'panomosity/runner'
 require 'panomosity/utils'
+require 'panomosity/errors'
 require 'panomosity/version'
 require 'pathname'
 require 'fileutils'
 require 'optparse'
+require 'json'
 
 module Panomosity
   def self.parse(arguments)
@@ -79,6 +83,33 @@ module Panomosity
         exit
       end
 
+      parser.on('--regional-distance-similarities-count COUNT', Integer, 'Set the minimum amount of regional control point counts for determining similar neighborhoods (default: 3)') do |count|
+        options[:regional_distance_similarities_count] = count
+      end
+
+      parser.on('--max-reduction-attempts COUNT', Integer, 'Set the max reduction attempts when removing neighborhood outliers (default: 2)') do |count|
+        options[:max_reduction_attempts] = count
+      end
+
+      desc = <<~DESC
+        Set distances to use when determining neighborhood region size in pairs
+        Use JSON e.g. '{"x1": 150, "x2": 30}'
+        Defaults:
+          Vertical pair is x is 10% of image width and y is 100px 
+          Horizontal pair is x is 100px and y is 10% of image height
+      DESC
+      parser.on('--distances [DISTANCE_JSON]', desc) do |distances|
+        options[:distances] = JSON.parse(distances) rescue nil
+      end
+
+      parser.on('--distances-horizontal [DISTANCE_JSON]', 'Same as above but only affects horizontal image pairs') do |distances|
+        options[:regional_distance_similarities_count] = JSON.parse(distances) rescue nil
+      end
+
+      parser.on('--distances-vertical [DISTANCE_JSON]', 'Same as above but only affects vertical image pairs') do |distances|
+        options[:regional_distance_similarities_count] = JSON.parse(distances) rescue nil
+      end
+
       parser.on('-h', '--help', 'Display this screen') do
         puts parser
         exit

data/lib/panomosity/control_point.rb

@@ -196,8 +196,8 @@ module Panomosity
       "#{to_s.sub(/\n/, '')} dist #{dist.round(4)} pixel_dist #{px.round(4)},#{py.round(4)},#{pdist.round(4)} pixel_r_dist #{prx.round(4)},#{pry.round(4)},#{prdist.round(4)} conn_type #{conn_type}"
     end
 
-    def attributes
-      @attributes
+    def attributes(raw: false)
+      @attributes.keep_if { |k, _| raw || !%i(raw i1 i2).include?(k) }
     end
   end
 end

data/lib/panomosity/errors.rb

@@ -0,0 +1,3 @@
+module Panomosity
+  class NoSimilarNeighborhoodsError < StandardError; end
+end

data/lib/panomosity/generalized_neighborhood.rb

@@ -0,0 +1,278 @@
+require 'panomosity/utils'
+
+module Panomosity
+  class GeneralizedNeighborhood
+    include Panomosity::Utils
+
+    CONTROL_POINT_ATTRIBUTES = [:dist_avg, :dist_std, :x_avg, :x_std, :y_avg, :y_std]
+    ATTRIBUTES = [:center, :scope, :control_points, :count] + CONTROL_POINT_ATTRIBUTES
+    DEFAULT_DISTANCE = 100
+
+    attr_reader :measure, :options
+    attr_accessor *ATTRIBUTES
+
+    class << self
+      include Panomosity::Utils
+
+      attr_accessor :options, :neighborhoods, :horizontal_similar_neighborhoods, :vertical_similar_neighborhoods, :horizontal_neighborhoods_by_similar_neighborhood, :vertical_neighborhoods_by_similar_neighborhood
+
+      def logger
+        @logger ||= Panomosity.logger
+      end
+
+      def attributes
+        { name: name }
+      end
+
+      def horizontal
+        horizontal_neighborhoods_by_similar_neighborhood
+      end
+
+      def vertical
+        vertical_neighborhoods_by_similar_neighborhood
+      end
+
+      def calculate_all(panorama:, options: {})
+        @neighborhoods = []
+        @options = options
+
+        Pair.create_pairs_from_panorama(panorama)
+        calculate_from_pairs
+        calculate_from_neighborhoods(type: :horizontal)
+        calculate_from_neighborhoods(type: :vertical)
+
+        @neighborhoods
+      end
+
+      def calculate_from_pairs
+        logger.debug 'calculating neighborhoods from pairs'
+        Pair.all.each do |pair|
+          control_points = pair.control_points.select(&:not_generated?)
+          control_points.each do |control_point|
+            base_params = { center: control_point, scope: pair, options: options }
+            position_params = { measure: { type: :position } }
+            neighborhood = calculate_neighborhood(**base_params.merge(position_params))
+            pair.generalized_neighborhoods << neighborhood
+            distance_params = { measure: { type: :distance, distances: { x1: neighborhood.dist_std } } }
+            distance_neighborhood = calculate_neighborhood(**base_params.merge(distance_params))
+            neighborhood.reference = distance_neighborhood
+            pair.generalized_neighborhoods << distance_neighborhood
+          end
+        end
+      end
+
+      def calculate_from_neighborhoods(type:)
+        count = options[:regional_distance_similarities_count] || 3
+        attempts = options[:max_reduction_attempts] || 2
+
+        # calculates similar neighborhoods based on the regional control point distances by pair
+        calculate_similar_neighborhoods(type: type, count: count)
+        if neighborhoods.empty?
+          logger.warn 'total neighborhoods came up empty, neighborhood default count to 2'
+          calculate_similar_neighborhoods(type: type, count: 2)
+          raise 'still could not find neighborhoods' if neighborhoods.empty?
+        end
+
+        std_outlier_reduction(type: type, max_reduction_attempts: attempts)
+        calculate_neighborhoods_by_similar_neighborhood(type: type)
+      end
+
+      def calculate_neighborhood(center:, scope:, options:, measure: {})
+        neighborhood = new(center: center, scope: scope, options: options)
+        neighborhood.update_measure(measure)
+        @neighborhoods << neighborhood.calculate
+        neighborhood
+      end
+
+      def similar_neighborhoods(type: :horizontal)
+        type == :horizontal ? @horizontal_similar_neighborhoods : @vertical_similar_neighborhoods
+      end
+
+      def similar_neighborhoods!(type: :horizontal)
+        neighborhoods = similar_neighborhoods(type: type)
+
+        if neighborhoods.nil? || neighborhoods.empty?
+          error = "No similar #{type} neighborhoods found"
+          raise NoSimilarNeighborhoodsError, error
+        else
+          neighborhoods
+        end
+      end
+
+      def neighborhoods_by_similar_neighborhood(type: :horizontal)
+        type == :horizontal ? @horizontal_neighborhoods_by_similar_neighborhood : @vertical_neighborhoods_by_similar_neighborhood
+      end
+
+      def calculate_similar_neighborhoods(type: :horizontal, count: 3)
+        similar_neighborhoods = neighborhoods.select(&:measure_position?).select(&:"#{type}?").select do |neighborhood|
+          neighborhood.scope.similar_neighborhoods << neighborhood if neighborhood.reference.count >= count
+        end
+        self.send(:"#{type}_similar_neighborhoods=", similar_neighborhoods)
+      end
+
+      def std_outlier_reduction(type: :horizontal, max_reduction_attempts: 2, reduction_attempts: 0)
+        return if reduction_attempts >= max_reduction_attempts
+
+        logger.debug "twice reducing #{type} neighborhood std outliers"
+
+        neighborhoods = similar_neighborhoods!(type: type)
+        std_dist_of_neighborhoods = neighborhoods.map(&:dist_std)
+        avg, std = *calculate_average_and_std(values: std_dist_of_neighborhoods)
+
+        similar_neighborhoods!(type: type).select! { |n| (avg - n.dist_std).abs <= std }
+        std_outlier_reduction(type: type, max_reduction_attempts: max_reduction_attempts, reduction_attempts: reduction_attempts + 1)
+      end
+
+      def calculate_neighborhoods_by_similar_neighborhood(type: :horizontal)
+        instance_variable_set("@#{type}_neighborhoods_by_similar_neighborhood", [])
+        similar_neighborhoods(type: type).each do |neighborhood|
+          base_params = { center: neighborhood, scope: self, options: options }
+          distance_params = { measure: { type: :distance, distances: { x1: neighborhood.dist_std } } }
+          neighborhoods_by_similar_neighborhood(type: type) << calculate_neighborhood(**base_params.merge(distance_params))
+        end
+
+        neighborhoods_by_similar_neighborhood(type: type).sort_by! { |n| -n.count }
+        neighborhoods_by_similar_neighborhood(type: type).max_by(5) { |n| n.count }.each do |n|
+          logger.debug "#{n.dist_avg} #{n.dist_std} #{n.count} x#{n.x_avg} y#{n.y_avg}"
+        end
+      end
+    end
+
+    def initialize(center:, scope:, options: {})
+      @center = center
+      @scope = scope
+      @options = options
+      @measure = Measure.new
+    end
+
+    def id
+      @id
+    end
+
+    def id=(id)
+      @id = id
+    end
+
+    def reference
+      @reference
+    end
+
+    def reference=(reference)
+      @reference = reference
+    end
+
+    def pair_scope?
+      scope.class.name == 'Panomosity::Pair'
+    end
+
+    def neighborhood_scope?
+      scope.class.name == self.class.name
+    end
+
+    def measure_position?
+      measure.position?
+    end
+
+    def measure_distance?
+      measure.distance?
+    end
+
+    def horizontal?
+      pair_scope? ? scope.horizontal? : center.horizontal?
+    end
+
+    def vertical?
+      pair_scope? ? scope.vertical? : center.vertical?
+    end
+
+    def type
+      horizontal? ? :horizontal : :vertical
+    end
+
+    def update_measure(params)
+      measure.update(params)
+      set_distance_defaults
+      set_measure_defaults
+    end
+
+    def distances_from_options(type: :horizontal)
+      if type == :both
+        options.fetch(:distances, {})
+      else
+        options[:distances]&.fetch(type, {}) || {}
+      end
+    end
+
+    def set_distance_defaults
+      return unless measure_position?
+
+      measure.update_distances(distances_from_options(type: :both))
+
+      if scope.horizontal?
+        measure.update_distances(x2: (scope.first_image.h * 0.1).round)
+        measure.update_distances(distances_from_options(type: :horizontal))
+      else
+        measure.update_distances(x1: (scope.first_image.w * 0.1).round)
+        measure.update_distances(distances_from_options(type: :vertical))
+      end
+
+      measure.distances[:x1] ||= DEFAULT_DISTANCE
+      measure.distances[:x2] ||= DEFAULT_DISTANCE
+    end
+
+    def set_measure_defaults
+      center_values = if measure_position?
+        { x1: center.x1, x2: center.y1 }
+      elsif pair_scope?
+        { x1: center.pdist }
+      else
+        { x1: center.dist_avg }
+      end
+
+      measure.update(center: center_values)
+    end
+
+    def calculate
+      if pair_scope?
+        elements = scope.control_points.select(&:not_generated?)
+
+        @control_points = if measure_position?
+          elements.select { |cp| measure.includes?(cp.x1, cp.y1) }
+        else
+          elements.select { |cp| measure.includes?(cp.pdist) }
+        end
+      else
+        @neighborhoods = scope.similar_neighborhoods(type: center.type).select { |n| measure.includes?(n.dist_avg) }
+        distance_neighborhoods = @neighborhoods.map(&:reference)
+        @control_points = distance_neighborhoods.map(&:control_points).flatten.uniq(&:raw)
+      end
+
+      @dist_avg, @dist_std = *calculate_average_and_std(values: control_points.map(&:pdist), ignore_empty: true)
+      @x_avg, @x_std = *calculate_average_and_std(values: control_points.map(&:px), ignore_empty: true)
+      @y_avg, @y_std = *calculate_average_and_std(values: control_points.map(&:py), ignore_empty: true)
+      @count = control_points.count
+
+      self
+    end
+
+    def attributes
+      attributes = CONTROL_POINT_ATTRIBUTES.reduce({}) { |h, k| h.merge!({ k => self.send(k) }) }
+      measure_attributes = measure.attributes.reduce({}) do |hash, (key, value)|
+        hash["measure_#{key}"] = value
+        hash
+      end
+      if pair_scope?
+        center_id = center[:id]
+        scope_id = [scope.control_points.first.n1, scope.control_points.first.n2]
+        scope_name = 'pair'
+      else
+        center_id = center.id
+        scope_id = nil
+        scope_name = 'neighborhood'
+      end
+      attributes.merge!(measure_attributes)
+      attributes.merge!(id: id, center: center_id, scope_id: scope_id, scope_name: scope_name, type: type, control_points: control_points.map{|c| c[:id]})
+      attributes
+    end
+  end
+end

data/lib/panomosity/image.rb

@@ -1,4 +1,4 @@
-# Exmaple of an image line
+# Example of an image line
 # i w2448 h3264 f0 v=0 Ra=0 Rb=0 Rc=0 Rd=0 Re=0 Eev6.433 Er1 Eb1 r0 p0 y0 TrX0.88957152 TrY0.79560269 TrZ1 Tpy0 Tpp0 j0 a=0 b=0 c=0 d=0 e=0 g0 t0 Va=0 Vb=0 Vc=0 Vd=0 Vx=0 Vy=0  Vm5 n"WZ8ppTx9PtcxASB3hbeeuS6Z"\n
 
 module Panomosity
@@ -167,5 +167,11 @@ module Panomosity
       return unless name =~ /c(\d+)_r(\d+)\.jpg/
       name.scan(/c(\d+)_r(\d+)\.jpg/).flatten.map(&:to_i).last
     end
+
+    def attributes
+      attributes = @attributes.keep_if { |k, _| !%i(raw).include?(k) }
+      attributes.merge!(column: column || 0, row: row || 0)
+      attributes
+    end
   end
 end

data/lib/panomosity/measure.rb

@@ -0,0 +1,56 @@
+require 'panomosity/utils'
+
+module Panomosity
+  class Measure
+    include Panomosity::Utils
+
+    ATTRIBUTES = [:type, :center, :distances]
+    DEFAULT_DISTANCE = 100
+
+    attr_reader :options
+    attr_accessor *ATTRIBUTES
+
+    def initialize(attributes = {})
+      @attributes = attributes
+      @attributes[:center] ||= {}
+      @attributes[:distances] ||= {}
+    end
+
+    def position?
+      type == :position
+    end
+
+    def distance?
+      type == :distance
+    end
+
+    def type
+      @attributes[:type]
+    end
+
+    def center
+      @attributes[:center]
+    end
+
+    def distances
+      @attributes[:distances]
+    end
+
+    def update(attributes = {})
+
+      @attributes.merge!(attributes)
+    end
+
+    def update_distances(values = {})
+      @attributes[:distances].merge!(values)
+    end
+
+    def includes?(*values)
+      center.values.zip(values, distances.values).all? { |center, value, distance| (center - value).abs <= distance }
+    end
+
+    def attributes
+      { type: type, center: center.values, distances: distances.values }
+    end
+  end
+end

data/lib/panomosity/optimisation_variable.rb

@@ -59,12 +59,12 @@ module Panomosity
     end
 
     def to_s
-      line_values = @@attributes.map { |attribute| "#{attribute}#{self.send(attribute)}" }
-      "v #{line_values.join(' ')}\n"
+      line_values = @@attributes.map { |attribute| "#{attribute}#{self.send(attribute)}" if self.send(attribute) }
+      "v #{line_values.compact.join(' ')}\n"
     end
 
     def attributes
-      @attributes.keep_if { |k,_| !%i(raw).include?(k) }
+      @attributes.keep_if { |k, _| !%i(raw).include?(k) }
     end
   end
 end

data/lib/panomosity/optimizer.rb

@@ -30,17 +30,14 @@ module Panomosity
         logger.info "applying custom values of xh_avg: #{xh_avg}, yh_avg: #{yh_avg}, xv_avg: #{xv_avg}, yv_avg: #{yv_avg}"
       end
 
-      unless xh_avg && yh_avg && xv_avg && yv_avg
-        Pair.calculate_neighborhoods(panorama)
-        Pair.calculate_neighborhood_groups
-      end
+      panorama.calculate_neighborhoods unless xh_avg && yh_avg && xv_avg && yv_avg
 
       ds = images.map(&:d).uniq.sort
       es = images.map(&:e).uniq.sort
 
       # get the average error for the best neighborhood group
-      x_avg = xh_avg || NeighborhoodGroup.horizontal.first.rx_avg
-      y_avg = yv_avg || NeighborhoodGroup.vertical.first.ry_avg
+      x_avg = xh_avg || GeneralizedNeighborhood.horizontal.first.x_avg
+      y_avg = yv_avg || GeneralizedNeighborhood.vertical.first.y_avg
 
       # start horizontally
       d_map = {}
@@ -57,8 +54,8 @@ module Panomosity
       logger.debug "created e_map #{e_map}"
 
       # add in the other offset
-      x_avg = xv_avg || NeighborhoodGroup.vertical.first.rx_avg
-      y_avg = yh_avg || NeighborhoodGroup.horizontal.first.ry_avg
+      x_avg = xv_avg || GeneralizedNeighborhood.vertical.first.x_avg
+      y_avg = yh_avg || GeneralizedNeighborhood.horizontal.first.y_avg
 
       de_map = {}
       d_map.each_with_index do |(dk,dv),di|
@@ -142,10 +139,9 @@ module Panomosity
     end
 
     def calculate_average_distance
-      Pair.calculate_neighborhoods(panorama)
-      Pair.calculate_neighborhood_groups
-      horizontal_distances = NeighborhoodGroup.horizontal[0..4].map(&:prdist_avg)
-      vertical_distances = NeighborhoodGroup.vertical[0..4].map(&:prdist_avg)
+      panorama.calculate_neighborhoods
+      horizontal_distances = GeneralizedNeighborhood.horizontal[0..4].map(&:dist_avg)
+      vertical_distances = GeneralizedNeighborhood.vertical[0..4].map(&:dist_avg)
       calculate_average(values: horizontal_distances + vertical_distances)
     end
 
@@ -161,9 +157,8 @@ module Panomosity
       images.each do |image|
         image.r = 0.0
       end
-      Pair.calculate_neighborhoods(panorama)
-      Pair.calculate_neighborhood_groups
-      cps = NeighborhoodGroup.horizontal.first.control_points
+      panorama.calculate_neighborhoods
+      cps = GeneralizedNeighborhood.horizontal.first.control_points
       avg, std = *calculate_average_and_std(values: cps.map(&:roll))
       # 0.1 degrees of std
       while std > 0.1

data/lib/panomosity/pair.rb

@@ -5,7 +5,7 @@ module Panomosity
     include Panomosity::Utils
     extend Panomosity::Utils
 
-    attr_accessor :pair, :control_points, :neighborhoods, :type
+    attr_accessor :pair, :control_points, :neighborhoods, :type, :generalized_neighborhoods, :similar_neighborhoods
 
     class << self
       attr_accessor :panorama, :logger
@@ -22,8 +22,8 @@ module Panomosity
         @pairs
       end
 
-      def good_control_points_to_keep(count: 3)
-        @pairs.map { |pair| pair.good_control_points_to_keep(count: count) }.flatten.uniq(&:raw)
+      def select_control_points_with_regional_distance_similarities
+        @pairs.map(&:select_control_points_with_regional_distance_similarities).flatten.uniq(&:raw)
       end
 
       def unconnected
@@ -44,6 +44,9 @@ module Panomosity
       rows = images.map(&:row).uniq.sort
 
       @pairs = []
+      @horizontal_pairs = []
+      @vertical_pairs = []
+      
       # horizontal pair creation
       rows.each do |row|
         columns.each do |column|
@@ -52,7 +55,9 @@ module Panomosity
           image_2 = images.find { |i| i.row == row && i.column == column.next }
           next if @panorama.calibration? && (image_1.nil? || image_2.nil?)
           control_points = @panorama.control_points.select { |cp| [cp.n1, cp.n2].sort == [image_1.id, image_2.id].sort }
-          @pairs << Pair.new([image_1, image_2].sort_by(&:id), control_points: control_points, type: :horizontal)
+          pair = Pair.new([image_1, image_2].sort_by(&:id), control_points: control_points, type: :horizontal)
+          @pairs << pair
+          @horizontal_pairs << pair
         end
       end
 
@@ -64,7 +69,9 @@ module Panomosity
           image_2 = images.find { |i| i.column == column && i.row == row.next }
           next if @panorama.calibration? && (image_1.nil? || image_2.nil?)
           control_points = @panorama.control_points.select { |cp| [cp.n1, cp.n2].sort == [image_1.id, image_2.id].sort }
-          @pairs << Pair.new([image_1, image_2].sort_by(&:id), control_points: control_points, type: :vertical)
+          pair = Pair.new([image_1, image_2].sort_by(&:id), control_points: control_points, type: :vertical)
+          @pairs << pair
+          @vertical_pairs << pair
         end
       end
     end
@@ -150,6 +157,8 @@ module Panomosity
       @pair = pair
       @control_points = control_points
       @neighborhoods = []
+      @generalized_neighborhoods = []
+      @similar_neighborhoods = []
       @type = type
     end
 
@@ -157,6 +166,10 @@ module Panomosity
       pair.map(&:id).to_s.gsub(' ', '')
     end
 
+    def ==(other)
+      to_s == other.to_s
+    end
+
     def info
       "#{to_s}(#{type}) image_1 d,e: #{pair.first.d},#{pair.first.e} | image_2 d,e: #{pair.last.d},#{pair.last.e}"
     end
@@ -218,6 +231,27 @@ module Panomosity
       end
     end
 
+    # Distance neighborhoods
+    def similar_control_points
+      @similar_control_points ||= similar_neighborhoods.map(&:reference).map(&:control_points).flatten.uniq(&:raw)
+    end
+
+    def select_control_points_with_regional_distance_similarities
+      # Keep all our control points if we have less than 10
+      if control_points.count >= 10
+        ratio = similar_control_points.count.to_f / control_points.count
+        if ratio < 0.2
+          Panomosity.logger.warn "#{to_s} keeping less than 20% (#{(ratio * 100).round(4)}%) of #{control_points.count} control points. Reverting and keeping all control points"
+          control_points
+        else
+          similar_control_points
+        end
+      else
+        Panomosity.logger.debug "Skipping pair #{to_s} since it has fewer than 10 control points"
+        control_points
+      end
+    end
+
     def best_neighborhood
       @best_neighborhood ||= @neighborhoods.max_by { |n| n.control_points.count }
     end
@@ -225,5 +259,17 @@ module Panomosity
     def control_points_of_best_neighborhood
       best_neighborhood ? best_neighborhood.control_points : []
     end
+
+    def attributes
+      x_avg, x_std = *calculate_average_and_std(values: control_points.map(&:px), ignore_empty: true)
+      y_avg, y_std = *calculate_average_and_std(values: control_points.map(&:py), ignore_empty: true)
+      dist_avg, dist_std = *calculate_average_and_std(values: control_points.map(&:pdist), ignore_empty: true)
+      i1 = control_points.first.n1
+      i2 = control_points.first.n2
+      {
+        id: [i1, i2], n: i1, N: i2, count: control_points.count, type: type,
+        x_avg: x_avg, x_std: x_std, y_avg: y_avg, y_std: y_std, dist_avg: dist_avg, dist_std: dist_std
+      }
+    end
   end
 end

data/lib/panomosity/panorama.rb

@@ -18,22 +18,25 @@ module Panomosity
       @logger = Panomosity.logger
     end
 
+    def calculate_neighborhoods
+      GeneralizedNeighborhood.calculate_all(panorama: self, options: options)
+    end
+
     def clean_control_points
-      if calibration?
-        Pair.calculate_neighborhoods(self, distance: 30)
-      else
-        Pair.calculate_neighborhoods(self, distance: 100)
-      end
-      control_points_to_keep = Pair.good_control_points_to_keep(count: 2)
+      options.merge!(distances: { x1: 30, x2: 30 }) if calibration? && !options[:distances].nil?
+      options.merge!(regional_distance_similarities_count: 2) unless options[:regional_distance_similarities_count]
+      calculate_neighborhoods
+
+      control_points_to_keep = Pair.select_control_points_with_regional_distance_similarities
       bad_control_points = control_points.reject { |cp| control_points_to_keep.map(&:raw).include?(cp.raw) }
       # far_control_points = control_points.select { |cp| cp.prdist > 50 }
       control_points_to_clean = bad_control_points.uniq(&:raw)
 
       # log warnings
       control_point_ratio = control_points_to_clean.count.to_f / control_points.count
-      logger.warn "Removing more than 30% (#{(control_point_ratio*100).round(4)}%) of control points. May potentially cause issues." if control_point_ratio >= 0.3
+      logger.warn "Removing more than 30% (#{(control_point_ratio * 100).round(4)}%) of control points. May potentially cause issues." if control_point_ratio >= 0.3
       control_point_pair_ratio = Pair.without_enough_control_points(ignore_connected: true).count.to_f / Pair.all.count
-      logger.warn "More than 50% (#{(control_point_pair_ratio*100).round(4)}%) of pairs have fewer than 3 control points. May potentially cause issues." if control_point_pair_ratio >= 0.5
+      logger.warn "More than 50% (#{(control_point_pair_ratio * 100).round(4)}%) of pairs have fewer than 3 control points. May potentially cause issues." if control_point_pair_ratio >= 0.5
 
       control_points_to_clean
     end
@@ -41,12 +44,8 @@ module Panomosity
     def fix_unconnected_image_pairs
       logger.info 'finding unconnected image pairs'
 
-      if calibration?
-        Pair.calculate_neighborhoods(self, distance: 30)
-      else
-        Pair.calculate_neighborhoods(self, distance: 100)
-      end
-      Pair.calculate_neighborhood_groups
+      options.merge!(distances: { x1: 30, x2: 30 }) if calibration? && !options[:distances].nil?
+      calculate_neighborhoods
 
       unconnected_image_pairs = Pair.unconnected
       logger.debug unconnected_image_pairs.map { |i| { type: i.type, pair: i.pair.map(&:id) } }
@@ -92,19 +91,19 @@ module Panomosity
     def generate_control_points(pair: nil, bad_control_point: nil, message: '')
       if pair
         if pair.horizontal?
-          group = NeighborhoodGroup.horizontal.first
+          group = GeneralizedNeighborhood.horizontal.first
         else
-          group = NeighborhoodGroup.vertical.first
+          group = GeneralizedNeighborhood.vertical.first
         end
       else
         if bad_control_point.conn_type == :horizontal
-          group = NeighborhoodGroup.horizontal.first
+          group = GeneralizedNeighborhood.horizontal.first
         else
-          group = NeighborhoodGroup.vertical.first
+          group = GeneralizedNeighborhood.vertical.first
         end
       end
 
-      control_point = ControlPoint.new(group.center.center.attributes)
+      control_point = ControlPoint.new(group.center.center.attributes(raw: true))
 
       if pair
         control_point[:n] = pair.first_image.id
@@ -153,8 +152,7 @@ module Panomosity
     end
 
     def diagnose
-      Pair.calculate_neighborhoods(self)
-      Pair.calculate_neighborhood_groups
+      calculate_neighborhoods
 
       recommendations = []
       messages = []
@@ -184,7 +182,7 @@ module Panomosity
       end
 
       # neighborhood group tests
-      group_count = NeighborhoodGroup.horizontal.count
+      group_count = GeneralizedNeighborhood.horizontal.count
       if group_count < 5
         message = <<~MESSAGE
           Total number of horizontal neighborhood groups is #{group_count} which is very low.
@@ -194,7 +192,7 @@ module Panomosity
         messages << message
       end
 
-      group_std_avg = calculate_average(values: NeighborhoodGroup.horizontal[0..4].map(&:prdist_std))
+      group_std_avg = calculate_average(values: GeneralizedNeighborhood.horizontal[0..4].map(&:dist_std))
       if group_std_avg > 1.0
         message = <<~MESSAGE
           The standard deviation of distances in the top 5 horizontal neighborhood groups is #{group_std_avg} which is high.
@@ -207,7 +205,7 @@ module Panomosity
         messages << message
       end
 
-      group_control_points = NeighborhoodGroup.horizontal.first.control_points.count
+      group_control_points = GeneralizedNeighborhood.horizontal.first.control_points.count
       total_control_points = Pair.horizontal.map(&:control_points).flatten.uniq(&:raw).count
       group_control_point_ratio = group_control_points.to_f / total_control_points
       if group_control_point_ratio < 0.2
@@ -222,7 +220,7 @@ module Panomosity
         recommendations << 'horizontal'
       end
 
-      group_count = NeighborhoodGroup.vertical.count
+      group_count = GeneralizedNeighborhood.vertical.count
       if group_count < 5
         message = <<~MESSAGE
           Total number of vertical neighborhood groups is #{group_count} which is very low.
@@ -232,7 +230,7 @@ module Panomosity
         messages << message
       end
 
-      group_std_avg = calculate_average(values: NeighborhoodGroup.vertical[0..4].map(&:prdist_std))
+      group_std_avg = calculate_average(values: GeneralizedNeighborhood.vertical[0..4].map(&:dist_std))
       if group_std_avg > 1.0
         message = <<~MESSAGE
           The standard deviation of distances in the top 5 vertical neighborhood groups is #{group_std_avg} which is high.
@@ -245,7 +243,7 @@ module Panomosity
         messages << message
       end
 
-      group_control_points = NeighborhoodGroup.vertical.first.control_points.count
+      group_control_points = GeneralizedNeighborhood.vertical.first.control_points.count
       total_control_points = Pair.vertical.map(&:control_points).flatten.uniq(&:raw).count
       group_control_point_ratio = group_control_points.to_f / total_control_points
       if group_control_point_ratio < 0.2
@@ -274,11 +272,24 @@ module Panomosity
           delta_d: delta_d,
           delta_e: delta_e,
           roll: roll,
-          horizontal: NeighborhoodGroup.horizontal.first.serialize,
-          vertical: NeighborhoodGroup.vertical.first.serialize
+          horizontal: GeneralizedNeighborhood.horizontal.first.attributes,
+          vertical: GeneralizedNeighborhood.vertical.first.attributes
         }
       }
+    rescue => error
+      message = "Got error #{error.message} when calculating neighborhoods. Recommending fallback"
 
+      logger.error message
+      error.backtrace.each { |line| logger.error line }
+
+      recommendations = %w(horizontal vertical)
+
+      diagnostic_report = {
+        messages: messages,
+        recommendations: recommendations,
+        data: {}
+      }
+    ensure
       File.open('diagnostic_report.json', 'w+') { |f| f.puts diagnostic_report.to_json }
 
       if recommendations.empty?
@@ -301,17 +312,13 @@ module Panomosity
       calibration_report = JSON.parse(File.read(filename))
 
       if @options[:report_type] == 'position'
-        if calibration?
-          Pair.calculate_neighborhoods(self, distance: 30)
-        else
-          Pair.calculate_neighborhoods(self, distance: 100)
-        end
-        Pair.calculate_neighborhood_groups
+        options.merge!(distances: { x1: 30, x2: 30 }) if calibration? && !options[:distances].nil?
+        calculate_neighborhoods
 
-        xh_avg = NeighborhoodGroup.horizontal.first.x_avg
-        yh_avg = NeighborhoodGroup.horizontal.first.y_avg
-        xv_avg = NeighborhoodGroup.vertical.first.x_avg
-        yv_avg = NeighborhoodGroup.vertical.first.y_avg
+        xh_avg = GeneralizedNeighborhood.horizontal.first.x_avg
+        yh_avg = GeneralizedNeighborhood.horizontal.first.y_avg
+        xv_avg = GeneralizedNeighborhood.vertical.first.x_avg
+        yv_avg = GeneralizedNeighborhood.vertical.first.y_avg
         calibration_report['position'] = {
           xh_avg: xh_avg,
           yh_avg: yh_avg,
@@ -327,7 +334,39 @@ module Panomosity
     end
 
     def calibration?
-      !!@input.split(/\n/).find { |line| '#panomosity calibration true' }
+      !!@input.split(/\n/).find { |line| line == '#panomosity calibration true' }
+    end
+
+    def save_file(filename)
+      logger.info "saving file #{filename}"
+
+      lines = @input.each_line.map do |line|
+        objects = [images, variable, control_points, optimisation_variables].flatten
+        object = objects.find { |object| object.raw == line }
+        object&.to_s || line
+      end.compact
+
+      File.open(filename, 'w') { |f| lines.each { |line| f.puts line } }
+    end
+
+    def attributes
+      calculate_neighborhoods
+      control_points = self.control_points.dup
+      control_points.each_with_index { |cp, i| cp[:id] = i }
+      neighborhoods = GeneralizedNeighborhood.neighborhoods.dup
+      neighborhoods.each_with_index { |n, i| n.id = i }
+      types = %i(horizontal vertical)
+      similar_neighborhoods = types.map { |type| GeneralizedNeighborhood.similar_neighborhoods(type: type) }.flatten
+      neighborhoods_by_similar_neighborhood = types.map { |type| GeneralizedNeighborhood.neighborhoods_by_similar_neighborhood(type: type) }.flatten
+      {
+        images: images.map(&:attributes),
+        variable: variable.attributes,
+        control_points: control_points.map(&:attributes),
+        optimisation_variables: optimisation_variables.map(&:attributes),
+        pairs: Pair.all.map(&:attributes),
+        similar_neighborhoods: similar_neighborhoods.map(&:attributes),
+        neighborhoods_by_similar_neighborhood: neighborhoods_by_similar_neighborhood.map(&:attributes)
+      }
     end
   end
 end

data/lib/panomosity/panorama_variable.rb

@@ -60,5 +60,9 @@ module Panomosity
       line_values = @@attributes.map { |attribute| "#{attribute}#{self.send(attribute)}" }
       "p #{line_values.join(' ')}\n"
     end
+
+    def attributes
+      @attributes.keep_if { |k, _| !%i(raw).include?(k) }
+    end
   end
 end

data/lib/panomosity/runner.rb

@@ -1,14 +1,16 @@
 require 'logger'
 require 'json'
 require 'csv'
+require 'panomosity/xlsx_writer'
 
 module Panomosity
   class Runner
     include Panomosity::Utils
+    include Panomosity::XLSXWriter
 
     attr_reader :logger
 
-    AVAILABLE_COMMANDS = %w(
+    AVAILABLE_COMMANDS = %w[
       add_calibration_flag
       apply_calibration_values
       check_position_changes
@@ -20,12 +22,14 @@ module Panomosity
       crop_centers
       diagnose
       export_control_point_csv
+      export_panorama_attributes_to_csv
       fix_conversion_errors
       fix_unconnected_image_pairs
       generate_border_line_control_points
       get_columns_and_rows
       get_control_point_info
       get_neighborhood_info
+      get_panorama_attributes
       import_control_point_csv
       merge_image_parameters
       nona_grid
@@ -33,7 +37,7 @@ module Panomosity
       prepare_for_pr0ntools
       remove_anchor_variables
       standardize_roll
-    )
+    ]
 
     def initialize(options)
       @options = options
@@ -226,12 +230,19 @@ module Panomosity
       logger.debug "cropping to #{percent}%"
       unless @options[:without_cropping]
         images = Image.parse(@input_file)
-        images.each do |image|
+
+        crop_center_commands = images.map do |image|
           width_offset = (image.w.to_f * (1 - scale_factor) / 2).round
           height_offset = (image.h.to_f * (1 - scale_factor) / 2).round
           logger.debug "cropping #{image.name}"
-          `convert #{image.name} -crop "#{percent}%x#{percent}%+#{width_offset}+#{height_offset}" #{image.name}`
+
+          "convert #{image.name} -crop '#{percent}%x#{percent}%+#{width_offset}+#{height_offset}' #{image.name}"
         end
+
+        crop_centers_command = crop_center_commands.join("\n")
+        parallel_command = "cat <<'EOF' | parallel \n%s\nEOF"
+        crop_centers_command = parallel_command % crop_centers_command
+        `#{crop_centers_command}`
       end
 
       # Since all images have been cropped, we need to change d,e params to move images based on how much was cropped
@@ -295,10 +306,10 @@ module Panomosity
     def export_control_point_csv
       logger.info 'exporting control point csv'
       panorama = Panorama.new(@input_file, @options)
-      Pair.calculate_neighborhoods(panorama, distance: 30)
+      panorama.calculate_neighborhoods
 
       filename = 'control_points.csv'
-      headers = %w(image_1_name image_2_name image_1_id image_2_id type width height d1 e1 d2 e2 x1 y1 x2 y2 rx ry r)
+      headers = %w[image_1_name image_2_name image_1_id image_2_id type width height d1 e1 d2 e2 x1 y1 x2 y2 rx ry r]
       CSV.open(filename, 'w+') do |csv|
         csv << headers
         panorama.control_points.each do |cp|
@@ -307,7 +318,7 @@ module Panomosity
                   cp.i2.d, cp.i2.e, cp.x1, cp.y1, cp.x2, cp.y2, cp.px, cp.py, cp.pdist]
         end
       end
-      
+
       logger.info 'Done. Check for control_points.csv'
     end
 
@@ -481,6 +492,78 @@ module Panomosity
       panorama.get_neighborhood_info
     end
 
+    def get_panorama_attributes
+      # for printing out the JSON
+      logger.level = Logger::UNKNOWN
+      logger.info 'getting panorama attributes'
+      panorama = Panorama.new(@input_file, @options)
+      puts panorama.attributes.to_json
+    end
+
+    def export_panorama_attributes_to_csv
+      panorama = Panorama.new(@input_file, @options)
+      panorama.calculate_neighborhoods
+
+      logger.info "Making #{Dir.pwd}/#{@input}_attributes Directory"
+
+      dir = "#{Dir.pwd}/#{@input}_attributes"
+      Dir.mkdir(dir) unless File.exists?(dir)
+
+      logger.info "Done. Moving into #{dir} Directory"
+
+      Dir.chdir "#{@input}_attributes"
+
+      logger.info "Making control_points.csv"
+
+      filename = 'control_points.csv'
+      headers = %w[image_1_name image_2_name image_1_id image_2_id type width height d1 e1 d2 e2 x1 y1 x2 y2 rx ry r]
+      CSV.open(filename, 'w+') do |csv|
+        csv << headers
+        panorama.control_points.each do |cp|
+          pair = Pair.all.find { |p| p.control_points.include?(cp) }
+          csv << [cp.i1.name, cp.i2.name, cp.i1.id, cp.i2.id, pair.type, cp.i1.w, cp.i1.h, cp.i1.d, cp.i1.e,
+                  cp.i2.d, cp.i2.e, cp.x1, cp.y1, cp.x2, cp.y2, cp.px, cp.py, cp.pdist]
+        end
+      end
+
+      logger.info "Done. Check for #{filename}"
+
+      attributes = panorama.attributes
+
+      logger.info "Making cp_pairs.csv"
+
+      filename = 'cp_pairs.csv'
+      headers = attributes[:pairs].first.keys
+      CSV.open(filename,"w") do |csv|
+        csv << headers
+        attributes[:pairs].each do |pair|
+          csv << pair.values
+        end
+      end
+
+      logger.info "Done. Check for #{filename}"
+
+      logger.info "Making neighborhood_pairs.csv"
+
+      filename = 'neighborhood_pairs.csv'
+      headers = ['data_type', attributes[:similar_neighborhoods].first.keys].flatten
+      CSV.open(filename,'w') do |csv|
+        csv << headers
+        attributes[:similar_neighborhoods].each do |neighborhood|
+          csv << ['similar_neighborhoods', neighborhood.values].flatten(1)
+        end
+        attributes[:neighborhoods_by_similar_neighborhood].each do |neighborhood|
+          csv << ['neighborhoods_by_similar_neighborhood', neighborhood.values].flatten(1)
+        end
+      end
+
+      logger.info "Done. Check for #{filename}"
+
+      csv_to_xlsx(['control_points.csv', 'cp_pairs.csv', 'neighborhood_pairs.csv'], 'attributes')
+
+      Dir.chdir '..'
+    end
+
     def import_control_point_csv
       logger.info 'importing control point csv'
       panorama = Panorama.new(@input_file, @options)
@@ -565,6 +648,7 @@ module Panomosity
         save_file
         logger.debug "running nona #{@output}"
         output = "nona --save-intermediate-images --intermediate-suffix=intermediate -v -m TIFF_m --seam=blend #{@output} -o #{res}_res_stitch_section_c#{column.to_s.rjust(5, '0')}_"
+        # output = "nona-mask #@output -m TIFF_m -o #{res}_res_stitch_section_c#{column.to_s.rjust(5, '0')}_"
         logger.debug output
         output
       end
@@ -601,12 +685,17 @@ module Panomosity
       logger.info 'preparing for pr0ntools'
 
       images = Image.parse(@input_file)
+      # ds --> X (cols)
+      # es --> Y (rows)
       ds = images.map(&:d).uniq.sort
       es = images.map(&:e).uniq.sort
 
-      # within two pixels of each other
-      d_groups = ds.map { |d| ds.select { |n| n.between?(d-2, d+2) }.sort }
-      e_groups = es.map { |e| es.select { |n| n.between?(e-2, e+2) }.sort }
+      # split ds and es into groups within two pixels of each other,
+      # required because d/e can vary a few pixels based on system precision
+      # each group of d/e represents an acceptable pixel range
+      # for images to be considered part of a col/row
+      d_groups = ds.map { |d| ds.select { |n| n.between?(d-2, d+2) }.sort }.uniq
+      e_groups = es.map { |e| es.select { |n| n.between?(e-2, e+2) }.sort }.uniq
 
       fov = images.map(&:v).uniq.find { |v| v != 0.0 }
       images.each do |i|
@@ -632,7 +721,12 @@ module Panomosity
       FileUtils.mkdir_p(new_dir)
 
       logger.info 'renaming files to format cX_rY.jpg'
-      images.each { |image| FileUtils.cp(image[:original_n], new_dir + image.n) }
+      images.each do |image|
+        original_name = image[:original_n]
+        new_name = new_dir + image.n
+        FileUtils.cp(original_name, new_name)
+        logger.debug "copied #{original_name} to #{new_name}"
+      end
 
       logger.info 'creating scan.json'
       image_first_col = images.find { |image| image.n == 'c0_r0.jpg' }

data/lib/panomosity/utils.rb

@@ -31,4 +31,4 @@ module Panomosity
       values
     end
   end
-end
+end

data/lib/panomosity/version.rb

@@ -1,3 +1,3 @@
 module Panomosity
-  VERSION = '0.1.36'
+  VERSION = '0.1.41'
 end

data/lib/panomosity/xlsx_writer.rb

@@ -0,0 +1,40 @@
+require 'write_xlsx'
+
+module Panomosity
+  module XLSXWriter
+    SEPERATE_CSV_COLUMNS = /,(?![^\[]*\])/.freeze
+
+
+    def csv_to_xlsx(input_files, output_file)
+      input_files = Array(input_files)
+
+      logger.info "Creating #{output_file}.xlsx"
+
+      workbook = WriteXLSX.new("#{output_file}.xlsx")
+
+      input_files.each do |csv|
+        logger.info "Creating #{csv} Worksheet"
+
+        worksheet = workbook.add_worksheet(csv)
+
+        seperate_csv_columns = /,(?![^\[]*\])/
+
+        File.open(csv, "r") do |f|
+          f.each_line do |rows|
+            cells = rows.split(SEPERATE_CSV_COLUMNS)
+
+            cells.each_with_index do |cell, column|
+              row = f.lineno - 1
+              data = cell.tr_s('"', '').strip
+              worksheet.write(row, column, data)
+            end
+          end
+        end
+      end
+
+      workbook.close
+
+      logger.info "Done. Check for #{output_file}.xlsx"
+    end
+  end
+end

data/panomosity.gemspec

@@ -6,7 +6,7 @@ require 'panomosity/version'
 Gem::Specification.new do |spec|
   spec.name          = 'panomosity'
   spec.version       = Panomosity::VERSION
-  spec.authors       = ['Oliver Garcia']
+  spec.authors       = ['Oliver Garcia', 'Joshua Stowers']
   spec.email         = ['ogarci5@gmail.com']
 
   spec.summary       = %q{Wrapper for the PTO file parsing needed for PanoTools.}
@@ -24,4 +24,6 @@ Gem::Specification.new do |spec|
   spec.add_development_dependency 'bundler', '~> 1.16'
   spec.add_development_dependency 'rake', '~> 10.0'
   spec.add_development_dependency 'rspec', '~> 3.0'
+  spec.add_runtime_dependency 'write_xlsx'
+
 end

metadata

@@ -1,14 +1,15 @@
 --- !ruby/object:Gem::Specification
 name: panomosity
 version: !ruby/object:Gem::Version
-  version: 0.1.36
+  version: 0.1.41
 platform: ruby
 authors:
 - Oliver Garcia
+- Joshua Stowers
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2019-06-20 00:00:00.000000000 Z
+date: 2020-10-30 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler
@@ -52,6 +53,20 @@ dependencies:
     - - "~>"
       - !ruby/object:Gem::Version
         version: '3.0'
+- !ruby/object:Gem::Dependency
+  name: write_xlsx
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '0'
 description: Custom scripts to help with PTO parsing and different strategies.
 email:
 - ogarci5@gmail.com
@@ -77,7 +92,10 @@ files:
 - exe/panomosity
 - lib/panomosity.rb
 - lib/panomosity/control_point.rb
+- lib/panomosity/errors.rb
+- lib/panomosity/generalized_neighborhood.rb
 - lib/panomosity/image.rb
+- lib/panomosity/measure.rb
 - lib/panomosity/neighborhood.rb
 - lib/panomosity/neighborhood_group.rb
 - lib/panomosity/optimisation_variable.rb
@@ -88,6 +106,7 @@ files:
 - lib/panomosity/runner.rb
 - lib/panomosity/utils.rb
 - lib/panomosity/version.rb
+- lib/panomosity/xlsx_writer.rb
 - panomosity.gemspec
 homepage: https://github.com/elevatesystems/panomosity
 licenses:
@@ -108,8 +127,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubyforge_project: 
-rubygems_version: 2.7.3
+rubygems_version: 3.0.8
 signing_key: 
 specification_version: 4
 summary: Wrapper for the PTO file parsing needed for PanoTools.