这里介绍一种图像修复办法,基本原理是,设定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
处理效果如下: