Mam plik kształtów krajów świata, pobrany z here. mogę wykreślić ją w R stosującWykreślanie mapy świata w rzucie prostokątnym daje "punkty nieskończone"
countries <- readOGR("shp","TM_WORLD_BORDERS-0.3",encoding="UTF-8",stringsAsFactors=F)
par(mar=c(0,0,0,0),bg=rgb(0.3,0.4,1))
plot(countries,col=rgb(1,0.8,0.4))
Teraz chcę wykreślić ją w projekcji ortograficznej (Ziemia widziana z kosmosu), więc staram
countries <- spTransform(countries,CRS("+proj=ortho +lat_0=-10 +lon_0=-60"))
Grałem też z x_0 i y_0 parametry (jak stwierdzono here), ale zawsze pojawia się błąd:
non finite transformation detected:
[1] 45.08332 39.76804 Inf Inf
Erro em .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, :
failure in Polygons 3 Polygon 1 points 1
Além disso: Mensagens de aviso perdidas:
In .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, :
108 projected point(s) not finite
czasami w 3 wielokąta, czasami w 7. Skąd pochodzą te "Inf"? Muszę zmienić dowolny parametr? Chcę wykreślić mapę jak ten
ale skupione przede Ameryce Południowej. Dzięki za pomoc!
Próbujesz wyświetlać wszystkie współrzędne od -180 do 180 i od -90 do 90. Ale nie jest to możliwe przy rodzaju projekcji, której próbujesz użyć. Spróbuj najpierw przyciąć do regionu, który chcesz wyświetlić. I zdecydowanie polecam przeczytać trochę literatury na temat prognoz geograficznych. –
Zobacz także: http://en.wikipedia.org/wiki/Orthographic_projection_in_cartography#Mathematics dla zasady kciuka, gdzie uprawiać. – plannapus
Dzięki, Pascal. Chcę narysować cały pół-świat, tak jak powyższy obrazek. Pomyślałem, że będę potrzebował tylko środkowej długości/długości, a algorytm będzie przecinał w razie potrzeby. – Rodrigo