boost/gil
2
0
Fork 0
mirror of https://github.com/boostorg/gil.git synced 2026-01-24 18:02:17 +00:00
Code Issues Projects Releases Wiki Activity
Files
7f4861e20093f92dbdf7fabcfb6873c63b676124
gil/develop/doc/html/reference/numeric_8hpp_source.html
Travis CI 7b3321edbb Deploy to GitHub Pages: 5c52dff369
2019-08-21 06:56:40 +00:00

74 lines
30 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.8.11"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>Generic Image Library: numeric.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.8.11 -->
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main&#160;Page</span></a></li>
<li><a href="pages.html"><span>Related&#160;Pages</span></a></li>
<li><a href="modules.html"><span>Modules</span></a></li>
<li><a href="annotated.html"><span>Classes</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
</ul>
</div>
<div id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File&#160;List</span></a></li>
</ul>
</div>
<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">numeric.hpp</div> </div>
</div><!--header-->
<div class="contents">
<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_NUMERIC_HPP</span></div><div class="line"><a name="l00009"></a><span class="lineno"> 9</span>&#160;<span class="preprocessor">#define BOOST_GIL_IMAGE_PROCESSING_NUMERIC_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/detail/math.hpp&gt;</span></div><div class="line"><a name="l00012"></a><span class="lineno"> 12</span>&#160;<span class="preprocessor">#include &lt;cmath&gt;</span></div><div class="line"><a name="l00013"></a><span class="lineno"> 13</span>&#160;<span class="preprocessor">#include &lt;boost/gil/image_view.hpp&gt;</span></div><div class="line"><a name="l00014"></a><span class="lineno"> 14</span>&#160;<span class="preprocessor">#include &lt;boost/gil/typedefs.hpp&gt;</span></div><div class="line"><a name="l00015"></a><span class="lineno"> 15</span>&#160;</div><div class="line"><a name="l00016"></a><span class="lineno"> 16</span>&#160;<span class="keyword">namespace </span><a class="code" href="namespaceboost.html">boost</a> { <span class="keyword">namespace </span>gil {</div><div class="line"><a name="l00017"></a><span class="lineno"> 17</span>&#160;</div><div class="line"><a name="l00029"></a><span class="lineno"> 29</span>&#160;<span class="keyword">inline</span> <span class="keywordtype">double</span> normalized_sinc(<span class="keywordtype">double</span> x)</div><div class="line"><a name="l00030"></a><span class="lineno"> 30</span>&#160;{</div><div class="line"><a name="l00031"></a><span class="lineno"> 31</span>&#160; <span class="keywordflow">return</span> std::sin(x * boost::gil::pi) / (x * boost::gil::pi);</div><div class="line"><a name="l00032"></a><span class="lineno"> 32</span>&#160;}</div><div class="line"><a name="l00033"></a><span class="lineno"> 33</span>&#160;</div><div class="line"><a name="l00041"></a><span class="lineno"><a class="line" href="group___image_processing_math.html#gaaa4583e6d012a8b49ad1738cec762651"> 41</a></span>&#160;<span class="keyword">inline</span> <span class="keywordtype">double</span> <a class="code" href="group___image_processing_math.html#gaaa4583e6d012a8b49ad1738cec762651">lanczos</a>(<span class="keywordtype">double</span> x, std::ptrdiff_t a)</div><div class="line"><a name="l00042"></a><span class="lineno"> 42</span>&#160;{</div><div class="line"><a name="l00043"></a><span class="lineno"> 43</span>&#160; <span class="comment">// means == but &lt;= avoids compiler warning</span></div><div class="line"><a name="l00044"></a><span class="lineno"> 44</span>&#160; <span class="keywordflow">if</span> (0 &lt;= x &amp;&amp; x &lt;= 0)</div><div class="line"><a name="l00045"></a><span class="lineno"> 45</span>&#160; <span class="keywordflow">return</span> 1;</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="keywordflow">if</span> (-a &lt; x &amp;&amp; x &lt; a)</div><div class="line"><a name="l00048"></a><span class="lineno"> 48</span>&#160; <span class="keywordflow">return</span> normalized_sinc(x) / normalized_sinc(x / static_cast&lt;double&gt;(a));</div><div class="line"><a name="l00049"></a><span class="lineno"> 49</span>&#160;</div><div class="line"><a name="l00050"></a><span class="lineno"> 50</span>&#160; <span class="keywordflow">return</span> 0;</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;</div><div class="line"><a name="l00053"></a><span class="lineno"> 53</span>&#160;<span class="keyword">inline</span> <span class="keywordtype">void</span> compute_tensor_entries(</div><div class="line"><a name="l00054"></a><span class="lineno"> 54</span>&#160; boost::gil::gray16_view_t dx,</div><div class="line"><a name="l00055"></a><span class="lineno"> 55</span>&#160; boost::gil::gray16_view_t dy,</div><div class="line"><a name="l00056"></a><span class="lineno"> 56</span>&#160; boost::gil::gray32f_view_t m11,</div><div class="line"><a name="l00057"></a><span class="lineno"> 57</span>&#160; boost::gil::gray32f_view_t m12_21,</div><div class="line"><a name="l00058"></a><span class="lineno"> 58</span>&#160; boost::gil::gray32f_view_t m22)</div><div class="line"><a name="l00059"></a><span class="lineno"> 59</span>&#160;{</div><div class="line"><a name="l00060"></a><span class="lineno"> 60</span>&#160; <span class="keywordflow">for</span> (std::ptrdiff_t y = 0; y &lt; dx.height(); ++y) {</div><div class="line"><a name="l00061"></a><span class="lineno"> 61</span>&#160; <span class="keywordflow">for</span> (std::ptrdiff_t x = 0; x &lt; dx.width(); ++x) {</div><div class="line"><a name="l00062"></a><span class="lineno"> 62</span>&#160; <span class="keyword">auto</span> dx_value = dx(x, y);</div><div class="line"><a name="l00063"></a><span class="lineno"> 63</span>&#160; <span class="keyword">auto</span> dy_value = dy(x, y);</div><div class="line"><a name="l00064"></a><span class="lineno"> 64</span>&#160; m11(x, y) = dx_value * dx_value;</div><div class="line"><a name="l00065"></a><span class="lineno"> 65</span>&#160; m12_21(x, y) = dx_value * dy_value;</div><div class="line"><a name="l00066"></a><span class="lineno"> 66</span>&#160; m22(x, y) = dy_value * dy_value;</div><div class="line"><a name="l00067"></a><span class="lineno"> 67</span>&#160; }</div><div class="line"><a name="l00068"></a><span class="lineno"> 68</span>&#160; }</div><div class="line"><a name="l00069"></a><span class="lineno"> 69</span>&#160;}</div><div class="line"><a name="l00070"></a><span class="lineno"> 70</span>&#160;</div><div class="line"><a name="l00080"></a><span class="lineno"> 80</span>&#160;<span class="keyword">template</span> &lt;<span class="keyword">typename</span> GradientView, <span class="keyword">typename</span> OutputView&gt;</div><div class="line"><a name="l00081"></a><span class="lineno"><a class="line" href="group___image_processing_math.html#ga5b662d207dd54aba778c0f8c3697c025"> 81</a></span>&#160;<span class="keyword">inline</span> <span class="keywordtype">void</span> <a class="code" href="group___image_processing_math.html#ga5b662d207dd54aba778c0f8c3697c025">compute_hessian_entries</a>(</div><div class="line"><a name="l00082"></a><span class="lineno"> 82</span>&#160; GradientView dx,</div><div class="line"><a name="l00083"></a><span class="lineno"> 83</span>&#160; GradientView dy,</div><div class="line"><a name="l00084"></a><span class="lineno"> 84</span>&#160; OutputView ddxx,</div><div class="line"><a name="l00085"></a><span class="lineno"> 85</span>&#160; OutputView dxdy,</div><div class="line"><a name="l00086"></a><span class="lineno"> 86</span>&#160; OutputView ddyy)</div><div class="line"><a name="l00087"></a><span class="lineno"> 87</span>&#160;{</div><div class="line"><a name="l00088"></a><span class="lineno"> 88</span>&#160; <span class="keyword">using</span> x_coord_t = <span class="keyword">typename</span> OutputView::x_coord_t;</div><div class="line"><a name="l00089"></a><span class="lineno"> 89</span>&#160; <span class="keyword">using</span> y_coord_t = <span class="keyword">typename</span> OutputView::y_coord_t;</div><div class="line"><a name="l00090"></a><span class="lineno"> 90</span>&#160; <span class="keyword">using</span> pixel_t = <span class="keyword">typename</span> std::remove_reference&lt;decltype(ddxx(0, 0))&gt;::type;</div><div class="line"><a name="l00091"></a><span class="lineno"> 91</span>&#160; <span class="keyword">using</span> channel_t = <span class="keyword">typename</span> std::remove_reference&lt;</div><div class="line"><a name="l00092"></a><span class="lineno"> 92</span>&#160; decltype(</div><div class="line"><a name="l00093"></a><span class="lineno"> 93</span>&#160; std::declval&lt;pixel_t&gt;().at(</div><div class="line"><a name="l00094"></a><span class="lineno"> 94</span>&#160; std::integral_constant&lt;int, 0&gt;{}</div><div class="line"><a name="l00095"></a><span class="lineno"> 95</span>&#160; )</div><div class="line"><a name="l00096"></a><span class="lineno"> 96</span>&#160; )</div><div class="line"><a name="l00097"></a><span class="lineno"> 97</span>&#160; &gt;::type;</div><div class="line"><a name="l00098"></a><span class="lineno"> 98</span>&#160;</div><div class="line"><a name="l00099"></a><span class="lineno"> 99</span>&#160; constexpr <span class="keywordtype">double</span> x_kernel[3][3] =</div><div class="line"><a name="l00100"></a><span class="lineno"> 100</span>&#160; {</div><div class="line"><a name="l00101"></a><span class="lineno"> 101</span>&#160; {1, 0, -1},</div><div class="line"><a name="l00102"></a><span class="lineno"> 102</span>&#160; {2, 0, -2},</div><div class="line"><a name="l00103"></a><span class="lineno"> 103</span>&#160; {1, 0, -1}</div><div class="line"><a name="l00104"></a><span class="lineno"> 104</span>&#160; };</div><div class="line"><a name="l00105"></a><span class="lineno"> 105</span>&#160; constexpr <span class="keywordtype">double</span> y_kernel[3][3] =</div><div class="line"><a name="l00106"></a><span class="lineno"> 106</span>&#160; {</div><div class="line"><a name="l00107"></a><span class="lineno"> 107</span>&#160; {1, 2, 1},</div><div class="line"><a name="l00108"></a><span class="lineno"> 108</span>&#160; {0, 0, 0},</div><div class="line"><a name="l00109"></a><span class="lineno"> 109</span>&#160; {-1, -2, -1}</div><div class="line"><a name="l00110"></a><span class="lineno"> 110</span>&#160; };</div><div class="line"><a name="l00111"></a><span class="lineno"> 111</span>&#160; constexpr <span class="keyword">auto</span> chosen_channel = std::integral_constant&lt;int, 0&gt;{};</div><div class="line"><a name="l00112"></a><span class="lineno"> 112</span>&#160; <span class="keywordflow">for</span> (y_coord_t y = 1; y &lt; ddxx.height() - 1; ++y)</div><div class="line"><a name="l00113"></a><span class="lineno"> 113</span>&#160; {</div><div class="line"><a name="l00114"></a><span class="lineno"> 114</span>&#160; <span class="keywordflow">for</span> (x_coord_t x = 1; x &lt; ddxx.width() - 1; ++x)</div><div class="line"><a name="l00115"></a><span class="lineno"> 115</span>&#160; {</div><div class="line"><a name="l00116"></a><span class="lineno"> 116</span>&#160; pixel_t ddxx_i;</div><div class="line"><a name="l00117"></a><span class="lineno"> 117</span>&#160; static_transform(ddxx_i, ddxx_i,</div><div class="line"><a name="l00118"></a><span class="lineno"> 118</span>&#160; [](channel_t) { <span class="keywordflow">return</span> <span class="keyword">static_cast&lt;</span>channel_t<span class="keyword">&gt;</span>(0); });</div><div class="line"><a name="l00119"></a><span class="lineno"> 119</span>&#160; pixel_t dxdy_i;</div><div class="line"><a name="l00120"></a><span class="lineno"> 120</span>&#160; static_transform(dxdy_i, dxdy_i,</div><div class="line"><a name="l00121"></a><span class="lineno"> 121</span>&#160; [](channel_t) { <span class="keywordflow">return</span> <span class="keyword">static_cast&lt;</span>channel_t<span class="keyword">&gt;</span>(0); });</div><div class="line"><a name="l00122"></a><span class="lineno"> 122</span>&#160; pixel_t ddyy_i;</div><div class="line"><a name="l00123"></a><span class="lineno"> 123</span>&#160; static_transform(ddyy_i, ddyy_i,</div><div class="line"><a name="l00124"></a><span class="lineno"> 124</span>&#160; [](channel_t) { <span class="keywordflow">return</span> <span class="keyword">static_cast&lt;</span>channel_t<span class="keyword">&gt;</span>(0); });</div><div class="line"><a name="l00125"></a><span class="lineno"> 125</span>&#160; <span class="keywordflow">for</span> (y_coord_t y_filter = 0; y_filter &lt; 2; ++y_filter)</div><div class="line"><a name="l00126"></a><span class="lineno"> 126</span>&#160; {</div><div class="line"><a name="l00127"></a><span class="lineno"> 127</span>&#160; <span class="keywordflow">for</span> (x_coord_t x_filter = 0; x_filter &lt; 2; ++x_filter)</div><div class="line"><a name="l00128"></a><span class="lineno"> 128</span>&#160; {</div><div class="line"><a name="l00129"></a><span class="lineno"> 129</span>&#160; <span class="keyword">auto</span> adjusted_y = y + y_filter - 1;</div><div class="line"><a name="l00130"></a><span class="lineno"> 130</span>&#160; <span class="keyword">auto</span> adjusted_x = x + x_filter - 1;</div><div class="line"><a name="l00131"></a><span class="lineno"> 131</span>&#160; ddxx_i.at(std::integral_constant&lt;int, 0&gt;{}) +=</div><div class="line"><a name="l00132"></a><span class="lineno"> 132</span>&#160; dx(adjusted_x, adjusted_y).at(chosen_channel)</div><div class="line"><a name="l00133"></a><span class="lineno"> 133</span>&#160; * x_kernel[y_filter][x_filter];</div><div class="line"><a name="l00134"></a><span class="lineno"> 134</span>&#160; dxdy_i.at(std::integral_constant&lt;int, 0&gt;{}) +=</div><div class="line"><a name="l00135"></a><span class="lineno"> 135</span>&#160; dx(adjusted_x, adjusted_y).at(chosen_channel)</div><div class="line"><a name="l00136"></a><span class="lineno"> 136</span>&#160; * y_kernel[y_filter][x_filter];</div><div class="line"><a name="l00137"></a><span class="lineno"> 137</span>&#160; ddyy_i.at(std::integral_constant&lt;int, 0&gt;{}) +=</div><div class="line"><a name="l00138"></a><span class="lineno"> 138</span>&#160; dy(adjusted_x, adjusted_y).at(chosen_channel)</div><div class="line"><a name="l00139"></a><span class="lineno"> 139</span>&#160; * y_kernel[y_filter][x_filter];</div><div class="line"><a name="l00140"></a><span class="lineno"> 140</span>&#160; }</div><div class="line"><a name="l00141"></a><span class="lineno"> 141</span>&#160; }</div><div class="line"><a name="l00142"></a><span class="lineno"> 142</span>&#160; ddxx(x, y) = ddxx_i;</div><div class="line"><a name="l00143"></a><span class="lineno"> 143</span>&#160; dxdy(x, y) = dxdy_i;</div><div class="line"><a name="l00144"></a><span class="lineno"> 144</span>&#160; ddyy(x, y) = ddyy_i;</div><div class="line"><a name="l00145"></a><span class="lineno"> 145</span>&#160; }</div><div class="line"><a name="l00146"></a><span class="lineno"> 146</span>&#160; }</div><div class="line"><a name="l00147"></a><span class="lineno"> 147</span>&#160;}</div><div class="line"><a name="l00148"></a><span class="lineno"> 148</span>&#160;</div><div class="line"><a name="l00155"></a><span class="lineno"><a class="line" href="group___image_processing_math.html#ga598481b60bf8ad39a0dce57d65d1820a"> 155</a></span>&#160;<span class="keyword">inline</span> <span class="keywordtype">void</span> <a class="code" href="group___image_processing_math.html#ga598481b60bf8ad39a0dce57d65d1820a">generate_normalized_mean</a>(boost::gil::gray32f_view_t dst)</div><div class="line"><a name="l00156"></a><span class="lineno"> 156</span>&#160;{</div><div class="line"><a name="l00157"></a><span class="lineno"> 157</span>&#160; <span class="keywordflow">if</span> (dst.width() != dst.height() || dst.width() % 2 != 1)</div><div class="line"><a name="l00158"></a><span class="lineno"> 158</span>&#160; <span class="keywordflow">throw</span> std::invalid_argument(<span class="stringliteral">&quot;kernel dimensions should be odd and equal&quot;</span>);</div><div class="line"><a name="l00159"></a><span class="lineno"> 159</span>&#160; <span class="keyword">const</span> <span class="keywordtype">float</span> entry = 1.0f / <span class="keyword">static_cast&lt;</span><span class="keywordtype">float</span><span class="keyword">&gt;</span>(dst.size());</div><div class="line"><a name="l00160"></a><span class="lineno"> 160</span>&#160;</div><div class="line"><a name="l00161"></a><span class="lineno"> 161</span>&#160; <span class="keywordflow">for</span> (<span class="keyword">auto</span>&amp; <a class="code" href="structboost_1_1gil_1_1pixel.html">pixel</a>: dst) {</div><div class="line"><a name="l00162"></a><span class="lineno"> 162</span>&#160; <a class="code" href="structboost_1_1gil_1_1pixel.html">pixel</a>.at(std::integral_constant&lt;int, 0&gt;{}) = entry;</div><div class="line"><a name="l00163"></a><span class="lineno"> 163</span>&#160; }</div><div class="line"><a name="l00164"></a><span class="lineno"> 164</span>&#160;}</div><div class="line"><a name="l00165"></a><span class="lineno"> 165</span>&#160;</div><div class="line"><a name="l00170"></a><span class="lineno"><a class="line" href="group___image_processing_math.html#gaf09a84c1960713b6a8df3e1d78b1b733"> 170</a></span>&#160;<span class="keyword">inline</span> <span class="keywordtype">void</span> <a class="code" href="group___image_processing_math.html#gaf09a84c1960713b6a8df3e1d78b1b733">generate_unnormalized_mean</a>(boost::gil::gray32f_view_t dst)</div><div class="line"><a name="l00171"></a><span class="lineno"> 171</span>&#160;{</div><div class="line"><a name="l00172"></a><span class="lineno"> 172</span>&#160; <span class="keywordflow">if</span> (dst.width() != dst.height() || dst.width() % 2 != 1)</div><div class="line"><a name="l00173"></a><span class="lineno"> 173</span>&#160; <span class="keywordflow">throw</span> std::invalid_argument(<span class="stringliteral">&quot;kernel dimensions should be odd and equal&quot;</span>);</div><div class="line"><a name="l00174"></a><span class="lineno"> 174</span>&#160;</div><div class="line"><a name="l00175"></a><span class="lineno"> 175</span>&#160; <span class="keywordflow">for</span> (<span class="keyword">auto</span>&amp; <a class="code" href="structboost_1_1gil_1_1pixel.html">pixel</a>: dst) {</div><div class="line"><a name="l00176"></a><span class="lineno"> 176</span>&#160; <a class="code" href="structboost_1_1gil_1_1pixel.html">pixel</a>.at(std::integral_constant&lt;int, 0&gt;{}) = 1.0f;</div><div class="line"><a name="l00177"></a><span class="lineno"> 177</span>&#160; }</div><div class="line"><a name="l00178"></a><span class="lineno"> 178</span>&#160;}</div><div class="line"><a name="l00179"></a><span class="lineno"> 179</span>&#160;</div><div class="line"><a name="l00185"></a><span class="lineno"><a class="line" href="group___image_processing_math.html#gae1bd55f00abfa4bad83560e35992f7f4"> 185</a></span>&#160;<span class="keyword">inline</span> <span class="keywordtype">void</span> <a class="code" href="group___image_processing_math.html#gae1bd55f00abfa4bad83560e35992f7f4">generate_gaussian_kernel</a>(boost::gil::gray32f_view_t dst, <span class="keywordtype">double</span> sigma)</div><div class="line"><a name="l00186"></a><span class="lineno"> 186</span>&#160;{</div><div class="line"><a name="l00187"></a><span class="lineno"> 187</span>&#160; <span class="keywordflow">if</span> (dst.width() != dst.height() || dst.width() % 2 != 1)</div><div class="line"><a name="l00188"></a><span class="lineno"> 188</span>&#160; <span class="keywordflow">throw</span> std::invalid_argument(<span class="stringliteral">&quot;kernel dimensions should be odd and equal&quot;</span>);</div><div class="line"><a name="l00189"></a><span class="lineno"> 189</span>&#160;</div><div class="line"><a name="l00190"></a><span class="lineno"> 190</span>&#160; <span class="keyword">const</span> <span class="keywordtype">double</span> denominator = 2 * boost::gil::pi * sigma * sigma;</div><div class="line"><a name="l00191"></a><span class="lineno"> 191</span>&#160; <span class="keyword">const</span> <span class="keyword">auto</span> middle = <a class="code" href="classboost_1_1gil_1_1point.html">boost::gil::point_t</a>(dst.width() / 2, dst.height() / 2);</div><div class="line"><a name="l00192"></a><span class="lineno"> 192</span>&#160; <span class="keywordflow">for</span> (boost::gil::gray32f_view_t::coord_t y = 0; y &lt; dst.height(); ++y)</div><div class="line"><a name="l00193"></a><span class="lineno"> 193</span>&#160; {</div><div class="line"><a name="l00194"></a><span class="lineno"> 194</span>&#160; <span class="keywordflow">for</span> (boost::gil::gray32f_view_t::coord_t x = 0; x &lt; dst.width(); ++x)</div><div class="line"><a name="l00195"></a><span class="lineno"> 195</span>&#160; {</div><div class="line"><a name="l00196"></a><span class="lineno"> 196</span>&#160; <span class="keyword">const</span> <span class="keyword">auto</span> delta_x = std::abs(middle.x - x);</div><div class="line"><a name="l00197"></a><span class="lineno"> 197</span>&#160; <span class="keyword">const</span> <span class="keyword">auto</span> delta_y = std::abs(middle.y - y);</div><div class="line"><a name="l00198"></a><span class="lineno"> 198</span>&#160; <span class="keyword">const</span> <span class="keywordtype">double</span> power = (delta_x * delta_x + delta_y * delta_y) / (2 * sigma * sigma);</div><div class="line"><a name="l00199"></a><span class="lineno"> 199</span>&#160; <span class="keyword">const</span> <span class="keywordtype">double</span> nominator = std::exp(-power);</div><div class="line"><a name="l00200"></a><span class="lineno"> 200</span>&#160; <span class="keyword">const</span> <span class="keywordtype">float</span> value = nominator / denominator;</div><div class="line"><a name="l00201"></a><span class="lineno"> 201</span>&#160; dst(x, y).at(std::integral_constant&lt;int, 0&gt;{}) = value;</div><div class="line"><a name="l00202"></a><span class="lineno"> 202</span>&#160; }</div><div class="line"><a name="l00203"></a><span class="lineno"> 203</span>&#160; }</div><div class="line"><a name="l00204"></a><span class="lineno"> 204</span>&#160;}</div><div class="line"><a name="l00205"></a><span class="lineno"> 205</span>&#160;</div><div class="line"><a name="l00206"></a><span class="lineno"> 206</span>&#160;}} <span class="comment">// namespace boost::gil</span></div><div class="line"><a name="l00207"></a><span class="lineno"> 207</span>&#160;</div><div class="line"><a name="l00208"></a><span class="lineno"> 208</span>&#160;<span class="preprocessor">#endif</span></div><div class="ttc" id="namespaceboost_html"><div class="ttname"><a href="namespaceboost.html">boost</a></div><div class="ttdef"><b>Definition:</b> algorithm.hpp:30</div></div>
<div class="ttc" id="structboost_1_1gil_1_1pixel_html"><div class="ttname"><a href="structboost_1_1gil_1_1pixel.html">boost::gil::pixel</a></div><div class="ttdoc">Represents a pixel value (a container of channels). Models: HomogeneousColorBaseValueConcept, PixelValueConcept, HomogeneousPixelBasedConcept. </div><div class="ttdef"><b>Definition:</b> metafunctions.hpp:23</div></div>
<div class="ttc" id="group___image_processing_math_html_ga5b662d207dd54aba778c0f8c3697c025"><div class="ttname"><a href="group___image_processing_math.html#ga5b662d207dd54aba778c0f8c3697c025">boost::gil::compute_hessian_entries</a></div><div class="ttdeci">void compute_hessian_entries(GradientView dx, GradientView dy, OutputView ddxx, OutputView dxdy, OutputView ddyy)</div><div class="ttdoc">Compute xy gradient, and second order x and y gradientsHessian matrix is defined as a matrix of parti...</div><div class="ttdef"><b>Definition:</b> numeric.hpp:81</div></div>
<div class="ttc" id="group___image_processing_math_html_gaaa4583e6d012a8b49ad1738cec762651"><div class="ttname"><a href="group___image_processing_math.html#gaaa4583e6d012a8b49ad1738cec762651">boost::gil::lanczos</a></div><div class="ttdeci">double lanczos(double x, std::ptrdiff_t a)</div><div class="ttdoc">Lanczos response at point xLanczos response is defined as: x == 0: 1 -a &lt; x &amp;&amp; x &lt; a: 0 otherwise: no...</div><div class="ttdef"><b>Definition:</b> numeric.hpp:41</div></div>
<div class="ttc" id="group___image_processing_math_html_ga598481b60bf8ad39a0dce57d65d1820a"><div class="ttname"><a href="group___image_processing_math.html#ga598481b60bf8ad39a0dce57d65d1820a">boost::gil::generate_normalized_mean</a></div><div class="ttdeci">void generate_normalized_mean(boost::gil::gray32f_view_t dst)</div><div class="ttdoc">Generate mean kernelFills supplied view with normalized mean in which all entries will be equal to...</div><div class="ttdef"><b>Definition:</b> numeric.hpp:155</div></div>
<div class="ttc" id="group___image_processing_math_html_gaf09a84c1960713b6a8df3e1d78b1b733"><div class="ttname"><a href="group___image_processing_math.html#gaf09a84c1960713b6a8df3e1d78b1b733">boost::gil::generate_unnormalized_mean</a></div><div class="ttdeci">void generate_unnormalized_mean(boost::gil::gray32f_view_t dst)</div><div class="ttdoc">Generate kernel with all 1sFills supplied view with 1s (ones) </div><div class="ttdef"><b>Definition:</b> numeric.hpp:170</div></div>
<div class="ttc" id="classboost_1_1gil_1_1point_html"><div class="ttname"><a href="classboost_1_1gil_1_1point.html">boost::gil::point&lt; std::ptrdiff_t &gt;</a></div></div>
<div class="ttc" id="group___image_processing_math_html_gae1bd55f00abfa4bad83560e35992f7f4"><div class="ttname"><a href="group___image_processing_math.html#gae1bd55f00abfa4bad83560e35992f7f4">boost::gil::generate_gaussian_kernel</a></div><div class="ttdeci">void generate_gaussian_kernel(boost::gil::gray32f_view_t dst, double sigma)</div><div class="ttdoc">Generate Gaussian kernelFills supplied view with values taken from Gaussian distribution. See https://en.wikipedia.org/wiki/Gaussian_blur. </div><div class="ttdef"><b>Definition:</b> numeric.hpp:185</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.8.11
</small>
</address>
</body>
</html>
Reference in New Issue View Git Blame Copy Permalink