SIFT+RANSAC识别物体，并标注边界

#include “sift.h“
#include “my_function.h“
#include 
#include 
#include “math.h“

//#include
using namespace std;



class CvModelEstimator2
{
public:
    CvModelEstimator2（int _modelPoints CvSize _modelSize int _maxBasicSolutions）;
    virtual ~CvModelEstimator2（）;

    virtual int runKernel（ const CvMat* m1 const CvMat* m2 CvMat* model ）=0;
    virtual bool runLMeDS（ const CvMat* m1 const CvMat* m2 CvMat* model
                           CvMat* mask double confidence=0.99 int maxIters=2000 ）;
    virtual bool runRANSAC（ const CvMat* m1 const CvMat* m2 CvMat* model
                            CvMat* mask double threshold
                            double confidence=0.99 int maxIters=2000 ）;
    virtual bool refine（ const CvMat* const CvMat* CvMat* int ） { return true; }
    virtual void setSeed（ int64 seed ）;

protected:
    virtual void computeReprojError（ const CvMat* m1 const CvMat* m2
                                     const CvMat* model CvMat* error ） = 0;
    virtual int findInliers（ const CvMat* m1 const CvMat* m2
                             const CvMat* model CvMat* error
                             CvMat* mask double threshold ）;
    virtual bool getSubset（ const CvMat* m1 const CvMat* m2
                            CvMat* ms1 CvMat* ms2 int maxAttempts=1000 ）;
    virtual bool checkSubset（ const CvMat* ms1 int count ）;

    CvRNG rng;
    int modelPoints;
    CvSize modelSize;
    int maxBasicSolutions;
    bool checkPartialSubsets;
};


class CvFMEstimator : public CvModelEstimator2
{
public:
    CvFMEstimator（ int _modelPoints ）;

    virtual int runKernel（ const CvMat* m1 const CvMat* m2 CvMat* model ）;
    virtual int run7Point（ const CvMat* m1 const CvMat* m2 CvMat* model ）;
    virtual int run8Point（ const CvMat* m1 const CvMat* m2 CvMat* model ）;
protected:
    virtual void computeReprojError（ const CvMat* m1 const CvMat* m2
                                     const CvMat* model CvMat* error ）;
};



CV_IMPL void cvConvertPointsHomogeneous（ const CvMat* src CvMat* dst ）
{
    Ptr temp denom;

    int i s_count s_dims d_count d_dims;
    CvMat _src _dst _ones;
    CvMat* ones = 0;

    if（ !CV_IS_MAT（src） ）
        CV_Error（ !src ? CV_StsNullPtr : CV_StsBadArg
        “The input parameter is not a valid matrix“ ）;

    if（ !CV_IS_MAT（dst） ）
        CV_Error（ !dst ? CV_StsNullPtr : CV_StsBadArg
        “The output parameter is not a valid matrix“ ）;

    if（ src == dst || src->data.ptr == dst->data.ptr ）
    {
        if（ src != dst && （!CV_ARE_TYPES_EQ（src dst） || !CV_ARE_SIZES_EQ（srcdst）） ）
            CV_Error（ CV_StsBadArg “Invalid inplace operation“ ）;
        return;
    }

    if（ src->rows > src->cols ）
    {
        if（ !（（src->cols > 1） ^ （CV_MAT_CN（src->type） > 1）） ）
            CV_Error（ CV_StsBadSize “Either the number of channels or columns