0001 data=u1; 0002 datais=isnan(data); 0003 cc=size(data); 0004 aa=cc(1); 0005 bb=cc(2); 0006 for nn=1:aa; 0007 for mm=2:bb-1; 0008 if datais(nn,mm)==1; 0009 data(nn,mm)=(data(nn,mm-1)+data(nn,mm+1))/2; 0010 end 0011 end 0012 end 0013 u1=data; 0014 clear data datais cc aa bb nn mm 0015 0016 data=v1; 0017 datais=isnan(data); 0018 cc=size(data); 0019 aa=cc(1); 0020 bb=cc(2); 0021 for nn=1:aa; 0022 for mm=2:bb-1; 0023 if datais(nn,mm)==1; 0024 data(nn,mm)=(data(nn,mm-1)+data(nn,mm+1))/2; 0025 end 0026 end 0027 end 0028 v1=data; 0029 clear data datais cc aa bb nn mm 0030 0031 data=w1; 0032 datais=isnan(data); 0033 cc=size(data); 0034 aa=cc(1); 0035 bb=cc(2); 0036 for nn=1:aa; 0037 for mm=2:bb-1; 0038 if datais(nn,mm)==1; 0039 data(nn,mm)=(data(nn,mm-1)+data(nn,mm+1))/2; 0040 end 0041 end 0042 end 0043 w1=data; 0044 clear data datais cc aa bb nn mm 0045 disp('done nanint')