আমি ডেনিশ পৌরসভাগুলির জন্য একটি শেফফিলের কয়েকটি সেন্ট্রয়েডগুলি সন্ধান করার চেষ্টা করছি এবং পরে তাদের মধ্যে ড্রাইভিংয়ের সময়টি সন্ধান করব। আমি ব্যবহার Rএর readShapePolyথেকে ফাংশন maptoolsসঙ্গে মিলিত gCentroidথেকে ফাংশন rgeos, এবং সবকিছু কাজ করে। তবে আমি স্প্যাটালপয়েন্টস পাই
SpatialPoints:
x y
1 571860.7 6225016
Coordinate Reference System (CRS) arguments: NA
যা স্পষ্টতই এমন কিছু নয় যা আমি গুগলে ভ্রমণের সময়গুলি ধরতে ব্যবহার করতে পারি। আমি এই সংখ্যাগুলি দ্রাঘিমাংশ-অক্ষাংশে রূপান্তর করার জন্য একটি উপায় খুঁজছি, তবে কীভাবে তা জানা নেই।
যখন আমি ব্যবহার ডেটা পড়তে 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
