Project (lat, lon) to rotated planar kilometres

Equirectangular projection centred at (lat_c, lon_c), then rotated rot_deg_cw degrees clockwise (default 17.5). Returns kilometres east-of-centre (after rotation) and kilometres north-of-centre (after rotation).

Usage

morie_tps_project_xy(
  lat,
  lon,
  rot_deg_cw = .MORIE_TPS_ROT_DEG_CW,
  lat_c = .MORIE_TPS_LAT_C,
  lon_c = .MORIE_TPS_LON_C
)

Arguments

lat

Numeric vector of latitudes (WGS84 degrees).

lon

Numeric vector of longitudes (WGS84 degrees).

rot_deg_cw

Clockwise rotation in degrees (default 17.5).

lat_c, lon_c

Centre-point of the projection (default downtown Toronto: 43.70 N, 79.40 W).

Value

A named list with numeric vectors x (km east of centre) and y (km north of centre).

Details

Clockwise convention: positive rot_deg_cw rotates the map so a line that previously sloped up-right slopes less (or down-right).