HtmIndex.hxx

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// CVS: $Revision: 4695 $, $Date: 2008-04-09 16:18:42 +0100 (Wed, 09 Apr 2008) $, $Author: EckhardSutorius $
//-----------------------------------------------------------------------------
#ifndef HTMINDEX_H
#define HTMINDEX_H

#include "htm.h"

#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)
  {
    // Calculate HTM indices
    for (int rowNo = row1; rowNo < row2; ++rowNo)
    {
      double ra, dec;
      data.value(raColno, rowNo, ra);
      data.value(decColno, rowNo, dec);

      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.

---