tpeqd: Two point equidistant and tilted perspective projections

In mapmiscTerra: Utilities for Producing Maps

Two point equidistant and tilted perspective projections

Description

Defines map projection

Usage

tpeqd(x, offset=c(0,0), axis='enu')
tpers(x,  hKm = 100*1000, tilt = -10, offset=c(0,0), axis='enu')

Arguments

x	A SpatialPoints* object of length 2 or a matrix with two columns.
hKm	Height veiwing the Earth from
tilt	Viewing angle
offset	2 coordinates to define the origin
axis	defaults to east, north, up. 'swu' would rotateo 90 degrees

Details

A coordinate reference system is returned

Value

Caracther string representing a crs.

References

https://en.wikipedia.org/wiki/Two-point_equidistant_projection/ https://proj.org/operations/projections/tpers.html

See Also

crs,project

Examples

tpeqd(rbind(c(0,0), c(10,50)))

data('netherlands')
nldCities = terra::unwrap(nldCities)

tcrs = tpeqd(x=nldCities[1:2,])

nldT = project(nldCities, tcrs)


map.new(nldT)
text(nldT,labels=nldT$name)
scaleBar(nldT, 'topright')

mapmiscTerra documentation built on July 9, 2023, 3:07 p.m.