Calculates coordinate mappings for inla.mesh projections.

inla.mesh.map(loc,
              projection = c("default", "longlat",
                             "longsinlat", "mollweide"),
              inverse = TRUE)

## Compute sensible default map axis limits
inla.mesh.map.lim(loc = NULL,
                  projection = c("default", "longlat",
                                 "longsinlat", "mollweide"))

Arguments

loc

Coordinates to be mapped.

projection

The projection type.

inverse

If TRUE, loc are map coordinates and coordinates in the mesh domain are calculated. If FALSE, loc are coordinates in the mesh domain and the forward map projection is calculated.

Value

For inla.mesh.map.lim, a list:

xlim

X axis limits in the map domain

ylim

Y axis limits in the map domain

No attempt is made to find minimal limits for partial spherical domains.

Author

Finn Lindgren finn.lindgren@gmail.com