filters_8cpp-source.html
来自「显著区域检测。求的图像中感兴趣区域的位置」· HTML 代码 · 共 297 行 · 第 1/2 页
HTML
297 行
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"><html><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8"><title>SaliencyToolbox: Filters.cpp Source File</title><link href="doxygen.css" rel="stylesheet" type="text/css"><link href="tabs.css" rel="stylesheet" type="text/css"></head><body><!-- Generated by Doxygen 1.5.2 --><div class="tabs"> <ul> <li><a href="index.html"><span>Main Page</span></a></li> <li><a href="annotated.html"><span>Data Structures</span></a></li> <li class="current"><a href="files.html"><span>Files</span></a></li> </ul></div><div class="tabs"> <ul> <li><a href="files.html"><span>File List</span></a></li> <li><a href="globals.html"><span>Globals</span></a></li> </ul></div><h1>Filters.cpp</h1><a href="Filters_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <a name="l00004"></a>00004 <span class="comment">// This file is part of the SaliencyToolbox - Copyright (C) 2006-2007</span><a name="l00005"></a>00005 <span class="comment">// by Dirk B. Walther and the California Institute of Technology.</span><a name="l00006"></a>00006 <span class="comment">// See the enclosed LICENSE.TXT document for the license agreement. </span><a name="l00007"></a>00007 <span class="comment">// More information about this project is available at: </span><a name="l00008"></a>00008 <span class="comment">// http://www.saliencytoolbox.net</span><a name="l00009"></a>00009 <a name="l00010"></a>00010 <span class="preprocessor">#include "<a class="code" href="mexLog_8h.html">mexLog.h</a>"</span><a name="l00011"></a>00011 <span class="preprocessor">#include "<a class="code" href="Image_8h.html">Image.h</a>"</span><a name="l00012"></a>00012 <span class="preprocessor">#include "<a class="code" href="Filters_8h.html">Filters.h</a>"</span><a name="l00013"></a>00013 <a name="l00014"></a>00014 <span class="preprocessor">#include <algorithm></span><a name="l00015"></a>00015 <span class="preprocessor">#include <cmath></span><a name="l00016"></a>00016 <span class="preprocessor">#include <limits></span><a name="l00017"></a>00017 <a name="l00018"></a>00018 <span class="comment">// ######################################################################</span><a name="l00019"></a>00019 <span class="comment">// kernel: 1 5 10 10 5 1</span><a name="l00020"></a><a class="code" href="Filters_8h.html#36928c645d464b6dfed7dfd50a5e2d28">00020</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#36928c645d464b6dfed7dfd50a5e2d28" title="Low-pass with a 6x6 separable kernel in Y and decimate in Y by factor 2.">lowPass6yDecY</a>(<span class="keyword">const</span> <a class="code" href="classImage.html" title="Provides access and processing functions for an image.">Image</a>& src)<a name="l00021"></a>00021 {<a name="l00022"></a>00022 <a class="code" href="mexLog_8h.html#092bd7bb2cb7fd9f483b6995cee61bc0">ASSERT</a>(<a class="code" href="Filters_8h.html#7c74fd3e38a3d834e5f48766f2896ed7">ecxStr</a>);<a name="l00023"></a>00023 <span class="keyword">const</span> <span class="keywordtype">int</span> w = src.<a class="code" href="classImage.html#36bcae4838228d574738249a87dc4464" title="Returns the width of the Image.">getWidth</a>(), hs = src.<a class="code" href="classImage.html#1cd7587f88b6932b5269051d1ba08ff9" title="Returns the height of the Image.">getHeight</a>();<a name="l00024"></a>00024 <span class="keyword">const</span> <span class="keywordtype">float</span> ecw = <a class="code" href="Filters_8h.html#a96a52e9dc841680c001410089a46154">ecx</a> / w;<a name="l00025"></a>00025 <span class="keywordtype">int</span> hr = hs / 2;<a name="l00026"></a>00026 <span class="keywordflow">if</span> (hr == 0) hr = 1;<a name="l00027"></a>00027 <a name="l00028"></a>00028 <a class="code" href="classImage.html" title="Provides access and processing functions for an image.">Image</a> result(w,hr);<a name="l00029"></a>00029 <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="l00030"></a>00030 <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="l00031"></a>00031 <a name="l00032"></a>00032 <span class="keywordflow">if</span> (hs <= 1)<a name="l00033"></a>00033 result = src;<a name="l00034"></a>00034 <span class="keywordflow">else</span> <span class="keywordflow">if</span> (hs == 2)<a name="l00035"></a>00035 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> x = 0; x < w; ++x)<a name="l00036"></a>00036 {<a name="l00037"></a>00037 <span class="comment">// use kernel [1 1]^T / 2</span><a name="l00038"></a>00038 *rptr++ = (sptr[0] + sptr[1]) / 2.0;<a name="l00039"></a>00039 sptr += 2;<a name="l00040"></a>00040 }<a name="l00041"></a>00041 <span class="keywordflow">else</span> <span class="keywordflow">if</span> (hs == 3)<a name="l00042"></a>00042 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> x = 0; x < w; ++x)<a name="l00043"></a>00043 {<a name="l00044"></a>00044 <span class="comment">// use kernel [1 2 1]^T / 4</span><a name="l00045"></a>00045 *rptr++ = (sptr[0] + sptr[1] * 2.0 + sptr[2]) / 4.0;<a name="l00046"></a>00046 sptr += 3;<a name="l00047"></a>00047 }<a name="l00048"></a>00048 <span class="keywordflow">else</span> <span class="comment">// general case with hs >= 4</span><a name="l00049"></a>00049 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> x = 0; x < w; ++x)<a name="l00050"></a>00050 {<a name="l00051"></a>00051 <span class="comment">// top most point - use kernel [10 10 5 1]^T / 26</span><a name="l00052"></a>00052 *rptr++ = ((sptr[0] + sptr[1]) * 10.0 + <a name="l00053"></a>00053 sptr[2] * 5.0 + sptr[3]) / 26.0;<a name="l00054"></a>00054 <span class="comment">//++sptr;</span><a name="l00055"></a>00055 <a name="l00056"></a>00056 <span class="comment">// general case</span><a name="l00057"></a>00057 <span class="keywordtype">int</span> y;<a name="l00058"></a>00058 <span class="keywordflow">for</span> (y = 0; y < (hs - 5); y += 2)<a name="l00059"></a>00059 {<a name="l00060"></a>00060 <span class="comment">// use kernel [1 5 10 10 5 1]^T / 32</span><a name="l00061"></a>00061 *rptr++ = ((sptr[1] + sptr[4]) * 5.0 +<a name="l00062"></a>00062 (sptr[2] + sptr[3]) * 10.0 +<a name="l00063"></a>00063 (sptr[0] + sptr[5])) / 32.0;<a name="l00064"></a>00064 sptr += 2;<a name="l00065"></a>00065 }<a name="l00066"></a>00066 <a name="l00067"></a>00067 <span class="comment">// find out how to treat the bottom most point</span><a name="l00068"></a>00068 <span class="keywordflow">if</span> (y == (hs - 5))<a name="l00069"></a>00069 {<a name="l00070"></a>00070 <span class="comment">// use kernel [1 5 10 10 5]^T / 31</span><a name="l00071"></a>00071 *rptr++ = ((sptr[1] + sptr[4]) * 5.0 +<a name="l00072"></a>00072 (sptr[2] + sptr[3]) * 10.0 +<a name="l00073"></a>00073 sptr[0]) / 31.0;<a name="l00074"></a>00074 sptr += 5;<a name="l00075"></a>00075 }<a name="l00076"></a>00076 <span class="keywordflow">else</span><a name="l00077"></a>00077 {<a name="l00078"></a>00078 <span class="comment">// use kernel [1 5 10 10]^T / 26</span><a name="l00079"></a>00079 *rptr++ = ( sptr[0] + sptr[1] * 5.0 +<a name="l00080"></a>00080 (sptr[2] + sptr[3]) * 10.0) / 26.0;<a name="l00081"></a>00081 sptr += 4;<a name="l00082"></a>00082 }<a name="l00083"></a>00083 }<a name="l00084"></a>00084 <a name="l00085"></a>00085 <span class="keywordflow">return</span> result; <a name="l00086"></a>00086 }<a name="l00087"></a>00087 <a name="l00088"></a>00088 <span class="comment">// ######################################################################</span><a name="l00089"></a>00089 <span class="comment">// kernel: 1 5 10 10 5 1</span><a name="l00090"></a><a class="code" href="Filters_8h.html#fd05058fad5d4735eef3352b93dcce30">00090</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#fd05058fad5d4735eef3352b93dcce30" title="Low-pass with a 6x6 separable kernel in X and decimate in X by factor 2.">lowPass6xDecX</a>(<span class="keyword">const</span> <a class="code" href="classImage.html" title="Provides access and processing functions for an image.">Image</a>& src)<a name="l00091"></a>00091 {<a name="l00092"></a>00092 <a class="code" href="mexLog_8h.html#092bd7bb2cb7fd9f483b6995cee61bc0">ASSERT</a>(<a class="code" href="Filters_8h.html#7c74fd3e38a3d834e5f48766f2896ed7">ecxStr</a>);<a name="l00093"></a>00093 <span class="keyword">const</span> <span class="keywordtype">int</span> ws = src.<a class="code" href="classImage.html#36bcae4838228d574738249a87dc4464" title="Returns the width of the Image.">getWidth</a>(), h = src.<a class="code" href="classImage.html#1cd7587f88b6932b5269051d1ba08ff9" title="Returns the height of the Image.">getHeight</a>();<a name="l00094"></a>00094 <span class="keyword">const</span> <span class="keywordtype">float</span> ecw = <a class="code" href="Filters_8h.html#a96a52e9dc841680c001410089a46154">ecx</a> / ws;<a name="l00095"></a>00095 <span class="keyword">const</span> <span class="keywordtype">int</span> h2 = h * 2, h3 = h * 3, h4 = h * 4, h5 = h * 5;<a name="l00096"></a>00096 <span class="keywordtype">int</span> wr = ws / 2;<a name="l00097"></a>00097 <span class="keywordflow">if</span> (wr == 0) wr = 1;<a name="l00098"></a>00098 <a name="l00099"></a>00099 <a class="code" href="classImage.html" title="Provides access and processing functions for an image.">Image</a> result(wr,h);<a name="l00100"></a>00100 <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="l00101"></a>00101 <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="l00102"></a>00102 <a name="l00103"></a>00103 <span class="keywordflow">if</span> (ws <= 1)<a name="l00104"></a>00104 result = src;<a name="l00105"></a>00105 <span class="keywordflow">else</span> <span class="keywordflow">if</span> (ws == 2)<a name="l00106"></a>00106 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> y = 0; y < h; ++y)<a name="l00107"></a>00107 {<a name="l00108"></a>00108 <span class="comment">// use kernel [1 1] / 2</span><a name="l00109"></a>00109 *rptr++ = (sptr[0] + sptr[h]) / 2.0;<a name="l00110"></a>00110 ++sptr;<a name="l00111"></a>00111 }<a name="l00112"></a>00112 <span class="keywordflow">else</span> <span class="keywordflow">if</span> (ws == 3)<a name="l00113"></a>00113 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> y = 0; y < h; ++y)<a name="l00114"></a>00114 {<a name="l00115"></a>00115 <span class="comment">// use kernel [1 2 1] / 4</span><a name="l00116"></a>00116 *rptr++ = (sptr[0] + sptr[h] * 2.0 + sptr[h2]) / 4.0;<a name="l00117"></a>00117 ++sptr;<a name="l00118"></a>00118 }<a name="l00119"></a>00119 <span class="keywordflow">else</span> <span class="comment">// general case for ws >= 4</span><a name="l00120"></a>00120 {<a name="l00121"></a>00121 <span class="comment">// left most point - use kernel [10 10 5 1] / 26</span><a name="l00122"></a>00122 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> y = 0; y < h; ++y)<a name="l00123"></a>00123 {<a name="l00124"></a>00124 *rptr++ = ((sptr[0] + sptr[h]) * 10.0 + <a name="l00125"></a>00125 sptr[h2] * 5.0 + sptr[h3]) / 26.0;<a name="l00126"></a>00126 ++sptr;<a name="l00127"></a>00127 }<a name="l00128"></a>00128 sptr -= h;<a name="l00129"></a>00129 <a name="l00130"></a>00130 <span class="comment">// general case</span><a name="l00131"></a>00131 <span class="keywordtype">int</span> x;⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?