geodesic: To convert from Cartesian coordinate to Geographic...

In diegoalarc/GeodesiCL: Geometric Geodesy Functions

To convert from Cartesian coordinate to Geographic coordinate.

Description

With this function it is possible to convert from Cartesian coordinate to Geographic coordinate and obtain the decimal precision that you assign.

Usage

geodesic(a, XYZ_df, digits = 4)

Arguments

a	Selection of Ellipsoid.
XYZ_df	Sexagesimal longitude and latitude as dataframe.
digits	Number of digits are rounded to. DEFAULT: 4

Value

data.frame with the data in the following order: "Pt", "Lat", "Lon", "H".

Note

create data frame of epsg codes by epsg <- rgdal::make_EPSG()

References

https://github.com/OSGeo/PROJ & https://github.com/cran/rgdal

Examples

# Point name
Pto <- "St1"

# Cartesian data
X <- 1711591.78090565
Y <- -5060304.1659587
Z <- -3473256.69328603

# Pto, X, Y and Z as data.frame
XYZ_df <- as.data.frame(cbind(Pto, X, Y, Z))

# To know the ellipsoids and the order open the Ellipsoids in the package
# and look for it number
Ellip <- Ellipsoids
#View(Ellip)

# We choose the number 5 which is GRS80
value <- geodesic(5, XYZ_df, digits = 4)
print(value)

diegoalarc/GeodesiCL documentation built on July 24, 2022, 9:14 p.m.