akariya

[ECE592] HW6 q1 Wcompress

Nov 24th, 2019
203
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 4.33 KB | None | 0 0
  1. clear();
  2. img = rgb2gray(imread("wolves.png"));
  3. crop = imcrop(img, [461 27 511 511]);
  4. [cr,bpp] = wcompress('c',crop,'mask.wtc','ezw','maxloop',12);
  5. Xc = wcompress('u','mask.wtc');
  6. delete('mask.wtc')
  7. colormap(pink(255))
  8. subplot(1,2,1); imshow(crop);  title('Original image')
  9. axis square
  10. subplot(1,2,2); imshow(uint8(Xc)); title('Compressed image')
  11. axis square
  12. D = double(abs(double(crop)-double(Xc)).^2);
  13. mse = sum(D(:))/numel(crop)
  14.  
  15. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  16. clear();
  17. img = rgb2gray(imread("wolves.png"));
  18. crop = imcrop(img, [461 27 511 511]);
  19. [cr,bpp] = wcompress('c',crop,'mask.wtc','spiht','maxloop',12);
  20. Xc = wcompress('u','mask.wtc');
  21. delete('mask.wtc')
  22. colormap(pink(255))
  23. subplot(1,2,1); imshow(crop);  title('Original image')
  24. axis square
  25. subplot(1,2,2); imshow(uint8(Xc)); title('Compressed image')
  26. axis square
  27. D = double(abs(double(crop)-double(Xc)).^2);
  28. mse = sum(D(:))/numel(crop)
  29.  
  30. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  31. clear();
  32. img = rgb2gray(imread("wolves.png"));
  33. crop = imcrop(img, [461 27 511 511]);
  34. [cr,bpp] = wcompress('c',crop,'mask.wtc','stw','maxloop',12);
  35. Xc = wcompress('u','mask.wtc');
  36. delete('mask.wtc')
  37. colormap(pink(255))
  38. subplot(1,2,1); imshow(crop);  title('Original image')
  39. axis square
  40. subplot(1,2,2); imshow(uint8(Xc)); title('Compressed image')
  41. axis square
  42. D = double(abs(double(crop)-double(Xc)).^2);
  43. mse = sum(D(:))/numel(crop)
  44.  
  45. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  46. clear();
  47. img = rgb2gray(imread("wolves.png"));
  48. crop = imcrop(img, [461 27 511 511]);
  49. [cr,bpp] = wcompress('c',crop,'mask.wtc','wdr','maxloop',12);
  50. Xc = wcompress('u','mask.wtc');
  51. delete('mask.wtc')
  52. colormap(pink(255))
  53. subplot(1,2,1); imshow(crop);  title('Original image')
  54. axis square
  55. subplot(1,2,2); imshow(uint8(Xc)); title('Compressed image')
  56. axis square
  57. D = double(abs(double(crop)-double(Xc)).^2);
  58. mse = sum(D(:))/numel(crop)
  59.  
  60. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  61. clear();
  62. img = rgb2gray(imread("wolves.png"));
  63. crop = imcrop(img, [461 27 511 511]);
  64. [cr,bpp] = wcompress('c',crop,'mask.wtc','aswdr','maxloop',12);
  65. Xc = wcompress('u','mask.wtc');
  66. delete('mask.wtc')
  67. colormap(pink(255))
  68. subplot(1,2,1); imshow(crop);  title('Original image')
  69. axis square
  70. subplot(1,2,2); imshow(uint8(Xc)); title('Compressed image')
  71. axis square
  72. D = double(abs(double(crop)-double(Xc)).^2);
  73. mse = sum(D(:))/numel(crop)
  74.  
  75. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  76. clear();
  77. img = rgb2gray(imread("wolves.png"));
  78. crop = imcrop(img, [461 27 511 511]);
  79. [cr,bpp] = wcompress('c',crop,'mask.wtc','spiht_3d','maxloop',12);
  80. Xc = wcompress('u','mask.wtc');
  81. delete('mask.wtc')
  82. colormap(pink(255))
  83. subplot(1,2,1); imshow(crop);  title('Original image')
  84. axis square
  85. subplot(1,2,2); imshow(uint8(Xc)); title('Compressed image')
  86. axis square
  87. D = double(abs(double(crop)-double(Xc)).^2);
  88. mse = sum(D(:))/numel(crop)
  89.  
  90. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  91. clear();
  92. img = rgb2gray(imread("wolves.png"));
  93. crop = imcrop(img, [461 27 511 511]);
  94. [cr,bpp] = wcompress('c',crop,'mask.wtc','lvl_mmc','maxloop',12);
  95. Xc = wcompress('u','mask.wtc');
  96. delete('mask.wtc')
  97. colormap(pink(255))
  98. subplot(1,2,1); imshow(crop);  title('Original image')
  99. axis square
  100. subplot(1,2,2); imshow(uint8(Xc)); title('Compressed image')
  101. axis square
  102. D = double(abs(double(crop)-double(Xc)).^2);
  103. mse = sum(D(:))/numel(crop)
  104.  
  105. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  106. clear();
  107. img = rgb2gray(imread("wolves.png"));
  108. crop = imcrop(img, [461 27 511 511]);
  109. [cr,bpp] = wcompress('c',crop,'mask.wtc','gbl_mmc_f','maxloop',12);
  110. Xc = wcompress('u','mask.wtc');
  111. delete('mask.wtc')
  112. colormap(pink(255))
  113. subplot(1,2,1); imshow(crop);  title('Original image')
  114. axis square
  115. subplot(1,2,2); imshow(uint8(Xc)); title('Compressed image')
  116. axis square
  117. D = double(abs(double(crop)-double(Xc)).^2);
  118. mse = sum(D(:))/numel(crop)
  119.  
  120. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  121. clear();
  122. img = rgb2gray(imread("wolves.png"));
  123. crop = imcrop(img, [461 27 511 511]);
  124. [cr,bpp] = wcompress('c',crop,'mask.wtc','gbl_mmc_h','maxloop',12);
  125. Xc = wcompress('u','mask.wtc');
  126. delete('mask.wtc')
  127. colormap(pink(255))
  128. subplot(1,2,1); imshow(crop);  title('Original image')
  129. axis square
  130. subplot(1,2,2); imshow(uint8(Xc)); title('Compressed image')
  131. axis square
  132. D = double(abs(double(crop)-double(Xc)).^2);
  133. mse = sum(D(:))/numel(crop)
  134.  
  135. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Advertisement
Add Comment
Please, Sign In to add comment