Skip to contents

These functions still attempt to do their job, but will be removed in a future version.

Usage

fm_spTransform(x, ...)

# Default S3 method
fm_spTransform(x, crs0 = NULL, crs1 = NULL, passthrough = FALSE, ...)

# S3 method for class 'SpatialPoints'
fm_spTransform(x, CRSobj, passthrough = FALSE, ...)

# S3 method for class 'SpatialPointsDataFrame'
fm_spTransform(x, CRSobj, passthrough = FALSE, ...)

fm_has_PROJ6()

fm_as_sp_crs(x, ...)

fm_sp_get_crs(x)

fm_sp2segment(...)

Arguments

x

A sp::Spatial object

...

Potential additional arguments

crs0

The source sp::CRS or inla.CRS object

crs1

The target sp::CRS or inla.CRS object

passthrough

Default is FALSE. Setting to TRUE allows objects with no CRS information to be passed through without transformation.

CRSobj

The target sp::CRS or inla.CRS object

Value

A CRS object, or NULL if no valid CRS identified

Details

This function is a convenience method to workaround PROJ4/PROJ6 differences, and the lack of a crs extraction method for Spatial objects. For newer code, use fm_crs() instead, that returns crs objects, and use fm_CRS() to extract/construct/convert to old style sp::CRS objects.

Functions

  • fm_spTransform(): [Deprecated] (See fm_transform() instead) Handle transformation of various inla objects according to coordinate reference systems of sp::CRS or INLA::inla.CRS class.

  • fm_spTransform(default): The default method handles low level transformation of raw coordinates.

  • fm_has_PROJ6(): Old checker for PROJ6.

  • fm_as_sp_crs(): Wrapper for fm_CRS() sp::Spatial and sp::CRS objects.

  • fm_sp_get_crs(): Wrapper for CRS(projargs) (PROJ4) and CRS(wkt) for sp::Spatial objects.

  • fm_sp2segment(): [Deprecated] in favour of fm_as_segm()

See also

Author

Finn Lindgren finn.lindgren@gmail.com

Examples

if (fm_safe_sp()) {
  s <- sp::SpatialPoints(matrix(1:6, 3, 2), proj4string = fm_CRS("sphere"))
  fm_CRS(s)
}
#> Coordinate Reference System:
#> Deprecated Proj.4 representation: +proj=geocent +R=1 +units=m +no_defs 
#> WKT2 2019 representation:
#> GEODCRS["unknown",
#>     DATUM["unknown",
#>         ELLIPSOID["unknown",1,0,
#>             LENGTHUNIT["metre",1,
#>                 ID["EPSG",9001]]]],
#>     PRIMEM["Reference meridian",0,
#>         ANGLEUNIT["degree",0.0174532925199433,
#>             ID["EPSG",9122]]],
#>     CS[Cartesian,3],
#>         AXIS["(X)",geocentricX,
#>             ORDER[1],
#>             LENGTHUNIT["metre",1,
#>                 ID["EPSG",9001]]],
#>         AXIS["(Y)",geocentricY,
#>             ORDER[2],
#>             LENGTHUNIT["metre",1,
#>                 ID["EPSG",9001]]],
#>         AXIS["(Z)",geocentricZ,
#>             ORDER[3],
#>             LENGTHUNIT["metre",1,
#>                 ID["EPSG",9001]]]]