Hi, Harald--
If you make a projection like this:
>>> proj = pf.h.proj(0,'Density')
then proj has most of the quantities you want.
>>> proj['Density']
will give you a flat array of the density in the 2d projection.
>>> proj.keys()
['pdy', 'pdx', 'weight_field', 'px', 'py', 'Density']
here, 'pdx' and 'pdy' are the half-widths of the zones, 'px' and 'py' are the positions of the zone centers
It also will get projections that you haven't already done. For instance, proj['x-velocity'] will make a projection of the x-velocity, then return the array.
I don't know if you can stick it into BinnedProfiles, but since the arrays are 2d, I always just make my own binning loop, I find it fast enough even on pretty big datasets.
I hope that helps,
d.