【问题标题】:RasterLayer in RR中的光栅层
【发布时间】:2014-03-31 19:33:18
【问题描述】:

我有一个包含多个 .adf 文件的 GIS 项目。由于我是 ArcMap 的新手,所以我不知道如何使用它们。我设法使用光栅库将它们读入 r

r <- raster("w001001.adf") 

当查看 sum(r) 时,我得到:

sum(r)
class       : RasterLayer 
dimensions  : 1641, 1358, 2228478  (nrow, ncol, ncell)
resolution  : 0.008333333, 0.008333333  (x, y)
extent      : -47.63578, -36.31911, -20.93489, -7.259889  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
data source : C:\Users\balal\Desktop\blabla\w001001.adf 
names       : w001001 
values      : 1, 15  (min, max)

我试图找到有关单元大小和每个单元的土地利用的信息(我假设土地利用与值相关,因为我有 15 个土地利用类别。)并将这些数据放入我的表格中可以写成.csv。这可能与R吗?还是我在看错误的文件?

抱歉,我不太清楚,但我对 GIS 几乎一无所知。谢谢!

【问题讨论】:

  • 您将受益于首先阅读光栅插图。
  • 分辨率单位取决于坐标系,在我看来它们就像十进制度 (en.wikipedia.org/wiki/Decimal_degrees)。将十进制度转换为距离有点棘手,因为它会随纬度而变化,但你所拥有的大约是 1 x 1 公里。关于问题的第二部分,您确定要将所有 200 万像素值写入 .csv 吗?
  • 查看 ?ratify 以引入与像素整数匹配的类别值。也 ?area 将给出 lonlat 数据的近似大小,但您应该使用 projection ()

标签: r gis arcmap


【解决方案1】:

您始终可以获取矩阵格式的栅格文件(单元格值)中的信息:

r.matrix<-as.matrix(r)

每个单元格的大小由下式给出:

resolution  : 0.008333333, 0.008333333  (x, y)

【讨论】:

    猜你喜欢
    • 2021-03-12
    • 1970-01-01
    • 2016-06-30
    • 1970-01-01
    • 2023-03-17
    • 1970-01-01
    • 2019-04-22
    • 2021-09-16
    相关资源
    最近更新 更多