我是在Quentin Read's blog之后做的。我的代码如下：

df <- tigris::counties(state = NULL, cb = T, class = 'sf')

df <- df %>% filter(!STATEFP %in% c('60', '66', '69', '72', '78'))

crs_lambert <- "+proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs"

df <- df %>%
  st_transform(crs = crs_lambert)

alaska <- df %>% filter(STATE_NAME %in% 'Alaska')
alaska_g <- st_geometry(alaska)
alaska_centroid <- st_centroid(st_union(alaska_g))

rot <- function(a) matrix(c(cos(a), sin(a), -sin(a), cos(a)), 2, 2)

alaska_trans <- (alaska_g - alaska_centroid) * rot(-39 * pi/180) / 2.3 + alaska_centroid + c(1000000, -5000000)

alaska <- alaska %>% st_set_geometry(alaska_trans) %>% st_set_crs(st_crs(df))

hawaii <- df %>% filter(STATE_NAME %in% 'Hawaii')

hawaii_g <- st_geometry(hawaii)
hawaii_centroid <- st_centroid(st_union(hawaii_g))

hawaii_trans <- (hawaii_g - hawaii_centroid) * rot(-35 * pi/180) + hawaii_centroid + c(5200000, -1400000)
hawaii <- hawaii %>% st_set_geometry(hawaii_trans) %>% st_set_crs(st_crs(df))

tnc_map_final <- df %>%
  filter(!STATE_NAME %in% c('Alaska', 'Hawaii')) %>%
  rbind(alaska) %>%
  rbind(hawaii)

ggplot(tnc_map_final) + 
  geom_sf()

结果，我得到了下面的Map。