reduct.theodolite: Data reduction for theodolite measurements using the...
In skyscapeR: Data Analysis and Visualization for Skyscape Archaeology

Description Usage Arguments References See Also Examples

This function calculates the true azimuth of a structure measured with a theodolite using the sun-sight technique.

reduct.theodolite(
  loc,
  az,
  date,
  time,
  tz,
  az.sun = 0,
  limb,
  alt,
  name,
  ID,
  HWT.ID
)

`loc`	Location, either a skyscapeR.object or a vector containing the latitude, longitude and elevation of location, in this order.
`az`	Array of azimuths. Use `ten` to convert to decimal point format if necessary.
`date`	Date of measurements as a string in the format: 'YYYY/MM/DD'
`time`	Time of sun-sight measurement in the format: 'HH:MM:SS'
`tz`	Timezone of input wither as a known acronym (e.g. "GMT", "CET") or a string with continent followed by country capital (e.g. "Europe/London").
`az.sun`	(Optional) Measured azimuth of the sun. Defaults to zero.
`limb`	(Optional) Measured limb of the sun. Options are left, right. If missing the center of the sun will be used for calculations.
`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.

Ruggles, C.L.N. (1999). Astronomy in Prehistoric Britain and Ireland. Yale University Press.

sunAz, ten

lat <- ten(35,50,37.8)
lon <- ten(14,34,6.4)
elev <- 100
az <- c( ten(298,24,10), ten(302,20,40))
az.sun <- ten(327,29,50)
date <- "2016/02/20"
time <- "11:07:17"

data <- reduct.theodolite(c(lat,lon,elev), az, date , time, tz= "Europe/Malta", az.sun)

# Declination will be automatically calculated if the altitude is also given:
data <- reduct.theodolite(c(lat,lon,elev), az, date , time, tz= "Europe/Malta", az.sun, alt=c(2,5))

# Alternatively, the altitude can be automatically retrieved from a horizon profile:
hor <- downloadHWT('HIFVTBGK')
data <- reduct.theodolite(hor, az, date, time, tz= "Europe/Malta", az.sun)