gps_pvt 0.8.2 → 0.8.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: b87c77fbb3844eb0fd835ef079b8a379708a6fe24e3de4ceeb6934206573f489
-  data.tar.gz: a280c1aa82bed41319e13069907940a062a073a0b01c95eefa3d3d8c1b32c517
+  metadata.gz: e2f23806309c4244499a91877edba1a5dfacb1a9108596527dffad650c9cd9ee
+  data.tar.gz: dc33e31c739a79dc7e7c28ad0d6f404b65788499dfa2f53f9b93edae7b589e7e
 SHA512:
-  metadata.gz: 761d7e99868e04d93adb48fac0345ff71467068657c544efe857ed4f6206ed8d11698871ad7dc73c232f4bd51094cf13c6752d96fc17307caccb4a2eab8dbe03
-  data.tar.gz: fd70c1d59dbf3885e7a2286726a9c16e0cfbec3097a105bbffde59e3a65ee22eb09f1d0b4032b371faf2d1695b074b0752f37763ef6b5716dd4214fed4566755
+  metadata.gz: 2a3ec484712ccbd06b99a08a672ef20980d68e5f1670b07b093f1d6e4d0480dce974cbe8a181ecb0e8d6074cc60c1c122778113e962625f57b41a7a85b61fc4a
+  data.tar.gz: 5368619df05a35b6295702e3ad64958aab84e92652b01747f31b5bb98974243ac56a5b4fc25bde1075109bd4d586eac14bf64e1c2ac6b5a31d9bc990ece92fd7

data/README.md

@@ -27,7 +27,7 @@ For Windows users, this gem requires Devkit because of native compilation.
 
 ## Usage
 
-### For user who just generate PVT solution
+### For user who just wants to generate PVT solution
 An attached executable is useful. After installation, type
 
     $ gps_pvt file_or_URI(s)
@@ -36,12 +36,12 @@ The format of file is automatically determined with its extension, such as .ubx
 
 | specification | recoginized as |
 ----|----
-| --rinex_nav=file_or_URI | [RINEX](https://www.igs.org/wg/rinex/#documents-formats) navigation file |
-| --rinex_obs=file_or_URI | [RINEX](https://www.igs.org/wg/rinex/#documents-formats) observation file |
-| --ubx=file_or_URI | [U-blox](https://www.u-blox.com/) dedicated format |
+| <a name=opt_rinex_nav>--rinex_nav=file_or_URI</a> | [RINEX](https://www.igs.org/wg/rinex/#documents-formats) navigation file |
+| <a name=opt_rinex_obs>--rinex_obs=file_or_URI</a> | [RINEX](https://www.igs.org/wg/rinex/#documents-formats) observation file |
+| <a name=opt_ubx>--ubx=file_or_URI</a> | [U-blox](https://www.u-blox.com/) dedicated format |
 | --sp3=file_or_URI | [Standard Product 3 Orbit Format](https://files.igs.org/pub/data/format/sp3c.txt) (supported gps_pvt version >= 0.6.0) |
 | --antex=file_or_URI | [Antenna Exchange Format](https://igs.org/wg/antenna#files) (supported gps_pvt version >= 0.6.0) |
-| --rinex_clk=file_or_URI | [RINEX clock](https://files.igs.org/pub/data/format/rinex_clock304.txt) file  (supported gps_pvt version >= 0.7.0) |
+| --rinex_clk=file_or_URI | [RINEX clock](https://files.igs.org/pub/data/format/rinex_clock304.txt) file (supported gps_pvt version >= 0.7.0) |
 
 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
 
@@ -49,13 +49,13 @@ Since version 0.2.0, SBAS and QZSS are supported in addition to GPS. Since versi
 
 Additionally, the following command options *--key=value* are available.
 
-| key | value | comment | version |
+| key | value | comment | since |
 ----|----|----|----
 | base_station | 3 \* (numeric+coordinate) | base position used for relative ENU position calculation. XYZ, NEU formats are acceptable. *ex1) --base_station=0X,0Y,0Z*, *ex2) --base_station=12.34N,56.789E,0U* | v0.1.7 |
 | 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 | | based on observation, automatically load ephemeris which is previously broadcasted from satellite and currently published online | v0.5.0 |
+| <a name=opt_online_ephemeris>online_ephemeris</a> | URL string | based on observation, ephemeris which is previously broadcasted from satellite and currently published online will automatically be loaded. If value is not given, the default source "ftp://gssc.esa.int/gnss/data/daily/%Y/brdc/BRDC00IGS_R_%Y%j0000_01D_MN.rnx.gz" is used. The value string is converted with [strftime](https://docs.ruby-lang.org/en/master/strftime_formatting_rdoc.html) before actual use. | v0.8.1 |
 
 ### For developer
 
@@ -148,6 +148,21 @@ receiver.solver.correction = { # provide by using a Hash
 }
 ```
 
+## Additional utilities
+
+### [to_ubx](exe/to_ubx)
+
+Utility to convert observation into u-blox ubx format and dump standard input. After installation of gps_pvt, to type
+
+    $ to_ubx file_or_URI(s) (options) > out.ubx
+
+saves resultant into out.ubx by using redirection. The shared options with gps_pvt executable are [rinex_obs](#opt_rinex_obs), [rinex_nav](#opt_rinex_nav), [ubx](#opt_ubx), and [online_ephemeris](#opt_online_ephemeris). In addition, the following options are available.
+
+| key | value | comment | since |
+----|----|----|----
+| ubx_rawx |  | Change output packet types to UBX-RAWX from its default UBX-RAW. | v0.8.1 |
+| broadcast_data |  | In addition to observation, ephemeris is inserted by using UBX-SFRB packets. If ubx_rawx option is specified, UBX-SFRBX is used instead of UBX-SFRB. | v0.8.1 |
+
 ## Development
 
 After checking out the repo, run `bin/setup` to install dependencies. Then, run `rake` to build library and run the tests. You can also run `bin/console` for an interactive prompt that will allow you to experiment.

data/exe/gps_pvt

@@ -59,7 +59,7 @@ options.reject!{|opt|
     misc_options[opt[0]] = t
     true
   when :online_ephemeris
-    misc_options[opt[0]] = opt[1]
+    (misc_options[opt[0]] ||= []) << opt[1]
     true
   else
     false

data/exe/to_ubx

@@ -40,6 +40,7 @@ files.collect!{|fname, ftype|
   ftype ||= case fname
   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_obs)=#{fname}"
   end
@@ -55,9 +56,12 @@ files.collect!{|fname, ftype|
 
 options.reject!{|opt|
   case opt[0]
-  when :ubx_rawx, :online_ephemeris, :broadcast_data, :eph_interval
+  when :ubx_rawx, :broadcast_data, :eph_interval
     misc_options[opt[0]] = opt[1]
     true
+  when :online_ephemeris
+    (misc_options[opt[0]] ||= []) << opt[1]
+    true
   else
     false
   end
@@ -95,14 +99,37 @@ files.each{|fname, ftype|
 # other files
 files.each{|fname, ftype|
   case ftype
+  when :ubx; rcv.parse_ubx(fname){}
   when :rinex_obs; rcv.parse_rinex_obs(fname){}
   end
 }
 
 obs.sort!{|a, b| a[0] <=> b[0]} # Sort by measurement time
 
-# Add dummy before actual first observation in order to notify time to solver
-obs.unshift([obs[0][0] - 1, GPS_PVT::GPS::Measurement::new])
+# Time packet for solution of leap seconds
+gen_gpstime = proc{
+  tmpl = [0xB5, 0x62, 0x01, 0x20, 16, 0]
+  gpst = GPS_PVT::GPS::Time
+  t_next = gpst::new(*gpst::leap_second_events.find{|wn, sec, leap| leap == 1}[0..1])
+  proc{|t_meas, meas|
+    next nil if t_meas < t_next
+    t_next = t_meas + 60 # 1 min. interval
+    ubx = tmpl.clone
+    t_sec = t_meas.seconds
+    t_msec = (t_sec * 1E3).round
+    t_nsec = ((t_sec * 1E3 - t_msec) * 1E6).round
+    leap = rcv.leap_seconds(t_meas)
+    ubx += [
+      t_msec, # ITOW  ms  GPS Millisecond time of Week
+      t_nsec, # Frac  ns  Nanoseconds remainder of rounded ms above, range -500000 .. 500000
+      t_meas.week, # week  -  GPS week (GPS time)
+      leap || 0, # LeapS  s  Leap Seconds (GPS-UTC)
+      leap ? 0x07 : 0x03, # validity bit field (0x01=ToW, 0x02=WN, 0x04=UTC)
+      10000, # TAcc  ns  Time Accuracy Estimate
+    ].pack("Vl<vcCV").unpack("C*")
+    GPS_PVT::UBX::update(ubx + [0, 0]).pack("C*")
+  }
+}.call
 
 gen_sfrb, gen_sfrbx = proc{
   cache = {}
@@ -335,6 +362,7 @@ gen_rawx = proc{|t_meas, meas| # Convert to RXM-RAWX(0x15)
 }
 
 gen_list = []
+gen_list << gen_gpstime unless misc_options[:ubx_rawx]
 gen_list << (misc_options[:ubx_rawx] ? gen_sfrbx : gen_sfrb) if misc_options[:broadcast_data]
 gen_list << (misc_options[:ubx_rawx] ? gen_rawx : gen_raw)
 STDOUT.binmode

data/ext/gps_pvt/GPS/GPS_wrap.cxx

@@ -1866,56 +1866,55 @@ int SWIG_Ruby_arity( VALUE proc, int minimal )
 #define SWIGTYPE_p_GPS_SpaceNodeT_double_t__Ionospheric_UTC_Parameters swig_types[16]
 #define SWIGTYPE_p_GPS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris swig_types[17]
 #define SWIGTYPE_p_GPS_TimeT_double_t swig_types[18]
-#define SWIGTYPE_p_GPS_TimeT_double_t__leap_second_event_t swig_types[19]
-#define SWIGTYPE_p_GPS_TimeT_double_t__leap_year_prop_res_t swig_types[20]
-#define SWIGTYPE_p_GPS_User_PVTT_double_t swig_types[21]
-#define SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t swig_types[22]
-#define SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t swig_types[23]
-#define SWIGTYPE_p_MatrixViewBaseT_t swig_types[24]
-#define SWIGTYPE_p_MatrixViewFilterT_MatrixViewBaseT_t_t swig_types[25]
-#define SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t swig_types[26]
-#define SWIGTYPE_p_PushableData swig_types[27]
-#define SWIGTYPE_p_RINEX_ClockT_double_t swig_types[28]
-#define SWIGTYPE_p_RINEX_ObservationT_double_t swig_types[29]
-#define SWIGTYPE_p_SBAS_EphemerisT_double_t swig_types[30]
-#define SWIGTYPE_p_SBAS_SolverOptionsT_double_t swig_types[31]
-#define SWIGTYPE_p_SBAS_SpaceNodeT_double_t swig_types[32]
-#define SWIGTYPE_p_SBAS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris swig_types[33]
-#define SWIGTYPE_p_SP3T_double_t swig_types[34]
-#define SWIGTYPE_p_System_ENUT_double_WGS84_t swig_types[35]
-#define SWIGTYPE_p_System_LLHT_double_WGS84_t swig_types[36]
-#define SWIGTYPE_p_System_XYZT_double_WGS84_t swig_types[37]
-#define SWIGTYPE_p_available_satellites_t swig_types[38]
-#define SWIGTYPE_p_char swig_types[39]
-#define SWIGTYPE_p_double swig_types[40]
-#define SWIGTYPE_p_enu_t swig_types[41]
-#define SWIGTYPE_p_eph_t swig_types[42]
-#define SWIGTYPE_p_float_t swig_types[43]
-#define SWIGTYPE_p_gps_space_node_t swig_types[44]
-#define SWIGTYPE_p_gps_time_t swig_types[45]
-#define SWIGTYPE_p_int swig_types[46]
-#define SWIGTYPE_p_int_t swig_types[47]
-#define SWIGTYPE_p_llh_t swig_types[48]
-#define SWIGTYPE_p_p_GPS_Solver_RAIM_LSRT_double_GPS_Solver_Base_DebugT_double_GPS_Solver_BaseT_double_t_t_t__user_pvt_t__detection_t swig_types[49]
-#define SWIGTYPE_p_p_GPS_Solver_RAIM_LSRT_double_GPS_Solver_Base_DebugT_double_GPS_Solver_BaseT_double_t_t_t__user_pvt_t__exclusion_t swig_types[50]
-#define SWIGTYPE_p_p_double swig_types[51]
-#define SWIGTYPE_p_range_correction_list_t swig_types[52]
-#define SWIGTYPE_p_s16_t swig_types[53]
-#define SWIGTYPE_p_s32_t swig_types[54]
-#define SWIGTYPE_p_s8_t swig_types[55]
-#define SWIGTYPE_p_satellites_t swig_types[56]
-#define SWIGTYPE_p_self_t swig_types[57]
-#define SWIGTYPE_p_std__vectorT_std__pairT_int_SBAS_SpaceNodeT_double_t__Satellite_const_p_t_t swig_types[58]
-#define SWIGTYPE_p_super_t swig_types[59]
-#define SWIGTYPE_p_swig__GC_VALUE swig_types[60]
-#define SWIGTYPE_p_u16_t swig_types[61]
-#define SWIGTYPE_p_u32_t swig_types[62]
-#define SWIGTYPE_p_u8_t swig_types[63]
-#define SWIGTYPE_p_uint_t swig_types[64]
-#define SWIGTYPE_p_unsigned_int swig_types[65]
-#define SWIGTYPE_p_xyz_t swig_types[66]
-static swig_type_info *swig_types[68];
-static swig_module_info swig_module = {swig_types, 67, 0, 0, 0, 0};
+#define SWIGTYPE_p_GPS_TimeT_double_t__leap_year_prop_res_t swig_types[19]
+#define SWIGTYPE_p_GPS_User_PVTT_double_t swig_types[20]
+#define SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t swig_types[21]
+#define SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t swig_types[22]
+#define SWIGTYPE_p_MatrixViewBaseT_t swig_types[23]
+#define SWIGTYPE_p_MatrixViewFilterT_MatrixViewBaseT_t_t swig_types[24]
+#define SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t swig_types[25]
+#define SWIGTYPE_p_PushableData swig_types[26]
+#define SWIGTYPE_p_RINEX_ClockT_double_t swig_types[27]
+#define SWIGTYPE_p_RINEX_ObservationT_double_t swig_types[28]
+#define SWIGTYPE_p_SBAS_EphemerisT_double_t swig_types[29]
+#define SWIGTYPE_p_SBAS_SolverOptionsT_double_t swig_types[30]
+#define SWIGTYPE_p_SBAS_SpaceNodeT_double_t swig_types[31]
+#define SWIGTYPE_p_SBAS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris swig_types[32]
+#define SWIGTYPE_p_SP3T_double_t swig_types[33]
+#define SWIGTYPE_p_System_ENUT_double_WGS84_t swig_types[34]
+#define SWIGTYPE_p_System_LLHT_double_WGS84_t swig_types[35]
+#define SWIGTYPE_p_System_XYZT_double_WGS84_t swig_types[36]
+#define SWIGTYPE_p_available_satellites_t swig_types[37]
+#define SWIGTYPE_p_char swig_types[38]
+#define SWIGTYPE_p_double swig_types[39]
+#define SWIGTYPE_p_enu_t swig_types[40]
+#define SWIGTYPE_p_eph_t swig_types[41]
+#define SWIGTYPE_p_float_t swig_types[42]
+#define SWIGTYPE_p_gps_space_node_t swig_types[43]
+#define SWIGTYPE_p_gps_time_t swig_types[44]
+#define SWIGTYPE_p_int swig_types[45]
+#define SWIGTYPE_p_int_t swig_types[46]
+#define SWIGTYPE_p_llh_t swig_types[47]
+#define SWIGTYPE_p_p_GPS_Solver_RAIM_LSRT_double_GPS_Solver_Base_DebugT_double_GPS_Solver_BaseT_double_t_t_t__user_pvt_t__detection_t swig_types[48]
+#define SWIGTYPE_p_p_GPS_Solver_RAIM_LSRT_double_GPS_Solver_Base_DebugT_double_GPS_Solver_BaseT_double_t_t_t__user_pvt_t__exclusion_t swig_types[49]
+#define SWIGTYPE_p_p_double swig_types[50]
+#define SWIGTYPE_p_range_correction_list_t swig_types[51]
+#define SWIGTYPE_p_s16_t swig_types[52]
+#define SWIGTYPE_p_s32_t swig_types[53]
+#define SWIGTYPE_p_s8_t swig_types[54]
+#define SWIGTYPE_p_satellites_t swig_types[55]
+#define SWIGTYPE_p_self_t swig_types[56]
+#define SWIGTYPE_p_std__vectorT_std__pairT_int_SBAS_SpaceNodeT_double_t__Satellite_const_p_t_t swig_types[57]
+#define SWIGTYPE_p_super_t swig_types[58]
+#define SWIGTYPE_p_swig__GC_VALUE swig_types[59]
+#define SWIGTYPE_p_u16_t swig_types[60]
+#define SWIGTYPE_p_u32_t swig_types[61]
+#define SWIGTYPE_p_u8_t swig_types[62]
+#define SWIGTYPE_p_uint_t swig_types[63]
+#define SWIGTYPE_p_unsigned_int swig_types[64]
+#define SWIGTYPE_p_xyz_t swig_types[65]
+static swig_type_info *swig_types[67];
+static swig_module_info swig_module = {swig_types, 66, 0, 0, 0, 0};
 #define SWIG_TypeQuery(name) SWIG_TypeQueryModule(&swig_module, &swig_module, name)
 #define SWIG_MangledTypeQuery(name) SWIG_MangledTypeQueryModule(&swig_module, &swig_module, name)
 
@@ -6654,7 +6653,16 @@ SWIGINTERN VALUE
 _wrap_Time_leap_second_events_get(VALUE self) {
   VALUE _val;
   
-  _val = SWIG_NewPointerObj(SWIG_as_voidptr(GPS_Time< double >::leap_second_events), SWIGTYPE_p_GPS_TimeT_double_t__leap_second_event_t,  0 );
+  {
+    _val = rb_ary_new2(sizeof(GPS_Time< double >::leap_second_events) / sizeof(GPS_Time< double >::leap_second_events[0]));
+    for(std::size_t i(0); i < sizeof(GPS_Time< double >::leap_second_events) / sizeof(GPS_Time< double >::leap_second_events[0]); ++i){
+      rb_ary_push(_val,
+        rb_ary_new3(3,
+          SWIG_From_int  (GPS_Time< double >::leap_second_events[i].corrected.week),
+          swig::from(GPS_Time< double >::leap_second_events[i].corrected.seconds),
+          SWIG_From_int  (GPS_Time< double >::leap_second_events[i].leap_seconds)));
+    }
+  }
   return _val;
 }
 
@@ -27470,7 +27478,6 @@ static swig_type_info _swigt__p_GPS_SpaceNodeT_double_t = {"_p_GPS_SpaceNodeT_do
 static swig_type_info _swigt__p_GPS_SpaceNodeT_double_t__Ionospheric_UTC_Parameters = {"_p_GPS_SpaceNodeT_double_t__Ionospheric_UTC_Parameters", "GPS_SpaceNode< double >::Ionospheric_UTC_Parameters *", 0, 0, (void*)0, 0};
 static swig_type_info _swigt__p_GPS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris = {"_p_GPS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris", "GPS_SpaceNode< double >::SatelliteProperties::Ephemeris *", 0, 0, (void*)0, 0};
 static swig_type_info _swigt__p_GPS_TimeT_double_t = {"_p_GPS_TimeT_double_t", "GPS_SpaceNode< double >::gps_time_t *|SBAS_SpaceNode< double >::gps_time_t *|GPS_Time< GLONASS_SpaceNode< double >::float_t > *|GPS_Time< GPS_Time< double >::float_t > *|GPS_Time< double > *", 0, 0, (void*)0, 0};
-static swig_type_info _swigt__p_GPS_TimeT_double_t__leap_second_event_t = {"_p_GPS_TimeT_double_t__leap_second_event_t", "GPS_Time< double >::leap_second_event_t *", 0, 0, (void*)0, 0};
 static swig_type_info _swigt__p_GPS_TimeT_double_t__leap_year_prop_res_t = {"_p_GPS_TimeT_double_t__leap_year_prop_res_t", "GPS_Time< double >::leap_year_prop_res_t *", 0, 0, (void*)0, 0};
 static swig_type_info _swigt__p_GPS_User_PVTT_double_t = {"_p_GPS_User_PVTT_double_t", "GPS_User_PVT< double > *", 0, 0, (void*)0, 0};
 static swig_type_info _swigt__p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t = {"_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t", "Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *|Matrix< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > *", 0, 0, (void*)0, 0};
@@ -27539,7 +27546,6 @@ static swig_type_info *swig_type_initial[] = {
   &_swigt__p_GPS_SpaceNodeT_double_t__Ionospheric_UTC_Parameters,
   &_swigt__p_GPS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris,
   &_swigt__p_GPS_TimeT_double_t,
-  &_swigt__p_GPS_TimeT_double_t__leap_second_event_t,
   &_swigt__p_GPS_TimeT_double_t__leap_year_prop_res_t,
   &_swigt__p_GPS_User_PVTT_double_t,
   &_swigt__p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t,
@@ -27608,7 +27614,6 @@ static swig_cast_info _swigc__p_GPS_SpaceNodeT_double_t[] = {  {&_swigt__p_GPS_S
 static swig_cast_info _swigc__p_GPS_SpaceNodeT_double_t__Ionospheric_UTC_Parameters[] = {  {&_swigt__p_GPS_Ionospheric_UTC_ParametersT_double_t, _p_GPS_Ionospheric_UTC_ParametersT_double_tTo_p_GPS_SpaceNodeT_double_t__Ionospheric_UTC_Parameters, 0, 0},  {&_swigt__p_GPS_SpaceNodeT_double_t__Ionospheric_UTC_Parameters, 0, 0, 0},{0, 0, 0, 0}};
 static swig_cast_info _swigc__p_GPS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris[] = {  {&_swigt__p_GPS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris, 0, 0, 0},  {&_swigt__p_GPS_EphemerisT_double_t, _p_GPS_EphemerisT_double_tTo_p_GPS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris, 0, 0},{0, 0, 0, 0}};
 static swig_cast_info _swigc__p_GPS_TimeT_double_t[] = {  {&_swigt__p_GPS_TimeT_double_t, 0, 0, 0},{0, 0, 0, 0}};
-static swig_cast_info _swigc__p_GPS_TimeT_double_t__leap_second_event_t[] = {  {&_swigt__p_GPS_TimeT_double_t__leap_second_event_t, 0, 0, 0},{0, 0, 0, 0}};
 static swig_cast_info _swigc__p_GPS_TimeT_double_t__leap_year_prop_res_t[] = {  {&_swigt__p_GPS_TimeT_double_t__leap_year_prop_res_t, 0, 0, 0},{0, 0, 0, 0}};
 static swig_cast_info _swigc__p_GPS_User_PVTT_double_t[] = {  {&_swigt__p_GPS_User_PVTT_double_t, 0, 0, 0},{0, 0, 0, 0}};
 static swig_cast_info _swigc__p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t[] = {  {&_swigt__p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0, 0, 0},{0, 0, 0, 0}};
@@ -27677,7 +27682,6 @@ static swig_cast_info *swig_cast_initial[] = {
   _swigc__p_GPS_SpaceNodeT_double_t__Ionospheric_UTC_Parameters,
   _swigc__p_GPS_SpaceNodeT_double_t__SatelliteProperties__Ephemeris,
   _swigc__p_GPS_TimeT_double_t,
-  _swigc__p_GPS_TimeT_double_t__leap_second_event_t,
   _swigc__p_GPS_TimeT_double_t__leap_year_prop_res_t,
   _swigc__p_GPS_User_PVTT_double_t,
   _swigc__p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t,

data/ext/ninja-scan-light/tool/swig/GPS.i

@@ -158,6 +158,16 @@ static std::string inspect_str(const VALUE &v){
   %typemap(typecheck, precedence=SWIG_TYPECHECK_POINTER) const std::tm & {
     $1 = (TYPE($input) == T_ARRAY) ? 1 : 0;
   }
+  %typemap(varout,fragment="SWIG_Traits_frag(FloatT)") leap_second_event_t [] {
+    $result = rb_ary_new2(sizeof($1) / sizeof($1[0]));
+    for(std::size_t i(0); i < sizeof($1) / sizeof($1[0]); ++i){
+      rb_ary_push($result,
+          rb_ary_new3(3,
+            SWIG_From(int)($1[i].corrected.week),
+            swig::from($1[i].corrected.seconds),
+            SWIG_From(int)($1[i].leap_seconds)));
+    }
+  }
 #endif
   %ignore canonicalize();
   %ignore GPS_Time(const int &_week, const float_t &_seconds);

data/ext/ninja-scan-light/tool/swig/spec/GPS_spec.rb

@@ -759,6 +759,10 @@ __RINEX_CLK_TEXT__
       puts "Measurement time: #{t_meas.to_a} (a.k.a #{"%d/%d/%d %02d:%02d:%02d UTC"%[*t_meas.c_tm]})"
       expect(t_meas.c_tm).to eq([2015, 6, 15, 23, 53, 33])
       expect(GPS::Time::new(0, t_meas.serialize)).to eq(t_meas)
+      expect(t_meas.leap_seconds).to eq(sn.iono_utc.delta_t_LS)
+      expect(GPS::Time::leap_second_events.select{|wn, sec, leap|
+        t_meas >= GPS::Time::new(wn, sec)
+      }[0][2]).to eq(sn.iono_utc.delta_t_LS)
       
       sn.update_all_ephemeris(t_meas)
       

data/lib/gps_pvt/receiver/extension.rb

@@ -29,15 +29,24 @@ class Receiver
     return (eph.valid?(t) ? eph : nil)
   end
 
-  def attach_online_ephemeris(uri_template = nil)
-    if (!uri_template) || (uri_template =~ /^\s*$/) then
-      uri_template = "ftp://gssc.esa.int/gnss/data/daily/%Y/brdc/BRDC00IGS_R_%Y%j0000_01D_MN.rnx.gz"
-    end
+  def attach_online_ephemeris(uri_template = [nil])
+    uri_template = uri_template.collect{|v|
+      if (!v) || (v =~ /^\s*$/) then
+        "ftp://gssc.esa.int/gnss/data/daily/%Y/brdc/BRDC00IGS_R_%Y%j0000_01D_MN.rnx.gz"
+      else
+        v
+      end
+    }.uniq
     loader = proc{|t_meas|
       utc = Time::utc(*t_meas.c_tm)
-      uri = URI::parse(utc.strftime(uri_template))
-      self.parse_rinex_nav(uri)
-      uri
+      uri_template.each{|v|
+        uri = URI::parse(utc.strftime(v))
+        begin
+          self.parse_rinex_nav(uri)
+        rescue Net::FTPError, Net::HTTPExceptions => e
+          $stderr.puts "Skip to read due to %s (%s)"%[e.inspect.gsub(/[\r\n]/, ' '), uri]
+        end
+      }
     }
     run_orig = self.method(:run)
     eph_list = {}

data/lib/gps_pvt/version.rb

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

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: gps_pvt
 version: !ruby/object:Gem::Version
-  version: 0.8.2
+  version: 0.8.3
 platform: ruby
 authors:
 - fenrir(M.Naruoka)
 autorequire: 
 bindir: exe
 cert_chain: []
-date: 2022-11-20 00:00:00.000000000 Z
+date: 2022-11-30 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rubyserial