image - stationary wavelet transform (MATLAB) -
anyone please explain being done following code. code performs image fusion using stationary wavelet transform.
%image decomposition using discrete stationary wavelet transform [a1l1,h1l1,v1l1,d1l1] = swt2(im1,1,'sym2'); [a2l1,h2l1,v2l1,d2l1] = swt2(im2,1,'sym2'); [a1l2,h1l2,v1l2,d1l2] = swt2(a1l1,1,'sym2'); [a2l2,h2l2,v2l2,d2l2] = swt2(a2l1,1,'sym2'); % fusion @ level2 afl2 = 0.5*(a1l2+a2l2); **what these equations ?** d = (abs(h1l2)-abs(h2l2))>=0; hfl2 = d.*h1l2 + (~d).*h2l2; d = (abs(v1l2)-abs(v2l2))>=0; vfl2 = d.*v1l2 + (~d).*v2l2; d = (abs(d1l2)-abs(d2l2))>=0; dfl2 = d.*d1l2 + (~d).*d2l2; % fusion @ level1 d = (abs(h1l1)-abs(h2l1))>=0; hfl1 = d.*h1l1 + (~d).*h2l1; d = (abs(v1l1)-abs(v2l1))>=0; vfl1 = d.*v1l1 + (~d).*v2l1; d = (abs(d1l1)-abs(d2l1))>=0; dfl1 = d.*d1l1 + (~d).*d2l1; % fused image afl1 = iswt2(afl2,hfl2,vfl2,dfl2,'sym2'); imf = iswt2(afl1,hfl1,vfl1,dfl1,'sym2');
here afl2, hfl2, vfl2, dfl2 @ fusion level 2 are
approximation coefficients horizontal detail coefficients vertical detail coefficients diagonal detail coefficients
also same @ next level , respective mathematical implementations according concept.
it important read concept documents once can understand implementation easily, can find info's following link, can directly move block diagram explains concept , physical implementation:
http://ijeetc.com/ijeetcadmin/upload/ijeetc_50e52508758cf.pdf
Comments
Post a Comment