Open esheldon opened 8 months ago
something like this in general
def make_image(self, pars): nrow, ncol = self.obs.image.shape wcs = self.obs.jacobian.get_galsim_wcs() model = self.make_model(pars) gsim = model.drawImage(nx=ncol, ny=nrow, wcs=wcs) return gsim.array
Note for some fitters, there can be multiple bands and can be convolved by the PSF
something like this in general