Daha çox

Robinson proyeksiyasına çevrildikdən sonra xəyal əsərləri

Robinson proyeksiyasına çevrildikdən sonra xəyal əsərləri


Robinson proyeksiyasından istifadə edərək bir xəritə yaratmaq üçün ggplot istifadə edirəm. Nəhayət, ayrı -ayrı ölkələri planlayacağam, amma qitə sərhədlərini bir fon qatı olaraq təqdim etməliyəm (çünki bəzi ölkələr planlanmayacaq). Verilərimin iki qatını (qitələr + ölkələr) biri ilə uğurla qura bilərəmmapprojquraşdırılmış proyeksiyalar:

Bununla birlikdə, şəkil faylının CRS -ni Robinsona çevirəndə qitələrdə qəribə əsərlər əldə edirəm. İşdə minimal bir iş nümunəsi:

kitabxana (dplyr) kitabxana (ggplot2) kitabxana (rgdal) # Təbii Yer şəklini əldə edin download.file (url = "http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/110m/cultural/ne_110m_admin_0_countries") .zip "," ne_110m_admin_0_countries.zip "," auto ") unzip (" ne_110m_admin_0_countries.zip ") file.remove (" ne_110m_admin_0_countries.zip ") # Müntəzəm məlumat dünyasını yükləyin və möhkəmləndirin <--cadmin" " ) qitələr.regular <- qala (dünya, bölgə = "qitə") # Robinsona çevir və dünyanı möhkəmləndir.robinson <- spTransform (dünya, CRS ("+proj = robin")) qitələr.robinson <- qala (dünya. robinson, region = "qitə") # Bütün qitələri qurmaq üçün qitələr <- data_frame (id = unikal (xarakter kimi (qitələr , aes (map_id = id)) + geom_map (xəritə = qitələr. nizamlı, color = "grey50", fill = "ağ", ölçüsü = 0.4) + genişləndirmə hədləri (x = qitələr. nizamlı $ uzun, y = qitələr. nizamlı $) lat) + coord_equal () map.regular

# Burada qəribə əsərlər var map.robinson <- ggplot (data = all.continents, aes (map_id = id)) + geom_map (map = continents.robinson, color = "grey50", fill = "white", size = 0.4) + expand_limits (x = qitələr.robinson $ uzunluğunda, y = qitələr.robinson $ lat) + coord_equal () xəritə.robinson

Bu qəribə əsərləri əldə etmədən bu şəkil fayllarını Robinson proyeksiyasına çevirməyin ən yaxşı yolu nədir?

YENİLƏNİB: Iant ♦ ın Natural Earth torpaq şəkillərindən istifadə etmək təklifinə əsasən, indi hər şey yaxşı işləyir. Düşünürəm ki, siyasi məlumatları çox iş görməyə məcbur etdim.

# Təbii Yer əldə et * torpaq * shapefiles download.file (url = "http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/110m/physical/ne_110m_land.zip", "ne_110m_land.zip" , "avtomatik") açın ("ne_110m_land.zip") file.remove ("ne_110m_land.zip") # Daimi məlumat dünyasını yükləyin və möhkəmləndirin world.land <- readOGR (".", "ne_110m_land") land.robinson <- spTransform (dünya.vətən, CRS ("+proj = robin")) qitələr.sadə <- qala (torpaq.robinson) # Hamısını qurmaq üçün qitələr. ))) # Hər şey indi işləyir map.fixed <- ggplot (data = all.continents, aes (map_id = id)) + geom_map (map = qitələr.sadə, color = "grey50", fill = "ağ", size = 0.4) + expand_limits (x = qitələr. Sadə $ uzun, y = qitələr. Sadə $ lat) + coord_equal () xəritə.


Ölkələr üçün NaturalEarth məlumat səhifəsinə baxaraq, göllü və gölsüz olduğunuzu düşünməyimə səbəb olan göllərlə sərhəd seçimini təklif edirik (ancaq bunu söyləmək çətindir). Nə etdiyinizə görə, fiziki dəstədən Sahil Xətti və ya Torpaqla daha yaxşı ola bilərsiniz.


Videoya baxın: Ruby on Rails by Leila Hofer