R зміцнюють, викликаючи розрив багатокутників


10

У мене виникають проблеми з побудовою моїх просторових даних за допомогою ggplot2. Карта виглядає чудово, коли планується за допомогою spplot, тому я припускаю, що розрив відбувається на стадії укріплення.

Код такий:

#install the packages
library(rgdal)
library(mapproj)
library(raster)
library(rgeos)
library(ggplot2)
library(plyr)

if (!require(gpclib)) install.packages("gpclib", type="source")
gpclibPermit()

setwd("C:/Users/My Documents")

#read in laa to regional mapping
#must aggregate to higher level regions as data is provided at this higher level
laa_region_mapping <- read.csv("laa_region.csv", header = TRUE)

#read in LAA polygons
laa_polygons <- readOGR("ctyua_ew_generalised_WGS84.json", "OGRGeoJSON")

#merge by laa to add region column to polygon data
laa_polygons_with_region_data <- merge(laa_polygons, laa_region_mapping,
                                by.x = "CTYUA13NM", by.y = "LAA",
                                all.x = TRUE, all.y = TRUE)

# aggregate laa polygons by the 21 regions (aggregate by regoin_code)
region_polygons <- raster::aggregate(laa_polygons_with_region_data, "region_code")

Сукупність спрацювала, як це видно з спплотом (зауважте: я знайшов, як об’єднати регіони з цього посту в SE: Об’єднайте просторові багатокутники за кодом в R )

#plot the resulting polygons using spplot
spplot(region_polygons)

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

Але коли я закріплюю просторові дані, щоб я міг використовувати ggplot, по краях спостерігається розрив.

#fortify and merge to create the data frame ggplot will show on the map
region_polygons@data$id <- rownames(region_polygons@data)
region_polygons.points <- fortify(region_polygons, region = "id")

# plot the fortified df using ggplot
ggplot(data = region_polygons.points, aes(x= long, y = lat, group = id, fill=id)) + geom_polygon()

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

Як я можу зупинити це сльозотечу?

Я розглядав подібні відповіді на SE, але відповіді дозволяють припустити, що розривання відбувається під час злиття ( Що є причиною "розриву" багатокутників (артефактів) за допомогою R, ggplot та geom_polygon? ). Я думаю, що моє розрив відбувається на стадії укріплення, оскільки сплот перед фортифікацією виглядає добре.


спочатку вам потрібно узагальнити свій набір даних, щоб усунути проблему (програма не може обробити стільки вершин)
Mapperz

Відповіді:


15

Ви повинні використовувати group=groupв aesкартографуванні. Ось відтворюваний приклад вашої проблеми:

library("ggplot2")
library("raster")

x <- getData('GADM', country='GBR', level=2)
y <- fortify(x, region="NAME_2")
head(y)
#     long   lat order  hole piece      group       id
# 1 -2.049 57.23     1 FALSE     1 Aberdeen.1 Aberdeen
# 2 -2.049 57.23     2 FALSE     1 Aberdeen.1 Aberdeen
# 3 -2.049 57.23     3 FALSE     1 Aberdeen.1 Aberdeen
# 4 -2.050 57.23     4 FALSE     1 Aberdeen.1 Aberdeen
# 5 -2.050 57.23     5 FALSE     1 Aberdeen.1 Aberdeen
# 6 -2.050 57.23     6 FALSE     1 Aberdeen.1 Aberdeen


# wrong group aesthetic
ggplot(data=y, aes(y=lat, x=long, group=id, fill=id)) +
  geom_polygon() + 
  guides(fill=FALSE)

wronge aes-group

# fixed plot
ggplot(data=y, aes(y=lat, x=long, group=group, fill=id)) +
  geom_polygon() +
  guides(fill=FALSE)

фіксований сюжет


1
Дякую @rcs, це саме проблема. Це зараз виправлено.
annievic

4
Коротке пояснення: У наборі даних з укріпленими просторовими багатокутниками idє ідентифікатор функції та groupідентифікатор окремих кілець (островів, дірок тощо). Отже, якщо ви малюєте з idгрупою, вона малює всі шматочки вашої функції як одне кільце, отже, і «рветься», коли стрибає між островами.
Spacedman

Для нащадків я залишу ще одну пропозицію, оскільки бачу це часто ... Поширеною причиною розриву багатокутника є несортовані дані. Якщо вказати правильну groupестетику це не вирішує (що не стосується цього конкретного прикладу), спробувати, y <- y[order(y$order),]ймовірно, буде. orderКолонка створюється fortifyфункцією з цієї самої причини.
dmp
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.