libs/gil/doc/html/reference/harris_8hpp_source.html

   1 <!-- HTML header for doxygen 1.8.13-->
   2 <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "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/xhtml;charset=UTF-8"/>
   6     <meta http-equiv="X-UA-Compatible" content="IE=9"/>
   7     <meta name="generator" content="Doxygen 1.8.15"/>
   8     <meta name="viewport" content="width=device-width, initial-scale=1"/>
   9     <title>Generic Image Library: harris.hpp Source File</title>
  10     <link href="tabs.css" rel="stylesheet" type="text/css"/>
  11     <script type="text/javascript" src="jquery.js"></script>
  12     <script type="text/javascript" src="dynsections.js"></script>
  13     <link href="doxygen.css" rel="stylesheet" type="text/css" />
  14     <link href="doxygen-boost.css" rel="stylesheet" type="text/css"/>
  15   </head>
  16   <body>
  17     <div class="boost-header">
  18       <table border="0" cellpadding="7" cellspacing="0" width="100%" summary="header">
  19         <tr>
  20           <td valign="top" width="300">
  21             <h3><a href="../index.html"><img alt="Boost GIL" src="../_static/gil.png" border="0"></a></h3>
  22           </td>
  23           <td ><h1 align="center"><a href="../index.html"></a></h1></td>
  24           <td></td>
  25         </tr>
  26       </table>
  27     </div>
  28     <hr/>
  29     <div id="top"><!-- do not remove this div, it is closed by doxygen! -->
  30 <!-- Generated by Doxygen 1.8.15 -->
  31 <script type="text/javascript" src="menudata.js"></script>
  32 <script type="text/javascript" src="menu.js"></script>
  33 <script type="text/javascript">
  34 /* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
  35 $(function() {
  36   initMenu('',false,false,'search.php','Search');
  37 });
  38 /* @license-end */</script>
  39 <div id="main-nav"></div>
  40 <div id="nav-path" class="navpath">
  41   <ul>
  42 <li class="navelem"><a class="el" href="dir_d44c64559bbebec7f509842c48db8b23.html">include</a></li><li class="navelem"><a class="el" href="dir_1878a3f4746a95c6aad317458cc7ef80.html">boost</a></li><li class="navelem"><a class="el" href="dir_df4750f408086f9b9c1b5ee4251365ff.html">gil</a></li><li class="navelem"><a class="el" href="dir_29d7d1879aced8e7726ea9a11f3c7333.html">image_processing</a></li>  </ul>
  43 </div>
  44 </div><!-- top -->
  45 <div class="header">
  46   <div class="headertitle">
  47 <div class="title">harris.hpp</div>  </div>
  48 </div><!--header-->
  49 <div class="contents">
  50 <div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno">    1</span>&#160;<span class="comment">//</span></div><div class="line"><a name="l00002"></a><span class="lineno">    2</span>&#160;<span class="comment">// Copyright 2019 Olzhas Zhumabek &lt;anonymous.from.applecity@gmail.com&gt;</span></div><div class="line"><a name="l00003"></a><span class="lineno">    3</span>&#160;<span class="comment">//</span></div><div class="line"><a name="l00004"></a><span class="lineno">    4</span>&#160;<span class="comment">// Use, modification and distribution are subject to the Boost Software License,</span></div><div class="line"><a name="l00005"></a><span class="lineno">    5</span>&#160;<span class="comment">// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at</span></div><div class="line"><a name="l00006"></a><span class="lineno">    6</span>&#160;<span class="comment">// http://www.boost.org/LICENSE_1_0.txt)</span></div><div class="line"><a name="l00007"></a><span class="lineno">    7</span>&#160;<span class="comment">//</span></div><div class="line"><a name="l00008"></a><span class="lineno">    8</span>&#160;<span class="preprocessor">#ifndef BOOST_GIL_IMAGE_PROCESSING_HARRIS_HPP</span></div><div class="line"><a name="l00009"></a><span class="lineno">    9</span>&#160;<span class="preprocessor">#define BOOST_GIL_IMAGE_PROCESSING_HARRIS_HPP</span></div><div class="line"><a name="l00010"></a><span class="lineno">   10</span>&#160;</div><div class="line"><a name="l00011"></a><span class="lineno">   11</span>&#160;<span class="preprocessor">#include &lt;boost/gil/image_view.hpp&gt;</span></div><div class="line"><a name="l00012"></a><span class="lineno">   12</span>&#160;<span class="preprocessor">#include &lt;boost/gil/typedefs.hpp&gt;</span></div><div class="line"><a name="l00013"></a><span class="lineno">   13</span>&#160;<span class="preprocessor">#include &lt;boost/gil/extension/numeric/kernel.hpp&gt;</span></div><div class="line"><a name="l00014"></a><span class="lineno">   14</span>&#160;</div><div class="line"><a name="l00015"></a><span class="lineno">   15</span>&#160;<span class="keyword">namespace </span>boost { <span class="keyword">namespace </span>gil {</div><div class="line"><a name="l00021"></a><span class="lineno">   21</span>&#160;</div><div class="line"><a name="l00022"></a><span class="lineno">   22</span>&#160;</div><div class="line"><a name="l00033"></a><span class="lineno">   33</span>&#160;<span class="keyword">template</span> &lt;<span class="keyword">typename</span> T, <span class="keyword">typename</span> Allocator&gt;</div><div class="line"><a name="l00034"></a><span class="lineno"><a class="line" href="group___corner_detection_algorithms.html#gab73f324a539c60eae52a8017feab44d8">   34</a></span>&#160;<span class="keywordtype">void</span> <a class="code" href="group___corner_detection_algorithms.html#gab73f324a539c60eae52a8017feab44d8">compute_harris_responses</a>(</div><div class="line"><a name="l00035"></a><span class="lineno">   35</span>&#160;    boost::gil::gray32f_view_t m11,</div><div class="line"><a name="l00036"></a><span class="lineno">   36</span>&#160;    boost::gil::gray32f_view_t m12_21,</div><div class="line"><a name="l00037"></a><span class="lineno">   37</span>&#160;    boost::gil::gray32f_view_t m22,</div><div class="line"><a name="l00038"></a><span class="lineno">   38</span>&#160;    boost::gil::detail::kernel_2d&lt;T, Allocator&gt; weights,</div><div class="line"><a name="l00039"></a><span class="lineno">   39</span>&#160;    <span class="keywordtype">float</span> k,</div><div class="line"><a name="l00040"></a><span class="lineno">   40</span>&#160;    boost::gil::gray32f_view_t harris_response)</div><div class="line"><a name="l00041"></a><span class="lineno">   41</span>&#160;{</div><div class="line"><a name="l00042"></a><span class="lineno">   42</span>&#160;    <span class="keywordflow">if</span> (m11.dimensions() != m12_21.dimensions() || m12_21.dimensions() != m22.dimensions()) {</div><div class="line"><a name="l00043"></a><span class="lineno">   43</span>&#160;        <span class="keywordflow">throw</span> std::invalid_argument(<span class="stringliteral">&quot;m prefixed arguments must represent&quot;</span></div><div class="line"><a name="l00044"></a><span class="lineno">   44</span>&#160;            <span class="stringliteral">&quot; tensor from the same image&quot;</span>);</div><div class="line"><a name="l00045"></a><span class="lineno">   45</span>&#160;    }</div><div class="line"><a name="l00046"></a><span class="lineno">   46</span>&#160;</div><div class="line"><a name="l00047"></a><span class="lineno">   47</span>&#160;    <span class="keyword">auto</span> <span class="keyword">const</span> window_length = weights.size();</div><div class="line"><a name="l00048"></a><span class="lineno">   48</span>&#160;    <span class="keyword">auto</span> <span class="keyword">const</span> width = m11.width();</div><div class="line"><a name="l00049"></a><span class="lineno">   49</span>&#160;    <span class="keyword">auto</span> <span class="keyword">const</span> height = m11.height();</div><div class="line"><a name="l00050"></a><span class="lineno">   50</span>&#160;    <span class="keyword">auto</span> <span class="keyword">const</span> half_length = window_length / 2;</div><div class="line"><a name="l00051"></a><span class="lineno">   51</span>&#160;</div><div class="line"><a name="l00052"></a><span class="lineno">   52</span>&#160;    <span class="keywordflow">for</span> (<span class="keyword">auto</span> y = half_length; y &lt; height - half_length; ++y)</div><div class="line"><a name="l00053"></a><span class="lineno">   53</span>&#160;    {</div><div class="line"><a name="l00054"></a><span class="lineno">   54</span>&#160;        <span class="keywordflow">for</span> (<span class="keyword">auto</span> x = half_length; x &lt; width - half_length; ++x)</div><div class="line"><a name="l00055"></a><span class="lineno">   55</span>&#160;        {</div><div class="line"><a name="l00056"></a><span class="lineno">   56</span>&#160;            <span class="keywordtype">float</span> ddxx = 0;</div><div class="line"><a name="l00057"></a><span class="lineno">   57</span>&#160;            <span class="keywordtype">float</span> dxdy = 0;</div><div class="line"><a name="l00058"></a><span class="lineno">   58</span>&#160;            <span class="keywordtype">float</span> ddyy = 0;</div><div class="line"><a name="l00059"></a><span class="lineno">   59</span>&#160;            <span class="keywordflow">for</span> (gil::gray32f_view_t::coord_t y_kernel = 0;</div><div class="line"><a name="l00060"></a><span class="lineno">   60</span>&#160;                y_kernel &lt; window_length;</div><div class="line"><a name="l00061"></a><span class="lineno">   61</span>&#160;                ++y_kernel) {</div><div class="line"><a name="l00062"></a><span class="lineno">   62</span>&#160;                <span class="keywordflow">for</span> (gil::gray32f_view_t::coord_t x_kernel = 0;</div><div class="line"><a name="l00063"></a><span class="lineno">   63</span>&#160;                    x_kernel &lt; window_length;</div><div class="line"><a name="l00064"></a><span class="lineno">   64</span>&#160;                    ++x_kernel) {</div><div class="line"><a name="l00065"></a><span class="lineno">   65</span>&#160;                    ddxx += m11(x + x_kernel - half_length, y + y_kernel - half_length)</div><div class="line"><a name="l00066"></a><span class="lineno">   66</span>&#160;                        .at(std::integral_constant&lt;int, 0&gt;{}) * weights.at(x_kernel, y_kernel);</div><div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;                    dxdy += m12_21(x + x_kernel - half_length, y + y_kernel - half_length)</div><div class="line"><a name="l00068"></a><span class="lineno">   68</span>&#160;                        .at(std::integral_constant&lt;int, 0&gt;{}) * weights.at(x_kernel, y_kernel);</div><div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160;                    ddyy += m22(x + x_kernel - half_length, y + y_kernel - half_length)</div><div class="line"><a name="l00070"></a><span class="lineno">   70</span>&#160;                        .at(std::integral_constant&lt;int, 0&gt;{}) * weights.at(x_kernel, y_kernel);</div><div class="line"><a name="l00071"></a><span class="lineno">   71</span>&#160;                }</div><div class="line"><a name="l00072"></a><span class="lineno">   72</span>&#160;            }</div><div class="line"><a name="l00073"></a><span class="lineno">   73</span>&#160;            <span class="keyword">auto</span> det = (ddxx * ddyy) - dxdy * dxdy;</div><div class="line"><a name="l00074"></a><span class="lineno">   74</span>&#160;            <span class="keyword">auto</span> trace = ddxx + ddyy;</div><div class="line"><a name="l00075"></a><span class="lineno">   75</span>&#160;            <span class="keyword">auto</span> harris_value = det - k * trace * trace;</div><div class="line"><a name="l00076"></a><span class="lineno">   76</span>&#160;            harris_response(x, y).at(std::integral_constant&lt;int, 0&gt;{}) = harris_value;</div><div class="line"><a name="l00077"></a><span class="lineno">   77</span>&#160;        }</div><div class="line"><a name="l00078"></a><span class="lineno">   78</span>&#160;    }</div><div class="line"><a name="l00079"></a><span class="lineno">   79</span>&#160;}</div><div class="line"><a name="l00080"></a><span class="lineno">   80</span>&#160;</div><div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160;}} <span class="comment">//namespace boost::gil</span></div><div class="line"><a name="l00082"></a><span class="lineno">   82</span>&#160;<span class="preprocessor">#endif</span></div><div class="ttc" id="group___corner_detection_algorithms_html_gab73f324a539c60eae52a8017feab44d8"><div class="ttname"><a href="group___corner_detection_algorithms.html#gab73f324a539c60eae52a8017feab44d8">boost::gil::compute_harris_responses</a></div><div class="ttdeci">void compute_harris_responses(boost::gil::gray32f_view_t m11, boost::gil::gray32f_view_t m12_21, boost::gil::gray32f_view_t m22, boost::gil::detail::kernel_2d&lt; T, Allocator &gt; weights, float k, boost::gil::gray32f_view_t harris_response)</div><div class="ttdoc">function to record Harris responsesThis algorithm computes Harris responses for structure tensor repr...</div><div class="ttdef"><b>Definition:</b> harris.hpp:34</div></div>
  51 </div><!-- fragment --></div><!-- contents -->
  52      <!-- HTML footer for doxygen 1.8.13-->
  53      <!-- start footer part -->
  54      <hr class="footer"/>
  55      <address class="footer">
  56       <small>
  57         Generated by &#160;<a href="http://www.doxygen.org/index.html">doxygen</a> 1.8.15
  58       </small>
  59     </address>
  60   </body>
  61 </html>