medianw.m

基于＂trace transform＂的识别２Ｄ图像的Ｃ＋＋源程序

%weighted median
function mW=medianW(A,W)
[Y,I]=sort(A);
V=W(I);
U=cumsum(V);  hlf=U(end)/2;
ind=min(find(U>=hlf));
if U(ind)>hlf | ind==length(A)
       mW=Y(ind);
else mW=(Y(ind)+Y(ind+1))/2; 
end