أكثر

مضلعات منطقة اللون وفقًا لتسمية المجموعة؟

مضلعات منطقة اللون وفقًا لتسمية المجموعة؟


أريد أن أقدم البيانات على خريطة لبلدنا.

لقد استخدمت R لإظهار المجموعات في بلد ما! هنا الخطوات التي اتخذتها:

setwd ('D: / r / Cluster2') القناة <- odbcConnectExcel ('cluster.xls') data <- sqlFetch (channel، 'clust9') y9 <- data.frame (inf = data $ infest، faible = data $ faible، moyen = data $ moyen، fort = data $ fort، lon = data $ Lon، lat = data $ Lat) مكتبة (حفرية) d = earth.dist (y9) km <- kmeans (d، Centers = 5) hc <- hclust (d) clust <- cutree (hc، k = 5) set.seed (123) مؤامرة (hc) y9 $ clust <- cutree (hc، k = 5) map.AL <- readOGR (dsn = " D: / r / الكتلة / shp "، layer =" ALG_boundaries ") map.df <- حصن (map.AL) ggplot (map.df) + geom_path (aes (x = long ، y = lat ، group = group) ) + geom_point (data = y9، aes (x = long، y = lat، color = factor (clust))، size = 4) + scale_color_discrete ("Cluster") + coord_fixed ()

تظهر على الخريطة نقاطًا فقط للمناطق التي لها نفس الخصائص. ومع ذلك ، ما أريد أن أعرضه هو كل المنطقة الملونة وليس مجرد نقطة.


يمكنك إنشاء عمود جديد في "data.frame" البيانات المكانية ثم تعيين تسمية المجموعة الخاصة بكل مقاطعة. بعد ذلك ، استخدم الطبقات في ggplotمضلعمعscale_fill_manualلمعالجة كل مجموعة بلون معين.

هذا مثال واحد:

# استيراد ملفات الأشكال تتطلب (rgdal) # قراءة ملف شكلي بحدود المقاطعات الجزائرية (تنزيل على: http://www.diva-gis.org/gdata) alg_provinces = readOGR (dsn = "C:…  DZA_adm" ، layer = "DZA_adm1") تتطلب (ggplot2) #Get rownames from "SpatialPolygonsDataFrame" object، slot "data" alg_provinces @ data $ id = rownames (alg_provinces @ data) #Transform object of class "SpatialPolygonsDataFrame" in "data.fpoint = حصن (alg_provinces، region = "id") تتطلب (plyr) #Join geometries (ملف shp.) معلومات مع معلومات مشتقة من جدول السمات (ملف .dbf) ، في كائن من الفئة "data.frame" alg.df = انضم (alg.points، alg_provinces @ data، by = "id") # هذا هو الجزء الذي سيتعين عليك فيه تعديل الكود الخاص بك. #Here قمت بتعيين كل مجموعة يدويًا إلى مقاطعتها المحددة. # أتمنى أن تستمتع بحقيقة أنني حاولت تقديم مثال لي يتناسب مع مثلك. الكتلة = data.frame (id = alg_provinces @ data $ id، alg_provinces @ data $ NAME_1، الكتلة = c ("0"، "1"، "2"، "1"، "3"، "2"، "0 "،" 1 "،" 0 "،" 3 "،" 3 "،" 3 "،" 1 "،" 1 "،" 3 "،" 5 "،" 1 "،" 0 "،" 3 "، "0" ، "5" ، "0" ، "3" ، "1" ، "2" ، "5" ، "1" ، "5" ، "1" ، "3" ، "2" ، "1 "،" 0 "،" 3 "،" 5 "،" 1 "،" 2 "،" 1 "،" 3 "،" 4 "،" 0 "،" 4 "،" 5 "،" 0 "، "3"، "1"، "0"، "4")) # دمج تسميات المجموعة مع البيانات المكانية alg = دمج (alg.df ، الكتلة ، بواسطة = "id") #Plot مع ggplot باستخدام الطبقات "geom_polygon "،" geom_path "وتعيين لون المجموعة المحدد باستخدام" Scale_fill_manual "ggplot (alg) + geom_polygon (aes (x = long ، y = lat ، group = NAME_1 ، fill = الكتلة)) + scale_fill_manual (القيم = c (" 0 "=" أبيض "،" 1 "=" أحمر "،" 2 "=" أصفر "،" 3 "=" أخضر "،" 4 "=" أزرق "،" 5 "=" أرجواني ")) + geom_path ( aes (x = long ، y = lat ، group = NAME_1)) + format_equal () + theme_bw () + xlab ("Longitude") + ylab ("Latitude")


الشبكات الواسعة للتنقل السكاني المتكرر في جزر ساموا وآثارها على انتقال الأمراض المعدية

تنطبق بعض قيود الوصول على البيانات التي تقوم عليها النتائج. يوجد في ساموا الأمريكية عدد قليل جدًا من السكان ، ومن المحتمل أن تسمح البيانات ذات المرجعية الجغرافية عالية الدقة بتحديد الأفراد والأسر ، وتنتهك السرية. قد تكون مجموعات البيانات التي تم إنشاؤها أثناء و / أو تحليلها أثناء الدراسة الحالية متاحة من المؤلف المقابل بناءً على طلب معقول يوفر استيفاء قيود السرية.


شاهد الفيديو: أسماء قروبات لطيفهه