gpx2exif 0.0.4 → 0.1.0

data/Gemfile

@@ -1,7 +1,7 @@
 source "http://rubygems.org"
 
 gem 'nokogiri'
-gem 'mini_exiftool', git: 'git://github.com/akwiatkowski/mini_exiftool.git'
+gem 'mini_exiftool'
 gem 'builder'
 
 group :development do

data/Gemfile.lock

@@ -1,9 +1,3 @@
-GIT
-  remote: git://github.com/akwiatkowski/mini_exiftool.git
-  revision: e2ca8e0c3de9190bd08520d6c4204d510963eb54
-  specs:
-    mini_exiftool (1.5.1)
-
 GEM
   remote: http://rubygems.org/
   specs:
@@ -14,8 +8,9 @@ GEM
       bundler (~> 1.0)
       git (>= 1.2.5)
       rake
+    mini_exiftool (1.5.1)
     multi_json (1.3.6)
-    nokogiri (1.5.3)
+    nokogiri (1.5.5)
     rake (0.9.2.2)
     rspec (2.3.0)
       rspec-core (~> 2.3.0)
@@ -37,7 +32,7 @@ DEPENDENCIES
   builder
   bundler (~> 1.0.0)
   jeweler (~> 1.6.4)
-  mini_exiftool!
+  mini_exiftool
   nokogiri
   rspec (~> 2.3.0)
   simplecov

data/README.md

@@ -65,6 +65,34 @@ How to use it
 
   generate_garmin_waypoints -y samples/sample_yaml_pois.yml -o file.gpx
 
+
+Render track with OpenStreetMap
+---------------------
+
+You can "convert" your tracks to images using this command.
+
+How to use it
+-------------
+
+1. Please check if you have installed RMagick gem.
+
+2. Run command.
+
+  gpx2png -g <input GPX file> -s <image size, format: WIDTHxHEIGHT> -o <output PPNG file>
+
+  Example:
+
+  gpx2png -g spec/fixtures/sample.gpx -s 800x600 -o map.png
+
+3. You can specify zoom.
+
+  gpx2png -g <input GPX file> -z <zoom, best results if between 9 and 15, max 18> -o <output PPNG file>
+
+  Example:
+
+  gpx2png -g spec/fixtures/sample.gpx -z 11 -o map.png
+
+
 Contributing to gpx2xif
 -------------------------------
 

data/Rakefile

@@ -22,7 +22,7 @@ Jeweler::Tasks.new do |gem|
   gem.email = "bobikx@poczta.fm"
   gem.authors = ["Aleksander Kwiatkowski"]
   # dependencies defined in Gemfile
-  gem.executables = ['geotag_all_images', 'geotag_simulate']
+  gem.executables = ['geotag_all_images', 'geotag_simulate', 'generate_garmin_waypoints', 'gpx2png']
 
   gem.files = FileList[
     "[A-Z]*", "{bin,generators,lib,test}/**/*"

data/VERSION

@@ -1 +1 @@
-0.0.4
+0.1.0

data/bin/generate_garmin_waypoints

@@ -22,7 +22,7 @@ OptionParser.new do |opts|
   end
 end.parse!
 
-g = GarminUtils::WaypointListGenerator.new
+g = GpxUtils::WaypointsExporter.new
 if options[:yaml]
   g.add_yaml_file(options[:yaml])
 end

data/bin/geotag_all_images

@@ -1,7 +1,7 @@
 #!/usr/bin/env ruby
 
 require 'rubygems'
-require 'gpx2exif'
+require 'geotagger/geotagger'
 
 puts "Are you sure? It is evil script which probably eat photos of your dog and family. Uppercase 'yes' and enter if you want to continue."
 str = gets
@@ -12,7 +12,7 @@ time_offset = time_offset.to_i
 
 exit(0) unless str.strip == 'YES'
 
-g = Gpx2exif::GeoManager.new
+g = Geotagger::Geotagger.new(verbose: true)
 g.add_all_files(time_offset)
 g.match_up
 g.save!

data/bin/geotag_simulate

@@ -1,7 +1,7 @@
 #!/usr/bin/env ruby
 
 require 'rubygems'
-require 'gpx2exif'
+require 'geotagger/geotagger'
 
 #puts "Are you sure? It is evil script which probably eat photos of your dog and family. Uppercase 'yes' and enter if you want to continue."
 #str = gets
@@ -11,7 +11,7 @@ puts "Do you want to add offset to image time? Default is 0 seconds."
 time_offset = gets
 time_offset = time_offset.to_i
 
-g = Gpx2exif::GeoManager.new
+g = Geotagger::Geotagger.new(verbose: true)
 g.add_all_files(time_offset)
 g.match_up
 g.simulate

data/bin/gpx2png

@@ -0,0 +1,57 @@
+#!/usr/bin/env ruby
+
+require 'rubygems'
+require 'gpx2exif'
+require 'gpx2png/osm'
+require 'optparse'
+
+options = { }
+OptionParser.new do |opts|
+  opts.banner = "Usage: gpx2png [options]"
+
+  opts.on("-g", "--gpx FILE", "Input GPX file") do |v|
+    options[:gpx] = v
+  end
+  opts.on("-z", "--zoom ZOOM", "Set zoom") do |v|
+    options[:zoom] = v
+  end
+  opts.on("-s", "--size WIDTHxHEIGHT", "Set output image size (better result)") do |v|
+    options[:size] = v
+  end
+  opts.on("-o", "--output FILE", "Output image file") do |v|
+    options[:output_file] = v
+  end
+end.parse!
+
+unless options[:gpx]
+  puts "Input GPX file needed"
+  @fail = true
+end
+unless options[:output_file]
+  puts "Output image file needed"
+  @fail = true
+end
+if options[:zoom].nil? and options[:size].nil?
+  puts "Zoom or image size needed"
+  @fail = true
+end
+
+unless @fail
+  g = GpxUtils::TrackImporter.new
+  g.add_file(options[:gpx])
+
+  e = Gpx2png::Osm.new
+  e.coords = g.coords
+
+  if options[:size]
+    # constant size
+    options[:size] =~ /(\d+)x(\d+)/
+    e.fixed_size($1.to_i, $1.to_i)
+  else
+    # constant zoom
+    e.zoom = options[:zoom].to_i
+    e.renderer_options = {crop_enabled: true}
+  end
+
+  e.save(options[:output_file])
+end

data/lib/geotagger.rb

@@ -0,0 +1,9 @@
+$:.unshift(File.dirname(__FILE__))
+
+require 'gpx_utils/waypoints_exporter'
+require 'gpx_utils/waypoints_importer'
+
+require 'gpx_utils/track_importer'
+
+module Geotagger
+end

data/lib/{gpx2exif → geotagger}/exif_editor.rb

@@ -3,12 +3,13 @@ require 'mini_exiftool'
 
 $:.unshift(File.dirname(__FILE__))
 
-module Gpx2exif
+module Geotagger
   class ExifEditor
 
-    def initialize
+    def initialize(options = {})
       @images = Array.new
       @global_time_offset = 0
+      @verbose = options[:verbose]
     end
 
     attr_reader :images
@@ -20,7 +21,7 @@ module Gpx2exif
         :time => get_photo_time(path) + time_offset + @global_time_offset
       }
       @images << i
-      puts "Added file #{path}, time #{i[:time]}"
+      puts "Added file #{path}, time #{i[:time]}" if @verbose
     end
 
     def get_photo_time(path)
@@ -54,14 +55,13 @@ module Gpx2exif
       photo.save
 
       photo2 = MiniExiftool.new path
-      puts " - coord saved lat #{photo2['GPSLatitude']} lon #{photo2['GPSLongitude']}"
+      puts " - coord saved lat #{photo2['GPSLatitude']} lon #{photo2['GPSLongitude']}" if @verbose
 
       # exiftool -GPSMapDatum="WGS-84" -gps:GPSLatitude="34,57,57"
       # -gps:GPSLatitudeRef="N" -gps:GPSLongitude="83,17,59" -gps:GPSLongitudeRef="W"
       # -gps:GPSAltitudeRef="0" -GPSAltitude=1426 -gps:GPSMeasureMode=2 -City="RabunBald"
       # -State="North Carolina" -Country="USA" ~/Desktop/RabunBaldSummit_NC.jpg
-
     end
 
   end
-end
+end

data/lib/{gpx2exif/geo_manager.rb → geotagger/geotagger.rb}

@@ -1,15 +1,21 @@
 require 'rubygems'
+require 'gpx_utils'
+require 'geotagger/exif_editor'
+require 'geotagger/track_importer'
 
 $:.unshift(File.dirname(__FILE__))
 
-module Gpx2exif
-  class GeoManager
+module Geotagger
+  class Geotagger
 
-    def initialize
+    def initialize(options = {})
+      @verbose = options[:verbose]
       @ee = ExifEditor.new
-      @gp = GpxParser.new
+      @ti = TrackImporter.new
+      @ti.verbose = @verbose
     end
 
+    # Add all GPX and images with
     def add_all_files(time_offset = 0)
       # add all GPX
       Dir.glob("**/*.GPX", File::FNM_CASEFOLD).each do |f|
@@ -26,7 +32,7 @@ module Gpx2exif
     end
 
     def add_gpx_file(path)
-      @gp.add_file(path)
+      @ti.add_file(path)
     end
 
     def add_image(path, time_offset = 0)
@@ -35,10 +41,10 @@ module Gpx2exif
 
     def match_up
       @ee.images.each do |i|
-        puts "* searching for #{i[:path]}, time #{i[:time]}"
-        i[:coord] = @gp.find_by_time(i[:time])
+        puts "* searching for #{i[:path]}, time #{i[:time]}" if @verbose
+        i[:coord] = @ti.find_by_time(i[:time])
         if i[:coord].nil?
-          puts " - not found"
+          puts " - not found" if @verbose
         end
       end
 
@@ -48,7 +54,7 @@ module Gpx2exif
     def save!
       @ee.images.each do |i|
         if not i[:coord].nil?
-          puts "! saving for #{i[:path]}"
+          puts "! saving for #{i[:path]}" if @verbose
           @ee.set_photo_coords_internal(i)
         end
 
@@ -57,7 +63,7 @@ module Gpx2exif
 
     def simulate
       to_process = @ee.images.select { |i| not i[:coord].nil? }
-      puts "Result: to update #{to_process.size} from #{@ee.images.size}"
+      puts "Result: to update #{to_process.size} from #{@ee.images.size}" if @verbose
     end
 
   end

data/lib/geotagger/track_importer.rb

@@ -0,0 +1,33 @@
+require 'rubygems'
+require 'nokogiri'
+
+$:.unshift(File.dirname(__FILE__))
+
+# Simple parsing GPX file
+module Geotagger
+  class TrackImporter < GpxUtils::TrackImporter
+
+    THRESHOLD = 5*60
+
+    attr_accessor :verbose
+
+    # Only import valid coords
+    def self.coord_valid?(lat, lon, elevation, time)
+      return true if lat and lon and time
+      return false
+    end
+
+    def find_by_time(time)
+      selected_coords = @coords.select { |c| (c[:time].localtime - time.localtime).abs < THRESHOLD }
+      selected_coords = selected_coords.sort { |a, b| (a[:time].localtime - time.localtime).abs <=> (b[:time].localtime - time.localtime).abs }
+      puts " - found #{selected_coords.size} coords within #{THRESHOLD}s from image time" if @verbose
+      if selected_coords.size > 0
+        puts " - best is #{selected_coords.first[:time].localtime}, time offset #{selected_coords.first[:time].localtime - time.localtime}" if @verbose
+        puts " - lat #{selected_coords.first[:lat]} lon #{selected_coords.first[:lon]}" if @verbose
+      end
+
+      return selected_coords.first
+    end
+
+  end
+end

data/lib/gpx2exif.rb

@@ -1,10 +1,6 @@
 $:.unshift(File.dirname(__FILE__))
 
-require 'garmin_utils'
+require 'gpx_utils'
 
-require 'gpx2exif/gpx_parser'
-require 'gpx2exif/exif_editor'
-require 'gpx2exif/geo_manager'
-
-module Gpx2exif
-end
+#module Gpx2exif
+#end

data/lib/gpx2png/base.rb

@@ -0,0 +1,30 @@
+require 'rubygems'
+
+$:.unshift(File.dirname(__FILE__))
+
+module Gpx2png
+  class Base
+
+    def initialize
+      @coords = Array.new
+      @zoom = 9
+      @verbose = true
+    end
+
+    def add(lat, lon)
+      @coords << { lat: lat, lon: lon }
+    end
+
+    attr_accessor :zoom, :color, :coords
+
+    # Some math stuff
+    def self.rad2deg(rad)
+      return rad * 180.0 / Math::PI
+    end
+
+    def self.deg2rad(deg)
+      return deg * Math::PI / 180.0
+    end
+
+  end
+end

data/lib/gpx2png/gpx2png.rb

@@ -0,0 +1,210 @@
+require 'rubygems'
+require 'chunky_png'
+require 'net/http'
+require "uri"
+
+
+$:.unshift(File.dirname(__FILE__))
+
+module Gpx2png
+  class Gpx2png
+
+    TILE_WIDTH = 256
+    TILE_HEIGHT = 256
+
+    def initialize
+      @coords = Array.new
+      @zoom = 9
+      @color = ChunkyPNG::Color.from_hex('#FF0000')
+    end
+
+    def add(lat, lon)
+      @coords << { lat: lat, lon: lon }
+    end
+
+    attr_accessor :zoom, :color, :coords
+
+    def dev
+      zoom = 15
+      @coords.collect { |c|
+        {
+          url: self.class.url(zoom, [c[:lat], c[:lon]]),
+          tile: self.class.convert(zoom, [c[:lat], c[:lon]]),
+          return: self.class.reverse_convert(zoom,
+                                             self.class.convert(zoom, [c[:lat], c[:lon]])
+          ),
+          point: self.class.point_on_image(zoom, [c[:lat], c[:lon]])
+        }
+      }
+    end
+
+    def to_png(filename)
+      download_and_join_tiles
+      @full_image.save(filename)
+      filename
+    end
+
+    # http://wiki.openstreetmap.org/wiki/Slippy_map_tilenames#X_and_Y
+    def self.convert(zoom, coord)
+      lat_deg, lon_deg = coord
+      lat_rad = deg2rad(lat_deg)
+      x = (((lon_deg + 180) / 360) * (2 ** zoom)).floor
+      y = ((1 - Math.log(Math.tan(lat_rad) + 1 / Math.cos(lat_rad)) / Math::PI) /2 * (2 ** zoom)).floor
+
+      return [x, y]
+    end
+
+    def self.url_convert(zoom, coord, server = 'b.')
+      x, y = convert(zoom, coord)
+      url(zoom, [x, y], server)
+    end
+
+    def self.url(zoom, coord, server = 'b.')
+      x, y = coord
+      url = "http://#{server}tile.openstreetmap.org\/#{zoom}\/#{x}\/#{y}.png"
+      return url
+    end
+
+    # top-left corner
+    def self.reverse_convert(zoom, coord)
+      x, y = coord
+      n = 2 ** zoom
+      lon_deg = x.to_f / n.to_f * 360.0 - 180.0
+      lat_deg = rad2deg(Math.atan(Math.sinh(Math::PI * (1.to_f - 2.to_f * y.to_f / n.to_f))))
+      return [lat_deg, lon_deg]
+    end
+
+    # return where you should put point on tile
+    def self.point_on_image(zoom, geo_coord)
+      osm_tile_coord = convert(zoom, geo_coord)
+      top_left_corner = reverse_convert(zoom, osm_tile_coord)
+      bottom_right_corner = reverse_convert(zoom, [
+        osm_tile_coord[0] + 1, osm_tile_coord[1] + 1
+      ])
+
+      # some line y = ax + b math
+
+      x_geo = geo_coord[1]
+      # offset
+      x_offset = x_geo - top_left_corner[1]
+      # scale
+      x_distance = (bottom_right_corner[1] - top_left_corner[1])
+      x = (TILE_WIDTH.to_f * (x_offset / x_distance)).round
+
+      y_geo = geo_coord[0]
+      # offset
+      y_offset = y_geo - top_left_corner[0]
+      # scale
+      y_distance = (bottom_right_corner[0] - top_left_corner[0])
+      y = (TILE_HEIGHT.to_f * (y_offset / y_distance)).round
+
+      return { osm_title_coord: osm_tile_coord, pixel_offset: [x, y] }
+    end
+
+
+    attr_reader :lat_min, :lat_max, :lon_min, :lon_max
+    attr_reader :tile_x_distance, :tile_y_distance
+
+    def download_and_join_tiles
+      @lat_min = @coords.collect { |c| c[:lat] }.min
+      @lat_max = @coords.collect { |c| c[:lat] }.max
+      @lon_min = @coords.collect { |c| c[:lon] }.min
+      @lon_max = @coords.collect { |c| c[:lon] }.max
+
+      @border_tiles = [
+        self.class.convert(@zoom, [@lat_min, @lon_min]),
+        self.class.convert(@zoom, [@lat_max, @lon_max])
+      ]
+
+      @tile_x_range = (@border_tiles[0][0])..(@border_tiles[1][0])
+      @tile_y_range = (@border_tiles[1][1])..(@border_tiles[0][1])
+
+      # new image
+      @full_image_x = (1 + @tile_x_range.max - @tile_x_range.min) * TILE_WIDTH
+      @full_image_y = (1 + @tile_y_range.max - @tile_y_range.min) * TILE_HEIGHT
+      puts "Output image dimension #{@full_image_x}x#{@full_image_y}"
+      @full_image = ChunkyPNG::Image.new(
+        @full_image_x,
+        @full_image_y,
+        ChunkyPNG::Color::WHITE
+      )
+
+      # {:x, :y, :blob}
+      @images = Array.new
+
+      @tile_x_range.each do |x|
+        @tile_y_range.each do |y|
+          url = self.class.url(@zoom, [x,y])
+
+          # blob time
+          uri = URI.parse(url)
+          response = Net::HTTP.get_response(uri)
+          blob = response.body
+          image = ChunkyPNG::Image.from_blob(blob)
+
+          @images << {
+            url: url,
+            image: image,
+            x: x,
+            y: y
+          }
+
+          # compose image
+          x_offset = (x - @tile_x_range.min) * TILE_WIDTH
+          y_offset = (y - @tile_y_range.min) * TILE_HEIGHT
+          @full_image.compose!(
+            image,
+            x_offset,
+            y_offset
+          )
+          
+          puts "processed #{x - @tile_x_range.min}x#{y - @tile_y_range.min} (max #{@tile_x_range.max - @tile_x_range.min}x#{@tile_y_range.max - @tile_y_range.min})"
+        end
+      end
+
+      # sweet, image is joined
+
+      # add some coords to the map
+      (1...@coords.size).each do |i|
+        lat_from = @coords[i-1][:lat]
+        lon_from = @coords[i-1][:lon]
+
+        lat_to = @coords[i][:lat]
+        lon_to = @coords[i][:lon]
+
+        point_from = self.class.point_on_image(@zoom, [lat_from, lon_from])
+        point_to = self.class.point_on_image(@zoom, [lat_to, lon_to])
+        # { osm_title_coord: osm_tile_coord, pixel_offset: [x, y] }
+
+        # first point
+        bitmap_xa = (point_from[:osm_title_coord][0] - @tile_x_range.min) * TILE_WIDTH + point_from[:pixel_offset][0]
+        bitmap_ya = (point_from[:osm_title_coord][1] - @tile_y_range.min) * TILE_HEIGHT + point_from[:pixel_offset][1]
+        bitmap_xb = (point_to[:osm_title_coord][0] - @tile_x_range.min) * TILE_WIDTH + point_to[:pixel_offset][0]
+        bitmap_yb = (point_to[:osm_title_coord][1] - @tile_y_range.min) * TILE_HEIGHT + point_to[:pixel_offset][1]
+
+        @full_image.line(
+          bitmap_xa, bitmap_ya,
+          bitmap_xb, bitmap_yb,
+          @color
+        )
+
+      end
+    end
+
+    def expand_map
+      # TODO expand min and max ranges
+    end
+
+
+    # Some math stuff
+    def self.rad2deg(rad)
+      return rad * 180.0 / Math::PI
+    end
+
+    def self.deg2rad(deg)
+      return deg * Math::PI / 180.0
+    end
+
+
+  end
+end