开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 图形图像处理 > 查看源码
KeyPoint.cpp
资源名称：OPENCV_SIFT_VC6.rar　[点击查看]
上传用户：yli86818
上传日期：2014-07-15
资源大小：273k
文件大小：2k
源码类别：

图形图像处理

开发平台：

Visual C++

KeyPoint.cpp：源码内容
  1. // KeyPoint.cpp: implementation of the CKeyPoint class.
  2. //
  3. //////////////////////////////////////////////////////////////////////
  5. #include "stdafx.h"
  6. #include "TestSIFT.h"
  7. #include "KeyPoint.h"
  9. #ifdef _DEBUG
  10. #undef THIS_FILE
  11. static char THIS_FILE[]=__FILE__;
  12. #define new DEBUG_NEW
  13. #endif
  15. //////////////////////////////////////////////////////////////////////
  16. // Construction/Destruction
  17. //////////////////////////////////////////////////////////////////////
  19. CKeyPoint:: CKeyPoint (IplImage* gaussImg, double x, double y, double imgScale,
  20. double kpScale, double orientation)
  21. {
  22. m_gaussImg = gaussImg;
  23. m_x = x;
  24. m_y = y;
  25. m_imgScale = imgScale;
  26. m_kpScale = kpScale;
  27. m_orientation = orientation;
  28. m_featureVec = NULL;
  29. }
  31. CKeyPoint::~CKeyPoint()
  32. {
  33. if (m_featureVec != NULL)
  34. delete m_featureVec;
  35. }
  37. void CKeyPoint::CreateVector (int xDim, int yDim, int oDim)
  38. {
  39. m_xDim = xDim;
  40. m_yDim = yDim;
  41. m_oDim = oDim;
  42. m_featureVec = new double[yDim * xDim * oDim];
  43. memset (m_featureVec,0,sizeof(double)*yDim*xDim*oDim);
  44. }
  47. void CKeyPoint::FVSet(int xI, int yI, int oI, double value)
  48. {
  49. // row major
  50. m_featureVec[(yI * m_xDim * m_oDim) + (xI * m_oDim) + oI] = value;
  51. }
  53. double CKeyPoint::FVGet(int xI, int yI, int oI)
  54. {
  55. // row major
  56. return m_featureVec[(yI * m_xDim * m_oDim) + (xI * m_oDim) + oI];
  57. }
  59. // Threshhold and normalize feature vector.
  60. // Note that the feature vector as a whole is normalized (Lowe's paper is
  61. // a bit unclear at that point).
  62. void CKeyPoint::CapAndNormalizeFV (double fvGradHicap)
  63. {
  64. // Straight normalization
  65. int len = m_xDim * m_yDim * m_oDim;
  66. double norm = 0.0;
  67. for (int i = 0 ; i < len ; i++)
  68. norm += m_featureVec[i]*m_featureVec[i];
  70. norm = sqrt (norm);
  71. if (norm == 0.0)
  72. return;
  74. for (i = 0 ; i < len; i++)
  75. m_featureVec[i] /= norm;
  77. // Hicap after normalization
  78. for (i = 0 ; i < len; i++) 
  79. if (m_featureVec[i] > fvGradHicap) 
  80. m_featureVec[i] = fvGradHicap;
  82. // Renormalize again
  83. norm = 0;
  84. for (i = 0 ; i < len ; i++)
  85. norm += m_featureVec[i]*m_featureVec[i];
  86. norm = sqrt (norm);
  88. for (i = 0 ; i < len; i++)
  89. m_featureVec[i] /= norm;
  90. }