ortho.proj {GEOmap}

Orthogonal Map Projection

Description

Orthogonal Map Projection

Usage

ortho.proj(lat, lon, lon0, lat1, R)

Arguments

lat

latitude, degrees

lon

longitude, degrees

lon0

view origin longitude, degrees

lat1

view origin latitude, degrees

R

Radius of sphere, default=1

Details

Assumes spherical globe. This function is not part of the normal GEOmap plotting routines.

Value

list

x

x, coordinate in units of R

y

y, coordinate in units of R

Author(s)

Jonathan M. Lees<[e-mail]>

References

Snyder, John P., Map Projections- a working manual, USGS, Professional Paper, 1987.

See Also

GLOBE.ORTH, setPROJ, projtype

Examples

olat = 0
         olon = 0

          tlat = 23
         tlon = 30
R = 1
ortho.proj(tlat, tlon, olon, olat, R)