從緯度/經度坐標中減去數據

Question

我有2個看起來像這樣的數據文件：

Model Data
long  lat   count
96.25 18.75     4
78.75 21.25     3
86.75 23.25     7
91.25 33.75     10

Observation Data
long  lat   count
96.75 25.75    10
86.75 23.25     7
78.75 21.25    11
95.25 30.25     5

我試圖減去匹配的經緯度組合（模型數據-觀測數據）的計數，以使第一個組合78.75和21.25的差值為-8。 任何不匹配要減去的經/緯度點都將被0減去或從中減去。

我已經嘗試過if語句來匹配減法點：

if (modeldata$long == obsdata$long & modeldata$lat == obsdata$lat) {
      obsdata$difference <- modeldata$count - obsdata$count
 }

但是，除非匹配點碰巧落在同一行內，否則這只會按順序減去行，而不是按匹配點。

我也收到以下警告：

警告信息：

1：在modeldata $ long == obsdata $ long中：較長的對象長度不是較短的對象長度的倍數

2：在modeldata $ lat == obsdata $ lat中：較長的對象長度不是較短的對象長度的倍數

3：如果if（modeldata $ long == obsdata $ long＆modeldata $ lat ==：條件的長度> 1，並且僅使用第一個元素

任何幫助將不勝感激！

Answer 1

您可以合並坐標，為NA加0並減去。

mdl <- read.table(text = "long  lat   count
96.25 18.75     4
78.75 21.25     3
86.75 23.25     7
91.25 33.75     10", header = TRUE)

obs <- read.table(text = "long  lat   count
96.75 25.75    10
                  86.75 23.25     7
                  78.75 21.25    11
                  95.25 30.25     5", header = TRUE)

xy <- merge(mdl, obs, by = c("long", "lat"), all.x = TRUE)

xy[is.na(xy)] <- 0

xy$diff <- xy$count.x - xy$count.y
xy

   long   lat count.x count.y diff
1 78.75 21.25       3      11   -8
2 86.75 23.25       7       7    0
3 91.25 33.75      10       0   10
4 96.25 18.75       4       0    4

Answer 2

您可以使用data.table加入並更新

library(data.table)

## reading your supplied data
# dt_model <- fread(
#   'long  lat   count
# 96.25 18.75     4
#   78.75 21.25     3
#   86.75 23.25     7
#   91.25 33.75     10'
# )
# 
# 
# dt_obs <- fread(
# "long  lat   count
#   96.75 25.75    10
#   86.75 23.25     7
#   78.75 21.25    11
#   95.25 30.25     5"
# )

setDT(dt_model)
setDT(dt_obs)

## this join & update will update the `dt_model`. 
dt_model[ 
    dt_obs
    , on = c("long", "lat")
    , count := count - i.count
    ]

dt_model
#     long   lat count
# 1: 96.25 18.75     4
# 2: 78.75 21.25    -8
# 3: 86.75 23.25     0
# 4: 91.25 33.75    10

注意一個明顯的警告，即加入坐標（浮點數/小數）可能並不總是能給出正確的答案

Answer 3

這是dplyr的選項

library(dplyr)
left_join(mdl, obs, by = c("long", "lat")) %>%
   transmute(long, lat, count = count.x - replace(count.y, is.na(count.y), 0))
#   long   lat count
#1 96.25 18.75     4
#2 78.75 21.25    -8
#3 86.75 23.25     0
#4 91.25 33.75    10