% im1 = imread('cameraman.tif');
% figure,imshow(im1);
% d1 = mean2(im1);
% d11 = std2(im1);
% d1,d11;
% 
% I = imread('cameraman.tif');
% imhist(I);
% [h, cx] = imhist(I);
% h = log10(h);
% figure, plot(cx, h);
% 
% close all;
% im = imread('cameraman.tif');
% imshow(im);
% img_fig = gcf;
% figure;
% imhist(im);
% [h, b] = imhist(img);
% h = log10(h);
% figure, plot(b, h);
% figure(img_fig);
% improfile('bilinear');

% RGB = imread('autumn.tif');
% imshow(RGB);
% HSV = rgb2hsv(RGB);
% HSV(:,:,3) = imadjust(HSV(:,:,3), [0.02 0.68], [0 1], 0.7);
% RGB = hsv2rgb(HSV);
% figure, imshow(RGB);

% close all;
% I = imread('concordorthophoto.png');
% imshow(I);
% figure, imhist(I);
% I = imadjust(I, [10 200] / 255, [ ], 1);
% figure, imshow(I);
% figure, imhist(I);

% close all;
% im = imread('cameraman.tif');
% im = im2double(im);
% h1 = [-1, -1, -1, -1, 8, -1, -1, -1, -1];
% im1 = filter2(h1, im);
% figure, imshow(im), title('1');
% figure, imshow(im1), title('2');

% close all;
% im = imread('cameraman.tif');
% imshow(im);
% I = medfilt2(im);
% figure, imshow(I);
% I = medfilt2(I);
% figure, imshow(I);

% script;
% echo off all;
% close all;
% I = imread('cameraman.tif');
% figure,imshow(I);
% ID = imnoise(I, 'gaussian', 0, 0.01);
% figure, imshow(ID);
% ID = imnoise(I, 'salt & pepper', 0.05);
% figure, imshow(ID);
% ID = imnoise(I, 'speckle', 0.08);
% figure, imshow(ID);

% echo off all;
% close all;
% I = imread('cameraman.tif');
% imshow(I);
% J = imnoise(I, 'gaussian', 0, 0.005);
% figure, imshow(J);
% K = wiener(J, [5,5]);
% figure, imshow(K);

% script;
% echo off all;
% close all;
% I = imread('cameraman.tif');
% imshow(I);
% Id = im2double(I);
% h = fftshift( fft2(Id) );
% figure, imshow(mat2gray(log10(abs(h))));
% 
% Io = ifft2(fftshift(h));
% figure, imshow(mat2gray(Io));
% clear;

% echo off all;
% close all;
% I = imread('liftingbody.png');
% imshow(I);
% S = qtdecomp(I, .27);
% M = full(S);
% figure, imshow(M);
% blocks = repmat(uint8(0), size(S));
% for dim = [512 256 128 64 32 16 8 4 2 1];
%     numblocks = length(find(S == dim));
%     if (numblocks > 0)
%         values = repmat(uint8(1), [dim dim numblocks]);
%         values(2:dim, 2:dim,:) = 0;
%         blocks = qtsetblk(blocks, S, dim, values);
%     end
% end
% blocks(end, 1:end) = 1;
% blocks(1:end, end) = 1;
% figure, imshow(blocks, []);

% echo off all;
% close all;
% I = imread('pillsetc.png');
% figure,imshow(I);
% I = im2bw(I,0.5);
% I = im2double(I);
% L = bwlabel(I);
% RGB = label2rgb(L);
% figure, imshow(RGB);
% stats = regionprops(L, 'all');
% Area = [stats.Area];

% echo off all;
% close all;
% I = imread('pears.png');
% %imshow(I);
% BW1 = im2bw(I, 0.5);
% %figure, imshow(BW1);
% BW2 = bwmorph(BW1, 'erode', 9);
% %figure, imshow(BW2);
% BW1 = BW1 & BW2;
% %figure, imshow(BW1);
% L = bwlabel(BW1);
% RGB = label2rgb(L);
% figure, imshow(RGB);

im2 = imread('rizitas.bmp');
im1 = imread('11.bmp');
im3 = imread('issue.bmp');
im4 = imread('44.bmp');
x = CombineWithAlpha(im1, im2, im3, im4, 0.38);
clear;