R新手求助：导入ASCII文件生成栅格时行列与范围设置异常

核心错误分析 #

• as.matrix参数无效：as.matrix()函数不支持nrow/ncol参数，你指定的401行、701列不会生效，矩阵行列数完全由输入的df1结构决定，这是栅格行列数未改变的直接原因。
• 栅格行顺序不匹配：多数ASCII地理栅格数据按**北→南（从上到下）的顺序存储，但rast()默认按南→北（从下到上）**读取矩阵，会导致地图错位或翻转。
• 矩阵构造逻辑错误：未强制按指定行列数重构数据，仅通过as.matrix转换无法改变原有结构。

修正代码 #

nrows = 401
ncols = 701
# 读取制表符分隔的ASCII数据
df1 <- read.table("test.txt", header = FALSE, sep = "\t")
# 强制将一维数据转为指定行列的矩阵（按列填充，匹配多数ASCII栅格存储逻辑）
r.mat <- matrix(df1$V1, nrow = nrows, ncol = ncols, byrow = FALSE)
# 反转矩阵行，适配栅格北→南的存储顺序
r.mat <- r.mat[nrow(r.mat):1, ]
# 创建栅格并指定范围（格式为xmin, xmax, ymin, ymax）
r <- rast(r.mat, extent = c(19497.75, 720497.8, 219499.8, 620499.8))
plot(r)

关键说明 #

• 若原始ASCII数据是按行优先存储，需将matrix()中的byrow参数改为TRUE。
• 确认extent参数顺序严格遵循xmin, xmax, ymin, ymax，若原始数据范围定义相反，需对应调整数值顺序。

内容的提问来源于stack exchange，提问作者M_BOKU