tpeqd {mapmisc}R Documentation

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

An object of class CRS.

References

http://en.wikipedia.org/wiki/Two-point_equidistant_projection http://proj4.org/projections/tpers.html

See Also

CRS,spTransform

Examples


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

data('netherlands')

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

if(require('rgdal', quietly=TRUE)) {
nldT = spTransform(nldCities, tcrs)
projection(nldT)

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

}


[Package mapmisc version 1.7.2 Index]