图像修复Criminisi算法

基于纹理合成的图像修复算法，其中最基础的就是Criminisi算法，该文件包含MATLAB代码，能处理彩色图像，有测试图片及最后的效果图。代码没有问题，可以直接使用。

/*
 * A best exemplar finder.  Scans over the entire image （using a
 * sliding window） and finds the exemplar which minimizes the sum
 * squared error （SSE） over the to-be-filled pixels in the target
 * patch. 
 *
 * @author Sooraj Bhat
 */
#include “mex.h“
#include 
#include “matrix.h“

void bestexemplarhelper（const int mm const int nn const int m const int n 
			const double *img const double *Ip 
			const mxLogical *toFill const mxLogical *sourceRegion         
			double *best） 
{
  register int ijiijjii2jj2MNIJndxndx2mn=m*nmmnn=mm*nn;
  double patchErr=0.0err=0.0bestErr=1000000000.0;

  /* foreach patch */
  N=nn-n+1;  M=mm-m+1;
  for （j=1; j<=N; ++j） 
  {
    J=j+n-1;
    for （i=1; i<=M; ++i） 
    {
      I=i+m-1;
      /*** Calculate patch error ***/
      /* foreach pixel in the current patch */
      for （jj=jjj2=1; jj<=J; ++jj++jj2） 
      {
	for （ii=iii2=1; ii<=I; ++ii++ii2） 
    {
	  ndx=ii-1+mm*（jj-1）;
	  if （!sourceRegion[ndx]）
	    goto skipPatch;
	  ndx2=ii2-1+m*（jj2-1）;
	  if （!toFill[ndx2]） 
      {
	    err=img[ndx      ] - Ip[ndx2    ]; patchErr += err*err;
	    err=img[ndx+=mmnn] - Ip[ndx2+=mn]; patchErr += err*err;
	    err=img[ndx+=mmnn] - Ip[ndx2+=mn]; patchErr += err*err;
	  }
	}
      }
      /*** Update ***/
      if （patchErr < bestErr） 
      {
	bestErr = patchErr; 
	best[0] = i; best[1] = I;
	best[2] = j; best[3] = J;
      }
      /*** Reset ***/
    skipPatch:
      patchErr = 0.0; 
    }
  }
}

/* best = bestexemplarhelper（mmnnmnimgIptoFillsourceRegion）; */
void mexFunction（int nlhsmxArray *plhs[]int nrhsconst mxArray *prhs[]） 
{
  int mmnnmn;
  double *img*Ip*best;
  mxLogical *toFill*sourceRegion;

  /* Extract the inputs */
  mm = （int）mxGetScalar（prhs[0]）;/*The value of the first real （nonimaginary） element of the mxArray.*/
  nn = （int）mxGetScalar（prhs[1]）;
  m  = （int）mxGetScalar（prhs[2]）;
  n  = （int）mxGetScalar（prhs[3]）;
  img = mxGetPr（prhs[4]）;/*The address of the first element of the real data*/
  Ip  = mxGetPr（prhs[5]）;
  toFill = mxGetLogicals（prhs[6]）;/*mxGetLogicals（）The address of the first logical element in the mxArray. The result is unspecified if the mxArray is not a logical array.*/
  sourceRegion = mxGetLogicals（prhs[7]）;
  
  /* Setup the output */
  plhs[0] = mxCreateDoubleMatrix（41mxREAL）;
  best = mxGetPr（plhs[0]）;
  best[0]=best[1]=best[2]=best[3]=0.0;

  /* Do the actual work */
  bestexemplarhelper（mmnnmnimgIptoFillsourceRegionbest）;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      162139  2015-03-26 10:20  20.png
     文件      136031  2015-03-26 10:20  201.png
     文件       37938  2015-03-29 20:15  benggee_7.jpg
     文件        2382  2015-04-06 14:14  bestexemplarhelper.asv
     文件        2503  2015-04-30 11:56  bestexemplarhelper.c
     文件        8536  2013-03-12 13:12  bestexemplarhelper.mexmaci64
     文件        6656  2015-03-26 10:46  bestexemplarhelper.mexw32
     文件         293  2015-04-30 21:37  Evaluate.m
     文件        8928  2015-05-03 23:46  inpaint7.asv
     文件        8928  2015-05-03 23:49  inpaint7.m
     文件         486  2015-04-11 14:39  plotall.m
     文件        1063  2004-11-29 23:04  README.txt