How to use custom PSF models

[Firestar-Reimu]

I have a SDSS image and I want to subtract the stars from a galaxy

SDSS PSF model is like this:

$${\rm PSF} = \exp{\left(-\frac{1}{2}(r/\sigma_1)^2\right)} + b \exp{\left(-\frac{1}{2}(r/\sigma_2)^2\right)} + p_0 \left(1+ \frac{1}{\beta}\left(\frac{r}{\sigma_p}\right)^2\right)^{-\beta/2}$$

https://www.sdss4.org/dr17/imaging/images/#PSFData

And I have all values of parameters.

But photutils seems to have no function to use custom PSF models.