myHOGpicture.m 777 B

12345678910111213141516171819202122232425262728
  1. function im = myHOGpicture(w, bsx, bsy )
  2. % myHOGpicture(w, bsx, bxy )
  3. % Make picture of positive HOG weights.
  4. % construct a "glyph" for each orientaion
  5. bim1 = zeros(bsy, bsx);
  6. bim1(:,round(bsx/2):round(bsx/2)+1) = 1;
  7. no = 9;
  8. bim = zeros([size(bim1) no]);
  9. bim(:,:,1) = bim1;
  10. for i = 2:no,
  11. bim(:,:,i) = imrotate(bim1, -(i-1)*(180/no), 'crop');
  12. end
  13. % make pictures of positive weights bs adding up weighted glyphs
  14. s = size(w);
  15. w(w < 0) = 0;
  16. im = zeros(bsy*s(1), bsx*s(2));
  17. for i = 1:s(1),
  18. iis = (i-1)*bsy+1:i*bsy;
  19. for j = 1:s(2),
  20. jjs = (j-1)*bsx+1:j*bsx;
  21. for k = 1:no,
  22. im(iis,jjs) = im(iis,jjs) + bim(:,:,k) * w(i,j,k);
  23. end
  24. end
  25. end
  26. end