usSequenceReader.h

/****************************************************************************
 *
 * This file is part of the ustk software.
 * Copyright (C) 2016 - 2017 by Inria. All rights reserved.
 *
 * This software is free software; you can redistribute it and/or
 * modify it under the terms of the GNU General Public License
 * ("GPL") version 2 as published by the Free Software Foundation.
 * See the file LICENSE.txt at the root directory of this source
 * distribution for additional information about the GNU GPL.
 *
 * For using ustk with software that can not be combined with the GNU
 * GPL, please contact Inria about acquiring a ViSP Professional
 * Edition License.
 *
 * This software was developed at:
 * Inria Rennes - Bretagne Atlantique
 * Campus Universitaire de Beaulieu
 * 35042 Rennes Cedex
 * France
 *
 * If you have questions regarding the use of this file, please contact
 * Inria at ustk@inria.fr
 *
 * This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
 * WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
 *
 * Authors:
 * Marc Pouliquen
 *
 *****************************************************************************/
 
#ifndef __usSequenceReader_h_
#define __usSequenceReader_h_
 
#include <algorithm>
#include <cstring>
#include <iostream>
#include <vector>
 //#include <unistd.h>
 
#include <visp3/core/vpConfig.h>
 
#ifdef VISP_HAVE_XML2
 
#include <visp3/core/vpException.h>
#include <visp3/core/vpImage.h>
#include <visp3/io/vpImageIo.h>
 
#include <visp3/ustk_core/usImageIo.h>
#include <visp3/ustk_core/usImageSettingsXmlParser.h>
 
template <class ImageType> class usSequenceReader
{
private:
  ImageType m_frame;
 
  double m_frameRate;
 
  long m_firstFrame;
  bool m_firstFrameIsSet;
 
  long m_lastFrame;
  bool m_lastFrameIsSet;
 
  long m_frameCount;
 
  std::string m_sequenceFileName;
  std::string m_genericImageFileName;
  bool m_fileNameIsSet;
 
  bool is_open;
 
  bool m_enableLoopCycling;
  int loopIncrement; // -1 or +1 depending on witch side we acquire
 
  usImageSettingsXmlParser m_xmlParser;
 
  std::vector<uint64_t> m_timestamps;
 
public:
  usSequenceReader();
 
  virtual ~usSequenceReader();
 
  // get images in grabber style
  void acquire(ImageType &image);
  void acquire(ImageType &image, uint64_t &timestamp);
 
  inline bool end()
  {
    if (!is_open) {
      open(m_frame);
      m_frameCount--;
    }
  //
    if (m_enableLoopCycling) {
      if (loopIncrement == 1) {
        if (m_frameCount > m_lastFrame) {
          m_frameCount -= 2;
          loopIncrement = -1;
          return false;
        }
      }
      else if (loopIncrement == -1) {
        if (m_frameCount < m_firstFrame) {
          m_frameCount += 2;
          loopIncrement = 1;
          return false;
        }
      }
    }
    else if (m_frameCount > m_lastFrame || m_frameCount < m_firstFrame)
      return true;
    return false;
  }
 
// get image by its number in the sequence
  void getFrame(ImageType &image, int index);
 
  // get image by its number in the sequence
  ImageType &getFrame(int index);
 
  // total number of frames in the sequence
  long getFrameCount();
 
  // attributes getters/setters
  double getFrameRate() const { return m_frameRate; }
  long getImageNumber() const { return m_frameCount; }
 
  // timestamps vector getter
  std::vector<uint64_t> getSequenceTimestamps() const { return m_timestamps; }
 
  // get the xml parser : usefull to acess all the image settings contained in it
  usImageSettingsXmlParser getXmlParser();
 
  void open(ImageType &image, uint64_t &timestamp);
  void open(ImageType &image);
 
  void setFirstFrameIndex(long firstIndex);
  void setLastFrameIndex(long lastIndex);
  void setLoopCycling(bool activateLoopCycling);
  void setSequenceFileName(const std::string &sequenceFileName);
};
 
/****************************************************************************
* Template implementations.
****************************************************************************/
 
template <class ImageType>
usSequenceReader<ImageType>::usSequenceReader()
  : m_frame(), m_frameRate(0.0), m_firstFrame(0), m_firstFrameIsSet(false), m_lastFrame(0), m_lastFrameIsSet(false),
  m_frameCount(0), m_sequenceFileName(""), m_genericImageFileName(""), m_fileNameIsSet(false), is_open(false),
  m_enableLoopCycling(false), loopIncrement(1)
{ }
 
template <class ImageType> usSequenceReader<ImageType>::~usSequenceReader() { }
 
template <class ImageType> void usSequenceReader<ImageType>::setSequenceFileName(const std::string &sequenceFileName)
{
  m_sequenceFileName = sequenceFileName;
  m_fileNameIsSet = true;
}
 
template <class ImageType> void usSequenceReader<ImageType>::setFirstFrameIndex(long firstIndex)
{
  m_firstFrame = firstIndex;
  m_firstFrameIsSet = true;
}
 
template <class ImageType> void usSequenceReader<ImageType>::setLastFrameIndex(long lastIndex)
{
  m_lastFrame = lastIndex;
  m_lastFrameIsSet = true;
}
 
template <class ImageType> void usSequenceReader<ImageType>::open(ImageType &image)
{
  (void)image;
  throw(vpException(vpException::notImplementedError));
}
 
template <class ImageType> void usSequenceReader<ImageType>::open(ImageType &image, uint64_t &timestamp)
{
  (void)image;
  (void)timestamp;
  throw(vpException(vpException::notImplementedError));
}
 
template <> inline void usSequenceReader<usImagePreScan2D<unsigned char> >::open(usImagePreScan2D<unsigned char> &image)
{
  if (!m_fileNameIsSet) {
    throw(vpException(vpException::badValue, "Sequence settings file name not set"));
  }
 
  m_xmlParser.parse(m_sequenceFileName);
  if (m_xmlParser.getImageType() != us::PRESCAN_2D)
    throw(vpException(vpException::badValue), "trying to open a non-usImagePreScan2D image !");
 
  setFirstFrameIndex(m_xmlParser.getSequenceStartNumber());
  setLastFrameIndex(m_xmlParser.getSequenceStopNumber());
  m_frameRate = m_xmlParser.getSequenceFrameRate();
  m_genericImageFileName = m_xmlParser.getImageFileName();
 
  // saving the settings for all the pre-scan sequence
  m_frame.setTransducerRadius(m_xmlParser.getTransducerSettings().getTransducerRadius());
  m_frame.setScanLinePitch(m_xmlParser.getTransducerSettings().getScanLinePitch());
  m_frame.setScanLineNumber(m_xmlParser.getTransducerSettings().getScanLineNumber());
  m_frame.setTransducerConvexity(m_xmlParser.getTransducerSettings().isTransducerConvex());
  m_frame.setAxialResolution(m_xmlParser.getAxialResolution());
  m_frame.setSamplingFrequency(m_xmlParser.getTransducerSettings().getSamplingFrequency());
  m_frame.setTransmitFrequency(m_xmlParser.getTransducerSettings().getTransmitFrequency());
 
  // Reading image
  char buffer[FILENAME_MAX];
  snprintf(buffer, FILENAME_MAX, m_genericImageFileName.c_str(), m_firstFrame);
  std::string imageFileName;
  // check timestamp
  std::vector<std::string> splitName =
    vpIoTools::splitChain(vpIoTools::splitChain(m_genericImageFileName, std::string("/")).back(), std::string("."));
 
  std::string parentName = vpIoTools::getParent(m_sequenceFileName);
  if (!parentName.empty()) {
    parentName = parentName + vpIoTools::path("/");
  }
  std::vector<std::string> splitNameDirs = vpIoTools::splitChain(m_genericImageFileName, std::string("/"));
  if (splitNameDirs.size() > 1) { // images are in a subdirectory of the one conatining the xml file
    for (unsigned int i = 0; i < splitNameDirs.size() - 1; i++)
      parentName += (splitNameDirs.at(i) + vpIoTools::path("/"));
  }
  if (splitName.size() == 2) { // no timestamp : image0002.png for example
    imageFileName = parentName + vpIoTools::splitChain(std::string(buffer), std::string("/")).back();
  }
  else if (splitName.size() == 3) { // timestamp included : image0002.156464063.png for example
    std::vector<std::string> dirFiles;
 
    if (vpIoTools::checkDirectory(vpIoTools::getParent(m_sequenceFileName))) { // correct path
      dirFiles = vpIoTools::getDirFiles(parentName);
 
      if (dirFiles.size() == (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() +
                                            2)) { // case of xml file in same directory
        for (unsigned int i = 0; i < dirFiles.size(); i++) {
          if (vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1) == std::string("xml"))
            dirFiles.erase(dirFiles.begin() + i);
        }
      }
      if (dirFiles.size() !=
          (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() + 1))
        throw(vpException(vpException::fatalError, "For imgage sequnces with timeStamps, the directory must contain "
                          "only the entire image sequence (no additionnal files allowed)"));
 
      // getting the all the timestamps of the sequence
      unsigned int i = 0;
      while (i < dirFiles.size()) {
        uint64_t timestamp_temp;
        std::istringstream(vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1)) >> timestamp_temp;
        m_timestamps.push_back(timestamp_temp);
        i++;
      }
      imageFileName = parentName + dirFiles.front();
    }
    else { // path not correct
      throw(vpException(vpException::fatalError), "Sequence filename incorrect !");
    }
  }
 
  vpImageIo::read(image, imageFileName);
 
  m_frame.setDepth((image.getHeight() - 1) * m_frame.getAxialResolution());
 
  // workaround to prevent the resize by setImagePreScanSettings (having a scanline number at 0 if not precised in the
  // xml)
  int scanlineNumber = image.getWidth();
  image.setImagePreScanSettings(m_frame);
  image.setScanLineNumber(scanlineNumber);
 
  m_frameCount = m_firstFrame + 1;
  is_open = true;
}
 
template <>
inline void usSequenceReader<usImagePreScan2D<unsigned char> >::open(usImagePreScan2D<unsigned char> &image,
                                                                     uint64_t &timestamp)
{
  if (!m_fileNameIsSet) {
    throw(vpException(vpException::badValue, "Sequence settings file name not set"));
  }
 
  m_xmlParser.parse(m_sequenceFileName);
  if (m_xmlParser.getImageType() != us::PRESCAN_2D)
    throw(vpException(vpException::badValue), "trying to open a non-usImagePreScan2D image !");
 
  setFirstFrameIndex(m_xmlParser.getSequenceStartNumber());
  setLastFrameIndex(m_xmlParser.getSequenceStopNumber());
  m_frameRate = m_xmlParser.getSequenceFrameRate();
  m_genericImageFileName = m_xmlParser.getImageFileName();
 
  // saving the settings for all the pre-scan sequence
  m_frame.setTransducerRadius(m_xmlParser.getTransducerSettings().getTransducerRadius());
  m_frame.setScanLinePitch(m_xmlParser.getTransducerSettings().getScanLinePitch());
  m_frame.setScanLineNumber(m_xmlParser.getTransducerSettings().getScanLineNumber());
  m_frame.setTransducerConvexity(m_xmlParser.getTransducerSettings().isTransducerConvex());
  m_frame.setAxialResolution(m_xmlParser.getAxialResolution());
  m_frame.setSamplingFrequency(m_xmlParser.getTransducerSettings().getSamplingFrequency());
  m_frame.setTransmitFrequency(m_xmlParser.getTransducerSettings().getTransmitFrequency());
 
  // Reading image
  char buffer[FILENAME_MAX];
  snprintf(buffer, FILENAME_MAX, m_genericImageFileName.c_str(), m_firstFrame);
  std::string imageFileName;
  // check timestamp
  std::vector<std::string> splitName =
    vpIoTools::splitChain(vpIoTools::splitChain(m_genericImageFileName, std::string("/")).back(), std::string("."));
 
  std::string parentName = vpIoTools::getParent(m_sequenceFileName);
  if (!parentName.empty()) {
    parentName = parentName + vpIoTools::path("/");
  }
  std::vector<std::string> splitNameDirs = vpIoTools::splitChain(m_genericImageFileName, std::string("/"));
  if (splitNameDirs.size() > 1) { // images are in a subdirectory of the one conatining the xml file
    for (unsigned int i = 0; i < splitNameDirs.size() - 1; i++)
      parentName += (splitNameDirs.at(i) + vpIoTools::path("/"));
  }
  if (splitName.size() == 2) { // no timestamp : image0002.png for example
    imageFileName = parentName + vpIoTools::splitChain(std::string(buffer), std::string("/")).back();
    timestamp = 0;
    m_timestamps.clear();
    m_timestamps.push_back(timestamp);
  }
  else if (splitName.size() == 3) { // timestamp included : image0002.156464063.png for example
    std::vector<std::string> dirFiles;
 
    if (vpIoTools::checkDirectory(vpIoTools::getParent(m_sequenceFileName))) { // correct path
      dirFiles = vpIoTools::getDirFiles(parentName);
 
      if (dirFiles.size() == (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() +
                                            2)) { // case of xml file in same directory
        for (unsigned int i = 0; i < dirFiles.size(); i++) {
          if (vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1) == std::string("xml"))
            dirFiles.erase(dirFiles.begin() + i);
        }
      }
      if (dirFiles.size() !=
          (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() + 1))
        throw(vpException(vpException::fatalError, "For imgage sequnces with timeStamps, the directory must contain "
                          "only the entire image sequence (no additionnal files allowed)"));
 
      // getting the all the timestamps of the sequence
      unsigned int i = 0;
      while (i < dirFiles.size()) {
        uint64_t timestamp_temp;
        std::istringstream(vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1)) >> timestamp_temp;
        m_timestamps.push_back(timestamp_temp);
        i++;
      }
      std::istringstream(vpIoTools::splitChain(dirFiles.at(0), std::string(".")).at(1)) >> timestamp;
 
      imageFileName = parentName + dirFiles.front();
    }
    else { // path not correct
      throw(vpException(vpException::fatalError), "Sequence filename incorrect !");
    }
  }
 
  vpImageIo::read(image, imageFileName);
 
  m_frame.setDepth((image.getHeight() - 1) * m_frame.getAxialResolution());
 
  // workaround to prevent the resize by setImagePreScanSettings (having a scanline number at 0 if not precised in the
  // xml)
  int scanlineNumber = image.getWidth();
  image.setImagePreScanSettings(m_frame);
  image.setScanLineNumber(scanlineNumber);
 
  m_frameCount = m_firstFrame + 1;
  is_open = true;
}
 
template <>
inline void usSequenceReader<usImagePostScan2D<unsigned char> >::open(usImagePostScan2D<unsigned char> &image)
{
  if (!m_fileNameIsSet) {
    throw(vpException(vpException::badValue, "Sequence settings file name not set"));
  }
 
  m_xmlParser.parse(m_sequenceFileName);
  if (m_xmlParser.getImageType() != us::POSTSCAN_2D)
    throw(vpException(vpException::badValue), "trying to open a non-usImagePostScan2D image !");
 
  setFirstFrameIndex(m_xmlParser.getSequenceStartNumber());
  setLastFrameIndex(m_xmlParser.getSequenceStopNumber());
  m_frameRate = m_xmlParser.getSequenceFrameRate();
  m_genericImageFileName = m_xmlParser.getImageFileName();
 
  // saving the settings for all the post scan sequence
  m_frame.setTransducerRadius(m_xmlParser.getTransducerSettings().getTransducerRadius());
  m_frame.setScanLinePitch(m_xmlParser.getTransducerSettings().getScanLinePitch());
  m_frame.setScanLineNumber(m_xmlParser.getTransducerSettings().getScanLineNumber());
  m_frame.setTransducerConvexity(m_xmlParser.getTransducerSettings().isTransducerConvex());
  m_frame.setWidthResolution(m_xmlParser.getWidthResolution());
  m_frame.setHeightResolution(m_xmlParser.getHeightResolution());
  m_frame.setSamplingFrequency(m_xmlParser.getTransducerSettings().getSamplingFrequency());
  m_frame.setTransmitFrequency(m_xmlParser.getTransducerSettings().getTransmitFrequency());
 
  // Reading image
  char buffer[FILENAME_MAX];
  snprintf(buffer, FILENAME_MAX, m_genericImageFileName.c_str(), m_firstFrame);
  std::string imageFileName;
  // check timestamp
  std::vector<std::string> splitName =
    vpIoTools::splitChain(vpIoTools::splitChain(m_genericImageFileName, std::string("/")).back(), std::string("."));
 
  std::string parentName = vpIoTools::getParent(m_sequenceFileName);
  if (!parentName.empty()) {
    parentName = parentName + vpIoTools::path("/");
  }
  std::vector<std::string> splitNameDirs = vpIoTools::splitChain(m_genericImageFileName, std::string("/"));
  if (splitNameDirs.size() > 1) { // images are in a subdirectory of the one conatining the xml file
    for (unsigned int i = 0; i < splitNameDirs.size() - 1; i++)
      parentName += (splitNameDirs.at(i) + vpIoTools::path("/"));
  }
  if (splitName.size() == 2) { // no timestamp : image0002.png for example
    imageFileName = parentName + vpIoTools::splitChain(std::string(buffer), std::string("/")).back();
  }
  else if (splitName.size() == 3) { // timestamp included : image0002.156464063.png for example
    std::vector<std::string> dirFiles;
 
    if (vpIoTools::checkDirectory(vpIoTools::getParent(m_sequenceFileName))) { // correct path
      dirFiles = vpIoTools::getDirFiles(parentName);
 
      if (dirFiles.size() == (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() +
                                            2)) { // case of xml file in same directory
        for (unsigned int i = 0; i < dirFiles.size(); i++) {
          if (vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1) == std::string("xml"))
            dirFiles.erase(dirFiles.begin() + i);
        }
      }
 
      if (dirFiles.size() !=
          (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() + 1))
        throw(vpException(vpException::fatalError, "For imgage sequnces with timeStamps, the directory must contain "
                          "only the entire image sequence (no additionnal files allowed)"));
 
      // getting the all the timestamps of the sequence
      unsigned int i = 0;
      while (i < dirFiles.size()) {
        uint64_t timestamp_temp;
        std::istringstream(vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1)) >> timestamp_temp;
        m_timestamps.push_back(timestamp_temp);
        i++;
      }
 
      imageFileName = parentName + dirFiles.front();
    }
    else { // path not correct
      throw(vpException(vpException::fatalError), "Sequence filename incorrect !");
    }
  }
 
  vpImageIo::read(image, imageFileName);
  m_frame.setDepth(image.getHeight() - 1 * m_frame.getHeightResolution());
 
  image.setTransducerRadius(m_frame.getTransducerRadius());
  image.setScanLinePitch(m_frame.getScanLinePitch());
  image.setScanLineNumber(m_frame.getScanLineNumber());
  image.setTransducerConvexity(m_frame.isTransducerConvex());
  image.setWidthResolution(m_frame.getWidthResolution());
  image.setHeightResolution(m_frame.getHeightResolution());
  image.setSamplingFrequency(m_frame.getSamplingFrequency());
  image.setTransmitFrequency(m_frame.getTransmitFrequency());
  image.setDepth(m_frame.getDepth());
 
  m_frameCount = m_firstFrame + 1;
  is_open = true;
}
 
template <>
inline void usSequenceReader<usImagePostScan2D<unsigned char> >::open(usImagePostScan2D<unsigned char> &image,
                                                                      uint64_t &timestamp)
{
  if (!m_fileNameIsSet) {
    throw(vpException(vpException::badValue, "Sequence settings file name not set"));
  }
 
  m_xmlParser.parse(m_sequenceFileName);
  if (m_xmlParser.getImageType() != us::POSTSCAN_2D)
    throw(vpException(vpException::badValue), "trying to open a non-usImagePostScan2D image !");
 
  setFirstFrameIndex(m_xmlParser.getSequenceStartNumber());
  setLastFrameIndex(m_xmlParser.getSequenceStopNumber());
  m_frameRate = m_xmlParser.getSequenceFrameRate();
  m_genericImageFileName = m_xmlParser.getImageFileName();
 
  // saving the settings for all the post scan sequence
  m_frame.setTransducerRadius(m_xmlParser.getTransducerSettings().getTransducerRadius());
  m_frame.setScanLinePitch(m_xmlParser.getTransducerSettings().getScanLinePitch());
  m_frame.setScanLineNumber(m_xmlParser.getTransducerSettings().getScanLineNumber());
  m_frame.setTransducerConvexity(m_xmlParser.getTransducerSettings().isTransducerConvex());
  m_frame.setWidthResolution(m_xmlParser.getWidthResolution());
  m_frame.setHeightResolution(m_xmlParser.getHeightResolution());
  m_frame.setSamplingFrequency(m_xmlParser.getTransducerSettings().getSamplingFrequency());
  m_frame.setTransmitFrequency(m_xmlParser.getTransducerSettings().getTransmitFrequency());
 
  // Reading image
  char buffer[FILENAME_MAX];
  snprintf(buffer, FILENAME_MAX, m_genericImageFileName.c_str(), m_firstFrame);
  std::string imageFileName;
  // check timestamp
  std::vector<std::string> splitName =
    vpIoTools::splitChain(vpIoTools::splitChain(m_genericImageFileName, std::string("/")).back(), std::string("."));
 
  std::string parentName = vpIoTools::getParent(m_sequenceFileName);
  if (!parentName.empty()) {
    parentName = parentName + vpIoTools::path("/");
  }
  std::vector<std::string> splitNameDirs = vpIoTools::splitChain(m_genericImageFileName, std::string("/"));
  if (splitNameDirs.size() > 1) { // images are in a subdirectory of the one conatining the xml file
    for (unsigned int i = 0; i < splitNameDirs.size() - 1; i++)
      parentName += (splitNameDirs.at(i) + vpIoTools::path("/"));
  }
  if (splitName.size() == 2) { // no timestamp : image0002.png for example
    imageFileName = parentName + vpIoTools::splitChain(std::string(buffer), std::string("/")).back();
    timestamp = 0;
  }
  else if (splitName.size() == 3) { // timestamp included : image0002.156464063.png for example
    std::istringstream(splitName.at(1)) >> timestamp;
    std::vector<std::string> dirFiles;
 
    if (vpIoTools::checkDirectory(vpIoTools::getParent(m_sequenceFileName))) { // correct path
      dirFiles = vpIoTools::getDirFiles(parentName);
 
      if (dirFiles.size() == (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() +
                                            2)) { // case of xml file in same directory
        for (unsigned int i = 0; i < dirFiles.size(); i++) {
          if (vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1) == std::string("xml"))
            dirFiles.erase(dirFiles.begin() + i);
        }
      }
 
      std::istringstream(vpIoTools::splitChain(dirFiles.at(0), std::string(".")).at(1)) >> timestamp;
 
      if (dirFiles.size() !=
          (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() + 1))
        throw(vpException(vpException::fatalError, "For imgage sequnces with timeStamps, the directory must contain "
                          "only the entire image sequence (no additionnal files allowed)"));
 
      // getting the all the timestamps of the sequence
      unsigned int i = 0;
      while (i < dirFiles.size()) {
        uint64_t timestamp_temp;
        std::istringstream(vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1)) >> timestamp_temp;
        m_timestamps.push_back(timestamp_temp);
        i++;
      }
 
      imageFileName = parentName + dirFiles.front();
    }
    else { // path not correct
      throw(vpException(vpException::fatalError), "Sequence filename incorrect !");
    }
  }
  vpImageIo::read(image, imageFileName);
  m_frame.setDepth(image.getHeight() - 1 * m_frame.getHeightResolution());
 
  image.setTransducerRadius(m_frame.getTransducerRadius());
  image.setScanLinePitch(m_frame.getScanLinePitch());
  image.setScanLineNumber(m_frame.getScanLineNumber());
  image.setTransducerConvexity(m_frame.isTransducerConvex());
  image.setWidthResolution(m_frame.getWidthResolution());
  image.setHeightResolution(m_frame.getHeightResolution());
  image.setSamplingFrequency(m_frame.getSamplingFrequency());
  image.setTransmitFrequency(m_frame.getTransmitFrequency());
  image.setDepth(m_frame.getDepth());
 
  m_frameCount = m_firstFrame + 1;
  is_open = true;
}
 
template <class ImageType> void usSequenceReader<ImageType>::acquire(ImageType &image)
{
  if (!is_open) {
    this->open(image);
    return;
  }
 
// Reading image
  char buffer[FILENAME_MAX];
  snprintf(buffer, FILENAME_MAX, m_genericImageFileName.c_str(), m_frameCount);
  std::string imageFileName;
  // check timestamp
  std::vector<std::string> splitName =
    vpIoTools::splitChain(vpIoTools::splitChain(m_genericImageFileName, std::string("/")).back(), std::string("."));
 
  std::string parentName = vpIoTools::getParent(m_sequenceFileName);
  if (!parentName.empty()) {
    parentName = parentName + vpIoTools::path("/");
  }
  std::vector<std::string> splitNameDirs = vpIoTools::splitChain(m_genericImageFileName, std::string("/"));
  if (splitNameDirs.size() > 1) { // images are in a subdirectory of the one containing the xml file
    for (unsigned int i = 0; i < splitNameDirs.size() - 1; i++)
      parentName += (splitNameDirs.at(i) + vpIoTools::path("/"));
  }
  if (splitName.size() == 2) { // no timestamp : image0002.png for example
    imageFileName = parentName + vpIoTools::splitChain(std::string(buffer), std::string("/")).back();
  }
  else if (splitName.size() == 3) { // timestamp included : image0002.156464063.png for example
    std::vector<std::string> dirFiles;
 
    if (vpIoTools::checkDirectory(vpIoTools::getParent(m_sequenceFileName))) { // correct path
      dirFiles = vpIoTools::getDirFiles(parentName);
 
      if (dirFiles.size() == (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() +
                                            2)) { // case of xml file in same directory
        for (unsigned int i = 0; i < dirFiles.size(); i++) {
          if (vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1) == std::string("xml"))
            dirFiles.erase(dirFiles.begin() + i);
        }
      }
 
      if (dirFiles.size() !=
          (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() + 1))
        throw(vpException(vpException::fatalError, "For imgage sequnces with timeStamps, the directory must contain "
                          "only the entire image sequence (no additionnal files allowed)"));
 
      imageFileName = parentName + dirFiles.at(m_frameCount);
    }
    else { // path not correct
      throw(vpException(vpException::fatalError), "Sequence filename incorrect !");
    }
  }
  vpImageIo::read(image, imageFileName);
  image.setTransducerSettings(m_frame);
  image.setScanLineNumber(image.getWidth());
  image.setDepth(m_frame.getDepth());
 
  m_frameCount += loopIncrement;
}
 
template <class ImageType> void usSequenceReader<ImageType>::acquire(ImageType &image, uint64_t &timestamp)
{
  if (!is_open) {
    this->open(image, timestamp);
    return;
  }
// Reading image
  char buffer[FILENAME_MAX];
  snprintf(buffer, FILENAME_MAX, m_genericImageFileName.c_str(), m_frameCount);
  std::string imageFileName;
  // check timestamp
  std::vector<std::string> splitName =
    vpIoTools::splitChain(vpIoTools::splitChain(m_genericImageFileName, std::string("/")).back(), std::string("."));
 
  std::string parentName = vpIoTools::getParent(m_sequenceFileName);
  if (!parentName.empty()) {
    parentName = parentName + vpIoTools::path("/");
  }
  std::vector<std::string> splitNameDirs = vpIoTools::splitChain(m_genericImageFileName, std::string("/"));
  if (splitNameDirs.size() > 1) { // images are in a subdirectory of the one conatining the xml file
    for (unsigned int i = 0; i < splitNameDirs.size() - 1; i++)
      parentName += (splitNameDirs.at(i) + vpIoTools::path("/"));
  }
  if (splitName.size() == 2) { // no timestamp : image0002.png for example
    imageFileName = parentName + vpIoTools::splitChain(std::string(buffer), std::string("/")).back();
    timestamp = 0;
  }
  else if (splitName.size() == 3) { // timestamp included : image0002.156464063.png for example
 
    std::vector<std::string> dirFiles;
 
    if (vpIoTools::checkDirectory(vpIoTools::getParent(m_sequenceFileName))) { // correct path
      dirFiles = vpIoTools::getDirFiles(parentName);
 
      if (dirFiles.size() == (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() +
                                            2)) { // case of xml file in same directory
        for (unsigned int i = 0; i < dirFiles.size(); i++) {
          if (vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1) == std::string("xml"))
            dirFiles.erase(dirFiles.begin() + i);
        }
      }
 
      std::istringstream(vpIoTools::splitChain(dirFiles.at(m_frameCount), std::string(".")).at(1)) >> timestamp;
 
      if (dirFiles.size() !=
          (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() + 1))
        throw(vpException(vpException::fatalError, "For imgage sequnces with timeStamps, the directory must contain "
                          "only the entire image sequence (no additionnal files allowed)"));
 
      imageFileName = parentName + dirFiles.at(m_frameCount);
    }
    else { // path not correct
      throw(vpException(vpException::fatalError), "Sequence filename incorrect !");
    }
  }
 
  vpImageIo::read(image, imageFileName);
  image.setTransducerSettings(m_frame);
  image.setScanLineNumber(image.getWidth());
  image.setDepth(m_frame.getDepth());
 
  m_frameCount += loopIncrement;
}
 
template <>
inline void usSequenceReader<usImagePreScan2D<unsigned char> >::acquire(usImagePreScan2D<unsigned char> &image)
{
  if (!is_open) {
    this->open(image);
    return;
  }
// Reading image
  char buffer[FILENAME_MAX];
  snprintf(buffer, FILENAME_MAX, m_genericImageFileName.c_str(), m_frameCount);
  std::string imageFileName;
  // check timestamp
  std::vector<std::string> splitName =
    vpIoTools::splitChain(vpIoTools::splitChain(m_genericImageFileName, std::string("/")).back(), std::string("."));
 
  std::string parentName = vpIoTools::getParent(m_sequenceFileName);
  if (!parentName.empty()) {
    parentName = parentName + vpIoTools::path("/");
  }
  std::vector<std::string> splitNameDirs = vpIoTools::splitChain(m_genericImageFileName, std::string("/"));
  if (splitNameDirs.size() > 1) { // images are in a subdirectory of the one conatining the xml file
    for (unsigned int i = 0; i < splitNameDirs.size() - 1; i++)
      parentName += (splitNameDirs.at(i) + vpIoTools::path("/"));
  }
  if (splitName.size() == 2) { // no timestamp : image0002.png for example
    imageFileName = parentName + vpIoTools::splitChain(std::string(buffer), std::string("/")).back();
  }
  else if (splitName.size() == 3) { // timestamp included : image0002.156464063.png for example
    std::vector<std::string> dirFiles;
 
    if (vpIoTools::checkDirectory(vpIoTools::getParent(m_sequenceFileName))) { // correct path
      dirFiles = vpIoTools::getDirFiles(parentName);
 
      if (dirFiles.size() == (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() +
                                            2)) { // case of xml file in same directory
        for (unsigned int i = 0; i < dirFiles.size(); i++) {
          if (vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1) == std::string("xml"))
            dirFiles.erase(dirFiles.begin() + i);
        }
      }
 
      if (dirFiles.size() !=
          (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() + 1))
        throw(vpException(vpException::fatalError, "For imgage sequnces with timeStamps, the directory must contain "
                          "only the entire image sequence (no additionnal files allowed)"));
 
      imageFileName = parentName + dirFiles.at(m_frameCount);
    }
    else { // path not correct
      throw(vpException(vpException::fatalError), "Sequence filename incorrect !");
    }
  }
 
  vpImageIo::read(image, imageFileName);
  image.setTransducerSettings(m_frame);
  image.setScanLineNumber(image.getWidth());
  image.setDepth(m_frame.getDepth());
  image.setAxialResolution(m_frame.getAxialResolution());
  image.setSamplingFrequency(m_frame.getSamplingFrequency());
  image.setTransmitFrequency(m_frame.getTransmitFrequency());
 
  m_frameCount += loopIncrement;
}
 
template <>
inline void usSequenceReader<usImagePreScan2D<unsigned char> >::acquire(usImagePreScan2D<unsigned char> &image,
                                                                        uint64_t &timestamp)
{
  if (!is_open) {
    this->open(image, timestamp);
    return;
  }
// Reading image
  char buffer[FILENAME_MAX];
  snprintf(buffer, FILENAME_MAX, m_genericImageFileName.c_str(), m_frameCount);
  std::string imageFileName;
  // check timestamp
  std::vector<std::string> splitName =
    vpIoTools::splitChain(vpIoTools::splitChain(m_genericImageFileName, std::string("/")).back(), std::string("."));
 
  std::string parentName = vpIoTools::getParent(m_sequenceFileName);
  if (!parentName.empty()) {
    parentName = parentName + vpIoTools::path("/");
  }
  std::vector<std::string> splitNameDirs = vpIoTools::splitChain(m_genericImageFileName, std::string("/"));
  if (splitNameDirs.size() > 1) { // images are in a subdirectory of the one conatining the xml file
    for (unsigned int i = 0; i < splitNameDirs.size() - 1; i++)
      parentName += (splitNameDirs.at(i) + vpIoTools::path("/"));
  }
  if (splitName.size() == 2) { // no timestamp : image0002.png for example
    imageFileName = parentName + vpIoTools::splitChain(std::string(buffer), std::string("/")).back();
    timestamp = 0;
  }
  else if (splitName.size() == 3) { // timestamp included : image0002.156464063.png for example
    std::vector<std::string> dirFiles;
 
    if (vpIoTools::checkDirectory(vpIoTools::getParent(m_sequenceFileName))) { // correct path
      dirFiles = vpIoTools::getDirFiles(parentName);
 
      std::sort(dirFiles.begin(), dirFiles.end());
      if (dirFiles.size() == (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() +
                                            2)) { // case of xml file in same directory
        for (unsigned int i = 0; i < dirFiles.size(); i++) {
          if (vpIoTools::splitChain(dirFiles.at(i), std::string(".")).at(1) == std::string("xml"))
            dirFiles.erase(dirFiles.begin() + i);
        }
      }
      std::istringstream(vpIoTools::splitChain(dirFiles.at(m_frameCount), std::string(".")).at(1)) >> timestamp;
 
      if (dirFiles.size() !=
          (unsigned int)(m_xmlParser.getSequenceStopNumber() - m_xmlParser.getSequenceStartNumber() + 1))
        throw(vpException(vpException::fatalError, "For imgage sequnces with timeStamps, the directory must contain "
                          "only the entire image sequence (no additionnal files allowed)"));
 
      imageFileName = parentName + dirFiles.at(m_frameCount);
    }
    else { // path not correct
      throw(vpException(vpException::fatalError), "Sequence filename incorrect !");
    }
  }
 
  vpImageIo::read(image, imageFileName);
  image.setTransducerSettings(m_frame);
  image.setScanLineNumber(image.getWidth());
  image.setDepth(m_frame.getDepth());
  image.setAxialResolution(m_frame.getAxialResolution());
  image.setSamplingFrequency(m_frame.getSamplingFrequency());
  image.setTransmitFrequency(m_frame.getTransmitFrequency());
 
  m_frameCount += loopIncrement;
}
 
template <class ImageType> void usSequenceReader<ImageType>::getFrame(ImageType &image, int index)
{
  if (!is_open) {
    open(image);
    return;
  }
  if (index < m_firstFrame || index > m_lastFrame) {
    throw(vpException(vpException::badValue, "position out of range"));
  }
 
// Reading image
  char buffer[FILENAME_MAX];
  snprintf(buffer, FILENAME_MAX, m_genericImageFileName.c_str(), index);
  std::string imageFileName = buffer;
 
  vpImageIo::read(image, imageFileName);
  image.setImageSettins(m_frame);
}
 
template <class ImageType> ImageType &usSequenceReader<ImageType>::getFrame(int index)
{
  ImageType image;
  getFrame(image, index);
  return image;
}
 
template <class ImageType> void usSequenceReader<ImageType>::setLoopCycling(bool activateLoopCycling)
{
  m_enableLoopCycling = activateLoopCycling;
}
 
template <class ImageType> long usSequenceReader<ImageType>::getFrameCount()
{
  if (is_open)
    return m_lastFrame - m_firstFrame + 1;
  throw(vpException(vpException::fatalError, "Cannot compute the frame number of the sequence before opening it !"));
}
 
template <class ImageType> usImageSettingsXmlParser usSequenceReader<ImageType>::getXmlParser()
{
  if (is_open)
    return m_xmlParser;
  throw(vpException(vpException::fatalError, "Sequence not opened, xml parser is empty !"));
}
 
#endif
#endif // __usSequenceReader_h_