基于R中的坐标合并df 10 x 10 km栅格网格

q35jwt9p  于 2023-04-27  发布在  其他
关注(0)|答案(1)|浏览(82)

1.我希望创建一个df,其中包含10 x10公里的网格在一个指定的区域(坐标边界)在R。
1.然后,我希望从新的df中分配一个网格正方形到一个现有的df,其中包含纬度和经度,年份,存在/不存在和其他数据当前df的样本:
| 纵向|纬度|年份|在场|设置|
| --------------|--------------|--------------|--------------|--------------|
| 67.71150|67.71150|二〇一七年|Y|1|
| 67.72150|67.72150|二〇一七年|N|二|
因此,新的df理想情况下看起来像这样:
| 纵向|纬度|年份|在场|设置|网格|
| --------------|--------------|--------------|--------------|--------------|--------------|
| 67.71150|67.71150|二〇一七年|Y|1|二|
| 67.72150|67.72150|二〇一七年|N|二|1|

wz3gfoph

wz3gfoph1#

你可以做这样的事

library(terra)
r <- rast(res=1/12, xmin=0, xmax=90, ymin=0, ymax=90) 
df$cell <- cellFromXY(r, df[, ("LONGITUDE", "LATITUDE")])

请注意,您谈论的是10 km(平面)分辨率,但您的坐标是角坐标(经度/纬度)。您不能直接匹配不同的坐标参考系统。我创建的栅格也是经度/纬度,它的5分钟分辨率约为5 km(高纬度除外)。您也可以创建平面栅格,并将坐标转换为相同的crs。

相关问题