verticalSection: Vertical sections

In unfoldr: Stereological Unfolding for Spheroidal Particles

Description

Compute vertical section profiles of a spheroid system

Usage

verticalSection(S, d, n = c(0, 1, 0), intern = FALSE)

Arguments

S	list of spheroids, see simPoissonSystem
d	distance of the intersecting plane from the origin of the box
n	normal vector which defines the interecting vertical plane
intern	logical, FALSE (default), return all section profiles otherwise only those which have their centers inside the correspondig intersection window

Details

The function intersects a spheroid system by a plane defined by the normal vector n either equal to c(0,1,0) (default) or c(1,0,0), which is called a vertical section. Depending on the type of spheroid (either "prolate or "oblate") the returned semi-axis lengths are those corresponding to the minor semi-axis or, respectively, major semi-axis in the way these are required for unfolding.

Value

list of sizes A, shape factors S and (vertical) angles alpha of section profiles in the plane w.r.t the 'z' axis between [0,π/2].

Author(s)

M. Baaske

Examples

 
 box <- list("xrange"=c(0,5),"yrange"=c(0,5),"zrange"=c(0,5))
 
 # (exact) bivariate size-shape (isotropic) orientation distribution (spheroids)
 theta <- list("size"=list("mx"=-2.5,"my"=0.5, "sdx"=0.35,"sdy"=0.25,"rho"=0.15),
		"orientation"=list("kappa"=1))

 S <- simPoissonSystem(theta,lam=100,size="rbinorm",box=box,
  type="prolate",perfect=TRUE,pl=1)

 sp <- verticalSection(S,d=2.5,n=c(0,1,0),intern=TRUE)
 summary(sp$alpha)
 

Example output

Size/Shape:  mx=-2.500000, sdx=0.350000, my=0.500000, sdy=0.250000, rho=0.150000 
Cumulative sum of probabilities: 0.902526, 0.997042, 0.999971, 1.000000 

Spheroid simulation with `rbinorm` (perfect=1): 
Mean number: 100.000000 (exact simulation: 138.500169) 
Number of spheroids: 13914 
Set label 'N'. 


Done. Simulated 13914 objects. 
Getting plane indices: [0 2 ] 
    Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
0.003743 0.432724 0.837595 0.812773 1.168389 1.569282 

unfoldr documentation built on Sept. 25, 2021, 1:07 a.m.