Saya mencoba untuk menemukan beberapa centroid dari shapefile untuk kota Denmark dan kemudian menemukan waktu mengemudi di antara mereka. Saya menggunakan R
's readShapePoly
fungsi dari maptools
dikombinasikan dengan gCentroid
fungsi dari rgeos
, dan semuanya bekerja. Namun, saya mendapatkan titik spasial seperti
SpatialPoints:
x y
1 571860.7 6225016
Coordinate Reference System (CRS) arguments: NA
Yang jelas bukan sesuatu yang bisa saya gunakan di Google untuk mengambil waktu perjalanan. Saya mencari cara untuk mengubah angka-angka ini menjadi bujur-lintang, tetapi tidak tahu caranya.
Ketika saya membaca data menggunakan readOGR
dari rgdal
perpustakaan saya mendapatkan koordinat yang sama tetapi ini memberi tahu saya berikut tentang apa yang saya asumsikan adalah proyeksi (tetapi koordinatnya sama)
Slot "proj4string":
CRS arguments:
+proj=utm +zone=32 +ellps=intl +units=m +no_defs
Contoh yang dapat direproduksi: Saya telah memasukkan data untuk contoh di sini: https://github.com/sebastianbarfort/shapefiles
Ini seharusnya mereproduksi masalah:
library(maptools)
library(rgdal)
library(rgeos)
map = readShapePoly("~/Downloads/shapefiles-master/kommuner1983.shp")
centroid = gCentroid(map)
centroid
sumber