gps_pvt 0.4.1 → 0.5.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: ffde7da5c82fa8eecc03dc1fb71c630ff95d8c1b016aefe4da64b63f8ab9e928
-  data.tar.gz: 8a56df6b69e8f91cb888aa1135846b6db75a053bbd833e50c62bf0f4937a7157
+  metadata.gz: 83dff6013480dd4df4faacde170a901d712a605093fa006d0d4542c58f43ac45
+  data.tar.gz: 33728f6f9c4969f04ed528831e14977f3f90846b0f740ce0644f06de18ac7322
 SHA512:
-  metadata.gz: 2889b55a9f925b9952ac1f29900bce98d8e0ab000418d7cecdd0405a2d81b84e70659ef4da271bb924231b7dcd80d8a550937478fbdb9f63ec0b38c0ea6d3061
-  data.tar.gz: '08d73a274ee81613d1d44129cfbc1a9732cc886ef5ee48d3eb70113e7d700244ad5ea977abae2d3294d86ac142380c05ba7ccf6b266327f3d001f80be584b0dc'
+  metadata.gz: fd13f81c590a85ffeed54fe62dfb45858373cf687b08e5dc37bf45805f54398c50ad71f03e512a28015d214239e189c0581317d703cb4af99f4f9d258d162d12
+  data.tar.gz: 64df6a70a7b02f94e6f94c878ada85bc293871e5478698a5fdc285f77418fd2d27ecf9aece5b42ca197d6a3dc79ab9294f81b06411b83ca008f3428011245b83

data/README.md

@@ -32,7 +32,7 @@ An attached executable is useful. After installation, type
 
     $ gps_pvt RINEX_or_UBX_file(s)
 
-The format of RINEX_or_UBX_file is automatically determined with its extention, such as .ubx will be treated as UBX format. If you want to specify the file format, instead of RINEX_or_UBX_file(s), use the following arguments:
+The format of RINEX_or_UBX_file is automatically determined with its extension, such as .ubx will be treated as UBX format. A gz compressed file can be specified directly, and URI such as https://... is also acceptable since version 0.5.0. If you want to specify the file format, instead of RINEX_or_UBX_file(s), use the following arguments:
 
     --rinex_nav=filename
     --rinex_obs=filename
@@ -50,6 +50,7 @@ Additionally, the following command options *--key=value* are available.
 | elevation_mask_deg | numeric | satellite elevation mask specified in degrees. *ex) --elevation_mask_deg=10* | v0.3.0 |
 | start_time | time string | start time to perform solution. GPS, UTC and other formats are supported. *ex1) --start_time=1234:5678* represents 5678 seconds in 1234 GPS week, *ex2) --start_time="2000-01-01 00:00:00 UTC"* is in UTC format. | v0.3.3 |
 | end_time | time string | end time to perform solution. Its format is the same as start_time. | v0.3.3 |
+| online_ephemeris | | automatically load ephemeris published online based on observation | v0.5.0 |
 
 ### For developer
 

data/exe/gps_pvt

@@ -1,16 +1,19 @@
 #!/usr/bin/env ruby
 
 require 'gps_pvt'
+require 'uri'
 
 # runnable quick example to solve PVT by using RINEX NAV/OBS or u-blox ubx
 
 $stderr.puts <<__STRING__
 Usage: #{__FILE__} GPS_file1 GPS_file2 ...
 As GPS_file, rinex_nav(*.YYn, *.YYh, *.YYq, *.YYg), rinex_obs(*.YYo), and ubx(*.ubx) format are currently supported.
+(YY = last two digit of year)
 File format is automatically determined based on its extention described in above parentheses.
 If you want to specify its format manually, --rinex_(nav|obs)=file_name or --ubx=file_name are available.
 Supported RINEX versions are 2 and 3.
-Note: YY = last two digit of year.
+A file having additional ".gz" extension is recognized as a file compressed by zlib.
+Major URL such as http(s)://... is acceptable as an input file name. 
 __STRING__
 
 options = []
@@ -54,6 +57,9 @@ options.reject!{|opt|
     end
     misc_options[opt[0]] = t
     true
+  when :online_ephemeris
+    misc_options[opt[0]] = opt[1]
+    true
   else
     false
   end
@@ -61,19 +67,45 @@ options.reject!{|opt|
 
 # Check file existence and extension
 files.collect!{|fname, ftype|
-  raise "File not found: #{fname}" unless File::exist?(fname)
   ftype ||= case fname
-  when /\.\d{2}[nhqg]$/; :rinex_nav
-  when /\.\d{2}o$/; :rinex_obs
+  when /\.\d{2}[nhqg](?:\.gz)?$/; :rinex_nav
+  when /\.\d{2}o(?:\.gz)?$/; :rinex_obs
   when /\.ubx$/; :ubx
   else
     raise "Format cannot be guessed, use --(format, ex. rinex_nav)=#{fname}"
   end
+  fname = proc{
+    next fname if File::exist?(fname)
+    if uri = URI::parse(fname) and !uri.instance_of?(URI::Generic) then
+      next uri
+    end
+    raise "File not found: #{fname}"
+  }.call
   [fname, ftype]
 }
 
 rcv = GPS_PVT::Receiver::new(options)
 
+proc{|src|
+  next unless src
+  loader = proc{|t_meas|
+    utc = Time::utc(*t_meas.c_tm)
+    y, yday = [:year, :yday].collect{|f| utc.send(f)}
+    uri = URI::parse(
+        "ftp://gssc.esa.int/gnss/data/daily/" + 
+          "%04d/brdc/BRDC00IGS_R_%04d%03d0000_01D_MN.rnx.gz"%[y, y, yday])
+    rcv.parse_rinex_nav(uri)
+    uri
+  }
+  run_orig = rcv.method(:run)
+  eph_list = {}
+  rcv.define_singleton_method(:run){|meas, t_meas, *args|
+    w_d = [t_meas.week, (t_meas.seconds.to_i / 86400)]
+    eph_list[w_d] ||= loader.call(t_meas)
+    run_orig.call(meas, t_meas, *args)
+  }
+}.call(misc_options[:online_ephemeris])
+
 proc{
   run_orig = rcv.method(:run)
   t_start, t_end = [nil, nil]

data/lib/gps_pvt/receiver.rb

@@ -4,6 +4,7 @@ Receiver class to be an top level interface to a user
 =end
 
 require_relative 'GPS'
+require_relative 'util'
 
 module GPS_PVT
 class Receiver
@@ -547,22 +548,24 @@ class Receiver
     $stderr.puts ", found packets are %s"%[ubx_kind.inspect]
   end
   
-  def parse_rinex_nav(fname)
+  def parse_rinex_nav(src)
+    fname = Util::get_txt(src)
     items = [
       @solver.gps_space_node,
       @solver.sbas_space_node,
       @solver.glonass_space_node,
     ].inject(0){|res, sn|
       loaded_items = sn.send(:read, fname)
-      raise "Format error! (Not RINEX) #{fname}" if loaded_items < 0
+      raise "Format error! (Not RINEX) #{src}" if loaded_items < 0
       res + loaded_items
     }
-    $stderr.puts "Read RINEX NAV file (%s): %d items."%[fname, items]
+    $stderr.puts "Read RINEX NAV file (%s): %d items."%[src, items]
   end
   
-  def parse_rinex_obs(fname, &b)
+  def parse_rinex_obs(src, &b)
+    fname = Util::get_txt(src)
     after_run = b || proc{|pvt| puts pvt.to_s if pvt}
-    $stderr.print "Reading RINEX observation file (%s)"%[fname]
+    $stderr.print "Reading RINEX observation file (%s)"%[src]
     types = nil
     glonass_freq = nil
     count = 0

data/lib/gps_pvt/util.rb

@@ -0,0 +1,32 @@
+require 'open-uri'
+require 'tempfile'
+require 'uri'
+require 'zlib'
+
+module GPS_PVT
+module Util
+  class << self
+    def inflate(src)
+      Zlib::GzipReader.send(*(src.kind_of?(IO) ? [:new, src] : [:open, src]))
+    end
+    def get_txt(fname_or_uri)
+      is_uri = fname_or_uri.kind_of?(URI)
+      ((is_uri && (RUBY_VERSION >= "2.5.0")) ? URI : Kernel) \
+          .send(:open, fname_or_uri){|src|
+        is_gz = (src.content_type =~ /gzip/) if is_uri
+        is_gz ||= (fname_or_uri.to_s =~ /\.gz$/)
+        is_file = src.kind_of?(File) || src.kind_of?(Tempfile)
+
+        return src.path if ((!is_gz) and is_file)
+
+        Tempfile::open(File::basename($0, '.*')){|dst|
+          dst.binmode
+          dst.write((is_gz ? inflate(is_file ? src.path : src) : src).read)
+          dst.rewind
+          dst.path
+        }
+      }
+    end
+  end
+end
+end

data/lib/gps_pvt/version.rb

@@ -1,5 +1,5 @@
 # frozen_string_literal: true
 
 module GPS_PVT
-  VERSION = "0.4.1"
+  VERSION = "0.5.0"
 end

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: gps_pvt
 version: !ruby/object:Gem::Version
-  version: 0.4.1
+  version: 0.5.0
 platform: ruby
 authors:
 - fenrir(M.Naruoka)
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2022-03-30 00:00:00.000000000 Z
+date: 2022-04-11 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake
@@ -93,11 +93,12 @@ files:
 - ext/ninja-scan-light/tool/swig/spec/GPS_spec.rb
 - ext/ninja-scan-light/tool/swig/spec/SylphideMath_spec.rb
 - gps_pvt.gemspec
+- gps_pvt.rbs
 - lib/gps_pvt.rb
 - lib/gps_pvt/receiver.rb
 - lib/gps_pvt/ubx.rb
+- lib/gps_pvt/util.rb
 - lib/gps_pvt/version.rb
-- sig/gps_pvt.rbs
 homepage: https://github.com/fenrir-naru/gps_pvt
 licenses: []
 metadata: