将 GRASS 光栅图读入 R

问题描述

我使用命令 readRAST (rgrass7) 将地图导入 R。现在，它不起作用。我尝试了几个选项，结果为 0。 更新了所有 R 包。安装 GRASS 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

解决方法

rgrass7 的开发者在此回答 github.com/rsbivand/rgrass7/issues/27 并建议现在解决此问题的方法。将在 rgrass7 的未来版本中修复