gps_pvt 0.3.3 → 0.5.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: b91245b6a851b08ef5ec4b23e6be27cdd1150eaf6f7036249bf971c105bd920d
-  data.tar.gz: 75e5a367d8e1cc0d0e8cc48dae63b3b9b6a72fc86eeeea76bc6d1562ba171f02
+  metadata.gz: 83dff6013480dd4df4faacde170a901d712a605093fa006d0d4542c58f43ac45
+  data.tar.gz: 33728f6f9c4969f04ed528831e14977f3f90846b0f740ce0644f06de18ac7322
 SHA512:
-  metadata.gz: 3763faded2f4a5ace56f3a5cfde25311792ea06bda6ce5eb987fa10e1add08d9eb505bdb697fefb84b95a212072a824c3c134d298f4adea13d4ea31edfbeb831
-  data.tar.gz: 76b55b2da8d816a4bd2be024191ac0b706e40624035f979411f9192d78e030f2ee10ab4ce186160e8de2006b175728db593b11c25020f734ed9ac22544f0d2fb
+  metadata.gz: fd13f81c590a85ffeed54fe62dfb45858373cf687b08e5dc37bf45805f54398c50ad71f03e512a28015d214239e189c0581317d703cb4af99f4f9d258d162d12
+  data.tar.gz: 64df6a70a7b02f94e6f94c878ada85bc293871e5478698a5fdc285f77418fd2d27ecf9aece5b42ca197d6a3dc79ab9294f81b06411b83ca008f3428011245b83

data/README.md

@@ -32,15 +32,15 @@ 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
     --ubx=filename
 
-From version 0.2.0, SBAS and QZSS are supported in addition to GPS. QZSS ranging is activated in default, however, SBAS is just utilized for ionospheric correction. If you want to activate SBAS ranging, "--with=(SBAS PRN number, ex. 137)" option is used with gps_pvt executable like
+Since version 0.2.0, SBAS and QZSS are supported in addition to GPS. Since version 0.4.0, GLONASS is also available. QZSS ranging is activated in default, however, SBAS is just utilized for ionospheric correction. GLONASS is also turned off by default. If you want to activate SBAS or GLONASS ranging, "--with=(system or PRN)" options are used with gps_pvt executable like
 
-    $ gps_pvt --with=137 RINEX_or_UBX_file(s)
+    $ gps_pvt --with=137 --with=GLONASS RINEX_or_UBX_file(s)
 
 Additionally, the following command options *--key=value* are available.
 
@@ -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/Rakefile

@@ -29,9 +29,11 @@ namespace :git do
           # same as #{repo}/.git/info/sparse-checkout
           "sparse-checkout set" + (<<-__SPARSE_PATTERNS__).lines.collect{|str| str.chomp.gsub(/^ */, ' ')}.join,
             /tool/param/
+            /tool/algorithm/integral.h
             /tool/navigation/GPS*
             /tool/navigation/SBAS*
             /tool/navigation/QZSS*
+            /tool/navigation/GLONASS*
             /tool/navigation/coordinate.h
             /tool/navigation/EGM.h
             /tool/navigation/MagneticField.h

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), rinex_obs(*.YYo), and ubx(*.ubx) format are currently supported.
+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}[nhq]$/; :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]