Title: | A simple interface to the PROJ.4 cartographic projections library |
---|---|
Description: | A simple interface to lat/long projection and datum transformation of the PROJ.4 cartographic projections library. It allows transformation of geographic coordinates from one projection and/or datum to another. |
Authors: | Simon Urbanek <[email protected]> |
Maintainer: | Simon Urbanek <[email protected]> |
License: | GPL-2 |
Version: | 1.0-14 |
Built: | 2024-11-03 05:20:46 UTC |
Source: | https://github.com/s-u/proj4 |
Projection of lat/long coordinates or its inverse.
project(xy, proj, inverse = FALSE, degrees = TRUE, silent = FALSE, ellps.default="sphere")
project(xy, proj, inverse = FALSE, degrees = TRUE, silent = FALSE, ellps.default="sphere")
xy |
input (list, matrix or 2d-array) - see details below. |
proj |
projection definition |
inverse |
if |
degrees |
if |
silent |
if set to |
ellps.default |
default ellipsoid that will be added if no datum
or ellipsoid parameter is specified in |
The input can be a list of two or more vectors (if the list contains
more than two entries, only first two entries are used and a warning
is issued), a two-dimensional matrix or array (the number of columns
or rows must be exactly two) or a vector of the length 2. For a 2x2
input the columns are taken as x
and y
. If the input
is a list then the result will be a list with the entries named
x
and y
, otherwise the result is a matrix with two
columns.
When the list form is used, inputs are recycled with a warning when necessary.
proj
specifies the target (or source if inverse)
projection. The format can be either a single (unnamed) string which
contains all parameters:
"+proj=lcc +lat_1=33 +lat_2=45 +lat_0=39 +lon_0=-96"
or an unnamed vector of complete individual parameters:
c("+proj=lcc","+lat_1=33","+lat_2=45","+lat_0=39","+lon_0=-96")
or a named vector or list that will be composed into parameters:
list(proj="lcc", lat_1=33, lat_2=45, lat_0=39, lon_0=-96)
if degrees
is TRUE
then the latitude and longitude are
expected to be in degrees, if FALSE
then in radians.
A two column matrix or list of coordinates. If the input was a list then the output will be a list, otherwise a matrix.
## this is just very simple, because we don't want to depend on ## maps package, so we can't show more useful stuff.. data(state) s <- project(state.center, "+proj=merc") plot(s, type='n', asp=1) text(s,, state.abb)
## this is just very simple, because we don't want to depend on ## maps package, so we can't show more useful stuff.. data(state) s <- project(state.center, "+proj=merc") plot(s, type='n', asp=1) text(s,, state.abb)
Transformation of geographics coordinates from one projection to another, using PROJ.4 library.
ptransform(data, src.proj, dst.proj, silent=TRUE)
ptransform(data, src.proj, dst.proj, silent=TRUE)
data |
input (list, matrix or 2d-array) - see details below. |
src.proj |
description of the source projection |
dst.proj |
description of the destination projection |
silent |
if set to |
The data can be a list of vectors (if the list contains more than
three entries, only first three entries are used and a warning is
issued), a two-dimensional matrix or array (the number of columns
or rows must be at most three). If the input is a list then the result
will be a list with the entries named x
, y
and z
,
otherwise the result is a matrix with three columns.
When the list form is used, inputs are recycled with a warning when necessary. All unspecified coordinates are set to zero.
Note that the project specification must include an ellipsoid
or a datum. Add ellps='sphere'
to obtain the same result as
older PROJ.4 versions with no datum specification.
Datum files must be installed in order to be able to perform datum
shifts (on Windows they should be located in c:\proj
) - see
PROJ.4 website for the download of datum files.
A matrix with three columns or list with projected coordinates. If the input was a list then the output will be a list, otherwise a matrix.
Simon Urbanek
data(state) sc <- cbind(state.center$x, state.center$y) ## this is essentially the same as project except ## that the default lat/long input is in radians tr <- ptransform(sc/180*pi, '+proj=latlong +ellps=sphere', '+proj=merc +ellps=sphere') ## we can compare it with the project result res <- project(sc, c(proj="merc")) ## ptransform has z coordinate which is 0 for this projection summary(tr - cbind(res, 0))
data(state) sc <- cbind(state.center$x, state.center$y) ## this is essentially the same as project except ## that the default lat/long input is in radians tr <- ptransform(sc/180*pi, '+proj=latlong +ellps=sphere', '+proj=merc +ellps=sphere') ## we can compare it with the project result res <- project(sc, c(proj="merc")) ## ptransform has z coordinate which is 0 for this projection summary(tr - cbind(res, 0))