buildhistexps.m

基于matlab写的对人体的分割程序

function [fgP imHist] = buildHistExps(im,resp)
%[fgPs imHistvq] = buildHistExp(im,resp)
%Returns multiple histogram probabilities
%Builds a cumulative histogram
%Set flags inside for 16 bin histogram and regression
%Right now, regression flag is broken

DO_16 = 1;
REGRESS = 0;

numTypes = size(resp,3);

if DO_16,
    imHist = imvq16(im);
else
    imHist = imvq(im);
end

siz = size(im(:,:,1));
N = prod(siz);
[aa,bb,imHist] = unique(imHist(:));
numBin = length(aa);
numTypes = size(resp,3);
imHist = imHist(:);
resp = reshape(resp,N,numTypes);
resp = resp./repmat(sum(resp),N,1);
fgP = zeros(numBin,numTypes);
%fgP = resp
%bgP = 1-fgP;

for k = 1:numBin,
   inds = find(imHist == k);
   fg = sum(resp(inds,:));
   bg = (length(inds) - fg)/(prod(siz) - 1);
   fgP(k,:) = fg./(fg + bg);
end

%Reshape back

DEBUG = 0;
if DEBUG,
    p = 1;
    imagesc(reshape(fgP(imHist,p),siz));
    %keyboard;
end