Description Usage Value See Also Examples
A helper function returning a CRS
object giving the OSGB36 coordinate reference.
1 | osgb36()
|
An object of class CRS
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 | osgb36()
# CRS arguments:
# +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.999601272 +x_0=400000 +y_0=-100000
# +datum=OSGB36 +units=m +no_defs +ellps=airy
# +towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894
# create a raster
r <- raster(nrows = 10, ncols = 10,
xmn = 0, xmx = 10)
# by default wgs84
projection(r)
# [1] "+proj=longlat +datum=WGS84"
# change it to OSGB36
projection(r) <- osgb36()
projection(r)
# [1] "+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.999601272 +x_0=400000\n +y_0=-100000
# +datum=OSGB36 +units=m +no_defs +ellps=airy
# +towgs84=446.448,-125.157,542.060,0.1502,0.2470,0.8421,-20.4894"
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.