... Neu T khong xac dinh%trong doi so vao thi T duoc set = 0. Ca hai gradient ngo ra duoc xac dinh%trong tam [0 1]%dec la loai mat na if (ndims(f)~=3)| (size(f,3)~=3) error('Input image must ... (ndims(f)~=3)| (size(f,3)~=3) error('Input image must be RGB');end%Tinh toan vi phan theo phuong x va y cua ba anh thanh phan su dung toan tu%Mat nash=fspecial(dec);sv=sh';Rx=imfilter(double(f(:,:,1)),sh,'replicate');Ry=imfilter(double(f(:,:,1)),sv,'replicate');Gx=imfilter(double(f(:,:,2)),sh,'replicate');Gy=imfilter(double(f(:,:,2)),sv,'replicate');Bx=imfilter(double(f(:,:,3)),sh,'replicate');By=imfilter(double(f(:,:,3)),sv,'replicate');%Tinh ... Bx.*By;A=0.5*(atan(2*gxy./(gxx-gyy+eps)));G1=0.5*((gxx+gyy)+(gxx-gyy).*cos(2*A)+2*gxy.*sin(2*A));%Bay gio lap lai voi goc + pi/2. Sau do chon gia tri cuc dai tai moi diem.A = A+pi/2;G2=0.5*((gxx+gyy)+(gxx-gyy).*cos(2*A)+2*gxy.*sin(2*A));G1=G1.^0.5;G2=G2.^0.5;%Tinh VG...