2021-12-14 在R中 提取坐标点对应的栅格值

1# 获取点坐标 坐标格式：

保存csv

2#代码 

library(raster)

library(rgdal)

#导入坐标文件

coor= read.csv(choose.files(),header = T)

#导入地图

raster=raster(choose.files())

plot(coor$LATITUDE ,coor$LONGITUDE , main="Point_ID")

#将坐标数据集更改为 SpatialPoints 格式

locationsC <- SpatialPoints(cbind(coor$LATITUDE ,coor$LONGITUDE),

                            proj4string=CRS("+init=epsg:4326"))

cord.UTM <- spTransform(locationsC, crs(raster))

#提取成功

res <- extract(raster, cord.UTM, method='simple')