reduct.theodolite: Data reduction for theodolite measurements using the...

Description Usage Arguments References See Also Examples

View source: R/data-reduction.R

Description

This function calculates the true azimuth of a structure measured with a theodolite using the sunsight technique.

Usage

1
2
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 and longitude 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 (eg. "GMT", "CET") or a string with continent followed by country capital (eg. "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 centre 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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
lat <- ten(35,50,37.8)
lon <- ten(14,34,6.4)
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), 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), 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('UFXERSLQ')
data <- reduct.theodolite(hor, az, date, time, tz= "Europe/Malta", az.sun)

f-silva-archaeo/skyscapeR documentation built on Jan. 15, 2019, 3:22 a.m.