define_new_grid.Rd#' This function defines a new floating grid coordinate reference system in the Oblique Mercator projection for a new floating grid based on the coordinates for a primary datum.
define_new_grid(pt_0_geo, pt_0_grid = c(0, 0), geo_crs = 4326, name = NULL)
| pt_0_geo | A numeric vector of the geographic coordinates
( |
|---|---|
| pt_0_grid | A numeric vector of the geographic coordinates
( |
| geo_crs | The coordinate reference system of the geographic coordinates,
defined using any method used by |
Object of class crs, which is a list with elements input
(length-1 character) and wkt (length-1 character).
This function requires the latitude and longitude coordinates (in WGS84 decimal degrees) of the primary datum for a new floating grid.
#> Coordinate Reference System: #> User input: +proj=omerc +lat_0=37.385 +lonc=-108.51 +alpha=0 +gamma=0 +k_0=1 +x_0=500 +y_0=500 #> wkt: #> PROJCRS["Great New Grid", #> BASEGEOGCRS["Great New Grid", #> DATUM["World Geodetic System 1984", #> ELLIPSOID["WGS 84",6378137,298.257223563, #> LENGTHUNIT["metre",1]], #> ID["EPSG",6326]], #> PRIMEM["Greenwich",0, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8901]]], #> CONVERSION["Great New Grid", #> METHOD["Hotine Oblique Mercator (variant B)", #> ID["EPSG",9815]], #> PARAMETER["Latitude of projection centre",37.385, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8811]], #> PARAMETER["Longitude of projection centre",-108.51, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8812]], #> PARAMETER["Azimuth of initial line",0, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8813]], #> PARAMETER["Angle from Rectified to Skew Grid",0, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8814]], #> PARAMETER["Scale factor on initial line",1, #> SCALEUNIT["unity",1], #> ID["EPSG",8815]], #> PARAMETER["Easting at projection centre",500, #> LENGTHUNIT["metre",1], #> ID["EPSG",8816]], #> PARAMETER["Northing at projection centre",500, #> LENGTHUNIT["metre",1], #> ID["EPSG",8817]]], #> CS[Cartesian,2], #> AXIS["(E)",east, #> ORDER[1], #> LENGTHUNIT["metre",1, #> ID["EPSG",9001]]], #> AXIS["(N)",north, #> ORDER[2], #> LENGTHUNIT["metre",1, #> ID["EPSG",9001]]]]# Or define from a different coordinate reference system: define_new_grid( pt_0_geo = c(720446.8, 4140492.5), pt_0_grid = c(500, 500), geo_crs = 6341, name = "Great New Grid" )#> Coordinate Reference System: #> User input: +proj=omerc +lat_0=37.384999622921 +lonc=-108.509999731967 +alpha=0 +gamma=0 +k_0=1 +x_0=500 +y_0=500 #> wkt: #> PROJCRS["Great New Grid", #> BASEGEOGCRS["Great New Grid", #> DATUM["World Geodetic System 1984", #> ELLIPSOID["WGS 84",6378137,298.257223563, #> LENGTHUNIT["metre",1]], #> ID["EPSG",6326]], #> PRIMEM["Greenwich",0, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8901]]], #> CONVERSION["Great New Grid", #> METHOD["Hotine Oblique Mercator (variant B)", #> ID["EPSG",9815]], #> PARAMETER["Latitude of projection centre",37.384999622921, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8811]], #> PARAMETER["Longitude of projection centre",-108.509999731967, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8812]], #> PARAMETER["Azimuth of initial line",0, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8813]], #> PARAMETER["Angle from Rectified to Skew Grid",0, #> ANGLEUNIT["degree",0.0174532925199433], #> ID["EPSG",8814]], #> PARAMETER["Scale factor on initial line",1, #> SCALEUNIT["unity",1], #> ID["EPSG",8815]], #> PARAMETER["Easting at projection centre",500, #> LENGTHUNIT["metre",1], #> ID["EPSG",8816]], #> PARAMETER["Northing at projection centre",500, #> LENGTHUNIT["metre",1], #> ID["EPSG",8817]]], #> CS[Cartesian,2], #> AXIS["(E)",east, #> ORDER[1], #> LENGTHUNIT["metre",1, #> ID["EPSG",9001]]], #> AXIS["(N)",north, #> ORDER[2], #> LENGTHUNIT["metre",1, #> ID["EPSG",9001]]]]