predict_sectorial: Predict sectorial PSF fitted model

In AlexCast/apsfs: Monte Carlo code for simulation of atmospheric Point Spread Functions and spatially-resolved Spherical Albedo Functions

Predict sectorial PSF fitted model

Description

Solves a fitted model from fit_sectorial for the PSF of the radial and azimuthal cummulative PSF at requested radius and azimuth, possibly including a third predictor variable.

Usage

predict_sectorial(
  r,
  a = NULL,
  fit,
  type = c("psf", "dpsf", "cpsf"),
  tpred = NULL
)

Arguments

r	The radial distances (km) at which the model should be evaluated.
a	The azimuthal distances (rad) at which the model should be evaluated.
fit	A model fit from fit_sectorial.
type	Type of prediction: 'psf', 'dpsf', or 'cumpsf'. See Details.
tpred	The level of the third predictor at which the model should be evaluated.

Details

If type = cumpsf, the model fit to the cumulative PSF will be evaluated at the desired positions. If type = dpsf, the density PSF is returned (1/km2). If type == psf, quadrature is used to calculate the average density PSF of the sector and the average is scaled by the area of the sector. This should equal the Monte Carlo results. Note that in this case, r will be sorted and the returned values are for the mid points of the input vector of positions, so will have a length of length(r) - 1. Default is to return the PSF.

Value

A matrix with the PSF, density PSF or cumulative PSF.

See Also

fit_sectorial, fit_annular, predict_annular

AlexCast/apsfs documentation built on Feb. 1, 2024, 9:48 p.m.