基于薄板样条的图像配准MATLAB代码

基于薄板样条的图像配准MATLAB代码，已测试可以正常运行

function out = idwMvInterp（imgw map maxhw p ）
% Description:
% Fill holes using Inverse Distance Weighting interpolation
%
% Inputs:
% imgw - input image
% map - Map of the canvas with 0 indicating holes and 1 indicating pixel
% maxhw - Radius of inverse weighted interpolation
% p - power for inverse weighted interpolation
%
% Output:
% out - interpolated image
%
% Author: Fitzgerald J Archibald
% Date: 23-Apr-09

outH  = size（imgw1）;
outW  = size（imgw2）;
out = imgw;

[yi_arr xi_arr] = find（map==0）; % Find locations needing fill
if isempty（yi_arr） == false
    color = size（imgw3）;
    for ix = 1:length（yi_arr）

        xi = xi_arr（ix）;
        yi = yi_arr（ix）;
        
        % Find min window which has non-hole neighbors
        yixL=max（yi-maxhw1）;
        yixU=min（yi+maxhwoutH）;
        xixL=max（xi-maxhw1）;
        xixU=min（xi+maxhwoutW）;

        % use inverse distance weighting filter for filling
        mapw = map（yixL:yixUxixL:xixU）;
        if isempty（find（mapw 1）） == false
            wk = compWk（mapw xi-xixL+1 yi-yixL+1 p）; % compute weights for interpolation
            for colIx = 1:color
                out（yixicolIx） = idw（imgw（yixL:yixU xixL:xixU colIx） mapw wk）; % interpolation
            end
        end
    end
end

return;

% compute wk
function wk = compWk（map cx cy p）

[hw] = size（map）;
[xy] = meshgrid（1:h1:w）;
y=（y-cx）.^2;
x=（x-cy）.^2;
d2=x+y; % square of distance
wk = 1./（d2‘.^（p/2））;

return;

% Inverse distance weighting
function out = idw（in map wk）

num=sum（double（in（find（map）））.*wk（find（map）））; % weight the available pixel values among the neighbors
den=sum（wk（find（map）））; % sum of contributing weights within radius r

out = num / den;

return

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-06-27 11:22  tpsWarp\
     目录           0  2015-06-27 10:54  tpsWarp\data\
     文件      110614  2014-02-12 12:48  tpsWarp\data\0505_02.jpg
     文件     2083861  2014-02-12 12:48  tpsWarp\data\MOV03798.MPG
     目录           0  2015-06-27 10:54  tpsWarp\doc\
     文件      316851  2014-02-12 12:48  tpsWarp\doc\readme.pdf
     文件        1800  2014-02-12 12:48  tpsWarp\idwMvInterp.m
     文件        2199  2014-02-12 12:48  tpsWarp\interp2d.m
     文件        2185  2014-02-12 12:48  tpsWarp\map.mat
     文件        1400  2014-02-12 12:48  tpsWarp\nearestInterp.m
     文件         399  2014-02-12 12:48  tpsWarp\tpsDemoLandmark.mat
     文件        3893  2015-06-29 09:03  tpsWarp\tpswarp.m
     文件        2245  2015-06-27 11:02  tpsWarp\tpsWarpDemo.m