1km วงกลมรอบจุด lat-long- ในหลายสถานที่ในโลก

ฉันมีคะแนน lat-long-points หลายร้อยจุดกระจายอยู่ทั่วโลกและต้องสร้างวงกลมรูปหลายเหลี่ยมรอบแต่ละจุดด้วยรัศมี 1,000 เมตร ฉันเข้าใจว่าจุดต่าง ๆ จะต้องถูกฉายจากองศา (ความยาว lat) เป็นบางสิ่งบางอย่างด้วยหน่วยเมตร แต่จะทำอย่างไรได้โดยไม่ต้องค้นหาด้วยตนเองและกำหนดเขต UTM สำหรับแต่ละจุด

นี่คือความกังวลสำหรับจุดแรกในฟินแลนด์

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"))

แต่ด้วยจุดที่สอง (แคนาดา) มันไม่ทำงาน (เพราะผิดโซน UTM)

the.points.projected <- spTransform(the.points.sp[2, ], CRS( "+init=epsg:32635" ))

วิธีนี้สามารถทำได้โดยไม่ต้องรับด้วยตนเองและระบุจุดโซน UTM ต่อจุด? ฉันไม่มีข้อมูลต่อจุดมากกว่าความยาว lat

ปรับปรุง:

การใช้และการรวมคำตอบที่ยอดเยี่ยมจากทั้ง AndreJ และ Mike T นี่คือรหัสสำหรับทั้งเวอร์ชันและพล็อต พวกเขาต่างกันในทศนิยมลำดับที่ 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

(ไม่แน่ใจว่าจะนำพล็อตเข้าสู่การอัพเดทได้อย่างไร)

คล้ายกับ @AndreJ แต่ใช้การประมาณการแบบgnomicแบบไดนามิกฉันหมายถึงการฉายภาพแบบอะซิมุ ธ ที่มีระยะเวลาเท่ากันแบบไดนามิกเพื่อความแม่นยำยิ่งขึ้น การฉายภาพ AEQ ที่มีศูนย์กลางอยู่ที่แต่ละจุดจะคาดการณ์ระยะทางเท่ากันในทุกทิศทางเช่นวงกลมบัฟเฟอร์ (เครื่องฉายภาพ Mercator จะมีการบิดเบือนในทิศเหนือและทิศตะวันออกเนื่องจากมันล้อมรอบด้านข้างของทรงกระบอก)

ดังนั้นสำหรับจุดแรกของคุณรอบ ๆ ฟินแลนด์สตริง PROJ.4 จะมีลักษณะดังนี้:

+proj=aeqd +lat_0=63.293001 +lon_0=27.472918 +x_0=0 +y_0=0

ดังนั้นคุณสามารถสร้างฟังก์ชั่น R เพื่อสร้างการฉายแบบไดนามิกนี้:

aeqd.buffer <- function(p, r)
{
    stopifnot(length(p) == 1)
    aeqd <- sprintf("+proj=aeqd +lat_0=%s +lon_0=%s +x_0=0 +y_0=0",
                    p@coords[[2]], p@coords[[1]])
    projected <- spTransform(p, CRS(aeqd))
    buffered <- gBuffer(projected, width=r, byid=TRUE)
    spTransform(buffered, p@proj4string)
}

จากนั้นทำสิ่งนี้สำหรับแคนาดา (รายการ 2):

aeqd.buffer(the.points.sp[2,], 1000)

แทนที่จะค้นหาโซน UTM ที่ถูกต้องคุณสามารถสร้างการฉายภาพตามขวางแบบกำหนดเองสำหรับทุกๆจุดด้วย

+proj=tmerc +lat_0=.... +lon_0=... +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

วาดวงกลมในเส้นโครงนั้น พิกัดจุดสุดยอดของวงกลมที่คาดการณ์จะเหมือนกันเสมอดังนั้นคุณต้องสร้างเพียงครั้งเดียว สำหรับสิ่งต่อไปนี้เพียงแค่กำหนด CRS ที่กำหนดเองใหม่ให้กับพวกเขา

Reproject วงกลมเป็น EPSG: 4326 เพื่อใช้งานต่อไป

ภายในระยะ 1,000 เมตรวงกลมจะเกือบแน่นอน ถ้าไม่ได้ (หรือวงกลมขนาดใหญ่) ให้ใช้แทนaeqdtmerc

จะทำอย่างไรถ้าคุณใช้วิธีสร้าง 1,000 เมตรใน EPSG: 4326 รอบ ๆ แต่ละจุดของคุณ ถ้าต้องการแปลง EPSG: 4326 เป็นระบบพิกัดอื่น ๆ ข้อดีของการฉายจุดคือคุณไม่ต้องกังวลกับความโค้งของโลกด้วย EPSG: 4326