mirror of
https://github.com/boostorg/gil.git
synced 2026-01-23 05:32:12 +00:00
175 lines
19 KiB
HTML
175 lines
19 KiB
HTML
<!-- 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: histogram_equalization.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&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&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">histogram_equalization.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 2020 Debabrata Mandal <mandaldebabrata123@gmail.com></span></div>
|
|
<div class="line"><a id="l00003" name="l00003"></a><span class="lineno"> 3</span><span class="comment">//</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><span class="comment">//</span></div>
|
|
<div class="line"><a id="l00008" name="l00008"></a><span class="lineno"> 8</span><span class="preprocessor">#ifndef BOOST_GIL_IMAGE_PROCESSING_HISTOGRAM_EQUALIZATION_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_HISTOGRAM_EQUALIZATION_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 <boost/gil/histogram.hpp></span></div>
|
|
<div class="line"><a id="l00012" name="l00012"></a><span class="lineno"> 12</span><span class="preprocessor">#include <boost/gil/image.hpp></span></div>
|
|
<div class="line"><a id="l00013" name="l00013"></a><span class="lineno"> 13</span> </div>
|
|
<div class="line"><a id="l00014" name="l00014"></a><span class="lineno"> 14</span><span class="preprocessor">#include <cmath></span></div>
|
|
<div class="line"><a id="l00015" name="l00015"></a><span class="lineno"> 15</span><span class="preprocessor">#include <map></span></div>
|
|
<div class="line"><a id="l00016" name="l00016"></a><span class="lineno"> 16</span><span class="preprocessor">#include <vector></span></div>
|
|
<div class="line"><a id="l00017" name="l00017"></a><span class="lineno"> 17</span> </div>
|
|
<div class="line"><a id="l00018" name="l00018"></a><span class="lineno"> 18</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="l00019" name="l00019"></a><span class="lineno"> 19</span> </div>
|
|
<div class="line"><a id="l00035" name="l00035"></a><span class="lineno"> 35</span> </div>
|
|
<div class="line"><a id="l00043" name="l00043"></a><span class="lineno"> 43</span><span class="keyword">template</span> <<span class="keyword">typename</span> SrcKeyType></div>
|
|
<div class="line"><a id="l00044" name="l00044"></a><span class="lineno"> 44</span><span class="keyword">auto</span> histogram_equalization(histogram<SrcKeyType> <span class="keyword">const</span>& src_hist)</div>
|
|
<div class="line"><a id="l00045" name="l00045"></a><span class="lineno"> 45</span> -> std::map<SrcKeyType, SrcKeyType></div>
|
|
<div class="line"><a id="l00046" name="l00046"></a><span class="lineno"> 46</span>{</div>
|
|
<div class="line"><a id="l00047" name="l00047"></a><span class="lineno"> 47</span> histogram<SrcKeyType> dst_hist;</div>
|
|
<div class="line"><a id="l00048" name="l00048"></a><span class="lineno"> 48</span> <span class="keywordflow">return</span> histogram_equalization(src_hist, dst_hist);</div>
|
|
<div class="line"><a id="l00049" name="l00049"></a><span class="lineno"> 49</span>}</div>
|
|
<div class="line"><a id="l00050" name="l00050"></a><span class="lineno"> 50</span> </div>
|
|
<div class="line"><a id="l00061" name="l00061"></a><span class="lineno"> 61</span><span class="keyword">template</span> <<span class="keyword">typename</span> SrcKeyType, <span class="keyword">typename</span> DstKeyType></div>
|
|
<div class="line"><a id="l00062" name="l00062"></a><span class="lineno"> 62</span><span class="keyword">auto</span> histogram_equalization(histogram<SrcKeyType> <span class="keyword">const</span>& src_hist, histogram<DstKeyType>& dst_hist)</div>
|
|
<div class="line"><a id="l00063" name="l00063"></a><span class="lineno"> 63</span> -> std::map<SrcKeyType, DstKeyType></div>
|
|
<div class="line"><a id="l00064" name="l00064"></a><span class="lineno"> 64</span>{</div>
|
|
<div class="line"><a id="l00065" name="l00065"></a><span class="lineno"> 65</span> <span class="keyword">static_assert</span>(</div>
|
|
<div class="line"><a id="l00066" name="l00066"></a><span class="lineno"> 66</span> std::is_integral<SrcKeyType>::value &&</div>
|
|
<div class="line"><a id="l00067" name="l00067"></a><span class="lineno"> 67</span> std::is_integral<DstKeyType>::value,</div>
|
|
<div class="line"><a id="l00068" name="l00068"></a><span class="lineno"> 68</span> <span class="stringliteral">"Source and destination histogram types are not appropriate"</span>);</div>
|
|
<div class="line"><a id="l00069" name="l00069"></a><span class="lineno"> 69</span> </div>
|
|
<div class="line"><a id="l00070" name="l00070"></a><span class="lineno"> 70</span> <span class="keyword">using </span>value_t = <span class="keyword">typename</span> histogram<SrcKeyType>::value_type;</div>
|
|
<div class="line"><a id="l00071" name="l00071"></a><span class="lineno"> 71</span> dst_hist.clear();</div>
|
|
<div class="line"><a id="l00072" name="l00072"></a><span class="lineno"> 72</span> <span class="keywordtype">double</span> sum = src_hist.sum();</div>
|
|
<div class="line"><a id="l00073" name="l00073"></a><span class="lineno"> 73</span> SrcKeyType min_key = (std::numeric_limits<DstKeyType>::min)();</div>
|
|
<div class="line"><a id="l00074" name="l00074"></a><span class="lineno"> 74</span> SrcKeyType max_key = (std::numeric_limits<DstKeyType>::max)();</div>
|
|
<div class="line"><a id="l00075" name="l00075"></a><span class="lineno"> 75</span> <span class="keyword">auto</span> cumltv_srchist = cumulative_histogram(src_hist);</div>
|
|
<div class="line"><a id="l00076" name="l00076"></a><span class="lineno"> 76</span> std::map<SrcKeyType, DstKeyType> color_map;</div>
|
|
<div class="line"><a id="l00077" name="l00077"></a><span class="lineno"> 77</span> std::for_each(cumltv_srchist.begin(), cumltv_srchist.end(), [&](value_t <span class="keyword">const</span>& v) {</div>
|
|
<div class="line"><a id="l00078" name="l00078"></a><span class="lineno"> 78</span> DstKeyType trnsfrmd_key =</div>
|
|
<div class="line"><a id="l00079" name="l00079"></a><span class="lineno"> 79</span> static_cast<DstKeyType>((v.second * (max_key - min_key)) / sum + min_key);</div>
|
|
<div class="line"><a id="l00080" name="l00080"></a><span class="lineno"> 80</span> color_map[std::get<0>(v.first)] = trnsfrmd_key;</div>
|
|
<div class="line"><a id="l00081" name="l00081"></a><span class="lineno"> 81</span> });</div>
|
|
<div class="line"><a id="l00082" name="l00082"></a><span class="lineno"> 82</span> std::for_each(src_hist.begin(), src_hist.end(), [&](value_t <span class="keyword">const</span>& v) {</div>
|
|
<div class="line"><a id="l00083" name="l00083"></a><span class="lineno"> 83</span> dst_hist[color_map[std::get<0>(v.first)]] += v.second;</div>
|
|
<div class="line"><a id="l00084" name="l00084"></a><span class="lineno"> 84</span> });</div>
|
|
<div class="line"><a id="l00085" name="l00085"></a><span class="lineno"> 85</span> <span class="keywordflow">return</span> color_map;</div>
|
|
<div class="line"><a id="l00086" name="l00086"></a><span class="lineno"> 86</span>}</div>
|
|
<div class="line"><a id="l00087" name="l00087"></a><span class="lineno"> 87</span> </div>
|
|
<div class="line"><a id="l00098" name="l00098"></a><span class="lineno"> 98</span><span class="keyword">template</span> <<span class="keyword">typename</span> SrcView, <span class="keyword">typename</span> DstView></div>
|
|
<div class="line"><a id="l00099" name="l00099"></a><span class="lineno"> 99</span><span class="keywordtype">void</span> histogram_equalization(</div>
|
|
<div class="line"><a id="l00100" name="l00100"></a><span class="lineno"> 100</span> SrcView <span class="keyword">const</span>& src_view,</div>
|
|
<div class="line"><a id="l00101" name="l00101"></a><span class="lineno"> 101</span> DstView <span class="keyword">const</span>& dst_view,</div>
|
|
<div class="line"><a id="l00102" name="l00102"></a><span class="lineno"> 102</span> std::size_t bin_width = 1,</div>
|
|
<div class="line"><a id="l00103" name="l00103"></a><span class="lineno"> 103</span> <span class="keywordtype">bool</span> mask = <span class="keyword">false</span>,</div>
|
|
<div class="line"><a id="l00104" name="l00104"></a><span class="lineno"> 104</span> std::vector<std::vector<bool>> src_mask = {})</div>
|
|
<div class="line"><a id="l00105" name="l00105"></a><span class="lineno"> 105</span>{</div>
|
|
<div class="line"><a id="l00106" name="l00106"></a><span class="lineno"> 106</span> gil_function_requires<ImageViewConcept<SrcView>>();</div>
|
|
<div class="line"><a id="l00107" name="l00107"></a><span class="lineno"> 107</span> gil_function_requires<MutableImageViewConcept<DstView>>();</div>
|
|
<div class="line"><a id="l00108" name="l00108"></a><span class="lineno"> 108</span> </div>
|
|
<div class="line"><a id="l00109" name="l00109"></a><span class="lineno"> 109</span> <span class="keyword">static_assert</span>(</div>
|
|
<div class="line"><a id="l00110" name="l00110"></a><span class="lineno"> 110</span> color_spaces_are_compatible<</div>
|
|
<div class="line"><a id="l00111" name="l00111"></a><span class="lineno"> 111</span> <span class="keyword">typename</span> color_space_type<SrcView>::type,</div>
|
|
<div class="line"><a id="l00112" name="l00112"></a><span class="lineno"> 112</span> <span class="keyword">typename</span> color_space_type<DstView>::type>::value,</div>
|
|
<div class="line"><a id="l00113" name="l00113"></a><span class="lineno"> 113</span> <span class="stringliteral">"Source and destination views must have same color space"</span>);</div>
|
|
<div class="line"><a id="l00114" name="l00114"></a><span class="lineno"> 114</span> </div>
|
|
<div class="line"><a id="l00115" name="l00115"></a><span class="lineno"> 115</span> <span class="comment">// Defining channel type</span></div>
|
|
<div class="line"><a id="l00116" name="l00116"></a><span class="lineno"> 116</span> <span class="keyword">using </span>source_channel_t = <span class="keyword">typename</span> channel_type<SrcView>::type;</div>
|
|
<div class="line"><a id="l00117" name="l00117"></a><span class="lineno"> 117</span> <span class="keyword">using </span>dst_channel_t = <span class="keyword">typename</span> channel_type<DstView>::type;</div>
|
|
<div class="line"><a id="l00118" name="l00118"></a><span class="lineno"> 118</span> <span class="keyword">using </span>coord_t = <span class="keyword">typename</span> SrcView::x_coord_t;</div>
|
|
<div class="line"><a id="l00119" name="l00119"></a><span class="lineno"> 119</span> </div>
|
|
<div class="line"><a id="l00120" name="l00120"></a><span class="lineno"> 120</span> std::size_t <span class="keyword">const</span> channels = num_channels<SrcView>::value;</div>
|
|
<div class="line"><a id="l00121" name="l00121"></a><span class="lineno"> 121</span> coord_t <span class="keyword">const</span> width = src_view.width();</div>
|
|
<div class="line"><a id="l00122" name="l00122"></a><span class="lineno"> 122</span> coord_t <span class="keyword">const</span> height = src_view.height();</div>
|
|
<div class="line"><a id="l00123" name="l00123"></a><span class="lineno"> 123</span> std::size_t pixel_max = (std::numeric_limits<dst_channel_t>::max)();</div>
|
|
<div class="line"><a id="l00124" name="l00124"></a><span class="lineno"> 124</span> std::size_t pixel_min = (std::numeric_limits<dst_channel_t>::min)();</div>
|
|
<div class="line"><a id="l00125" name="l00125"></a><span class="lineno"> 125</span> </div>
|
|
<div class="line"><a id="l00126" name="l00126"></a><span class="lineno"> 126</span> <span class="keywordflow">for</span> (std::size_t i = 0; i < channels; i++)</div>
|
|
<div class="line"><a id="l00127" name="l00127"></a><span class="lineno"> 127</span> {</div>
|
|
<div class="line"><a id="l00128" name="l00128"></a><span class="lineno"> 128</span> histogram<source_channel_t> h;</div>
|
|
<div class="line"><a id="l00129" name="l00129"></a><span class="lineno"> 129</span> fill_histogram(nth_channel_view(src_view, i), h, bin_width, <span class="keyword">false</span>, <span class="keyword">false</span>, mask, src_mask);</div>
|
|
<div class="line"><a id="l00130" name="l00130"></a><span class="lineno"> 130</span> h.normalize();</div>
|
|
<div class="line"><a id="l00131" name="l00131"></a><span class="lineno"> 131</span> <span class="keyword">auto</span> h2 = cumulative_histogram(h);</div>
|
|
<div class="line"><a id="l00132" name="l00132"></a><span class="lineno"> 132</span> <span class="keywordflow">for</span> (std::ptrdiff_t src_y = 0; src_y < height; ++src_y)</div>
|
|
<div class="line"><a id="l00133" name="l00133"></a><span class="lineno"> 133</span> {</div>
|
|
<div class="line"><a id="l00134" name="l00134"></a><span class="lineno"> 134</span> <span class="keyword">auto</span> src_it = nth_channel_view(src_view, i).row_begin(src_y);</div>
|
|
<div class="line"><a id="l00135" name="l00135"></a><span class="lineno"> 135</span> <span class="keyword">auto</span> dst_it = nth_channel_view(dst_view, i).row_begin(src_y);</div>
|
|
<div class="line"><a id="l00136" name="l00136"></a><span class="lineno"> 136</span> <span class="keywordflow">for</span> (std::ptrdiff_t src_x = 0; src_x < width; ++src_x)</div>
|
|
<div class="line"><a id="l00137" name="l00137"></a><span class="lineno"> 137</span> {</div>
|
|
<div class="line"><a id="l00138" name="l00138"></a><span class="lineno"> 138</span> <span class="keywordflow">if</span> (mask && !src_mask[src_y][src_x])</div>
|
|
<div class="line"><a id="l00139" name="l00139"></a><span class="lineno"> 139</span> dst_it[src_x][0] = channel_convert<dst_channel_t>(src_it[src_x][0]);</div>
|
|
<div class="line"><a id="l00140" name="l00140"></a><span class="lineno"> 140</span> <span class="keywordflow">else</span></div>
|
|
<div class="line"><a id="l00141" name="l00141"></a><span class="lineno"> 141</span> dst_it[src_x][0] = <span class="keyword">static_cast<</span>dst_channel_t<span class="keyword">></span>(</div>
|
|
<div class="line"><a id="l00142" name="l00142"></a><span class="lineno"> 142</span> h2[src_it[src_x][0]] * (pixel_max - pixel_min) + pixel_min);</div>
|
|
<div class="line"><a id="l00143" name="l00143"></a><span class="lineno"> 143</span> }</div>
|
|
<div class="line"><a id="l00144" name="l00144"></a><span class="lineno"> 144</span> }</div>
|
|
<div class="line"><a id="l00145" name="l00145"></a><span class="lineno"> 145</span> }</div>
|
|
<div class="line"><a id="l00146" name="l00146"></a><span class="lineno"> 146</span>}</div>
|
|
<div class="line"><a id="l00147" name="l00147"></a><span class="lineno"> 147</span> </div>
|
|
<div class="line"><a id="l00148" name="l00148"></a><span class="lineno"> 148</span>}} <span class="comment">//namespace boost::gil</span></div>
|
|
<div class="line"><a id="l00149" name="l00149"></a><span class="lineno"> 149</span> </div>
|
|
<div class="line"><a id="l00150" name="l00150"></a><span class="lineno"> 150</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  <a href="http://www.doxygen.org/index.html">doxygen</a> 1.9.8
|
|
</small>
|
|
</address>
|
|
</body>
|
|
</html>
|