RDT(input,Lframe)

function RDtransform = RDT(input,Lframe)
[m,n]=size(input);
ini=1;


v=1;
while n>v 
    v=v*2;
end
n=v;

N=log2(n);
for j=1:2
 if j==2
    in=n/Lframe;
    out=2*n/Lframe;
 else
     in=1;
     out=n/Lframe;
 end
for i=in:out
    x(i,:)=input(j,ini:Lframe+ini);
    ini=ini+Lframe;
end
end

for ini=1:out
    sum(ini)=0;
for k=0:N-4
    a=2^k;
    z=(Lframe/a)-2;
    for t=1:z
        a1=a*(t-1);
        a2=a*t;
        a3=a*(t+1);
        sum(ini)=sum(ini)+(1/4)*abs(x(ini,a1)+x(ini,a3)-2*x(ini,a2));
    end
    S(ini,k+1)=(a/Lframe)*sum(ini);
end
end

for k=1:N-3
    OUTPUT(k)=0;
    for ini=1:out
        OUTPUT(k)=OUTPUT(k)+(S(ini,k))^2;
    end
    RDtransform(k)=sqrt(OUTPUT(k));
end
return