WIDXline() [1/2]

std::string sff::WIDXline	(	const sff::WID2 &	wid2,
		const bool &	debug
	)

write WID2 information in extended format

write WID information in extended format.

This function breaks the SFF format definition. It writes the WID2 class data to a line with identifier WIDX rather than WID2. The function is implemented within this library to provide transparent reading of WID-data. Writing will be used outside of the library only.

Reasons to implement this function are:

1. The GSE WID2 format does not use microseconds for the time of the first sample. The information present in sff::WID2 however is precise to the level of microseconds and is also used in this way when processing e.g. MiniSEED data. It is thus desirable to read and write this information with the full available content.
2. The GSE WID2 format restricts the sampling to be defined in terms of a sampling rate. This is inappropriate to extremely long period data with e.g. 1 Minute sampling interval (i.e. 60s). Round-off errors will lead to time residuals in these cases of more than a second within one day.

See also

WIDXline(const std::string& line)

Definition at line 115 of file widXio.cc.

References sff::WID2::auxid, sff::WID2::calib, sff::WID2::calper, sff::WID2::channel, sff::WID2::date, sff::WID2::dt, sff::WID2::hang, sff::WID2::instype, sff::WID2::nsamples, sff::helper::nsignificantdigits(), sff::WID2::station, sff::WID2::vang, and WIDXID.

Referenced by sff::WID2::read(), and test_extended_functions().

  {
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): entered function" <<
                                      std::endl; 
    }
    std::ostringstream os;
    os << sff::WIDXID;
    os << " " << wid2.date.hierarchicalstring();
    os << " ";
    os.width(5);
    os.setf(std::ios_base::left, std::ios_base::adjustfield);
    os << wid2.station.substr(0,5);
    os << " ";
    os.width(3);
    os.setf(std::ios_base::left, std::ios_base::adjustfield);
    os << wid2.channel.substr(0,3);
    os << " ";
    os.width(4);
    os.setf(std::ios_base::left, std::ios_base::adjustfield);
    os << wid2.auxid.substr(0,4);
    os << " ";
    os.width(6);
    os.setf(std::ios_base::left, std::ios_base::adjustfield);
    os << wid2.instype.substr(0,6);
    os << " " << wid2.nsamples;
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): partial result:" <<
                                      os.str() << std::endl; 
    }
    double rate=1./wid2.dt;
    double value=wid2.dt;
    if (sff::helper::nsignificantdigits(rate) < 4)
    {
      os << " r ";
      value=rate;
    }
    else
    {
      os << " i ";
    }
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): partial result:" <<
                                      os.str() << std::endl; 
    }
    os.setf(std::ios_base::scientific, std::ios_base::floatfield);
    int ndigits=sff::helper::nsignificantdigits(value, debug);
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): ndigits:" <<
                                      ndigits << std::endl; 
    }
    os.precision(ndigits);
    os.width(ndigits+6);
    os << value;
    os << " ";
    os.precision(sff::helper::nsignificantdigits(wid2.calib, debug));
    os << wid2.calib;
    os << " ";
    os.precision(sff::helper::nsignificantdigits(wid2.calper, debug));
    os << wid2.calper;
    os << " ";
    os.precision(sff::helper::nsignificantdigits(wid2.hang, debug));
    os << wid2.hang;
    os << " ";
    os.precision(sff::helper::nsignificantdigits(wid2.vang, debug));
    os << wid2.vang;
    if (debug) 
    { 
      std::cerr << "DEBUG (WIDXline): final result:" <<
                                      os.str() << std::endl; 
    }
    return(os.str());
  } // std::string WIDXline(const sff::WID2& wid2)

Here is the call graph for this function:

Here is the caller graph for this function: