Image Component Library (ICL)
|
00001 /******************************************************************** 00002 ** Image Component Library (ICL) ** 00003 ** ** 00004 ** Copyright (C) 2006-2013 CITEC, University of Bielefeld ** 00005 ** Neuroinformatics Group ** 00006 ** Website: www.iclcv.org and ** 00007 ** http://opensource.cit-ec.de/projects/icl ** 00008 ** ** 00009 ** File : ICLMath/src/ICLMath/LeastSquareModelFitting2D.h ** 00010 ** Module : ICLMath ** 00011 ** Authors: Christof Elbrechter ** 00012 ** ** 00013 ** ** 00014 ** GNU LESSER GENERAL PUBLIC LICENSE ** 00015 ** This file may be used under the terms of the GNU Lesser General ** 00016 ** Public License version 3.0 as published by the ** 00017 ** ** 00018 ** Free Software Foundation and appearing in the file LICENSE.GPL ** 00019 ** included in the packaging of this file. Please review the ** 00020 ** following information to ensure the license requirements will ** 00021 ** be met: http://www.gnu.org/licenses/lgpl-3.0.txt ** 00022 ** ** 00023 ** The development of this software was supported by the ** 00024 ** Excellence Cluster EXC 277 Cognitive Interaction Technology. ** 00025 ** The Excellence Cluster EXC 277 is a grant of the Deutsche ** 00026 ** Forschungsgemeinschaft (DFG) in the context of the German ** 00027 ** Excellence Initiative. ** 00028 ** ** 00029 ********************************************************************/ 00030 00031 #pragma once 00032 00033 #include <ICLMath/LeastSquareModelFitting.h> 00034 00035 namespace icl{ 00036 namespace math{ 00037 00039 00041 class LeastSquareModelFitting2D : public LeastSquareModelFitting<double,Point32f>{ 00043 typedef LeastSquareModelFitting<double,Point32f> Super; 00044 public: 00046 LeastSquareModelFitting2D(){} 00047 00049 LeastSquareModelFitting2D(int modelDim, DesignMatrixGen gen, 00050 DynMatrix<double> *constraintMatrix = 0): 00051 Super(modelDim,gen,constraintMatrix){} 00052 00054 00055 static inline void line_gen(const Point32f &p, double *d){ 00056 d[0] = p.x; 00057 d[1] = p.y; 00058 d[2] = 1; 00059 } 00060 00062 00063 static inline void circle_gen(const Point32f &p, double *d){ 00064 d[0] = sqr(p.x) + sqr(p.y); 00065 d[1] = p.x; 00066 d[2] = p.y; 00067 d[3] = 1; 00068 } 00069 00071 00072 static inline void restr_ellipse_gen(const Point32f &p, double *d){ 00073 d[0] = sqr(p.x); 00074 d[1] = sqr(p.y); 00075 d[2] = p.x; 00076 d[3] = p.y; 00077 d[4] = 1; 00078 } 00079 00081 00082 static inline void ellipse_gen(const Point32f &p, double *d){ 00083 d[0] = sqr(p.x); 00084 d[1] = p.x * p.y; 00085 d[2] = sqr(p.y); 00086 d[3] = p.x; 00087 d[4] = p.y; 00088 d[5] = 1; 00089 } 00090 00091 inline std::vector<double> fit(const std::vector<Point32f> &points){ 00092 return Super::fit(points); 00093 } 00094 inline icl64f getError(const std::vector<double> &model,const Point32f &p) { 00095 return Super::getError(model,p); 00096 } 00097 }; 00098 00099 00100 } // namespace math 00101 } 00102