Ermitteln der Maschenblocknummer unter Verwendung der Koordinaten R

Question

Ich versuche die Netzblocknummer (MeshblockN) für die angegebenen Koordinaten zu finden. Also habe ich zwei Datensätze. 1. Datensatz A mit eindeutiger ID und Koordinaten. 2. Formdatei mit Mesh-Blocknummer und weiteren Details - Datensatz B.

Was ich bisher gemacht habe, geladen die Shape-Datei - Datensatz B in R mit readOGR() Funktion. Und lud den Datensatz A, der lat und long hat. Dann angewendet coordinates() Funktion auf Datensatz A. Jetzt habe ich die zwei erforderlichen Eingaben für pts.poly <- point.in.poly(A,B). Aber ich bin immer diesen Fehler -

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 

Kann mir jemand helfen, gehe ich vor? Ich bin mir nicht sicher, was daran falsch ist.

Answer 1

Da niemand so weit geantwortet hat, aktualisiere ich meine eigene Version, löste ich es über Funktion in SP-Paket. over(A,B)