mapmisc (version 1.7.0)

tpeqd: 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

Value

An object of class CRS.

Details

A coordinate reference system is returned

References

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

See Also

CRS,spTransform

Examples

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

}

# }