วิธีสร้าง Linestrings จากจุดอย่างมีประสิทธิภาพ?

Dec 30 2020

ฉันมี geom POINT ในกรอบข้อมูลสองเฟรมที่แยกจากกัน สิ่งที่ฉันต้องการทำคือการเชื่อมต่อจุดด้วยเส้น (ในภายหลังบนแผนที่) ดังนั้นฉันจึงต้องการสร้าง Linestring สำหรับแต่ละคู่ของจุดจากกรอบข้อมูลเหล่านั้น ฉันทำแบบนี้:

coordsCust <- table %>%
  st_as_sf(coords = c("lonCust","latCust"), crs = 4326)

coordsApp <- table %>%
  st_as_sf(coords = c("lonApp","latApp"), crs = 4326) %>%
  st_geometry()

และ Linestring:

lines <- st_sfc(mapply(function(a,b){
  st_cast(st_union(a,b),"LINESTRING")}, 
  coordsCust$geometry, coordsApp$geometry, SIMPLIFY=FALSE))

รหัสนี้ใช้งานได้ฉันสามารถสร้าง Linestrings สำหรับจุดแต่ละคู่ทีละแถว:

LINESTRING (14.035 51.65182, 14.33418 53.53346)
LINESTRING (20.42767 49.98073, 16.62978 52.31037)
LINESTRING (20.18762 50.03337, 16.62978 52.31037)
LINESTRING (19.04625 49.79234, 16.62978 52.31037)
LINESTRING (21.35808 50.92382, 16.62978 52.31037)

ปัญหาคือสำหรับ 30,000 แถวโซลูชันนี้ทำงานได้ช้ามาก - ประมาณ 21 วินาที มีวิธีอื่นในการสร้าง linestrings จากจุดหรือไม่? สิ่งที่ทำงานได้เร็วขึ้นมาก? ฉันค้นหาวิธีแก้ปัญหาบางอย่างบนเว็บ แต่ก็ไร้ผล ฉันได้อ่านบางอย่างเกี่ยวกับการแปลง sf เป็นเมทริกซ์และการใช้งานpmapแต่ไม่รู้ว่าจะใช้งานได้อย่างไรที่นี่

UPDATE:ถ้าฉันต้องการใช้ฟังก์ชัน sfheaders :: sf_linestring ฉันต้องเข้าร่วมรูปทรงเรขาคณิตจากทั้งสองชุดข้อมูล ฉันทำแบบนี้:

df <- cbind(coordsCust,coordsApp)

และกรอบข้อมูลสุดท้าย (ฉันแสดงส่วนที่สำคัญที่สุดของมัน) แสดงอยู่ด้านล่าง:

ขออภัย sf_linestring ทำงานไม่ถูกต้องบนดาต้าเฟรมนี้ ฉันต้องการสร้าง linestring ระหว่าง POINT สำหรับแต่ละแถวแยกกันดังที่แสดงบนหน้าจอ

คำตอบ

2 SymbolixAU Dec 30 2020 at 04:49

หากไม่มีชุดข้อมูล exmaple ก็ยากที่จะตอบคำถามของคุณได้อย่างสมบูรณ์ แต่ถ้าคุณสามารถทำให้ data.frame ของคุณอยู่ในรูปแบบ 'ยาว' คุณsfheadersสามารถทำได้ในทันที