代码之家 › 专栏 › 技术社区 › GCGM

在R中创建交互式绘图

shapefile raster plot ggplot2 r

GCGM · 技术社区 · 6 年前

我画了一个光栅

plot(NDVI[[4]])

[![在此处输入图像描述][1][1]

之后,我添加了第二个包含多边形限制的层:

plot(fields, add=TRUE)

我的目标是用户点击图像来选择其中的一些多边形(比如3个)。这些单击用于标识这些多边形,这些多边形稍后将用于导出它们所代表区域内的平均光栅值。

到目前为止,我一直在更新一个包含点的shapefile,但我想让它更具交互性

points<-readOGR("Points_crops.shp") 
fields<-readOGR("Boundaries.shp")
fields_sub <- fields[!is.na(sp::over(fields, sp::geometry(points))), ]
NDVI_mean<-lapply(NDVI, FUN=function (NDVI) {data.frame(mean=extract(NDVI,fields_sub,fun=mean))})

1 回复 | 直到 6 年前

GCGM 6 年前

对于那些感兴趣的人,我已经使用 click spatialPoints 设定 proj4string (与参考光栅相同)。

  points<-click(NDVI[[4]], n=5, xy=TRUE, show=TRUE)
  points$value<-NULL
  points<-SpatialPoints(points, proj4string = crs(S2_stack_crop[[2]]))

一旦我有了这些要点 空间点 我可以继续下一步

推荐文章

cezar Dennis Kioko · 使用python从shapefile中删除记录

7 年前

ShaiNe Ram · 如何使用JAVA代码读取形状(.shp)文件并显示其值

7 年前

Alex Blasco · 如何获得形状文件中多边形的所有点。shp使用GeoTools?

7 年前

tda · 使用shapefile[duplicate]在Python中查找点的封闭多边形

7 年前

Jereme · 如何将形状文件导入SQL Server?

7 年前

IsisDorus · 将Lat/Lon点映射到R中的形状文件

7 年前

maurobio · 使用pyshp创建直线段的形状文件?

7 年前

Chris · 计算多个多边形之间共享边界的长度

8 年前

27titanik · 如何将shp文件中的数据与数据帧中包含的数据进行匹配?

9 年前

Fabio Lamanna · python-创建等间距空间网格的shapefile

10 年前