added sample density compensation from Pipe and Menon

This is what we use around here.  I've verified it on a phantom for a radial pattern and it works very well.