RANSAC算法用MATLAB写的代码

用MATLAB写的RANSAC，用来进行消除误匹配的，可以提高匹配质量。

function [results options] = RANSAC（X options）

% [results options] = RANSAC（X options）
%
% DESC:
% estimate the vector of parameters Theta using RANSAC （see source [1]
% [2]）
%
% VERSION:
% 1.1.5
%
% INPUT:
%
% X                 = input data. The data id provided as a matrix that has
%                     dimesnsions 2dxN where d is the data dimensionality
%                     and N is the number of elements
%
% options           = structure containing the following fields:
%
%   sigma               = noise std
%   P_inlier            = Chi squared probability threshold for inliers
%                         （i.e. the probability that an point whose squared
%                          error is less than T_noise_squared is an inlier）
%                         （default = 0.99）
%   T_noise_squared     = Error threshold （overrides sigma）
%   epsilon             = False Alarm Rate （i.e. the probability we never
%                         pick a good minimal sample set） （default = 1e-3）
%   Ps                  = sampling probability （ 1 x size（X 2） ）
%                         （default: uniform i.e. Ps is empty）
%   ind_tabu            = logical array indicating the elements that should
%                         not be considered to construct the MSS （default
%                         is empty）
%   validateMSS_fun     = function that validates a MSS
%                         Should be in the form of:
%
%                         flag = validateMSS_foo（X s）
%
%   validateTheta_fun   = function that validates a parameter vector
%                         Should be in the form of:
%
%                         flag = validateTheta_foo（X Theta s）
%
%   est_fun             = function that estimates Theta.
%                         Should be in the form of:
%
%                         [Theta k] = estimate_foo（X s）
%
%   man_fun             = function that returns the residual error.
%                         Should be in the form of:
%
%                         [E T_noise_squared] = man_fun（Theta X）
%
%   mode                = algorithm flavour
%                         ‘RANSAC‘  -> Fischler & Bolles
%                         ‘MSAC‘    -> Torr & Zisserman
%
%
%   max_iters           = maximum number of iterations  （default = inf）
%   min_iters           = minimum number of iterations  （default = 0）
%   max_no_updates      = maximum number of iterations with no updates
%                         （default = inf）
%   fix_seed            = true to fix the seed of the random number
%                         generator so that the results on the same data
%                         set are repeatable （default = false）
%   reestimate          = true to resestimate the parameter vector using
%                         all the detected inliers
%                         （default = false）
%   verbose             = true for verbose output
%                         （default = true）
%   notify_iters        = if verbose output is on then print some
%                     

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-08-08 06:39  RANSAC\
     目录           0  2011-08-08 06:39  RANSAC\Common\
     文件        1363  2011-08-08 06:39  RANSAC\Common\chi2inv_LUT.m
     文件       33667  2010-08-19 02:40  RANSAC\Common\chi2inv_LUT.mat
     目录           0  2011-08-08 06:39  RANSAC\Common\chi2inv_mathworks\
     文件         620  2011-08-08 06:39  RANSAC\Common\chi2inv_mathworks\generate_chi2inv_LUT.m
     文件        1129  2011-08-08 06:39  RANSAC\Common\get_consensus_set.m
     文件        2348  2011-08-08 06:39  RANSAC\Common\get_consensus_set_rank.m
     文件         642  2011-08-08 06:39  RANSAC\Common\get_iter_RANSAC.m
     文件        3290  2011-08-08 06:39  RANSAC\Common\get_minimal_sample_set.m
     文件         784  2011-08-08 06:39  RANSAC\Common\get_q_RANSAC.m
     文件        1105  2011-08-08 06:39  RANSAC\Common\get_rand.m
     文件        1873  2011-08-08 06:39  RANSAC\Common\get_rand_prob.m
     文件        5330  2011-08-08 06:39  RANSAC\Common\stabilize.m
     文件        7639  2010-08-19 02:40  RANSAC\COPYING.LESSER.txt
     文件       35147  2010-08-19 02:40  RANSAC\COPYING.txt
     目录           0  2011-08-08 06:39  RANSAC\Docs\
     文件     1272258  2011-08-08 06:39  RANSAC\Docs\RANSAC4Dummies.pdf
     目录           0  2011-08-08 06:39  RANSAC\Examples\
     文件       18087  2010-08-19 02:40  RANSAC\Examples\LineData.mat
     文件        3212  2011-08-08 06:39  RANSAC\Examples\test_RANSAC_affine.m
     文件        3004  2011-08-08 06:39  RANSAC\Examples\test_RANSAC_homography.m
     文件        2884  2011-08-08 06:39  RANSAC\Examples\test_RANSAC_line_01.m
     文件        1485  2011-08-08 06:39  RANSAC\Examples\test_RANSAC_line_02.m
     文件        2622  2011-08-08 06:39  RANSAC\Examples\test_RANSAC_plane.m
     文件        3158  2011-08-08 06:39  RANSAC\Examples\test_RANSAC_RST.m
     文件          11  2010-08-19 02:40  RANSAC\Examples\VERSION.txt
     目录           0  2011-08-08 06:39  RANSAC\Models\
     目录           0  2011-08-08 06:39  RANSAC\Models\Affine\
     文件        2628  2011-08-08 06:39  RANSAC\Models\Affine\affineLS.m
     文件        2054  2011-08-08 06:39  RANSAC\Models\Affine\error_affine.m
............此处省略36个文件信息