如何将R模型应用于一组光栅和矢量预测器图层

如何将R模型应用于一组光栅和矢量预测器图层,r,vector,spatial,raster,grass,R,Vector,Spatial,Raster,Grass,我有一个R模型,它预测了大约6个预测值。我有一组与预测器相对应的空间图层,一些作为栅格,一些作为矢量要素类,我希望通过将模型应用于每个输入栅格单元来生成预测值的栅格 (我在编程ArcPy和Python方面很有经验,但对于R来说就不那么有经验了,因为我没有使用Grass) 在R中,我是否应该遍历输出光栅的所有单元格,找到所有相关的预测值,然后调用模型并将结果写入输出光栅?或者有没有更有效的方法,比如使用光栅堆栈 哪些R包对此有用?我应该用草吗 其中两个输入光栅的大小和对齐方式相同,而第三个光栅的单

我有一个R模型,它预测了大约6个预测值。我有一组与预测器相对应的空间图层,一些作为栅格,一些作为矢量要素类,我希望通过将模型应用于每个输入栅格单元来生成预测值的栅格

(我在编程ArcPy和Python方面很有经验,但对于R来说就不那么有经验了,因为我没有使用Grass)

  • 在R中,我是否应该遍历输出光栅的所有单元格,找到所有相关的预测值,然后调用模型并将结果写入输出光栅?或者有没有更有效的方法,比如使用光栅堆栈

  • 哪些R包对此有用?我应该用草吗

  • 其中两个输入光栅的大小和对齐方式相同,而第三个光栅的单元大小较大。在开始循环之前,是否需要以某种方式转换第三个光栅以匹配其他光栅?我该怎么做

  • 作为向量特征类的预测器层由python脚本生成,我希望从同一个R脚本调用该脚本。我可以将它们作为光栅获得(但这将是一个缓慢的过程),或者:

    • 我可以在R-how中将它们转换为光栅吗
    • 我可以照原样用吗?如何从位置的矢量要素类中获取值(对应于光栅单元的中心)

  • 最后,R是最好的地方吗?还是草?还是Python?

    光栅软件包将是您解决该问题的首选方案

  • 使用
    光栅化
    将矢量转换为光栅。为此,您可以访问以下QA链接

  • 将所有光栅转换为相同的投影和像素分辨率。请访问以下网站链接

  • 将所有光栅堆叠在下面

  • 将包含纬度、经度和因变量的点向量图层覆盖在光栅堆栈上,并提取光栅值

  • 使用相应的R包应用任何回归模型

  • 使用
    raster::predict
    函数将模型应用于光栅图层范围,使用光栅堆栈作为自变量,如中所示


  • 这是一个非常广泛的问题,有很多组成部分-你应该在R中找到一个关于空间数据的教程。如果你在任何特定的问题上卡住了,那么回到这里,问一个更集中的问题。