filters_8cpp-source.html

来自「显著区域检测。求的图像中感兴趣区域的位置」· HTML 代码 · 共 297 行 · 第 1/2 页

<a name="l00132"></a>00132       <span class="keywordflow">for</span> (x = 0; x &lt; (ws - 5); x += 2)
<a name="l00133"></a>00133         {
<a name="l00134"></a>00134           <span class="keywordflow">for</span> (<span class="keywordtype">int</span> y = 0; y &lt; h; ++y)
<a name="l00135"></a>00135             {
<a name="l00136"></a>00136               <span class="comment">// use kernel [1 5 10 10 5 1] / 32</span>
<a name="l00137"></a>00137               *rptr++ = ((sptr[h]  + sptr[h4])  *  5.0 +
<a name="l00138"></a>00138                          (sptr[h2] + sptr[h3])  * 10.0 +
<a name="l00139"></a>00139                          (sptr[0]  + sptr[h5])) / 32.0;
<a name="l00140"></a>00140               ++sptr;
<a name="l00141"></a>00141             }
<a name="l00142"></a>00142           sptr += h;
<a name="l00143"></a>00143         }
<a name="l00144"></a>00144         
<a name="l00145"></a>00145       <span class="comment">// find out how to treat the right most point</span>
<a name="l00146"></a>00146       <span class="keywordflow">if</span> (x == (ws - 5))
<a name="l00147"></a>00147         <span class="keywordflow">for</span> (<span class="keywordtype">int</span> y = 0; y &lt; h; ++y)
<a name="l00148"></a>00148           {
<a name="l00149"></a>00149             <span class="comment">// use kernel [1 5 10 10 5] / 31</span>
<a name="l00150"></a>00150             *rptr++ = ((sptr[h]  + sptr[h4])  *  5.0 +
<a name="l00151"></a>00151                        (sptr[h2] + sptr[h3])  * 10.0 +
<a name="l00152"></a>00152                         sptr[0]) / 31.0;
<a name="l00153"></a>00153             ++sptr;
<a name="l00154"></a>00154           }
<a name="l00155"></a>00155       <span class="keywordflow">else</span>
<a name="l00156"></a>00156         <span class="keywordflow">for</span> (<span class="keywordtype">int</span> y = 0; y &lt; h; ++y)
<a name="l00157"></a>00157           {
<a name="l00158"></a>00158             <span class="comment">// use kernel [1 5 10 10] / 26</span>
<a name="l00159"></a>00159             *rptr++ = ( sptr[0]  + sptr[h]   * 5.0 + 
<a name="l00160"></a>00160                        (sptr[h2] + sptr[h3]) * 10.0) / 26.0;
<a name="l00161"></a>00161             ++sptr;
<a name="l00162"></a>00162           }
<a name="l00163"></a>00163     }
<a name="l00164"></a>00164   <span class="keywordflow">return</span> result;
<a name="l00165"></a>00165 }
<a name="l00166"></a>00166 
<a name="l00167"></a>00167 <span class="comment">// ######################################################################</span>
<a name="l00168"></a><a class="code" href="Filters_8h.html#40cfbd12584dfa6a575c85c97a03aa18">00168</a> <a class="code" href="classImage.html" title="Provides access and processing functions for an image.">Image</a> <a class="code" href="Filters_8cpp.html#40cfbd12584dfa6a575c85c97a03aa18" title="2d convolution that avoids bleeding energy over the edge.">conv2PreserveEnergy</a>(<span class="keyword">const</span> <a class="code" href="classImage.html" title="Provides access and processing functions for an image.">Image</a>&amp; src, <span class="keyword">const</span> <a class="code" href="classImage.html" title="Provides access and processing functions for an image.">Image</a>&amp; f)
<a name="l00169"></a>00169 {
<a name="l00170"></a>00170   <a class="code" href="mexLog_8h.html#092bd7bb2cb7fd9f483b6995cee61bc0">ASSERT</a>(src.<a class="code" href="classImage.html#b87bf4b103cb5a9cf090f8729f362335" title="Returns whether the Image is initialized.">isInitialized</a>());
<a name="l00171"></a>00171   <a class="code" href="mexLog_8h.html#092bd7bb2cb7fd9f483b6995cee61bc0">ASSERT</a>(<a class="code" href="Filters_8h.html#7c74fd3e38a3d834e5f48766f2896ed7">ecxStr</a>);
<a name="l00172"></a>00172 
<a name="l00173"></a>00173   <span class="keyword">const</span> <span class="keywordtype">int</span> sw = src.<a class="code" href="classImage.html#36bcae4838228d574738249a87dc4464" title="Returns the width of the Image.">getWidth</a>();
<a name="l00174"></a>00174   <span class="keyword">const</span> <span class="keywordtype">int</span> sh = src.<a class="code" href="classImage.html#1cd7587f88b6932b5269051d1ba08ff9" title="Returns the height of the Image.">getHeight</a>();
<a name="l00175"></a>00175   <span class="keyword">const</span> <span class="keywordtype">float</span> ecw = <a class="code" href="Filters_8h.html#a96a52e9dc841680c001410089a46154">ecx</a> / sw;
<a name="l00176"></a>00176 
<a name="l00177"></a>00177   <a class="code" href="classImage.html#f9f1ae298f5d8f8b44e4e34a72f61c2f" title="Read-only iterator.">Image::const_iterator</a> filter = f.<a class="code" href="classImage.html#e0666b33daa7ba55e3348d4fff02356f" title="Returns a read-only iterator to the beginning of the Image data.">begin</a>();
<a name="l00178"></a>00178   <span class="keyword">const</span> <span class="keywordtype">int</span> fw = f.<a class="code" href="classImage.html#36bcae4838228d574738249a87dc4464" title="Returns the width of the Image.">getWidth</a>();
<a name="l00179"></a>00179   <span class="keyword">const</span> <span class="keywordtype">int</span> fh = f.<a class="code" href="classImage.html#1cd7587f88b6932b5269051d1ba08ff9" title="Returns the height of the Image.">getHeight</a>();
<a name="l00180"></a>00180 
<a name="l00181"></a>00181   <a class="code" href="mexLog_8h.html#092bd7bb2cb7fd9f483b6995cee61bc0">ASSERT</a>((fw &amp; 1) &amp;&amp; (fh &amp; 1));
<a name="l00182"></a>00182 
<a name="l00183"></a>00183   <a class="code" href="classImage.html" title="Provides access and processing functions for an image.">Image</a> result(sw, sh);
<a name="l00184"></a>00184   <a class="code" href="classImage.html#f9f1ae298f5d8f8b44e4e34a72f61c2f" title="Read-only iterator.">Image::const_iterator</a> sptr = src.<a class="code" href="classImage.html#e0666b33daa7ba55e3348d4fff02356f" title="Returns a read-only iterator to the beginning of the Image data.">begin</a>();
<a name="l00185"></a>00185   <a class="code" href="classImage.html#01d14564f4b61f97d305e293716933cf" title="Read/write iterator.">Image::iterator</a> rptr = result.<a class="code" href="classImage.html#9987bf538358f25470af7d628e97694a" title="Returns a read/write iterator to the beginning of the Image data.">beginw</a>();
<a name="l00186"></a>00186 
<a name="l00187"></a>00187   <span class="keyword">const</span> <span class="keywordtype">int</span> fend = fw * fh - 1;
<a name="l00188"></a>00188   <span class="keyword">const</span> <span class="keywordtype">int</span> fw2 = (fw - 1) / 2;
<a name="l00189"></a>00189   <span class="keyword">const</span> <span class="keywordtype">int</span> fh2 = (fh - 1) / 2;
<a name="l00190"></a>00190 
<a name="l00191"></a>00191   <span class="keyword">const</span> <span class="keywordtype">int</span> sSkip = sh - fh;
<a name="l00192"></a>00192 
<a name="l00193"></a>00193   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> rx = 0; rx &lt; sw; ++rx)
<a name="l00194"></a>00194     {
<a name="l00195"></a>00195       <span class="comment">// Determine if we're safely inside the image in the x-direction:</span>
<a name="l00196"></a>00196       <span class="keyword">const</span> <span class="keywordtype">bool</span> isXclean = ((rx &gt;= fw2) &amp;&amp; (rx &lt;  (sw - fw2)));
<a name="l00197"></a>00197       <span class="keywordflow">for</span> (<span class="keywordtype">int</span> ry = 0; ry &lt; sh; ++ry, ++rptr)
<a name="l00198"></a>00198         {
<a name="l00199"></a>00199           <span class="comment">// Determine if we're safely inside the image in the y-direction:</span>
<a name="l00200"></a>00200           <span class="keyword">const</span> <span class="keywordtype">bool</span> isYclean = ((ry &gt;= fh2) &amp;&amp; (ry &lt;  (sh - fh2)));
<a name="l00201"></a>00201           
<a name="l00202"></a>00202           <span class="keywordflow">if</span> (isXclean &amp;&amp; isYclean)
<a name="l00203"></a>00203             {
<a name="l00204"></a>00204               <span class="comment">// well inside the image: use straight-forward convolution</span>
<a name="l00205"></a>00205               <span class="keywordtype">float</span> rval = 0.0f;
<a name="l00206"></a>00206               <a class="code" href="classImage.html#f9f1ae298f5d8f8b44e4e34a72f61c2f" title="Read-only iterator.">Image::const_iterator</a> fptr = filter+fend;
<a name="l00207"></a>00207               <a class="code" href="classImage.html#f9f1ae298f5d8f8b44e4e34a72f61c2f" title="Read-only iterator.">Image::const_iterator</a> sptr2 = sptr + sh*(rx-fw2) + ry - fh2;
<a name="l00208"></a>00208 
<a name="l00209"></a>00209               <span class="keywordflow">for</span> (<span class="keywordtype">int</span> fx = 0; fx &lt; fw; ++fx)
<a name="l00210"></a>00210                 {
<a name="l00211"></a>00211                   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> fy = 0; fy &lt; fh; ++fy)
<a name="l00212"></a>00212                     rval += (*sptr2++) * (*fptr--);
<a name="l00213"></a>00213 
<a name="l00214"></a>00214                   sptr2 += sSkip;
<a name="l00215"></a>00215                 }
<a name="l00216"></a>00216               *rptr = rval;
<a name="l00217"></a>00217               <span class="keywordflow">continue</span>;
<a name="l00218"></a>00218             }
<a name="l00219"></a>00219           <span class="keywordflow">else</span>
<a name="l00220"></a>00220             {
<a name="l00221"></a>00221               <span class="comment">// near the border: assume that the value missing pixels</span>
<a name="l00222"></a>00222               <span class="comment">// is equal to the average over the present pixels</span>
<a name="l00223"></a>00223               <span class="comment">// to minimize artifacts at the edge.</span>
<a name="l00224"></a>00224 
<a name="l00225"></a>00225               <span class="keywordtype">float</span> rval = 0.0f;
<a name="l00226"></a>00226               <span class="keywordtype">float</span> sSum = 0.0f;
<a name="l00227"></a>00227               <span class="keywordtype">int</span> sCount = 0;
<a name="l00228"></a>00228               <span class="keywordtype">float</span> fsum_skipped = 0.0f;
<a name="l00229"></a>00229 
<a name="l00230"></a>00230               <span class="keywordflow">for</span> (<span class="keywordtype">int</span> fx = 0; fx &lt; fw; ++fx)
<a name="l00231"></a>00231                 {
<a name="l00232"></a>00232                   <span class="keyword">const</span> <span class="keywordtype">int</span> sx = rx + fx - fw2;
<a name="l00233"></a>00233                   <span class="keywordflow">if</span> (sx &gt;= 0 &amp;&amp; sx &lt; sw)
<a name="l00234"></a>00234                     {
<a name="l00235"></a>00235                       <span class="keywordflow">for</span> (<span class="keywordtype">int</span> fy = 0; fy &lt; fh; ++fy)
<a name="l00236"></a>00236                         {
<a name="l00237"></a>00237                           <span class="keyword">const</span> <span class="keywordtype">float</span> fil = filter[fend - fx * fh - fy];
<a name="l00238"></a>00238                           <span class="keyword">const</span> <span class="keywordtype">int</span> sy = ry + fy - fh2;
<a name="l00239"></a>00239                           <span class="keywordflow">if</span> (sy &gt;= 0 &amp;&amp; sy &lt; sh)
<a name="l00240"></a>00240                             {
<a name="l00241"></a>00241                               <span class="comment">// here we have the pixel present: </span>
<a name="l00242"></a>00242                               <span class="comment">// store it and count it to compute the mean</span>
<a name="l00243"></a>00243                               <span class="keyword">const</span> <span class="keywordtype">float</span> sVal = sptr[sx * sh + sy];
<a name="l00244"></a>00244                               rval += sVal * fil;
<a name="l00245"></a>00245                               sSum += sVal;
<a name="l00246"></a>00246                               ++sCount;
<a name="l00247"></a>00247                             }
<a name="l00248"></a>00248                           <span class="keywordflow">else</span>
<a name="l00249"></a>00249                             {
<a name="l00250"></a>00250                               <span class="comment">// this is accumulataing the filter values </span>
<a name="l00251"></a>00251                               <span class="comment">// over missing pixels</span>
<a name="l00252"></a>00252                               fsum_skipped += fil;
<a name="l00253"></a>00253                             }
<a name="l00254"></a>00254                         }
<a name="l00255"></a>00255                     }
<a name="l00256"></a>00256                   <span class="keywordflow">else</span>
<a name="l00257"></a>00257                     {
<a name="l00258"></a>00258                       <span class="comment">// this is accumulating the filter values</span>
<a name="l00259"></a>00259                       <span class="comment">// over an entire column of mixxing pixels</span>
<a name="l00260"></a>00260                       <span class="keywordflow">for</span> (<span class="keywordtype">int</span> fy = 0; fy &lt; fh; ++fy)
<a name="l00261"></a>00261                         fsum_skipped += filter[fend - fx * fh - fy];
<a name="l00262"></a>00262                     }
<a name="l00263"></a>00263                 }
<a name="l00264"></a>00264               <span class="comment">// compute the average over the present pixels</span>
<a name="l00265"></a>00265               <span class="keyword">const</span> <span class="keywordtype">float</span> sAvg = sSum / sCount;
<a name="l00266"></a>00266               
<a name="l00267"></a>00267               <span class="comment">// add the that average x accumulated filter values to the result</span>
<a name="l00268"></a>00268               *rptr = rval + (fsum_skipped * sAvg);
<a name="l00269"></a>00269             }
<a name="l00270"></a>00270         }
<a name="l00271"></a>00271     }
<a name="l00272"></a>00272   <span class="keywordflow">return</span> result;
<a name="l00273"></a>00273 }
</pre></div><hr size="1"><address style="text-align: right;"><small>Generated on Fri Sep 7 13:09:49 2007 for SaliencyToolbox by&nbsp;
<a href="http://www.doxygen.org/index.html">
<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.2 </small></address>
</body>
</html>