readfile()

File datrw::tsoft::readfile	(	std::istream &	is,
		const ReaderConfig &	rc
	)

read complete file

function to read a complete TSOFT file

Definition at line 48 of file tsoftsffcontainer.cc.

References datrw::tsoft::SFFchannelid::auxid, datrw::tsoft::SFFchannelid::channel, channelinfofree(), datrw::tsoft::Channeldata::chinfo(), datrw::tsoft::TSOFTchannelid::datatype, datrw::tsoft::Datatrace::date(), datrw::tsoft::TSOFTfile::dc(), datrw::tsoft::TSOFTfile::free(), datrw::tsoft::TSOFTchannelid::instrument, datrw::tsoft::SFFchannelid::instrument, datrw::tsoft::Datatrace::interval(), datrw::tsoft::TSOFTchannelid::location, datrw::tsoft::File::Mfree, datrw::tsoft::Datatrace::Mseries, datrw::tsoft::File::Mtraces, datrw::tsoft::TSOFTfile::nchannels(), datrw::tsoft::Channeldata::ntraces(), datrw::tsoft::Datatrace::series(), datrw::tsoft::SFFchannelid::station, datrw::tsoft::Channelinfo::thedatatype(), datrw::tsoft::Channelinfo::theinstrument(), datrw::tsoft::Channelinfo::thelocation(), datrw::tsoft::Channeldata::trace(), and translate().

Referenced by datrw::itsoftstream::read().

    {
      TSOFTfile infile(is, rc);
      File outfile;
      outfile.Mfree.append(infile.free());
      const Datacontainer& dc=infile.dc();

      // count traces
      int ntraces=0;
      for (int ich=0; ich<infile.nchannels(); ++ich)
      {
        const Channeldata& cd=dc[ich];
        /*
        std::cout << "ich " << ich << 
          " has " << cd.nsequences() << " sequences" << std::endl;
          */
        ntraces += cd.ntraces();
      }
      // std::cout << "found " << ntraces << " traces" << std::endl;

      // reserve space for n traces
      outfile.Mtraces=Ttraceseries(ntraces);

      // cycle trough all input traces
      aff::Iterator<Ttraceseries> I(outfile.Mtraces);
      for (int ich=0; ich<infile.nchannels(); ++ich)
      {
        const Channeldata& cd=dc[ich];
        const Channelinfo& ci=cd.chinfo();
        sff::FREE channelfree;
        channelfree.append(channelinfofree(ci));
        for (int itr=0; itr<cd.ntraces(); ++itr)
        {
          I->Mfree=channelfree;
          const Datatrace& ds=cd.trace(itr);

          // extract samples
          I->Mseries=ds.series();

          // extract header
          I->Mwid2.date=ds.date();
          I->Mwid2.dt=libtime::time2double(ds.interval());
          I->Mwid2.nsamples=I->Mseries.size();

          TSOFTchannelid tci;
          tci.location=ci.thelocation();
          tci.instrument=ci.theinstrument();
          tci.datatype=ci.thedatatype();
          /*
          std::cout << "translate: " 
            << tci.location << ":" 
            << tci.instrument << ":"
            << tci.datatype << ":" << std::endl;
            */
          SFFchannelid sci=translate(tci);

          I->Mwid2.station=sci.station;
          I->Mwid2.channel=sci.channel;
          I->Mwid2.instype=sci.instrument;
          I->Mwid2.auxid=sci.auxid;

          ++I;
        }
      }

      return(outfile);
    } // File readfile(const std::istream& is)

Here is the call graph for this function:

Here is the caller graph for this function: