मैं ggplotवेल्स में प्रशासनिक क्षेत्रों में जनसंख्या के कोरोप्लेथ मानचित्र के निर्माण की प्रक्रिया में हूं । मैंने आयुध सर्वेक्षण से बाउंड्री-लाइन डेटा डाउनलोड किया है और सही शेपफाइल (समुदाय_वर्ड_ग्रीन.शप) प्रतीत होता है। आर का उपयोग करते हुए, मुझे शेपफाइल में पढ़ने के रूप में मिला है।
require(maptools)
shape <- readShapePoly(wards)
str(shape)जो मुझे यह आशाजनक आउटपुट देता है:
Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots
..@ data :'data.frame': 1690 obs. of 4 variables:
.. ..$ NAME : Factor w/ 1507 levels "Abbey Cwmhir",..: 969 90 111 200 441 477 1455 249 255 305 ...
.. ..$ DESCRIPTIO: Factor w/ 4 levels "COMMUNITY","COMMUNITY WARD",..: 2 2 1 1 2 2 2 2 1 1 ...
.. ..$ COMMUNITY : Factor w/ 858 levels "Abbey Cwmhir",..: 67 67 81 128 152 152 152 152 157 190 ...
.. ..$ FILE_NAME : Factor w/ 23 levels "ABERTAWE_-_SWANSEA",..: 1 1 1 1 1 1 1 1 1 1 ...
.. ..- attr(*, "data_types")= chr [1:4] "C" "C" "C" "C"
..@ polygons :List of 1690
.. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots
.. .. .. ..@ Polygons :List of 1
.. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots
.. .. .. .. .. .. ..@ labpt : num [1:2] 259009 188524
.. .. .. .. .. .. ..@ area : num 1923892
.. .. .. .. .. .. ..@ hole : logi FALSE
.. .. .. .. .. .. ..@ ringDir: int 1
.. .. .. .. .. .. ..@ coords : num [1:1629, 1:2] 259413 259420 259427 259427 259432 ...
.. .. .. ..@ plotOrder: int 1
.. .. .. ..@ labpt : num [1:2] 259009 188524
.. .. .. ..@ ID : chr "0"
.. .. .. ..@ area : num 1923892अब अगर मैं ऐसा करता हूं:
bar <- fortify(shape, region = "NAME")मुझे एक अच्छा डेटा फ्रेम मिलता है, जो barमुझे उम्मीद के मुताबिक बहुत अच्छा लगता है:
> str(bar)
'data.frame': 4744053 obs. of 7 variables:
$ long : num 302962 302970 302974 303013 303015 ...
$ lat : num 280066 280076 280078 280097 280105 ...
$ order: int 1 2 3 4 5 6 7 8 9 10 ...
$ hole : logi FALSE FALSE FALSE FALSE FALSE FALSE ...
$ piece: Factor w/ 29 levels "1","2","3","4",..: 1 1 1 1 1 1 1 1 1 1 ...
$ group: Factor w/ 1762 levels "Abbey Cwmhir.1",..: 1 1 1 1 1 1 1 1 1 1 ...
$ id : chr "Abbey Cwmhir" "Abbey Cwmhir" "Abbey Cwmhir" "Abbey Cwmhir" ..हालांकि, यह एक बड़ा डेटा फ्रेम है और ggplotइसे प्रदर्शित करने की कोशिश करते समय कश से निकल जाता है। वास्तव में मैं केवल एक समय में एक क्षेत्र को देखना चाहता हूं। ऐसा लगता FILE_NAMEहै कि आकार वस्तु का कारक वह है जो मैं चाहता हूं क्योंकि यह ज्यादातर काउंटियों और प्रमुख सम्मेलनों से मेल खाती है।
> unique(shape@data$FILE_NAME)
[1] ABERTAWE_-_SWANSEA
[2] BLAENAU_GWENT_-_BLAENAU_GWENT
[3] BRO_MORGANNWG_-_THE_VALE_OF_GLAMORGAN
[4] CAERDYDD_-_CARDIFF
[5] CAERFFILI_-_CAERPHILLY
[6] CASNEWYDD_-_NEWPORT
[7] CASTELL-NEDD_PORT_TALBOT_-_NEATH_PORT_TALBOT
[8] CONWY_-_CONWY
[9] GWYNEDD_-_GWYNEDD
[10] MERTHYR_TUDFUL_-_MERTHYR_TYDFIL
[11] PEN-Y-BONT_AR_OGWR_-_BRIDGEND
[12] POWYS_-_POWYS
[13] RHONDDA_CYNON_TAF_-_RHONDDA_CYNON_TAFF
[14] SIR BENFRO - PEMBROKESHIRE
[15] SIR_BENFRO_-_PEMBROKESHIRE
[16] SIR_CEREDIGION_-_CEREDIGION
[17] SIR_DDINBYCH_-_DENBIGHSHIRE
[18] SIR_FYNWY_-_MONMOUTHSHIRE
[19] SIR_GAERFYRDDIN_-_CARMARTHENSHIRE
[20] SIR_YNYS_MON_-_ISLE_OF_ANGLESEY
[21] SIR_Y_FFLINT_-_FLINTSHIRE
[22] TOR-FAEN_-_TORFAEN
[23] WRECSAM_-_WREXHAM
23 Levels: ABERTAWE_-_SWANSEA ... WRECSAM_-_WREXHAMप्रश्न। मैं shapeआकृति से निकाले गए ऑब्जेक्ट से डेटा का केवल सबसेट कैसे चुन सकता हूं ? उदाहरण के लिए, केवल POWYS_-_POWYSभागों? अगर मैं किसी भी तरह FILE_NAMEसे बनाए गए डेटा फ्रेम में शामिल कर fortifyसकता हूं तो मैं आसानी से barडेटा फ्रेम को कम कर सकता हूं लेकिन मुझे नहीं पता कि यह कैसे करना है। या fortifyवस्तु के केवल हिस्सों को निकालने के लिए उपयोग करने का एक तरीका है ?