STARbar: Calculate displayed and projected leaf areas of a 3D plant
In YplantQMC: Plant architectural analysis with Yplant and QuasiMC.

Description Usage Arguments Details Author(s) References See Also Examples

From a 3D plant, calculates the area of all leaves facing the viewing angle (from a given azimuth and altitude), or 'projected area', and the area of the leaves that is visible ('displayed area'). By default, repeats calculations from many viewing angles to estimate the hemi-spherically averaged STAR, or \overline{STAR}.

## S3 method for class 'plant3d'
STARbar(object, 
	method=c("gridtracer","exact","QuasiMC","slowquasimc"),
	integration=c("Turtlesky","Yplant","Turtle244"), 
	progressbar=TRUE, 
	returnldr=FALSE, 
	quiet=FALSE,
	npside=50, 
	silhouette=TRUE, 
	azimuth=NA, 
	altitude=NA, ...)
## S3 method for class 'plant3dlist'
STARbar(object, 
	quiet=FALSE, ...)

`object`	An object of class 'plant3d', see `constructplant`
`method`	The method to calculate the displayed area. See Details.
`integration`	The integration method to calculate the average STAR over the hemisphere. See Details.
`progressbar`	If TRUE (default), displays a graphical progress bar.
`npside`	For `method = "gridtracer"`, the number of grid cells per side for raytracing.
`returnldr`	If TRUE, returns a dataframe with results per viewing angle.
`quiet`	If TRUE, prints no progress to the screen.
`silhouette`	If TRUE, also calculates the 2D convex hull around the projected plant (see `Silhouette`).
`azimuth, altitude`	Azimuth and altitude of the viewing direction (Optional, in Radians).
`...`	Further arguments are ignored, for now.

This function calculates the displayed area (DA) and projected area (PA) of the entire plant, either over some specified viewing angle(s), or (by default) the hemispherically average values. From these averages, the \overline{STAR} is estimated (see Duursma et al. 2012).

The \overline{STAR} can also be calculated with the summary.plant3d function.

There are four methods for the calculation of the displayed and projected area from a given viewing angle:

gridtracer: A regular grid is placed over the projected plant. Displayed area is calculated from the number of intersecting grid points, and the grid size. The argument npside sets the number of grid points on the shorter side of the projected plant. Slow for large plants, but good estimates are still obtained with low values of npside.
exact: A polygon method to calculate the exact displayed area. Useful for testing, but very slow.
QuasiMC: Runs QuasiMC over the entire hemisphere - the QuasiMC model does the averaging. Relatively slow for small plants, but fast for large plants.
slowquasimc: Do not use this method (for testing only). Runs QuasiMC once per viewing angle, then averages the results. Painfully slow.

There are three integration methods. Note that if method = "QuasiMC", the integration method is ignored.

Turtlesky: Integrates over 58 points, that are placed approx. uniformly over the hemisphere.
Yplant: Integrates over 160 angles, distributed over the hemisphere as in the original Yplant. Weighing is applied over these angles to yield the average STAR.
Turtle244: As Turtlesky, but uses 244 points (slow method).

Remko Duursma

Duursma, R.A., D.S. Falster, F. Valladares, F.J. Sterck, R.W. Pearcy, C.H. Lusk, K.M. Sendall, M. Nordenstahl, N.C. Houter, B.J. Atwell, N. Kelly, J.W.G. Kelly, M. Liberloo, D.T. Tissue, B.E. Medlyn and D.S. Ellsworth. 2012. Light interception efficiency explained by two simple variables: a test using a diversity of small- to medium-sized woody plants. New Phytologist. 193:397-408.

projectplant,constructplant,Silhouette, summary.plant3d

# Get STARbar for the built-in Toona plant:
# (Settings are fast, not accurate)
STARbar(toona, method="gridtracer", npside=15)

## Not run: 
# For exact STAR, use:
STARbar(toona, method="exact")

# To produce an LDR file (as in the original Yplant), use:
clidstar <- STARbar(toona, method="gridtracer", npside=30, integration="Yplant", returnldr=T)
write.table(clidstar$ldr, "toona.LDR")

## End(Not run)