Advertisement
yaramohamed1

HSV

Apr 10th, 2015
244
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 4.53 KB | None | 0 0
  1. fileID = fopen('result1.txt','w');
  2. totalNumber=294;
  3. prefix_image1='test';
  4. fileformat1='.jpg';
  5. for nums=1:totalNumber
  6.  A = imread(strcat(prefix_image1,num2str(nums),fileformat1));
  7.  A = rgb2hsv(A);
  8. hsv1 = A(:, :, 1);
  9. hsv2 = A(:, :, 2);
  10. hsv3 = A(:, :, 3);
  11. Hn1 = imhist(hsv1)./numel(hsv1);
  12. Hn2 = imhist(hsv2)./numel(hsv2);
  13. Hn3 = imhist(hsv3)./numel(hsv3);
  14. minDifference=1000000;
  15. %class1=0;
  16.  
  17.  
  18. NumberOfimages=60;  
  19. prefix_image='coast';
  20. fileformat='.jpg';
  21. for num=1:NumberOfimages
  22.   B = imread(strcat(prefix_image,num2str(num),fileformat));
  23.   B = rgb2hsv(B);  
  24.   hsv11 = B(:, :, 1);
  25.   hsv22 = B(:, :, 2);
  26.   hsv33 = B(:, :, 3);
  27.   Hn11 = imhist(hsv11)./numel(hsv11);
  28.   Hn22 = imhist(hsv22)./numel(hsv22);
  29.   Hn33 = imhist(hsv33)./numel(hsv33);
  30.   difference=(Hn1-Hn11);  
  31.   d = size(difference);
  32.   f1 = sum((Hn1 - Hn11).^2);
  33.   f2=f1./d;
  34.   %-------------------------------
  35.   difference2=(Hn2-Hn22);
  36.   d = size(difference2);
  37.   f3 = sum((Hn2 - Hn22).^2);
  38.   f4=f3./d;
  39.   %-----------------------------------
  40.   difference3=(Hn3-Hn33);
  41.   d = size(difference3);
  42.   f5 = sum((Hn3 - Hn3).^2);
  43.   f6=f5./d;
  44.   %------------------------------------------
  45.   final=(f2+f4+f6)/3;
  46.   if (final <minDifference)
  47.   minDifference=final;
  48.   class1=strcat(prefix_image,num2str(num),fileformat);
  49.   className='Coast';
  50.   end
  51. end
  52.  
  53.  
  54. NumberOfimages=52;
  55. prefix_image='forest';
  56. fileformat='.jpg';
  57. for num=1:NumberOfimages
  58.   B = imread(strcat(prefix_image,num2str(num),fileformat));
  59.   B = rgb2hsv(B);  
  60.   hsv11 = B(:, :, 1);
  61.   hsv22 = B(:, :, 2);
  62.   hsv33 = B(:, :, 3);
  63.   Hn11 = imhist(hsv11)./numel(hsv11);
  64.   Hn22 = imhist(hsv22)./numel(hsv22);
  65.   Hn33 = imhist(hsv33)./numel(hsv33);
  66.   difference=(Hn1-Hn11);  
  67.   d = size(difference);
  68.   f1 = sum((Hn1 - Hn11).^2);
  69.   f2=f1./d;
  70.   %-------------------------------
  71.   difference2=(Hn2-Hn22);
  72.   d = size(difference2);
  73.   f3 = sum((Hn2 - Hn22).^2);
  74.   f4=f3./d;
  75.   %-----------------------------------
  76.   difference3=(Hn3-Hn33);
  77.   d = size(difference3);
  78.   f5 = sum((Hn3 - Hn3).^2);
  79.   f6=f5./d;
  80.   %------------------------------------------
  81.   final=(f2+f4+f6)/3;
  82.   if (final <minDifference)
  83.   minDifference=final;
  84.   class1=strcat(prefix_image,num2str(num),fileformat);
  85.   className='Forest';
  86.   end
  87. end
  88.  
  89. NumberOfimages=40;
  90. prefix_image='highway';
  91. fileformat='.jpg';
  92. for num=1:NumberOfimages
  93.   B = imread(strcat(prefix_image,num2str(num),fileformat));
  94.   B = rgb2hsv(B);  
  95.   hsv11 = B(:, :, 1);
  96.   hsv22 = B(:, :, 2);
  97.   hsv33 = B(:, :, 3);
  98.   Hn11 = imhist(hsv11)./numel(hsv11);
  99.   Hn22 = imhist(hsv22)./numel(hsv22);
  100.   Hn33 = imhist(hsv33)./numel(hsv33);
  101.   difference=(Hn1-Hn11);  
  102.   d = size(difference);
  103.   f1 = sum((Hn1 - Hn11).^2);
  104.   f2=f1./d;
  105.   %-------------------------------
  106.   difference2=(Hn2-Hn22);
  107.   d = size(difference2);
  108.   f3 = sum((Hn2 - Hn22).^2);
  109.   f4=f3./d;
  110.   %-----------------------------------
  111.   difference3=(Hn3-Hn33);
  112.   d = size(difference3);
  113.   f5 = sum((Hn3 - Hn3).^2);
  114.   f6=f5./d;
  115.   %------------------------------------------
  116.   final=(f2+f4+f6)/3;
  117.   if final <minDifference
  118.   minDifference=final;
  119.   class1=strcat(prefix_image,num2str(num),fileformat);
  120.   className='Highway';
  121.   end
  122. end
  123.  
  124.  
  125. NumberOfimages=46;
  126. prefix_image='insidecity';
  127. fileformat='.jpg';
  128. for num=1:NumberOfimages
  129.   B = imread(strcat(prefix_image,num2str(num),fileformat));
  130.   B = rgb2hsv(B);  
  131.   hsv11 = B(:, :, 1);
  132.   hsv22 = B(:, :, 2);
  133.   hsv33 = B(:, :, 3);
  134.   Hn11 = imhist(hsv11)./numel(hsv11);
  135.   Hn22 = imhist(hsv22)./numel(hsv22);
  136.   Hn33 = imhist(hsv33)./numel(hsv33);
  137.   difference=(Hn1-Hn11);  
  138.   d = size(difference);
  139.   f1 = sum((Hn1 - Hn11).^2);
  140.   f2=f1./d;
  141.   %-------------------------------
  142.   difference2=(Hn2-Hn22);
  143.   d = size(difference2);
  144.   f3 = sum((Hn2 - Hn22).^2);
  145.   f4=f3./d;
  146.   %-----------------------------------
  147.   difference3=(Hn3-Hn33);
  148.   d = size(difference3);
  149.   f5 = sum((Hn3 - Hn3).^2);
  150.   f6=f5./d;
  151.   %------------------------------------------
  152.   final=(f2+f4+f6)/3;
  153.   if final <minDifference
  154.   minDifference=final;
  155.   class1=strcat(prefix_image,num2str(num),fileformat);
  156.   className='Insidecity';  
  157.   end
  158. end
  159.    testImage=strcat(prefix_image1,num2str(nums));
  160.    fprintf('Image Name : %s',testImage);
  161.    fprintf('Class Type : %s\n',className);
  162.    fprintf('Class Image Name : %s\n',class1);
  163.    fprintf('---------------------------------\n');
  164.    fprintf(fileID,'Image Name : %s',testImage);
  165.    fprintf(fileID,'Class Type : %s\n',className);
  166.    fprintf(fileID,'Image Name : %s',class1);
  167.    fprintf(fileID,'---------------------------------\n');
  168. end
  169.  
  170. fclose(fileID);
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement