The two binning routines made available bin _only_, they do not normalise to the number of pixels in the bins. In order to achieve this, adapt the weighted binning routine by changing line 39 to:
Ibin(bini)=sum(Isort(qi(1):qi(2)).*weightfact)./sum(weightfact);
I have not looked at a similar solution for the equal intensity binning routine yet.
Leave a Reply