จะแนบ data.frame อย่างง่ายกับ SpatialPolygonDataFrame ใน R ได้อย่างไร

ให้ df = data frame, sp = วัตถุรูปหลายเหลี่ยมเชิงพื้นที่ และโดย = ชื่อหรือหมายเลขคอลัมน์ของคอลัมน์ทั่วไป จากนั้นคุณสามารถรวม data frame เข้ากับ sp object ได้โดยใช้โค้ดบรรทัดต่อไปนี้

sp@data = data.frame(sp@data, df[match(sp@data[,by], df[,by]),])

นี่คือวิธีการทำงานของโค้ด ฟังก์ชันการจับคู่ภายในจะจัดแนวคอลัมน์เพื่อรักษาลำดับไว้ ดังนั้นเมื่อเรารวมมันเข้ากับ sp@data ลำดับจะถูกรักษาไว้อย่างถูกต้อง การตรวจสอบอย่างรวดเร็วเพื่อดูว่าโค้ดใช้งานได้หรือไม่คือการตรวจสอบทั้งสองคอลัมน์ที่สอดคล้องกับคอลัมน์ทั่วไปและดูว่ามันเหมือนกันหรือไม่ (คอลัมน์ทั่วไปมีการทำซ้ำและง่ายต่อการลบสำเนา แต่ฉันเก็บไว้เหมือนเดิม เช็คดีๆ)

มันง่ายอย่างนี้:

require(sp) # the trick is that this package must be loaded!

oo <- merge(SPDF,DF, by="QDGC")

ฉันทดสอบด้วยตัวเองแล้ว แต่จะใช้งานได้ก็ต่อเมื่อคุณใช้ merge from package sp นี่เป็นค่าเริ่มต้นเมื่อโหลดแพ็คเกจ sp ฟังก์ชัน merge จะถูกโอเวอร์โหลด และใช้ sp::merge หากอาร์กิวเมนต์แรกเป็นโครงสร้างเชิงพื้นที่

การผสานสามารถสร้างดาต้าเฟรมที่มีแถวมากกว่าต้นฉบับได้หากไม่มีการแมป 1-1 อย่างง่ายของทั้งสองดาต้าเฟรม ในกรณีนี้ จะต้องคัดลอกเรขาคณิตทั้งหมดและสร้างรูปหลายเหลี่ยมหลายอัน ซึ่งอาจไม่ใช่เรื่องดี

หากคุณมี dataframe ซึ่งมีจำนวนแถวเท่ากันกับ SpatialPointsDataFrame คุณสามารถแทนที่ช่อง @data ได้โดยตรง

library(sp)
example(overlay) # to get the srdf object
srdf@data
spplot(srdf)
srdf@data=data.frame(x=runif(3),xx=rep(0,3))
spplot(srdf)

หากคุณนับจำนวนแถวผิด:

srdf@data=data.frame(x=runif(2),xx=rep(0,2))
spplot(srdf)
Error in data.frame(..., check.names = FALSE) : 
  arguments imply differing number of rows: 3, 2

บางทีฟังก์ชัน joinCountryData2Map ในแพ็คเกจ rworldmap อาจเป็นแรงบันดาลใจได้ (แต่ผมอาจจะผิดเหมือนครั้งที่แล้วนะครับ)

อีกวิธีหนึ่งคือการใช้ฟังก์ชัน append_data จากแพ็คเกจ tmaptools . มันถูกเรียกด้วยข้อโต้แย้งเหล่านี้:

append_data(shp, data, key.shp = NULL, key.data = NULL,
  ignore.duplicates = FALSE, ignore.na = FALSE,
  fixed.order = is.null(key.data) && is.null(key.shp))

น่าเสียดายที่มันถูกเรียกว่า append เนื่องจากฉันเข้าใจว่าการผนวกมากกว่านั้นมีความหมายว่า rbind และเราต้องการมีบางอย่างเช่น join หรือ merge ที่นี่

หากเพิกเฉยต่อข้อเท็จจริงดังกล่าว ฟังก์ชันจะมีประโยชน์มากในการทำให้แน่ใจว่าคุณได้การรวมที่ถูกต้อง และหากมีบางแถวปรากฏเพียงด้านเดียวของการรวม จากเอกสาร:

ภายใต้ความครอบคลุม (รายการรูปร่างที่ไม่สอดคล้องกับบันทึกข้อมูล) ความครอบคลุมมากเกินไป (บันทึกข้อมูลที่ไม่สอดคล้องกับรายการรูปร่างตามลำดับ) รวมถึงการมีอยู่ของค่าคีย์ที่ซ้ำกันจะถูกตรวจสอบและรายงานโดยอัตโนมัติผ่านข้อความคอนโซล ด้วย under_coverage และ over_coverage ค่าคีย์ความครอบคลุมด้านล่างและด้านบนจากการเรียก append_data ครั้งล่าสุดสามารถดึงข้อมูลได้

หากเป็นเชปไฟล์สองไฟล์ที่จำเป็นต้องรวมเป็นออบเจ็กต์เดียว ให้ใช้ rbind()

เมื่อใช้ rbind() โปรดตรวจสอบให้แน่ใจว่าอาร์กิวเมนต์ทั้งสองที่คุณใช้คือ SpatialDataFrames คุณสามารถตรวจสอบได้โดยใช้ class(sf) หากไม่ใช่ dataframe ให้ใช้ st_as_sf() เพื่อแปลงเป็น SpatialDataFrame ก่อนที่จะทำการ rbind

หมายเหตุ : คุณยังสามารถใช้สิ่งนี้เพื่อต่อท้าย NULLs ได้ โดยเฉพาะเมื่อคุณใช้ผลลัพธ์จากการวนซ้ำและต้องการสะสมผลลัพธ์