Nothing
profitMakePointSource=function(xcen,ycen,mag=0,magzero=0,
modellist=list(sersic=list(mag=0,re=1,nser=0.5,axrat=1,ang=0)),
psf=NULL,image=matrix(0,25,25),finesample=1L, add=FALSE, plot=FALSE,
returnfine=FALSE, ...)
{
profitCheckIsPositiveInteger(finesample)
haspsfimg = !is.null(psf)
haspsfmodel = !is.null(modellist)
stopifnot(xor(haspsfimg,haspsfmodel))
hasimage = !is.null(image)
scale=10^(-0.4*(mag-magzero))
dimimg = dim(image)
if(missing(xcen)) xcen = dimimg[1]/2
if(missing(ycen)) ycen = dimimg[2]/2
pixlocs = c(1,1)
if(haspsfimg)
{
dimpsf = dim(psf)
remapgrid=expand.grid(seq(-dimpsf[1]/2,dimpsf[1]/2,len=dimpsf[1]),seq(-dimpsf[2]/2,dimpsf[2]/2,len=dimpsf[2]))
offs=(c(xcen,ycen)-(dimpsf%%2)/2)%%1
pixlocs = pixlocs + floor(c(xcen,ycen)-offs-dimpsf/2)
remapgrid[,1]=remapgrid[,1]-offs[1]
remapgrid[,2]=remapgrid[,2]-offs[2]
remappsf=matrix(profitInterp2d(remapgrid[,1],remapgrid[,2],psf)[,3],dim(psf)[1],dim(psf)[2])
remappsf=remappsf/sum(remappsf)
output=remappsf*scale
} else {
stopifnot(hasimage)
for(comp in names(modellist))
{
compmag = modellist[[comp]]$mag
stopifnot(!is.null(compmag))
modellist[[comp]]$xcen = rep(xcen,length(compmag))
modellist[[comp]]$ycen = rep(ycen,length(compmag))
}
output = profitMakeModel(modellist,dim=dimimg, finesample=finesample, returnfine = returnfine)$z*scale
}
if(add) output=profitAddMats(image,output,pixlocs)
if(plot){
magimage(output, ...)
}
return=output
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.