histogram.html 13 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177
  1. <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
  2. "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
  3. <html xmlns="http://www.w3.org/1999/xhtml">
  4. <head>
  5. <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
  6. <title>Tutorial: Histogram - Boost.GIL documentation</title>
  7. <link rel="stylesheet" href="../_static/pygments.css" type="text/css" />
  8. <link rel="stylesheet" href="../_static/style.css" type="text/css" />
  9. <script type="text/javascript">
  10. var DOCUMENTATION_OPTIONS = {
  11. URL_ROOT: '../',
  12. VERSION: '',
  13. COLLAPSE_MODINDEX: false,
  14. FILE_SUFFIX: '.html'
  15. };
  16. </script>
  17. <script type="text/javascript" src="../_static/jquery.js"></script>
  18. <script type="text/javascript" src="../_static/underscore.js"></script>
  19. <script type="text/javascript" src="../_static/doctools.js"></script>
  20. <link rel="index" title="Index" href="../genindex.html" />
  21. <link rel="search" title="Search" href="../search.html" />
  22. <link rel="top" title="Boost.GIL documentation" href="../index.html" />
  23. <link rel="next" title="Tutorial: Image Gradient" href="gradient.html" />
  24. <link rel="prev" title="Tutorial: Video Lecture" href="video.html" />
  25. </head>
  26. <body>
  27. <div class="header">
  28. <table border="0" cellpadding="7" cellspacing="0" width="100%" summary=
  29. "header">
  30. <tr>
  31. <td valign="top" width="300">
  32. <h3><a href="../index.html"><img
  33. alt="C++ Boost" src="../_static/gil.png" border="0"></a></h3>
  34. </td>
  35. <td >
  36. <h1 align="center"><a href="../index.html"></a></h1>
  37. </td>
  38. <td>
  39. <div id="searchbox" style="display: none">
  40. <form class="search" action="../search.html" method="get">
  41. <input type="text" name="q" size="18" />
  42. <input type="submit" value="Search" />
  43. <input type="hidden" name="check_keywords" value="yes" />
  44. <input type="hidden" name="area" value="default" />
  45. </form>
  46. </div>
  47. <script type="text/javascript">$('#searchbox').show(0);</script>
  48. </td>
  49. </tr>
  50. </table>
  51. </div>
  52. <hr/>
  53. <div class="content">
  54. <div class="navbar" style="text-align:right;">
  55. <a class="prev" title="Tutorial: Video Lecture" href="video.html"><img src="../_static/prev.png" alt="prev"/></a>
  56. <a class="next" title="Tutorial: Image Gradient" href="gradient.html"><img src="../_static/next.png" alt="next"/></a>
  57. </div>
  58. <div class="section" id="tutorial-histogram">
  59. <h1>Tutorial: Histogram</h1>
  60. <div class="contents local topic" id="contents">
  61. <ul class="simple">
  62. <li><a class="reference internal" href="#original-implementation" id="id1">Original implementation</a></li>
  63. <li><a class="reference internal" href="#gil-implementation" id="id2">GIL implementation</a></li>
  64. </ul>
  65. </div>
  66. <p>This is a short tutorial presenting an example of a very simple sample of code
  67. from an existing code base that calculates histogram of an image.
  68. Next, the program is rewritten using GIL featres.</p>
  69. <div class="section" id="original-implementation">
  70. <h2><a class="toc-backref" href="#id1">Original implementation</a></h2>
  71. <p>Actual code from a commercial software product that computes the luminosity
  72. histogram (variable names have been changed and unrelated parts removed):</p>
  73. <div class="highlight-cpp"><div class="highlight"><pre><span class="kt">void</span> <span class="nf">luminosity_hist</span><span class="p">(</span>
  74. <span class="n">std</span><span class="o">::</span><span class="kt">uint8_t</span> <span class="k">const</span><span class="o">*</span> <span class="n">r</span><span class="p">,</span> <span class="n">std</span><span class="o">::</span><span class="kt">uint8_t</span> <span class="k">const</span><span class="o">*</span> <span class="n">g</span><span class="p">,</span> <span class="n">std</span><span class="o">::</span><span class="kt">uint8_t</span> <span class="k">const</span><span class="o">*</span> <span class="n">b</span><span class="p">,</span>
  75. <span class="kt">int</span> <span class="n">rows</span><span class="p">,</span> <span class="kt">int</span> <span class="n">cols</span><span class="p">,</span> <span class="kt">int</span> <span class="n">sRowBytes</span><span class="p">,</span> <span class="n">Histogram</span><span class="o">*</span> <span class="n">hist</span><span class="p">)</span>
  76. <span class="p">{</span>
  77. <span class="k">for</span> <span class="p">(</span><span class="kt">int</span> <span class="n">r</span> <span class="o">=</span> <span class="mi">0</span><span class="p">;</span> <span class="n">r</span> <span class="o">&lt;</span> <span class="n">rows</span><span class="p">;</span> <span class="n">r</span><span class="o">++</span><span class="p">)</span>
  78. <span class="p">{</span>
  79. <span class="k">for</span> <span class="p">(</span><span class="kt">int</span> <span class="n">c</span> <span class="o">=</span> <span class="mi">0</span><span class="p">;</span> <span class="n">c</span> <span class="o">&lt;</span> <span class="n">cols</span><span class="p">;</span> <span class="n">c</span><span class="o">++</span><span class="p">)</span>
  80. <span class="p">{</span>
  81. <span class="kt">int</span> <span class="n">v</span> <span class="o">=</span> <span class="n">RGBToGray</span><span class="p">(</span><span class="n">r</span><span class="p">[</span><span class="n">c</span><span class="p">],</span> <span class="n">g</span><span class="p">[</span><span class="n">c</span><span class="p">],</span> <span class="n">b</span><span class="p">[</span><span class="n">c</span><span class="p">]);</span> <span class="c1">// call internal function or macro</span>
  82. <span class="p">(</span><span class="o">*</span><span class="n">hist</span><span class="p">)[</span><span class="n">v</span><span class="p">]</span><span class="o">++</span><span class="p">;</span>
  83. <span class="p">}</span>
  84. <span class="n">r</span> <span class="o">+=</span> <span class="n">sRowBytes</span><span class="p">;</span>
  85. <span class="n">g</span> <span class="o">+=</span> <span class="n">sRowBytes</span><span class="p">;</span>
  86. <span class="n">b</span> <span class="o">+=</span> <span class="n">sRowBytes</span><span class="p">;</span>
  87. <span class="p">}</span>
  88. <span class="p">}</span>
  89. </pre></div>
  90. </div>
  91. <p>Let&#8217;s consider the following issues of the implementation above:</p>
  92. <ul class="simple">
  93. <li>Works only for RGB (duplicate versions exist for other color spaces)</li>
  94. <li>Works only for 8-bit images (duplicate versions exist)</li>
  95. <li>Works only for planar images</li>
  96. </ul>
  97. </div>
  98. <div class="section" id="gil-implementation">
  99. <h2><a class="toc-backref" href="#id2">GIL implementation</a></h2>
  100. <div class="highlight-cpp"><div class="highlight"><pre><span class="k">template</span> <span class="o">&lt;</span><span class="k">typename</span> <span class="n">GrayView</span><span class="p">,</span> <span class="k">typename</span> <span class="n">R</span><span class="o">&gt;</span>
  101. <span class="kt">void</span> <span class="n">grayimage_histogram</span><span class="p">(</span><span class="n">GrayView</span><span class="o">&amp;</span> <span class="n">img</span><span class="p">,</span> <span class="n">R</span><span class="o">&amp;</span> <span class="n">hist</span><span class="p">)</span>
  102. <span class="p">{</span>
  103. <span class="k">for</span> <span class="p">(</span><span class="k">typename</span> <span class="n">GrayView</span><span class="o">::</span><span class="n">iterator</span> <span class="n">it</span><span class="o">=</span><span class="n">img</span><span class="p">.</span><span class="n">begin</span><span class="p">();</span> <span class="n">it</span><span class="o">!=</span><span class="n">img</span><span class="p">.</span><span class="n">end</span><span class="p">();</span> <span class="o">++</span><span class="n">it</span><span class="p">)</span>
  104. <span class="o">++</span><span class="n">hist</span><span class="p">[</span><span class="o">*</span><span class="n">it</span><span class="p">];</span>
  105. <span class="p">}</span>
  106. <span class="k">template</span> <span class="o">&lt;</span><span class="k">typename</span> <span class="n">View</span><span class="p">,</span> <span class="k">typename</span> <span class="n">R</span><span class="o">&gt;</span>
  107. <span class="kt">void</span> <span class="n">luminosity8bit_hist</span><span class="p">(</span><span class="n">View</span><span class="o">&amp;</span> <span class="n">img</span><span class="p">,</span> <span class="n">R</span><span class="o">&amp;</span> <span class="n">hist</span><span class="p">)</span>
  108. <span class="p">{</span>
  109. <span class="n">grayimage_histogram</span><span class="p">(</span><span class="n">color_converted_view</span><span class="o">&lt;</span><span class="n">gray8_pixel_t</span><span class="o">&gt;</span><span class="p">(</span><span class="n">img</span><span class="p">),</span><span class="n">hist</span><span class="p">);</span>
  110. <span class="p">}</span>
  111. </pre></div>
  112. </div>
  113. <p>Using the Boost.Lambda library (or C++11 lambda) features it can written
  114. even simpler:</p>
  115. <div class="highlight-cpp"><div class="highlight"><pre><span class="k">using</span> <span class="n">boost</span><span class="o">::</span><span class="n">lambda</span><span class="p">;</span>
  116. <span class="k">template</span> <span class="o">&lt;</span><span class="k">typename</span> <span class="n">GrayView</span><span class="p">,</span> <span class="k">typename</span> <span class="n">R</span><span class="o">&gt;</span>
  117. <span class="kt">void</span> <span class="n">grayimage_histogram</span><span class="p">(</span><span class="n">GrayView</span><span class="o">&amp;</span> <span class="n">img</span><span class="p">,</span> <span class="n">R</span><span class="o">&amp;</span> <span class="n">hist</span><span class="p">)</span>
  118. <span class="p">{</span>
  119. <span class="n">for_each_pixel</span><span class="p">(</span><span class="n">img</span><span class="p">,</span> <span class="o">++</span><span class="n">var</span><span class="p">(</span><span class="n">hist</span><span class="p">)[</span><span class="n">_1</span><span class="p">]);</span>
  120. <span class="p">}</span>
  121. </pre></div>
  122. </div>
  123. <p>Let&#8217;s consider the following advantages of the GIL version:</p>
  124. <ul class="simple">
  125. <li>Works with any supported channel depth, color space, channel ordering
  126. (RGB vs BGR), and row alignment policy.</li>
  127. <li>Works for both planar and interleaved images.</li>
  128. <li>Works with new color spaces, channel depths and image types that can be
  129. provided in future extensions of GIL</li>
  130. <li>The second version is as efficient as the hand-coded version</li>
  131. </ul>
  132. <p>Shortly, it is also very flexible.</p>
  133. <p>For example, to compute the histogram of the second channel of the top left
  134. quadrant of the image, taking every other row and column, call:</p>
  135. <div class="highlight-cpp"><div class="highlight"><pre><span class="n">grayimage_histogram</span><span class="p">(</span>
  136. <span class="n">nth_channel_view</span><span class="p">(</span>
  137. <span class="n">subsampled_view</span><span class="p">(</span>
  138. <span class="n">subimage_view</span><span class="p">(</span><span class="n">img</span><span class="p">,</span>
  139. <span class="mi">0</span><span class="p">,</span><span class="mi">0</span><span class="p">,</span> <span class="n">img</span><span class="p">.</span><span class="n">width</span><span class="p">()</span> <span class="o">/</span> <span class="mi">2</span><span class="p">,</span> <span class="n">img</span><span class="p">.</span><span class="n">height</span><span class="p">()</span> <span class="o">/</span> <span class="mi">2</span><span class="p">),</span> <span class="c1">// upper left quadrant</span>
  140. <span class="mi">2</span><span class="p">,</span> <span class="mi">2</span> <span class="c1">// skip every other row and column</span>
  141. <span class="p">),</span>
  142. <span class="mi">1</span> <span class="c1">// index of the second channel (for example, green for RGB)</span>
  143. <span class="p">),</span>
  144. <span class="n">hist</span>
  145. <span class="p">);</span>
  146. </pre></div>
  147. </div>
  148. <p>Since GIL operates on the source pixels of <code class="docutils literal"><span class="pre">img</span></code> object directly, no extra
  149. memory is allocated and no images are copied.</p>
  150. </div>
  151. </div>
  152. <div class="navbar" style="text-align:right;">
  153. <a class="prev" title="Tutorial: Video Lecture" href="video.html"><img src="../_static/prev.png" alt="prev"/></a>
  154. <a class="next" title="Tutorial: Image Gradient" href="gradient.html"><img src="../_static/next.png" alt="next"/></a>
  155. </div>
  156. </div>
  157. <div class="footer" role="contentinfo">
  158. Last updated on 2019-12-10 00:12:10.
  159. Created using <a href="http://sphinx-doc.org/">Sphinx</a> 1.5.6.
  160. </div>
  161. </body>
  162. </html>