一种莫名的图像修复算法

这里介绍一种图像修复办法，基本原理是，设定2个阈值，一大一小，在阈值范围内的数据不做处理，在阈值之外的进行变分处理。
变分公式是：

同理，Wn,Ww,Ws,具体见代码；



计算取其加权平均值，即处理后的数据

Matlba代码如下：

InputData=importdata('039高度.txt');%输入数据
InputData=reshape(InputData,128,128)';
[MinInput,MaxInput]=FindRange(InputData);%数据最大最小值
highest=0.2*MaxInput;%变分范围
lowest=0.9*MinInput;
K=1000;
[MM,NN]=size(InputData);
Binaryzation=zeros(MM,NN);
OutputData=InputData;
for i=1:MM
    for j=1:NN
        if(InputData(i,j)<highest && InputData(i,j)>lowest)
            Binaryzation(i,j)=1;
        end
    end
end
for u=1:K
    for i=1:MM
        for j=1:NN
            if Binaryzation(i,j)==0
                if j>1&&j<NN
                    E=j+1;
                    W=j-1; 
                elseif j==1
                    E=j+1;
                    W=j+1;  
                else
                    E=j-1;
                    W=j-1;
                end
                if i>1&&i<MM
                    S=i+1;
                    N=i-1;
                elseif i==1
                    S=i+1;
                    N=i+1;
                else
                    S=i-1;
                    N=i-1;
                end
                Io=InputData(i,j);
                Ie=InputData(i,E);
                Iw=InputData(i,W);
                In=InputData(N,j);
                Is=InputData(S,j);

                Ine=InputData(N,E);
                Inw=InputData(N,W);
                Ise=InputData(S,E);
                Isw=InputData(S,W);

                We=1+(Ie-Io)^2+(Ine+In-Ise-Is)^2/16;                            
                We=1/sqrt(We);

                Ww=1+(Iw-Io)^2+(Inw+In-Is-Isw)^2/16;
                Ww=1/sqrt(Ww);

                Ws=1+(Is-Io)^2+(Ise+Ie-Iw-Isw)^2/16;
                Ws=1/sqrt(Ws);

                Wn=1+(In-Io)^2+(Ine+Ie-Iw-Inw)^2/16;
                Wn=1/sqrt(Wn);

                OutputData(i,j)=(Ie*We+Iw*Ww+In*Wn+Is*Ws)/(We+Ww+Ws+Wn);

            end
        end
    end
end

处理效果如下：