Learn R Programming

sphunif (version 1.4.3)

rhea: Rhea craters from Hirata (2016)

Description

Craters on Rhea from Hirata (2016).

Usage

rhea

Arguments

Format

A data frame with 3596 rows and 4 variables:

name

name of the crater (if named).

diameter

diameter of the crater (in km).

theta

longitude angle \(\theta \in [0, 2\pi)\) of the crater center.

phi

latitude angle \(\phi \in [-\pi/2, \pi/2]\) of the crater center.

Details

The \((\theta, \phi)\) angles are such their associated planetocentric coordinates are: $$(\cos(\phi) \cos(\theta), \cos(\phi) \sin(\theta), \sin(\phi))',$$ with \((0, 0, 1)'\) denoting the north pole.

The script performing the data preprocessing is available at rhea.R.

References

Hirata, N. (2016) Differential impact cratering of Saturn's satellites by heliocentric impactors. Journal of Geophysical Research: Planets, 121:111--117. tools:::Rd_expr_doi("10.1002/2015JE004940")

Examples

Run this code
# Load data
data("rhea")

# Add Cartesian coordinates
rhea$X <- cbind(cos(rhea$theta) * cos(rhea$phi),
                sin(rhea$theta) * cos(rhea$phi),
                sin(rhea$phi))

# Tests
unif_test(data = rhea$X[rhea$diam > 15 & rhea$diam < 20, ],
          type = c("PCvM", "PAD", "PRt"), p_value = "asymp")

Run the code above in your browser using DataLab