将经度/纬度坐标转换为兰伯特等角圆锥投影

Question

我不是地图和坐标系方面的专家。我需要将经度和纬度坐标转换为 LCC（兰伯特保角圆锥投影）。我有一个需要在地图中绘制的城市坐标列表。问题是地图的投影是"+proj=lcc +lat_1=43 +lat_2=62 +lat_0=30 +lon_0=10 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs"。

如何将经纬度坐标转换为地图的投影？

在带有城市和坐标的tibble 下方：

cities <- tibble(city.name = c('Lisbon', 'Barcelona', 'Brussels'),
                 lon = c(-9.139337,2.173404,4.351710),
                 lat = c(38.72225,41.38506,50.85034))

如何根据地图的投影将这些坐标转换成坐标？

"+proj=lcc +lat_1=43 +lat_2=62 +lat_0=30 +lon_0=10 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs"

谢谢！

Answer 1

您可以使用sf 包，其中包含sf::st_transform()。这会将坐标投影到 sf 对象中。

library(tidyverse)
library(sf)

cities <- tibble(city.name = c('Lisbon', 'Barcelona', 'Brussels'),
                 lon = c(-9.139337,2.173404,4.351710),
                 lat = c(38.72225,41.38506,50.85034))

city_proj <- "+proj=lcc +lat_1=43 +lat_2=62 +lat_0=30 +lon_0=10 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs"

cities <- st_as_sf(cities, coords = c(2, 3))
st_crs(cities) <- 4326
cities <- st_transform(cities, crs = city_proj)

ggplot(cities) + 
  geom_sf()

编辑：将初始 CRS 更改为标准 CRS。这使得转换更有意义。