BoundingBoxRefine.cpp

Go to the documentation of this file.

/*
 * Copyright (c) 2004 CSIRO ICT Centre
 *
 * $Id: BoundingBoxRefine.cpp 587 2004-12-03 15:06:33Z nch $
 */

/* BoundingBoxRefine is a Refine that reports
 * only those candidate pairs who are close enough to each other than
 * the bounding box around the angular threshold circle centred on one
 * point contains the other point. */

#include <iostream>
#include "String.h"

#include "BoundingBoxRefine.h"
#include "Object.h"
#include "ObjectPairConsumer.h"


String BoundingBoxRefine::s_name("Bounding Box Refine");


BoundingBoxRefine::BoundingBoxRefine(ObjectPairConsumer * mPPCons,
                                     ObjectPairConsumer * uPPCons)
    : Refine(mPPCons, uPPCons)
{}

BoundingBoxRefine::~BoundingBoxRefine()
{}

/* checks whether the bounding boxes of the two points intersect.
 * If so, then a match is passed through for reporting or further
 * refinement.  If not, the pair is rejected. */
bool BoundingBoxRefine::refine(Object const * redObject,
                               Object const * blueObject)
{
  bool result = false;

  double radius = Object::computeDistanceBound(redObject->getSD(),
                                              blueObject->getSD());

  // if the distance between the two is less than the sum of their bounds,
  // then reject the pair
  if (std::abs(redObject->getDec() - blueObject->getDec()) > radius)
  {
    reportNoMatch(redObject, blueObject);
  }
  else
  {
    // for the second parameter, either point will suffice.  The constraint
    // will be slightly stronger if we picker the point furthest from the
    // equator, but it's not really worth the effort.
    double raBound = Object::computeRACorrection(radius, redObject->getDec());

    double distance = blueObject->getRa() - redObject->getRa();
    if (distance < 0)
      distance = 360 + distance;
    if (distance > 180)
      distance = 360 - distance;
    if (distance > raBound)
    {
      reportNoMatch(redObject, blueObject);
    }
    else
    {
      result = reportMatch(redObject, blueObject);
    }
  }

  return result;
}