boost/gil
2
0
Fork 0
mirror of https://github.com/boostorg/gil.git synced 2026-01-25 06:12:22 +00:00
Code Issues Projects Releases Wiki Activity
Files
43e8aef5bd0bd032f628401be21fbdf232e82b58
gil/develop/doc/html/reference/hessian_8hpp_source.html
github-actions[bot] 3b73200ee1 deploy: 49d9276706
2025-04-04 16:28:54 +00:00

139 lines
14 KiB
HTML
Raw Blame History

<!-- HTML header for doxygen 1.8.13-->
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.9.8"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>Generic Image Library: hessian.hpp Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
<link href="doxygen-boost.css" rel="stylesheet" type="text/css"/>
</head>
<body>
<div class="boost-header">
<table border="0" cellpadding="7" cellspacing="0" width="100%" summary="header">
<tr>
<td valign="top" width="300">
<h3><a href="../index.html"><img alt="Boost GIL" src="../_static/gil.png" border="0"></a></h3>
</td>
<td ><h1 align="center"><a href="../index.html"></a></h1></td>
<td></td>
</tr>
</table>
</div>
<hr/>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<!-- Generated by Doxygen 1.9.8 -->
<script type="text/javascript" src="menudata.js"></script>
<script type="text/javascript" src="menu.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&amp;dn=expat.txt MIT */
$(function() {
initMenu('',false,false,'search.php','Search');
});
/* @license-end */
</script>
<div id="main-nav"></div>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&amp;dn=expat.txt MIT */
$(document).ready(function() { init_codefold(0); });
/* @license-end */
</script>
<div id="nav-path" class="navpath">
<ul>
<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>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle"><div class="title">hessian.hpp</div></div>
</div><!--header-->
<div class="contents">
<div class="fragment"><div class="line"><a id="l00001" name="l00001"></a><span class="lineno"> 1</span><span class="comment">//</span></div>
<div class="line"><a id="l00002" name="l00002"></a><span class="lineno"> 2</span><span class="comment">// Copyright 2019 Olzhas Zhumabek &lt;anonymous.from.applecity@gmail.com&gt;</span></div>
<div class="line"><a id="l00003" name="l00003"></a><span class="lineno"> 3</span><span class="comment">// Copyright 2021 Scramjet911 &lt;36035352+Scramjet911@users.noreply.github.com&gt;</span></div>
<div class="line"><a id="l00004" name="l00004"></a><span class="lineno"> 4</span><span class="comment">// Use, modification and distribution are subject to the Boost Software License,</span></div>
<div class="line"><a id="l00005" name="l00005"></a><span class="lineno"> 5</span><span class="comment">// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at</span></div>
<div class="line"><a id="l00006" name="l00006"></a><span class="lineno"> 6</span><span class="comment">// http://www.boost.org/LICENSE_1_0.txt)</span></div>
<div class="line"><a id="l00007" name="l00007"></a><span class="lineno"> 7</span> </div>
<div class="line"><a id="l00008" name="l00008"></a><span class="lineno"> 8</span><span class="preprocessor">#ifndef BOOST_GIL_IMAGE_PROCESSING_HESSIAN_HPP</span></div>
<div class="line"><a id="l00009" name="l00009"></a><span class="lineno"> 9</span><span class="preprocessor">#define BOOST_GIL_IMAGE_PROCESSING_HESSIAN_HPP</span></div>
<div class="line"><a id="l00010" name="l00010"></a><span class="lineno"> 10</span> </div>
<div class="line"><a id="l00011" name="l00011"></a><span class="lineno"> 11</span><span class="preprocessor">#include &lt;boost/gil/image_view.hpp&gt;</span></div>
<div class="line"><a id="l00012" name="l00012"></a><span class="lineno"> 12</span><span class="preprocessor">#include &lt;boost/gil/typedefs.hpp&gt;</span></div>
<div class="line"><a id="l00013" name="l00013"></a><span class="lineno"> 13</span><span class="preprocessor">#include &lt;boost/gil/image_processing/kernel.hpp&gt;</span></div>
<div class="line"><a id="l00014" name="l00014"></a><span class="lineno"> 14</span><span class="preprocessor">#include &lt;stdexcept&gt;</span></div>
<div class="line"><a id="l00015" name="l00015"></a><span class="lineno"> 15</span> </div>
<div class="line"><a id="l00016" name="l00016"></a><span class="lineno"> 16</span><span class="keyword">namespace </span><a class="code hl_namespace" href="namespaceboost.html">boost</a> { <span class="keyword">namespace </span>gil {</div>
<div class="line"><a id="l00017" name="l00017"></a><span class="lineno"> 17</span> </div>
<div class="line"><a id="l00026" name="l00026"></a><span class="lineno"> 26</span><span class="keyword">template</span> &lt;<span class="keyword">typename</span> GradientView, <span class="keyword">typename</span> T, <span class="keyword">typename</span> Allocator, <span class="keyword">typename</span> OutputView&gt;</div>
<div class="line"><a id="l00027" name="l00027"></a><span class="lineno"> 27</span><span class="keyword">inline</span> <span class="keywordtype">void</span> compute_hessian_responses(</div>
<div class="line"><a id="l00028" name="l00028"></a><span class="lineno"> 28</span> GradientView ddxx,</div>
<div class="line"><a id="l00029" name="l00029"></a><span class="lineno"> 29</span> GradientView dxdy,</div>
<div class="line"><a id="l00030" name="l00030"></a><span class="lineno"> 30</span> GradientView ddyy,</div>
<div class="line"><a id="l00031" name="l00031"></a><span class="lineno"> 31</span> <span class="keyword">const</span> detail::kernel_2d&lt;T, Allocator&gt;&amp; weights,</div>
<div class="line"><a id="l00032" name="l00032"></a><span class="lineno"> 32</span> OutputView dst)</div>
<div class="line"><a id="l00033" name="l00033"></a><span class="lineno"> 33</span>{</div>
<div class="line"><a id="l00034" name="l00034"></a><span class="lineno"> 34</span> <span class="keywordflow">if</span> (ddxx.dimensions() != ddyy.dimensions()</div>
<div class="line"><a id="l00035" name="l00035"></a><span class="lineno"> 35</span> || ddyy.dimensions() != dxdy.dimensions()</div>
<div class="line"><a id="l00036" name="l00036"></a><span class="lineno"> 36</span> || dxdy.dimensions() != dst.dimensions()</div>
<div class="line"><a id="l00037" name="l00037"></a><span class="lineno"> 37</span> || weights.center_x() != weights.center_y())</div>
<div class="line"><a id="l00038" name="l00038"></a><span class="lineno"> 38</span> {</div>
<div class="line"><a id="l00039" name="l00039"></a><span class="lineno"> 39</span> <span class="keywordflow">throw</span> std::invalid_argument(<span class="stringliteral">&quot;dimensions of views are not the same&quot;</span></div>
<div class="line"><a id="l00040" name="l00040"></a><span class="lineno"> 40</span> <span class="stringliteral">&quot; or weights don&#39;t have equal width and height&quot;</span></div>
<div class="line"><a id="l00041" name="l00041"></a><span class="lineno"> 41</span> <span class="stringliteral">&quot; or weights&#39; dimensions are not odd&quot;</span>);</div>
<div class="line"><a id="l00042" name="l00042"></a><span class="lineno"> 42</span> }</div>
<div class="line"><a id="l00043" name="l00043"></a><span class="lineno"> 43</span> <span class="comment">// Use pixel type of output, as values will be written to output</span></div>
<div class="line"><a id="l00044" name="l00044"></a><span class="lineno"> 44</span> <span class="keyword">using </span>pixel_t = <span class="keyword">typename</span> std::remove_reference&lt;decltype(std::declval&lt;OutputView&gt;()(0, 0))&gt;::type;</div>
<div class="line"><a id="l00045" name="l00045"></a><span class="lineno"> 45</span> </div>
<div class="line"><a id="l00046" name="l00046"></a><span class="lineno"> 46</span> <span class="keyword">using </span>channel_t = <span class="keyword">typename</span> std::remove_reference</div>
<div class="line"><a id="l00047" name="l00047"></a><span class="lineno"> 47</span> &lt;</div>
<div class="line"><a id="l00048" name="l00048"></a><span class="lineno"> 48</span> <span class="keyword">decltype</span>(std::declval&lt;pixel_t&gt;().at(std::integral_constant&lt;int, 0&gt;{}))</div>
<div class="line"><a id="l00049" name="l00049"></a><span class="lineno"> 49</span> &gt;::type;</div>
<div class="line"><a id="l00050" name="l00050"></a><span class="lineno"> 50</span> </div>
<div class="line"><a id="l00051" name="l00051"></a><span class="lineno"> 51</span> </div>
<div class="line"><a id="l00052" name="l00052"></a><span class="lineno"> 52</span> <span class="keyword">auto</span> center = weights.center_y();</div>
<div class="line"><a id="l00053" name="l00053"></a><span class="lineno"> 53</span> <span class="keywordflow">for</span> (<span class="keyword">auto</span> y = center; y &lt; dst.height() - center; ++y)</div>
<div class="line"><a id="l00054" name="l00054"></a><span class="lineno"> 54</span> {</div>
<div class="line"><a id="l00055" name="l00055"></a><span class="lineno"> 55</span> <span class="keywordflow">for</span> (<span class="keyword">auto</span> x = center; x &lt; dst.width() - center; ++x)</div>
<div class="line"><a id="l00056" name="l00056"></a><span class="lineno"> 56</span> {</div>
<div class="line"><a id="l00057" name="l00057"></a><span class="lineno"> 57</span> <span class="keyword">auto</span> ddxx_i = channel_t();</div>
<div class="line"><a id="l00058" name="l00058"></a><span class="lineno"> 58</span> <span class="keyword">auto</span> ddyy_i = channel_t();</div>
<div class="line"><a id="l00059" name="l00059"></a><span class="lineno"> 59</span> <span class="keyword">auto</span> dxdy_i = channel_t();</div>
<div class="line"><a id="l00060" name="l00060"></a><span class="lineno"> 60</span> <span class="keywordflow">for</span> (<span class="keyword">typename</span> OutputView::coord_t w_y = 0; w_y &lt; static_cast&lt;std::ptrdiff_t&gt;(weights.size()); ++w_y)</div>
<div class="line"><a id="l00061" name="l00061"></a><span class="lineno"> 61</span> {</div>
<div class="line"><a id="l00062" name="l00062"></a><span class="lineno"> 62</span> <span class="keywordflow">for</span> (<span class="keyword">typename</span> OutputView::coord_t w_x = 0; w_x &lt; static_cast&lt;std::ptrdiff_t&gt;(weights.size()); ++w_x)</div>
<div class="line"><a id="l00063" name="l00063"></a><span class="lineno"> 63</span> {</div>
<div class="line"><a id="l00064" name="l00064"></a><span class="lineno"> 64</span> ddxx_i += ddxx(x + w_x - center, y + w_y - center)</div>
<div class="line"><a id="l00065" name="l00065"></a><span class="lineno"> 65</span> .at(std::integral_constant&lt;int, 0&gt;{}) * weights.at(w_x, w_y);</div>
<div class="line"><a id="l00066" name="l00066"></a><span class="lineno"> 66</span> ddyy_i += ddyy(x + w_x - center, y + w_y - center)</div>
<div class="line"><a id="l00067" name="l00067"></a><span class="lineno"> 67</span> .at(std::integral_constant&lt;int, 0&gt;{}) * weights.at(w_x, w_y);</div>
<div class="line"><a id="l00068" name="l00068"></a><span class="lineno"> 68</span> dxdy_i += dxdy(x + w_x - center, y + w_y - center)</div>
<div class="line"><a id="l00069" name="l00069"></a><span class="lineno"> 69</span> .at(std::integral_constant&lt;int, 0&gt;{}) * weights.at(w_x, w_y);</div>
<div class="line"><a id="l00070" name="l00070"></a><span class="lineno"> 70</span> }</div>
<div class="line"><a id="l00071" name="l00071"></a><span class="lineno"> 71</span> }</div>
<div class="line"><a id="l00072" name="l00072"></a><span class="lineno"> 72</span> <span class="keyword">auto</span> determinant = ddxx_i * ddyy_i - dxdy_i * dxdy_i;</div>
<div class="line"><a id="l00073" name="l00073"></a><span class="lineno"> 73</span> dst(x, y).at(std::integral_constant&lt;int, 0&gt;{}) = determinant;</div>
<div class="line"><a id="l00074" name="l00074"></a><span class="lineno"> 74</span> }</div>
<div class="line"><a id="l00075" name="l00075"></a><span class="lineno"> 75</span> }</div>
<div class="line"><a id="l00076" name="l00076"></a><span class="lineno"> 76</span>}</div>
<div class="line"><a id="l00077" name="l00077"></a><span class="lineno"> 77</span> </div>
<div class="line"><a id="l00078" name="l00078"></a><span class="lineno"> 78</span>}} <span class="comment">// namespace boost::gil</span></div>
<div class="line"><a id="l00079" name="l00079"></a><span class="lineno"> 79</span> </div>
<div class="line"><a id="l00080" name="l00080"></a><span class="lineno"> 80</span><span class="preprocessor">#endif</span></div>
<div class="ttc" id="anamespaceboost_html"><div class="ttname"><a href="namespaceboost.html">boost</a></div><div class="ttdoc">defined(BOOST_NO_CXX17_HDR_MEMORY_RESOURCE)</div><div class="ttdef"><b>Definition</b> algorithm.hpp:36</div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- HTML footer for doxygen 1.8.13-->
<!-- start footer part -->
<hr class="footer"/>
<address class="footer">
<small>
Generated by &#160;<a href="http://www.doxygen.org/index.html">doxygen</a> 1.9.8
</small>
</address>
</body>
</html>
Reference in New Issue View Git Blame Copy Permalink