skyscapeR (version 1.0.0)

reduct.theodolite: Data reduction for theodolite measurements using the sun-sight method

Description

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

Usage

reduct.theodolite(
  loc,
  az,
  date,
  time,
  tz,
  az.sun = 0,
  limb,
  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.

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.

References

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

See Also

sunAz, ten

Examples

# NOT RUN {
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)
# }