基于最小生成树的全局优化立体匹配方法

基于最小生成树的全局优化立体匹配方法，全局优化，效果非常好，跑赢局部优化方法，效果刚刚的。并且已经进行了simd优化

/*
 * ctmf.c - Constant-time median filtering
 * Copyright （C） 2006  Simon Perreault
 *
 * Reference: S. Perreault and P. Hébert “Median Filtering in Constant Time“
 * IEEE Transactions on Image Processing September 2007.
 *
 * This program has been obtained from http://nomis80.org/ctmf.html. No patent
 * covers this program although it is subject to the following license:
 *
 *  This program is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU General Public License as published by
 *  the Free Software Foundation either version 3 of the License or
 *  （at your option） any later version.
 *
 *  This program is distributed in the hope that it will be useful
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU General Public License for more details.
 *
 *  You should have received a copy of the GNU General Public License
 *  along with this program.  If not see .
 *
 * Contact:
 *  Laboratoire de vision et systèmes numériques
 *  Pavillon Adrien-Pouliot
 *  Université Laval
 *  Sainte-Foy Québec Canada
 *  G1K 7P4
 *
 *  perreaul@gel.ulaval.ca
 */

/* Standard C includes */
#include 
#include 
#include 
#include 
#define __MMX__
/* Type declarations */
#ifdef _MSC_VER
#include setsd.h>
typedef UINT8 uint8_t;
typedef UINT16 uint16_t;
typedef UINT32 uint32_t;
#pragma warning（ disable: 4799 ）
#else
#include 
#endif

/* Intrinsic declarations */
#if defined（__SSE2__） || defined（__MMX__）
#if defined（__SSE2__）
#include 
#elif defined（__MMX__）
#include 
#endif
#if defined（__GNUC__）
#include 
#elif defined（_MSC_VER）
#include 
#endif
#elif defined（__ALTIVEC__）
#include 
#elif defined（__AVX__）
#include 
#endif

/* Compiler peculiarities */
#if defined（__GNUC__）
#include 
#define inline __inline__
#define align（x） __attribute__ （（aligned （x）））
#elif defined（_MSC_VER）
#define inline __inline
#define align（x） __declspec（align（x））
#else
#define inline
#define align（x）
#endif

#ifndef MIN
#define MIN（ab） （（a） > （b） ? （b） : （a））
#endif

#ifndef MAX
#define MAX（ab） （（a） < （b） ? （b） : （a））
#endif

/**
 * This structure represents a two-tier histogram. The first tier （known as the
 * “coarse“ level） is 4 bit wide and the second tier （known as the “fine“ level）
 * is 8 bit wide. Pixels inserted in the fine level also get inserted into the
 * coarse bucket designated by the 4 MSBs of the fine bucket value.
 *
 * The structure is aligned on 16 bytes which is a prerequisite for SIMD
 * instructions. Each bucket is 16 bit wide which means that extra care must be
 * taken to prevent overflow.
 */
typedef struct align（16）
{
    uint16_t coarse[16];
    uint16_t fine[16][16];
} 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    ..A..H.     29696  2016-01-22 15:03  A Non-Local Cost Aggregation Method for Stereo Matching\.vs\qx_nonlocal_cost_aggregation\v14\.suo

     文件     168765  2012-05-09 19:53  A Non-Local Cost Aggregation Method for Stereo Matching\cones\disparity.pgm

     文件     506288  2012-05-09 18:03  A Non-Local Cost Aggregation Method for Stereo Matching\cones\left.ppm

     文件     506288  2012-05-09 18:03  A Non-Local Cost Aggregation Method for Stereo Matching\cones\right.ppm

     文件      27648  2012-05-09 19:51  A Non-Local Cost Aggregation Method for Stereo Matching\cones\Thumbs.db

     文件      16948  2016-01-21 17:35  A Non-Local Cost Aggregation Method for Stereo Matching\ctmf.c

     文件        265  2012-02-22 22:21  A Non-Local Cost Aggregation Method for Stereo Matching\ctmf.h

     文件       5862  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\cl.command.1.tlog

     文件     108754  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\CL.read.1.tlog

     文件       5190  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\CL.write.1.tlog

     文件      17490  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\ctmf.obj

     文件      46946  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\example.obj

     文件     331814  2012-05-09 18:02  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\left.ppm

     文件          2  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\link-cvtres.read.1.tlog

     文件          2  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\link-cvtres.write.1.tlog

     文件          2  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\link.1876-cvtres.read.1.tlog

     文件          2  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\link.1876-cvtres.write.1.tlog

     文件          2  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\link.1876.read.1.tlog

     文件          2  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\link.1876.write.1.tlog

     文件       3648  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\link.command.1.tlog

     文件       6712  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\link.read.1.tlog

     文件       2640  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\link.write.1.tlog

     文件        580  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\mt.command.1.tlog

     文件        542  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\mt.read.1.tlog

     文件        542  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\mt.write.1.tlog

     文件     116334  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\qx_basic.obj

     文件      88487  2016-01-22 15:07  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\qx_mst_kruskals_image.obj

     文件       5736  2016-01-22 14:52  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\qx_nonlo.7BE6E10D.tlog\CL.command.1.tlog

     文件     132244  2016-01-22 14:52  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\qx_nonlo.7BE6E10D.tlog\CL.read.1.tlog

     文件       5748  2016-01-22 14:52  A Non-Local Cost Aggregation Method for Stereo Matching\Debug\qx_nonlo.7BE6E10D.tlog\CL.write.1.tlog

............此处省略165个文件信息