GEOmap (version 2.5-11)

ortho.proj: Orthogonal Map Projection

Description

Orthogonal Map Projection

Usage

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

Value

list

x

x, coordinate in units of R

y

y, coordinate in units of R

Arguments

lat

latitude, degrees

lon

longitude, degrees

lon0

view origin longitude, degrees

lat1

view origin latitude, degrees

R

Radius of sphere, default=1

Author

Jonathan M. Lees<[e-mail]>

Details

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

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)