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