Advertisement
yaramohamed1

M2

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