将 GRASS 栅格地图读入 R答案

【问题标题】：reading GRASS raster maps into R将 GRASS 栅格地图读入 R
【发布时间】：2021-07-05 04:52:27
【问题描述】：

我使用命令 readRAST (rgrass7) 将地图导入 R。现在，它不起作用。我尝试了几个选项，结果为 0。更新了所有 R 包。安装了草 7.8.5。在 GRASS 中创建 Rgui。获取

    elev <- readRAST("elevation", ignore.stderr=TRUE)
    Error in sp::CRS(getLocationProj()) : NA

我仍然可以毫无问题地通过 readVECT 获取矢量图。

有什么办法可以解决吗？如果没有，目前将 GRASS 栅格地图导入 R 的方法是什么？

谢谢。

#================

R version 4.0.5 (2021-03-31) -- "Shake and Throw"
Copyright (C) 2021 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64 (64-bit)

During startup - Warning messages:
1: Setting LC_CTYPE=en_AU.cp1252 failed
2: Setting LC_COLLATE=en_AU.cp1252 failed
3: Setting LC_TIME=en_AU.cp1252 failed
4: Setting LC_MONETARY=en_AU.cp1252 failed

    library(sp)
    library(rgrass7)
    Loading required package: XML
    GRASS GIS interface loaded with GRASS version: GRASS 7.8.5 (2020)
    and location: nc_spm_08_grass7
    library(rgdal)
    rgdal: version: 1.5-23, (SVN revision 1121)
    Geospatial Data Abstraction Library extensions to R successfully loaded
    Loaded GDAL runtime: GDAL 3.2.1, released 2020/12/29
    Path to GDAL shared files: C:/Users/andrew/Documents/R/win-library/4.0/rgdal/gdal
    GDAL binary built with GEOS: TRUE
    Loaded PROJ runtime: Rel. 7.2.1, January 1st, 2021, [PJ_VERSION: 721]
    Path to PROJ shared files: C:/Users/andrew/Documents/R/win-library/4.0/rgdal/proj
    PROJ CDN enabled: FALSE
    Linking to sp version:1.4-5
    To mute warnings of possible GDAL/OSR exportToProj4() degradation,
    use options("rgdal_show_exportToProj4_warnings"="none") before loading rgdal.
    Overwritten PROJ_LIB was C:/Users/andrew/Documents/R/win-library/4.0/rgdal/proj
    gmeta()
    gisdbase C:\andrew\GISDATABASE
    location nc_spm_08_grass7
    mapset PERMANENT
    rows 280
    columns 767
    north 307000
    south 27000
    west 151000
    east 918000
    nsres 1000
    ewres 1000
    use_sp()

    sessionInfo()
    R version 4.0.5 (2021-03-31)
    Platform: x86_64-w64-mingw32/x64 (64-bit)
    Running under: Windows 10 x64 (build 18363)


    Matrix products: default
    
    locale:
    [1] LC_COLLATE=C LC_CTYPE=English_Australia.1252
    [3] LC_MONETARY=C LC_NUMERIC=C
    [5] LC_TIME=C
    
    attached base packages:
    [1] stats graphics grDevices utils datasets methods base
    
    other attached packages:
    [1] rgdal_1.5-23 rgrass7_0.2-6 XML_3.99-0.6 sp_1.4-5
    
    loaded via a namespace (and not attached):
    [1] compiler_4.0.5 grid_4.0.5 lattice_0.20-41
    > elev <- readRAST("elevation", ignore.stderr=TRUE)
    Error in sp::CRS(getLocationProj()) : NA

【问题讨论】：

您可以使用raster 包来读取像raster("filename.extension") 这样的栅格。
问题不在于只读取光栅文件名.extension，而是关于正确读取作为 GRASS 一部分的光栅文件（包含所有包含的信息）。

标签： r raster grass

【解决方案1】：

rgrass7 的开发者在此处github.com/rsbivand/rgrass7/issues/27 进行了回答，并提出了现在解决此问题的方法。它会在 rgrass7 的未来版本中修复

【讨论】：