আমার শত শত ল্যাট-লম্বা-পয়েন্ট রয়েছে যা সারা বিশ্বে ছড়িয়ে আছে এবং অবশ্যই একেবারে 1000 মিটার ব্যাসার্ধের সাথে অবশ্যই তাদের প্রত্যেকের চারদিকে বৃত্ত-বহুভুজ তৈরি করতে হবে। আমি বুঝতে পারি পয়েন্টগুলি প্রথমে ডিগ্রি (ল্যাট লম্বা) থেকে মিটার-ইউনিটগুলির সাথে কোনও কিছুর দিকে প্রক্ষেপণ করা উচিত, তবে প্রতিটি পয়েন্টের জন্য ইউটিএম-জোনগুলি ম্যানুয়ালি অনুসন্ধান এবং সংজ্ঞা না দিয়ে কীভাবে এটি করা যেতে পারে?
এখানে ফিনল্যান্ডের প্রথম পয়েন্টের জন্য একটি মেগাও রয়েছে।
library(sp)
library(rgdal)
library(rgeos)
the.points.latlong <- data.frame(
Country=c("Finland", "Canada", "Tanzania", "Bolivia", "France"),
lat=c(63.293001, 54.239631, -2.855123, -13.795272, 48.603949),
long=c(27.472918, -90.476303, 34.679950, -65.691146, 4.533465))
the.points.sp <- SpatialPointsDataFrame(the.points.latlong[, c("long", "lat")], data.frame(ID=seq(1:nrow(the.points.latlong))), proj4string=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))
the.points.projected <- spTransform(the.points.sp[1, ], CRS( "+init=epsg:32635" )) # Only first point (Finland)
the.circles.projected <- gBuffer(the.points.projected, width=1000, byid=TRUE)
plot(the.circles.projected)
points(the.points.projected)
the.circles.sp <- spTransform(the.circles.projected, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))তবে দ্বিতীয় পয়েন্ট (কানাডা) দিয়ে এটি কাজ করে না (কারণ ভুল ইউটিএম-অঞ্চল)।
the.points.projected <- spTransform(the.points.sp[2, ], CRS( "+init=epsg:32635" ))ম্যানুয়ালি ইউটিএম-জোন পয়েন্ট প্রতি বিন্দু না পেয়ে এবং নির্দিষ্ট না করে কীভাবে এটি করা যায়? আমার কাছে লম্বা লম্বা চেয়ে পয়েন্ট প্রতি আর কোনও তথ্য নেই।
হালনাগাদ:
AndreJ এবং মাইক টি উভয়ের দুর্দান্ত উত্তরগুলি ব্যবহার এবং একত্রিত করার জন্য, এখানে উভয় সংস্করণ এবং প্লটগুলির কোড। তারা 4 র্থ দশমিক বা তার চেয়ে পৃথক, তবে উভয়ই খুব ভাল উত্তর!
gnomic.buffer <- function(p, r) {
stopifnot(length(p) == 1)
gnom <- sprintf("+proj=gnom +lat_0=%s +lon_0=%s +x_0=0 +y_0=0",
p@coords[[2]], p@coords[[1]])
projected <- spTransform(p, CRS(gnom))
buffered <- gBuffer(projected, width=r, byid=TRUE)
spTransform(buffered, p@proj4string)
}
custom.buffer <- function(p, r) {
stopifnot(length(p) == 1)
cust <- sprintf("+proj=tmerc +lat_0=%s +lon_0=%s +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs",
p@coords[[2]], p@coords[[1]])
projected <- spTransform(p, CRS(cust))
buffered <- gBuffer(projected, width=r, byid=TRUE)
spTransform(buffered, p@proj4string)
}
test.1 <- gnomic.buffer(the.points.sp[2,], 1000)
test.2 <- custom.buffer(the.points.sp[2,], 1000)
library(ggplot2)
test.1.f <- fortify(test.1)
test.2.f <- fortify(test.2)
test.1.f$transf <- "gnomic"
test.2.f$transf <- "custom"
test.3.f <- rbind(test.1.f, test.2.f)
p <- ggplot(test.3.f, aes(x=long, y=lat, group=transf))
p <- p + geom_path()
p <- p + facet_wrap(~transf)
p(কীভাবে আপডেটটিতে প্লট পাবেন তা নিশ্চিত নন)।