Make it easier to specify how the PSF profile should be centered on an image about piff HOT 4 CLOSED

works for me

from piff.

from piff.

Can't the user just say (int(x), int(y)) to centre in a pixel?

Yes, that is the current workaround.

im = psf.draw(x=int(x+0.5), y=int(y+0.5))

But it feels slightly clunky compared to

im = psf.draw(x,y, center=True)

I don't understand the center=pos option.

Well, you might not want the image you are drawing onto to have the same meanings of (x,y) as the original image where the interpolation is defined. For instance, you might have the image on which you are rendering the PSF be centered at (0,0) for convenience of some measuring code. Furthermore, you might be concerned about some pixelization effects in your code, so maybe you don't always want the PSF centered in the center of the central pixel. Maybe you want to dither around by a random sub-pixel offset relative to the center. Then you would do

cenx, ceny = np.random.rand(2) - 0.5
center = galsim.PositionD(cenx,ceny)
psf.draw(x,y, image=target_image, center=center)
measure_psf(target_image, center=center)

Anyway, if you don't have a use for it yourself, you don't need to use it (natch). But in some DES work we've been doing, there was some confusion on how Piff was working here, so this is intended to make it easier for the user to be clear about what they are asking Piff to do wrt the center of the rendering.

from piff.

FYI, I'll mostly be leveraging the work I did on this GalSim PR where I added a center parameter to the drawImage command. This will encompass how the center=pos option will work. The other two options described above will just define two reasonable values for what that pos might be.

from piff.