Mapa mapy Pacyfiku z wypełnionymi krajami

Question

Próbuję wykreślić mapę Pacyfiku za pomocą World2Hires w bibliotece mapproj R, ale przy próbie wypełnienia krajów jest dziwny błąd. Jak to naprawić?

library(maps) 
library(mapproj) 
library(mapdata) 
map("world2Hires", 
    xlim=c(120, 260), 
    ylim=c(-60, 40), 
    boundary=TRUE, 
    interior=TRUE, 
    fill=TRUE, 
    col="gray30", 
) 
map.axes() 

Oto wynik:

Answer 1

Problem wydaje się być z małym podzbiorem obszarów, które powodują zawijanie. Od pewnego próbnego i błędu zapisywania oryginalnego połączenia , takiego jak mapnames <- map(...), a następnie przekazywania podzbiorów tej listy do argumentu regions= w nowym wywołaniu, mogłem uniknąć owijania wypełnień. Np .:

library(maps) 
library(mapproj) 
library(mapdata) 
map("world2Hires", regions=mapnames$names[c(1:7,14:641)], 
    xlim=c(120, 260), 
    ylim=c(-60, 40), 
    boundary=TRUE, 
    interior=TRUE, 
    fill=TRUE 
) 
map.axes() 

Co do dokładniejszej lub sensownego rozwiązania, aby temu zapobiec, jestem zakłopotany. Odtwarzanie z opcją wrap= nic nie pomaga, podobnie jak w przypadku innych opcji. Na marginesie ten problem nie pojawia się przy użyciu bazy danych "world", ale pojawia się tylko dla "world2" i "world2Hires".

Answer 2

Odpowiedź @thelatemail to najlepsze i najprostsze rozwiązanie, jakie widziałem w tym problemie. Aby było to bardziej uniwersalne, lepiej jest usunąć wielokąty według nazwy. Dzieje się tak dlatego, że w zależności od limitów, które dajesz pierwsze wywołanie map(), indeksy nazw wieloboków mogą być różne.

library(maps) 
library(mapproj) 
library(mapdata) 

mapnames <- map("world2Hires", xlim=c(120, 260), ylim=c(-60, 40), 
       fill=TRUE, plot=FALSE) 

mapnames2 <- map("world2Hires", xlim=c(100, 200), ylim=c(-20, 60), 
       fill=TRUE, plot=FALSE) 

mapnames$names[10] 
[1] "Mali" 
mapnames2$names[10] 
[1] "Thailand" 

Istnieje 8 krajów, które są przecięte na południku zerowym: Wielka Brytania, Francja, Hiszpania, Algieria, Mali, Burkina Faso, Ghana, Togo. Łącząc te nazwy krajów z mapnames$names można usunąć wielokątów niezależnie od oryginalnego stopniu:

remove <- c("UK:Great Britain", "France", "Spain", "Algeria", "Mali", 
      "Burkina Faso", "Ghana", "Togo") 

map("world2Hires", regions=mapnames$names[!(mapnames$names %in% remove)], 
    xlim=c(120, 260), 
    ylim=c(-60, 40), 
    boundary=TRUE, 
    interior=TRUE, 
    fill=TRUE 
) 
map.axes() 

Można również użyć Grepl(), ale dlatego, wielokąty są nazwane heirarchically, można usunąć niektóre podrzędne wielokątów narodów w pytaniu. Na przykład mapnames$names[grepl("UK", mapnames$names)] zwraca 34 dopasowania.

Zaproponowałbym to jako edycję, ale nie mam jeszcze uprawnień.