A = imread('test.png');
redChannel = A(:, :, 1);
greenChannel = A(:, :, 2);
blueChannel = A(:, :, 3);
HnRed1 = imhist(redChannel)./numel(redChannel);
HnGreen1 = imhist(greenChannel)./numel(greenChannel);
HnBlue1 = imhist(blueChannel)./numel(blueChannel);
minDifference=1000000;
class='';


NumberOfimages=60;
prefix_image='coast';
fileformat='.jpg';
for num=1:NumberOfimages
  B = imread(strcat(prefix_image,num2str(num),fileformat));
  redChannel1 = B(:, :, 1);
  greenChannel1 = B(:, :, 2);
  blueChannel1 = B(:, :, 3);
  HnRed2 = imhist(redChannel1)./numel(redChannel1);
  HnGreen2 = imhist(greenChannel1)./numel(greenChannel1);
  HnBlue2 = imhist(blueChannel1)./numel(blueChannel1);
  difference=(HnRed1-HnRed2);
  d = size(difference);
  f1 = sum((HnRed1 - HnRed2).^2);
  f2=f1./d;
  %-------------------------------
  difference2=(HnGreen1-HnGreen2);
  d = size(difference2);
  f3 = sum((HnGreen1 - HnGreen2).^2);
  f4=f3./d;
  %-----------------------------------
  difference3=(HnBlue1-HnBlue2);
  d = size(difference3);
  f5 = sum((HnBlue1 - HnBlue2).^2);
  f6=f5./d;
  %------------------------------------------
  final=(f2+f4+f6)/3;
  if (final < minDifference)
  minDifference=final;
  class=strcat(prefix_image,num2str(num),fileformat)
  end
end


NumberOfimages=52;
prefix_image='forest';
fileformat='.jpg';
for num=1:NumberOfimages
  B = imread(strcat(prefix_image,num2str(num),fileformat));
  redChannel1 = B(:, :, 1);
  greenChannel1 = B(:, :, 2);
  blueChannel1 = B(:, :, 3);
  HnRed2 = imhist(redChannel1)./numel(redChannel1);
  HnGreen2 = imhist(greenChannel1)./numel(greenChannel1);
  HnBlue2 = imhist(blueChannel1)./numel(blueChannel1);
  difference=(HnRed1-HnRed2);
  d = size(difference)
  f1 = sum((HnRed1 - HnRed2).^2);
  f2=f1./d;
  %-------------------------------
  difference2=(HnGreen1-HnGreen2);
  d = size(difference2)
  f3 = sum((HnGreen1 - HnGreen2).^2);
  f4=f3./d;
  %-----------------------------------
  difference3=(HnBlue1-HnBlue2);
  d = size(difference3)
  f5 = sum((HnBlue1 - HnBlue2).^2);
  f6=f5./d;
  %------------------------------------------
  final=(f2+f4+f6)./3;
if (final < minDifference)
  minDifference=final;
  class=strcat(prefix_image,num2str(num),fileformat)
  end
end


NumberOfimages=40;
prefix_image='highway';
fileformat='.jpg';
for num=1:NumberOfimages
  B = imread(strcat(prefix_image,num2str(num),fileformat));
  redChannel1 = B(:, :, 1);
  greenChannel1 = B(:, :, 2);
  blueChannel1 = B(:, :, 3);
  HnRed2 = imhist(redChannel1)./numel(redChannel1);
  HnGreen2 = imhist(greenChannel1)./numel(greenChannel1);
  HnBlue2 = imhist(blueChannel1)./numel(blueChannel1);
  difference=(HnRed1-HnRed2);
  d = size(difference)
  f1 = sum((HnRed1 - HnRed2).^2);
  f2=f1./d;
  %-------------------------------
  difference2=(HnGreen1-HnGreen2);
  d = size(difference2)
  f3 = sum((HnGreen1 - HnGreen2).^2);
  f4=f3./d;
  %-----------------------------------
  difference3=(HnBlue1-HnBlue2);
  d = size(difference3)
  f5 = sum((HnBlue1 - HnBlue2).^2);
  f6=f5./d;
  %------------------------------------------
  final=(f2+f4+f6)./3;
  if (final < minDifference)
  minDifference=final;
  class=strcat(prefix_image,num2str(num),fileformat)
end



NumberOfimages=46;
prefix_image='insidecity';
fileformat='.jpg';
for num=1:NumberOfimages
  B = imread(strcat(prefix_image,num2str(num),fileformat));
  redChannel1 = B(:, :, 1);
  greenChannel1 = B(:, :, 2);
  blueChannel1 = B(:, :, 3);
  HnRed2 = imhist(redChannel1)./numel(redChannel1);
  HnGreen2 = imhist(greenChannel1)./numel(greenChannel1);
  HnBlue2 = imhist(blueChannel1)./numel(blueChannel1);
  difference=(HnRed1-HnRed2);
  d = size(difference)
  f1 = sum((HnRed1 - HnRed2).^2);
  f2=f1./d;
  %-------------------------------
  difference2=(HnGreen1-HnGreen2);
  d = size(difference2)
  f3 = sum((HnGreen1 - HnGreen2).^2);
  f4=f3./d;
  %-----------------------------------
  difference3=(HnBlue1-HnBlue2);
  d = size(difference3)
  f5 = sum((HnBlue1 - HnBlue2).^2);
  f6=f5./d;
  %------------------------------------------
  final=(f2+f4+f6)./3;
   if (final < minDifference)
  minDifference=final;
  class=strcat(prefix_image,num2str(num),fileformat)
end
 Disp(minDifference)
 Disp(class)
imshow('class');