使用terra包绘制矢量数据时出现下标越界错误且无坐标轴和图例

1. 先确认矢量数据的属性结构 #

先运行以下代码检查属性表的基本情况，确保第1个属性确实存在且有效：

# 查看属性列的前几行数据
head(values(v))
# 统计属性列的数量
ncol(v)

如果ncol(v)返回0，说明你的矢量数据根本没有属性列，plot(v,1)报错就是因为没有可调用的第1个属性。这种情况要检查原始shp文件的dbf附属文件是否完好，或者重新导出带属性字段的矢量文件。

2. 解决绘图无坐标轴和图例的问题 #

强制显示坐标轴 #

默认调用属性绘图时可能会隐藏坐标轴，手动添加axes=TRUE参数即可：

plot(v, 1, axes=TRUE)

修复图例显示问题 #

plg参数需要传入完整的配置列表，只写位置字符串无效，正确写法如下：

plot(v, 1, axes=TRUE, plg=list(pos="topleft"))

如果还是没图例，先检查第1个属性的取值是否全部相同——用unique(v[[1]])查看，若所有要素属性值一致，图例自然不会生成。

3. 重置绘图设备 #

有时候绘图设备的异常状态会导致显示bug，重置后再尝试绘图：

dev.off()
plot(v, 1, axes=TRUE, plg=list(pos="topleft"))

4. 更新terra包版本 #

旧版本的terra可能存在绘图逻辑的bug，更新到最新版本再测试：

install.packages("terra")
library(terra)

内容的提问来源于stack exchange，提问作者we need a Mat. Stat.