def __init__(self, array, pixel_scale, renormalize=False, **kwargs):
<EOL>super().__init__(array=array, pixel_scale=pixel_scale)<EOL>if renormalize:<EOL><INDENT>self[:, :] = np.divide(self, np.sum(self))<EOL><DEDENT>
Class storing a 2D Point Spread Function (PSF), including its blurring kernel.

Parameters
----------
array : ndarray
    The 2d PSF blurring kernel.
renormalize : bool
    Renormalize the PSF such that he sum of kernel values total 1.0?