أكثر

القطع الأثرية الوهمية بعد التحويل إلى إسقاط روبنسون

القطع الأثرية الوهمية بعد التحويل إلى إسقاط روبنسون


أنا أستخدم ggplot لإنشاء خريطة باستخدام إسقاط Robinson. في النهاية ، سأقوم بتخطيط البلدان الفردية ، لكني بحاجة إلى توفير حدود القارة كطبقة خلفية (حيث لن يتم رسم بعض البلدان). يمكنني رسم طبقتين من بياناتي (القارات + البلدان) بنجاح باستخدام واحدة منمابروجالإسقاطات المضمنة في:

ومع ذلك ، عندما أقوم بتحويل CRS الخاص بملف الشكل إلى Robinson ، أحصل على قطع أثرية غريبة في القارات. إليك مثال بسيط للعمل:

مكتبة (dplyr) مكتبة (ggplot2) (rgdal) # Get Natural Earth shapefiles 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 ") # تحميل وتحصين عالم البيانات العادي <- readOGR (" ._10 "،" ne_1_1 ) القارات. روبنسون ، المنطقة = "القارة") # القارات لرسم جميع القارات <- data_frame (id = unique (as.character (القارات.إنتظام $ id))) # هذا يعمل map. Regular <- ggplot (data = all.continents) ، aes (map_id = id)) + geom_map (الخريطة = القارات.المنتظمة ، اللون = "grey50" ، التعبئة = "الأبيض" ، الحجم = 0.4) + expand_limits (x = القارات. منتظم $ long ، y = القارات. منتظم $ خط الطول) + تنسيق_المساواة () خريطة

# يحتوي هذا على خريطة مصنوعات غريبة.روبنسون <- ggplot (data = all.continents، aes (map_id = id)) + geom_map (map = Continental.robinson، color = "grey50"، fill = "white"، size = 0.4) + expand_limits (س = القارات.روبنسون دولار طويل ، ص = قارات.روبنسون دولار لات) + منسق_المساواة () map.robinson

ما هي أفضل طريقة لتحويل ملفات الأشكال هذه إلى إسقاط روبنسون دون الحصول على كل القطع الأثرية الغريبة؟

تحديث: باتباع اقتراح iant ♦ لاستخدام أشكال الأرض الطبيعية للأرض ، كل شيء يعمل بشكل جيد الآن. أعتقد أنني كنت أحاول جعل البيانات السياسية تقوم بالكثير من العمل.

# Get Natural Earth * land * shapefiles download.file (url = "http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/110m/physical/ne_110m_land.zip"، "ne_110m_land.zip" ، "auto") unzip ("ne_110m_land.zip") file.remove ("ne_110m_land.zip") # تحميل وتحصين البيانات العادية world.land <- readOGR ("."، "ne_110m_land") land.robinson <- spTransform (world.land، CRS ("+ proj = robin")) القارات. بسيط <- تحصين (land.robinson) # القارات لرسم جميع القارات <- data_frame (id = unique (as.character (القارات.بسيط $ id ))) # كل شيء يعمل الآن. 0.4) + expand_limits (س = القارات.بسيط بالدولار طويل ، ص = القارات.بسيط دولار لات) + منسق_المساواة () خريطة ثابتة


يتيح لك النظر إلى صفحة بيانات NaturalEarth للبلدان اختيار الحدود مع البحيرات وبدونها مما يقودني إلى الاعتقاد بأن لديك تلك التي بها بحيرات (ولكن من الصعب تحديد ذلك). لما تفعله قد تكون أفضل حالًا مع الساحل أو الأرض من المجموعة المادية.


شاهد الفيديو: أكثر القطع الأثرية غموضا والتي حيرت العلماء ليومنا هذا