Mulai sekarangMulai gratis

Mengonversi dari satu CRS/proyeksi ke yang lain

Proses mengonversi dari satu CRS atau proyeksi ke lainnya ditangani oleh metode spTransform() dalam paket rgdal. spTransform() memiliki metode untuk semua objek sp termasuk SpatialPolygonsDataFrame, tetapi tidak berfungsi pada objek raster. Ini karena transformasi raster sedikit lebih rumit; grid persegi panjang yang ditransformasikan tidak lagi berbentuk persegi panjang. Anda dapat melihat ?raster::projectRaster jika Anda ingin tahu tentang cara mentransformasikan raster.

Transformasi itu sederhana. Argumen pertama untuk spTransform(), x, adalah objek spasial yang akan ditransformasikan dan argumen kedua, CRS, adalah spesifikasi CRS yang diinginkan. CRS dapat ditentukan dengan string PROJ4, yang bisa Anda susun secara manual, tetapi jauh lebih mudah mengambilnya dari objek yang sudah ada (misalnya dengan fungsi proj4string()).

Saatnya membuat kedua himpunan data poligon Anda berada dalam CRS yang sama.

Latihan ini merupakan bagian dari kursus

Memvisualisasikan Data Geospasial di R

Lihat Kursus

Instruksi latihan

  • Transformasikan neighborhoods agar memiliki CRS yang sama dengan nyc_tracts menggunakan spTransform() dengan argumen CRS disetel ke proj4string(nyc_tracts).
  • Verifikasi transformasi dengan melihat head() dari coordinates(neighborhoods).
  • Periksa apakah himpunan data kini selaras dengan memplot neighborhoods, lalu memplot nyc_tracts dengan add = TRUE dan col = "red", dan terakhir memplot water dengan add = TRUE dan col = "blue".

Latihan interaktif langsung praktik

Cobalah latihan ini dengan melengkapi kode contoh ini.

library(sp)
library(raster)

# Use spTransform on neighborhoods: neighborhoods



# head() on coordinates() of neighborhoods


# Plot neighborhoods, nyc_tracts and water
Edit dan Jalankan Kode