HtmIndex.hxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// CVS: $Revision: 7503 $, $Date: 2010-09-07 16:34:15 +0100 (Tue, 07 Sep 2010) $, $Author: EckhardSutorius $
//-----------------------------------------------------------------------------
#ifndef HTMINDEX_H
#define HTMINDEX_H

#include "htm.h"
#include <limits>

#include "Constants.h"
#include "Converter.h"
#include "DataMethod.hxx"
#include "Logger.h"
#include "StringOps.h"
//-----------------------------------------------------------------------------
template<typename DataType>
class HtmIndex : public DataMethod<DataType>
{
public:
  HtmIndex(int htmLevel)
  {
    DataMethod<DataType>::methodName = "HtmIndex";
    Logger log;
    log.addMessage("Making HTM indicies with level = "
                   + StringOps::NumToString(htmLevel));
    mHtmPtr = new SpatialIndex(htmLevel);
  }

  virtual ~HtmIndex() { delete mHtmPtr; }

  void setSource(IntMap& tbls, IntMap& cols, StringMap units)
  {
    raColno  = cols[0];
    decColno = cols[1];
    Converter converter;
    mDecFact = converter.getFactor(units[0], "Degrees");
    mRAFact  = converter.getFactor(units[1], "Degrees");
  }

  void setTarget(IntMap& cols, StringMap units)
  {
    htmColno = cols[0];
  }

  void doit(TableData<DataType>& data, int row1, int row2)
  {
    Logger log;
    // Calculate HTM indices
    for (int rowNo = row1; rowNo < row2; ++rowNo)
    {
      double ra, dec;
      data.value(raColno, rowNo, ra);
      data.value(decColno, rowNo, dec);

      // Assign to correct columns
      if (ra != ra or std::isinf(ra) or
          dec != dec  or std::isinf(dec))
      {
        string mess = "WCS error in htmID ";
        mess += StringOps::NumToString(htmColno) + " ";
        mess += StringOps::NumToString(rowNo) + " : ";
        mess += "value replaced with dbl_max";
        log.addError(mess);
        data.assign(htmColno, rowNo, 0.0);
      } else {
        data.assign(htmColno, rowNo, static_cast<double>(
                      mHtmPtr->idByPoint(ra * mRAFact, dec * mDecFact)));
      }
    }
  }

  // No data from raw data included
  void doit2(TableData<DataType>& rawData, TableData<DataType>& data, int row1, int row2) { }

private:
  int htmColno;
  int raColno;
  int decColno;
  double mDecFact;
  double mRAFact;
  SpatialIndex* mHtmPtr;
};
//-----------------------------------------------------------------------------
#endif
//-----------------------------------------------------------------------------
// Change log:
//
//  3-Jun-2004,  IAB: Original version.
// 22-Mar-2007,  RSC: Robustified against dodgy WCS values.
//  2-May-2007,  RSC: Reverted to original design of allowing dodgy WCS values.
//  7-Apr-2008, ETWS: Upgraded to use new detection table layout.