%
% RealImage.m
%
% Generate real images for GW and PR algorithm tests
%
% Hui Li
% 3/7/2000
%	

gamma = 1.495;

imST = imread('../image/StanfordTowerLarge.jpg','jpg');
imwrite(imST, 'ST1.jpg','jpg');

imST1 = double(imST)./256;
imST1 = imST1.^(1/gamma);

imshow(imST1);
imwrite(imST1, 'ST1GammaInv.jpg','jpg');


% construct bright spot on the image

imST2 = imST;

imST2(369:400,497:528,1) = zeros(32,32);
imST2(369:400,497:528,2) = ones(32,32).*255;
imST2(369:400,497:528,3) = ones(32,32).*250;

imwrite(imST2, 'ST2.jpg','jpg');

imST2 = double(imST2)./256;
imST2 = imST2.^(1/gamma);

imshow(imST2);
imwrite(imST2, 'ST2GammaInv.jpg','jpg');


