Itti 的算法代码

Itti是Koch的学生，现已由他的显著性视觉模型扬名天下，附件即是他的代码

/*!@file Filters.cpp various filter functions.
 */

// This file is part of the SaliencyToolbox - Copyright （C） 2006-2007
// by Dirk B. Walther and the California Institute of Technology.
// See the enclosed LICENSE.TXT document for the license agreement. 
// More information about this project is available at: 
// http://www.saliencytoolbox.net

#include “mexLog.h“
#include “Image.h“
#include “Filters.h“

#include 
#include 
#include 

// ######################################################################
// kernel: 1 5 10 10 5 1
Image lowPass6yDecY（const Image& src）
{
  ASSERT（ecxStr）;
  const int w = src.getWidth（） hs = src.getHeight（）;
  const float ecw = ecx / w;
  int hr = hs / 2;
  if （hr == 0） hr = 1;
  
  Image result（whr）;
  Image::iterator rptr = result.beginw（）;
  Image::const_iterator sptr = src.begin（）;
  
  if （hs <= 1）
    result = src;
  else if （hs == 2）
    for （int x = 0; x < w; ++x）
      {
        // use kernel [1 1]^T / 2
        *rptr++ = （sptr[0] + sptr[1]） / 2.0;
        sptr += 2;
      }
  else if （hs == 3）
    for （int x = 0; x < w; ++x）
      {
        // use kernel [1 2 1]^T / 4
        *rptr++ = （sptr[0] + sptr[1] * 2.0 + sptr[2]） / 4.0;
        sptr += 3;
      }
  else // general case with hs >= 4
    for （int x = 0; x < w; ++x）
      {
        // top most point - use kernel [10 10 5 1]^T / 26
        *rptr++ = （（sptr[0] + sptr[1]） * 10.0 + 
                    sptr[2] * 5.0 + sptr[3]） / 26.0;
        //++sptr;
        
        // general case
        int y;
        for （y = 0; y < （hs - 5）; y += 2）
          {
            // use kernel [1 5 10 10 5 1]^T / 32
            *rptr++ = （（sptr[1] + sptr[4]）  *  5.0 +
                       （sptr[2] + sptr[3]）  * 10.0 +
                       （sptr[0] + sptr[5]）） / 32.0;
            sptr += 2;
          }
        
        // find out how to treat the bottom most point
        if （y == （hs - 5））
          {
            // use kernel [1 5 10 10 5]^T / 31
             *rptr++ = （（sptr[1] + sptr[4]）  *  5.0 +
                        （sptr[2] + sptr[3]）  * 10.0 +
                         sptr[0]）            / 31.0;
            sptr += 5;
          }
        else
          {
            // use kernel [1 5 10 10]^T / 26
            *rptr++ = （ sptr[0] + sptr[1]  *  5.0 +
                       （sptr[2] + sptr[3]） * 10.0） / 26.0;
            sptr += 4;
          }
        }
            
   return result;      
}

// ######################################################################
// kernel: 1 5 10 10 5 1
Image lowPass6xDecX（const Image& src）
{
  ASSERT（ecxStr）;
  const int ws = src.getWidth（） h = src.getHeight（）;
  const float ecw = ecx / ws;
  const int h2 = h * 2 h3 = h * 3 h4 = h * 4 h5 = h * 5;
  int wr = ws / 2;
  if （wr == 0） wr = 1;
  
  Image result（wrh）;
  Image::iterator rptr = result.beginw（）;
  Image::const_iterator sptr = src.begin（）;

  if （ws <= 1）
    result = src;
  else if （ws == 2）
    for （int y = 0; y < h; ++y）
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2007-09-07 15:02  SaliencyToolbox\bin\
     文件     1515874  2007-07-20 16:46  SaliencyToolbox\bin\mexConv2PreserveEnergy.dll
     文件       33750  2007-07-20 16:51  SaliencyToolbox\bin\mexConv2PreserveEnergy.mexa64
     文件       27960  2007-07-20 16:50  SaliencyToolbox\bin\mexConv2PreserveEnergy.mexglx
     文件       34460  2007-07-20 16:48  SaliencyToolbox\bin\mexConv2PreserveEnergy.mexmac
     文件       33888  2007-07-20 16:48  SaliencyToolbox\bin\mexConv2PreserveEnergy.mexmaci
     文件     1515868  2007-07-20 16:46  SaliencyToolbox\bin\mexGaussianSubsample.dll
     文件       33744  2007-07-20 16:51  SaliencyToolbox\bin\mexGaussianSubsample.mexa64
     文件       27922  2007-07-20 16:50  SaliencyToolbox\bin\mexGaussianSubsample.mexglx
     文件       34456  2007-07-20 16:48  SaliencyToolbox\bin\mexGaussianSubsample.mexmac
     文件       33884  2007-07-20 16:48  SaliencyToolbox\bin\mexGaussianSubsample.mexmaci
     文件     1515813  2007-07-20 16:46  SaliencyToolbox\bin\mexLocalMaxima.dll
     文件       33702  2007-07-20 16:51  SaliencyToolbox\bin\mexLocalMaxima.mexa64
     文件       27952  2007-07-20 16:50  SaliencyToolbox\bin\mexLocalMaxima.mexglx
     文件       34444  2007-07-20 16:48  SaliencyToolbox\bin\mexLocalMaxima.mexmac
     文件       33872  2007-07-20 16:48  SaliencyToolbox\bin\mexLocalMaxima.mexmaci
     文件        2227  2007-09-07 15:00  SaliencyToolbox\CHANGELOG
     目录           0  2007-09-07 15:02  SaliencyToolbox\doc\
     文件       12612  2007-09-07 14:47  SaliencyToolbox\doc\index.html
     目录           0  2007-09-07 15:02  SaliencyToolbox\doc\mdoc\
     文件         273  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\alpha.png
     文件         327  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\c++.png
     文件         252  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\c.png
     文件         214  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\demoicon.gif
     文件         133  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\down.png
     文件         265  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\fortran.png
     文件         255  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\hp.png
     文件        7968  2007-09-07 14:43  SaliencyToolbox\doc\mdoc\index.html
     文件         136  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\left.png
     文件         272  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\linux.png
     文件        1002  2007-06-19 09:35  SaliencyToolbox\doc\mdoc\m2html.css
............此处省略396个文件信息