Blobber.cpp

Go to the documentation of this file.

#include "Blobber.h"

#include <vector>

#include "Cube.h"
#include "IException.h"
#include "Progress.h"
#include "SpecialPixel.h"
#include "Table.h"
#include "TableField.h"

using std::vector;

namespace Isis {
  Blobber::Blobber() : _blobname("_undefined_"), _fieldname("_undefined_"),
    _name("Blob") { }
  Blobber::Blobber(const QString &blobname, const QString &fieldname,
                   const QString &name) : _blobname(blobname),
    _fieldname(fieldname), _name(name) {
  }

  Blobber::Blobber(Cube &cube, const QString &blobname,
                   const QString &fieldname,
                   const QString &name) :
    _blobname(blobname),
    _fieldname(fieldname),
    _name(name) {
    load(cube);
  }

  Blobber Blobber::deepcopy() const {
    Blobber myblob = *this;
    myblob._buf = _buf.copy();
    return (myblob);
  }

  void Blobber::load(const QString &filename) {
    Cube cube;
    cube.open(filename);
    load(cube);
    return;
  }


  void Blobber::load(Cube &cube) {
    Table tbl(getBlobName());
    cube.read(tbl);
    TableField data = tbl[0][getFieldName()];
    if (data.isDouble()) {
      loadDouble(tbl);
    }
    else if (data.isInteger()) {
      loadInteger(tbl);
    }
    else {
      QString msg = "Field type for " + getFieldName() +
                   " is not double or integer";
      throw IException(IException::Programmer, msg, _FILEINFO_);
    }
  }


  void Blobber::loadDouble(Table &tbl) {
    int nlines = tbl.Records();
    int nsamps = tbl[0][getFieldName()].size();
    BlobBuf pixels(nlines, nsamps);
    for (int i = 0 ; i < nlines ; i++) {
      vector<double> d = tbl[i][getFieldName()];
      for (unsigned int j = 0 ; j < d.size() ; j++) {
        pixels[i][j] = d[j];
      }
    }
    _buf = pixels;
  }

  void Blobber::loadInteger(Table &tbl) {
    int nlines = tbl.Records();
    int nsamps = tbl[0][getFieldName()].size();
    BlobBuf pixels(nlines, nsamps);
    for (int i = 0 ; i < nlines ; i++) {
      vector<int> d = tbl[i][getFieldName()];
      for (unsigned int j = 0 ; j < d.size(); j++) {
        pixels[i][j] = int2ToDouble(d[j]);
      }
    }
    _buf = pixels;
  }

  double Blobber::int2ToDouble(int value) const {
    if (value == NULL2) return NULL8;
    else if (value == LOW_REPR_SAT2) return LOW_REPR_SAT8;
    else if (value == LOW_INSTR_SAT2) return LOW_INSTR_SAT8;
    else if (value == HIGH_INSTR_SAT2) return HIGH_INSTR_SAT8;
    else if (value == HIGH_REPR_SAT2) return HIGH_REPR_SAT8;
    else return value;

  }


  double Blobber::int2ToDouble(unsigned int value) const {
    if (value == NULLUI4) return NULL8;
    else if (value == LOW_REPR_SATUI4) return LOW_REPR_SAT8;
    else if (value == LOW_INSTR_SATUI4) return LOW_INSTR_SAT8;
    else if (value == HIGH_INSTR_SATUI4) return HIGH_INSTR_SAT8;
    else if (value == HIGH_REPR_SATUI4) return HIGH_REPR_SAT8;
    else return value;

  }




}  //  end namespace Isis