R/profitMakeGaussianPSF.R

In ProFit: Fit Projected 2D Profiles to Galaxy Images

profitMakeGaussianPSF=function(fwhm=3, dim=c(25,25), trim=1-pi/4, plot=FALSE, ...){
  if(length(dim)==1){dim=rep(dim,2)}
  kernel=profitMakeModel(modellist=list(sersic=list(xcen=dim[1]/2, ycen=dim[2]/2, mag=0, re=fwhm/2, nser=0.5, axrat=1, ang=0)), dim=dim)$z
  if(trim>0){
    cut=floor(trim*prod(dim)/8)*8/prod(dim)
    kernel[kernel<quantile(kernel,cut)]=0
  }
  kernel/sum(kernel)
  
  if(plot){
	  magimage(kernel, ...)
  }
  
  return=kernel
}