BRISK图像配准算法

本算法是BRISK特征提取进而特征描述完成图像配准的MATLAB代码

function Iout=affine_warp（IinMmode）
% Affine transformation function （Rotation Translation Resize）
% This function transforms a volume with a 3x3 transformation matrix 
%
% Iout=affine_warp（IinMinvmode）
%
% inputs
%   Iin: The input image
%   Minv: The （inverse） 3x3 transformation matrix
%   mode: If 0: linear interpolation and outside pixels set to nearest pixel
%            1: linear interpolation and outside pixels set to zero
%            （cubic interpolation only support by compiled mex file）
%            2: cubic interpolation and outsite pixels set to nearest pixel
%            3: cubic interpolation and outside pixels set to zero
%
% output
%   Iout: The transformed image
%
% example
%   % Read image
%   I=im2double（imread（‘lenag2.png‘））
%   % Make a transformation matrix
%   M=make_transformation_matrix（[2 3]2[1.0 1.1]）;
%   % Transform the image
%   Iout=affine_warp（IM0）
%   % Show the image
%   figure imshow（Iout）;
%
% Function is written by D.Kroon University of Twente （February 2009）
  
% Make all xy indices
[xy]=ndgrid（0:size（Iin1）-10:size（Iin2）-1）;

% Calculate center of the image
% mean= size（Iin）/2;
% Make center of the image coordinates 00
%xd=x-mean（1）; 
%yd=y-mean（2）;
xd=x;
yd=y;

% Calculate the Transformed coordinates
Tlocalx = mean（1） + M（11） * xd + M（12） *yd + M（13） * 1;
Tlocaly = mean（2） + M（21） * xd + M（22） *yd + M（23） * 1;

switch（mode）
	case 0
		Interpolation=‘bilinear‘;
		Boundary=‘replicate‘;
	case 1
		Interpolation=‘bilinear‘;
		Boundary=‘zero‘;
	case 2
		Interpolation=‘bicubic‘;
		Boundary=‘replicate‘;
	otherwise
		Interpolation=‘bicubic‘;
		Boundary=‘zero‘;
end
Iout=image_interpolation（IinTlocalxTlocalyInterpolationBoundary）;

function Iout = image_interpolation（IinTlocalxTlocalyInterpolationBoundaryImageSize）
% This function is used to transform an 2D image in a backwards way with an
% transformation image.
%
%   Iout = image_interpolation（IinTlocalxTlocalyInterpolationBoundaryImageSize）
%
% inputs
%	   Iin : 2D greyscale or color input image
%	   TlocalxTlocaly : （Backwards） Transformation images for all image pixels
%	   Interpolation:
%       ‘nearest‘    - nearest-neighbor interpolation
%       ‘bilinear‘   - bilinear interpolation
%       ‘bicubic‘    - cubic interpolation; the default method
%      Boundary:
%       ‘zero‘       - outside input image are implicilty assumed to be zero
%       ‘replicate‘  - Input array values outside the bounds of the array
%                      are assumed to equal the nearest array border value
%	（optional）
%	   ImageSize:    - Size of output image
% outputs
%  	   Iout : The transformed image
%
% Function is written by D.Kroon University of Twente （September 2010）

if（~isa（Iin‘double‘）） Iin=double（Iin）; end
if（nargin<6） ImageSize=[size（Iin1） size（Iin2）]; end
if（ndims（Iin）==2） lo=1; else lo=3; end

switch（lower（Interpolation））
    cas

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

     文件       9721  2010-09-06 13:04  QRC_BRISK\affine_warp.m

     文件      19029  2011-09-22 17:24  QRC_BRISK\brisk.cpp

     文件     109568  2011-09-22 17:24  QRC_BRISK\brisk.mexw32

     文件       2972  2011-09-22 12:14  QRC_BRISK\brisk_interface.h

     文件     424843  2014-12-24 10:43  QRC_BRISK\data\da.png

     文件     239755  2014-11-06 15:40  QRC_BRISK\data\lena1.png

     文件     211716  2014-11-06 15:40  QRC_BRISK\data\lena2.png

     文件      72037  2014-12-24 10:43  QRC_BRISK\data\xiao.png

     文件       4970  2015-01-03 21:19  QRC_BRISK\demo.asv

     文件     137216  2011-09-22 17:24  QRC_BRISK\demo.exe

     文件       5023  2015-01-03 21:48  QRC_BRISK\demo.m

     文件        363  2015-01-03 21:44  QRC_BRISK\filter.asv

     文件        381  2015-01-03 21:46  QRC_BRISK\filter.m

     文件     370688  2011-09-22 12:14  QRC_BRISK\opencv_calib3d220.dll

     文件    2010624  2011-09-22 12:14  QRC_BRISK\opencv_core220.dll

     文件     637440  2011-09-22 12:14  QRC_BRISK\opencv_features2d220.dll

     文件      22016  2011-09-22 12:14  QRC_BRISK\opencv_flann220.dll

     文件     776192  2011-09-22 12:14  QRC_BRISK\opencv_highgui220.dll

     文件    1242112  2011-09-22 12:14  QRC_BRISK\opencv_imgproc220.dll

     文件       4055  2014-11-26 15:52  QRC_BRISK\QRC_ransac.m

     文件       2681  2014-12-07 18:45  QRC_BRISK\SIFTMATCH.asv

     文件       2693  2014-12-13 18:23  QRC_BRISK\SIFTMATCH.m

     目录          0  2015-01-03 20:25  QRC_BRISK\data

     目录          0  2015-01-03 21:24  QRC_BRISK

----------- ---------  ---------- -----  ----

              6306095                    24