R语言绘制研究区概况图

要绘制研究区的概况图、土地利用图和DEM（数字高程模型）土壤类型图，你可以使用R语言中的`raster`包和`sp`包来完成。

library(raster)
library(sp)
library(RColorBrewer)

# 读取研究区边界数据
study_area <- readOGR(dsn = "路径", layer = "图层名称")  # 根据实际情况替换路径和图层名称

# 读取土地利用数据
land_use <- raster("路径/土地利用.tif")  # 根据实际情况替换路径和土地利用文件名

# 读取DEM数据
dem <- raster("路径/DEM.tif")  # 根据实际情况替换路径和DEM文件名

# 设置图形设备
dev.new(width = 10, height = 8)  # 自定义图形设备的宽度和高度

# 绘制研究区概况图
plot(study_area, main = "研究区概况图")

# 绘制土地利用图
colormap <- brewer.pal(n = 8, name = "Set3")  # 使用颜色方案设置颜色映射
plot(land_use, col = colormap, main = "土地利用图")

# 绘制DEM土壤类型图
colormap <- terrain.colors(10)  # 使用颜色方案设置颜色映射
plot(dem, col = colormap, main = "DEM土壤类型图")

# 添加图例
par(oma = c(0, 0, 0, 3))  # 调整外边距，给图例留出空间
legend(x = "right", legend = "图例", fill = colormap)

# 关闭图形设备
dev.off()
 

需要将代码中的`路径`、`图层名称`、`土地利用.tif`和`DEM.tif`替换为实际的文件和路径名称。同时，可以根据需要选择适合的颜色方案（使用`brewer.pal`和`terrain.colors`函数），并在绘制图像时进行相应的调整。

这段代码会首先读取研究区的边界数据（以矢量格式存储），然后读取土地利用数据和DEM数据（以栅格格式存储）。接下来，它会创建一个新的图形设备，并分别绘制研究区概况图。