Ошибочный сюжет Колорадо

Я пробую картографический пакет R. Пришлось потрудиться, чтобы найти шейп-файл штата США, который хотел бы работать с картографическими материалами — многие из них казались слишком большими и т. д. Казалось, все идет хорошо, но в штате Колорадо есть ошибки.

library(cartography)
library(sf)
library(RColorBrewer)
library(maps)
library(ggplot2)

rm(list = ls())


# USA shape file
states <- st_as_sf(map("state", plot = F, fill = TRUE))

#seems to plot correctly here
#ggplot(states) + geom_sf(aes(fill = ID))

usa <- st_transform(states, 
                    CRS("+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=37.5 +lon_0=-96"))

# still seems to plot fine
#ggplot(usa) + geom_sf(aes(fill = ID))

usa <- st_buffer(usa, dist=0)

datamap <- usa

datamap$randoVar <- sample(1:3, length(datamap$ID), replace = T)

datamap_pencil <- getPencilLayer(
  x = datamap, 
  buffer = 500,
  size = 400, 
  lefthanded = F
)

plot(st_geometry(usa), col = "white", border = "black", bg = "lightblue1")

typoLayer(
  x = datamap_pencil,
  var="randoVar",
  col = c("aquamarine4", "yellow3","#3c5cb0"),
  lwd = .7,
  legend.values.order = 1:3,
  legend.pos = "bottomleft",
  legend.title.txt = "",
  add = TRUE
)

labelLayer(x = datamap, txt = "ID", 
           cex = 0.9, halo = TRUE, r = 0.15)

Я впервые заметил, потому что, когда я попытался объединить файл данных и выполнить заливку с помощью этой функции, в Colorado появилось сообщение «Нет данных». Аналогичным образом приведенный выше код указывает на то, что гемометрия состояния или идентификатор отключены. Я недостаточно знаю ГИС, чтобы понять, почему. Мне пришлось изменить проекцию CRS, чтобы я мог буферизовать файл карты (getPencilLayer продолжал выдавать ошибку самопересечения, которая, кажется, распространена при отображении R).

введите здесь описание изображения

Любые идеи о том, что делать?


r gis cartography
person James Holland    schedule 07.02.2020    source источник
comment
По какой-то причине st_buffer превращает Колорадо в очень маленький треугольник. Вы можете увидеть это с помощью ggplot(usa[usa$ID=='colorado',]) + geom_sf(). Я не знаю, почему это происходит, но это немного изолирует проблему.   -  person Kent Johnson    schedule 07.02.2020

Ответы (2)


arrow_upward
1
arrow_downward

Ну, в итоге я исправил с помощью шейп-файла из переписи населения США.

https://www2.census.gov/geo/tiger/TIGER2017/STATE/< /а>

states <- st_read("#mypath#/tl_2017_us_state/tl_2017_us_state.shp")
states <- states[!(states$NAME %in% 
                     c("Commonwealth of the Northern Mariana Islands", "United States Virgin Islands", 
                       "Puerto Rico", "American Samoa", "Hawaii", "Guam", "Alaska")
), ]

usa <- st_transform(states, 
                     CRS("+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=37.5 +lon_0=-96"))

usa <- st_buffer(usa, dist=0)

Не уверен, как настроить получение геодаты с карты ("состояние"...), но это сработало для меня.

Использовал эти настройки буфера и размера (позже после слияния данных)

datamap_pencil <- getPencilLayer(
  x = datamap, 
  buffer = 500,
  size = 400, 
  lefthanded = F
)
person James Holland    schedule 07.02.2020

arrow_upward
0
arrow_downward

Уже ответили здесь: https://gis.stackexchange.com/a/351910/142200

Проблема в том, что исходный объект карты недействителен

library(cartography)
library(sf)
library(RColorBrewer)
library(maps)
library(ggplot2)

rm(list = ls())
# USA shape file
states <- st_as_sf(map("state", plot = F, fill = TRUE))

usa <- st_transform(states,
                    "+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=37.5 +lon_0=-96")


datamap <- usa
# Check validity----
st_is_valid(datamap)
#>  [1]  TRUE  TRUE  TRUE FALSE FALSE  TRUE  TRUE  TRUE  TRUE  TRUE  TRUE  TRUE
#> [13]  TRUE  TRUE  TRUE FALSE  TRUE  TRUE  TRUE  TRUE  TRUE  TRUE  TRUE  TRUE
#> [25]  TRUE  TRUE  TRUE  TRUE  TRUE  TRUE FALSE  TRUE  TRUE  TRUE  TRUE  TRUE
#> [37]  TRUE  TRUE  TRUE FALSE FALSE FALSE  TRUE  TRUE FALSE  TRUE  TRUE  TRUE
#> [49]  TRUE

#Make valid
library(lwgeom)
datamap<-st_make_valid(datamap)
st_is_valid(datamap)
#>  [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
#> [16] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
#> [31] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
#> [46] TRUE TRUE TRUE TRUE

# Start cartography

datamap$randoVar <- sample(1:3, length(datamap$ID), replace = T)

datamap_pencil <- getPencilLayer(
  x = datamap, 
  buffer = 500,
  size = 400, 
  lefthanded = F
)
plot(st_geometry(usa), col = "white", border = "black", bg = "lightblue1")
typoLayer(
  x = datamap_pencil,
  var="randoVar",
  col = c("aquamarine4", "yellow3","#3c5cb0"),
  lwd = .7,
  legend.values.order = 1:3,
  legend.pos = "bottomleft",
  legend.title.txt = "",
  add = TRUE
)

введите здесь описание изображения

Создано 25 февраля 2020 г. с помощью пакета reprex (v0.3.0)

person dieghernan    schedule 25.02.2020