R：使用 maptools::elide 缩放阿拉斯加

Question

我正在构建一个州的 shapefile，其中阿拉斯加和夏威夷被表示为德克萨斯州南部的某个地方，以便于制作说明性地图。使用 maptools 包和来自 https://rud.is/b/2014/11/16/moving-the-earth-well-alaska-hawaii-with-r/ 的一些代码，我已经能够使用来自 TIGER 的形状来做到这一点。

但是，现在我运行遇到了麻烦，因为我想将城市添加到我的地图中。将我的形状制作成 sp 对象，然后使用 maptools::elide 以前在阿拉斯加工作得很好，但是使用比例省略在点集合上的工作方式不同，所以我的城市在错误的地方结束：

library(maptools)
library(sf)
library(tmap)
library(tidyverse)
library(tidygeocoder)

  
ak_city_sf <-
  tribble(~city_name, ~city_search_string,
          "Juneau", "Juneau, Alaska, United States",
          "Anchorage", "Anchorage, Alaska, United States",
          "Utqiagvik", "Utqiagvik, Alaska, United States",
          "Scammon Bay", "Scammon Bay, Alaska, United States") %>%
  geocode(city_search_string, method = 'osm', lat = latitude , long = longitude) %>%
  st_as_sf(coords = c("longitude","latitude"))

st_crs(ak_city_sf) <- 4326

ak_city_sf <-
  ak_city_sf %>%
  st_transform(2163)
    
ak_state_sf <-
  tigris::states(cb = T) %>%
  filter(STUSPS == "AK") %>%
  st_transform(2163)

# before transformation, everything looks fine...

tm_shape(ak_state_sf) +
  tm_borders() +
  tm_shape(ak_city_sf) +
  tm_dots(size = .1) +
  tm_text("city_name",
          size = .5)

SCALE_FACTOR <- 10000


ak_state_sf_scaled <-
  ak_state_sf %>%
  as("Spatial") %>%
  elide(scale = SCALE_FACTOR) %>%
  st_as_sf()

st_crs(ak_state_sf_scaled) <- 2163

ak_city_sf_scaled <-
  ak_city_sf %>%
  as("Spatial") %>%
  elide(scale = SCALE_FACTOR) %>%
  st_as_sf()

st_crs(ak_city_sf_scaled) <- 2163
# after scaling, things don't look so good

tm_shape(ak_state_sf_scaled) +
  tm_borders() +
  tm_shape(ak_city_sf_scaled) +
  tm_dots(size = .1) +
  tm_text("city_name",
          size = .5)

maptools::elide 似乎是执行此类操作的最佳命令（即使它迫使我转换为 sp 对象）。规模的文档对我来说意义不大。（我不认为我可以将它们组合成一个对象，因为它们是城市的点和州的多边形）。我怎样才能像缩放状态一样缩放点？

Answer 1

要缩放或旋转两个单独的几何图形以便将它们映射在一起，有必要定义一个质心，您的州几何图形和城市几何图形都将围绕该质心旋转或缩放。这具有一定的逻辑意义，因为除非您定义（隐式或显式）缩放或旋转的中心，否则缩放或旋转点集合（在本例中代表城市）不会发生。

一旦你定义了一个公共质心（在我的例子中，我只是使用了我正在转换的状态的质心），你可以使用此处显示的仿射变换： https://geocompr.robinlovelace.net/geometric-operations.html#affine-transformations.

R：使用 maptools::elide 缩放阿拉斯加

R: scaling Alaska using maptools::elide

r

gis

maptools