reduct.compass: Data reduction for compass measurements

In skyscapeR: Data Analysis and Visualization for Skyscape Archaeology

Description

This function calculates the true azimuth of a structure measured with a compass.

Usage

reduct.compass(loc, mag.az, date, magdec, alt, name, ID, HWT.ID)

Arguments

loc	Location, either a skyscapeR.object or a vector containing the latitude, longitude and elevation of location, in this order.
mag.az	Array of magnetic azimuth measurements.
date	(Optional) Date of measurements as a string in the format: 'YYYY/MM/DD'. Only necessary if magdec is not given.
magdec	(Optional) Magnetic declination, if known.
alt	(Optional) Altitude, necessary for automatic declination calculation. If missing and loc is a skyscapeR.horizon object then the altitude will be automatically read from the horizon profile.
name	(Optional) Names or labels to identify each measurement.
ID	(Optional) IDs or codes to identify each measurement.
HWT.ID	(Optional) HeyWhatsThat IDs relating to a previously generated horizon profile for measurement.

See Also

mag.dec, az2dec, hor2alt

Examples

loc <- c(35,-7, 100)
mag.az <- c(89.5, 105, 109.5)
data <- reduct.compass(loc, mag.az, "2016/04/02")

# Declination will be automatically calculated if the altitude is also given:
data <- reduct.compass(loc, mag.az, "2016/04/02", alt=c(1,2,0))

# Alternatively, the altitude can be automatically retrieved from a horizon profile:
hor <- downloadHWT('HIFVTBGK')
data <- reduct.compass(hor, mag.az, "2016/04/02")

skyscapeR documentation built on Oct. 30, 2021, 1:06 a.m.