【问题标题】:Export multiple rasters from GRASS从 GRASS 导出多个栅格
【发布时间】:2014-08-08 16:16:34
【问题描述】:

我正在使用 GRASS 和 R 来分析一些高光谱数据。 我想将我的光谱层导出到一个与 ENVI 兼容的 .img 文件。 到目前为止,我已经使用 g.mlist 选择了合适的层,现在我使用 i.group 从这些层创建组,我可以使用 r.out.gdal 导出。

我看到 i.group 接受一个输入参数,即栅格图层名称列表 - 我有大约 100 个,所以我正在寻找一种方法来避免单独输入它们!

我的问题本质上是,我已经将名称从 g.mlist 导出到一个文本文件 - 有没有办法在输入参数中使用这个文件?

这可能是一个非常基本的问题,并且是由于我的编程经验不足。

我想使用类似的东西:

system("i.group group=MSP1 input=/home/Documents/MSP/list.txt")

除了 /home/Documents/MSP/list.txt 不是栅格地图。

或者,如果有人知道将多个图像层导出到一个文件中的更好方法,那也很棒。

【问题讨论】:

    标签: r spatial grass


    【解决方案1】:

    我没有像我发布的问题那样尝试从我的 Rscript 中调用 GRASS,而是从单独的脚本中调用了 i.group,它似乎可以工作。 所以我现在的电话是:

    Rscript MSP.R
    sh ExportENVI.sh
    

    ExportENVI 脚本是:

    #!/bin/bash
    #Exports corrected layers to single image file for ENVI
    
    group_list=$(g.mlist type=rast pattern=msp* separator=,)
    echo $group_list
    i.group group=MSP1 input=$group_list
    r.out.gdal input=MSP1 output=/home/Documents/MSP/group_output.img format=HFA
    r.out.gdal input=MSP1 output=/home/Documents/MSP/group_output.img format=ENVI
    exit 0
    

    我不确定 ENVI 的正确格式,如果我需要在这里同时做 .img 和 .hdr,但我会继续工作。

    【讨论】:

      猜你喜欢
      • 2021-07-05
      • 2017-05-02
      • 1970-01-01
      • 1970-01-01
      • 2015-09-13
      • 2019-04-02
      • 2019-01-19
      • 1970-01-01
      • 2017-01-26
      相关资源
      最近更新 更多