new_dwt.m

it is used to find traffic

function  [Image_DWT_Fn,watermarked] = new_dwt(a);

Level=1;

% W_img  = imread('mit.jpg');  % image acquisition
% W_img  = rgb2gray(W_img);  
% for i=201:256
%     for j=158:256
%         W_img(i,j)= 0;
%     end
% end


W_img=imread('mitlogo.jpg');
W_img=rgb2gray(W_img);
W1=imresize(W_img,[64 64]);

% figure;
% imshow(W1);

for i=1:64
    for j=1:64
        if(W1(i,j)<128)
            W1(i,j)=0;
        else
            W1(i,j)=255;
        end
    end
end

% figure('Name','Watermarking image');
%      imshow(uint8(W1));

[Image_Input_DWT_Fn,Image_DWT_Fn,LL_1,LL_2,HL_1,HL_2,LH_1,LH_2,HH_1,HH_2] = DWT_Level_Fn(Level);

[m,n]=size(HL_1);

%W_img


for i=1:m
    for j=1:n
        HL_1_W(i,j)=double(int8(HL_1(i,j))+(int8(a*W_img(i,j))));
    end
end

[watermarked] = IDWT_Level_Fn(LL_1,LL_2,HL_1_W,HL_2,LH_1,LH_2,HH_1,HH_2,Level);

% figure;
% imshow(uint8(watermarked));

% [Error_Location_Fn Error_Fn Error_Value_Max_Fn Square_Error_Fn RMS_Error_Fn PSNR_Fn PSNR_Infinity_Fn] = Error_PSNR_Calculation(Image_Input_DWT_Fn,watermarked);
% Square_Error_Fn
% Error_Value_Max_Fn
% RMS_Error_Fn
% PSNR_Fn
% PSNR_dB_Fn = 10*log10(PSNR_Fn)
%