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')
