निर्देशांक को readShapePoly से R में परिवर्तित कर लंबे-लम्बे निर्देशांक में परिवर्तित करें

Question

मैं डेनिश नगर पालिकाओं के लिए एक आकृति के कुछ केन्द्रक खोजने की कोशिश कर रहा हूं और बाद में उनके बीच ड्राइविंग समय का पता लगाता हूं। मैं उपयोग Rके readShapePolyसे समारोह maptoolsके साथ संयुक्त gCentroidरूप से समारोह rgeos, और सब कुछ काम करता है। हालांकि, मुझे इस तरह के स्थानिक अंक मिलते हैं

SpatialPoints:
     x       y
1 571860.7 6225016
Coordinate Reference System (CRS) arguments: NA 

जो स्पष्ट रूप से ऐसा कुछ नहीं है जिसका उपयोग मैं Google में यात्रा के समय को हड़पने के लिए कर सकता हूं। मैं इन संख्याओं को देशांतर-अक्षांश में बदलने का मार्ग खोज रहा हूं, लेकिन पता नहीं कैसे।

जब मैं लाइब्रेरी readOGRसे उपयोग किए गए डेटा को पढ़ता हूं तो मुझे rgdalसमान निर्देशांक मिलते हैं लेकिन यह मुझे निम्नलिखित के बारे में बताता है कि मैं क्या मान रहा हूं वह प्रक्षेपण है (लेकिन निर्देशांक समान हैं)

Slot "proj4string":
CRS arguments:
+proj=utm +zone=32 +ellps=intl +units=m +no_defs

उदाहरण के लिए: मैंने यहाँ उदाहरण के लिए डेटा रखा है: https://github.com/sebastianbarfort/shapefiles

इस समस्या को पुन: पेश करना चाहिए:

library(maptools)
library(rgdal)
library(rgeos)

map = readShapePoly("~/Downloads/shapefiles-master/kommuner1983.shp")
centroid = gCentroid(map)
centroid

Answer 1

spTransformनिर्देशांक को WGS84 में बदलने के लिए उपयोग करें :

library("rgdal")
library("rgeos")

map <- readOGR(".", "kommuner1983")
map_wgs84 <- spTransform(map, CRS("+proj=longlat +datum=WGS84"))
plot(map_wgs84, axes=TRUE)

gCentroid(map_wgs84)
# SpatialPoints:
#       x     y
# 1 10.05 55.96
# Coordinate Reference System (CRS) arguments: +proj=longlat +datum=WGS84
# +ellps=WGS84 +towgs84=0,0,0 

rgdal::readOGRप्रक्षेपण जानकारी को स्वचालित रूप से पढ़ने में सक्षम है। maptoolsकार्य न तो प्रक्षेपण जानकारी को पढ़ते हैं और न ही लिखते हैं, इन विवरणों को मैन्युअल रूप से प्रबंधित करने के लिए इसे छोड़ दें।

Comments

ये वही है जिसे मैं ढूंढ रहा था।

— sBarfort

@Barfort यदि उत्तर आपके प्रश्न को हल करता है तो आप इसे स्वीकार करना चाहते हैं। इस तरह दूसरों को पता चल जाएगा कि यह सही उत्तर है।

— RK

आपके पास जेंकॉइड के रूप में समान नियंत्रण नहीं है, लेकिन SpatialPolygons वस्तु में केन्द्रक को पकड़ता है। आप उन्हें एक मैट्रिक्स के रूप में पुनः प्राप्त कर सकते हैं: निर्देशांक (map_wgs84)

— जेफरी इवांस