Mercurial > hg > smallbox
annotate util/sparco utils/thumbFromOp.m @ 247:ecce33192fcc tip
Added tag ver_2.1 for changeset cef4500b936f
author | luisf <luis.figueira@eecs.qmul.ac.uk> |
---|---|
date | Wed, 31 Oct 2012 12:24:44 +0000 |
parents | 62f20b91d870 |
children |
rev | line source |
---|---|
ivan@77 | 1 function P = thumbFromOp(op,m,n,sm,sn,grayscale) |
ivan@77 | 2 % Output matrix P, of size m x n, sampled at |
ivan@77 | 3 % sm x sn upper top |
ivan@77 | 4 |
ivan@77 | 5 % Copyright 2008, Ewout van den Berg and Michael P. Friedlander |
ivan@77 | 6 % http://www.cs.ubc.ca/labs/scl/sparco |
ivan@77 | 7 % $Id: thumbFromOp.m 1040 2008-06-26 20:29:02Z ewout78 $ |
ivan@77 | 8 |
ivan@77 | 9 if nargin < 6, grayscale = 0; end |
ivan@77 | 10 |
ivan@77 | 11 info = op([],0); |
ivan@77 | 12 |
ivan@77 | 13 sm = min(info{1},sm); |
ivan@77 | 14 sn = min(info{2},sn); |
ivan@77 | 15 M = zeros(sm,sn); |
ivan@77 | 16 |
ivan@77 | 17 for i=1:sn |
ivan@77 | 18 v = zeros(info{2},1); v(i) = 1; |
ivan@77 | 19 w = real(op(v,1)); |
ivan@77 | 20 |
ivan@77 | 21 M(:,i) = w(1:sm); |
ivan@77 | 22 end |
ivan@77 | 23 |
ivan@77 | 24 mn = min(min(M)); |
ivan@77 | 25 mx = max(max(M)); |
ivan@77 | 26 M = (M - mn) / (mx-mn); |
ivan@77 | 27 |
ivan@77 | 28 idxm = floor(linspace(1,sm+1,m+1)); idxm = idxm(1:end-1); |
ivan@77 | 29 idxn = floor(linspace(1,sn+1,n+1)); idxn = idxn(1:end-1); |
ivan@77 | 30 |
ivan@77 | 31 if grayscale |
ivan@77 | 32 P = 1-M(idxm,idxn); |
ivan@77 | 33 else |
ivan@77 | 34 clrmap = hsv; |
ivan@77 | 35 M = 1 + round(M * (length(clrmap)-1)); |
ivan@77 | 36 P = zeros(m,n,3); |
ivan@77 | 37 for j1=1:m |
ivan@77 | 38 for j2=1:n |
ivan@77 | 39 P(j1,j2,:) = clrmap(M(idxm(j1),idxn(j2)),:); |
ivan@77 | 40 end |
ivan@77 | 41 end |
ivan@77 | 42 end |