clear?all;
close?all;
ima?=??double（imread（‘barbara.png‘））;
hD?=?3;
hW?=?10;
D=2*hD;
W=2*hW+1;
D1=D+1;

sigma=20;
%?h?=6*sigma;
%?h=h^2;
h?=?quantile_gauss（D1?0.88）?*?（4*?sigma^2）?.*?D1.^2;
randn（‘seed‘0）;
ima_nse=ima+sigma*randn（size（ima））;
imshow（ima_nse[]）;
[widthheight]?=?size（ima_nse）;
ima_nse2?=?padarray（ima_nse[hD?hD]‘symmetric‘）;?
ima_nse3?=?padarray（ima_nse[hW+hD?hW+hD]‘symmetric‘）;?
output?=zeros（widthheight）;
Z?=zeros（widthheight）;

k=1;
w?=?cell（W*W）;
for?m?=?-hW:hW
????for?n?=?-hW:hW
????????if（m==0?&&?n==0）?continue;?end;????%保證搜索窗比較窗不重疊

????????Sd?=?（ima_nse2-ima_nse3（hW+1+m:hW+m+width+DhW+1+n:hW+n+height+D））.^2;
????????Sd?=?cumsum