Get EPSG codes or PROJ.4 codes for projections from the State Plane Coordinate System.
Usage
state_plane(state, plane_id = NULL, type = c("epsg", "proj4"))
Arguments
- state
The postal code for the state.
- plane_id
The state plane geographic zone identifier. A
NULL
value will return a projection for the entire state.- type
The type of output to return: either an EPSG code or PROJ4 string.
See also
For documentation of the underlying State Plane Coordinate System
projection data frame, see state_proj
.
Examples
if (require(USAboundariesData)) {
state_plane(state = "MA", type = "epsg")
state_plane(state = "MA", type = "proj4")
state_plane(state = "MA", plane_id = "island", type = "epsg")
state_plane(state = "MA", plane_id = "island", type = "proj4")
# Show the difference made by a state plane projection
if (require(sf)) {
va <- us_states(states = "VA", resolution = "high")
plot(st_geometry(va), graticule = TRUE)
va <- st_transform(va, state_plane("VA"))
plot(st_geometry(va), graticule = TRUE)
}
}
#> Loading required package: USAboundariesData
#> Loading required package: sf
#> Linking to GEOS 3.10.2, GDAL 3.4.1, PROJ 8.2.1; sf_use_s2() is TRUE