matlab torr3D reconstruction 工具箱

matlab三维重建工具箱，欢迎初学机器视觉的朋友学习。

/*
patch_match.c     implements ...

  takes two images and does sum of differences of intensities at windows specified in each

  2001 written by Phil Torr
  Microsoft Research Cambridge
*/

//if correlation exceeds jump out then stop
#include 
#include 
#include “mex.h“

void mexFunction （
				  int		nlhs		   /* number of expected outputs */
				  mxArray	**plhs 	   /* matrix pointer array returning outputs */
				  int		nrhs		   /* number of inputs */
				  const mxArray **prhs	   /* matrix pointer array for inputs */
				  ） {
	int       width height i j l k border =1;
	double    *out *im1 *im2;
	bool is_max;
	int x1y1x2y2x1iy1ix2iy2ihalf_sizejump_out;
	double corr = 0.0;
	
	/* parameter checks */
	if （（nrhs < 7） || （nlhs != 1）） {
		mexErrMsgTxt （“Usage: Y = patch_match（im1im2x1y1x2y2half_sizejump_out）\n\n“）;
		return;
	}
	

	//not any checking here yet!!!!

	/* reading the parameters */
	height = mxGetM （prhs [0]）;
	width = mxGetN （prhs [0]）;
	im1 = （double *） mxGetPr （prhs [0]）;
	im2 = （double *） mxGetPr （prhs [1]）;
	x1 = （int）（mxGetScalar（prhs[2]））-1;
	y1 = （int）（mxGetScalar（prhs[3]））-1;
	x2 = （int）（mxGetScalar（prhs[4]））-1;
	y2 = （int）（mxGetScalar（prhs[5]））-1;
	half_size = （int）（mxGetScalar（prhs[6]））;
	jump_out = （double）（mxGetScalar（prhs[7]））;

//	printf（“x1 y1 x2 y2 %d %d %d %d\n“x1y1 x2y2）;
	
	//out = mxGetScalar（plhs [0]）;

	/* require memory for return */
	plhs [0] = mxCreateDoubleMatrix （1 1 mxREAL）;
	out = （double *） mxGetPr （plhs [0]）;
	

	// do correlation over a patch
	for （i = -half_size; i <= half_size; i++）
	{
		for （j = -half_size; j <= half_size; j++）
		{
			x1i = x1 + i;
			x2i = x2 + i;

			y1i = y1 + j;
			y2i = y2 + j;
			corr += fabs（im1 [y1i * height + x1i]- im2 [y2i * height + x2i]）;
//			printf（“im1 i j x1i y1i %d %d %d %d %lf \n“ ijx1iy1iim1 [x1i * height + y1i]）;
//			printf（“im2 i j x2i y2i %d %d %d %d %lf \n“ ijx2iy2iim2 [x2i * height + y2i]）;
//			printf（“corr i j  %d %d %lf \n“ ijcorr）;

		}
		if （corr > jump_out）
		{
			out[0] = corr;
			return;
		}

	}

	out[0] = corr;

//	printf（“corr \n %lf“ corr）;

	return;
}



			//	out[0] = im1 [x1 * height + y1]- im2 [x2 * height + y2];
//	out[0] = im1 [y1 * height + x1]- im2 [y2 * height + x2];
//	out[0] = im1 [y1 * width + x1]- im2 [y2 * width + x2];

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-02-14 22:56  torr 3D reconstruction\
     文件       40960  2018-02-14 22:55  torr 3D reconstruction\birch_match.dll
     文件       20480  2018-02-14 22:55  torr 3D reconstruction\bugfixes.doc
     文件      921654  2018-02-14 22:55  torr 3D reconstruction\j1.bmp
     文件      921654  2018-02-14 22:55  torr 3D reconstruction\j2.bmp
     文件        2426  2018-02-14 22:55  torr 3D reconstruction\patch_match.c
     文件       40960  2018-02-14 22:55  torr 3D reconstruction\patch_match.dll
     文件       28672  2018-02-14 22:55  torr 3D reconstruction\readme.doc
     文件         718  2018-02-14 22:55  torr 3D reconstruction\readme.txt
     文件       37040  2018-02-14 22:55  torr 3D reconstruction\stereo.fig
     文件       53122  2018-02-14 22:55  torr 3D reconstruction\tab.cor
     文件      128540  2018-02-14 22:55  torr 3D reconstruction\table.fmatches
     文件      105588  2018-02-14 22:55  torr 3D reconstruction\table.matches
     文件      783728  2018-02-14 22:55  torr 3D reconstruction\torrsam.ps
     文件        1240  2018-02-14 22:55  torr 3D reconstruction\torr_add_manual_matches.m
     文件        2949  2018-02-14 22:55  torr 3D reconstruction\torr_birch_match.c
     文件       40960  2018-02-14 22:55  torr 3D reconstruction\torr_birch_match.dll
     文件        2408  2018-02-14 22:55  torr 3D reconstruction\torr_calc_cubic_coefs.m
     文件       17800  2018-02-14 22:55  torr 3D reconstruction\torr_cal_menu.fig
     文件        7385  2018-02-14 22:55  torr 3D reconstruction\torr_cal_menu.m
     文件        3972  2018-02-14 22:55  torr 3D reconstruction\torr_charris.m
     文件        2932  2018-02-14 22:55  torr 3D reconstruction\torr_compare_epipoles.m
     文件        4297  2018-02-14 22:55  torr 3D reconstruction\torr_compF_sc.m
     文件        5555  2018-02-14 22:55  torr 3D reconstruction\torr_comp_sam_lin_script.m
     文件        3471  2018-02-14 22:55  torr 3D reconstruction\torr_corn_matcher.m
     文件        1602  2018-02-14 22:55  torr 3D reconstruction\torr_correctx4F.m
     文件         882  2018-02-14 22:55  torr 3D reconstruction\torr_cor_script.m
     文件       13776  2018-02-14 22:55  torr 3D reconstruction\torr_c_menu.fig
     文件        7496  2018-02-14 22:55  torr 3D reconstruction\torr_c_menu.m
     文件        1498  2018-02-14 22:55  torr 3D reconstruction\torr_display_epipolar.m
     文件        2528  2018-02-14 22:55  torr 3D reconstruction\torr_display_epipoles.m
............此处省略76个文件信息