用坐标R求网格块数
我试图为提供的坐标找到网格块编号(MeshblockN)。所以我有两个数据集。1.具有唯一id和坐标的数据集A。2.具有网格块编号和其他详细信息的形状文件-数据集B用坐标R求网格块数,r,geospatial,R,Geospatial,我试图为提供的坐标找到网格块编号(MeshblockN)。所以我有两个数据集。1.具有唯一id和坐标的数据集A。2.具有网格块编号和其他详细信息的形状文件-数据集B 到目前为止,我使用readOGR()函数将图形文件-数据集B加载到R中。并加载具有lat和long的数据集A。然后将coordinates()函数应用于数据集A。现在我有了pts.poly所需的两个输入,因为到目前为止没有人回答这个问题,我正在更新自己的版本,我使用sp包中的over函数解决了它超过(A,B) Error in da
到目前为止,我使用
readOGR()
函数将图形文件-数据集B加载到R中。并加载具有lat和long的数据集A。然后将coordinates()函数应用于数据集A。现在我有了pts.poly所需的两个输入,因为到目前为止没有人回答这个问题,我正在更新自己的版本,我使用sp包中的over函数解决了它<代码>超过(A,B)
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