座標を使用してメッシュブロック番号を検索するR

Question

指定された座標に対してメッシュブロック番号（MeshblockN）を検索しようとしています。だから私は2つのデータセットを持っています。ユニークなIDおよび座標を有するデータセットA。 2.メッシュブロック番号とその他の詳細を含む形状ファイル - データセットB.

これまでに行ったことは、readOGR()関数を使用してRにデータセットBをロードしました。そして、緯度と経度のあるデータセットAをロードしました。次に、coordinates（）関数をデータセットAに適用しました。今では、pts.poly <- point.in.poly(A,B)の2つの必須入力があります。しかし、私はこのエラーが発生しています -

Error in data.frame(z@data, stats::na.omit(sp::over(pts, polys))) : 
    arguments imply differing number of rows: 1, 0 

> head(A@data) 
    row_names unique_id 
1   1 2837061 
> head(A@coords) 
    longitude latitude 
1 170.5065 -45.86696 

> head(B@data, n = 1) 
    Meshblock MeshblockN AreaUnitCo AreaUnitNa  
0 MB 0352700 0352700  510402 Chelsea 



> head(B@polygons, n = 1) 
[[1]] 
An object of class "Polygons" 
Slot "Polygons": 
[[1]] 
An object of class "Polygon" 
Slot "labpt": 
[1] 174.72094 -36.82045 

Slot "area": 
[1] 2.858646e-05 

Slot "hole": 
[1] FALSE 

Slot "ringDir": 
[1] 1 

Slot "coords": 
      [,1]  [,2] 
    [1,] 174.7181 -36.81820 
    [2,] 174.7184 -36.81831 
    [3,] 174.7188 -36.81824 
    . 
    . 
[500,] 174.7241 -36.82205 
[ reached getOption("max.print") -- omitted 363 rows ] 



Slot "plotOrder": 
[1] 1 

Slot "labpt": 
[1] 174.72094 -36.82045 

Slot "ID": 
[1] "0" 

Slot "area": 
[1] 2.858646e-05 

誰でも手伝ってもらえますか？これについて何が間違っているか分かりません。

Answer 1

私は自分のバージョンを更新していますので、spパッケージのover関数を使って解決しました。 over(A,B)