usNeedleDetectionTools.h

/****************************************************************************
 *
 * This file is part of the ustk software.
 * Copyright (C) 2016 - 2017 by Inria. All rights reserved.
 *
 * This software is free software; you can redistribute it and/or
 * modify it under the terms of the GNU General Public License
 * ("GPL") version 2 as published by the Free Software Foundation.
 * See the file LICENSE.txt at the root directory of this source
 * distribution for additional information about the GNU GPL.
 *
 * For using ustk with software that can not be combined with the GNU
 * GPL, please contact Inria about acquiring a ViSP Professional
 * Edition License.
 *
 * This software was developed at:
 * Inria Rennes - Bretagne Atlantique
 * Campus Universitaire de Beaulieu
 * 35042 Rennes Cedex
 * France
 *
 * If you have questions regarding the use of this file, please contact
 * Inria at ustk@inria.fr
 *
 * This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
 * WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
 *
 * Authors:
 * Pierre Chatelain
 * Alexandre Krupa
 *
 *****************************************************************************/
 
#ifndef __usNeedleDetectionTools_h_
#define __usNeedleDetectionTools_h_
 
#include <algorithm>
#include <float.h>
#include <iostream>
#include <math.h>
#include <stdio.h>
#include <string.h>
#include <vector>
 
// ViSP
#include <visp3/core/vpColVector.h>
#include <visp3/core/vpList.h>
#include <visp3/core/vpMatrix.h>
 
// UsTK
#include <visp3/ustk_core/usConfig.h>
 
namespace usNeedleDetectionTools
{
double angle(double *p1, double *p2, double *q1, double *q2);
 
vpMatrix approximateCoordinates(vpMatrix X, vpMatrix MSS, unsigned int s);
 
void arithmeticMean(const vpMatrix points, double *mean, unsigned int npts, unsigned int d);
 
void computeQuantile(unsigned char *data, unsigned int dataSize, unsigned int nDesired, unsigned int &q,
                     unsigned int &nThresholded);
void computeQuantile(short *data, unsigned int dataSize, unsigned int nDesired, unsigned int &q,
                     unsigned int &nThresholded);
void computeQuantile(const unsigned int *data, unsigned int dataSize, unsigned int nDesired, unsigned int &q,
                     unsigned int &nThresholded);
 
void computeControlPoints(const vpMatrix &model, vpMatrix &mss, double **points, const unsigned int &nPoints);
 
double dist3(double *x, double *y);
 
bool findEntry(const vpMatrix &model, double *entry, unsigned int nPoints, const vpColVector &origin,
               const vpColVector &entryPlane, int VOI[6]);
 
vpColVector geometricMedian(const vpMatrix points, unsigned int npts, unsigned int d);
 
bool inside(const vpColVector point, int *VOI);
 
bool intersectionLinePlane(vpMatrix line, const vpColVector &plane, double offset, double &x, double &y, double &z);
 
void linearRegression(vpMatrix &points, double &x, double &y, double &z, double &u, double &v, double &w);
void linearRegression(vpMatrix &points, double &a, double &b, double &u, double &v, double &w);
 
short quantile(short *data, unsigned int num, unsigned int n);
 
double variance(vpMatrix points, unsigned int npts, unsigned int d);
}
 
#endif // __usNeedleDetectionTools_h_