BTC.rar_btc_matlab truncati_signature

clear all; close all; clc; p = imread('c:\snehal\pic.jpg'); p=imresize(p,[128,128]); subplot(2, 2, 1:2); imshow(p); title('Original'); [row col plane]=size(p); p=double(p); % calculating global avg of each block block=64; tiilx=row-block; tilly=col-block; totalcnt=1; factorx=0; factory=0; TotalBlock=(row/block)^2; TotalBlock; Gavg = zeros(TotalBlock,3); TotalBlock=1; Gavg(TotalBlock,1)=0; for x=1:block:tiilx+1 for y=1:block:tilly+1 for i=x:(x+block-1) for j=y:(block+y-1) Gavg(TotalBlock,1)=(p(i,j,1)+Gavg(TotalBlock,1)); Gavg(TotalBlock,2)=(p(i,j,2)+Gavg(TotalBlock,2)); Gavg(TotalBlock,3)=(p(i,j,3)+Gavg(TotalBlock,3)); end end TotalBlock=TotalBlock+1 end end TotalBlock=TotalBlock-1 divideFactor=block*block; for x=1:TotalBlock for i=1:3 Gavg(x,i)=round(Gavg(x,i)/divideFactor); end end disp(Gavg) % till here TotalBlock=(row/block)^2; info=zeros(TotalBlock,6); TotalBlock=1 np=zeros(row,col,3) %% Red component largecnt=0; scnt=0; tlarge=0; ts=0; Lavg=0; Savg=0; TotalBlock=1 for x=1:block:tiilx+1 for y=1:block:tilly+1 for i=x:(x+block-1) for j=y:(block+y-1) TotalBlock; if p(i,j,1) > Gavg(TotalBlock,1) largecnt= p(i,j,1)+largecnt; tlarge=tlarge+1; np(i,j,1)=1; else scnt=p(i,j,1)+scnt; ts=ts+1; np(i,j,1)=0; end end end tlarge=tlarge-1; ts=ts-1; Lavg=largecnt/tlarge; Savg=scnt/ts; info(TotalBlock,1)=Lavg; info(TotalBlock,2)=Savg; TotalBlock=TotalBlock+1 tlarge=0 ts=0 largecnt=0; scnt=0; end end TotalBlock=1 largecnt=0; scnt=0; tlarge=0; ts=0; Lavg=0; Savg=0; % Green component for x=1:block:tiilx+1 for y=1:block:tilly+1 for i=x:(x+block-1) for j=y:(block+y-1) TotalBlock; if p(i,j,2) > Gavg(TotalBlock,2) largecnt= p(i,j,2)+largecnt; tlarge=tlarge+1; np(i,j,2)=1; else scnt=p(i,j,2)+scnt; ts=ts+1; np(i,j,2)=0; end end end tlarge=tlarge-1; ts=ts-1; Lavg=largecnt/tlarge; Savg=scnt/ts; info(TotalBlock,3)=Lavg; info(TotalBlock,4)=Savg; TotalBlock=TotalBlock+1 tlarge=0 ts=0 largecnt=0; scnt=0; end end TotalBlock=1 largecnt=0; scnt=0; tlarge=0; ts=0; Lavg=0; Savg=0; % Blue component for x=1:block:tiilx+1 for y=1:block:tilly+1 for i=x:(x+block-1) for j=y:(block+y-1) TotalBlock; if p(i,j,3) > Gavg(TotalBlock,3) largecnt= p(i,j,3)+largecnt; tlarge=tlarge+1; np(i,j,3)=1; else scnt=p(i,j,3)+scnt; ts=ts+1; np(i,j,3)=0; end end end tlarge=tlarge-1; ts=ts-1; Lavg=largecnt/tlarge; Savg=scnt/ts; info(TotalBlock,5)=Lavg; info(TotalBlock,6)=Savg; TotalBlock=TotalBlock+1 tlarge=0 ts=0 largecnt=0; scnt=0; end end TotalBlock=1 disp(info); % Final Image for x=1:block:tiilx+1 for y=1:block:tilly+1 for i=x:(x+block-1) for j=y:(block+y-1) TotalBlock; if np(i,j,1) ==1 np(i,j,1)=info(TotalBlock,1); else np(i,j,1)=info(TotalBlock,2); end if np(i,j,2) ==1 np(i,j,2)=info(TotalBlock,3); else np(i,j,2)=info(TotalBlock,4); end if np(i,j,3) ==1 np(i,j,3)=info(TotalBlock,5); else np(i,j,3)=info(TotalBlock,6); end end end TotalBlock=TotalBlock+1 end end np=uint8(np); subplot(2, 2, 3:4); imshow(np); Title('Output Image');