my_lot.m

是属于清华大学出的书 数字信号处理及电子课件光盘内容

function I=My_Lot(X,p)
[N,K]=size(X);
x=double(X);
y1=zeros(N,K);
M=8;  L=2*M;
a=[0.13*pi 0.16*pi 0.13*pi];
P=fast_LOT(a);
y2=zeros(N+L,K);
x2=zeros(N+L,K);
x2(1:M,1:K)=x(M:-1:1,1:K);
x2((M+1):(N+M),1:K)=x';
x2((N+M+1):N+L,1:K)=x(N:-1:N-M+1,1:K);
x3=zeros(N+L,K);
for i=1:N+L-1
   if rem(i,8)==0
       for j=1:K
           if rem(j,8)==0
               y2(i-8+1:i,j-8+1:j)=P'*x2(i-8+1:i+8,j-8+1:j);
               y2(i-8+p+1:i,j-8+p+1:j)=0;
           end
       end
   end
end
for i=1:N+L-1
   if rem(i,8)==0
       for j=1:K
           if rem(j,8)==0
               Q=P*y2(i-8+1:i,j-8+1:j);
               x3(i-8+1:i+8,j-8+1:j)=x3(i-8+1:i+8,j-8+1:j)+Q;
           end
       end
   end
end
x4=x3((M+1):(N+M),1:K);
x5=x4';
x5=uint8(x5);
I=x5;