Menyeleksi subset neighborhoods
Anda tidak memerlukan semua wilayah lingkungan lain di New York, jadi Anda akan mengambil subset hanya lingkungan di New York County. Anda sudah tahu caranya!
neighborhoods adalah SpatialPolygonsDataFrame dan Anda telah mempelajari di Bab 2 cara mengambil subset berdasarkan kolom di slot data. Kuncinya adalah membuat vektor logika, lalu mengambil subset SpatialPolygonsDataFrame seperti mengambil subset data frame.
Bagaimana cara mengidentifikasi lingkungan yang tepat? Coba jalankan:
head(neighborhoods@data)
CountyFIPS adalah kode numerik yang mengidentifikasi county. Jika Anda dapat mengetahui kode untuk New York County, Anda dapat menyisakan hanya baris dengan nilai tersebut.
Latihan ini adalah bagian dari kursus
Memvisualisasikan Data Geospasial di R
Petunjuk latihan
- Objek
nyc_tracts_mergejuga memiliki kode county pada kolomCOUNTYFP. Cari nilaiunique()untuk menemukan kode New York County. - Ambil subset
neighborhoodsdengan menambahkan ekspresi logika yang menguji apakahneighborhoods$CountyFIPSmemiliki nilai yang benar. - Sunting plot Anda untuk menggunakan
manhat_hoodsalih-alih neighborhoods. - Tambahkan layer
tm_text(), memetakantextke"NTAName".
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
library(tmap)
# Find unique() nyc_tracts_merge$COUNTYFP
# Add logical expression to pull out New York County
manhat_hoods <- neighborhoods[___, ]
tm_shape(nyc_tracts_merge) +
tm_fill(col = "estimate") +
tm_shape(water) +
tm_fill(col = "grey90") +
# Edit to use manhat_hoods instead
tm_shape(neighborhoods) +
tm_borders() +
# Add a tm_text() layer