volInfo.cpp

 
#include <iostream>
#include <map>
#include <DGtal/base/Common.h>
#include <DGtal/io/readers/GenericReader.h>
#include <DGtal/helpers/StdDefs.h>
#include <DGtal/images/Image.h>
#include <DGtal/images/ImageContainerBySTLVector.h>
 
#include "CLI11.hpp"
 
 
using namespace std;
using namespace DGtal;
using namespace Z3i;
 
void missingParam ( std::string param )
{
  trace.error() <<" Parameter: "<<param<<" is required..";
  trace.info() <<std::endl;
  exit ( 1 );
}
 
 
int main(int argc, char**argv)
{
    // parse command line using CLI ----------------------------------------------
  CLI::App app;
  std::string inputFileName;
    
  app.description("Retreive information from vol file\n Basic usage: \n \tvolInfo  <volFileName>  ");
  app.add_option("-i,--input,1", inputFileName, "Input vol file." )->required()->check(CLI::ExistingFile);
  
  app.get_formatter()->column_width(40);
  CLI11_PARSE(app, argc, argv);
  // END parse command line using CLI ----------------------------------------------
 
  trace.beginBlock("Loading file");
  typedef ImageContainerBySTLVector<Z3i::Domain, unsigned char>  MyImageC;
 
  MyImageC  imageC = GenericReader<MyImageC>::import(inputFileName);
  
  trace.info()<<imageC<<std::endl;
  trace.info()<<"Scanning the values:"<<std::endl;
 
  std::map<unsigned char, size_t> values;
  for(auto v: imageC.range())
    values[ v ] ++;
  
  for(auto val: values)
    std::cout<<"\tvalue "<< (int)val.first<<": "<< val.second<<" voxels."<<std::endl;
  
  trace.endBlock();
  return 0;
}