fast_marching.hpp 4.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121
  1. /*M///////////////////////////////////////////////////////////////////////////////////////
  2. //
  3. // IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
  4. //
  5. // By downloading, copying, installing or using the software you agree to this license.
  6. // If you do not agree to this license, do not download, install,
  7. // copy or use the software.
  8. //
  9. //
  10. // License Agreement
  11. // For Open Source Computer Vision Library
  12. //
  13. // Copyright (C) 2000-2008, Intel Corporation, all rights reserved.
  14. // Copyright (C) 2009-2011, Willow Garage Inc., all rights reserved.
  15. // Third party copyrights are property of their respective owners.
  16. //
  17. // Redistribution and use in source and binary forms, with or without modification,
  18. // are permitted provided that the following conditions are met:
  19. //
  20. // * Redistribution's of source code must retain the above copyright notice,
  21. // this list of conditions and the following disclaimer.
  22. //
  23. // * Redistribution's in binary form must reproduce the above copyright notice,
  24. // this list of conditions and the following disclaimer in the documentation
  25. // and/or other materials provided with the distribution.
  26. //
  27. // * The name of the copyright holders may not be used to endorse or promote products
  28. // derived from this software without specific prior written permission.
  29. //
  30. // This software is provided by the copyright holders and contributors "as is" and
  31. // any express or implied warranties, including, but not limited to, the implied
  32. // warranties of merchantability and fitness for a particular purpose are disclaimed.
  33. // In no event shall the Intel Corporation or contributors be liable for any direct,
  34. // indirect, incidental, special, exemplary, or consequential damages
  35. // (including, but not limited to, procurement of substitute goods or services;
  36. // loss of use, data, or profits; or business interruption) however caused
  37. // and on any theory of liability, whether in contract, strict liability,
  38. // or tort (including negligence or otherwise) arising in any way out of
  39. // the use of this software, even if advised of the possibility of such damage.
  40. //
  41. //M*/
  42. #ifndef OPENCV_VIDEOSTAB_FAST_MARCHING_HPP
  43. #define OPENCV_VIDEOSTAB_FAST_MARCHING_HPP
  44. #include <cmath>
  45. #include <queue>
  46. #include <algorithm>
  47. #include "opencv2/core.hpp"
  48. namespace cv
  49. {
  50. namespace videostab
  51. {
  52. //! @addtogroup videostab_marching
  53. //! @{
  54. /** @brief Describes the Fast Marching Method implementation.
  55. See http://iwi.eldoc.ub.rug.nl/FILES/root/2004/JGraphToolsTelea/2004JGraphToolsTelea.pdf
  56. */
  57. class CV_EXPORTS FastMarchingMethod
  58. {
  59. public:
  60. FastMarchingMethod() : inf_(1e6f), size_(0) {}
  61. /** @brief Template method that runs the Fast Marching Method.
  62. @param mask Image mask. 0 value indicates that the pixel value must be inpainted, 255 indicates
  63. that the pixel value is known, other values aren't acceptable.
  64. @param inpaint Inpainting functor that overloads void operator ()(int x, int y).
  65. @return Inpainting functor.
  66. */
  67. template <typename Inpaint>
  68. Inpaint run(const Mat &mask, Inpaint inpaint);
  69. /**
  70. @return Distance map that's created during working of the method.
  71. */
  72. Mat distanceMap() const { return dist_; }
  73. private:
  74. enum { INSIDE = 0, BAND = 1, KNOWN = 255 };
  75. struct DXY
  76. {
  77. float dist;
  78. int x, y;
  79. DXY() : dist(0), x(0), y(0) {}
  80. DXY(float _dist, int _x, int _y) : dist(_dist), x(_x), y(_y) {}
  81. bool operator <(const DXY &dxy) const { return dist < dxy.dist; }
  82. };
  83. float solve(int x1, int y1, int x2, int y2) const;
  84. int& indexOf(const DXY &dxy) { return index_(dxy.y, dxy.x); }
  85. void heapUp(int idx);
  86. void heapDown(int idx);
  87. void heapAdd(const DXY &dxy);
  88. void heapRemoveMin();
  89. float inf_;
  90. cv::Mat_<uchar> flag_; // flag map
  91. cv::Mat_<float> dist_; // distance map
  92. cv::Mat_<int> index_; // index of point in the narrow band
  93. std::vector<DXY> narrowBand_; // narrow band heap
  94. int size_; // narrow band size
  95. };
  96. //! @}
  97. } // namespace videostab
  98. } // namespace cv
  99. #include "fast_marching_inl.hpp"
  100. #endif