Menggabungkan data dari CRS/proyeksi yang berbeda
Setiap objek spasial memiliki sistem referensi koordinat (CRS) yang terkait dengannya. Umumnya, ini ditetapkan saat data diimpor dan akan dibaca langsung dari berkas spasial. Inilah cara neighborhoods dan nyc_tracts memperoleh informasi sistem koordinatnya.
Baik paket sp maupun raster memiliki fungsi proj4string() yang mengembalikan CRS dari objek yang dipanggil.
Mencoba bekerja dengan data spasial menggunakan CRS yang berbeda itu mirip dengan bekerja dengan himpunan data yang satu dalam mil dan yang lain dalam kilometer. Keduanya mengukur hal yang sama, tetapi angkanya tidak dapat dibandingkan secara langsung.
Mari kita lihat kedua objek poligon kita.
Latihan ini adalah bagian dari kursus
Memvisualisasikan Data Geospasial di R
Petunjuk latihan
- Panggil
proj4string()padaneighborhoods, lalu ulangi padanyc_tracts. Verifikasi bahwa kedua string tersebut berbeda. - Lihat
head()daricoordinates()milikneighborhoodsdan ulangi untuknyc_tracts. Dapatkah Anda melihat masalahnya?nyc_tractsmemiliki koordinat x sekitar -70, tetapineighborhoodssekitar 1.000.000! - Plot
neighborhoods, lalu plotnyc_tractsdengancol = "red"danadd = TRUEuntuk menambahkannya di atasnya.
Latihan interaktif praktis
Cobalah latihan ini dengan menyelesaikan kode contoh berikut.
library(sp)
# proj4string() on nyc_tracts and neighborhoods
# coordinates() on nyc_tracts and neighborhoods
# plot() neighborhoods and nyc_tracts